JP2009092454A - Device and method for analyzing multilayer film - Google Patents

Device and method for analyzing multilayer film Download PDF

Info

Publication number
JP2009092454A
JP2009092454A JP2007261840A JP2007261840A JP2009092454A JP 2009092454 A JP2009092454 A JP 2009092454A JP 2007261840 A JP2007261840 A JP 2007261840A JP 2007261840 A JP2007261840 A JP 2007261840A JP 2009092454 A JP2009092454 A JP 2009092454A
Authority
JP
Japan
Prior art keywords
wave number
wavelength
sample
reflectance
layer structure
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.)
Granted
Application number
JP2007261840A
Other languages
Japanese (ja)
Other versions
JP4834847B2 (en
Inventor
Makoto Okawachi
真 大川内
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.)
Otsuka Electronics Co Ltd
Original Assignee
Otsuka Electronics Co Ltd
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 Otsuka Electronics Co Ltd filed Critical Otsuka Electronics Co Ltd
Priority to JP2007261840A priority Critical patent/JP4834847B2/en
Publication of JP2009092454A publication Critical patent/JP2009092454A/en
Application granted granted Critical
Publication of JP4834847B2 publication Critical patent/JP4834847B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a multilayer film analyzer and a multilayer film analyzing method capable of measuring the thickness of a multilayer film sample having wavelength dependence with higher accuracy. <P>SOLUTION: A wavenumber conversion section 721 calculates a wavenumber K<SB>1</SB>(λ) and a wavenumber-converted reflectivity R<SB>1</SB>'(λ)(=R(λ)/(1-R(λ))) successively for each wavelength stored in a buffer section 71, and outputs them to a buffer section 731. The buffer section 731 stores those wavenumbers K<SB>1</SB>(λ) and reflectivities R<SB>1</SB>'(λ) output successively from the conversion section 721, while associating them with each other. A Fourier conversion section 741 applies Fourier conversion to a wavenumber-converted reflectivity R<SB>1</SB>'(K<SB>1</SB>) stored in the buffer section 731 as to the wavenumber K<SB>1</SB>to calculate a power spectrum P<SB>1</SB>. A peak search section 751 searches for peaks appearing in the power spectrum P<SB>1</SB>calculated by the conversion section 741, and acquires a film thickness corresponding to each peak position. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

この発明は、多層膜解析装置および多層膜解析方法に関し、より特定的には膜厚をより高い精度で測定する技術に関する。   The present invention relates to a multilayer film analysis apparatus and a multilayer film analysis method, and more specifically to a technique for measuring a film thickness with higher accuracy.

従来から、機能性樹脂フィルムや半導体基板といった多層膜試料の膜厚を測定する方法として、反射率スペクトルを用いる方法が知られている。より具体的には、まず、白色光を多層膜試料に照射し、その反射光の反射率スペクトルを取得する。そして、取得された反射率スペクトルのデータに対して、ゼロフィリング処理やフィルタ処理などを行なった後、フーリエ変換を用いてパワースペクトルを取得する。さらに、このパワースペクトルに現れるピークから光学的膜厚nd(但し、nは膜の屈折率、dは膜厚)を算出し、この光学的膜厚ndを既知の屈折率nで除することにより、幾何学的な膜厚dを算出する。このような膜厚測定方法は、たとえば、特開2005−308394号公報(特許文献1)、特開2005−37315号公報(特許文献2)、特開平07−294220号公報(特許文献3)、特開平10−311708号公報(特許文献4)、特開2003−130614号公報(特許文献5)などに記載されている。   Conventionally, a method using a reflectance spectrum is known as a method for measuring the film thickness of a multilayer film sample such as a functional resin film or a semiconductor substrate. More specifically, first, the multilayer film sample is irradiated with white light, and the reflectance spectrum of the reflected light is acquired. And after performing a zero filling process, a filter process, etc. with respect to the data of the acquired reflectance spectrum, a power spectrum is acquired using a Fourier transform. Further, the optical film thickness nd (where n is the refractive index of the film and d is the film thickness) is calculated from the peak appearing in the power spectrum, and this optical film thickness nd is divided by the known refractive index n. The geometric film thickness d is calculated. Such film thickness measuring methods are disclosed in, for example, JP-A-2005-308394 (Patent Document 1), JP-A-2005-37315 (Patent Document 2), JP-A-07-294220 (Patent Document 3), It describes in Unexamined-Japanese-Patent No. 10-311708 (patent document 4), Unexamined-Japanese-Patent No. 2003-130614 (patent document 5), etc.

一般的に、薄膜を構成する物質(たとえば、誘電体など)は、その屈折率に波長依存性を有する。しかしながら、上述の先行技術に開示される方法では、多層膜試料を構成する薄膜の屈折率が波長依存性を持たないと仮定した上で、パワースペクトルを算出している。そのため、その屈折率に波長依存性をもつ多層膜試料を測定する場合には、この屈折率に起因した誤差が発生する。   In general, a substance constituting a thin film (for example, a dielectric) has wavelength dependency on its refractive index. However, in the method disclosed in the above prior art, the power spectrum is calculated on the assumption that the refractive index of the thin film constituting the multilayer film sample does not have wavelength dependency. Therefore, when measuring a multilayer film sample having a wavelength dependency on the refractive index, an error due to the refractive index occurs.

これに対して、特開2001−227916号公報(特許文献6)には、パワースペクトルに基づいて算出される光学的膜厚ndに対して、屈折率nの波長依存性を考慮して補正を行なう方法が開示されている。
特開2005−308394号公報 特開2005−37315号公報 特開平07−294220号公報 特開平10−311708号公報 特開2003−130614号公報 特開2001−227916号公報
On the other hand, Japanese Patent Laid-Open No. 2001-227916 (Patent Document 6) corrects the optical film thickness nd calculated based on the power spectrum in consideration of the wavelength dependence of the refractive index n. A method of performing is disclosed.
JP 2005-308394 A JP 2005-37315 A JP 07-294220 A Japanese Patent Laid-Open No. 10-311708 JP 2003-130614 A JP 2001-227916 A

しかしながら、特開2001−227916号公報(特許文献6)に開示される方法においても、その式中に現れるいくつかの項(たとえば、屈折率n,dn/dk,k=2π/λなど)は依然として波長依存性を有しているため、波長依然性に起因する誤差を完全に取り除くことはできないという問題があった。 However, even in the method disclosed in Japanese Patent Laid-Open No. 2001-227916 (Patent Document 6), some terms appearing in the equation (for example, refractive index n 2 , dn f / dk, k = 2π / λ, etc.) ) Still has wavelength dependency, and there is a problem that errors due to wavelength dependency cannot be completely removed.

この発明は、このような問題点を解決するためになされたものであって、その目的は、波長依存性を有する多層膜試料の膜厚をより高い精度を測定することが可能な多層膜解析装置および多層膜解析方法を提供することである。   The present invention has been made to solve such problems, and its purpose is to analyze a multilayer film that can measure the film thickness of a multilayer film sample having wavelength dependency with higher accuracy. An apparatus and a multilayer analysis method are provided.

この発明のある局面に従う多層膜解析装置は、基板上に少なくとも1つの膜が形成された試料に対して所定の波長範囲をもつ測定光を照射する光源と、試料で反射された光または試料を透過した光に基づいて、反射率または透過率の波長分布特性を取得する分光部と、波長分布特性における各波長とその波長における値との対応関係を、各波長についての波数と所定の関係式に従って算出される変換値との対応関係に変換する波数特性変換手段と、波数特性変換手段によってそれぞれ変換された対応関係からなる波数分布特性に対して、波数について周波数変換する周波数変換手段と、周波数変換後の波数分布特性に現れるピークに基づいて、試料における膜厚を算出する膜厚算出手段とを備える。   A multilayer film analyzer according to an aspect of the present invention includes a light source that irradiates measurement light having a predetermined wavelength range to a sample in which at least one film is formed on a substrate, and light or sample reflected by the sample. Based on the transmitted light, the spectroscopic unit that acquires the reflectance or transmittance wavelength distribution characteristic, the correspondence between each wavelength in the wavelength distribution characteristic and the value at that wavelength, the wave number for each wavelength and a predetermined relational expression A wave number characteristic converting means for converting to a corresponding relationship with a conversion value calculated according to the above, a frequency converting means for converting the frequency of the wave number with respect to the wave number distribution characteristic having a corresponding relationship converted by the wave number characteristic converting means, and a frequency And a film thickness calculating means for calculating the film thickness of the sample based on the peak appearing in the wave number distribution characteristics after the conversion.

好ましくは、波数特性変換手段は、予め入力される試料に含まれる膜の屈折率と対応の波長とに基づいて波数を算出する。   Preferably, the wave number characteristic conversion means calculates the wave number based on the refractive index of the film included in the sample inputted in advance and the corresponding wavelength.

好ましくは、関係式は、各波長における反射率または透過率の値を周波数変換手段での周波数変換に係る位相因子に対して線形化するための関数である。   Preferably, the relational expression is a function for linearizing the value of reflectance or transmittance at each wavelength with respect to a phase factor related to frequency conversion in the frequency conversion means.

さらに好ましくは、波長分布特性は、試料で反射された光に基づいて取得される反射率Rの波長分布特性であり、関係式は、R/(1−R)である。   More preferably, the wavelength distribution characteristic is a wavelength distribution characteristic of the reflectance R acquired based on the light reflected by the sample, and the relational expression is R / (1-R).

好ましくは、波数特性変換手段は、試料に複数の膜が形成されている場合に、複数の膜のそれぞれについて波数分布特性を算出し、膜厚算出手段は、複数の膜についての周波数変換後の波波数分布特性のそれぞれに現れるピークに基づいて、複数の膜の各々の膜厚を算出する。   Preferably, the wave number characteristic conversion unit calculates the wave number distribution characteristic for each of the plurality of films when the sample has a plurality of films, and the film thickness calculation unit calculates the frequency after frequency conversion for the plurality of films. Based on the peak appearing in each of the wave number distribution characteristics, the thickness of each of the plurality of films is calculated.

好ましくは、少なくとも屈折率および膜厚が規定された層構造情報を複数格納する層構造格納手段と、複数の層構造情報のうち少なくとも1つの層構造情報に基づいて、層構造候補を順次選択する選択手段と、層構造候補に対応する層構造情報に基づいて、層構造候補における反射率または透過率の波長分布特性を算出する波長分布特性算出手段と、層構造候補の波長分布特性から、波数特性変換手段および周波数変換手段と同一の処理によって周波数変換後の波数分布特性を算出する波数分布特性算出手段と、実測された波長分布特性から算出された周波数変換後の波数分布特性と、波数分布特性算出手段によって算出された周波数変換後の波数分布特性との間の相関値を算出する相関値算出手段と、相関値に基づいて、選択手段が順次選択する層構造候補のうち実測値に最も近いものを決定する決定手段とをさらに備える。   Preferably, at least one layer structure candidate is sequentially selected based on at least one layer structure information out of the plurality of layer structure information, and a layer structure storage means for storing a plurality of layer structure information in which at least a refractive index and a film thickness are defined. From the selection means, the wavelength distribution characteristic calculating means for calculating the wavelength distribution characteristic of the reflectance or transmittance in the layer structure candidate based on the layer structure information corresponding to the layer structure candidate, and the wave number from the wavelength distribution characteristic of the layer structure candidate Wave number distribution characteristic calculation means for calculating the wave number distribution characteristic after frequency conversion by the same processing as the characteristic conversion means and frequency conversion means, the wave number distribution characteristic after frequency conversion calculated from the actually measured wavelength distribution characteristic, and the wave number distribution Correlation value calculation means for calculating a correlation value between the frequency distribution characteristics after frequency conversion calculated by the characteristic calculation means, and the selection means sequentially selects based on the correlation value. Further comprising determination means for determining the closest match to the measured value of the layer structure candidates.

この発明の別の局面に従う多層膜解析方法は、基板上に少なくとも1つの膜が形成された試料に対して所定の波長範囲をもつ測定光を照射するステップと、試料で反射された光または試料を透過した光に基づいて、反射率または透過率の波長分布特性を取得するステップと、波長分布特性における各波長とその波長における値との対応関係を、各波長についての波数と所定の関係式に従って算出される変換値との対応関係に変換するステップと、対応関係に変換するステップにおいてそれぞれ変換された対応関係からなる波数分布特性に対して、波数について周波数変換するステップと、周波数変換後の波数分布特性に現れるピークに基づいて、試料における膜厚を算出するステップとを備える。   A multilayer film analysis method according to another aspect of the present invention includes a step of irradiating a sample having at least one film formed on a substrate with measurement light having a predetermined wavelength range, and light reflected from the sample or the sample. The step of obtaining the wavelength distribution characteristic of reflectance or transmittance based on the light transmitted through the light, the correspondence between each wavelength in the wavelength distribution characteristic and the value at that wavelength, the wave number for each wavelength and a predetermined relational expression Converting to the corresponding relationship with the conversion value calculated in accordance with the step, the step of converting the frequency of the wave number with respect to the wave number distribution characteristic composed of the corresponding relationship converted in the step of converting to the corresponding relationship, Calculating the film thickness of the sample based on the peak appearing in the wave number distribution characteristics.

この発明によれば、波長依存性を有する多層膜試料の膜厚をより高い精度を測定することができる。   According to this invention, it is possible to measure the film thickness of a multilayer film sample having wavelength dependency with higher accuracy.

この発明の実施の形態について、図面を参照しながら詳細に説明する。なお、図中の同一または相当部分については、同一符号を付してその説明は繰返さない。   Embodiments of the present invention will be described in detail with reference to the drawings. Note that the same or corresponding parts in the drawings are denoted by the same reference numerals and description thereof will not be repeated.

図1は、この発明の実施の形態に従う多層膜解析装置100の概略構成図である。
本実施の形態に従う多層膜解析装置100は、代表的に、多層膜試料で反射された反射光の波長分布特性(以下「スペクトル」とも称す。)に基づいて、多層膜試料における各層の膜厚測定や層構造の解析を行なう。なお、反射光のスペクトルに代えて、多層膜試料を透過した光のスペクトル(透過光スペクトル)を用いてもよい。
FIG. 1 is a schematic configuration diagram of a multilayer film analyzing apparatus 100 according to an embodiment of the present invention.
Multilayer film analysis apparatus 100 according to the present embodiment typically has a film thickness of each layer in the multilayer film sample based on the wavelength distribution characteristics (hereinafter also referred to as “spectrum”) of the reflected light reflected by the multilayer film sample. Perform measurement and analysis of layer structure. Instead of the spectrum of the reflected light, the spectrum of the light transmitted through the multilayer film sample (transmitted light spectrum) may be used.

本明細書において多層膜試料とは、基板上に少なくとも1つの膜が形成された試料を意味する。代表的な多層膜試料としては、半導体基板、ガラス基板、サファイア基板、石英基板、フィルムなどの表面に薄膜が形成(コーティング)された多層薄膜が挙げられる。より具体的には、薄膜形成されたガラス基板は、液晶ディスプレイ(LCD:Liquid Crystal Display)やプラズマディスプレイ(PDP:Plasma Display Panel)などのフラットパネルディスプレイ(FPD:Flat Panel Display)のディスプレイ部として使用される。また、薄膜形成されたサファイア基板は、窒化物半導体(GaN:Gallium Nitride)系のLED(Light Emitting Diode)やLD(Laser Diode)として使用される。また、薄膜形成された石英基板は、各種の光学フィルタや光学部品およびプロジェクション液晶などに使用される。なお、以下では、多層膜試料を単に「試料」とも記す。   In this specification, a multilayer film sample means a sample in which at least one film is formed on a substrate. Typical multilayer film samples include a multilayer thin film in which a thin film is formed (coated) on the surface of a semiconductor substrate, a glass substrate, a sapphire substrate, a quartz substrate, a film, or the like. More specifically, the thin-film-formed glass substrate is used as a display part of a flat panel display (FPD) such as a liquid crystal display (LCD) or a plasma display panel (PDP). Is done. The thin-film sapphire substrate is used as a nitride semiconductor (GaN: Gallium Nitride) LED (Light Emitting Diode) or LD (Laser Diode). The quartz substrate formed with a thin film is used for various optical filters, optical components, projection liquid crystals, and the like. Hereinafter, the multilayer film sample is also simply referred to as “sample”.

図1を参照して、多層膜解析装置100は、測定用光源10と、コリメートレンズ12と、カットフィルタ14と、結像レンズ16,36と、絞り部18と、ビームスプリッタ20,30と、観察用光源22と、光ファイバ24と、出射部26と、ピンホールミラー32と、軸変換ミラー34と、観察用カメラ38と、表示部39と、対物レンズ40とステージ50と、可動機構52と、分光測定部60と、データ処理部70とを備える。   Referring to FIG. 1, the multilayer film analyzing apparatus 100 includes a measurement light source 10, a collimating lens 12, a cut filter 14, imaging lenses 16 and 36, a diaphragm unit 18, beam splitters 20 and 30, Observation light source 22, optical fiber 24, emitting section 26, pinhole mirror 32, axis conversion mirror 34, observation camera 38, display section 39, objective lens 40, stage 50, and movable mechanism 52. And a spectroscopic measurement unit 60 and a data processing unit 70.

測定用光源10は、多層膜試料の反射光のスペクトルを取得するために、所定の波長範囲をもつ測定光を発生する光源であり、その波長幅が比較的広い白色光光源が用いられる。代表的に、測定用光源10は、重水素ランプ(Dランプ)やタングステンランプ、またはそれらの組合せからなる。 The measurement light source 10 is a light source that generates measurement light having a predetermined wavelength range in order to acquire a spectrum of reflected light from the multilayer film sample, and a white light source having a relatively wide wavelength range is used. Typically, the measurement light source 10 is composed of a deuterium lamp (D 2 lamp) or a tungsten lamp or a combination thereof.

コリメートレンズ12と、カットフィルタ14と、結像レンズ16と、絞り部18とは、測定用光源10とビームスプリッタ30とを結ぶ光軸AX2上に配置され、測定用光源10から出射される測定光を光学的に調整する。   The collimating lens 12, the cut filter 14, the imaging lens 16, and the aperture unit 18 are arranged on the optical axis AX <b> 2 connecting the measurement light source 10 and the beam splitter 30, and are measured from the measurement light source 10. Adjust light optically.

具体的には、コリメートレンズ12は、測定用光源10からの測定光が最初に入射する光学部品であり、拡散光線として伝播する測定光を屈折させて平行光線に変換する。コリメートレンズ12を通過した測定光はカットフィルタ14に入射する。カットフィルタ14は、測定光に含まれる不要な波長成分を遮断する。代表的に、カットフィルタ14は、ガラス基板などに蒸着された多層膜によって形成される。結像レンズ16は、測定光のビーム径を調整するために、カットフィルタ14を通過した測定光を平行光線から収束光線に変換する。結像レンズ16を通過した測定光は絞り部18に入射する。絞り部18は、測定光の光量を所定量に調整した上でビームスプリッタ30へ出射する。好ましくは、絞り部18は、結像レンズ16によって変換された測定光の結像位置に配置される。なお、絞り部18の絞り量は、試料に入射する測定光の被写界深度や必要な光強度などに応じて適宜設定される。   Specifically, the collimating lens 12 is an optical component on which the measurement light from the measurement light source 10 first enters, and refracts the measurement light propagating as a diffused light and converts it into parallel light. The measurement light that has passed through the collimating lens 12 enters the cut filter 14. The cut filter 14 blocks unnecessary wavelength components included in the measurement light. Typically, the cut filter 14 is formed of a multilayer film deposited on a glass substrate or the like. The imaging lens 16 converts the measurement light that has passed through the cut filter 14 from parallel light into convergent light in order to adjust the beam diameter of the measurement light. The measurement light that has passed through the imaging lens 16 is incident on the diaphragm 18. The diaphragm 18 adjusts the amount of measurement light to a predetermined amount and then emits the light to the beam splitter 30. Preferably, the diaphragm unit 18 is disposed at the imaging position of the measurement light converted by the imaging lens 16. The diaphragm amount of the diaphragm unit 18 is appropriately set according to the depth of field of the measurement light incident on the sample, the required light intensity, and the like.

一方、観察用光源22は、試料への焦点合わせや測定位置の確認に使用される観察光を生成する光源である。そして、観察用光源22が発生する観察光は、試料で反射可能な波長を含むように選択される。観察用光源22は、光ファイバ24を介して出射部26と接続されており、観察用光源22で生成された観察光は、光導波路である光ファイバ24を伝播した後に出射部26からビームスプリッタ20へ向けて出射される。   On the other hand, the observation light source 22 is a light source that generates observation light used for focusing on a sample and confirming a measurement position. The observation light generated by the observation light source 22 is selected so as to include a wavelength that can be reflected by the sample. The observation light source 22 is connected to the emission unit 26 via the optical fiber 24, and the observation light generated by the observation light source 22 propagates through the optical fiber 24, which is an optical waveguide, and then travels from the emission unit 26 to the beam splitter. It is emitted toward 20.

出射部26は、試料に所定の観察基準像が投射されるように、観察用光源22で生成された観察光の一部をマスクするマスク部26aを含む。この観察基準像は、その表面に何らの模様(パターン)も形成されていない試料(代表的に、透明なガラス基板など)に対しても、ユーザによる焦点合わせを容易化するためのものである。なお、レチクル像の形状はいずれであってもよいが、一例として同心円状や十字状のパターンなどを用いることができる。   The emission unit 26 includes a mask unit 26a that masks a part of the observation light generated by the observation light source 22 so that a predetermined observation reference image is projected onto the sample. This observation reference image is for facilitating user focusing even on a sample (typically a transparent glass substrate or the like) on which no pattern is formed. . Note that the shape of the reticle image may be any, but as an example, a concentric or cruciform pattern may be used.

すなわち、観察用光源22で生成された直後の観察光のビーム断面における光強度(光量)は略均一であるが、マスク部26aがこの観察光の一部をマスク(遮へい)することで、観察光は、そのビーム断面において光強度が略ゼロである領域(影領域)が形成される。この影領域が観察基準像として試料に投射される。   That is, the light intensity (light quantity) in the beam cross section of the observation light immediately after being generated by the observation light source 22 is substantially uniform, but the mask portion 26a masks (shields) a part of the observation light, thereby observing. The light has a region (shadow region) where the light intensity is substantially zero in the beam cross section. This shadow region is projected onto the sample as an observation reference image.

ステージ50は、試料を配置するための試料台であり、その配置面は平坦に形成される。このステージ50は、一例として機械的に連結された可動機構52によって、3方向(X方向・Y方向・Z方向)に自在に駆動される。可動機構52は、代表的に3軸分のサーボモータと、各サーボモータを駆動するためのサーボドライバとを含んで構成される。そして、可動機構52は、ユーザまたは図示しない制御装置などからのステージ位置指令に応答してステージ50を駆動する。このステージ50の駆動によって、試料と後述する対物レンズ40との間の位置関係が変更される。   The stage 50 is a sample stage for arranging samples, and the arrangement surface thereof is formed flat. The stage 50 is freely driven in three directions (X direction, Y direction, and Z direction) by a movable mechanism 52 that is mechanically coupled as an example. The movable mechanism 52 typically includes a servo motor for three axes and a servo driver for driving each servo motor. The movable mechanism 52 drives the stage 50 in response to a stage position command from a user or a control device (not shown). The positional relationship between the sample and an objective lens 40 described later is changed by driving the stage 50.

対物レンズ40と、ビームスプリッタ20と、ビームスプリッタ30と、ピンホールミラー32とは、ステージ50の平坦面に垂直な方向に延伸する光軸AX1上に配置される。   The objective lens 40, the beam splitter 20, the beam splitter 30, and the pinhole mirror 32 are disposed on the optical axis AX1 that extends in a direction perpendicular to the flat surface of the stage 50.

ビームスプリッタ30は、測定用光源10で生成される測定光を反射することで、その伝播方向を光軸AX1の紙面下向きに変換する。また、ビームスプリッタ30は、光軸AX1を紙面上向きに伝播する試料からの反射光を透過させる。代表的に、ビームスプリッタ30はハーフミラーで構成される。   The beam splitter 30 reflects the measurement light generated by the measurement light source 10, thereby converting the propagation direction of the optical axis AX1 downward in the drawing. The beam splitter 30 transmits the reflected light from the sample propagating along the optical axis AX1 upward in the drawing. Typically, the beam splitter 30 is composed of a half mirror.

一方、ビームスプリッタ20は、観察用光源22で生成される観察光を反射することで、その伝播方向を光軸AX1の紙面下向きに変換する。同時に、ビームスプリッタ20は、光軸AX1を紙面下向きに伝播するビームスプリッタ30で反射された測定光を透過させる。すなわち、ビームスプリッタ20は、測定用光源10から集光光学系である対物レンズ40までの光学経路上の所定位置において観察光を注入する光注入部として機能する。このビームスプリッタ20で合成された測定光と観察光とは、対物レンズ40に入射する。また、ビームスプリッタ20は、光軸AX1を紙面上向きに伝播する試料からの反射光を透過させる。代表的に、ビームスプリッタ20はハーフミラーで構成される。   On the other hand, the beam splitter 20 reflects the observation light generated by the observation light source 22, thereby converting the propagation direction of the optical axis AX <b> 1 downward in the drawing. At the same time, the beam splitter 20 transmits the measurement light reflected by the beam splitter 30 that propagates along the optical axis AX1 downward in the drawing. That is, the beam splitter 20 functions as a light injection unit that injects observation light at a predetermined position on the optical path from the measurement light source 10 to the objective lens 40 that is a condensing optical system. The measurement light and the observation light synthesized by the beam splitter 20 enter the objective lens 40. The beam splitter 20 transmits the reflected light from the sample propagating along the optical axis AX1 upward in the drawing. Typically, the beam splitter 20 is composed of a half mirror.

対物レンズ40は、光軸AX1を紙面下向きに伝播する測定光および観察光を集光するための集光光学系である。すなわち、対物レンズ40は、試料またはその近接した位置で結像するように測定光および観察光を収束させる。また、対物レンズ40は、所定の倍率(たとえば、10倍,20倍,30倍,40倍など)を有する拡大レンズであり、試料の光学特性を測定する領域を対物レンズ40に入射する光のビーム断面に比較してより微小化できる。   The objective lens 40 is a condensing optical system for condensing measurement light and observation light that propagates along the optical axis AX1 downward in the drawing. That is, the objective lens 40 converges the measurement light and the observation light so as to form an image at the sample or a position close thereto. The objective lens 40 is a magnifying lens having a predetermined magnification (for example, 10 times, 20 times, 30 times, 40 times, etc.), and a region for measuring the optical characteristics of the sample is incident on the objective lens 40. Compared to the beam cross section, it can be further miniaturized.

また、対物レンズ40から試料に入射した測定光および観察光は、試料で反射され、光軸AX1上を紙面上向きに伝播する。この反射光は、対物レンズ40に透過した後、ビームスプリッタ20および30を透過してピンホールミラー32まで到達する。   Further, the measurement light and observation light incident on the sample from the objective lens 40 are reflected by the sample and propagate upward on the optical axis AX1. The reflected light passes through the objective lens 40 and then passes through the beam splitters 20 and 30 to reach the pinhole mirror 32.

ピンホールミラー32は、試料で生じる反射光のうち、測定反射光と観察反射光とを分離する光分離部として機能する。具体的には、ピンホールミラー32は、光軸AX1を紙面上向きに伝播する試料からの反射光を反射する反射面を含み、その反射面と光軸AX1との交点を中心とする穴あき部(ピンホール)32aが形成されている。このピンホール32aの大きさは、測定用光源10からの測定光が試料で反射されて生じる測定反射光の、ピンホールミラー32の位置におけるビーム径に比較して小さくなるように形成される。また、このピンホール32aは、それぞれ測定光および観察光が試料で反射されて生じる測定反射光および観察反射光の結像位置と一致するように配置される。このような構成によって、試料で生じた反射光は、ピンホール32aを通過して分光測定部60に入射する。一方、反射光の残部は、その伝播方向を変換されて軸変換ミラー34へ入射する。   The pinhole mirror 32 functions as a light separation unit that separates measurement reflected light and observation reflected light among the reflected light generated in the sample. Specifically, the pinhole mirror 32 includes a reflective surface that reflects the reflected light from the sample propagating along the optical axis AX1 upward in the drawing, and has a hole with a center at the intersection of the reflective surface and the optical axis AX1. (Pinhole) 32a is formed. The size of the pinhole 32a is formed so as to be smaller than the beam diameter at the position of the pinhole mirror 32 of the measurement reflected light generated when the measurement light from the measurement light source 10 is reflected by the sample. Further, the pinhole 32a is disposed so as to coincide with the imaging positions of the measurement reflected light and the observation reflected light generated by reflecting the measurement light and the observation light on the sample, respectively. With such a configuration, the reflected light generated in the sample passes through the pinhole 32a and enters the spectroscopic measurement unit 60. On the other hand, the remaining part of the reflected light is converted in its propagation direction and enters the axis conversion mirror 34.

分光測定部60は、ピンホールミラー32を通過した測定反射光のスペクトルを測定し、その測定結果をデータ処理部70へ出力する。より詳細には、分光測定部60は、回折格子(グレーティング)62と、検出部64と、カットフィルタ66と、シャッタ68とを含む。   The spectroscopic measurement unit 60 measures the spectrum of the reflected measurement light that has passed through the pinhole mirror 32 and outputs the measurement result to the data processing unit 70. More specifically, the spectroscopic measurement unit 60 includes a diffraction grating (grating) 62, a detection unit 64, a cut filter 66, and a shutter 68.

カットフィルタ66と、シャッタ68と、回折格子62とは、光軸AX1上に配置される。カットフィルタ66は、ピンホールを通過して分光測定部60に入射する測定反射光に含まれる測定範囲外の波長成分を制限するための光学フィルタであり、特に測定範囲外の波長成分を遮断する。シャッタ68は、検出部64をリセットするときなどに、検出部64に入射する光を遮断するために使用される。シャッタ68は、代表的に電磁力によって駆動する機械式のシャッタからなる。   The cut filter 66, the shutter 68, and the diffraction grating 62 are disposed on the optical axis AX1. The cut filter 66 is an optical filter for limiting a wavelength component outside the measurement range included in the measurement reflected light that passes through the pinhole and enters the spectroscopic measurement unit 60, and particularly blocks the wavelength component outside the measurement range. . The shutter 68 is used to block light incident on the detection unit 64 when the detection unit 64 is reset. The shutter 68 is typically a mechanical shutter that is driven by electromagnetic force.

回折格子62は、入射する測定反射光を分光した上で、各分光波を検出部64へ導く。具体的には、回折格子62は、反射型の回折格子であり、所定の波長間隔毎の回折波が対応する各方向に反射するように構成される。このような構成を有する回折格子62に測定反射波が入射すると、含まれる各波長成分は対応する方向に反射されて、検出部64の所定の検出領域に入射する。回折格子62は、代表的にフラットフォーカス型球面グレーティングからなる。   The diffraction grating 62 separates the incident measurement reflected light and guides each spectral wave to the detection unit 64. Specifically, the diffraction grating 62 is a reflection type diffraction grating, and is configured so that a diffracted wave for each predetermined wavelength interval is reflected in each corresponding direction. When the measurement reflected wave is incident on the diffraction grating 62 having such a configuration, each wavelength component included is reflected in a corresponding direction and is incident on a predetermined detection region of the detection unit 64. The diffraction grating 62 is typically composed of a flat focus type spherical grating.

検出部64は、試料の反射率スペクトルを測定するために、回折格子62で分光された測定反射光に含まれる各波長成分の光強度に応じた電気信号を出力する。検出部64は、代表的にフォトダイオードなどの検出素子をアレイ状に配置したフォトダイオードアレイや、マトリックス状に配置されたCCD(Charged Coupled Device)などからなる。   The detection unit 64 outputs an electrical signal corresponding to the light intensity of each wavelength component included in the measurement reflected light that has been dispersed by the diffraction grating 62 in order to measure the reflectance spectrum of the sample. The detection unit 64 typically includes a photodiode array in which detection elements such as photodiodes are arranged in an array, a CCD (Charged Coupled Device) arranged in a matrix, and the like.

なお、回折格子62および検出部64は、光学特性の測定波長範囲および測定波長間隔などに応じて適宜設計される。   The diffraction grating 62 and the detection unit 64 are appropriately designed according to the measurement wavelength range and the measurement wavelength interval of the optical characteristics.

データ処理部70は、検出部64によって取得された反射率スペクトルに対して、本発明に係る特徴的な処理を行なうことで、試料の膜厚を測定する。さらに、データ処理部70は、試料の層構造を解析することも可能である。膜厚の測定方法および層構造の解析方法については後述する。そして、データ処理部70は、測定した試料の膜厚や層構造などを光学特性として出力する。   The data processing unit 70 measures the film thickness of the sample by performing characteristic processing according to the present invention on the reflectance spectrum acquired by the detection unit 64. Furthermore, the data processing unit 70 can also analyze the layer structure of the sample. A method for measuring the film thickness and a method for analyzing the layer structure will be described later. Then, the data processing unit 70 outputs the measured film thickness and layer structure of the sample as optical characteristics.

一方、ピンホールミラー32で反射された観測反射光は光軸AX3に沿って伝播し、軸変換ミラー34へ入射する。軸変換ミラー34は、観測反射光の伝播方向を光軸AX3から光軸AX4に変換する。すると、観測反射光は、光軸AX4に沿って伝播し、観察用カメラ38へ入射する。   On the other hand, the observation reflected light reflected by the pinhole mirror 32 propagates along the optical axis AX3 and enters the axis conversion mirror 34. The axis conversion mirror 34 converts the propagation direction of the observation reflected light from the optical axis AX3 to the optical axis AX4. Then, the observation reflected light propagates along the optical axis AX4 and enters the observation camera 38.

観察用カメラ38は、観察反射光によって得られる反射像を取得する撮像部であり、代表的にはCCD(Charged Coupled Device)やCMOS(Complementary Metal Oxide Semiconductor)センサなどからなる。なお、観察用カメラ38は、代表的に可視帯域に感度をもつものであり、所定の測定範囲に感度をもつ検出部64とは異なる感度特性をもつ場合が多い。そして、観察用カメラ38は、観察反射光によって得られる反射像に応じた映像信号を表示部39へ出力する。表示部39は、観察用カメラ38からの映像信号に基づいて反射像を画面上に表示する。ユーザは、この表示部39に表示される反射像を目視して、試料に対する焦点合わせや測定位置の確認などを行なう。表示部39は、代表的に液晶ディスプレイ(LCD)などからなる。なお、観察用カメラ38および表示部39に代えて、ユーザが反射像を直接的に目視できるファインダーを設けてもよい。   The observation camera 38 is an imaging unit that acquires a reflected image obtained by the observation reflected light, and typically includes a CCD (Charged Coupled Device), a CMOS (Complementary Metal Oxide Semiconductor) sensor, or the like. Note that the observation camera 38 typically has sensitivity in the visible band, and often has sensitivity characteristics different from those of the detection unit 64 having sensitivity in a predetermined measurement range. Then, the observation camera 38 outputs a video signal corresponding to the reflected image obtained by the observation reflected light to the display unit 39. The display unit 39 displays a reflected image on the screen based on the video signal from the observation camera 38. The user visually observes the reflected image displayed on the display unit 39 to perform focusing on the sample and confirmation of the measurement position. The display unit 39 typically includes a liquid crystal display (LCD). Instead of the observation camera 38 and the display unit 39, a finder that allows the user to directly view the reflected image may be provided.

(膜厚測定処理その1)
以下では、データ処理部70における試料の膜厚を測定する処理について説明する。
(Thickness measurement process 1)
Below, the process which measures the film thickness of the sample in the data processing part 70 is demonstrated.

図2は、代表的な多層膜試料の断面模式図である。
説明を簡単化するために、図2に示すような、空気や真空などの雰囲気層(添え字0)、測定対象とする薄膜(添え字1)、基板層(添え字2)の3層からなる代表的な多層膜試料を考える。各層における屈折率を添え字iを用いて、屈折率nと表す。互いに異なる屈折率をもつ層の界面では光の反射が生じるため、屈折率の異なるi層とi+1層との間の各境界面でのP偏光成分およびS偏光成分の振幅反射率(Fresnel係数)は次の式のように表わすことができる。
FIG. 2 is a schematic cross-sectional view of a typical multilayer film sample.
In order to simplify the explanation, from the three layers shown in FIG. 2, an atmosphere layer such as air or vacuum (subscript 0), a thin film to be measured (subscript 1), and a substrate layer (subscript 2). Consider a typical multilayer film sample. With index i the refractive index of each layer, expressed as the refractive index n i. Since reflection of light occurs at the interface between layers having different refractive indexes, the amplitude reflectance (Fresnel coefficient) of the P-polarized component and the S-polarized component at each interface between the i layer and i + 1 layer having different refractive indexes. Can be expressed as:

Figure 2009092454
Figure 2009092454

ここで、φは、i層における入射角である。この入射角φは、以下のようなSnellの法則によって、最上層の雰囲気層(0層)における入射角から計算できる。 Here, φ i is an incident angle in the i layer. This incident angle φ i can be calculated from the incident angle in the uppermost atmosphere layer (0 layer) according to Snell's law as follows.

sinφ=Nsinφ
光が干渉可能な厚さをもつ薄膜内では、上式で表される反射率で反射する光が薄膜内を何度も往復する。そのため、表面で直接反射した光と薄膜内を多重反射した後の光との間ではその光路長が異なるため、位相が互いに異なったものとなり、薄膜表面において光の干渉が生じる。このような、各薄膜内における光の干渉効果を示すために、i層の薄膜における光の位相角βを導入すると、以下のように表わすことができる。
N 0 sinφ 0 = N i sinφ i
In a thin film having a thickness that allows light to interfere, the light reflected at the reflectance expressed by the above formula reciprocates in the thin film many times. For this reason, the light path length is different between the light directly reflected on the surface and the light after multiple reflection in the thin film, so that the phases are different from each other, and light interference occurs on the thin film surface. In order to show the light interference effect in each thin film, when the phase angle β i of light in the i-layer thin film is introduced, it can be expressed as follows.

Figure 2009092454
Figure 2009092454

ここで、dはi層の厚さを示し、λは入射光の波長を示す。
より単純化するために、多層膜試料に対して垂直に光が照射される場合、すなわち入射角φ=0とすると、P偏光とS偏光との区別はなくなり、各層の界面における振幅反射率および薄膜の位相角βは以下のようになる。
Here, d i represents the thickness of the i layer, and λ represents the wavelength of incident light.
For further simplification, when light is irradiated perpendicularly to the multilayer film sample, that is, when the incident angle φ i = 0, there is no distinction between P-polarized light and S-polarized light, and the amplitude reflectance at the interface of each layer The phase angle β 1 of the thin film is as follows.

Figure 2009092454
Figure 2009092454

さらに、図2に示す3層系の多層膜試料における反射率Rは、以下のようになる。   Further, the reflectance R in the three-layer multilayer film sample shown in FIG. 2 is as follows.

Figure 2009092454
Figure 2009092454

上式において、位相角βについてのフーリエ変換を考えると、位相因子(Phase Factor)であるcos2βは反射率Rに対して非線形となる。そこで、この位相因子cos2βについて線形性を有する関数への変換を行う。一例として、この反射率Rを以下の式のように変換し、独自の変数である「波数変換反射率」R’を定義する。 In the above equation, when considering the Fourier transform for the phase angle β 1 , the phase factor cos 2β 1 is non-linear with respect to the reflectance R. Therefore, to convert to a function having a linearity for the phase factor cos2β 1. As an example, the reflectance R is converted as shown in the following equation, and a unique variable “wave number conversion reflectance” R ′ is defined.

Figure 2009092454
Figure 2009092454

この波数変換反射率R’は、位相因子cos2βについての1次式となり、線形性を有することになる。ここで、式中のRは波数変換反射率R’における切片であり、Rは波数変換反射率R’における傾きである。すなわち、この波数変換反射率R’は、各波長における反射率Rの値を周波数変換に係る位相因子cos2βに対して線形化するための関数である。なお、このような位相因子について線形化するための関数としては、1/(1−R)という関数を用いてもよい。 This wave number conversion reflectance R ′ is a linear expression for the phase factor cos 2β 1 and has linearity. Here, R a in the equation is an intercept in the wave number conversion reflectance R ′, and R b is an inclination in the wave number conversion reflectance R ′. That is, the wave number conversion reflectance R ′ is a function for linearizing the value of the reflectance R at each wavelength with respect to the phase factor cos 2β 1 related to frequency conversion. As a function for linearizing such a phase factor, a function 1 / (1-R) may be used.

したがって、対象とする薄膜における波数Kは以下のように定義できる。 Therefore, the wave number K 1 in the thin film of interest can be defined as follows.

Figure 2009092454
Figure 2009092454

ここで、薄膜中の波長λの光速度をsとし、真空中の波長λの光速度をcとすると、屈折率n=c/sで表される。また、薄膜中をx方向に進行する光によって生じる電磁波E(x,t)は、波数K,角周波数ω,位相δを用いて、E(x,t)=Eexp[j(ωt−Kx+δ)]と表される。すなわち、薄膜中の電磁波の伝搬特性は波数Kに依存する。これらの関係から、真空中において波長λをもつ光は、薄膜中ではその光速度が低下するため、波長もλからλ/nまで長くなることがわかる。このような波長分散現象を考慮して、波数変換反射率R’を以下のように定義する。 Here, when the light velocity of the wavelength λ in the thin film is s and the light velocity of the wavelength λ in the vacuum is c, the refractive index is represented by n 1 = c / s. An electromagnetic wave E (x, t) generated by light traveling in the x direction in the thin film uses E (x, t) = E 0 exp [j (ωt) using the wave number K 1 , the angular frequency ω, and the phase δ. −K 1 x + δ)]. That is, the propagation characteristic of the electromagnetic wave in the thin film depends on the wave number K. From these relationships, it can be seen that light having a wavelength λ in a vacuum has a longer wavelength from λ to λ / n 1 because the speed of light decreases in the thin film. In consideration of such a wavelength dispersion phenomenon, the wave number conversion reflectance R ′ is defined as follows.

Figure 2009092454
Figure 2009092454

この関係から、波数変換反射率R’を波数Kについてフーリエ変換すると、膜厚dに相当する周期成分にピークが現れることにより、このピーク位置を特定することで、膜厚dを算出することができる。 From this relationship, when the wave number conversion reflectance R 'Fourier transform wave number K, by the peaks appearing in the period component corresponding to the film thickness d 1, by identifying the peak position and calculates the film thickness d 1 be able to.

すなわち、本実施の形態に従う多層膜解析装置100は、測定される反射率スペクトルにおける各波長とその波長における反射率の値との対応関係を、各波長から算出される波数と上述の関係式に従って算出される波数変換反射率R’との対応関係に変換する。そして、波数Kを含む波数変換反射率R’の関数に対して、波数Kについて周波数変換し、この周波数変換後の特性に現れるピークに基づいて、対象となる試料における膜厚を算出する。   That is, multilayer analysis apparatus 100 according to the present embodiment determines the correspondence between each wavelength in the measured reflectance spectrum and the reflectance value at that wavelength according to the wave number calculated from each wavelength and the above-described relational expression. Conversion is made to a correspondence relationship with the calculated wave number conversion reflectance R ′. Then, the wave number K is frequency-converted with respect to the function of the wave number conversion reflectance R ′ including the wave number K, and the film thickness of the target sample is calculated based on the peak appearing in the characteristics after the frequency conversion.

ここで、RおよびRは、薄膜内における干渉現象とは無関係な値ではあるが、薄膜の屈折率nを含む各層の界面における振幅反射率に依存する。そのため、屈折率が波長分散をもつ場合には、その値は波長(すなわち、波数K)に依存する関数値となり、波数Kに関して一定値とはならない。そこで、フーリエ変換を⊃で表し、R’,R,R,cos2Kを波数Kでフーリエ変換した後の関数であるパワースペクトルをそれぞれP,P,P,Fとすると、以下の式が成立する。 Here, R a and R b are values irrelevant to the interference phenomenon in the thin film, but depend on the amplitude reflectivity at the interface of each layer including the refractive index n 1 of the thin film. Therefore, when the refractive index has chromatic dispersion, the value is a function value depending on the wavelength (that is, the wave number K), and does not become a constant value with respect to the wave number K. Therefore, when Fourier transform is represented by ⊃, and R ′, R a , R b , and cos 2K 1 d 1 are subjected to Fourier transform with a wave number K, the power spectra as functions are respectively represented by P, P a , P b , and F, respectively. The following equation holds.

Figure 2009092454
Figure 2009092454

式中のPにおける膜厚に依存する成分は相対的に小さく、かつパワースペクトルFとは独立のピークをもつので、パワースペクトルFに影響を与えない。 Component dependent on the film thickness at P a in the formula is relatively small, and since the power spectrum F with separated peaks, does not affect the power spectrum F.

一方、式中のPは、パワースペクトルFとコンボリューションされることにより、Pにおける膜厚成分がパワースペクトルFの膜厚成分に変調を加えることになる。しかしながら、Pは、薄膜内における干渉現象に無関係であり、隣接する2つの層における屈折率の波長依存性のみに影響を受けるため、波数Kに対するPの膜厚成分はFの膜厚成分に比較して無視できる程度に小さい。たとえば、Rが膜厚qの周期関数であるとし、そのフーリエ変換後のPがコンボリューションによりパワースペクトルFの膜厚成分dに変調を加えたとすると、スペクトルとして現れるピークは、「d−q」または「d+q」となるが、qの値が非常に小さいのでピーク位置dに対する影響は小さい。 On the other hand, P b in the equation is convolved with the power spectrum F, so that the film thickness component in P b modulates the film thickness component of the power spectrum F. However, P b is independent of the interference phenomena in the film, since it is affected only in the wavelength dependence of the refractive index in the two adjacent layers, the thickness component of P b for the wavenumber K is the thickness component of the F Is negligibly small compared to For example, assuming that R b is a periodic function of the film thickness q and P b after the Fourier transform is modulated on the film thickness component d of the power spectrum F by convolution, the peak appearing as the spectrum is “d− q ”or“ d + q ”, but since the value of q is very small, the influence on the peak position d is small.

さらに、フーリエ変換を行なう際には、測定対象の薄膜の最大膜厚を考慮して、ナイキストのサンプリング定理に従って、波数変換反射率R’に対して適切なサンプル間隔およびサンプル数でサンプリングが行なわれる。このようにサンプリングされた波数変換反射率R’に基づいて算出されたパワースペクトルの膜厚分解能rに対して、Pの膜厚成分qはより小さい可能性が高く(q<r)、膜厚dの測定結果にはほとんど影響を与えないといえる。 Further, when performing Fourier transform, sampling is performed at an appropriate sample interval and number of samples with respect to the wave number conversion reflectivity R ′ according to the Nyquist sampling theorem in consideration of the maximum film thickness of the thin film to be measured. . Thus for a film thickness resolution r of the power spectrum calculated based on the sampled-wavenumber conversion reflectance R ', thickness component q of P b is high is less than possibility (q <r), film It can be said that the measurement result of the thickness d is hardly affected.

このように、算出された反射率スペクトルを、薄膜における波長分散を考慮した波数についての関数に変換した上で、フーリエ変換を行なうことにより、薄膜の膜厚を正確に算出することができる。   Thus, the film thickness of the thin film can be accurately calculated by performing Fourier transform after converting the calculated reflectance spectrum into a function with respect to the wave number in consideration of wavelength dispersion in the thin film.

なお、上述の説明では、反射率スペクトルを用いる場合について例示したが、透過率スペクトルを用いてもよい。この場合には、測定された透過率をT、「波数変換透過率」をT’とすると、以下のような関係式で表される。   In the above description, the case where the reflectance spectrum is used is exemplified, but the transmittance spectrum may be used. In this case, when the measured transmittance is T, and the “wave number conversion transmittance” is T ′, the following relational expression is obtained.

Figure 2009092454
Figure 2009092454

透過率スペクトルを用いる場合においても、透過率Tは位相因子cos2βに対して非線形となる。そのため、上述したのと同様の理由から、位相因子cos2βについて線形性を有する波数変換透過率T’を採用する。上式によれば、波数変換透過率T’は、位相因子cos2βについての1次式となり、上述したのと同様の手順に従って、薄膜の膜厚を正確に算出することができる。すなわち、この波数変換透過率T’は、各波長における透過率Tの値を周波数変換に係る位相因子cos2βに対して線形化するための関数である。 In the case of using the transmission spectra, the transmittance T becomes non-linear with respect to the phase factor cos2β 1. Therefore, for the same reason as described above, the wave number conversion transmittance T ′ having linearity with respect to the phase factor cos 2β 1 is adopted. According to the above equation, the wave number conversion transmittance T ′ is a linear equation for the phase factor cos 2β 1 , and the film thickness of the thin film can be accurately calculated according to the same procedure as described above. That is, the wave number conversion transmittance T ′ is a function for linearizing the value of the transmittance T at each wavelength with respect to the phase factor cos 2β 1 related to frequency conversion.

(測定例1)
図3は、測定例1の試料から測定された反射率スペクトルの一例を示す図である。図3に示す試料は、石英基板上に5.0μmの窒化シリコンの膜を成膜したものである。
(Measurement Example 1)
FIG. 3 is a diagram illustrating an example of a reflectance spectrum measured from the sample of Measurement Example 1. The sample shown in FIG. 3 is obtained by forming a 5.0 μm silicon nitride film on a quartz substrate.

図4は、図3に示す反射率スペクトルから従来の方法を用いて膜厚を算出した場合の結果を示す図である。   FIG. 4 is a diagram showing the results when the film thickness is calculated from the reflectance spectrum shown in FIG. 3 using a conventional method.

図5は、図3に示す反射率スペクトルから本実施の形態に従う方法を用いて膜厚を算出した場合の結果を示す図である。   FIG. 5 is a diagram showing the results when the film thickness is calculated from the reflectance spectrum shown in FIG. 3 using the method according to the present embodiment.

図4を参照して、従来の方法では、図3に示す反射率スペクトルに対して、ゼロフィリング処理やフィルタ処理などの前処理を行なった後、フーリエ変換によってパワースペクトルを算出した。なお、本測定例1では、高速フーリエ変換(FFT:Fast Fourier Transform)を用いてフーリエ変換を行った。このように反射率スペクトルをそのままフーリエ変換して得られるパワースペクトルは、光学距離(optical distance)についての関数となる。そして、このパワースペクトル中のピーク位置に相当する光学距離を読み取り、この光学距離を、測定対象の試料である窒化シリコンのみかけの屈折率(たとえば、測定波長範囲における平均屈折率)で除することで、測定対象の薄膜の膜厚を算出した。図4に示すように、従来の方法を用いた場合には、本来5.0μmとなるべき値が「5.237μm」として比較的多くの誤差を伴って算出されていることがわかる。   Referring to FIG. 4, in the conventional method, the power spectrum is calculated by Fourier transform after pre-processing such as zero filling processing and filter processing is performed on the reflectance spectrum shown in FIG. In this measurement example 1, Fourier transform was performed using Fast Fourier Transform (FFT). Thus, the power spectrum obtained by directly Fourier-transforming the reflectance spectrum is a function with respect to the optical distance. Then, the optical distance corresponding to the peak position in the power spectrum is read, and this optical distance is divided by the apparent refractive index (for example, the average refractive index in the measurement wavelength range) of the silicon nitride that is the sample to be measured. Thus, the film thickness of the thin film to be measured was calculated. As shown in FIG. 4, when the conventional method is used, the value that should be 5.0 μm is calculated as “5.237 μm” with relatively many errors.

これに対して、図5を参照して、本実施の形態に従う方法では、図3に示す反射率スペクトルを波数変換反射率スペクトルに変換した後、従来の方法と同様に、フーリエ変換によってパワースペクトルを算出した。この波数変換反射率スペクトルをフーリエ変換して得られるパワースペクトルは、薄膜の幾何学的膜厚についての関数となる。そのため、図5に示すパワースペクトルに現れるピーク位置を特定することで、薄膜の膜厚を直接求めることができる。図5に示すように、本実施の形態に従う方法を用いた場合には、薄膜の本来の膜厚である「5.0μm」が算出されていることがわかる。   On the other hand, referring to FIG. 5, in the method according to the present embodiment, the reflectance spectrum shown in FIG. 3 is converted into a wavenumber-converted reflectance spectrum, and then the power spectrum is obtained by Fourier transform as in the conventional method. Was calculated. The power spectrum obtained by Fourier transforming this wave number conversion reflectance spectrum is a function of the geometric thickness of the thin film. Therefore, the film thickness of the thin film can be directly obtained by specifying the peak position appearing in the power spectrum shown in FIG. As shown in FIG. 5, when the method according to the present embodiment is used, it is understood that “5.0 μm”, which is the original film thickness of the thin film, is calculated.

このように、図4と図5とを比較すると、本実施の形態に従う方法によれば、従来法に比較して、より高い精度で膜厚を測定できることがわかる。   Thus, comparing FIG. 4 with FIG. 5, it can be seen that the film thickness can be measured with higher accuracy than the conventional method according to the method according to the present embodiment.

(膜厚測定処理その2)
上述の説明では、測定対象の薄膜が1層である場合について例示したが、複数の薄膜が存在している場合にも同様の手順でそれぞれの膜厚を測定することができる。
(Thickness measurement process 2)
In the above description, the case where the thin film to be measured is a single layer has been exemplified, but each film thickness can be measured in the same procedure even when a plurality of thin films are present.

一例として、図6に示すような2層の薄膜層が形成された試料を考える。第1層の屈折率をn、膜厚をdとし、第2層の屈折率をn21、膜厚をdとすると、波数変換反射率R’は以下のように表される。 As an example, consider a sample in which two thin film layers as shown in FIG. 6 are formed. When the refractive index of the first layer is n 1 , the film thickness is d 1 , the refractive index of the second layer is n 21 , and the film thickness is d 2 , the wavenumber conversion reflectance R ′ is expressed as follows.

Figure 2009092454
Figure 2009092454

ここで、第1層および第2層の膜厚dおよびdを算出する場合には、波数K,Kでそれぞれ変換した波数変換反射率R’(K),R’(K)を用いる。具体的には、以下のように表される。 Here, when calculating the film thicknesses d 1 and d 2 of the first layer and the second layer, wave number conversion reflectances R 1 ′ (K 1 ) and R 2 ′ converted by the wave numbers K 1 and K 2 , respectively. (K 2 ) is used. Specifically, it is expressed as follows.

Figure 2009092454
Figure 2009092454

これらの式中において、d’およびd’は正しい膜厚値ではないが、波数変換反射率R’(K)の第2項に相当するパワースペクトル中のピークから本来の膜厚dが求まり、波数変換反射率R’(K)の第3項に相当するパワースペクトル中のピークから本来の膜厚dが求まる。 In these equations, d 1 ′ and d 2 ′ are not correct film thickness values, but the original film thickness from the peak in the power spectrum corresponding to the second term of the wave number conversion reflectance R 1 ′ (K 1 ). d 1 is obtained, and the original film thickness d 2 is obtained from the peak in the power spectrum corresponding to the third term of the wave number conversion reflectance R 2 ′ (K 2 ).

なお、実際には、第1層および第2層は、その屈折率が近似していることが多く、両者の界面における反射率は、他の界面における反射率に比較して相対的に小さくなることが多い。その結果、波数変換反射率の関数に含まれるRやRに比較して、Rの値が小さくなり、パワースペクトルから、波数変換反射率R’(K)の第3項に相当するピークを識別することが困難である場合も多い。このような場合には、波数変換反射率R’(K)の第4項に相当するパワースペクトルのピーク位置(d’+d)と、波数変換反射率R’(K)の第2項に相当するパワースペクトルのピーク位置(d’)とを算出した上で、両者の差をとることで、膜厚dを算出することができる。 Actually, the refractive index of the first layer and the second layer is often approximate, and the reflectance at the interface between the two is relatively smaller than the reflectance at the other interface. There are many cases. As a result, the value of R c becomes smaller than R b and R d included in the function of wave number conversion reflectivity, and the third term of the wave number conversion reflectivity R 2 ′ (K 2 ) is obtained from the power spectrum. Often it is difficult to identify the corresponding peak. In such a case, 'the peak position of the power spectrum corresponding to the fourth term of the (K 2) (d 1' wavenumber conversion reflectance R 2 and + d 2), wavenumber conversion reflectance R 2 '(K 2) After calculating the peak position (d 1 ′) of the power spectrum corresponding to the second term, the film thickness d 2 can be calculated by taking the difference between the two .

(測定例2)
図7は、測定例2の試料から測定された反射率スペクトルの一例を示す図である。図7に示す試料は、基板上に、3.0μmの薄膜(第1層)と、2.0μmの薄膜(第2層)とを成膜したものである。なお、図7には、各薄膜および基板の透過率についても示す。
(Measurement example 2)
FIG. 7 is a diagram illustrating an example of the reflectance spectrum measured from the sample of Measurement Example 2. The sample shown in FIG. 7 is obtained by depositing a 3.0 μm thin film (first layer) and a 2.0 μm thin film (second layer) on a substrate. FIG. 7 also shows the transmittance of each thin film and substrate.

図8は、図7に示す反射率スペクトルから変換した第1層に係る波数変換反射率R’(K)についてのパワースペクトルである。図9は、図7に示す反射率スペクトルから変換した第2層に係る波数変換反射率R’(K)についてのパワースペクトルである。 FIG. 8 is a power spectrum for the wave number conversion reflectance R 1 ′ (K 1 ) of the first layer converted from the reflectance spectrum shown in FIG. FIG. 9 is a power spectrum for the wave number conversion reflectance R 2 ′ (K 2 ) according to the second layer converted from the reflectance spectrum shown in FIG.

図8を参照して、本来ならば波数変換反射率R’(K)の各項に対応する4つのピークが観測されるはずであるが、上述したように、薄膜間の反射率が相対的に低いことや、フーリエ変換処理におけるノイズカットのしきい値などを考慮すると、主として、波数変換反射率R’(K)の第2項および第4項に由来する2つのピークが現れているものと考えられる。この2つのピークのうち、より小さい膜厚に位置するピークに対応する膜厚を第1層の膜厚dとして求めることができる。すなわち、図8に示すように、膜厚が3μmの位置にピークが存在しており、このことから、第1層の膜厚d=3μmと求めることができる。 Referring to FIG. 8, four peaks corresponding to each term of wave number conversion reflectance R 1 ′ (K 1 ) should be observed, but as described above, the reflectance between thin films is Considering the relatively low value and the noise cut threshold value in the Fourier transform process, two peaks derived mainly from the second and fourth terms of the wave number transform reflectivity R 1 ′ (K 1 ) It is thought that it is appearing. Of the two peaks, it is possible to obtain a film thickness corresponding to the peak located at a smaller thickness as the thickness d 1 of the first layer. That is, as shown in FIG. 8, a peak exists at a position where the film thickness is 3 μm, and from this, it can be determined that the film thickness d 1 of the first layer is 3 μm.

次に、図9を参照して、波数変換反射率R’(K)についてのパワースペクトルにおいても、図8と同様に、主として波数変換反射率R’(K)の第2項および第4項に由来する2つのピークのみが現れている。ここで、波数変換反射率R’(K)についてのパワースペクトルを、図8に示すパワースペクトルから算出される第1層の膜厚に相当するd’だけシフトさせた上で、波数変換反射率R’(K)の第4項に由来するピークに対応する膜厚を第2層の膜厚dとして求めることができる。すなわち、図9に示すように、シフト後の第4項に由来するピークが2μmの位置に存在しており、このことから、第2層の膜厚d=2μmと求めることができる。 Next, with reference to FIG. 9, 'even in the power spectra for (K 2), similarly to FIG. 8, mainly wavenumber conversion reflectance R 2' wavenumber conversion reflectance R 2 the second term of (K 2) And only two peaks from the fourth term appear. Here, after shifting the power spectrum of the wave number conversion reflectance R 2 ′ (K 2 ) by d 1 ′ corresponding to the film thickness of the first layer calculated from the power spectrum shown in FIG. The film thickness corresponding to the peak derived from the fourth term of the conversion reflectance R 2 ′ (K 2 ) can be obtained as the film thickness d 2 of the second layer. That is, as shown in FIG. 9, the peak derived from the fourth term after the shift exists at a position of 2 μm, and from this, the film thickness d 2 of the second layer can be obtained as 2 μm.

上述と同様の手順によって、3層以上の多層膜試料についても、波数変換反射率または波数変換透過率を用いて、各層の膜厚を屈折率の波長依存性の影響を受けることなく算出することができる。
(データ処理部の構成)
図10は、この発明の実施の形態に従うデータ処理部70の概略のハードウェア構成を示す模式図である。
Calculate the film thickness of each layer without being affected by the wavelength dependence of the refractive index using the wave number conversion reflectance or wave number conversion transmittance for a multilayer film sample of three or more layers by the same procedure as described above. Can do.
(Configuration of data processing unit)
FIG. 10 is a schematic diagram showing a schematic hardware configuration of data processing unit 70 according to the embodiment of the present invention.

図10を参照して、データ処理部70は、代表的にコンピュータによって実現され、オペレーティングシステム(OS:Operating System)を含む各種プログラムを実行するCPU(Central Processing Unit)200と、CPU200でのプログラムの実行に必要なデータを一時的に記憶するメモリ部212と、CPU200で実行されるプログラムを不揮発的に記憶するハードディスク部(HDD:Hard Disk Drive)210とを含む。また、ハードディスク部210には、後述するような処理を実現するためのプログラムが予め記憶されており、このようなプログラムは、FDDドライブ216またはCD−ROMドライブ214によって、それぞれフレキシブルディスク216aまたはCD−ROM(Compact Disk-Read Only Memory)214aなどから読み取られる。   Referring to FIG. 10, data processing unit 70 is typically realized by a computer, and includes a CPU (Central Processing Unit) 200 that executes various programs including an operating system (OS), and a program executed by CPU 200. A memory unit 212 that temporarily stores data necessary for execution and a hard disk unit (HDD: Hard Disk Drive) 210 that stores a program executed by the CPU 200 in a nonvolatile manner are included. The hard disk unit 210 stores a program for realizing processing as will be described later, and such a program is stored in the flexible disk 216a or the CD-ROM by the FDD drive 216 or the CD-ROM drive 214, respectively. It is read from a ROM (Compact Disk-Read Only Memory) 214a or the like.

CPU200は、キーボードやマウスなどからなる入力部208を介してユーザなどからの指示を受取るとともに、プログラムの実行によって測定される測定結果などをディスプレイ部204へ出力する。   The CPU 200 receives an instruction from a user or the like via the input unit 208 including a keyboard and a mouse, and outputs a measurement result measured by executing the program to the display unit 204.

(制御構造)
図11は、この発明の実施の形態に従う膜厚測定処理を実現するための制御構造を示すブロック図である。図11に示すブロック図は、CPU200がハードディスク部210などの予め格納されたプログラムをメモリ部212などに読み出して実行することで実現される。
(Control structure)
FIG. 11 is a block diagram showing a control structure for realizing a film thickness measurement process according to the embodiment of the present invention. The block diagram shown in FIG. 11 is realized by the CPU 200 reading a program stored in advance such as the hard disk unit 210 into the memory unit 212 and executing the program.

図11を参照して、データ処理部70(図1)は、バッファ部71と、波数変換部721,722,…,72mと、バッファ部731,732,…,73mと、フーリエ変換部741,742,…,74mと、ピーク探索部751,752,…,75mと、総合算出部76とをその機能として含む。なお、図10には、m層の薄膜が形成された試料において、各層の膜厚を測定する汎用的な構成を示す。   11, data processing unit 70 (FIG. 1) includes buffer unit 71, wave number conversion units 721, 722,..., 72m, buffer units 731, 732,. , 74m, peak search units 751, 752,..., 75m, and a total calculation unit 76 are included as functions. FIG. 10 shows a general-purpose configuration for measuring the thickness of each layer in a sample in which m layers of thin films are formed.

バッファ部71は、分光測定部60から出力される反射率スペクトルR(λ)を一時的に格納する。より具体的には、分光測定部60からは所定の波長分解能毎に反射率の値が出力されるので、バッファ部71は、波長とその波長における反射率とを対応付けて格納する。   The buffer unit 71 temporarily stores the reflectance spectrum R (λ) output from the spectroscopic measurement unit 60. More specifically, since the reflectance value is output from the spectroscopic measurement unit 60 for each predetermined wavelength resolution, the buffer unit 71 stores the wavelength and the reflectance at that wavelength in association with each other.

波数変換部721は、ユーザなどによって予め入力される薄膜のうち第1層の屈折率n1に基づいて、バッファ部71に一時的に格納される反射率スペクトルを波数変換する。すなわち、波数変換部721は、反射率スペクトルにおける各波長とその波長における反射率との対応関係を、各波長についての波数Kと上述の関係式に従って算出される対応の波数変換反射率R’との対応関係に変換する。より具体的には、波数変換部721は、バッファ部71に格納される波長毎に、波数K(λ)および波数変換反射率R’(λ)(=R(λ)/(1−R(λ)))を順次算出し、バッファ部731へ出力する。 The wave number conversion unit 721 performs wave number conversion on the reflectance spectrum temporarily stored in the buffer unit 71 based on the refractive index n1 of the first layer among thin films input in advance by a user or the like. That is, the wave number conversion unit 721 determines the correspondence between each wavelength in the reflectance spectrum and the reflectance at that wavelength, and the corresponding wave number conversion reflectance R 1 calculated according to the wave number K 1 for each wavelength and the above relational expression. Convert to a correspondence with '. More specifically, the wave number conversion unit 721 generates a wave number K 1 (λ) and a wave number conversion reflectance R 1 ′ (λ) (= R (λ) / (1−) for each wavelength stored in the buffer unit 71. R (λ))) are sequentially calculated and output to the buffer unit 731.

バッファ部731は、波数変換部721から順次出力される波数K(λ)と波数変換反射率R’(λ)とを対応付けて格納する。すなわち、バッファ部731には、波数K(λ)に関する波数変換反射率の波数分布特性である波数変換反射率R’(K)が格納される。 The buffer unit 731 stores the wave number K 1 (λ) and the wave number conversion reflectance R 1 ′ (λ) sequentially output from the wave number conversion unit 721 in association with each other. That is, the buffer unit 731 stores the wave number conversion reflectance R 1 ′ (K 1 ) that is the wave number distribution characteristic of the wave number conversion reflectance with respect to the wave number K 1 (λ).

フーリエ変換部741は、バッファ部731に格納される波数変換反射率R’(K)を波数Kについてフーリエ変換を行って、パワースペクトルPを算出する。なお、フーリエ変換の方法としては、高速フーリエ変換(FFT)や離散コサイン変換(DCT:Discrete Cosine Transform)などを用いることができる。 The Fourier transform unit 741 calculates the power spectrum P 1 by performing a Fourier transform on the wave number conversion reflectance R 1 ′ (K 1 ) stored in the buffer unit 731 with respect to the wave number K 1 . In addition, as a method of Fourier transform, fast Fourier transform (FFT), discrete cosine transform (DCT: Discrete Cosine Transform), etc. can be used.

ピーク探索部751は、フーリエ変換部741によって算出されたパワースペクトルPの中に現れるピークを探索し、各ピーク位置に対応する膜厚を取得する。 Peak search unit 751 searches for a peak appearing in the power spectrum P 1 calculated by Fourier transform unit 741, it acquires a thickness corresponding to the peak position.

上述のように、波数変換部721と、バッファ部731と、フーリエ変換部741と、ピーク探索部751とは、第1層の薄膜に起因するパワースペクトルに現れるピークに対応する膜厚を取得する。   As described above, the wave number conversion unit 721, the buffer unit 731, the Fourier transform unit 741, and the peak search unit 751 acquire the film thickness corresponding to the peak appearing in the power spectrum caused by the thin film of the first layer. .

同様にして、波数変換部722と、バッファ部732と、フーリエ変換部742と、ピーク探索部752とは、第2層の薄膜に起因するパワースペクトルに現れるピークに対応する膜厚を取得する。さらに、波数変換部72mと、バッファ部73mと、フーリエ変換部74mと、ピーク探索部75mとは、第m層の薄膜に起因するパワースペクトルに現れるピークに対応する膜厚を取得する。   Similarly, the wave number conversion unit 722, the buffer unit 732, the Fourier transform unit 742, and the peak search unit 752 acquire the film thickness corresponding to the peak appearing in the power spectrum caused by the thin film of the second layer. Furthermore, the wave number converting unit 72m, the buffer unit 73m, the Fourier transforming unit 74m, and the peak searching unit 75m acquire the film thickness corresponding to the peak appearing in the power spectrum caused by the m-th layer thin film.

総合算出部76は、ピーク探索部751,752,…,75mから出力される各パワースペクトルに現れるピークに対応する膜厚に基づいて、各薄膜における膜厚を算出する。より具体的には、図8および図9に示すように、まず、ピーク探索部751から出力されるピークに対応する膜厚に基づいて第1層の膜厚が決定され、この第1層の膜厚とピーク探索部752から出力されるピークに対応する膜厚とに基づいて第2層の膜厚が決定される。以下同様にして、各層の膜厚が順次算出される。そして、総合算出部76は、算出した各層の膜厚の測定結果を、ユーザや図示しない上位のコンピュータなどへ出力する。   The total calculation unit 76 calculates the film thickness of each thin film based on the film thickness corresponding to the peak appearing in each power spectrum output from the peak search units 751, 752,. More specifically, as shown in FIG. 8 and FIG. 9, first, the film thickness of the first layer is determined based on the film thickness corresponding to the peak output from the peak search unit 751. The film thickness of the second layer is determined based on the film thickness and the film thickness corresponding to the peak output from the peak search unit 752. In the same manner, the thickness of each layer is sequentially calculated. Then, the total calculation unit 76 outputs the calculated measurement results of the thicknesses of the respective layers to the user, a host computer (not shown), or the like.

図11に示す制御構造と本願発明との対応関係については、波数変換部721,722,…,72mが「波数特性変換手段」に対応し、フーリエ変換部741,742,…,74mが「周波数変換手段」に対応し、総合算出部76が「膜厚算出手段」に対応する。   11, the wave number conversion units 721, 722,..., 72m correspond to “wave number characteristic conversion means”, and the Fourier conversion units 741, 742,. The total calculation unit 76 corresponds to “conversion means” and “film thickness calculation means”.

(処理手順)
図12は、この発明の実施の形態に従う膜厚測定に係る処理手順を示すフローチャートである。
(Processing procedure)
FIG. 12 is a flowchart showing a processing procedure related to film thickness measurement according to the embodiment of the present invention.

図12を参照して、まず、ユーザは、測定対象の試料をステージ50(図1)上に配置する(ステップS100)。また、ユーザは、測定対象の試料を構成する薄膜の屈折率をデータ処理部70へ入力する(ステップS102)。ユーザによる薄膜の屈折率の入力に応じて、データ処理部70のCPU200は、これらの値をメモリ部212などに格納する。   Referring to FIG. 12, first, the user places a sample to be measured on stage 50 (FIG. 1) (step S100). In addition, the user inputs the refractive index of the thin film constituting the sample to be measured to the data processing unit 70 (step S102). In response to the input of the refractive index of the thin film by the user, the CPU 200 of the data processing unit 70 stores these values in the memory unit 212 or the like.

ここで、観察用光源22からは観察光の照射が開始され、ユーザは、表示部39に表示される観察用カメラ38で撮影された反射像を参照しながら、可動機構52にステージ位置指令を与えて、測定対象の領域の調整や焦点合わせを行なう(ステップS104)。   Here, irradiation of observation light is started from the observation light source 22, and the user gives a stage position command to the movable mechanism 52 while referring to the reflection image taken by the observation camera 38 displayed on the display unit 39. Then, adjustment and focus adjustment of the area to be measured are performed (step S104).

その後、ユーザが測定開始指令を与えると、測定用光源10から測定光の発生が開始され、分光測定部60は、試料からの反射光から取得される反射率スペクトルをデータ処理部70へ出力する(ステップS106)。   Thereafter, when the user gives a measurement start command, generation of measurement light from the measurement light source 10 is started, and the spectroscopic measurement unit 60 outputs the reflectance spectrum acquired from the reflected light from the sample to the data processing unit 70. (Step S106).

データ処理部70のCPU200は、分光測定部60からの反射率スペクトルをメモリ部212などに一時的に格納する(ステップS108)。CPU200は、入力された屈折率の1つに基づいて、反射率スペクトルを波数変換する(ステップS110)。そして、CPU200は、この波数変換して算出される波数変換反射率をメモリ部212などに格納する(ステップS112)。さらに、CPU200は、算出された波数変換反射率に対してフーリエ変換を行って、パワースペクトルを算出する(ステップS114)。その後、CPU200は、算出したパワースペクトルに現れるピークおよびそのピークに対応する膜厚を取得する(ステップS116)。   The CPU 200 of the data processing unit 70 temporarily stores the reflectance spectrum from the spectroscopic measurement unit 60 in the memory unit 212 or the like (step S108). CPU 200 performs wave number conversion on the reflectance spectrum based on one of the input refractive indexes (step S110). Then, the CPU 200 stores the wave number conversion reflectance calculated by the wave number conversion in the memory unit 212 or the like (step S112). Further, the CPU 200 performs a Fourier transform on the calculated wave number conversion reflectivity to calculate a power spectrum (step S114). Thereafter, the CPU 200 acquires the peak appearing in the calculated power spectrum and the film thickness corresponding to the peak (step S116).

その後、CPU200は、入力されたすべての屈折率について、ピークおよびそのピークに対応する膜厚の取得処理が実行されたか否かを判断する(ステップS118)。入力されたすべての屈折率について、ピークおよびそのピークに対応する膜厚の取得処理が実行されていなければ(ステップS118においてNO)、残りの反射率の1つについて、ステップS110以下の処理が繰返し実行される。   Thereafter, the CPU 200 determines whether or not the peak and the film thickness acquisition process corresponding to the peak have been executed for all input refractive indexes (step S118). If the processing for obtaining the peak and the film thickness corresponding to the peak has not been executed for all input refractive indexes (NO in step S118), the processing from step S110 onward is repeated for one of the remaining reflectances. Executed.

一方、入力されたすべての屈折率について、ピークおよびそのピークに対応する膜厚の取得処理が実行されていれば(ステップS118においてYES)、CPU200は、各屈折率について算出されるピークおよびそのピークに対応する膜厚に基づいて、各薄膜における膜厚を算出する(ステップS120)。そして、CPU200は、算出した各薄膜における膜厚をディスプレイ部204やインターフェイス部206を介して図示しない上位のコンピュータなどへ出力する(ステップS122)。   On the other hand, if the processing for obtaining the peak and the film thickness corresponding to the peak has been executed for all input refractive indexes (YES in step S118), CPU 200 calculates the peak calculated for each refractive index and the peak. Based on the film thickness corresponding to, the film thickness in each thin film is calculated (step S120). Then, the CPU 200 outputs the calculated film thickness of each thin film to a host computer (not shown) or the like via the display unit 204 or the interface unit 206 (step S122).

以上のような一連の処理に従って、多層膜試料の各層の膜厚が測定できる。
(層構造同定処理)
上述したように、反射率スペクトルをフーリエ変換して得られるパワースペクトルには、多層膜試料の層構造に応じたピークが現れる。しかしながら、算出されるパワースペクトルは、その算出過程において、各種のフィルタ処理・打ち切り処理・アポダイゼーション(apodization)処理といった前処理や、ベースラインの設定に伴うDCオフセットピークといった擬似ピークの発生の影響を受ける。また、層間の間にも光学的な現象が生じ得る。そのため、パワースペクトルには、1つの層から1つのピークが現れるわけではなく、層の組合せに応じた複数のピークが出現する。そのため、層構造が予め既知の場合であっても、パワースペクトルに現れるピークがいずれの層の干渉に起因するものであるかをスペクトルの形状だけで判断することは困難な場合が多い。ましてや層構造が未知の場合には、スペクトルの形状から層構造を特定することは困難である。
The film thickness of each layer of the multilayer film sample can be measured according to the series of processes as described above.
(Layer structure identification process)
As described above, a peak corresponding to the layer structure of the multilayer film sample appears in the power spectrum obtained by Fourier transforming the reflectance spectrum. However, in the calculation process, the calculated power spectrum is affected by the generation of pseudo-peaks such as pre-processing such as various filter processing, truncation processing, and apodization processing, and DC offset peak associated with the baseline setting. . Also, an optical phenomenon can occur between the layers. Therefore, in the power spectrum, one peak does not appear from one layer, but a plurality of peaks appear according to the combination of layers. For this reason, even when the layer structure is known in advance, it is often difficult to determine from which spectrum the peak that appears in the power spectrum is due to the interference of which layer. If the layer structure is unknown, it is difficult to identify the layer structure from the shape of the spectrum.

そこで、本実施の形態に従う多層膜解析装置は、対象とする試料の層構造を同定する機能を備える。この層構造を同定する機能は、代表的にデータ処理部70に実装される。具体的には、予め規定された複数の層構造のうち所定数を抽出することで仮想の層構造を設定して、この仮想の層構造から反射率スペクトル(理論値)を計算し、さらにこの反射率スペクトルをフーリエ変換してパワースペクトルを算出する。同時に、実測値の反射率スペクトルから同じ処理手順でフーリエ変換してパワースペクトルを算出する。そして、これらのパワースペクトルの間で相関係数を算出し、この算出された相関係数に基づいて、仮想の層構造を順次変更していき、最も相関係数が高くなった組合せをその試料における層構造として決定する。   Therefore, the multilayer analysis apparatus according to the present embodiment has a function of identifying the layer structure of the target sample. The function of identifying this layer structure is typically implemented in the data processing unit 70. Specifically, a virtual layer structure is set by extracting a predetermined number from a plurality of predefined layer structures, and a reflectance spectrum (theoretical value) is calculated from the virtual layer structure. The power spectrum is calculated by Fourier transforming the reflectance spectrum. At the same time, a power spectrum is calculated from the reflectance spectrum of the actually measured value by Fourier transform in the same processing procedure. Then, a correlation coefficient is calculated between these power spectra, and the virtual layer structure is sequentially changed based on the calculated correlation coefficient, and the combination with the highest correlation coefficient is selected as the sample. Determined as the layer structure.

図13は、この発明の実施の形態に従う層構造同定処理を実現するための制御構造を示すブロック図である。図13に示すブロック図は、CPU200がハードディスク部210などの予め格納されたプログラムをメモリ部212などに読み出して実行することで実現される。   FIG. 13 is a block diagram showing a control structure for realizing the layer structure identification process according to the embodiment of the present invention. The block diagram shown in FIG. 13 is realized by the CPU 200 reading a program stored in advance such as the hard disk unit 210 into the memory unit 212 and executing the program.

図13を参照して、データ処理部70は、層構造データベース300と、組合せ部304と、反射率算出部306と、理論スペクトル算出部308と、パワースペクトル算出部310,320と、バッファ部318と、相関係数算出部330と、層構造決定部332とをその機能として含む。   Referring to FIG. 13, the data processing unit 70 includes a layer structure database 300, a combination unit 304, a reflectance calculation unit 306, a theoretical spectrum calculation unit 308, power spectrum calculation units 310 and 320, and a buffer unit 318. And a correlation coefficient calculation unit 330 and a layer structure determination unit 332 as functions thereof.

層構造データベース300は、複数の層構造情報302を予め格納する。この層構造情報302の各々は、各層について、少なくとも屈折率および膜厚の情報を含む。この層構造情報は、上述の膜厚測定処理によって算出された膜厚や既知の屈折率などに基づいて、ユーザが任意に設定することができる。   The layer structure database 300 stores a plurality of layer structure information 302 in advance. Each of the layer structure information 302 includes at least information of refractive index and film thickness for each layer. The layer structure information can be arbitrarily set by the user based on the film thickness calculated by the above-described film thickness measurement process, the known refractive index, and the like.

組合せ部304は、層構造データベース300に格納される複数の層構造情報302のうち、所定数の層構造情報302を組合せて、候補となる層構造を選択する。なお、組合せ部304が組合せる層の数については、予めユーザが設定するようにしてもよい。   The combination unit 304 selects a candidate layer structure by combining a predetermined number of layer structure information 302 among the plurality of layer structure information 302 stored in the layer structure database 300. Note that the number of layers to be combined by the combining unit 304 may be set in advance by the user.

反射率算出部306は、組合せ部304によって組合せられた候補に対応する層構造情報(膜厚および屈折率)に基づいて、当該層構造における反射率を解析的に算出する。なお、反射率に代えて絶対透過率を算出してもよい。   Based on the layer structure information (film thickness and refractive index) corresponding to the candidates combined by the combination unit 304, the reflectance calculation unit 306 analytically calculates the reflectance in the layer structure. Note that the absolute transmittance may be calculated instead of the reflectance.

この反射率(あるいは、屈折率)の算出については、上述したように、隣接する層界面におけるフレネル係数の計算、および各層の干渉による位相角の計算を順次行ない、薄膜内において無限に反射が繰返されるとして、無限等比級数の収束値(の2乗)として反射率(または、絶対透過率)を解析的に計算できる。   In calculating the reflectivity (or refractive index), as described above, the Fresnel coefficient at the adjacent layer interface and the phase angle due to interference of each layer are sequentially calculated, and reflection is repeated infinitely within the thin film. As a result, the reflectance (or absolute transmittance) can be analytically calculated as the convergence value (square) of an infinite geometric series.

理論スペクトル算出部308は、反射率算出部306によって算出された候補となる層構造の反射率の演算式に基づいて、各波長における強度値を順次算出して、候補となる層構造における反射率スペクトル(理論値)を算出する。なお、この反射率スペクトルの波長分解能は、測定される反射率スペクトルの分解能と同様とすることが好ましい。   The theoretical spectrum calculation unit 308 sequentially calculates the intensity value at each wavelength based on the calculation formula of the reflectance of the candidate layer structure calculated by the reflectance calculation unit 306, and reflects the reflectance in the candidate layer structure. A spectrum (theoretical value) is calculated. The wavelength resolution of the reflectance spectrum is preferably the same as the resolution of the measured reflectance spectrum.

パワースペクトル算出部310は、理論スペクトル算出部308で算出される反射率スペクトル(理論値)をフーリエ変換してパワースペクトル(理論値)を算出する。ここでパワースペクトル算出部310は、図11の波数変換部721と、バッファ部731と、フーリエ変換部741とによって行なわれるパワースペクトルの算出処理と同一の処理を実行する。これは、実測された反射率スペクトルから算出されるパワースペクトルと、理論上で算出されるパワースペクトルとを同一条件で比較するためである。   The power spectrum calculation unit 310 performs Fourier transform on the reflectance spectrum (theoretical value) calculated by the theoretical spectrum calculation unit 308 to calculate a power spectrum (theoretical value). Here, the power spectrum calculation unit 310 executes the same processing as the power spectrum calculation processing performed by the wave number conversion unit 721, the buffer unit 731, and the Fourier transform unit 741 in FIG. 11. This is because the power spectrum calculated from the actually measured reflectance spectrum and the theoretically calculated power spectrum are compared under the same conditions.

一方、バッファ部318は、分光測定部60で実測された反射率スペクトルを一時的に格納する。そして、パワースペクトル算出部320は、この反射率スペクトル(実測値)を上述したパワースペクトル算出部310と同一のパワースペクトル算出処理によって、パワースペクトル(実測値)を算出する。   On the other hand, the buffer unit 318 temporarily stores the reflectance spectrum actually measured by the spectroscopic measurement unit 60. Then, the power spectrum calculation unit 320 calculates the power spectrum (measured value) from the reflectance spectrum (measured value) by the same power spectrum calculation process as that of the power spectrum calculation unit 310 described above.

相関係数算出部330は、パワースペクトル算出部310によって算出されたパワースペクトル(理論値)と、パワースペクトル算出部320によって算出されたパワースペクトル(実測値)との間で相関係数を算出する。このような相関係数の算出の一例としては、パワースペクトル(実測値)をR(K)とし、パワースペクトル(理論値)をR(K)とすると、2つのパワースペクトルの相互相関関数Cmc=<R(K)R(K+κ)>に従って算出することができる。あるいは、Pearsonの線形相関係数を算出してもよい。 The correlation coefficient calculation unit 330 calculates a correlation coefficient between the power spectrum (theoretical value) calculated by the power spectrum calculation unit 310 and the power spectrum (actual measurement value) calculated by the power spectrum calculation unit 320. . As an example of calculation of such a correlation coefficient, when the power spectrum (actual value) is R m (K) and the power spectrum (theoretical value) is R c (K), the cross-correlation function of the two power spectra It can be calculated according to C mc = <R m (K) R c (K + κ)>. Alternatively, Pearson's linear correlation coefficient may be calculated.

層構造決定部332は、相関係数算出部330によって算出された相関係数とそのときの候補となる層構造とを対応付けて格納し、すべての組合せについての相関係数の算出が完了した後に、最も相関係数が高い組合せをその試料における層構造として決定する。   The layer structure determination unit 332 stores the correlation coefficient calculated by the correlation coefficient calculation unit 330 in association with the candidate layer structure at that time, and the calculation of correlation coefficients for all combinations is completed. Later, the combination with the highest correlation coefficient is determined as the layer structure in the sample.

図14は、この発明の実施の形態に従う層構造同定処理に係る処理手順を示すフローチャートである。   FIG. 14 is a flowchart showing a processing procedure related to the layer structure identification processing according to the embodiment of the present invention.

図14を参照して、まず、ユーザは、図12に示す膜厚測定に係る処理手順のステップS100〜S106と同様の手順に従って、試料からの反射光に基づいて反射率スペクトル(実測値)を取得する(ステップS200)。この反射率スペクトルがデータ処理部70へ入力されると、データ処理部70のCPU200は、この反射率スペクトルをメモリ部212などに一時的に格納する(ステップS202)。そして、CPU200は、図12に示す膜厚測定に係る処理手順のステップS110〜S114と同様の手順に従って、反射率スペクトル(実測値)からパワースペクトルを算出する(ステップS204)。   Referring to FIG. 14, first, the user follows a procedure similar to steps S100 to S106 of the processing procedure relating to the film thickness measurement shown in FIG. 12 on the basis of the reflectance spectrum (measured value) based on the reflected light from the sample. Obtain (step S200). When the reflectance spectrum is input to the data processing unit 70, the CPU 200 of the data processing unit 70 temporarily stores the reflectance spectrum in the memory unit 212 or the like (step S202). Then, the CPU 200 calculates a power spectrum from the reflectance spectrum (actual measurement value) according to the same procedure as steps S110 to S114 of the processing procedure relating to the film thickness measurement shown in FIG. 12 (step S204).

一方、CPU200は、ハードディスク部210(層構造データベース300)に格納される複数の層構造情報302のうち、所定数の層構造情報302を組合せて、候補となる層構造を選択する(ステップS206)。そして、CPU200は、候補となる層構造に対応する層構造情報(膜厚および屈折率)に基づいて、選択した層構造における反射率を解析的に算出する(ステップS208)。さらに、CPU200は、候補となる層構造の反射率の演算式に基づいて、選択した層構造における反射率スペクトル(理論値)を算出する(ステップS210)。さらに、CPU200は、図12に示す膜厚測定に係る処理手順のステップS110〜S114と同様の手順に従って、反射率スペクトル(理論値)からパワースペクトルを算出する(ステップS212)。   On the other hand, the CPU 200 selects a candidate layer structure by combining a predetermined number of layer structure information 302 among the plurality of layer structure information 302 stored in the hard disk unit 210 (layer structure database 300) (step S206). . Then, the CPU 200 analytically calculates the reflectance in the selected layer structure based on the layer structure information (film thickness and refractive index) corresponding to the candidate layer structure (step S208). Further, the CPU 200 calculates the reflectance spectrum (theoretical value) in the selected layer structure based on the calculation formula of the reflectance of the candidate layer structure (step S210). Further, the CPU 200 calculates a power spectrum from the reflectance spectrum (theoretical value) according to a procedure similar to steps S110 to S114 of the processing procedure relating to film thickness measurement shown in FIG. 12 (step S212).

その後、CPU200は、パワースペクトル(理論値)と、パワースペクトル(実測値)との間で相関係数を算出し(ステップS214)、この算出した相関係数と選択されている層構造(組合せ)とを対応付けて、メモリ部212などに格納する(ステップS216)。そして、CPU200は、ハードディスク部210(層構造データベース300)に格納される複数の層構造情報302について、未選択の組合せが存在するか否かを判断する(ステップS218)。未選択の組合せが存在する場合(ステップS218においてYES)には、ステップS206以下の処理が繰返される。   Thereafter, the CPU 200 calculates a correlation coefficient between the power spectrum (theoretical value) and the power spectrum (actual measurement value) (step S214), and the calculated correlation coefficient and the selected layer structure (combination). Are associated with each other and stored in the memory unit 212 or the like (step S216). Then, the CPU 200 determines whether or not an unselected combination exists for the plurality of layer structure information 302 stored in the hard disk unit 210 (layer structure database 300) (step S218). If there is an unselected combination (YES in step S218), the processes in and after step S206 are repeated.

これに対して、未選択の組合せが存在しない場合(ステップS218においてNO)には、CPU200は、メモリ部212に格納した相関係数のうち、最も高いものに対応する層構造を対象となる試料における層構造に決定する(ステップS220)。そして、CPU200は、決定した層構造をディスプレイ部204やインターフェイス部206を介して図示しない上位のコンピュータなどへ出力する(ステップS222)。   On the other hand, when there is no unselected combination (NO in step S218), the CPU 200 targets the layer structure corresponding to the highest one of the correlation coefficients stored in the memory unit 212. The layer structure is determined at (step S220). Then, the CPU 200 outputs the determined layer structure to an upper computer (not shown) or the like via the display unit 204 or the interface unit 206 (step S222).

以上のような一連の処理に従って、層構造の同定ができる。
図15は、任意の試料に対して測定された反射率スペクトルを示す測定例である。
The layer structure can be identified according to the series of processes as described above.
FIG. 15 is a measurement example showing a reflectance spectrum measured for an arbitrary sample.

図16は、この発明の実施の形態に従う層構造同定処理を用いて図15に示す反射率スペクトルに対して同定した結果の一例を示す図である。   FIG. 16 shows an example of the result of identification for the reflectance spectrum shown in FIG. 15 using the layer structure identification process according to the embodiment of the present invention.

図15および図16に示すように、本実施の形態に従う層構造同定処理によれば、多層膜試料の膜構造を高精度に同定することができる。   As shown in FIGS. 15 and 16, according to the layer structure identification process according to the present embodiment, the film structure of the multilayer film sample can be identified with high accuracy.

この発明の実施の形態によれば、多層膜試料から測定された反射率スペクトルからパワースペクトルを算出する過程において、多層膜試料における各波長の波長分散を考慮して、波数の関数である波数変換反射率を導入する。そして、パワースペクトルは、この波数変換反射率の関数に対して、波数についてフーリエ変換することで算出される。さらに、このパワースペクトルに現れるピークに基づいて対象となる試料における膜厚を算出する。このように、パワースペクトルの算出過程において、各波長の波長分散を考慮した独自の関数を採用することで、屈折率の波長依存性の影響を排除して、より高い精度で多層膜試料の膜厚を測定することができる。   According to the embodiment of the present invention, in the process of calculating the power spectrum from the reflectance spectrum measured from the multilayer film sample, the wave number conversion that is a function of the wave number is considered in consideration of the chromatic dispersion of each wavelength in the multilayer film sample. Introduce reflectance. The power spectrum is calculated by performing a Fourier transform on the wave number with respect to the function of the wave number conversion reflectance. Further, the film thickness of the target sample is calculated based on the peak appearing in the power spectrum. In this way, in the process of calculating the power spectrum, by adopting a unique function that takes into account the chromatic dispersion of each wavelength, the influence of the wavelength dependency of the refractive index is eliminated, and the film of the multilayer film sample with higher accuracy. Thickness can be measured.

また、この発明の実施の形態によれば、予め定められた複数の層構造情報から候補となる層構造を順次選択した上で、この候補となる層構造における反射率スペクトル(理論値)を算出し、さらにこの反射率スペクトルからパワースペクトル(理論値)を算出する。そして、実際に測定された反射率スペクトルから算出されるパワースペクトル(実測値)とパワースペクトル(理論値)との相関係数に基づいて、最も確からしい層構造が決定される。これにより、多数の層が積層された多層膜試料であっても、その層構造をより正確に特定できる。   Further, according to the embodiment of the present invention, a candidate layer structure is sequentially selected from a plurality of predetermined layer structure information, and a reflectance spectrum (theoretical value) in the candidate layer structure is calculated. Further, a power spectrum (theoretical value) is calculated from the reflectance spectrum. The most probable layer structure is determined based on the correlation coefficient between the power spectrum (measured value) calculated from the actually measured reflectance spectrum and the power spectrum (theoretical value). Thereby, even if it is a multilayer film sample in which many layers were laminated, the layer structure can be specified more correctly.

今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した説明ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。   The embodiment disclosed this time should be considered as illustrative in all points and not restrictive. The scope of the present invention is defined by the terms of the claims, rather than the description above, and is intended to include any modifications within the scope and meaning equivalent to the terms of the claims.

この発明の実施の形態に従う多層膜解析装置の概略構成図である。1 is a schematic configuration diagram of a multilayer film analysis apparatus according to an embodiment of the present invention. 代表的な多層膜試料の断面模式図である。It is a cross-sectional schematic diagram of a typical multilayer film sample. 測定例1の試料から測定された反射率スペクトルの一例を示す図である。It is a figure which shows an example of the reflectance spectrum measured from the sample of the measurement example 1. FIG. 図3に示す反射率スペクトルから従来の方法を用いて膜厚を算出した場合の結果を示す図である。It is a figure which shows the result at the time of calculating a film thickness using the conventional method from the reflectance spectrum shown in FIG. 図3に示す反射率スペクトルから本実施の形態に従う方法を用いて膜厚を算出した場合の結果を示す図である。It is a figure which shows the result at the time of calculating a film thickness from the reflectance spectrum shown in FIG. 3 using the method according to this Embodiment. 2層の薄膜層が形成された多層膜試料の断面模式図である。It is a cross-sectional schematic diagram of the multilayer film sample in which the two thin film layers were formed. 測定例2の試料から測定された反射率スペクトルの一例を示す図である。It is a figure which shows an example of the reflectance spectrum measured from the sample of the measurement example 2. FIG. 図7に示す反射率スペクトルから変換した第1層に係る波数変換反射率についてのパワースペクトルである。It is a power spectrum about the wave number conversion reflectance which concerns on the 1st layer converted from the reflectance spectrum shown in FIG. 図7に示す反射率スペクトルから変換した第2層に係る波数変換反射率についてのパワースペクトルである。It is a power spectrum about the wave number conversion reflectance which concerns on the 2nd layer converted from the reflectance spectrum shown in FIG. この発明の実施の形態に従うデータ処理部の概略のハードウェア構成を示す模式図である。It is a schematic diagram which shows the schematic hardware constitutions of the data processing part according to embodiment of this invention. この発明の実施の形態に従う膜厚測定処理を実現するための制御構造を示すブロック図である。It is a block diagram which shows the control structure for implement | achieving the film thickness measurement process according to embodiment of this invention. この発明の実施の形態に従う膜厚測定に係る処理手順を示すフローチャートである。It is a flowchart which shows the process sequence which concerns on the film thickness measurement according to embodiment of this invention. この発明の実施の形態に従う層構造同定処理を実現するための制御構造を示すブロック図である。It is a block diagram which shows the control structure for implement | achieving the layer structure identification process according to embodiment of this invention. この発明の実施の形態に従う層構造同定処理に係る処理手順を示すフローチャートである。It is a flowchart which shows the process sequence which concerns on the layer structure identification process according to embodiment of this invention. 任意の試料に対して測定された反射率スペクトルを示す測定例である。It is a measurement example which shows the reflectance spectrum measured with respect to arbitrary samples. この発明の実施の形態に従う層構造同定処理を用いて図15に示す反射率スペクトルに対して同定した結果の一例を示す図である。It is a figure which shows an example of the result identified with respect to the reflectance spectrum shown in FIG. 15 using the layer structure identification process according to embodiment of this invention.

符号の説明Explanation of symbols

10 測定用光源、12 コリメートレンズ、14,66 カットフィルタ、16,36 結像レンズ、18 絞り部、20,30 ビームスプリッタ、22 観察用光源、24 光ファイバ、26 出射部、26a マスク部、32 ピンホールミラー、32a ピンホール、34 軸変換ミラー、38 観察用カメラ、39 表示部、40 対物レンズ、50 ステージ、52 可動機構、60 分光測定部、62 回折格子、64 検出部、68 シャッタ、70 データ処理部、71,731,732,…,73m バッファ部、741,742,…,74m フーリエ変換部、751,752,…,75m ピーク探索部、76 総合算出部、100 多層膜解析装置、204 ディスプレイ部、206 インターフェイス部、208 入力部、210 ハードディスク部(HDD)、212 メモリ部、214 CD−ROMドライブ、214a CD−ROM、216 FDDドライブ、216a フレキシブルディスク、300 層構造データベース、302 層構造情報、304 組合せ部、306 反射率算出部、308 理論スペクトル算出部、310,320 パワースペクトル算出部、318 バッファ部、330 相関係数算出部、332 層構造決定部。   DESCRIPTION OF SYMBOLS 10 Measurement light source, 12 Collimate lens, 14,66 Cut filter, 16,36 Imaging lens, 18 Aperture part, 20, 30 Beam splitter, 22 Observation light source, 24 Optical fiber, 26 Output part, 26a Mask part, 32 Pinhole mirror, 32a Pinhole, 34 axis conversion mirror, 38 Observation camera, 39 Display unit, 40 Objective lens, 50 Stage, 52 Movable mechanism, 60 Spectrometer, 62 Diffraction grating, 64 Detector, 68 Shutter, 70 Data processing unit 71,731,732, ..., 73m Buffer unit, 741,742, ..., 74m Fourier transform unit, 751,752, ..., 75m Peak search unit, 76 Total calculation unit, 100 Multi-layer film analysis device, 204 Display unit, 206 interface unit, 208 input unit, 210 Hard disk unit (HDD), 212 memory unit, 214 CD-ROM drive, 214a CD-ROM, 216 FDD drive, 216a flexible disk, 300 layer structure database, 302 layer structure information, 304 combination unit, 306 reflectance calculation unit, 308 Theoretical spectrum calculation unit, 310, 320 Power spectrum calculation unit, 318 buffer unit, 330 correlation coefficient calculation unit, 332 layer structure determination unit.

Claims (7)

基板上に少なくとも1つの膜が形成された試料に対して所定の波長範囲をもつ測定光を照射する光源と、
前記試料で反射された光または前記試料を透過した光に基づいて、反射率または透過率の波長分布特性を取得する分光部と、
前記波長分布特性における各波長とその波長における反射率または透過率の値との対応関係を、各波長についての波数と所定の関係式に従って算出される変換値との対応関係に変換する波数特性変換手段と、
前記波数特性変換手段によってそれぞれ変換された対応関係からなる波数分布特性に対して、波数について周波数変換する周波数変換手段と、
周波数変換後の波数分布特性に現れるピークに基づいて、前記試料における膜厚を算出する膜厚算出手段とを備える、多層膜解析装置。
A light source that emits measurement light having a predetermined wavelength range to a sample having at least one film formed on a substrate;
A spectroscopic unit that acquires reflectance or transmittance wavelength distribution characteristics based on light reflected by the sample or light transmitted through the sample;
Wave number characteristic conversion for converting the correspondence between each wavelength in the wavelength distribution characteristic and the value of reflectance or transmittance at that wavelength into the correspondence between the wave number for each wavelength and the conversion value calculated according to a predetermined relational expression Means,
Frequency conversion means for frequency-converting the wave number with respect to the wave number distribution characteristics each having a corresponding relationship converted by the wave number characteristic conversion means;
A multilayer film analysis apparatus comprising: a film thickness calculation unit configured to calculate a film thickness of the sample based on a peak appearing in a wave number distribution characteristic after frequency conversion.
前記波数特性変換手段は、予め入力される前記試料に含まれる膜の屈折率と対応の波長とに基づいて前記波数を算出する、請求項1に記載の多層膜解析装置。   The multilayer film analyzer according to claim 1, wherein the wave number characteristic conversion unit calculates the wave number based on a refractive index of a film included in the sample input in advance and a corresponding wavelength. 前記関係式は、各波長における反射率または透過率の値を前記周波数変換手段での周波数変換に係る位相因子に対して線形化するための関数である、請求項1または2に記載の多層膜解析装置。   3. The multilayer film according to claim 1, wherein the relational expression is a function for linearizing a reflectance or transmittance value at each wavelength with respect to a phase factor related to frequency conversion by the frequency conversion unit. Analysis device. 前記波長分布特性は、前記試料で反射された光に基づいて取得される反射率Rの波長分布特性であり、
前記関係式は、R/(1−R)である、請求項3に記載の多層膜解析装置。
The wavelength distribution characteristic is a wavelength distribution characteristic of reflectance R acquired based on light reflected by the sample,
The multilayer analysis apparatus according to claim 3, wherein the relational expression is R / (1-R).
前記波数特性変換手段は、前記試料に複数の膜が形成されている場合に、前記複数の膜のそれぞれについて前記波数分布特性を算出し、
前記膜厚算出手段は、前記複数の膜についての周波数変換後の前記波波数分布特性のそれぞれに現れるピークに基づいて、前記複数の膜の各々の膜厚を算出する、請求項1〜4のいずれか1項に記載の多層膜解析装置。
The wave number characteristic conversion means calculates the wave number distribution characteristic for each of the plurality of films when a plurality of films are formed on the sample,
The film thickness calculation means calculates the film thickness of each of the plurality of films based on a peak appearing in each of the wave number distribution characteristics after frequency conversion for the plurality of films. The multilayer film analysis apparatus according to any one of the above.
少なくとも屈折率および膜厚が規定された層構造情報を複数格納する層構造格納手段と、
複数の前記層構造情報のうち少なくとも1つの層構造情報に基づいて、層構造候補を順次選択する選択手段と、
前記層構造候補に対応する前記層構造情報に基づいて、前記層構造候補における反射率または透過率の波長分布特性を算出する波長分布特性算出手段と、
前記層構造候補の前記波長分布特性から、前記波数特性変換手段および前記周波数変換手段と同一の処理によって周波数変換後の波数分布特性を算出する波数分布特性算出手段と、
実測された前記波長分布特性から算出された周波数変換後の波数分布特性と、前記波数分布特性算出手段によって算出された周波数変換後の波数分布特性との間の相関値を算出する相関値算出手段と、
前記相関値に基づいて、前記選択手段が順次選択する層構造候補のうち実測値に最も近いものを決定する決定手段とをさらに備える、請求項1〜5のいずれか1項に記載の多層膜解析装置。
Layer structure storage means for storing a plurality of layer structure information in which at least the refractive index and the film thickness are defined,
Selection means for sequentially selecting layer structure candidates based on at least one layer structure information among the plurality of layer structure information;
Based on the layer structure information corresponding to the layer structure candidate, wavelength distribution characteristic calculating means for calculating a wavelength distribution characteristic of reflectance or transmittance in the layer structure candidate;
From the wavelength distribution characteristic of the layer structure candidate, a wave number distribution characteristic calculating means for calculating a wave number distribution characteristic after frequency conversion by the same processing as the wave number characteristic converting means and the frequency converting means,
Correlation value calculating means for calculating a correlation value between the wave number distribution characteristic after frequency conversion calculated from the actually measured wavelength distribution characteristic and the wave number distribution characteristic after frequency conversion calculated by the wave number distribution characteristic calculating means. When,
6. The multilayer film according to claim 1, further comprising: a determination unit that determines a layer structure candidate that is sequentially selected by the selection unit based on the correlation value, and that determines the closest layer to the actual measurement value. Analysis device.
基板上に少なくとも1つの膜が形成された試料に対して所定の波長範囲をもつ測定光を照射するステップと、
前記試料で反射された光または前記試料を透過した光に基づいて、反射率または透過率の波長分布特性を取得するステップと、
前記波長分布特性における各波長とその波長における値との対応関係を、各波長についての波数と所定の関係式に従って算出される変換値との対応関係に変換するステップと、
前記対応関係に変換するステップにおいてそれぞれ変換された対応関係からなる波数分布特性に対して、波数について周波数変換するステップと、
周波数変換後の波数分布特性に現れるピークに基づいて、前記試料における膜厚を算出するステップとを備える、多層膜解析方法。
Irradiating a sample having at least one film formed on a substrate with measurement light having a predetermined wavelength range;
Obtaining a reflectance or transmittance wavelength distribution characteristic based on light reflected by the sample or transmitted through the sample;
Converting the correspondence between each wavelength in the wavelength distribution characteristic and the value at that wavelength into a correspondence between the wave number for each wavelength and a conversion value calculated according to a predetermined relational expression;
For the wave number distribution characteristic consisting of the corresponding relationship converted in the step of converting into the corresponding relationship, frequency converting the wave number;
And a step of calculating a film thickness of the sample based on a peak appearing in a wave number distribution characteristic after frequency conversion.
JP2007261840A 2007-10-05 2007-10-05 Multilayer film analysis apparatus and multilayer film analysis method Active JP4834847B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007261840A JP4834847B2 (en) 2007-10-05 2007-10-05 Multilayer film analysis apparatus and multilayer film analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007261840A JP4834847B2 (en) 2007-10-05 2007-10-05 Multilayer film analysis apparatus and multilayer film analysis method

Publications (2)

Publication Number Publication Date
JP2009092454A true JP2009092454A (en) 2009-04-30
JP4834847B2 JP4834847B2 (en) 2011-12-14

Family

ID=40664586

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007261840A Active JP4834847B2 (en) 2007-10-05 2007-10-05 Multilayer film analysis apparatus and multilayer film analysis method

Country Status (1)

Country Link
JP (1) JP4834847B2 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011095138A (en) * 2009-10-30 2011-05-12 Yokogawa Electric Corp Method for measuring film thickness of multilayer film and device therefor
JP2012154931A (en) * 2011-01-26 2012-08-16 Mitsutoyo Corp Method and device for execution of film thickness measurement using white light scanning interferometry
WO2016163352A1 (en) * 2015-04-08 2016-10-13 株式会社 荏原製作所 Film thickness measuring method, film thickness measuring device, polishing method, and polishing device
CN109163666A (en) * 2017-06-05 2019-01-08 大塚电子株式会社 Optical measuring device and measuring method
US10480997B1 (en) 2018-05-18 2019-11-19 Otsuka Electronics Co., Ltd. Optical measurement apparatus and optical measurement method
JP2020122680A (en) * 2019-01-29 2020-08-13 大塚電子株式会社 Optical measuring system and optical measurement method
JP2020176842A (en) * 2019-04-15 2020-10-29 大塚電子株式会社 Optical unit, optical measurement device, and optical measurement method
JP7121222B1 (en) * 2021-05-21 2022-08-17 浜松ホトニクス株式会社 Film thickness measuring device and film thickness measuring method
WO2022244309A1 (en) * 2021-05-21 2022-11-24 浜松ホトニクス株式会社 Film thickness measurement device and film thickness measurement method

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6487579B1 (en) 2018-01-09 2019-03-20 浜松ホトニクス株式会社 Film thickness measuring apparatus, film thickness measuring method, film thickness measuring program, and recording medium for recording film thickness measuring program

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH074922A (en) * 1993-06-21 1995-01-10 Jasco Corp Apparatus and method for measurement of film thickness of semiconductor multilayer thin film
JPH0755435A (en) * 1993-08-20 1995-03-03 Dainippon Screen Mfg Co Ltd Film thickness measuring method for multilayer film sample
JPH0763524A (en) * 1993-08-24 1995-03-10 Jasco Corp Method and device for evaluating thin film
JPH07159131A (en) * 1993-12-07 1995-06-23 Nkk Corp Measurement of ellipso-parameter and ellipsometer
JP2003075127A (en) * 2001-09-06 2003-03-12 Horiba Jobin Yvon Co Ltd Measuring method of ultrathin film and thin film
JP2005037315A (en) * 2003-07-18 2005-02-10 Toppan Printing Co Ltd Film thickness measuring method and film thickness measuring instrument for multilayer film
JP2006504092A (en) * 2002-10-28 2006-02-02 テヴエト プロセス コントロール テクノロジーズ リミテッド Method and apparatus for thickness analysis of complex layer structures
JP2006519993A (en) * 2003-03-06 2006-08-31 ザイゴ コーポレーション Profiling complex surface structures using scanning interferometry
JP2007155630A (en) * 2005-12-08 2007-06-21 Mitsubishi Chemical Engineering Corp Thickness measuring method and thickness measuring apparatus of laminated thin film

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH074922A (en) * 1993-06-21 1995-01-10 Jasco Corp Apparatus and method for measurement of film thickness of semiconductor multilayer thin film
JPH0755435A (en) * 1993-08-20 1995-03-03 Dainippon Screen Mfg Co Ltd Film thickness measuring method for multilayer film sample
JPH0763524A (en) * 1993-08-24 1995-03-10 Jasco Corp Method and device for evaluating thin film
JPH07159131A (en) * 1993-12-07 1995-06-23 Nkk Corp Measurement of ellipso-parameter and ellipsometer
JP2003075127A (en) * 2001-09-06 2003-03-12 Horiba Jobin Yvon Co Ltd Measuring method of ultrathin film and thin film
JP2006504092A (en) * 2002-10-28 2006-02-02 テヴエト プロセス コントロール テクノロジーズ リミテッド Method and apparatus for thickness analysis of complex layer structures
JP2006519993A (en) * 2003-03-06 2006-08-31 ザイゴ コーポレーション Profiling complex surface structures using scanning interferometry
JP2005037315A (en) * 2003-07-18 2005-02-10 Toppan Printing Co Ltd Film thickness measuring method and film thickness measuring instrument for multilayer film
JP2007155630A (en) * 2005-12-08 2007-06-21 Mitsubishi Chemical Engineering Corp Thickness measuring method and thickness measuring apparatus of laminated thin film

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011095138A (en) * 2009-10-30 2011-05-12 Yokogawa Electric Corp Method for measuring film thickness of multilayer film and device therefor
JP2012154931A (en) * 2011-01-26 2012-08-16 Mitsutoyo Corp Method and device for execution of film thickness measurement using white light scanning interferometry
TWI661174B (en) * 2015-04-08 2019-06-01 日商荏原製作所股份有限公司 Film thickness measurement method, film thickness measurement device, polishing method, and polishing device
WO2016163352A1 (en) * 2015-04-08 2016-10-13 株式会社 荏原製作所 Film thickness measuring method, film thickness measuring device, polishing method, and polishing device
JP2016200459A (en) * 2015-04-08 2016-12-01 株式会社荏原製作所 Film thickness measuring method, film thickness measuring device, polishing method, and polishing device
CN107429989A (en) * 2015-04-08 2017-12-01 株式会社荏原制作所 Film thickness measuring method, film thickness measuring device, Ginding process and lapping device
US20180130667A1 (en) * 2015-04-08 2018-05-10 Ebara Corporation Film thickness measuring method, film thickness measuring apparatus, polishing method, and polishing apparatus
US10256104B2 (en) 2015-04-08 2019-04-09 Ebara Corporation Film thickness measuring method, film thickness measuring apparatus, polishing method, and polishing apparatus
US10309767B2 (en) 2017-06-05 2019-06-04 Otsuka Electronics Co., Ltd. Optical measurement apparatus and optical measurement method
US10288412B2 (en) 2017-06-05 2019-05-14 Otsuka Electronics Co., Ltd. Optical measurement apparatus and optical measurement method
CN109163666A (en) * 2017-06-05 2019-01-08 大塚电子株式会社 Optical measuring device and measuring method
CN109163666B (en) * 2017-06-05 2021-06-11 大塚电子株式会社 Optical measuring device and optical measuring method
TWI759495B (en) * 2017-06-05 2022-04-01 日商大塚電子股份有限公司 Optical measurement apparatus and optical measurement method
US10480997B1 (en) 2018-05-18 2019-11-19 Otsuka Electronics Co., Ltd. Optical measurement apparatus and optical measurement method
KR20190132243A (en) 2018-05-18 2019-11-27 오츠카 일렉트로닉스 가부시키가이샤 Optical measuring device and optical measuring method
JP2020122680A (en) * 2019-01-29 2020-08-13 大塚電子株式会社 Optical measuring system and optical measurement method
JP7199093B2 (en) 2019-01-29 2023-01-05 大塚電子株式会社 Optical measurement system and optical measurement method
JP2020176842A (en) * 2019-04-15 2020-10-29 大塚電子株式会社 Optical unit, optical measurement device, and optical measurement method
JP7219463B2 (en) 2019-04-15 2023-02-08 大塚電子株式会社 Optical unit, optical measuring device and optical measuring method
JP7121222B1 (en) * 2021-05-21 2022-08-17 浜松ホトニクス株式会社 Film thickness measuring device and film thickness measuring method
WO2022244309A1 (en) * 2021-05-21 2022-11-24 浜松ホトニクス株式会社 Film thickness measurement device and film thickness measurement method

Also Published As

Publication number Publication date
JP4834847B2 (en) 2011-12-14

Similar Documents

Publication Publication Date Title
JP4834847B2 (en) Multilayer film analysis apparatus and multilayer film analysis method
JP5309359B2 (en) Film thickness measuring apparatus and film thickness measuring method
KR101841776B1 (en) Optical characteristic measuring apparatus and optical characteristic measuring method
KR20130018553A (en) Film thickness measurement apparatus
JP5172203B2 (en) Optical characteristic measuring apparatus and measuring method
US10365163B2 (en) Optical critical dimension metrology
US8947673B2 (en) Estimating thickness based on number of peaks between two peaks in scanning white light interferometry
TWI292030B (en) High density multi-channel detecting device
JP5172204B2 (en) Optical characteristic measuring apparatus and focus adjustment method
US20060098206A1 (en) Apparatus and method for measuring thickness and profile of transparent thin film using white-light interferometer
JP5172040B2 (en) Surface shape measuring method and surface shape measuring apparatus
JP5394317B2 (en) Rotationally symmetric aspherical shape measuring device
JPH07134007A (en) Measuring device for film thickness of thin-film by high space resolution
JP5517621B2 (en) High sensitivity spectrum analysis unit
JP2010002328A (en) Film thickness measuring instrument
US7292337B2 (en) Optical processor using detecting assembly and method using same
TWI793321B (en) Optical measuring device and optical measuring method
JP6725988B2 (en) Thickness measuring device and thickness measuring method
US20140029006A1 (en) Surface plasmon sensor and method of measuring refractive index
KR101388424B1 (en) Apparatus for measuring a thickness using digital light processing and method using the same
JP5028660B2 (en) Optical characteristic measuring apparatus and optical characteristic measuring method
JP2007155480A (en) Surface measurement apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110712

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110713

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110812

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: 20110830

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20110902

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20141007

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4834847

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250