JP7080992B2 - 欠陥の発生を推定するシステム、及びコンピューター可読媒体 - Google Patents

欠陥の発生を推定するシステム、及びコンピューター可読媒体 Download PDF

Info

Publication number
JP7080992B2
JP7080992B2 JP2020556468A JP2020556468A JP7080992B2 JP 7080992 B2 JP7080992 B2 JP 7080992B2 JP 2020556468 A JP2020556468 A JP 2020556468A JP 2020556468 A JP2020556468 A JP 2020556468A JP 7080992 B2 JP7080992 B2 JP 7080992B2
Authority
JP
Japan
Prior art keywords
edge
probability
defect
data
wafer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020556468A
Other languages
English (en)
Other versions
JPWO2020100180A1 (ja
Inventor
宏 福田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi High Tech Corp
Original Assignee
Hitachi High Tech Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi High Tech Corp filed Critical Hitachi High Tech Corp
Publication of JPWO2020100180A1 publication Critical patent/JPWO2020100180A1/ja
Application granted granted Critical
Publication of JP7080992B2 publication Critical patent/JP7080992B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/06Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and measuring the absorption
    • G01N23/18Investigating the presence of flaws defects or foreign matter
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70491Information management, e.g. software; Active and passive control, e.g. details of controlling exposure processes or exposure tool monitoring processes
    • G03F7/705Modelling or simulating from physical phenomena up to complete wafer processes or whole workflow in wafer productions
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70605Workpiece metrology
    • G03F7/70616Monitoring the printed patterns
    • G03F7/7065Defects, e.g. optical inspection of patterned layer for defects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/0006Industrial image inspection using a design-rule based approach
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01LSEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
    • H01L22/00Testing or measuring during manufacture or treatment; Reliability measurements, i.e. testing of parts without further processing to modify the parts as such; Structural arrangements therefor
    • H01L22/10Measuring as part of the manufacturing process
    • H01L22/12Measuring as part of the manufacturing process for structural parameters, e.g. thickness, line width, refractive index, temperature, warp, bond strength, defects, optical inspection, electrical measurement of structural dimensions, metallurgic measurement of diffusions
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01LSEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
    • H01L22/00Testing or measuring during manufacture or treatment; Reliability measurements, i.e. testing of parts without further processing to modify the parts as such; Structural arrangements therefor
    • H01L22/20Sequence of activities consisting of a plurality of measurements, corrections, marking or sorting steps
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/40Imaging
    • G01N2223/418Imaging electron microscope
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/611Specific applications or type of materials patterned objects; electronic devices
    • G01N2223/6116Specific applications or type of materials patterned objects; electronic devices semiconductor wafer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/646Specific applications or type of materials flaws, defects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/60Specific applications or type of materials
    • G01N2223/646Specific applications or type of materials flaws, defects
    • G01N2223/6462Specific applications or type of materials flaws, defects microdefects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/22Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
    • G01N23/225Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion
    • G01N23/2251Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • G06T2207/10061Microscopic image from scanning electron microscope
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20076Probabilistic image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30148Semiconductor; IC; Wafer

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computer Hardware Design (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Molecular Biology (AREA)
  • Evolutionary Computation (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
  • Testing Or Measuring Of Semiconductors Or The Like (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Telephonic Communication Services (AREA)
  • Train Traffic Observation, Control, And Security (AREA)
  • Maintenance And Management Of Digital Transmission (AREA)

Description

本開示は、欠陥の発生を推定するシステム、及びコンピューター可読媒体に係り、特に、確率的に非常に稀に発生する微細なパターン欠陥の発生確率を高精度に推定するシステム、及びコンピューター可読媒体に関する。
半導体デバイスの微細化に伴い、ナノメートルオーダーの幅を持つパターンを可視化できる走査電子顕微鏡(Scanning Electron Microscope:SEM)を用いた計測や検査の重要性が増大しつつある。非特許文献1には、確率的に非常に稀に発生する微細なパターン欠陥が説明されている。
Peter De Bisschop,"Stochastic effects in EUV lithography: random, local CD variability, and printing failures," Journal of Micro/Nanolithography, MEMS, and MOEMS,16(4),041013(2017)
確率的なパターン欠陥は、例えば膨大な数のパターンを露光したときに、ランダムに発生する欠陥である。その発生確率は極めて低いが、先端デバイス上には極めて多くの微細パターンが存在するため、発生確率が低い欠陥であっても、歩留まりに影響を与えることになる。また、確率的なパターン欠陥は極めて微細であるため、光学式検査装置での検出は困難であり、SEMのような微細な欠陥を大きく拡大して画像化する装置を用いて検査することが望まれる。しかしながら、高倍率である程、所定面積の検査を行うためには多くの時間を要することになる。
一方で、半導体デバイスの製造工程では、製造装置、材料、プロセス等の条件の変動(ウエハ異物、レジスト塗布現像面内分布変動、マスクや光学系の劣化)をいち早く捉え、これらの適切な調整や必要な処理を行うことが求められている。例えばレジストに光を投影してパターニングを行う投影露光装置のフォーカス条件や露光量等が変動すると、パターンの出来栄えが変わるため、試料を評価して、適切なタイミングで調整を行うことが求められるが、確率的なパターン欠陥の有無を検出するためには多くの時間を要し、一方で高頻度な製造装置の調整は半導体デバイスの生産量低下の要因となる。
非特許文献1には、確率的パターン欠陥についての紹介はあるが、確率的パターン欠陥のような極めて稀に発生する欠陥についての適切な評価法に関する開示がない。以下に、少ない検査点数で確率的パターン欠陥のような欠陥の発生を推定することを目的とするシステム、及びコンピューター可読媒体について説明する。
上記目的を達成するための一態様として、ウエハ上の複数の計測点の計測データに基づいて求められる、パターンのエッジが、第1の位置に存在する確率に関する第1のデータを、取得或いは発生、前記エッジが前記第1の位置にあるときに、当該第1の位置とは異なる第2の位置と前記第1の位置を含む領域を膜欠陥が覆う確率に関する第2のデータを取得或いは発生、前記第1のデータと前記第2のデータの積に基づいて、前記欠陥の発生確率を予測するシステム、及びコンピューター可読媒体を提案する。
上記構成によれば、少ない検査点数で確率的パターン欠陥のような欠陥の発生を推定することが可能となる。
確率欠陥(Stochastic defect)の発生確率予測システムの概要を示す図。 ウエハ上に形成されたパターンの近傍に生じる欠陥の種類を説明する図。 レジストパターンのエッジ存在確率、レジストパターン間の非パターン領域の一部に残存する膜欠陥の連続確率、及びパターン欠陥発生確率の関係を示す図。 連続膜欠陥確率P2と膜欠陥発生確率P3との関係を示す図。 加速条件下でのLCDU計測に基づいて、予め膜欠陥確率関数P2を求めておくことによって、確率欠陥の発生確率を予測するシステムの概要を示す図。 加速条件下での計測に基づいて、予め膜欠陥確率関数P2を求めておき、当該P2と実測により求められたエッジ残存確率P1から、欠陥確率P_defectを推定する工程を示す図。 欠陥発生確率推定システムの概要を示す図。 計測ツールの一種である走査電子顕微鏡の一例を示す図。 半導体製造工程と欠陥発生確率推定工程との関係を示すフローチャート。
昨今、波長13nmのEUV(Extreme Ultraviolet)光を用いた投影露光法により20nm以下の集積回路パターンをウエハ上に大量かつ高速に転写するEUVリソグラフィの開発が進んでいる。しかしながら、EUVリソグラフィを用いてパターンを微細化しようとしたとき、ストカスティック欠陥、又は確率的パターン欠陥の存在が歩留まり低下の要因になることが考えられる。すなわち、微細ラインアンドスペースパターンや密集ホール、ドット等、同一寸法・同一設計の膨大な数のパターンをEUV露光したとき、ライン間のブリッジ、ラインの断線、ホールやドットの消失等の致命的な欠陥がランダムな位置に発生する。
その発生確率は10-12~10-4程度と極めて低いが、先端デバイスでは300mm直径のウエハ上には、デバイスレイヤ毎に1012個以上の微細パターンが存在するので、仮に上記確率で欠陥が発生するとウエハ上に多数の欠陥が発生し歩留まりが低下してしまう。又、欠陥発生確率は、パターンの仕上がり寸法によって大きく変化する(例えば、1nmで1桁等)。このため、パターンの寸法管理の重要性が指摘されている。
確率的欠陥、又は確率的パターン欠陥の発生確率は極めて低いので、これを直接管理するには、膨大な数のパターンを検査する必要がある。一方、欠陥寸法は微細なので、光学的に検出することは難しく、電子線検査が望まれる。しかしながら、光学検査に比べて低速の電子線検査で、膨大な数のパターンを検査するには膨大な時間を要する。例えば、1ミクロン角の画像を時間当たり1万枚取得する場合、1mm角領域の全面を検査するには100時間を要する。画像サイズを10ミクロン角、時間当たり画像数を10万枚としても、1cm角より大きな集積回路チップやウエハ全面の検査を行うのは現実的に極めて困難であることがわかる。
一方、ウエハ上のある部分的な領域(例えば1mm角)に含まれる多数の同一設計のパターン(例えばメモリセルパターン等)の寸法を計測すると、その分布が正規分布から外れる。以下に、パターン寸法の分布と確率的パターン欠陥との関係を利用して、現実的に可能な数のパターン寸法計測結果から、確率的パターン欠陥の発生確率を高い精度で予測する例について説明する。
まず、EUV光等によるレジストの反応過程について説明する。レジスト膜にパターン状の光を照射すると、フッ化アルゴン(ArF)エキシマレーザー(波長193nm)等のDUV(Deep Ultraviolet)光ではレジスト膜中の分子結合が光子を吸収して光反応を生じる。又、EUV光では膜中の原子が光子を吸収して、等方的に光電子を放出する。光電子は周辺原子による弾性散乱により進行方向を変え、非弾性散乱によりエネルギを失うとともに2次電子を発生する。この現象は、光電子はエネルギが一定レベル以下となった時点で停止する。エネルギが5~10meVの2次電子により、分子が反応を生じる。当該反応はレジスト材料の反応原理により様々である。
現在代表的なレジスト系は化学増幅レジスト(CAR)で、酸発生剤(PAG)と複数個の極性変換基を含むマトリクスポリマーを構成要素とする。DUVでは光子、EUVでは2次電子のエネルギを吸収した酸発生剤(PAG)分子が酸を生成し、酸発生位置を中心とする酸拡散長の範囲内で酸が触媒として作用して、触媒連鎖長個(1個の酸が誘起する触媒反応数)のポリマー・分子の極性基の極性反転反応を生ずる。ポリマー内の極性反転基数がある閾値を超えるとポリマーの溶解性が反転して、現像プロセスによりパターンが形成される。
なお、40nm以下のパターニングにおいて、パターンのエッジに、図2(a)に示すようなナノメータレベルのランダムな凹凸(Line Edge Roughness:LER)が形成されることが確認されている。その主な原因は、上記光子数のショットノイズとレジスト材料の空間的な組成ばらつき、反応の確率的変動等である。
LERはフッ化アルゴンエキシマレーザーを用いた露光によって形成されるパターンと、EUV光を用いた露光によって形成されるパターンの双方に形成される一方、確率的パターン欠陥はフッ化アルゴンエキシマレーザーを用いた露光によって形成されるパターンでは確認されなかった。従って、後者は、20nm以下の領域に特有のレジスト材料の問題か、又は、波長の違い(光子数及び光子エネルギの違いによる反応過程の違い)によるものと考えられるが、その原因は必ずしも明らかでない。
ここで、主パターン形成と確率的パターン欠陥発生を以下のように定義する。設計上のレジストが残るべき領域(パターン部)に所定の厚さの膜が残り、かつ、レジストが除去されるべき領域(非パターン部)のレジストが完全に除去された状態をパターン形成と定義する。次に、上記非パターン部の一部で所期膜厚の所定の割合(例えば10~30%)以上のレジストが残存する場合を残存膜欠陥、上記パターン部の一部でレジスト膜が所期膜厚の所定の割合(例えば10~30%)以上消失する場合を消失膜欠陥、と定義する。
上記所定の割合はエッチングプロセスによって決まる。さらに、上記膜欠陥の領域が設計パターンの機能を阻害する範囲に広がり、引き続き行われるエッチングにおけるレジストマスク機能に支障をきたす場合を、パターン欠陥と定義する。例えば、残存膜欠陥が2つの独立した設計パターン間を結ぶ領域に広がった場合(ブリッジ欠陥、例えば図2(b))、微細開口パターンの全面を埋めてしまった場合(消失開口、図2(d)の中央のパターン)、膜消失欠陥が、線状の設計パターンを切断する場合(ブリーク欠陥、図2(c))、微細ドットパターンの一部が失われた場合(消失ドット、図2(e))、等がある。
なお、確率的欠陥の発生の推測に基づいて、例えば、露光装置の露光条件へのフィードバックの適切なタイミングを特定でき、結果として高頻度な露光条件調整を行うことなく、半導体製造を継続することができるが、これに限ることはなく、例えばマスク汚染を起因とする確率的欠陥の発生を推定することも可能である。マスク汚染の確率的欠陥の発生の推定によって、マスククリーニングを適切なタイミングで実施することが可能となる。更に、ウエハ裏面に付着する異物を原因とするフォーカスシフトに起因する確率的欠陥の推定に基づいて、ウエハ裏面のクリーニング処理のタイミングを適切に計画することもできる。また、本開示の推定法は、レジストプロセス(塗布現像トラック)の異常によるレジスト感度シフト等に起因する確率的欠陥の発生推定にも用いることが可能である。
更に、本開示の手法によれば、異常なウエハを検出することもでき、異常なウエハの排除にも適用することが可能となる。
発明者の検討によれば、上記ランダムに発生するパターン欠陥の発生確率は、上記主パターンのエッジ位置の存在確率と、上記エッジに接する領域に連続的に膜欠陥が発生する確率の2つの確率の積で表わされる。即ち、2つの事象の条件付き確率で表すことができる。そこで、以下に説明する実施例では、図1に示すように、SEM等の計測手段により上記パターンを計測して上記エッジ位置の存在確率を求め、物理化学シミュレーション又は経験的・統計モデルにより上記連続膜欠陥発生確率を求め、両者の積和によりストカスティックなパターン欠陥(以下パターン欠陥又は単に欠陥と記す)の発生確率を推定する。
ここで発明者の検討によれば、上記連続膜欠陥発生確率はエッジ位置の関数となる。従って、プロセス条件、ウエハ状態、マスク状態等が変化したとき、エッジ位置の存在確率分布が変化するとともに、エッジ位置の変化に応じて膜欠陥発生確率が変化して、欠陥発生確率が変化する。さらに、上記各種条件・状態変化に伴う膜欠陥発生確率のエッジ位置依存性の変化は小さいので、エッジ位置の関数としての確率関数の形を一定と近似する。
これにより、10の-12乗の確率で生じる欠陥を検出するために必要だった10の12乗個のパターン検査を、例えば10の-6乗の確率で生じるエッジ位置分布の異常を検知するために必要な6乗個のパターン計測で代替することができる。これにより現実的な時間内で上記欠陥発生の危険性を検知することが可能となる。
以下、簡単のため、x方向に強度変調した周期パターン(y方向に長手方向を持つラインアンドスペースパターン等)において、レジストパターン間の非パターン領域の一部を残存膜欠陥が覆う場合(図2(b)に相当する)を例に、図3を用いて説明する。発明者の検討によれば、この場合の確率的パターン欠陥の発生確率は、近似的に以下のように求めることができる。
まず、主パターンのエッジが位置x_edge(第1の位置)に存在する確率をP1(x_edge)とする。P1(x_edge)は、基準エッジ位置(図3(a)の設計エッジ位置)に対する相対的なエッジ位置座標をx_edgeとしたときの、複数の計測点の計測データから求められるx_edgeの頻度分布を示している。更に、エッジがx_edgeにあるときに、当該x_edgeと、欠陥評価点x(第2の位置)との間に連続して膜欠陥を生じる確率関数をP2(x、x_edge)とする。このように求められたP1(x_edge)とP2(x、x_edge)より、xにおけるパターン欠陥確率P_defect(x)を数1によって求める。
Figure 0007080992000001
即ち、パターン内のある位置とエッジを結ぶパターン欠陥発生確率は、エッジ出現確率P1(x_edge)と、エッジと上記位置を結ぶ領域を膜欠陥が覆う確率P2(x,x_edge)の2つの部分確率の積で表される。
次に、実際の製造工程における欠陥発生とその検出について説明する。まず前提として、露光プロセス、レジスト材料、露光マスク、ウエハの基準となる状態(以下プロセス状態と呼ぶ)において、パターン欠陥発生確率が所望のレベルに抑えられているとする。そこに、露光条件変動、マスク汚染、ウエハ状態変化(裏面異物によるフォーカス位置変動等)等によりプロセス状態が変化したとする。
このとき、上式の2つの要素、エッジ出現確率P1(x_edge)、及びエッジがx_edgeにあるときある位置x(例えば非パターン領域の中心等の欠陥基準点)を結ぶ領域を膜欠陥が覆う確率P2(x,x_edge)、各々の変化について説明する。
まず、実際のウエハ上で同一設計パターンに対して形成されるレジストパターンのエッジには前述のLERが存在し、図3(a)に示すようにエッジに沿ってエッジ位置x_edgeが変動する。そこで、基準プロセス条件において、例えば1本のラインパターンのエッジに沿って多数のエッジ点を計測することにより、P1(x_edge)を得る(図3(b))。P1(x_edge)は、上記プロセス状態変化に応じて、図3(b)の条件AからBのように変化する。
なお、上記分布変化の割合は平均値から離れた部分で大きくなる。例えば、正規分布やポアソン分布ではその一般的性質として、平均値又は標準偏差が変化したときの値の変化率は、平均値から離れるほど増大する。これは、必ずしも正規分布やポアソン分布では表されない実際の寸法(エッジ位置)分布においても観察される。
次に、エッジと欠陥評価点を結ぶ領域を膜欠陥が覆う確率P2(x,x_edge)の変化について述べる。P2(x,x_edge)は、いくつかの仮説モデルにより計算できるが、例えば、数2のように単位面積あたり膜欠陥が生じる確率P3(x´)の上記領域内における直積として表すことができる。
Figure 0007080992000002
後述するシミュレーションにより求めたP3を図4(a)に示す。露光条件が変化したときのP3自体の変化は、図中実線及び点線に示すように比較的小さい。単位面積あたりの膜欠陥発生確率は上記状態変化(例えば高々数%程度の露光量プロファイルの変化)にほぼ比例するからである。
一方、図4(a)に示したP3に対して上式で計算したP2のx、x_edge依存性を図4(b)及び図(c)に示す。図4(c)には、非パターン領域の中心(x=16nm)におけるP2のx―x_edge依存性を示す。P2の値は、|x―x_edge|に対して指数関数的に減少する。その他のモデルにおいても、およそ|x―x_edge|に対して指数関数的に減少する。このように、条件変化によるP2の変化は、x_edgeの変化を介する部分が支配的である。そこで、本例では基本的に、条件が変化したときのP2の変化をx_edgeの変化のみに依存すると近似する。
図3(b)(c)の結果より、欠陥確率は数1を用いて図3(d)のように求められる。わずかな(2nm程度の)平均エッジ位置の変化により、非パターン領域部の欠陥確率が2桁程度変化していることがわかる。
以上より、プロセス状態が変化すると、エッジ位置x_edgeの分布P1(x_edge)が変化するとともに、エッジと欠陥基準点を結ぶ領域を膜欠陥が覆う確率P2(x,x_edge)が指数関数的に変化する。本例では、パターン欠陥発生確率変化に及ぼす要因は、エッジ変動が支配的であると近似して、エッジ位置の分布については、実際のウエハ上で観察される計測値を用い、直接計測することのできない膜欠陥発生確率空間分布の積和部については、統計解析による推定値を用いて、演算式に基づいて欠陥発生確率を推定する。又は欠陥発生の危険性を予知する。
先に述べたとおり、欠陥発生確率は極めて小さい(例えば10-12)ので、実際に欠陥を検出するためには膨大な数のサンプル計測が必要で、現実的には困難であった。
以下に説明する実施例では、これを、各々エッジ位置のみの関数として表される2つのより大きな確率(例えば各々10-6)の積として表すことにより、1012点より遥かに少ない例えば10点程度の計測により求めることができる。100万点のエッジ位置の計測結果は、確立された計測ツールの1種であるCD-SEM(Critical Dimension-SEM)によって、例えば数百から数千枚の画像を撮像することにより得ることができ、その計測に要する時間も数分から数十分程度である。
なお、SEMによるエッジ位置の計測は、1本の走査線上でエッジ位置を検出してもよいが、実際には近接する複数の走査線に対して得られた信号の平均、又は複数の走査線の各々に対して得られたエッジ位置の平均を用いてもよい。但し、得られたエッジ位置が局所的なエッジ位置の変動を反映していることが望ましく、上記平均に用いるエッジ方向の長さは、最小パターン幅と同程度かそれ以下(例えば典型的には10~20nm)であることが望ましい。
なお、エッジ位置が平均値から離れたところ(例えばスペース幅が狭まったところ)では、状態変化によるエッジ出現確率の変化が大きく、かつ、|x―x_edge|(例えばエッジとスペース中心間の距離)が小さくP2の値が大きいので、欠陥発生確率に及ぼす影響が大きい。従って、エッジ位置分布のテール部の計測が特に重要である。
図8は、主にP1を求めるためのエッジ位置を計測し、計測データ(メトロロジーデータ)を出力するための計測ツールの1種であるCD-SEMの一例を示す図である。CD-SEMは走査電子顕微鏡と、検出器の出力に基づいて、パターンの寸法やエッジ位置(例えば設計エッジ位置との差異)を特定するための1以上のコンピューターシステムより構成されており、当該コンピューターシステムは、所定のコンピューター読み取り可読媒体に記憶されたプログラムを読み出し、後述するような欠陥確率推定処理を実行するように構成されている。コンピューターシステムは、走査電子顕微鏡と通信可能に構成されている。コンピューターシステムは、1以上の伝送媒体で走査電子顕微鏡に接続され、走査電子顕微鏡と離間した位置に設置されても良いし、走査電子顕微鏡のモジュールとするようにしても良い。
図8に例示する走査電子顕微鏡は、電子源801から引出電極802によって電子ビーム803が引き出され、不図示の加速電極により加速される。加速された電子ビーム803は、集束レンズの一形態であるコンデンサレンズ804により絞られた後、走査偏向器805により偏向される。これにより、電子ビーム803は、試料809上を一次元的又は二次元的に走査する。試料809に入射する電子ビーム803は、試料台808に内蔵された電極に、負電圧を印加することによって形成される減速電界により減速されると共に、対物レンズ806のレンズ作用により集束されて試料809の表面を照射される。試料室807内部は真空が保たれている。
試料809上の照射箇所からは電子810(二次電子、後方散乱電子等)が放出される。放出された電子810は、試料台808に内蔵された前記電極に印加された負電圧に基づく加速作用により、電子源801の方向に加速される。加速された電子810は変換電極812に衝突し、二次電子811を発生させる。変換電極812から放出された二次電子811は、検出器813により捕捉され、捕捉された二次電子量により検出器813の出力Iが変化する。この出力Iの変化に応じ、表示装置の輝度が変化する。例えば二次元像を形成する場合には、走査偏向器805への偏向信号と、検出器813の出力Iとを同期させ、走査領域の画像を形成する。
なお、図8に例示するSEMは、試料809から放出された電子810を変換電極812において二次電子811に一端変換して検出する例を示しているが、無論このような構成に限られることはなく、例えば加速された電子の軌道上に、電子倍像管や検出器の検出面を配置する構成を採用しても良い。制御装置814は、撮像レシピと呼ばれるSEMを制御するための動作プログラムに従って、上記SEMの各光学要素に必要な制御信号を供給する。
次に検出器813で検出された信号はA/D変換器815によってデジタル信号に変換され、画像処理部816に送られる。画像処理部816は複数の走査によって得られた信号をフレーム単位で積算することによって積算画像を生成する。
ここで、走査領域の1回の走査で得られる画像を1フレームの画像と呼ぶ。例えば、8フレームの画像を積算する場合、8回の2次元走査によって得られた信号を画素単位で加算平均処理を行うことによって、積算画像を生成する。同一走査領域を複数回走査して、走査毎に1フレームの画像を複数個生成して保存することもできる。
更に画像処理部816は、デジタル画像を一時記憶するための画像記憶媒体である画像メモリ818と、画像メモリ818に記憶された画像から特徴量(ラインやホールの幅の寸法値、ラフネス指標値、パターン形状を示す指標値、パターンの面積値、エッジ位置となる画素位置等)の算出を行うCPU817を有する。
さらにまた、各パターンの計測値や各画素の輝度値等を保存する記憶媒体819を有する。全体制御はワークステーション820によって行われる、必要な装置の操作、検出結果の確認等がグラフィカルユーザーインタフェース(以下、GUIと表記する)によって実現できるようになっている。また、画像メモリは、走査偏向器805に供給される走査信号に同期して、検出器の出力信号(試料から放出される電子量に比例する信号)を、対応するメモリ上のアドレス(x,y)に記憶するように構成されている。また、画像処理部816は、メモリに記憶された輝度値からラインプロファイルを生成し、閾値法等を用いてエッジ位置を特定し、エッジ間の寸法を測定する演算処理装置としても機能する。
以下、欠陥発生確率推定方法、欠陥発生確率推定を実行するシステム、及び欠陥発生確率推定をコンピューターシステムに実行させるプログラムが記憶されたコンピューター可読媒体について説明する。
以下、欠陥発生確率推定の一実施例について説明する。まず、CDSEMを用いて微細パターン付ウエハ上の100μm角程度の領域に一様に存在する周期パターンの局所寸法を多数点で計測して、各パターン中心に対するパターンエッジ位置x1の分布F(x1)を求める。ここでは、EUV露光装置と化学増幅系レジストを用いてウエハ上に形成した16nmラインアンドスペース(L/S)パターンに対して、SEMの視野(FOV)を1ミクロン角として、約250枚のSEM画像を取得する。各画像に含まれる約25本のラインの各々に対して、その中心から右側エッジまでの距離をエッジ方向5nm毎に5nm範囲で平均して計測する。これにより画像あたり約4000点のエッジ位置、合計100万点のエッジ位置を計測し、エッジ位置の分布、P1(x_edge)を求めた。
次に、計算機を用いて欠陥発生の部分確率を求める。エッジがx_edgeにあるときに、x_edgeとxを結ぶ領域に膜欠陥を生じる確率関数P2(x,x_edge)を、以下の2つの方法により求めた。
まず、第1の方法として、設計パターンのレイアウト情報及び露光条件から、設計パターン内位置xにおける(単位面積当たり)膜欠陥発生確率P3(x)を以下のようにして求める。モンテカルロ法を用いた物理化学シミュレーションにより、位置xにおけるレジスト材料ポリマー極性基の極性反転反応やポリマーの溶解性反転等のイベント密度n(局所体積当りイベント数)を計算し、位置xとこれらのイベント数密度nを確率変数とする第1の確率密度分布関数g1(x、n)を算出する。モンテカルロ法を用いたg1の算出法については、別の実施例において説明する。
g1は、所定の物理法則に基づくシミュレーションにより求めるが、x、nを関数とする平均・標準偏差を持つ正規分布、ポアソン分布等を用いてもよい。これらの分布関数の母数としては、上記物理化学シミュレーションの結果をフィッティングして求めることができる。
次にg1を所定のnの範囲で積分してxにおいて膜欠陥が発生する単位面積当たり確率P3(x)を求める。さらに、x=x_edge~xの間でP3(x)の直積を求め、これをP2(x、x_edge)する。この方法で求めるのは、ショットノイズが偶発的に隣接接続することにより生じる欠陥の確率である。
第2の方法では数3を用いて、P2を算出する。まず、上記パターンに対して光学像シミュレーションにより光強度分布I(x)を計算する。次に、座標r0で吸収された光子が、エッジ位置x_edgeと位置xを結ぶ領域のポリマーの溶解性を反転させる確率、P4(x,x_edge,r0)を求める。P4は、上記光子吸収により生じた光電子の散乱による軌跡に沿って、上記領域で所定の密度であり、2次電子が発生する確率で近似してもよい。これらの確率は、所定の物理法則に基づくモンテカルロシミュレーションにより求めるが、x,x_edge,nを確率変数とする様々な確率分布関数で近似してもよい。モンテカルロ法を用いた算出法については、別の実施例において説明する。エッジがx_edgeにあるときに、位置xとの間に連続膜欠陥を生じる確率関数P2を、数3を用いて求める。conv(A,B)はAとBの畳み込みである。
Figure 0007080992000003
P4に基づくP2の導出は、例えばウエハ上に形成されたレジスト膜に投影された光の光子により発生した光電子が、レジスト膜の座標x_edgeから座標xの間に近接して2次電子を発生する確率を求めることと、定義することができる。
第1の方法と第2の方法で求められたP2は、異なる要因に基づく膜欠陥を生じる確率関数であるため、線形結合等により、両者を足し合わせる。但し、一方の影響が支配的であり、他方の影響を無視できるのであれば、一方の方法のみからP2を算出するようにしても良い。
以上のようにして求められたP1(x_edge)、P2(x,x_edge)、及び数1を用いてxにおける欠陥確率P_defect(x)を求める。
上述の処理を、ウエハ内30チップに対してチップ内の1領域、及びウエハ内3チップに対してチップ内の10領域、合計60領域に対して行い、欠陥発生確率のウエハ面内、及びチップ内分布を求めた。ウエハエッジ部とチップの特定位置においてエッジ位置分布のテールに広がりが観察され、これに応じて欠陥発生確率が10-11から10-9に2桁劣化することが予測された。ウエハ上及びチップ上の当該位置の1mm角領域に対して高速電子線検査装置を用いて欠陥検査を行ったところ、ほぼ上で予測された密度の欠陥が検出された。
なお、本実施例では、ラインアンドスペースのような特定の一次元方向に長いパターンが周期的に配列された周期パターンを例にとって説明したが、孤立パターン、又はホールやドット等の2次元パターンの場合も同様にして適用可能である。
次に、上記物理化学シミュレーションを用いずに、計測結果の統計処理により数1のP2(x、x_edge)を求める方法について説明する。
露光量やマスク寸法等の露光・マスク条件を変えると、レジストパターンの平均仕上がり寸法が変化する。このときパターン欠陥発生確率を計測すると、欠陥確率は上記平均仕上がり寸法に対して極めて敏感に(指数関数的に)変化する。先に説明したように、本開示では、これを「局所的な」エッジ位置の発生確率と、それに伴う局所的な連続膜欠陥生成確率の変化の積として捉える。そこで、LCDU又はLER計測により求めた局所的なエッジ位置の発生確率P1(x_edge)と、同じサンプルに対して実際に計測した欠陥発生確率P_defectに対して、数4を満たすようなP2´(x_edge)を、例えば統計的フィッティングにより求める。
Figure 0007080992000004
特に、本実施例では、P2´(x_edge)を欠陥発生確率が高い加速条件下で求め、次にこれを通常のウエハ処理条件に適用する(図5)。ここでLCDU(Local Critical Dimension Uniformity)とは、局所的なフィールドにおけるCD値の均一性を示す指標値であり、局所的なフィールドの複数のパターンを計測した結果得られる値である。
まず、欠陥発生確率が高く、現実的に可能な計測点数内で0でない欠陥が現れるような複数(i=1~m)の露光・マスク条件下で、エッジ位置x_edgeの頻度分布(又は寸法分布)P1(x_edge,i)と、欠陥発生確率Pdi_observed(i)を計測する(図6(a)、(b)のAとBに相当)。次に、数5のδが最小となるように、P2´(x_edge)を決定する(図6(c))。
Figure 0007080992000005
次に、通常の露光・マスク条件下でエッジ位置頻度分布P1(x_edge)を計測し(図6(d)のC、D)、これに対して、上記加速条件下で求めたP2´(P2)を用いて、数1式より欠陥発生確率を推定する(図6(e)のC、D)。
なお、本開示によって求められる欠陥発生確率の計算結果は、寸法またはエッジ位置の頻度分布のテールが支配していると考えられる。
通常の露光条件では、現実的に可能な時間内で、(式1)の被積分関数のピークを与えるxにエッジ分布のテールが伸びるまで計測を続けることが困難な場合がある。この場合、分布のテールを適宜外挿する等してもよい。特に、P1(x_edge)は、分布の中心から離れるにつれ正規分布からずれるため、正規分布からずれ始めた時点で、図6(d)の曲線C、Dのように、正規分布に接してのびる指数分布等によりテールを外挿近似(曲線C、Dの点線部分)する等してもよい。
また、P2´(x_edge)は適当な関数型を仮定することが望ましい。その原型としては実施例1又は2で求めた物理化学シミュレーションに基づく方法や、観測された欠陥確率寸法依存性を寸法分布でデコンボリューションする等の方法が適用できる。
さらに、本実施例は実施例1で説明した方法と組み合わせることもできる。例えば、本実施例を、ある設計マスクパターンに対して適用して、P2´(x_edge)を得たとする。次に、マスクパターンの設計を変更した場合、変更後のマスクに対するP2´(x_edge)は、改めて求めなおす必要がある。
この場合、改めて本実施例を適用する代わりに、実施例1で示した方法により、変更前後の連続膜欠陥確率P2_old(x_edge)、P2_new(x_edge)を計算により求め、変更後のマスクに対するP2´を、数6を用いて求める。
Figure 0007080992000006
上記変更は、露光条件やその他プロセス条件を変えた場合も同様である。この場合、実施例1を単独で用いる場合と比べて、シミュレーションを厳密にキャリブレーションしないでも、様々な条件に対応できる。
なお、寸法ばらつき(LCDU又はLER)があると、パターン欠陥確率はその平均パターン寸法依存性と寸法分布のコンボリューションとなり、劣化(増大)する。しかしながら、上記計測された欠陥確率は、既にLCDU、LERにより局所的にエッジ位置が変動するパターン群の平均寸法に対して得られた結果である。即ち、平均寸法がある値のときの欠陥発生確率を論じており、寸法ばらつき分布との関係は明らかでない。本開示によれば、観察された寸法ばらつき(LCDU又はLER)とその変化から、欠陥発生確率とその変化を定量的に推定することができる。
図9は、具体的な欠陥発生確率P_defectの算出工程の一例を示すフローチャートである。図9の例では、投影露光装置の条件出しを行うときに用いられるFEM(Focus Exposure Matrix)ウエハ等から、フォーカス、或いは露光量が適正ではないが故に欠陥が頻発する露光条件(加速条件)で露光された複数のパターン(露光条件)を選択し、当該選択された露光条件のパターンの計測に基づいて、P2´(x_edge)を導出するためのP1(x_edge)とP_defectを導出する例について説明する。
まず、投影露光装置の適正な露光条件を特定するために、露光条件(フォーカス、露光量など)を変化させ、それぞれの露光条件にて露光を行う(ステップ901)。次に各露光条件で形成されたパターンについてCD-SEM等を用いて多点計測を行う(ステップ902)。多点計測の結果、半導体デバイスの設計データ、或いは設計データに基づいて定められる許容値に合致する計測値のパターンを形成する露光条件を選択することによって、量産時の露光条件を決定する(ステップ903)。
一方、露光条件を決める条件出し工程では、露光条件を変化させて露光を行い、様々な露光条件でパターンを形成するため、その中には、パターンが適正に形成されない(欠陥が含まれる)露光条件が含まれる。よって、その中から欠陥発生率の高い露光条件(加速条件)のパターンを選択し、露光条件ごとにP1とP_defectを導出(ステップ904、905)し、上述のような統計フィッティングにより、P2´を導出する(ステップ906)。例えば、10万点の計測点の内、欠陥数が100点のウエハ(或いはチップ)、及び欠陥数が50点のウエハのような欠陥確率が異なる複数の試料ごとに、P1とP_defectを導出し、これらの試料に共通するP2´を導出することで、後述する欠陥発生確率を導出に要するデータを予め用意する。
実際の半導体量産工程では、露光装置の条件出し工程で選択された露光条件を用いたパターニング(ステップ907)が行われ、プロセス管理のために、CD-SEM等を用いたパターン計測(ステップ908)、及び必要に応じた露光装置の調整が行われる。本開示は、確率的欠陥の発生を限られた計測点に対する計測に基づいて、推定することによって、計測に要する手間や時間を抑制するためのものであり、本開示の手法によれば、高頻度な露光装置の調整による生産性の低下を抑制と、高い歩留まりの維持の両立が可能となる。
図9に例示するように、量産工程で形成されたウエハに対する多点計測(例えば10万点)結果に基づいてP1(x_edge)を導出する(ステップ909)。次に、ステップ906で予め求められているP2´(x_edge)と、ステップ909で導出されたP1(x_edge)から、P_defectを導出する(ステップ910)ことによって、量産時の欠陥発生確率を導出する。
導出された欠陥発生確率と、予め設定されている所定条件(例えば1012の計測点中、n点「以下」の欠陥)を比較し、所定条件を満たす場合には、そのまま製造を継続する(ステップ913)。一方、所定条件を満たさない場合(n点以上の欠陥が発生する場合)は、原因の究明を行い、しかるべき対策を行う。例えば、露光装置の設定により対処可能な場合は、調整(ステップ912)を行い、製造を継続する(ステップ913)。また、原因が不明の場合には、当該ウエハの製造を中断する。
コンピューターシステムは、欠陥発生確率をCD-SEMに設けられた表示装置等に適時表示することで、オペレータは露光装置の適切な調整タイミングを把握することができる。また、欠陥発生確率が所定の条件を満たさなくなったときに、警報を発生することによって、オペレータは露光装置の適切な調整タイミングを知ることが可能となる。
なお、本実施例は教師付き機械学習により実装してもよい。即ち、例えばディープ・ニューラルネットワークやコンボリューショナル・ニューラルネットワーク等のニューラルネットワークを用いた機械学習法により、加速条件(i=1~m)を教師として、その寸法頻度分布関数P1(x_edge,i)を入力ベクトルとして、対応する条件に対する欠陥発生確率Pdi_observed(i)を出力とするニューラルネットワークを作成する。換言すれば、前記P1(x_edge)を入力、前記欠陥の発生確率を出力とする教師データを用いて学習されたパラメータを前記第2のデータとして中間層に備えた学習モデルを生成する。
しかる後に、上記ニューラルネットワークに、任意の寸法頻度分布関数P1(x_edge,i)を入力して、欠陥発生確率を予測する。即ち、構築された学習モデルに、P1(x_edge)を入力することによって、前記欠陥の発生確率を出力する。これにより、数5によることなく、連続膜欠陥確率P2(x_edge)は、自動的に上記ネットワーク中に埋め込まれる。予測結果は全数検査により適宜検証して、検証結果を用いてネットワークを補正する。
次に、欠陥発生確率推定システムの概要について説明する。図7に例示する欠陥発生確率推定システムは1以上のコンピューターシステムによって構成される。図7に例示するシステムには、CDSEM、電子線欠陥検査装置、計測結果データベース部、物理化学シミュレーター部、膜欠陥確率関数算出部、確率関数データベース部、パターン欠陥確率予測部、等のサブシステムが含まれている。
に例示したCDSEMは、測定対象パターン情報を入力して、多数パターンを計測した画像又は寸法データより、P1(x_edge)を求め、計測結果データベース部に保存する。なお、P1等はCDSEMに設けられたプロセッサによって演算し、その演算結果を、CDSEMに設けられたプロセッサで取得し、数1に例示したような演算を当該プロセッサで行うようにしても良いし、CDSEMとは別に設けられたプロセッサにて、CDSEMから多点計測結果(複数計測点のエッジ位置情報)を受け取り、当該プロセッサでP1等のデータ(第1のデータ)を演算により発生させるようにしても良い。
必要に応じて、物理化学シミュレーター部は、実施例1で述べた上記g1(x、n)、P3(x)、I(x)、P4(x、x_edge、0)、P2(x、x_edge)等を算出する。また必要に応じて、電子線欠陥検査装置は、対象パターン群のパターン欠陥数(欠陥確率)を計測する。連続膜欠陥確率関数算出部は、前記実施例2に記載の方法により上記物理化学シミュレーター部の計算結果を用いてP2(x、x_edge)を算出する。又は、連続膜欠陥確率関数算出部は、前記実施例2に記載の方法により上記P1(x_edge)と電子線欠陥検査装置により計測されたパターン欠陥数(欠陥確率)からP2(x、x_edge)を算出する。
確率関数データベース部は算出されたP2(x、x_edge)を保存する。パターン欠陥確率予測部は、上記P1(x_edge)とP2(x、x_edge)から欠陥確率P_defect(x)を算出する。以下、各サブシステムの詳細について説明する。なお、図7に記載の他のサブシステムについては後述する。
CDSEMで計測した画像、エッジ座標、寸法等のデータは、例えば計測結果データベースに保存される。物理化学シミュレーター部では、前記、g1(x、n)、P3(x)、I(x)、P4(x、x_edge、0)、P2(x、x_edge)等を算出する。以下算出方法について説明する。
例えば、前記実施例の第1の方法で用いるP2(x、x_edge)を、次のようにして計算する。まず、予測対象パターンの図形情報と露光条件・材料パラメータ等のプロセス情報を入力する。レジスト膜中の光学像強度分布を求め、光学像強度に比例した確率で光子の吸収イベントを生成させる。像強度分布は、平面方向はHopkins theory、膜深さ方向はLambert Beerの法則に従い近似計算したが、より厳密な方法を用いてもよい。EUV露光では、光子を吸収した原子が放出した光電子の散乱軌跡に沿って複数の2次電子が生成する。この過程を電子散乱のモンテカルロ法で計算した。2次電子発生位置で化学増幅系レジストの酸が発生したと近似して、酸の存在確率分布を、上記酸発生位置を中心として酸拡散長をぼけ量とするガウス分布と仮定する。酸の存在確率に比例する確率でレジストポリマーに含まれる複数極性基の極性反転反応イベントを生成させる。但し1個の酸あたりの極性反転数、単位体積あたりの極性反転数には上限を設定する。
次に、上記計算空間(レジスト膜)を、例えば1nm刻みの3次元立体格子で分割し、これによってできるボクセル(格子によってできる単位立方体)により、溶解性反転の最小単位であるレジストのマトリクスポリマーや分子を近似する。各ボクセルに含まれる各種イベント数をカウントし、ボクセルあたりの極性反転反応イベント数が一定の閾値を超えた場合、当該ボクセルの溶解性が反転すると判定する。さらに、平面方向の固定位置で、深さ方向に溶解性反転したボクセル数(ボクセルの溶解性反転イベント数)をカウントして、これをある閾値と比較することにより、局所的な(単位面積当たりの)膜の溶解性を判定する。例えば、ネガ型レジストの場合、深さ方向溶解性反転ボクセル数が第1の閾値を超えたとき膜が不溶化して局所的に残存膜欠陥が生じる、第2の閾値(>第1の閾値)を越えたとき膜全体が不溶化するとする。
上記レジスト空間をある一定の露光強度で露光したときの各ボクセルに含まれるイベント数は、モンテカルロシミュレーションの結果として統計的にばらつくので、各イベント数nを確率変数とする確率関数が求められる。1次元パターンの光学像においてレジスト膜内の光の減衰を無視すると、露光強度は座標xのみに依存するので、上記確率関数はxの関数となる。このようにして、深さ方向溶解性反転ボクセル数に対する確率関数g1に対して、例えばネガ型レジストにおける残存膜欠陥発生確率を数7により求めた。
Figure 0007080992000007
更に数2により、P2(x、x_edge)を計算した。同様にしてポジ型レジストにおける消失膜欠陥発生確率も同様の式により求められる。又、ネガ型レジストにおける消失膜欠陥発生確率、又はポジ型レジストにおける残存膜欠陥発生確率は、積分範囲をn<nc2に適宜変更する等により求められる。
また、実施例1の第2の方法で用いるP2(x、x_edge)については、次のようにして求める。1個の光電子の散乱軌跡上で隣接する2次電子発生は独立でない(相関がある)ので、連鎖の生じる確率自体を評価する。1個の光電子の散乱中に、至近距離(例えば1nm以下)で連続して発生した2次電子を2次電子鎖と定義する。
まず、モンテカルロ法を用いて多数の光電子の散乱軌跡に沿った2次電子発生イベントを発生させ、2次電子鎖を抽出する。2次電子鎖に対して光子吸収位置から最も遠い点をその終点、光子吸収位置から終点方向へ投影したときの電子鎖の長さを長さと定義し、終点までの距離(range)と長さ(length)を確率変数とする確立密度分布、P_SEstring (length, range)を求める。次に、これを用いて、r0=(x0、y0)で発生した光電子が、x_edgeにあるエッジからxまで延びる2次電子鎖を形成する確率密度分布P4(x、x_edge、r0)を計算する。
エッジがx_edgeにあるとき、x_edgeとxを結ぶ領域を欠陥が覆う確率P2(x、x_edge)は、発生する全ての光電子についてP4(x、x_edge、r0)足し合わせて求める。即ち、光子吸収イベント数の分布I(r0)とP4(x、x_edge、r0)を数8のように畳み込み積分することによって求められる。
Figure 0007080992000008
上記P2(x、x_edge)の算出は、連続膜欠陥確率関数算出部で行うとともに、物理化学シミュレーター部で算出したg1(x、n)、P3(x)、I(x)、P4(x、x_edge、r0)、P2(x、x_edge)等は、確率関数データベース部に保存する。又、連続膜欠陥確率関数算出部は、前記実施例2で説明した統計フィッティングによるP2(x、x_edge)の算出も行う。
欠陥確率予測部は、上記P1(x_edge)、P2(x、x_edge)を入力、パターン欠陥発生確率P_defect(x)を算出する。あわせて、これらの確率関数の表示部、欠陥数、欠陥確率検証結果等の入力を含む外部入出力部、これらを制御するためのGUIを有する。物理化学シミュレーターで行う計算の一部は予測システムで行ってもよい。g1(x、n)、P3(x)、I(x)、P4(x、x_edge、r0)、P2(x、x_edge)等は、予測システム内部でデータベースとして保存して、必要に応じて利用してもよい。
ここでは、システムの運用に関する実施例について述べる。実施例1に述べたように、本開示は、ウエハ又はチップ上であらかじめ決められた複数の部分領域に対して適用することにより、上記各領域のパターン欠陥発生確率の分布を予測する方法やシステムに関するものである。
一方、ウエハ又はチップ上の1部分領域の欠陥発生確率予測に必要なパターン寸法の計測点数は大きいので、ウエハ上多数の部分領域に対して欠陥発生確率を予測するには長時間を要する。そこで、ウエハ上の欠陥発生危険度の空間分布を別手段により予測し、危険度が高い領域に限定して上述のような欠陥発生確率を求めるようにしてもよい。上記別手段としては、露光装置等の設定露光・フォーカス分布、ウエハ形状分布計測結果、光学的寸法計測結果、光学的欠陥検査結果、過去の欠陥発生密度分布、等を用いることができる。
また、図7に示すように、上記確率予測システムの後段に、欠陥発生危険度判定部、警告システムを加えてもよい。上記危険度がある閾値を超えた場合、当該ウエハ及び後続ウエハの処理を中断し、当該箇所の欠陥発生確率が増大した原因を調査して、適宜プロセス条件等にフィードバックする。
又、上記危険度が一定以上の部分領域に対しては、図7に示すように、電子線欠陥検査装置等による全パターン検査を行い、上記部分領域内の欠陥数、欠陥密度を実測することが好ましい。さらに、予測精度検証部において予測精度の検証を行い、実測された欠陥密度が予測値とずれた場合、欠陥確率モデルにフィードバックする。フィードバックには例えばベイズ統計モデルを適用することができる。量産プロセスにおいて、これを継続することにより、自己学習により予測モデルが高精度化する。
801…電子源、802…引出電極、803…電子ビーム、804…コンデンサレンズ、805…走査偏向器、806…対物レンズ、807…試料室、808…試料台、809…試料、810…電子、811…二次電子、812…変換電極、813…検出器、814…制御装置、815…A/D変換器、816…画像処理部、817…CPU、818…画像メモリ、819…記憶媒体、820…ワークステーション

Claims (17)

  1. 計測ツールによって得られた計測データから、ウエハの欠陥の発生確率を推定するためのコンピューターで実行される方法を実施するためのコンピューターシステムで実行可能なプログラム命令を保存する非一時的コンピューター可読媒体であって、
    コンピューターで実行される方法は、
    ウエハ上の複数の計測点の計測データに基づいて求められる、パターンのエッジが、第1の位置に存在する確率に関する第1のデータを、取得或いは発生し、
    前記エッジが前記第1の位置にあるときに、当該第1の位置とは異なる第2の位置と前記第1の位置を含む領域を欠陥が覆う確率に関する第2のデータを取得或いは発生し、
    前記第1のデータと前記第2のデータの積に基づいて、前記欠陥の発生確率を予測する。
  2. 請求項1において、
    前記第1のデータは、前記第1の位置の設計パターンに対する相対的なエッジ位置座標をx_edgeとしたときの、前記複数の計測点の計測データから求められるx_edgeの頻度分布、P1(x_edge)である非一時的コンピューター可読媒体。
  3. 請求項1において、
    前記第2のデータは、前記エッジがx_edgeにあるときに、当該x_edgeと、前記第2の位置xとの間に連続して膜欠陥を生じる確率関数P2(x,x_edge)である非一時的コンピューター可読媒体。
  4. 請求項3において、
    前記ウエハにパターンを投影する露光装置によって得られる光強度分布から、単位面積当たりの局所的な膜欠陥確率P3を求め、当該P3について、前記x_edgeから前記第2の位置xの範囲の直積を求めることによって、前記確率関数P2を導出する非一時的コンピューター可読媒体。
  5. 請求項3において、
    前記ウエハにパターンを投影する露光装置によって、前記ウエハ上に形成されたレジスト膜に投影された光の光子により発生した光電子が、上記レジスト膜の前記_edgeから前記第2の位置xの間に近接して2次電子を発生する確率を求めることによって、前記確率関数P2を導出する非一時的コンピューター可読媒体。
  6. 請求項1において、
    前記欠陥発生確率の推定を行うウエハとは異なるウエハ上に形成された複数のパターンの計測に基づいて、前記第1のデータと、前記欠陥の発生確率を導出し、当該異なるウエハに関する第1のデータと、欠陥の発生分布に基づいて、前記第2のデータを導出する非一時的コンピューター可読媒体。
  7. 請求項6において、
    前記第1のデータは、前記第1の位置の設計パターンに対する相対的なエッジ位置座標をx_edgeとしたときの、前記複数の計測点の計測データから求められるx_edgeの頻度分布、P1(x_edge)である非一時的コンピューター可読媒体。
  8. 請求項において、
    前記第2のデータは、前記エッジがx_edgeにあるときに、当該x_edgeと、前記第2の位置との間に連続して膜欠陥を生じる確率関数P2(x,x_edge)である非一時的コンピューター可読媒体。
  9. 請求項6において、
    前記欠陥発生確率の推定を行うウエハとは異なるウエハは前記欠陥発生確率の推定を行うウエハとは、製造条件の異なるウエハである非一時的コンピューター可読媒体。
  10. 請求項において、
    前記P1(x_edge)を入力、前記欠陥の発生確率を出力とする教師データを用いて学習されたパラメータを前記第2のデータとして中間層に備えた学習モデルを生成し、当該学習モデルに、前記P1(x_edge)を入力することによって、前記欠陥の発生確率を出力する非一時的コンピューター可読媒体。
  11. ウエハ上の欠陥の発生確率を予測するために構成されたシステムであって、
    ウエハに対してビームを照射することによって得られる信号の検出に基づいて、出力を発生するように構成された計測ツールと、
    ウエハ上の複数の計測点の計測データに基づいて求められる、パターンのエッジが、第1の位置に存在する確率に関する第1のデータを、取得或いは発生し、
    前記エッジが前記第1の位置にあるときに、当該第1の位置とは異なる第2の位置と前記第1の位置との間の領域を膜欠陥が覆う確率に関する第2のデータを取得或いは発生し、
    前記第1のデータと前記第2のデータの積を求めるように構成されたコンピューターを含むシステム。
  12. 請求項11において、
    前記第1のデータは、前記第1の位置の設計パターンに対する相対的なエッジ位置座標をx_edgeとしたときの、前記複数の計測点の計測データから求められるx_edgeの頻度分布、P1(x_edge)であるシステム。
  13. 請求項11において、
    前記第2のデータは、前記エッジがx_edgeにあるときに、当該x_edgeと、前記第2の位置xとの間に連続して欠陥を生じる確率関数P2(x,x_edge)であるシステム。
  14. 請求項13において、
    前記コンピューターは、前記ウエハにパターンを投影する露光装置によって得られる光強度分布から、単位面積当たりの局所的な膜欠陥確率P3を求め、当該P3について、前記x_edgeから前記第2の位置xの範囲の直積を求めることによって、前記確率関数P2を導出するシステム。
  15. 請求項13において、
    前記コンピューターは、前記ウエハにパターンを投影する露光装置によって、前記ウエハ上に形成されたレジスト膜に投影された光の光子により発生した光電子が、上記レジスト膜の前記x_edgeから前記第2の位置xの間に近接して2次電子を発生する確率を求めることによって、前記確率関数P2を導出するシステム。
  16. 請求項11において、
    前記コンピューターは、前記第1のデータと前記欠陥の発生確率の実測値を入力として、前記第2のデータを出力するように構成されたシステム
  17. 請求項12において、
    前記コンピューターは、前記第1のデータを入力、前記欠陥の発生確率を出力とする教師データを用いて学習されたパラメータを前記第2のデータとして中間層に備えた学習モデルを生成し、当該学習モデルに、前記第1のデータを入力することによって、前記欠陥の発生確率を出力するシステム。
JP2020556468A 2018-11-12 2018-11-12 欠陥の発生を推定するシステム、及びコンピューター可読媒体 Active JP7080992B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2018/041781 WO2020100180A1 (ja) 2018-11-12 2018-11-12 欠陥の発生を推定するシステム、及びコンピューター可読媒体

Publications (2)

Publication Number Publication Date
JPWO2020100180A1 JPWO2020100180A1 (ja) 2021-10-07
JP7080992B2 true JP7080992B2 (ja) 2022-06-06

Family

ID=70731370

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020556468A Active JP7080992B2 (ja) 2018-11-12 2018-11-12 欠陥の発生を推定するシステム、及びコンピューター可読媒体

Country Status (6)

Country Link
US (2) US11747291B2 (ja)
JP (1) JP7080992B2 (ja)
KR (2) KR102563023B1 (ja)
CN (1) CN112889140B (ja)
TW (2) TWI724619B (ja)
WO (1) WO2020100180A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102563023B1 (ko) * 2018-11-12 2023-08-03 주식회사 히타치하이테크 결함의 발생을 추정하는 시스템, 및 컴퓨터 가독 매체
WO2020176270A1 (en) * 2019-02-25 2020-09-03 Applied Materials Israel Ltd. System and method for detecting rare stochastic defects
US11423529B2 (en) 2020-02-18 2022-08-23 Applied Materials Isreal Ltd. Determination of defect location for examination of a specimen
KR102384189B1 (ko) 2020-09-25 2022-04-08 충북대학교 산학협력단 머신러닝에 기초한 반도체 불량 시점 예측 장치 및 방법
US11699623B2 (en) * 2020-10-14 2023-07-11 Applied Materials, Inc. Systems and methods for analyzing defects in CVD films
KR20230028647A (ko) * 2021-08-20 2023-03-02 삼성전자주식회사 Euv 리소그래피의 결함 예측 방법 및 이를 이용한 반도체 소자의 제조 방법
CN115167019A (zh) * 2022-07-15 2022-10-11 苏州华星光电技术有限公司 曝光缺陷感知装置、显示终端及曝光缺陷感知方法
WO2024213339A1 (en) * 2023-04-11 2024-10-17 Asml Netherlands B.V. Method for efficient dynamic sampling plan generation and accurate probe die loss projection

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007140128A (ja) 2005-11-18 2007-06-07 Dainippon Printing Co Ltd 基板選択装置
JP2010034402A (ja) 2008-07-30 2010-02-12 Toshiba Corp パターン形状予測方法
JP2010243283A (ja) 2009-04-03 2010-10-28 Hitachi High-Technologies Corp 欠陥検査方法および欠陥検査装置
JP2013236087A (ja) 2012-04-19 2013-11-21 Applied Materials Israel Ltd トポグラフィック属性を利用した欠陥分類
JP2014240765A (ja) 2013-06-11 2014-12-25 株式会社日立ハイテクノロジーズ ラインパターンの形状評価方法及びその装置
JP2017096625A (ja) 2014-02-21 2017-06-01 株式会社日立ハイテクノロジーズ パターン測定装置、及びコンピュータープログラム
JP2018056327A (ja) 2016-09-29 2018-04-05 株式会社日立ハイテクノロジーズ パターン評価装置及びコンピュータープログラム

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3660763B2 (ja) * 1996-06-26 2005-06-15 株式会社日立製作所 被検査パターンの検査方法及び製造プロセス診断方法並びに半導体基板の製造方法
JP3808817B2 (ja) * 2002-09-05 2006-08-16 株式会社東芝 マスク欠陥検査方法、半導体装置の製造方法、マスク欠陥検査装置、欠陥影響度マップ作成方法およびプログラム
JP4993934B2 (ja) * 2006-03-31 2012-08-08 Hoya株式会社 パターン欠陥検査方法、フォトマスクの製造方法、及び表示デバイス用基板の製造方法
JP5330019B2 (ja) * 2009-02-18 2013-10-30 ルネサスエレクトロニクス株式会社 マスクパターンの検査方法およびマスクパターン検査装置
US8953869B2 (en) * 2012-06-14 2015-02-10 Kla-Tencor Corporation Apparatus and methods for inspecting extreme ultra violet reticles
JP2014130077A (ja) * 2012-12-28 2014-07-10 Hitachi High-Technologies Corp パターン形状評価方法、半導体装置の製造方法及びパターン形状評価装置
WO2015121127A1 (en) * 2014-02-11 2015-08-20 Asml Netherlands B.V. Model for calculating a stochastic variation in an arbitrary pattern
US10197908B2 (en) * 2016-06-21 2019-02-05 Lam Research Corporation Photoresist design layout pattern proximity correction through fast edge placement error prediction via a physics-based etch profile modeling framework
WO2018042581A1 (ja) * 2016-09-01 2018-03-08 株式会社 日立ハイテクノロジーズ パターン計測装置、及びコンピュータープログラム
CN110050230B (zh) * 2016-12-02 2021-06-11 Asml荷兰有限公司 用于估计随机变量的模型
US10474042B2 (en) * 2017-03-22 2019-11-12 Kla-Tencor Corporation Stochastically-aware metrology and fabrication
CN108335990B (zh) * 2018-01-31 2021-07-27 中国科学院微电子研究所 一种定位工艺缺陷的方法及装置
KR102563023B1 (ko) * 2018-11-12 2023-08-03 주식회사 히타치하이테크 결함의 발생을 추정하는 시스템, 및 컴퓨터 가독 매체

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007140128A (ja) 2005-11-18 2007-06-07 Dainippon Printing Co Ltd 基板選択装置
JP2010034402A (ja) 2008-07-30 2010-02-12 Toshiba Corp パターン形状予測方法
JP2010243283A (ja) 2009-04-03 2010-10-28 Hitachi High-Technologies Corp 欠陥検査方法および欠陥検査装置
JP2013236087A (ja) 2012-04-19 2013-11-21 Applied Materials Israel Ltd トポグラフィック属性を利用した欠陥分類
JP2014240765A (ja) 2013-06-11 2014-12-25 株式会社日立ハイテクノロジーズ ラインパターンの形状評価方法及びその装置
JP2017096625A (ja) 2014-02-21 2017-06-01 株式会社日立ハイテクノロジーズ パターン測定装置、及びコンピュータープログラム
JP2018056327A (ja) 2016-09-29 2018-04-05 株式会社日立ハイテクノロジーズ パターン評価装置及びコンピュータープログラム

Also Published As

Publication number Publication date
TWI807258B (zh) 2023-07-01
TW202018435A (zh) 2020-05-16
JPWO2020100180A1 (ja) 2021-10-07
CN112889140B (zh) 2023-06-02
US20230333033A1 (en) 2023-10-19
TW202135194A (zh) 2021-09-16
KR102453077B1 (ko) 2022-10-11
KR20210052531A (ko) 2021-05-10
KR102563023B1 (ko) 2023-08-03
US11747291B2 (en) 2023-09-05
CN112889140A (zh) 2021-06-01
KR20220140033A (ko) 2022-10-17
US12111272B2 (en) 2024-10-08
US20210396692A1 (en) 2021-12-23
WO2020100180A1 (ja) 2020-05-22
TWI724619B (zh) 2021-04-11

Similar Documents

Publication Publication Date Title
JP7080992B2 (ja) 欠陥の発生を推定するシステム、及びコンピューター可読媒体
US11802763B2 (en) Pattern measurement device and computer program
JP2020511003A (ja) オーバレイ計量データの確率論的挙動の影響の判別
KR101720992B1 (ko) 검사 장치 및 검사 장치 시스템
WO2019145278A1 (en) Pre-scan feature determination methods and systems
WO2018061135A1 (ja) パターン計測装置、及びコンピュータープログラム
US20170336717A1 (en) Exposure Condition Evaluation Device
JP2024508077A (ja) 確率的プロセスウィンドウの検出
US20220285226A1 (en) Apparatus and methods for three dimensional reticle defect smart repair
US20230019371A1 (en) Measurement Method, Measurement System, and Non-Transitory Computer Readable Medium
JP6581835B2 (ja) 半導体デバイスの評価条件設定方法、及び評価条件設定装置
TWI815419B (zh) 用於判定與微影製程相關之隨機度量之方法
TW201945853A (zh) 用於判定關於由微影製程形成之特徵之經校正尺寸參數值的方法及相關設備
US20240369944A1 (en) Method for determining a stochastic metric relating to a lithographic process
CN117255972A (zh) 用于确定与光刻过程有关的随机度量的方法
CN118302663A (zh) 筛选边缘安置均匀性晶片随机指标
JP2015115193A (ja) 走査型電子顕微鏡装置の評価方法、及び走査型電子顕微鏡装置の評価装置
JP2004040038A (ja) マスクの洗浄方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210331

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210331

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220510

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220525

R150 Certificate of patent or registration of utility model

Ref document number: 7080992

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150