JP2009194356A - Etching amount calculating method, storage medium, and etching amount calculating apparatus - Google Patents

Etching amount calculating method, storage medium, and etching amount calculating apparatus Download PDF

Info

Publication number
JP2009194356A
JP2009194356A JP2008192311A JP2008192311A JP2009194356A JP 2009194356 A JP2009194356 A JP 2009194356A JP 2008192311 A JP2008192311 A JP 2008192311A JP 2008192311 A JP2008192311 A JP 2008192311A JP 2009194356 A JP2009194356 A JP 2009194356A
Authority
JP
Japan
Prior art keywords
interference
period
etching amount
interference wave
waveform
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
JP2008192311A
Other languages
Japanese (ja)
Other versions
JP5026363B2 (en
Inventor
Susumu Saito
進 斉藤
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.)
Tokyo Electron Ltd
Original Assignee
Tokyo Electron 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 Tokyo Electron Ltd filed Critical Tokyo Electron Ltd
Priority to JP2008192311A priority Critical patent/JP5026363B2/en
Priority to CN2008101888359A priority patent/CN101488456B/en
Priority to KR1020090001096A priority patent/KR101073229B1/en
Priority to US12/353,457 priority patent/US8014891B2/en
Priority to TW098101548A priority patent/TWI508161B/en
Publication of JP2009194356A publication Critical patent/JP2009194356A/en
Application granted granted Critical
Publication of JP5026363B2 publication Critical patent/JP5026363B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Drying Of Semiconductors (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide an etching amount calculating method which can stably and accurately calculate the amount of etching even if a disturbance is added. <P>SOLUTION: In the etching of a wafer W which forms a trench 132 using a mask film 131, laser beam L1 is irradiated on the wafer W. Superposed interference wave is calculated by receiving the superposed interference light resulting from superposition of mask film interference light and trench interference light. A waveform of a window which defines the present timing T as the endpoint is extracted. Max entropy method is used for the waveform of the window, and the trench interference light is detected from frequency distribution obtained by performing frequency analysis. While shifting the endpoint of the window by only Δt, the calculation of the superposed interference wave, the extraction of the waveform of the window, the frequency analysis of the wave form of the window, and the detection of the trench interference period are repeated, and the detected trench interference periods are integrated and averaged at each repetition. The etching amount of the trench 132 is calculated based on the integrated and averaged trench interference periods. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

本発明は、エッチング量算出方法、記憶媒体及びエッチング量算出装置に関し、特にマスク膜を用いてウエハ上にトレンチやホール等の凹部を形成する際のエッチング量算出方法に関する。   The present invention relates to an etching amount calculation method, a storage medium, and an etching amount calculation device, and more particularly, to an etching amount calculation method when a recess such as a trench or a hole is formed on a wafer using a mask film.

半導体デバイスの製造工程では、ウエハにおいてマスク膜を用いて被エッチング層にトレンチやホールを形成するエッチングが行われている。エッチングではマスク膜で覆われていない部分の被エッチング層がプラズマによって物理的・化学的に削られていくが、トレンチの形成では該トレンチの深さを制御する必要がある。したがって、エッチング中にトレンチの深さ、すなわち、エッチング量を算出する必要があるが、従来、エッチング量の算出方法としては光の干渉を利用する方法が広く用いられている。   In the manufacturing process of a semiconductor device, etching for forming a trench or a hole in a layer to be etched is performed using a mask film on a wafer. In etching, a portion of the layer to be etched that is not covered with the mask film is physically and chemically removed by plasma. However, in forming a trench, it is necessary to control the depth of the trench. Therefore, it is necessary to calculate the depth of the trench, that is, the etching amount during etching, but conventionally, a method using light interference has been widely used as a method for calculating the etching amount.

図22は、エッチング中における光の干渉を説明するための図である。   FIG. 22 is a diagram for explaining light interference during etching.

図22において、被エッチング層130上に形成されたマスク膜131を有するウエハWにはエッチングによってトレンチ132が形成されているが、該ウエハWへレーザ光L1を照射すると、マスク膜131の表面からの反射光L、マスク膜131及び被エッチング層130の境界面からの反射光L及びトレンチ132の底からの反射光Lが生じる。 In FIG. 22, a trench 132 is formed by etching on a wafer W having a mask film 131 formed on the etching target layer 130. When the wafer W is irradiated with laser light L 1, the surface of the mask film 131 is exposed. Reflected light L 2 , reflected light L 3 from the interface between the mask film 131 and the etching target layer 130, and reflected light L 4 from the bottom of the trench 132 are generated.

反射光L〜Lを検出器で受光する際、図22に示すように、各反射光の光路長がマスク膜131の厚さやトレンチ132の深さ分だけ異なるため、検出器の受光面では各反射光の位相が異なり、干渉光(例えば、反射光L及び反射光Lの干渉光(以下「トレンチ干渉光」という。)や反射光L及び反射光Lの干渉光(以下「マスク膜干渉光」という。))が発生する。 When the reflected lights L 2 to L 4 are received by the detector, the optical path length of each reflected light differs by the thickness of the mask film 131 and the depth of the trench 132 as shown in FIG. Then, the phase of each reflected light is different, and interference light (for example, interference light of reflected light L 2 and reflected light L 4 (hereinafter referred to as “trench interference light”) or interference light of reflected light L 2 and reflected light L 3 ( Hereinafter, it is referred to as “mask film interference light”)).

そして、エッチング中においてトレンチ132の深さは時々刻々と変化するため、反射光Lと反射光Lの光路長差も時々刻々と変化して干渉光の強度が変化する、すなわち、反射光Lと反射光Lから干渉波(以下「トレンチ干渉波」という。)が発生する。干渉波の周期はトレンチ132の深さの変化速度(エッチングレート)によって決定されるため、干渉波の周期からエッチングレートを算出し、さらには算出されたエッチングレートとエッチング時間とからエッチング量(トレンチ132の深さ)を算出することができる。 Since the depth of the trench 132 changes every moment during the etching, the optical path length difference between the reflected light L 2 and the reflected light L 4 also changes every moment and the intensity of the interference light changes, that is, the reflected light. Interference waves (hereinafter referred to as “trench interference waves”) are generated from L 2 and reflected light L 4 . Since the period of the interference wave is determined by the change rate (etching rate) of the depth of the trench 132, the etching rate is calculated from the period of the interference wave, and further, the etching amount (trench) is calculated from the calculated etching rate and etching time. 132 depth) can be calculated.

ところで、エッチング中ではマスク膜131も微量ずつエッチングされて厚さが変化するため、反射光L及び反射光Lからも干渉波(以下「マスク膜干渉波」という。)が発生する。各干渉波は同一の検出器に検出されるため、該検出器が検出した干渉波は異なる周期を有する複数の干渉波が重畳されたもの(以下、「重畳干渉波」という。)となる(図23参照。)。 By the way, during etching, the mask film 131 is also etched by a small amount to change its thickness, so that interference waves (hereinafter referred to as “mask film interference waves”) are also generated from the reflected light L 2 and the reflected light L 3 . Since each interference wave is detected by the same detector, the interference wave detected by the detector is a superposition of a plurality of interference waves having different periods (hereinafter, referred to as “superimposed interference wave”) ( (See FIG. 23).

図23に示すような重畳干渉波からトレンチ132の深さ(被エッチング層130のエッチング量)を算出するには、重畳干渉波からトレンチ干渉波を分離する必要がある。   In order to calculate the depth of the trench 132 (the etching amount of the etching target layer 130) from the superimposed interference wave as shown in FIG. 23, it is necessary to separate the trench interference wave from the superimposed interference wave.

図23の重畳干渉波では短周期の干渉波と長周期の干渉波が比較的明確に分離可能である。ここで、エッチング中におけるトレンチ132の深さの変化速度がマスク膜131の厚さの変化速度よりも大きいため、トレンチ干渉波はマスク膜干渉波よりも周期が短い。したがって、図23の重畳干渉波における短周期の干渉波がトレンチ干渉波であり、短周期の干渉波における極値間の時間(図中の「Δt」)からトレンチ干渉波の周期を容易に算出することができる。   In the superimposed interference wave of FIG. 23, the short-cycle interference wave and the long-cycle interference wave can be separated relatively clearly. Here, since the rate of change of the depth of the trench 132 during etching is larger than the rate of change of the thickness of the mask film 131, the trench interference wave has a shorter period than the mask film interference wave. Therefore, the short-period interference wave in the superimposed interference wave in FIG. 23 is a trench interference wave, and the period of the trench interference wave is easily calculated from the time between the extreme values in the short-cycle interference wave (“Δt” in the figure). can do.

重畳干渉波から極値間の時間を読み取る方法では、重畳干渉波において短周期の干渉波と長周期の干渉波が比較的明確に分離可能である必要があるため、短周期の干渉波と長周期の干渉波が分離困難な重畳干渉波ではトレンチ干渉波の周期を算出することができない。また、短周期の干渉波における極値間ではトレンチ干渉波の周期が一定であると見なすため、算出されるトレンチ干渉波の周期(被エッチング層130のエッチングレート)は、図24に示すように、段々状になる。すなわち、重畳干渉波から極値間の時間を読み取る方法は分解能が低い。   In the method of reading the time between extreme values from the superimposed interference wave, it is necessary that the short-period interference wave and the long-period interference wave are relatively clearly separable in the superimposed interference wave. With a superimposed interference wave in which it is difficult to separate the interference wave of the period, the period of the trench interference wave cannot be calculated. Further, since the period of the trench interference wave is assumed to be constant between the extreme values in the short-period interference wave, the calculated period of the trench interference wave (etching rate of the etching target layer 130) is as shown in FIG. , Step by step. That is, the method of reading the time between extreme values from the superimposed interference wave has low resolution.

そこで、近年、重畳干渉波から極値間の時間を読み取ることなく周波数解析によってトレンチ干渉波の周期を算出する方法が開発されている。この方法では、重畳干渉波から周波数解析(例えば、高速フーリエ変換法)によって周波数分布(図26(A)参照。)を得て、該周波数分布からトレンチ干渉波の周期を検出する(例えば、特許文献1参照。)。
特開平2−71517号公報
Therefore, in recent years, a method of calculating the period of the trench interference wave by frequency analysis without reading the time between extreme values from the superimposed interference wave has been developed. In this method, the frequency distribution (see FIG. 26A) is obtained from the superimposed interference wave by frequency analysis (for example, fast Fourier transform method), and the period of the trench interference wave is detected from the frequency distribution (for example, patents). Reference 1).
Japanese Patent Laid-Open No. 2-71517

しかしながら、図25に示すように、重畳干渉波においてレーザ光源や検出器の異常(図中「I」)又はマスク膜干渉波とトレンチ干渉波の干渉による見かけの周期変化(図中「II」)等の外乱が重畳干渉波に加わることがある。上述した周波数解析を用いる方法ではエッチング中の全時間に亘る重畳干渉波を解析して得られた周波数分布からトレンチ干渉波の周期を検出するだけなので、外乱が重畳干渉波に加わった場合、本来存在しない干渉周期にピークが発生する等して周波数分布が不正確になり(図26(B)参照。)、その結果、エッチング量の算出を安定して正確に行うことができない。   However, as shown in FIG. 25, an apparent period change (“II” in the figure) due to an abnormality in the laser light source or detector (“I” in the figure) or interference between the mask film interference wave and the trench interference wave in the superimposed interference wave. Such disturbances may be added to the superimposed interference wave. In the method using the frequency analysis described above, the period of the trench interference wave is only detected from the frequency distribution obtained by analyzing the superimposed interference wave over the entire time during the etching. The frequency distribution becomes inaccurate because a peak occurs in a non-existing interference cycle (see FIG. 26B), and as a result, the etching amount cannot be calculated stably and accurately.

本発明の目的は、外乱が加わってもエッチング量の算出を安定して正確に行うことができるエッチング量算出方法、記憶媒体及びエッチング量算出装置を提供することにある。   An object of the present invention is to provide an etching amount calculation method, a storage medium, and an etching amount calculation apparatus capable of stably and accurately calculating an etching amount even when a disturbance is applied.

上記目的を達成するために、請求項1記載のエッチング量算出方法は、マスク膜を用いて凹部を形成する基板のエッチングにおいて前記凹部のエッチング量を算出するエッチング量算出方法であって、前記基板に光を照射する照射ステップと、少なくとも前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉光が他の干渉光に重畳された重畳干渉光を受光する受光ステップと、前記受光された重畳干渉光から重畳干渉波を算出する干渉波算出ステップと、前記重畳干渉波から所定期間の波形を抽出する波形抽出ステップと、前記抽出された波形に周波数解析を施す周波数解析ステップと、前記周波数解析によって得られた周波数分布から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出する干渉周期検出ステップと、前記所定期間を所定時間だけずらしつつ前記干渉波算出ステップ、前記波形抽出ステップ、前記周波数解析ステップ及び前記干渉周期検出ステップを繰り返し、繰り返しの度に前記検出された干渉波の周期を積算平均する積算平均ステップと、前記積算平均された干渉波の周期に基づいて前記凹部のエッチング量を算出するエッチング量算出ステップとを有することを特徴とする。   In order to achieve the above object, the etching amount calculation method according to claim 1 is an etching amount calculation method for calculating the etching amount of the recess in etching of the substrate in which the recess is formed using a mask film, An irradiation step for irradiating light, a light receiving step for receiving superimposed interference light in which interference light of at least reflected light from the mask film and reflected light from the bottom of the recess is superimposed on other interference light, and the light receiving An interference wave calculation step of calculating a superimposed interference wave from the superimposed interference light, a waveform extraction step of extracting a waveform of a predetermined period from the superimposed interference wave, a frequency analysis step of performing frequency analysis on the extracted waveform, An interference period for detecting the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess from the frequency distribution obtained by the frequency analysis And repeating the step of calculating the interference wave, the step of extracting the waveform, the step of analyzing the frequency, and the step of detecting the interference period while shifting the predetermined period by a predetermined time, and the period of the detected interference wave is repeated each time it is repeated. An integration averaging step for integrating and averaging, and an etching amount calculating step for calculating the etching amount of the concave portion based on the period of the integrated and averaged interference wave are provided.

請求項2記載のエッチング量算出方法は、請求項1記載のエッチング量算出方法において、前記所定期間は、前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波よりも周期が長い前記他の干渉光の波形の1周期よりも大きいことを特徴とする。   The etching amount calculation method according to claim 2 is the etching amount calculation method according to claim 1, wherein the predetermined period has a period longer than an interference wave of reflected light from the mask film and reflected light from the bottom of the recess. It is characterized by being longer than one period of the long waveform of the other interference light.

請求項3記載のエッチング量算出方法は、請求項1又は2記載のエッチング量算出方法において、前記他の干渉光の波形の周期が前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期よりも長い場合に、前記周波数解析ステップに先立って、前記重畳干渉波から抽出された所定期間の波形から前記他の干渉光の波形が占める部分の殆どを除去する解析前処理ステップをさらに有し、前記周波数解析ステップでは、前記他の干渉光の波形が占める部分の殆どが除去された波形に周波数解析を施すことを特徴とする。   The etching amount calculation method according to claim 3 is the etching amount calculation method according to claim 1 or 2, wherein the waveform period of the other interference light is reflected light from the mask film and reflected light from the bottom of the recess. Pre-analysis processing for removing most of the portion occupied by the waveform of the other interference light from the waveform of the predetermined period extracted from the superimposed interference wave prior to the frequency analysis step when the period of the interference wave is longer than The frequency analysis step further includes performing frequency analysis on the waveform from which most of the portion occupied by the waveform of the other interference light has been removed.

請求項4記載のエッチング量算出方法は、請求項3記載のエッチング量算出方法において、前記解析前処理ステップでは、前記抽出された波形から、該抽出された波形を二次多項式で近似した波形を除去することを特徴とする。   The etching amount calculation method according to claim 4 is the etching amount calculation method according to claim 3, wherein, in the pre-analysis processing step, a waveform obtained by approximating the extracted waveform with a second-order polynomial is extracted from the extracted waveform. It is characterized by removing.

請求項5記載のエッチング量算出方法は、請求項3又は4記載のエッチング量算出方法において、前記所定期間は前記他の干渉光の波形の1/4周期以下であることを特徴とする。   The etching amount calculation method according to claim 5 is the etching amount calculation method according to claim 3 or 4, wherein the predetermined period is equal to or less than ¼ period of the waveform of the other interference light.

請求項6記載のエッチング量算出方法は、請求項3乃至5のいずれか1項に記載のエッチング量算出方法において、前記基板の表面における前記凹部の開口率が0.5%以下であるか、又は前記凹部がディープトレンチであることを特徴とする。   The etching amount calculation method according to claim 6 is the etching amount calculation method according to any one of claims 3 to 5, wherein the aperture ratio of the recesses on the surface of the substrate is 0.5% or less. Alternatively, the concave portion is a deep trench.

請求項7記載のエッチング量算出方法は、請求項1乃至6のいずれか1項に記載のエッチング量算出方法において、前記周波数解析では最大エントロピー法を用いることを特徴とする。   The etching amount calculation method according to claim 7 is the etching amount calculation method according to any one of claims 1 to 6, wherein the frequency analysis uses a maximum entropy method.

請求項8記載のエッチング量算出方法は、請求項1乃至7のいずれか1項に記載のエッチング量算出方法において、前記周波数分布から検出された前記干渉波の周期が異常値に該当する場合には該干渉波の周期を除去する干渉周期修正ステップをさらに有することを特徴とする。   The etching amount calculation method according to claim 8 is the etching amount calculation method according to any one of claims 1 to 7, wherein a period of the interference wave detected from the frequency distribution corresponds to an abnormal value. Further includes an interference period correcting step for removing the period of the interference wave.

請求項9記載のエッチング量算出方法は、請求項8記載のエッチング量算出方法において、前記干渉周期修正ステップでは、前記異常値に該当する前記干渉波の周期が求められた前記所定期間の前の前記所定期間又は後の前記所定期間から求められる前記干渉波の周期を、前記異常値に該当する前記干渉波の周期が求められた前記所定期間の干渉波の周期とみなすことを特徴とする。   The etching amount calculation method according to claim 9 is the etching amount calculation method according to claim 8, wherein, in the interference period correction step, a period before the predetermined period in which the period of the interference wave corresponding to the abnormal value is obtained. The period of the interference wave obtained from the predetermined period or the subsequent predetermined period is regarded as the period of the interference wave of the predetermined period in which the period of the interference wave corresponding to the abnormal value is obtained.

請求項10記載のエッチング量算出方法は、請求項1乃至9のいずれか1項に記載のエッチング量算出方法において、前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を予め予測し、前記干渉周期検出ステップでは、前記周波数解析によって得られた周波数分布において、前記予測された周期近傍から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出することを特徴とする。   The etching amount calculation method according to claim 10 is the etching amount calculation method according to any one of claims 1 to 9, wherein interference waves of reflected light from the mask film and reflected light from the bottom of the recess are reduced. In the interference period detection step, in the frequency distribution obtained by the frequency analysis, an interference wave of reflected light from the mask film and reflected light from the bottom of the recess in the frequency distribution obtained by the frequency analysis is predicted. The period is detected.

請求項11記載のエッチング量算出方法は、請求項1乃至10のいずれか1項に記載のエッチング量算出方法において、前記他の干渉光は、前記マスク膜の表面からの反射光、並びに前記マスク膜及び前記基板の表面の境界面からの反射光の干渉光であることを特徴とする。   The etching amount calculation method according to claim 11 is the etching amount calculation method according to any one of claims 1 to 10, wherein the other interference light includes reflected light from the surface of the mask film, and the mask. The interference light is reflected light from the boundary surface between the film and the surface of the substrate.

上記目的を達成するために、請求項12記載の記憶媒体は、マスク膜を用いて凹部を形成する基板のエッチングにおいて前記凹部のエッチング量を算出するエッチング量算出方法をコンピュータに実行させるプログラムを格納するコンピュータで読み取り可能な記憶媒体であって、前記エッチング量算出方法は、前記基板に光を照射する照射ステップと、少なくとも前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉光が他の干渉光に重畳された重畳干渉光を受光する受光ステップと、前記受光された重畳干渉光から重畳干渉波を算出する干渉波算出ステップと、前記重畳干渉波から所定期間の波形を抽出する波形抽出ステップと、前記抽出された波形に周波数解析を施す周波数解析ステップと、前記周波数解析によって得られた周波数分布から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出する干渉周期検出ステップと、前記所定期間を所定時間だけずらしつつ前記干渉波算出ステップ、前記波形抽出ステップ、前記周波数解析ステップ及び前記干渉周期検出ステップを繰り返し、繰り返しの度に前記検出された干渉波の周期を積算平均する積算平均ステップと、前記積算平均された干渉波の周期に基づいて前記凹部のエッチング量を算出するエッチング量算出ステップとを有することを特徴とする。   In order to achieve the above object, a storage medium according to claim 12 stores a program for causing a computer to execute an etching amount calculation method for calculating an etching amount of the recess in etching a substrate on which the recess is formed using a mask film. The etching amount calculation method includes: an irradiation step of irradiating the substrate with light; and interference light of at least reflected light from the mask film and reflected light from the bottom of the recess Receiving a superimposed interference light superimposed on another interference light, an interference wave calculating step calculating a superimposed interference wave from the received superimposed interference light, and extracting a waveform of a predetermined period from the superimposed interference wave A waveform extracting step, a frequency analyzing step for performing frequency analysis on the extracted waveform, and a frequency analysis step obtained by the frequency analysis. An interference period detecting step for detecting a period of interference waves of reflected light from the mask film and reflected light from the bottom of the recess from a frequency distribution; the interference wave calculating step while shifting the predetermined period by a predetermined time; and the waveform The extraction step, the frequency analysis step, and the interference cycle detection step are repeated, and the integration average step of integrating and averaging the detected interference wave cycle each time it is repeated, and based on the integration averaged interference wave cycle And an etching amount calculating step for calculating the etching amount of the recess.

上記目的を達成するために、請求項13記載のエッチング量算出装置は、マスク膜を用いて凹部を形成する基板のエッチングにおいて前記凹部のエッチング量を算出するエッチング量算出装置において、前記基板に光を照射する照射部と、少なくとも前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉光が他の干渉光に重畳された重畳干渉光を受光する受光部と、前記受光された重畳干渉光から重畳干渉波を算出する干渉波算出部と、前記重畳干渉波から所定期間の波形を抽出する波形抽出部と、前記抽出された波形に周波数解析を施す周波数解析部と、前記周波数解析によって得られた周波数分布から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出する干渉周期検出部と、前記所定期間を所定時間だけずらしつつ前記重畳干渉波の算出、前記所定期間の波形の抽出、前記周波数解析及び前記干渉波の周期の検出を繰り返し、繰り返しの度に前記検出された干渉波の周期を積算平均する積算平均部と、前記積算平均された干渉波の周期に基づいて前記凹部のエッチング量を算出するエッチング量算出部とを備えることを特徴とする。   In order to achieve the above object, an etching amount calculation apparatus according to claim 13 is an etching amount calculation apparatus that calculates an etching amount of the recess in etching of the substrate on which the recess is formed using a mask film. An irradiating unit that irradiates the interference light of at least the reflected light from the mask film and the reflected light from the bottom of the recess, and the light receiving unit that receives the superimposed interference light superimposed on the other interference light; and An interference wave calculation unit that calculates a superimposed interference wave from the superimposed interference light, a waveform extraction unit that extracts a waveform of a predetermined period from the superimposed interference wave, a frequency analysis unit that performs frequency analysis on the extracted waveform, and the frequency An interference period detection unit for detecting a period of interference waves of reflected light from the mask film and reflected light from the bottom of the recess from a frequency distribution obtained by analysis; and the predetermined period is predetermined. Accumulation that repeats calculation of the superimposed interference wave, extraction of the waveform during the predetermined period, frequency analysis, and detection of the period of the interference wave while shifting the interval, and integrating and averaging the period of the detected interference wave at each repetition An average portion and an etching amount calculation unit that calculates an etching amount of the concave portion based on the cycle of the integrated and averaged interference wave are provided.

請求項1記載のエッチング量算出方法、請求項12記載の記憶媒体及び請求項13記載のエッチング量算出装置によれば、所定期間を所定時間だけずらしつつ重畳干渉波の算出、所定期間の波形の抽出、周波数解析、並びにマスク膜からの反射光及び凹部の底からの反射光の干渉波の周期の検出を繰り返し、繰り返しの度に検出された干渉波の周期を積算平均し、積算平均された干渉波の周期に基づいて凹部のエッチング量を算出する。したがって、例え、抽出された或る所定期間の波形に外乱が加わる場合であっても、当該或る所定期間の波形に基づいて検出された干渉波の周期は、他の所定期間の波形に基づいて検出された干渉波の周期と積算平均されるので、外乱が加わった所定期間の波形に基づいて検出された干渉波の周期が、積算平均される干渉波の周期に及ぼす影響を小さくすることができ、もって、外乱が加わってもエッチング量の算出を安定して正確に行うことができる。   According to the etching amount calculation method according to claim 1, the storage medium according to claim 12, and the etching amount calculation device according to claim 13, the superimposed interference wave is calculated while the predetermined period is shifted by a predetermined time, and the waveform of the predetermined period is calculated. Extraction, frequency analysis, and detection of the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess were repeated, and the period of the interference wave detected at each repetition was integrated and averaged. The etching amount of the recess is calculated based on the period of the interference wave. Therefore, even if a disturbance is applied to the extracted waveform of a certain predetermined period, the period of the interference wave detected based on the waveform of the certain predetermined period is based on the waveform of another predetermined period. Therefore, the influence of the period of the interference wave detected based on the waveform of the predetermined period with the disturbance on the period of the interference wave to be averaged is reduced. Therefore, the etching amount can be calculated stably and accurately even when a disturbance is applied.

請求項2記載のエッチング量算出方法によれば、上記所定期間は、マスク膜からの反射光及び凹部の底からの反射光の干渉波よりも周期が長い他の干渉光の波形の1周期よりも大きいので、所定期間の波形の周波数解析の信頼性を向上することができ、もって、エッチング量の算出をより正確に行うことができる。   According to the etching amount calculation method according to claim 2, the predetermined period is longer than one period of the waveform of the other interference light having a longer period than the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess. Therefore, the reliability of the frequency analysis of the waveform during the predetermined period can be improved, and the etching amount can be calculated more accurately.

請求項3記載のエッチング量算出方法によれば、他の干渉光の波形の周期がマスク膜からの反射光及び凹部の底からの反射光の干渉波の周期よりも長い場合に、周波数解析に先立って、重畳干渉波から抽出された所定期間の波形から他の干渉光の波形が占める部分の殆どが除去されるので、凹部の底からの反射光の光量が少ない場合であって、他の干渉光の波形が重畳干渉波の殆どを占める場合であっても、除去後の波形においてマスク膜からの反射光及び凹部の底からの反射光の干渉波が占める部分の割合を大きくすることができ、もって、周波数解析においてマスク膜からの反射光及び凹部の底からの反射光の干渉波の周期を正確に算出することができる。   According to the etching amount calculation method according to claim 3, when the period of the waveform of the other interference light is longer than the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess, the frequency analysis is performed. First, since most of the portion occupied by the waveform of the other interference light is removed from the waveform of the predetermined period extracted from the superimposed interference wave, the amount of reflected light from the bottom of the recess is small, Even when the waveform of the interference light occupies most of the superimposed interference wave, the ratio of the portion occupied by the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess in the waveform after removal may be increased. Therefore, the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess can be accurately calculated in the frequency analysis.

請求項4記載のエッチング量算出方法によれば、周波数解析に先立って抽出された波形から、該抽出された波形を二次多項式で近似した波形が除去される。凹部の底からの反射光の光量が少ない場合であって、他の干渉光の波形が重畳干渉波の殆どを占める場合、重畳干渉波の波形は他の干渉光の波形とほぼ等しくなるので、抽出された重畳干渉波を二次多項式で近似した波形も他の干渉光の波形とほぼ等しくなる。したがって、抽出された波形から他の干渉光の波形が占める部分の殆どを確実に除去することができる。   According to the etching amount calculation method of the fourth aspect, a waveform obtained by approximating the extracted waveform with a second order polynomial is removed from the waveform extracted prior to the frequency analysis. When the amount of reflected light from the bottom of the recess is small, and the waveform of the other interference light occupies most of the superimposed interference wave, the waveform of the superimposed interference wave is almost equal to the waveform of the other interference light, A waveform obtained by approximating the extracted superimposed interference wave with a second-order polynomial is also substantially equal to the waveform of the other interference light. Therefore, most of the portion occupied by the waveform of the other interference light can be reliably removed from the extracted waveform.

請求項5記載のエッチング量算出方法によれば、上記所定期間は他の干渉光の波形の1/4周期以下である。重畳干渉波の殆どを占める他の干渉光の波形は正弦波に近いので、他の干渉光の波形の1/4周期以下の部分を抽出すれば、該抽出された波形は二次多項式によって正確に近似することができる。これにより、抽出された波形から他の干渉光の波形が占める部分の殆どを正確に除去することができる。   According to the etching amount calculation method of the fifth aspect, the predetermined period is equal to or shorter than a quarter cycle of the waveform of the other interference light. Since the other interfering light waveform that occupies most of the superimposed interference wave is close to a sine wave, if the portion of the other interfering light waveform that is less than ¼ period is extracted, the extracted waveform is accurately expressed by a second order polynomial. Can be approximated. Thereby, most of the portion occupied by the waveform of the other interference light can be accurately removed from the extracted waveform.

請求項7記載のエッチング量算出方法によれば、周波数解析では最大エントロピー法を用いるので、所定期間の波形の数が少なくても周波数解析の信頼性を向上することができ、もって、エッチング量の算出をより正確に行うことができる。   According to the etching amount calculation method of claim 7, since the maximum entropy method is used in the frequency analysis, the reliability of the frequency analysis can be improved even if the number of waveforms in the predetermined period is small. Calculation can be performed more accurately.

請求項8記載のエッチング量算出方法によれば、周波数分布から検出された干渉波の周期が異常値に該当する場合には該干渉波の周期を除去するので、外乱が加わった所定期間の波形に基づいて検出された干渉波の周期の積算平均される干渉波の周期に対する影響を除去することができ、もって、外乱が加わってもエッチング量の算出をより安定して正確に行うことができる。   According to the etching amount calculation method according to claim 8, when the period of the interference wave detected from the frequency distribution corresponds to an abnormal value, the period of the interference wave is removed. Thus, the influence of the period of the interference wave detected based on the interference wave period on the integrated average can be removed, and the etching amount can be calculated more stably and accurately even when disturbance is applied. .

請求項9記載のエッチング量算出方法によれば、異常値に該当する干渉波の周期が求められた所定期間の前の所定期間又は後の所定期間から求められる干渉波の周期を、異常値に該当する干渉波の周期が求められた所定期間の干渉波の周期とみなすので、外乱が加わった所定期間の波形に基づいて検出された干渉光の影響を確実に除去することができる。   According to the etching amount calculation method according to claim 9, the period of the interference wave obtained from the predetermined period before or after the predetermined period from which the period of the interference wave corresponding to the abnormal value is obtained is set to the abnormal value. Since the period of the corresponding interference wave is regarded as the period of the interference wave obtained for the predetermined period, the influence of the interference light detected based on the waveform of the predetermined period to which the disturbance is added can be reliably removed.

請求項10記載のエッチング量算出方法によれば、マスク膜からの反射光及び凹部の底からの反射光の干渉波の周期を予め予測し、所定期間の波形の周波数解析によって得られた周波数分布において、予測された周期近傍から干渉波の周期を検出するので、干渉波の周期の検出を迅速に行うことができるとともに、異常値が検出されるのを抑制することができる。   According to the etching amount calculation method according to claim 10, the frequency distribution obtained by predicting in advance the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the concave portion, and analyzing the waveform in a predetermined period. Therefore, since the period of the interference wave is detected from the vicinity of the predicted period, it is possible to quickly detect the period of the interference wave and to suppress the detection of the abnormal value.

以下、本発明の実施の形態について図面を参照しながら説明する。   Hereinafter, embodiments of the present invention will be described with reference to the drawings.

まず、本発明の第1の実施の形態に係るエッチング量算出方法及びエッチング量算出方法が適用される基板処理装置について説明する。この基板処理装置は基板としての半導体ウエハ(以下、単に「ウエハ」という。)Wにプラズマを用いたエッチングを施すように構成されている。なお、ウエハWは、前述の図22に示すように、被エッチング層130と、該被エッチング層130上に所定のパターンで形成されたマスク膜131とを有する。   First, the etching amount calculation method and the substrate processing apparatus to which the etching amount calculation method according to the first embodiment of the present invention is applied will be described. This substrate processing apparatus is configured to perform etching using plasma on a semiconductor wafer (hereinafter simply referred to as “wafer”) W as a substrate. Note that the wafer W includes the etching target layer 130 and the mask film 131 formed in a predetermined pattern on the etching target layer 130 as shown in FIG.

図1は、本実施の形態に係るエッチング量算出方法が適用される基板処理装置の構成を概略的に示す断面図である。   FIG. 1 is a cross-sectional view schematically showing a configuration of a substrate processing apparatus to which the etching amount calculation method according to the present embodiment is applied.

図1において、基板処理装置10は、例えば、アルミニウム等の導電性材料からなる処理室11と、ウエハWを載置する載置台として処理室11内の底面に配設される下部電極12と、該下部電極12の上方に所定の間隔を隔てて配設されたシャワーヘッド13とを備える。   In FIG. 1, a substrate processing apparatus 10 includes, for example, a processing chamber 11 made of a conductive material such as aluminum, a lower electrode 12 disposed on a bottom surface in the processing chamber 11 as a mounting table on which a wafer W is mounted, And a shower head 13 disposed above the lower electrode 12 at a predetermined interval.

処理室11の下部には真空排気装置(図示しない)が接続された排気部14が接続され、下部電極12には整合器15を介して高周波電源16が接続され、シャワーヘッド13の内部のバッファ室17には処理ガス導入管18が接続され、該処理ガス導入管18には処理ガス供給装置19が接続されている。シャワーヘッド13は下部において、バッファ室17と、シャワーヘッド13及び下部電極12の間の空間である処理空間Sとを連通させる複数のガス穴20を有する。シャワーヘッド13は、処理ガス導入管18からバッファ室17に導入された処理ガスを、複数のガス穴20を介して処理空間Sに供給する。   An exhaust unit 14 connected to a vacuum exhaust device (not shown) is connected to the lower part of the processing chamber 11, and a high frequency power supply 16 is connected to the lower electrode 12 via a matching unit 15, and a buffer inside the shower head 13. A processing gas introduction pipe 18 is connected to the chamber 17, and a processing gas supply device 19 is connected to the processing gas introduction pipe 18. The shower head 13 has a plurality of gas holes 20 in the lower portion for communicating the buffer chamber 17 with the processing space S that is a space between the shower head 13 and the lower electrode 12. The shower head 13 supplies the processing gas introduced into the buffer chamber 17 from the processing gas introduction pipe 18 into the processing space S through the plurality of gas holes 20.

この基板処理装置10では、排気部14によって処理室11内を所定の真空度まで減圧した後、下部電極12から処理空間Sに高周波電圧を印加した状態で、シャワーヘッド13から処理空間Sへ処理ガスを供給し、処理空間Sにおいて処理ガスからプラズマを発生させる。該発生したプラズマは、ウエハWにおいてマスク膜131に覆われていない被エッチング層130に衝突・接触して該被エッチング層130をエッチングし、該被エッチング層130にトレンチ132(凹部)を形成する。   In the substrate processing apparatus 10, after the processing chamber 11 is depressurized to a predetermined degree of vacuum by the exhaust unit 14, the processing is performed from the shower head 13 to the processing space S in a state where a high frequency voltage is applied from the lower electrode 12 to the processing space S. Gas is supplied, and plasma is generated from the processing gas in the processing space S. The generated plasma collides and contacts the etching target layer 130 not covered with the mask film 131 on the wafer W to etch the etching target layer 130, thereby forming a trench 132 (recess) in the etching target layer 130. .

処理室11内のシャワーヘッド13には下部電極12に載置されたウエハWを上方から観測するためのモニタ装置21が配設されている。モニタ装置21は円筒状の部材からなり、シャワーヘッド13を貫通する。モニタ装置21の上端には石英ガラス等の透明体からなる窓部材22が設けられている。また、処理室11の上方にはモニタ装置21の上端と集光レンズ23を介して対向する光ファイバ24が配置されている。   The shower head 13 in the processing chamber 11 is provided with a monitor device 21 for observing the wafer W placed on the lower electrode 12 from above. The monitor device 21 is made of a cylindrical member and penetrates the shower head 13. A window member 22 made of a transparent material such as quartz glass is provided at the upper end of the monitor device 21. In addition, an optical fiber 24 is disposed above the processing chamber 11 so as to face the upper end of the monitor device 21 via a condenser lens 23.

光ファイバ24は、被エッチング層130のエッチング量を算出するエッチング量算出装置25に接続されている。エッチング量算出装置25は、光ファイバ24にそれぞれ接続されたレーザ光源26(照射部)及び検出器27(受光部)と、検出器27に接続された演算部28(干渉波算出部、波形抽出部、周波数解析部、干渉周期検出部、積算平均部、エッチング量算出部)とを備え、基板処理装置10のコントローラ29の制御下で作動する。レーザ光源26としては、例えば、半導体レーザが用いられる。また、検出器27としては、例えば、ホトマルやフォトダイオードが用いられる。なお、コントローラ29は、演算部28だけでなく、基板処理装置10の各構成要素、例えば、高周波電源16にも接続され、各構成要素の動作を制御する。   The optical fiber 24 is connected to an etching amount calculation device 25 that calculates the etching amount of the etching target layer 130. The etching amount calculation device 25 includes a laser light source 26 (irradiation unit) and a detector 27 (light receiving unit) connected to the optical fiber 24, and a calculation unit 28 (interference wave calculation unit, waveform extraction) connected to the detector 27. Unit, frequency analysis unit, interference period detection unit, integration average unit, etching amount calculation unit) and operate under the control of the controller 29 of the substrate processing apparatus 10. As the laser light source 26, for example, a semiconductor laser is used. As the detector 27, for example, a photomultiplier or a photodiode is used. The controller 29 is connected not only to the computing unit 28 but also to each component of the substrate processing apparatus 10, for example, the high-frequency power supply 16, and controls the operation of each component.

エッチング量算出装置25は、レーザ光源26からのレーザ光を光ファイバ24、集光レンズ23及びモニタ装置21を介して下部電極12上のウエハWに向けて照射するとともに、ウエハWからの反射光、即ち、トレンチ干渉光(マスク膜からの反射光及び凹部の底からの反射光の干渉光)やマスク膜干渉光(他の干渉光)が重畳された重畳干渉光を、光ファイバ24等を介して検出器27によって受光する。検出器27によって受光された重畳干渉光は電気信号に変換されて演算部28に送信される。   The etching amount calculation device 25 irradiates the laser light from the laser light source 26 toward the wafer W on the lower electrode 12 via the optical fiber 24, the condenser lens 23, and the monitor device 21, and also reflects the reflected light from the wafer W. That is, trench interference light (interference light of reflected light from the mask film and reflected light from the bottom of the recess) and superposed interference light on which mask film interference light (other interference light) is superimposed are transmitted through the optical fiber 24 and the like. The light is received by the detector 27. The superimposed interference light received by the detector 27 is converted into an electrical signal and transmitted to the calculation unit 28.

演算部28は受信した電気信号に基づいて重畳干渉光から重畳干渉波を算出する。また、演算部28は、算出された重畳干渉波に基づいて、後述する図6のエッチング量算出方法を実行してトレンチ132のエッチング量を算出する。   The calculation unit 28 calculates a superimposed interference wave from the superimposed interference light based on the received electrical signal. Further, the calculation unit 28 calculates the etching amount of the trench 132 by executing an etching amount calculation method of FIG. 6 described later based on the calculated superimposed interference wave.

ところで、エッチングではエッチングレートが一定ではなく、種々の要因(処理空間Sの圧力の変化や高周波電圧の乱れ)によって変化する。特に、アスペクト比の大きいトレンチ(例えば、ディープトレンチ)をエッチングにて形成する場合には、トレンチ132のエッチング量(エッチング深さ)が大きくなると、トレンチ132の入り口にデポが付着してトレンチ132内へのプラズマの進入が抑制されるため、エッチングレートが低下する(図2参照。)。また、図2に示すように、エッチングレートは微小変化を繰り返すため、トレンチ132のエッチング量を正確に算出するには小刻みにエッチングレートを算出する必要がある。   By the way, in etching, the etching rate is not constant, and varies depending on various factors (change in pressure in the processing space S and disturbance of high-frequency voltage). In particular, when a trench having a large aspect ratio (for example, a deep trench) is formed by etching, when the etching amount (etching depth) of the trench 132 is increased, deposits are attached to the entrance of the trench 132 and the trench 132 is filled. Since the entry of plasma into the substrate is suppressed, the etching rate is reduced (see FIG. 2). Further, as shown in FIG. 2, since the etching rate repeats minute changes, it is necessary to calculate the etching rate in small increments in order to accurately calculate the etching amount of the trench 132.

エッチングにおいて小刻みにエッチングレートを算出するには、通常、ウエハWからの反射光の波形を各タイミングにおいて微分すればよいが、上述したように、ウエハWからの反射光はトレンチ干渉光やマスク膜干渉光が重畳された重畳干渉光であるため、単純に反射光の波形を各タイミングにおいて微分してもトレンチ132のエッチングレートを正確に求めることができない。   In order to calculate the etching rate in small increments in etching, it is usually sufficient to differentiate the waveform of the reflected light from the wafer W at each timing. However, as described above, the reflected light from the wafer W is reflected by the trench interference light or the mask film. Since the interference light is superimposed interference light, the etching rate of the trench 132 cannot be obtained accurately even if the waveform of the reflected light is simply differentiated at each timing.

そこで、本実施の形態では周波数解析を行い、重畳干渉波からトレンチ干渉波の周期(以下、「トレンチ干渉周期」という。)を算出する。また、周波数解析では一定以上のデータ長、具体的には、解析対象の波形の1周期以上のデータ長を用いることが周波数解析の信頼性向上に寄与することから、あるタイミングのエッチングレートを求めるために、本実施の形態では、重畳干渉波から所定期間の波形を抽出し、該抽出された波形を周波数解析する。また、上述したように、エッチングレートはトレンチ干渉周期と関連があるため、本実施の形態では、抽出された波形から、まず、トレンチ干渉周期を求め、該トレンチ干渉波からエッチングレートを算出する。   Therefore, in the present embodiment, frequency analysis is performed, and the period of the trench interference wave (hereinafter referred to as “trench interference period”) is calculated from the superimposed interference wave. In addition, in the frequency analysis, using a data length of a certain value or more, specifically, using a data length of one cycle or more of the waveform to be analyzed contributes to improving the reliability of the frequency analysis. Therefore, an etching rate at a certain timing is obtained. Therefore, in this embodiment, a waveform for a predetermined period is extracted from the superimposed interference wave, and the extracted waveform is subjected to frequency analysis. As described above, since the etching rate is related to the trench interference period, in this embodiment, first, the trench interference period is obtained from the extracted waveform, and the etching rate is calculated from the trench interference wave.

図3は、本実施の形態に係るエッチング量算出方法における重畳干渉波からの所定期間の波形の抽出を説明するための図である。   FIG. 3 is a diagram for explaining extraction of a waveform for a predetermined period from the superimposed interference wave in the etching amount calculation method according to the present embodiment.

図3において、重畳干渉波30は約30秒周期で振動しているので、上記所定期間を30秒間に設定する。ここでは、タイミングAのトレンチ干渉周期を求めるために、タイミングAから30秒前までの期間31における重畳干渉波30の波形(図中の四角形で囲った部分の波形)を抽出する。なお、本実施の形態では上記所定期間を以下「窓」と称する。窓31は始点32及び終点33を有し、始点32はタイミングAから30秒前に該当し、終点33はタイミングAに該当する。   In FIG. 3, since the superimposed interference wave 30 vibrates at a cycle of about 30 seconds, the predetermined period is set to 30 seconds. Here, in order to obtain the trench interference period at the timing A, the waveform of the superimposed interference wave 30 in the period 31 from the timing A to 30 seconds before (the waveform of the portion surrounded by a square in the drawing) is extracted. In the present embodiment, the predetermined period is hereinafter referred to as “window”. The window 31 has a start point 32 and an end point 33. The start point 32 corresponds to 30 seconds before the timing A, and the end point 33 corresponds to the timing A.

そして、本実施の形態では、抽出された窓31の波形に周波数解析を施す。ここで、抽出された波形はせいぜい1周期分であるため、解析方法として最大エントロピー法を用いる。最大エントロピー法は観測現象の周波数分布を極めて短い測定時間の測定結果から高分解能で算出する方法であり(「科学計算のための波形データ処理」(CQ出版社,昭和61年4月30日初版発行)参照。)、解析対象の波形の数をさほど必要としないため、何周期分もの波形を必要とする高速フーリエ変換法よりも本実施の形態に係るエッチング量算出方法に適している。   In this embodiment, frequency analysis is performed on the extracted waveform of the window 31. Here, since the extracted waveform is at most one cycle, the maximum entropy method is used as an analysis method. The maximum entropy method is a method of calculating the frequency distribution of observation phenomena from the measurement results of extremely short measurement time with high resolution ("Waveform data processing for scientific calculation" (CQ Publishing Co., Ltd., April 30, 1986, first edition) Since the number of waveforms to be analyzed is not so large, it is more suitable for the etching amount calculation method according to the present embodiment than the fast Fourier transform method that requires many cycles of waveforms.

図4は、図3における窓の波形から最大エントロピー法を用いた周波数解析によって得られた周波数分布を示す図である。   FIG. 4 is a diagram showing a frequency distribution obtained by frequency analysis using the maximum entropy method from the waveform of the window in FIG.

重畳干渉波30には主としてトレンチ干渉波とマスク膜干渉波の2つが含まれているため、窓31の波形から得られる周波数分布には、図4に示すように、ピークを示す周波数(干渉波の周期)が主として2つ(図4では約0.012Hzと約0.037Hz)存在する。ここで、上述したように、トレンチ干渉波はマスク膜干渉波よりも周期が短い(周波数が高い)ため、約0.037Hzの周波数がトレンチ干渉周期に該当する。そこで、約0.037Hzの周波数をトレンチ干渉周期として検出する。このように、本実施の形態では、周波数解析によって得られる周波数分布には2つのピークが存在することが予見されることから、周波数解析に先立ってトレンチ干渉周期を予め予測し、得られた周波数分布において予測されたトレンチ干渉周期の近傍からトレンチ干渉周期を検出するのが好ましい。   Since the superimposed interference wave 30 mainly includes two of the trench interference wave and the mask film interference wave, the frequency distribution obtained from the waveform of the window 31 has a peak frequency (interference wave) as shown in FIG. There are mainly two periods (about 0.012 Hz and about 0.037 Hz in FIG. 4). Here, as described above, since the trench interference wave has a shorter cycle (higher frequency) than the mask film interference wave, a frequency of about 0.037 Hz corresponds to the trench interference cycle. Therefore, a frequency of about 0.037 Hz is detected as the trench interference period. Thus, in this embodiment, since it is predicted that there are two peaks in the frequency distribution obtained by the frequency analysis, the trench interference period is predicted in advance prior to the frequency analysis, and the obtained frequency is obtained. It is preferable to detect the trench interference period from the vicinity of the trench interference period predicted in the distribution.

窓31は、タイミングAから30秒前までの重畳干渉波30の波形を含むため、図4に示す周波数分布はタイミングAから30秒前までの重畳干渉波30における周波数分布となる。したがって、図4に示す周波数分布から検出されたトレンチ干渉周期は、タイミングAから30秒前までの重畳干渉波30におけるトレンチ干渉光の平均周期となるが、本実施の形態では、図4に示す周波数分布から検出されたトレンチ干渉周期を便宜的にタイミングAにおけるトレンチ干渉周期とみなす。なお、本実施の形態では、後述するように重畳干渉波30に対して複数の窓を設定し、各窓の波形の周波数分布から得られたトレンチ干渉周期を積算平均してトレンチ干渉周期の全体平均値を算出するため、窓31におけるトレンチ干渉光の平均周期をタイミングAにおけるトレンチ干渉周期とみなすことの弊害が解消される。   Since the window 31 includes the waveform of the superimposed interference wave 30 30 seconds before the timing A, the frequency distribution shown in FIG. 4 is the frequency distribution in the superimposed interference wave 30 30 seconds before the timing A. Therefore, the trench interference period detected from the frequency distribution shown in FIG. 4 is the average period of the trench interference light in the superimposed interference wave 30 from the timing A to 30 seconds before, but in this embodiment, it is shown in FIG. The trench interference period detected from the frequency distribution is regarded as a trench interference period at timing A for convenience. In the present embodiment, as will be described later, a plurality of windows are set for the superimposed interference wave 30, and the trench interference period obtained from the frequency distribution of the waveform of each window is integrated and averaged to obtain the entire trench interference period. Since the average value is calculated, the adverse effect of considering the average period of the trench interference light in the window 31 as the trench interference period at the timing A is eliminated.

また、本実施の形態に係るエッチング量算出方法では、最初の所定期間(始点32がエッチング開始時に該当し且つ終点33がエッチング開始時から30秒後に該当する窓31に対応する)経過時からエッチング量を算出するタイミングまでの全期間におけるエッチングレートの平均値を算出し、該エッチングレートの平均値からエッチング量を算出する。   In addition, in the etching amount calculation method according to the present embodiment, the etching is performed from the time when the first predetermined period (the start point 32 corresponds to the start of etching and the end point 33 corresponds to the corresponding window 31 30 seconds after the start of etching) has elapsed. The average value of the etching rate over the entire period until the timing for calculating the amount is calculated, and the etching amount is calculated from the average value of the etching rate.

図5は、本実施の形態に係るエッチング量算出方法におけるエッチングレートの平均値の算出方法を説明するための図であり、エッチング開始後80秒が経過した場合を示す。   FIG. 5 is a diagram for explaining a method for calculating an average value of etching rates in the etching amount calculating method according to the present embodiment, and shows a case where 80 seconds have elapsed after the start of etching.

図5では、重畳干渉波50に対してΔt(所定時間)だけずらされたn個の窓W(k=1〜n、nは自然数)が設定され、各窓Wにおいて周波数分布が求められ、各周波数分布からn個のトレンチ干渉周期f(k=1〜n、nは自然数)が検出される。 In FIG. 5, n windows W k (k = 1 to n, where n is a natural number) shifted by Δt (predetermined time) with respect to the superimposed interference wave 50 are set, and the frequency distribution is obtained in each window W k . N trench interference periods f k (k = 1 to n, n is a natural number) are detected from each frequency distribution.

次いで、n個のトレンチ干渉周期fが下記式(1)に基づいて積算平均されて Next, n trench interference periods f k are averaged based on the following formula (1).

Figure 2009194356
Figure 2009194356

トレンチ干渉周期faveがエッチング開始後80秒までのトレンチ干渉波の平均値として算出される。さらに測定波長(レーザ光源26からのレーザ光の波長)をλとすると下記式(2)からエッチング開始後80秒までのエッチングレートの平均値が算出される。エッチングレートの平均値 = fave×λ/2 … (2)
次いで、下記式(3)からエッチング開始後80秒までのエッチング量が算出される。エッチング量 = エッチングレートの平均値×エッチング時間 … (3)
本実施の形態に係るエッチング量算出方法において、或る窓W(tは1〜nのいずれかの自然数)における重畳干渉波50に外乱が加わると、該窓Wから求められたトレンチ干渉周期ftは異常値となるが、該トレンチ干渉周期ftは他の窓W(uは1〜nのいずれかであって、t以外の自然数)から求められたトレンチ干渉周期fと積算平均されるため、トレンチ干渉周期ftが積算平均されたトレンチ干渉周期faveに及ぼす影響は小さい。
The trench interference period f ave is calculated as the average value of the trench interference wave up to 80 seconds after the start of etching. Further, if the measurement wavelength (the wavelength of the laser beam from the laser light source 26) is λ, the average value of the etching rate up to 80 seconds after the start of etching is calculated from the following equation (2). Average value of etching rate = f ave × λ / 2 (2)
Next, the etching amount up to 80 seconds after the start of etching is calculated from the following formula (3). Etching amount = Average value of etching rate x Etching time (3)
In the etching amount calculating method according to the present embodiment, when (the t one of natural numbers 1 to n) one window W t disturbance is applied to the superposed interference wave 50 in the trench interference obtained from the window W t Although the period ft is an abnormal value, the trench interference period ft is an integrated average of the trench interference period f u obtained from another window W u (where u is any one of 1 to n and is a natural number other than t). Therefore , the influence of the trench interference period ft on the integrated and averaged trench interference period f ave is small.

次に、本実施の形態に係るエッチング量算出方法について説明する。   Next, an etching amount calculation method according to the present embodiment will be described.

図6は、本実施の形態に係るエッチング量算出方法を示すフローチャートである。   FIG. 6 is a flowchart showing an etching amount calculation method according to the present embodiment.

図6において、まず、基板処理装置10がウエハWの被エッチング層130のエッチングを開始した後、レーザ光源26がレーザ光L1を光ファイバ24、集光レンズ23及びモニタ装置21を介してウエハWに向けて照射し(ステップS61)(照射ステップ)、検出器27がウエハWからの反射光である重畳干渉光を光ファイバ24等を介して受光する(ステップS62)(受光ステップ)。   In FIG. 6, first, after the substrate processing apparatus 10 starts etching the etching target layer 130 of the wafer W, the laser light source 26 transmits the laser light L <b> 1 through the optical fiber 24, the condenser lens 23, and the monitor device 21. (Step S61) (irradiation step), and the detector 27 receives superimposed interference light, which is reflected light from the wafer W, via the optical fiber 24 or the like (step S62) (light reception step).

次いで、ステップS63において、演算部28は、現タイミングTが予め設定されたエッチングの終了時間を過ぎているか否かを判別し、エッチングの終了時間を過ぎている場合(ステップS63でYES)には本処理を終了し、エッチングの終了時間を過ぎていない場合(ステップS63でNO)には、演算部28は検出器27によって受光された重畳干渉光に基づいてエッチング開始時から現タイミングTまでの重畳干渉波を算出(更新)する(ステップS64)(干渉波算出ステップ)。   Next, in step S63, the calculation unit 28 determines whether or not the current timing T has passed a preset etching end time, and if the etching end time has passed (YES in step S63). When this process is finished and the etching end time has not passed (NO in step S63), the calculation unit 28 determines from the etching start to the current timing T based on the superimposed interference light received by the detector 27. The superimposed interference wave is calculated (updated) (step S64) (interference wave calculation step).

次いで、演算部28は、現タイミングTを終点とする窓の波形を抽出し(ステップS65)(波形抽出ステップ)、該抽出された窓の波形に最大エントロピー法を用いて周波数解析を施す(ステップS66)(周波数解析ステップ)。ここで窓の始点から終点までの時間はマスク膜干渉波の1周期よりも長く設定されている。   Next, the computing unit 28 extracts the waveform of the window whose end point is the current timing T (step S65) (waveform extraction step), and performs frequency analysis using the maximum entropy method on the extracted waveform of the window (step S65). S66) (frequency analysis step). Here, the time from the start point to the end point of the window is set longer than one period of the mask film interference wave.

その後、演算部28は、周波数解析によって得られた周波数分布において、予め予測されたトレンチ干渉周期近傍におけるピークを示す周波数を、現タイミングTにおけるトレンチ干渉周期として検出する(ステップS67)(干渉周期検出ステップ)。   Thereafter, the calculation unit 28 detects, as a trench interference period at the current timing T, a frequency indicating a peak near the trench interference period predicted in advance in the frequency distribution obtained by the frequency analysis (step S67) (interference period detection). Step).

次いで、演算部28は、今回検出された現タイミングTにおけるトレンチ干渉周期と、最初の所定期間経過時から現タイミングTまでの間に検出された各タイミングにおけるトレンチ干渉周期とを上記式(1)を用いて積算平均し(ステップS68)(積算平均ステップ)、積算平均されたトレンチ干渉周期を、エッチング開始時から現タイミングTまでの間における測定波長及びエッチング時間(エッチング開始時から現タイミングTまでの時間)に基づき、上記式(2),(3)を用いてトレンチ132のエッチング量に換算する(ステップS69)(エッチング量算出ステップ)。   Next, the calculation unit 28 calculates the trench interference period at the current timing T detected this time and the trench interference period at each timing detected between the elapse of the first predetermined period and the current timing T by the above equation (1). (Step S68) (integration average step), and the integrated and averaged trench interference period is measured from the etching start time to the current timing T and the etching time (from the etching start time to the current timing T). Is converted into the etching amount of the trench 132 using the above formulas (2) and (3) (step S69) (etching amount calculating step).

その後、演算部28は、現タイミングTにΔtを加えて現タイミングTを更新する、すなわち、窓の終点をΔtだけずらし(ステップS70)、ステップS63に戻る。   Thereafter, the calculation unit 28 adds Δt to the current timing T to update the current timing T, that is, shifts the window end point by Δt (step S70), and returns to step S63.

本実施の形態に係るエッチング量算出方法によれば、窓の終点をΔtだけずらしつつ重畳干渉波の算出、重畳干渉波からの窓の波形の抽出、周波数解析及び現タイミングTにおけるトレンチ干渉周期の検出を繰り返し、繰り返しの度に検出された現タイミングTにおけるトレンチ干渉周期と、最初の所定期間経過時から現タイミングTまでの間に検出された各タイミングにおけるトレンチ干渉周期とを積算平均し、該積算平均されたトレンチ干渉周期をトレンチ132のエッチング量に換算する。したがって、例え、抽出された或る窓の波形に外乱が加わる場合であっても、当該或る窓から求められた異常値であるトレンチ干渉周期は、他の窓から求められたトレンチ干渉周期と積算平均されるので、異常値であるトレンチ干渉周期が積算平均されたトレンチ干渉周期に及ぼす影響を小さくすることができ、もって、外乱が加わってもトレンチ132のエッチング量の算出を安定して正確に行うことができる。   According to the etching amount calculation method according to the present embodiment, the superimposed interference wave is calculated while the end point of the window is shifted by Δt, the window waveform is extracted from the superimposed interference wave, the frequency analysis, and the trench interference period at the current timing T. The detection is repeated, and the trench interference period at the current timing T detected at each repetition and the trench interference period at each timing detected between the elapse of the first predetermined period and the current timing T are integrated and averaged, The averaged trench interference period is converted into the etching amount of the trench 132. Therefore, even if a disturbance is added to the waveform of a certain window extracted, the trench interference period which is an abnormal value obtained from the certain window is the same as the trench interference period obtained from the other window. Since the integrated averaging is performed, the influence of the trench interference period, which is an abnormal value, on the integrated averaged trench interference period can be reduced, so that the calculation of the etching amount of the trench 132 can be performed stably and accurately even when a disturbance is applied. Can be done.

また、本実施の形態に係るエッチング量算出方法では、上記窓の始点から終点までの時間はステップS64で算出された重畳干渉波の1周期よりも大きいので、窓における重畳干渉波の周波数解析の信頼性を向上することができる。   In the etching amount calculation method according to the present embodiment, the time from the start point to the end point of the window is longer than one period of the superimposed interference wave calculated in step S64. Reliability can be improved.

さらに、本実施の形態に係るエッチング量算出方法では、周波数解析において最大エントロピー法を用いるので、窓の波形の数が少なくても周波数解析の信頼性を向上することができる。   Furthermore, since the maximum entropy method is used in the frequency analysis in the etching amount calculation method according to the present embodiment, the reliability of the frequency analysis can be improved even if the number of window waveforms is small.

また、本実施の形態に係るエッチング量算出方法では、トレンチ干渉周期を予め予測し、周波数解析によって得られた周波数分布において、予測されたトレンチ干渉周期近傍から現タイミングTにおけるトレンチ干渉周期を検出するので、トレンチ干渉周期の検出を迅速に行うことができるとともに、トレンチ干渉周期として異常値が検出されるのを抑制することができる。   In the etching amount calculation method according to the present embodiment, the trench interference period is predicted in advance, and the trench interference period at the current timing T is detected from the vicinity of the predicted trench interference period in the frequency distribution obtained by frequency analysis. Therefore, the detection of the trench interference period can be performed quickly, and the detection of an abnormal value as the trench interference period can be suppressed.

次に、本発明の第2の実施の形態に係るエッチング量算出方法について説明する。   Next, an etching amount calculation method according to the second embodiment of the present invention will be described.

本実施の形態は、その構成、作用が上述した第1の実施の形態と基本的に同じであるので、重複した構成、作用については説明を省略し、以下に異なる構成、作用についての説明を行う。   Since the configuration and operation of this embodiment are basically the same as those of the first embodiment described above, the description of the overlapping configuration and operation will be omitted, and the description of the different configuration and operation will be described below. Do.

図7は、本実施の形態に係るエッチング量算出方法を示すフローチャートである。   FIG. 7 is a flowchart showing an etching amount calculation method according to the present embodiment.

図7において、まず、ステップS61乃至S67が実行され、次いで、ステップS71において、演算部28が、ステップS67において現タイミングTにおけるトレンチ干渉周期として検出されたトレンチ干渉周期は異常値(例えば、ステップS66で得られた周波数分布における最大周波数又は最小周波数)に該当するか否かを判別する。   In FIG. 7, first, steps S61 to S67 are executed, and then, in step S71, the calculation unit 28 detects the trench interference period detected as the trench interference period at the current timing T in step S67 (for example, step S66). It is discriminated whether it corresponds to the maximum frequency or the minimum frequency in the frequency distribution obtained in (1).

ステップS71の判別の結果、現タイミングTにおけるトレンチ干渉周期として検出されたトレンチ干渉周期が異常値に該当しない場合(ステップS71でNO)には、そのままステップS68に進み、検出された干渉周期が異常値に該当する場合(ステップSでYES)には、現タイミングTにおけるトレンチ干渉周期として検出されたトレンチ干渉周期を除去するとともに、現タイミングTより1つ前のタイミングに対応する窓から検出されたトレンチ干渉周期を、現タイミングTにおけるトレンチ干渉周期として設定することにより、トレンチ干渉周期を修正する(ステップS72)(干渉周期修正ステップ)。   As a result of the determination in step S71, if the trench interference period detected as the trench interference period at the current timing T does not correspond to the abnormal value (NO in step S71), the process proceeds to step S68 as it is, and the detected interference period is abnormal. If the value corresponds to the value (YES in step S), the trench interference period detected as the trench interference period at the current timing T is removed, and it is detected from the window corresponding to the timing immediately before the current timing T. By setting the trench interference period as the trench interference period at the current timing T, the trench interference period is corrected (step S72) (interference period correcting step).

次いで、演算部28は、ステップS68乃至S70を実行する。   Next, the calculation unit 28 executes Steps S68 to S70.

本実施の形態に係るエッチング量算出方法によれば、現タイミングTにおけるトレンチ干渉周期として検出されたトレンチ干渉周期が異常値に該当する場合には該検出されたトレンチ干渉周期を除去するとともに、現タイミングTより1つ前のタイミングに対応する窓から検出されたトレンチ干渉周期を、現タイミングTにおけるトレンチ干渉周期として設定するので、異常値に該当するトレンチ干渉周期の積算平均されたトレンチ干渉周期に対する影響を除去することができ、もって、外乱が加わってもトレンチ132のエッチング量の算出をより安定して正確に行うことができる。   According to the etching amount calculation method according to the present embodiment, when the trench interference period detected as the trench interference period at the current timing T corresponds to an abnormal value, the detected trench interference period is removed and the current Since the trench interference period detected from the window corresponding to the timing before the timing T is set as the trench interference period at the current timing T, the trench interference period corresponding to the abnormal value is integrated with respect to the averaged trench interference period. The influence can be removed, so that the amount of etching of the trench 132 can be calculated more stably and accurately even when a disturbance is applied.

上述した本実施の形態に係るエッチング量算出方法では、検出されたトレンチ干渉周期が異常値に該当する場合には、現タイミングTより1つ前のタイミングに対応する窓から検出されたトレンチ干渉周期を、現タイミングTにおけるトレンチ干渉周期として設定したが、現タイミングTより1つ後のタイミングに対応する窓から検出されたトレンチ干渉周期を、現タイミングTにおけるトレンチ干渉周期として設定してもよい。   In the etching amount calculation method according to the present embodiment described above, when the detected trench interference period corresponds to an abnormal value, the trench interference period detected from the window corresponding to the timing immediately before the current timing T. Is set as the trench interference period at the current timing T. However, the trench interference period detected from the window corresponding to the timing immediately after the current timing T may be set as the trench interference period at the current timing T.

次に、本発明の第3の実施の形態に係るエッチング量算出方法について説明する。   Next, an etching amount calculation method according to the third embodiment of the present invention will be described.

本実施の形態は、その構成、作用が上述した第1の実施の形態と基本的に同じであり、抽出された窓の波形に周波数解析を施す前に、該窓の波形に前処理を施す点において異なるのみであるので、重複した構成、作用については説明を省略し、以下に異なる構成、作用についての説明を行う。   This embodiment is basically the same in configuration and operation as the first embodiment described above, and pre-processes the extracted waveform of the window before performing frequency analysis on the extracted waveform. Since only the points are different, the description of the duplicated configuration and action will be omitted, and the different configuration and action will be described below.

ウエハWの表面においてトレンチ132の開口部が占める割合の百分率(以下、「開口率」という。)が小さい場合、例えば、0.5%を下回る場合、トレンチ132の底からの反射光Lの絶対光量が少なくなるため、検出器27が受光する重畳干渉光においてトレンチ干渉光が占める部分の割合が小さくなる。 When the percentage of the ratio of the opening of the trench 132 to the surface of the wafer W (hereinafter referred to as “opening ratio”) is small, for example, below 0.5%, the reflected light L 4 from the bottom of the trench 132 is reduced. Since the absolute light amount is reduced, the ratio of the portion occupied by the trench interference light in the superimposed interference light received by the detector 27 is reduced.

図9は、開口率が変化した場合における重畳干渉波の一部の変化を示す図である。   FIG. 9 is a diagram illustrating a change in part of the superimposed interference wave when the aperture ratio changes.

図9に示すように、開口率が5%の場合、重畳干渉波では2種類の干渉波(マスク膜干渉波、トレンチ干渉波)が重畳されていることが明りょうであるが、開口率が0.5%の場合、反射光Lの絶対光量が少なくなり、重畳干渉波にはトレンチ干渉波の波形が殆ど現れない。このような現象は、ウエハWの表面においてホールの開口部が占める割合が小さい場合やトレンチ(又はホール)のアスペクト比が大きい場合(例えば、トレンチ132がディープトレンチである場合)においても、トレンチやホールの底からの反射光の絶対光量が少なくなるため、起こりうる。 As shown in FIG. 9, when the aperture ratio is 5%, it is clear that two types of interference waves (mask film interference wave and trench interference wave) are superimposed on the superimposed interference wave. If 0.5%, the absolute amount of the reflected light L 4 is reduced, almost no waveform of the trench interference wave to the superposition interference wave. Such a phenomenon occurs even when the ratio of the hole openings on the surface of the wafer W is small or when the aspect ratio of the trench (or hole) is large (for example, when the trench 132 is a deep trench). This can happen because the absolute amount of reflected light from the bottom of the hole is reduced.

重畳干渉波にトレンチ干渉波の波形が殆ど現れない場合(開口率が0.5%の場合)、該重畳干渉波から窓31の波形を抽出し、抽出された波形をそのまま周波数解析すると、得られる周波数分布においてトレンチ干渉波の周期(トレンチ干渉周期)のピークが小さくなる。   When the waveform of the trench interference wave hardly appears in the superimposed interference wave (when the aperture ratio is 0.5%), the waveform of the window 31 is extracted from the superimposed interference wave, and the extracted waveform is directly subjected to frequency analysis. In the frequency distribution obtained, the peak of the period of the trench interference wave (trench interference period) becomes small.

図10は、開口率が変化した場合における窓の波形の周波数分布の変化を示す図である。   FIG. 10 is a diagram illustrating a change in the frequency distribution of the waveform of the window when the aperture ratio changes.

図10に示すように、開口率が5%の場合、周波数分布において2つのピーク(トレンチ干渉波の周期(約0.8Hz)、マスク膜干渉波の周期(約0.1Hz))が明りょうに現れているが、開口率が0.5%の場合、周波数分布において1つのピーク(マスク膜干渉波の周期)のみが明りょうに現れてトレンチ干渉波の周期が殆ど現れていない。その結果、トレンチ干渉周期を正確に検出することができず、エッチレートを正確に算出することができない。   As shown in FIG. 10, when the aperture ratio is 5%, two peaks in the frequency distribution (the period of the trench interference wave (about 0.8 Hz) and the period of the mask film interference wave (about 0.1 Hz)) are clear. However, when the aperture ratio is 0.5%, only one peak (period of the mask film interference wave) clearly appears in the frequency distribution, and the period of the trench interference wave hardly appears. As a result, the trench interference period cannot be accurately detected, and the etch rate cannot be accurately calculated.

そこで、本実施の形態では、重畳干渉波から窓31によって抽出された波形を周波数解析する前に、重畳干渉波からマスク膜干渉波が占める部分の殆どを除去する。   Therefore, in this embodiment, most of the portion occupied by the mask film interference wave is removed from the superimposed interference wave before frequency analysis of the waveform extracted from the superimposed interference wave by the window 31.

図11は、開口率が0.5%の場合における重畳干渉波を示す図である。   FIG. 11 is a diagram showing a superimposed interference wave when the aperture ratio is 0.5%.

図11に示すように、開口率が0.5%の場合、図11に示すように、重畳干渉波には周期の短いトレンチ干渉波が殆ど現れないため、重畳干渉波はほぼマスク膜干渉波によって占められている。したがって、重畳干渉波を近似した波形はマスク膜干渉波の波形とほぼ等しくなる。そこで、本実施の形態では、重畳干渉波から該重畳干渉波を近似した波形を除去する。これにより、重畳干渉波からマスク膜干渉波が占める部分の殆どを除去することができる。   As shown in FIG. 11, when the aperture ratio is 0.5%, as shown in FIG. 11, since the trench interference wave having a short period does not appear in the superimposed interference wave, the superimposed interference wave is almost the mask film interference wave. Is occupied by. Therefore, the waveform approximating the superimposed interference wave is substantially equal to the waveform of the mask film interference wave. Therefore, in this embodiment, a waveform approximating the superimposed interference wave is removed from the superimposed interference wave. Thereby, most of the portion occupied by the mask film interference wave can be removed from the superimposed interference wave.

また、図11に示すように、マスク膜干渉波がほぼ占める重畳干渉波は正弦波に近く、正弦波の1/4周期以下の部分は二次多項式によって正確に近似することができる。そこで、本実施の形態では、重畳干渉波から窓31の波形を抽出する際、マスク膜干渉波の1/4周期以下の部分を抽出する。   Further, as shown in FIG. 11, the superimposed interference wave that is substantially occupied by the mask film interference wave is close to a sine wave, and the portion of the sine wave having a quarter period or less can be accurately approximated by a quadratic polynomial. Therefore, in the present embodiment, when extracting the waveform of the window 31 from the superimposed interference wave, a portion having a quarter period or less of the mask film interference wave is extracted.

すなわち、本実施の形態では、図12に示すように、マスク膜干渉波の1/4周期以下に相当するn個の窓W(k=1〜n、nは自然数)をそれぞれΔtだけずらして設定し、各窓Wの波形を抽出し、抽出した波形から該抽出した波形を二次多項式によって近似した波形を除去し、マスク膜干渉波が占める部分の殆どが除去された波形を得(図13)、除去後の波形を周波数解析する。これにより、図14に示すように、トレンチ干渉周期(約0.8Hz)のピークが明りょうに現れる周波数分布を得ることができる。なお、図14の周波数分布においてマスク膜干渉波の周期(約0.1Hz)のピークが現れないのは、窓Wの波形からマスク膜干渉波が占める部分の殆どが除去されたためである。 That is, in this embodiment, as shown in FIG. 12, n windows W k (k = 1 to n, n is a natural number) corresponding to a quarter period or less of the mask film interference wave are shifted by Δt. set Te, extracting a waveform of each window W k, extracted extract out waveform from a waveform to remove the waveform approximated by quadratic polynomial, to obtain a waveform most is removed portion where the mask film interference wave occupies (FIG. 13) The frequency of the waveform after removal is analyzed. Thereby, as shown in FIG. 14, a frequency distribution in which the peak of the trench interference period (about 0.8 Hz) clearly appears can be obtained. Incidentally, no peak appears the cycle (approximately 0.1 Hz) of the mask film interference wave in the frequency distribution of FIG. 14, because the most from the waveform of the window W k of the portion occupied by the mask film interference wave is removed.

図15は、本実施の形態に係るエッチング量算出方法を示すフローチャートである。なお、本実施の形態に係るエッチング量算出方法は、開口率が小さい場合、例えば、該開口率が0.5%を下回る場合に実行される。   FIG. 15 is a flowchart showing an etching amount calculation method according to the present embodiment. Note that the etching amount calculation method according to the present embodiment is executed when the aperture ratio is small, for example, when the aperture ratio is less than 0.5%.

図15において、まず、ステップS61乃至S64が実行され、次いで、演算部28は、重畳干渉波における、現タイミングTを終点とする1/4周期以下の部分を窓の波形として抽出する(ステップS65)(波形抽出ステップ)。   In FIG. 15, first, steps S61 to S64 are executed, and then the calculation unit 28 extracts, as a window waveform, a portion of the superimposed interference wave that is ¼ period or less with the current timing T as the end point (step S65). ) (Waveform extraction step).

次いで、演算部28は、抽出された窓の波形を二次多項式によって近似した波形(以下、単に「近似波形」という。)を算出し(ステップS151)、抽出された窓の波形から該算出された近似波形を除去し(ステップS152)(解析前処理ステップ)、マスク膜干渉波が占める部分の殆どが除去された波形を得、近似波形除去後の波形を周波数解析する(ステップS153)。   Next, the calculation unit 28 calculates a waveform (hereinafter simply referred to as “approximate waveform”) obtained by approximating the extracted window waveform by a second-order polynomial (step S151), and the calculation is performed from the extracted window waveform. The approximate waveform is removed (step S152) (pre-analysis processing step), a waveform from which most of the portion occupied by the mask film interference wave is removed is obtained, and the waveform after the approximate waveform is removed is subjected to frequency analysis (step S153).

次いで、演算部28は、ステップS67乃至S70を実行する。   Next, the calculation unit 28 executes Steps S67 to S70.

本実施の形態に係るエッチング量算出方法によれば、周波数解析に先立って抽出された窓の波形からマスク膜干渉波が占める部分の殆どが除去されるので、開口率が小さい場合、例えば、開口率が0.5%を下回る場合であっても、近似波形除去後の波形においてトレンチ干渉波が占める部分の割合を大きくすることができ、これにより、周波数解析によってトレンチ干渉周期のピークが明りょうに現れる周波数分布を得ることができる。その結果、トレンチ干渉波の周期を正確に算出することができる。なお、上述した本実施の形態に係るエッチング量算出方法は、マスク膜干渉光の波形の周期がトレンチ干渉光の干渉波の周期よりも長い場合にのみ用いることができる。   According to the etching amount calculation method according to the present embodiment, most of the portion occupied by the mask film interference wave is removed from the waveform of the window extracted prior to the frequency analysis. Even when the rate is less than 0.5%, it is possible to increase the proportion of the portion of the trench interference wave occupied in the waveform after the approximate waveform removal, and thus the peak of the trench interference period becomes clear by frequency analysis. Can be obtained. As a result, the period of the trench interference wave can be accurately calculated. Note that the etching amount calculation method according to the present embodiment described above can be used only when the period of the waveform of the mask film interference light is longer than the period of the interference wave of the trench interference light.

本実施の形態に係るエッチング量算出方法では、周波数解析に先立って抽出された窓の波形から、該抽出された窓の波形を二次多項式で近似した波形(近似波形)が除去される。開口率が小さい場合、重畳干渉波の波形はマスク膜干渉波とほぼ等しくなるので、近似波形もマスク膜干渉波とほぼ等しくなる。したがって、抽出された窓の波形からマスク膜干渉波が占める部分の殆どを確実に除去することができる。   In the etching amount calculation method according to the present embodiment, a waveform (approximate waveform) obtained by approximating the extracted window waveform with a quadratic polynomial is removed from the extracted window waveform prior to the frequency analysis. When the aperture ratio is small, the waveform of the superimposed interference wave is substantially equal to the mask film interference wave, and therefore the approximate waveform is also substantially equal to the mask film interference wave. Therefore, most of the portion occupied by the mask film interference wave can be reliably removed from the extracted waveform of the window.

また、本実施の形態に係るエッチング量算出方法では、マスク膜干渉波の1/4周期以下の部分を窓の波形として抽出する。重畳干渉波の殆どを占めるマスク膜干渉波は正弦波に近いので、マスク膜干渉波の1/4周期以下の部分を抽出すれば、該抽出された窓の波形は二次多項式によって正確に近似することができる。これにより、抽出された窓の波形からマスク膜干渉波が占める部分の殆どを正確に除去することができる。   Further, in the etching amount calculation method according to the present embodiment, a portion having a quarter period or less of the mask film interference wave is extracted as a window waveform. Since the mask film interference wave that occupies most of the superimposed interference wave is close to a sine wave, if the portion of the mask film interference wave that is less than ¼ period is extracted, the extracted waveform of the window can be accurately approximated by a quadratic polynomial. can do. Thereby, most of the portion occupied by the mask film interference wave can be accurately removed from the extracted waveform of the window.

なお、本実施の形態における開口率が小さい場合は、ウエハWの表面においてトレンチ132の開口部が占める割合が小さい場合だけでなく、ウエハWの表面においてホールの開口部が占める割合が小さい場合やトレンチ(又はホール)のアスペクト比が大きい場合も該当する。   Note that when the aperture ratio in the present embodiment is small, not only the ratio of the openings of the trenches 132 on the surface of the wafer W is small but also the ratio of the openings of the holes on the surface of the wafer W is small. This is also the case when the aspect ratio of the trench (or hole) is large.

上述した各実施の形態では、周波数解析において最大エントロピー法を用いたが、各窓内の干渉波形の数が多い場合には高速フーリエ変換法を用いてもよい。高速フーリエ変換法は最大エントロピー法よりも必要とする計算回数が少ないため、トレンチ132のエッチング量をより早く算出することができる。   In each of the above-described embodiments, the maximum entropy method is used in frequency analysis. However, when the number of interference waveforms in each window is large, a fast Fourier transform method may be used. Since the fast Fourier transform method requires fewer calculations than the maximum entropy method, the etching amount of the trench 132 can be calculated earlier.

また、上述した各実施の形態に係るエッチング量算出方法を用いて或るトレンチのエッチング量(エッチング深さ)を算出する場合において、マスク膜131がレーザ光を透過させるような膜であるとき、図8に示すように、算出されたエッチング量(図中「モニター深さ」)と実測されたエッチング量(図中「エッチング深さ」)との間に誤差が生じる場合がある。これは、重畳干渉光が主として反射光L及び反射光Lの干渉光ではなく、反射光L及び反射光Lの干渉光を含み、反射光Lの光路長変化にマスク膜131の厚さ変化だけでなくマスク膜131の屈折率も影響するためであると考えられる。 Further, when the etching amount (etching depth) of a certain trench is calculated using the etching amount calculation method according to each of the above-described embodiments, when the mask film 131 is a film that transmits laser light, As shown in FIG. 8, an error may occur between the calculated etching amount (“monitor depth” in the figure) and the actually measured etching amount (“etching depth” in the figure). This is because the superimposed interference light mainly includes the interference light of the reflected light L 3 and the reflected light L 4 , not the interference light of the reflected light L 2 and the reflected light L 4 , and the mask film 131 is affected by the change in the optical path length of the reflected light L 3. This is probably because not only the thickness change but also the refractive index of the mask film 131 affects.

算出されたエッチング量と実測されたエッチング量との間に誤差が生じる場合には、エッチング量算出に先立ち、試験用のウエハWを用いてトレンチ132のエッチング量(エッチングレート)を実測するとともに、上述した各実施の形態に係るエッチング量算出方法によってトレンチ132のエッチング量を算出し、実測されたエッチング量と算出されたエッチング量との回帰式等を求めるのが好ましい。そして、以後のエッチングにおいて、上述した各実施の形態に係るエッチング量算出方法によってトレンチ132のエッチング量を算出した後、該算出されたエッチング量を回帰式によって補正すればよい。   When an error occurs between the calculated etching amount and the actually measured etching amount, the etching amount (etching rate) of the trench 132 is measured using the test wafer W prior to the etching amount calculation, and It is preferable to calculate the etching amount of the trench 132 by the etching amount calculation method according to each of the above-described embodiments, and obtain a regression equation between the actually measured etching amount and the calculated etching amount. In the subsequent etching, after the etching amount of the trench 132 is calculated by the etching amount calculation method according to each of the embodiments described above, the calculated etching amount may be corrected by a regression equation.

上述した各実施の形態では、トレンチ132のエッチング量が算出されたが、図6、図7や図15のエッチング量算出方法を実行してホールのエッチング量を算出してもよい。   In each of the above-described embodiments, the etching amount of the trench 132 is calculated. However, the etching amount of holes may be calculated by executing the etching amount calculation method shown in FIGS.

また、本発明の目的は、上述した各実施の形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体を、コンピュータ(例えば、コントローラ29)に供給し、コンピュータのCPUが記憶媒体に格納されたプログラムコードを読み出して実行することによっても達成される。   Another object of the present invention is to supply a computer (for example, the controller 29) with a storage medium storing software program codes for realizing the functions of the above-described embodiments, and the CPU of the computer is stored in the storage medium. It is also achieved by reading and executing the program code.

この場合、記憶媒体から読み出されたプログラムコード自体が上述した各実施の形態の機能を実現することになり、プログラムコード及びそのプログラムコードを記憶した記憶媒体は本発明を構成することになる。   In this case, the program code itself read from the storage medium realizes the functions of the above-described embodiments, and the program code and the storage medium storing the program code constitute the present invention.

また、プログラムコードを供給するための記憶媒体としては、例えば、RAM、NV−RAM、フロッピー(登録商標)ディスク、ハードディスク、光磁気ディスク、CD−ROM、CD−R、CD−RW、DVD(DVD−ROM、DVD−RAM、DVD−RW、DVD+RW)等の光ディスク、磁気テープ、不揮発性のメモリカード、他のROM等の上記プログラムコードを記憶できるものであればよい。或いは、上記プログラムコードは、インターネット、商用ネットワーク、若しくはローカルエリアネットワーク等に接続される不図示の他のコンピュータやデータベース等からダウンロードすることによりコンピュータに供給されてもよい。   Examples of the storage medium for supplying the program code include RAM, NV-RAM, floppy (registered trademark) disk, hard disk, magneto-optical disk, CD-ROM, CD-R, CD-RW, DVD (DVD). -ROM, DVD-RAM, DVD-RW, DVD + RW) and other optical disks, magnetic tapes, non-volatile memory cards, other ROMs, etc., as long as they can store the program code. Alternatively, the program code may be supplied to the computer by downloading from another computer or database (not shown) connected to the Internet, a commercial network, a local area network, or the like.

また、コンピュータが読み出したプログラムコードを実行することにより、上記各実施の形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、CPU上で稼動しているOS(オペレーティングシステム)等が実際の処理の一部又は全部を行い、その処理によって上述した各実施の形態の機能が実現される場合も含まれる。   Further, by executing the program code read by the computer, not only the functions of the above-described embodiments are realized, but also an OS (operating system) running on the CPU based on the instruction of the program code. Includes a case where part or all of the actual processing is performed and the functions of the above-described embodiments are realized by the processing.

更に、記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書き込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるCPU等が実際の処理の一部又は全部を行い、その処理によって上述した各実施の形態の機能が実現される場合も含まれる。   Further, after the program code read from the storage medium is written in a memory provided in a function expansion board inserted into the computer or a function expansion unit connected to the computer, the function expansion is performed based on the instruction of the program code. This includes a case where the CPU or the like provided in the board or the function expansion unit performs part or all of the actual processing, and the functions of the above-described embodiments are realized by the processing.

上記プログラムコードの形態は、オブジェクトコード、インタプリタにより実行されるプログラムコード、OSに供給されるスクリプトデータ等の形態から成ってもよい。   The form of the program code may include an object code, a program code executed by an interpreter, script data supplied to the OS, and the like.

次に、本発明の実施例について説明する。   Next, examples of the present invention will be described.

実施例1
まず、シリコンからなる被エッチング層130上に酸化膜からなるレーザ光L1を透過させるマスク膜131が形成されたウエハWを準備し、基板処理装置10でエッチングにより被エッチング層130にディープトレンチ132を形成した。このときのエッチング条件は以下の通りであった。
実際のエッチングレート : 1200nm/分
選択比 : 10対1(被エッチング層130対マスク膜131)
開口比 : 0.05
測定波長(レーザ光L1の波長) : 300nm
サンプリングレート : 10Hz
ディープトレンチ132のエッチングにおいて、窓の始点から終点までの時間を30秒に設定し、図6のエッチング量算出方法を実行して各タイミングにおける積算平均されたトレンチ干渉周期を求め、該積算平均されたトレンチ干渉周期から各タイミングにおけるディープトレンチ132のエッチングレートを求めた。そして、求められたエッチングレートをグラフに示した(図16参照。)。
Example 1
First, a wafer W on which a mask film 131 that transmits a laser beam L1 made of an oxide film is formed on a layer to be etched 130 made of silicon is prepared, and a deep trench 132 is formed in the layer to be etched 130 by etching with the substrate processing apparatus 10. Formed. The etching conditions at this time were as follows.
Actual etching rate: 1200 nm / min Selection ratio: 10 to 1 (etched layer 130 to mask film 131)
Aperture ratio: 0.05
Measurement wavelength (wavelength of laser light L1): 300 nm
Sampling rate: 10Hz
In the etching of the deep trench 132, the time from the start point to the end point of the window is set to 30 seconds, and the etching amount calculation method of FIG. 6 is executed to obtain the integrated and averaged trench interference period at each timing. The etching rate of the deep trench 132 at each timing was obtained from the trench interference period. Then, the obtained etching rate is shown in a graph (see FIG. 16).

比較例1
また、上述したディープトレンチ132のエッチングにおいて、検出器27によってウエハWからの重畳干渉波を観測し、該重畳干渉波における短周期の干渉波を読み取り、該短周期の干渉波における各極値間の時間からトレンチ干渉周期を求め、該トレンチ干渉周期から各極値間におけるディープトレンチ132のエッチングレートを求めた。そして、求められたエッチングレートをグラフに示した(図16参照。)。
Comparative Example 1
Further, in the above-described etching of the deep trench 132, a superposition interference wave from the wafer W is observed by the detector 27, a short period interference wave in the superposition interference wave is read, and between each extreme value in the short period interference wave. From this time, the trench interference period was obtained, and the etching rate of the deep trench 132 between each extreme value was obtained from the trench interference period. Then, the obtained etching rate is shown in a graph (see FIG. 16).

図16のグラフより、実施例1のエッチングレートの方が、比較例1のエッチングレートよりも変動が小さく安定していることが分かった。   From the graph of FIG. 16, it was found that the etching rate of Example 1 was smaller and more stable than the etching rate of Comparative Example 1.

さらに、実施例1のエッチング量における誤差と比較例1のエッチング量における誤差を時系列的にグラフへ示した(図17参照。)。   Further, the error in the etching amount of Example 1 and the error in the etching amount of Comparative Example 1 are shown in a graph in time series (see FIG. 17).

図17のグラフより、実施例1のエッチング量の方が、比較例1のエッチング量よりも誤差が小さいことが分かった。これにより、図6のエッチング量算出方法はエッチング量の算出を正確に行えることが分かった。   From the graph of FIG. 17, it was found that the etching amount of Example 1 had a smaller error than the etching amount of Comparative Example 1. Accordingly, it was found that the etching amount calculation method of FIG. 6 can accurately calculate the etching amount.

実施例2
次に、シリコンからなる被エッチング層130上に酸化膜からなるマスク膜131が形成されたウエハWを準備し、基板処理装置10でエッチングにより被エッチング層130にシャロートレンチ132を形成した。このときのエッチング条件は以下の通りであった。
実際のエッチングレート : 360nm/分
選択比 : 10対1(被エッチング層130対マスク膜131)
開口比 : 0.2
測定波長(レーザ光L1の波長) : 300nm
サンプリングレート : 10Hz
シャロートレンチ132のエッチングにおいて、窓の始点から終点までの時間を25秒
に設定し、図6のエッチング量算出方法を実行して各タイミングにおけるシャロートレンチ132のエッチング量(エッチング深さ)を算出した。そして、算出されたエッチング量をグラフに示した(図18参照。)。
Example 2
Next, a wafer W having a mask film 131 made of an oxide film formed on the layer to be etched 130 made of silicon was prepared, and the shallow trench 132 was formed in the layer to be etched 130 by etching with the substrate processing apparatus 10. The etching conditions at this time were as follows.
Actual etching rate: 360 nm / min Selection ratio: 10 to 1 (etched layer 130 to mask film 131)
Aperture ratio: 0.2
Measurement wavelength (wavelength of laser light L1): 300 nm
Sampling rate: 10Hz
In the etching of the shallow trench 132, the time from the start point to the end point of the window was set to 25 seconds, and the etching amount calculation method of FIG. 6 was executed to calculate the etching amount (etching depth) of the shallow trench 132 at each timing. . Then, the calculated etching amount is shown in a graph (see FIG. 18).

比較例2
また、上述したシャロートレンチ132のエッチングにおいて、検出器27によってウエハWからの重畳干渉波を観測し、エッチング開始時から各タイミングまでの全ての重畳干渉波から周波数解析によって周波数分布を得て、該周波数分布に基づいてエッチング開始時から各タイミングまでの間におけるトレンチ干渉周期を求め、該トレンチ干渉周期から各タイミングにおけるシャロートレンチ132のエッチング量(エッチング深さ)を算出した。すなわち、図3に示す窓を用いずに重畳干渉波からエッチング量を算出した。そして、算出されたエッチング量をグラフに示した(図18参照。)。
Comparative Example 2
Further, in the above-described etching of the shallow trench 132, the superimposed interference wave from the wafer W is observed by the detector 27, the frequency distribution is obtained by frequency analysis from all the superimposed interference waves from the start of etching to each timing, Based on the frequency distribution, the trench interference period between the start of etching and each timing was obtained, and the etching amount (etching depth) of the shallow trench 132 at each timing was calculated from the trench interference period. That is, the etching amount was calculated from the superimposed interference wave without using the window shown in FIG. Then, the calculated etching amount is shown in a graph (see FIG. 18).

図18のグラフでは、比較例2のエッチング量のデータが乱れているが、これは重畳干渉波に外乱が加わったためと考えられる。一方、実施例2のエッチング量のデータは乱れていなかった。実施例2と比較例2は乱れが加わった同じ重畳干渉波から求められたエッチング量であるため、これにより、図6のエッチング量算出方法は重畳干渉波に外乱が加わってもエッチング量の算出を安定して正確に行うことができることが分かった。   In the graph of FIG. 18, the etching amount data of Comparative Example 2 is disturbed. This is considered to be because a disturbance is added to the superimposed interference wave. On the other hand, the etching amount data of Example 2 was not disturbed. Since Example 2 and Comparative Example 2 are etching amounts obtained from the same superimposed interference wave to which disturbance is added, the etching amount calculation method in FIG. 6 calculates the etching amount even when disturbance is added to the superimposed interference wave. Has been found to be stable and accurate.

実施例3
上述した実施例1,2とは他のウエハWにおけるトレンチ132のエッチングにおいて、図6のエッチング量算出方法を実行して各タイミングにおけるトレンチ132のエッチングレートを求めた。そして、求められたエッチングレートをグラフに示した(図19参照。)。
Example 3
In the etching of the trench 132 on the other wafer W from the first and second embodiments described above, the etching amount calculation method of FIG. 6 was executed to obtain the etching rate of the trench 132 at each timing. Then, the obtained etching rate is shown in a graph (see FIG. 19).

比較例3
また、実施例3と同じエッチングにおいて、最大エントロピー法ではなく高速フーリエ変換法を用いた以外は図6のエッチング量算出方法と同じ条件のエッチング量算出方法を実行して各タイミングにおけるトレンチ132のエッチングレートを求めた。そして、求められたエッチングレートをグラフに示した(図19参照。)。
Comparative Example 3
Further, in the same etching as in the third embodiment, the etching amount calculation method under the same conditions as the etching amount calculation method in FIG. 6 is executed except that the fast entropy method is used instead of the maximum entropy method, and the trench 132 is etched at each timing. I asked for the rate. Then, the obtained etching rate is shown in a graph (see FIG. 19).

図19のグラフより、実施例3のエッチングレートの方が、比較例3のエッチングレートよりも変動が小さく安定していることが分かった。これにより、最大エントロピー法を用いるとエッチング量の算出を安定して行えることが分かった。   From the graph of FIG. 19, it was found that the etching rate of Example 3 was smaller and more stable than the etching rate of Comparative Example 3. Thus, it was found that the etching amount can be calculated stably when the maximum entropy method is used.

実施例4
まず、開口率が5%のウエハWと開口率が0.5%のウエハWを準備し、各ウエハWの被エッチング層130をエッチングする際に、図15のエッチング量算出方法を実行して各エッチングレートを求めた。そして、求められた各エッチングレートをグラフに示した(図20参照。)。
Example 4
First, a wafer W with an aperture ratio of 5% and a wafer W with an aperture ratio of 0.5% are prepared, and the etching amount calculation method shown in FIG. Each etching rate was determined. The obtained etching rates are shown in a graph (see FIG. 20).

比較例4
実施例4と同様に、開口率が5%のウエハWと開口率が0.5%のウエハWを準備し、各ウエハWの被エッチング層130をエッチングする際に、図6のエッチング量算出方法を実行して各エッチングレートを求めた。そして、求められた各エッチングレートをグラフに示した(図21参照。)。
Comparative Example 4
Similar to the fourth embodiment, when the wafer W having an aperture ratio of 5% and the wafer W having an aperture ratio of 0.5% are prepared and the etching target layer 130 of each wafer W is etched, the etching amount calculation of FIG. The method was carried out to determine each etching rate. The obtained etching rates are shown in a graph (see FIG. 21).

図20及び図21のグラフを比較すると、図6のエッチング量算出方法を実行した場合、開口率が5%におけるエッチングレートは安定しているものの、開口率が0.5%におけるエッチングレートは安定していないのに対して、図15のエッチング量算出方法を実行した場合、開口率が5%におけるエッチングレート及び開口率が0.5%におけるエッチングレートのいずれも安定していることが分かった。これにより、抽出された窓の波形を周波数解析する前に該抽出された窓の波形から該窓の波形の近似波形を除去すると、開口率が小さい場合でも正確にエッチングレートを求めることができることが分かった。   Comparing the graphs of FIGS. 20 and 21, when the etching amount calculation method of FIG. 6 is executed, the etching rate is stable when the aperture ratio is 5%, but the etching rate is stable when the aperture ratio is 0.5%. On the other hand, when the etching amount calculation method of FIG. 15 was executed, it was found that both the etching rate at an aperture ratio of 5% and the etching rate at an aperture ratio of 0.5% were stable. . Thus, when the approximate waveform of the window waveform is removed from the extracted window waveform before frequency analysis of the extracted window waveform, the etching rate can be accurately obtained even when the aperture ratio is small. I understood.

本発明の第1の実施の形態に係るエッチング量算出方法が適用される基板処理装置の構成を概略的に示す断面図である。1 is a cross-sectional view schematically showing a configuration of a substrate processing apparatus to which an etching amount calculation method according to a first embodiment of the present invention is applied. トレンチのエッチングにおけるエッチングレートの低下を説明するための図である。It is a figure for demonstrating the fall of the etching rate in the etching of a trench. 本実施の形態に係るエッチング量算出方法における重畳干渉波からの所定期間の波形の抽出を説明するための図である。It is a figure for demonstrating extraction of the waveform of the predetermined period from the superimposed interference wave in the etching amount calculation method which concerns on this Embodiment. 図3における窓の波形から最大エントロピー法を用いた周波数解析によって得られた周波数分布を示す図である。It is a figure which shows the frequency distribution obtained by the frequency analysis using the maximum entropy method from the waveform of the window in FIG. 本実施の形態に係るエッチング量算出方法におけるエッチングレートの平均値の算出方法を説明するための図である。It is a figure for demonstrating the calculation method of the average value of the etching rate in the etching amount calculation method which concerns on this Embodiment. 本実施の形態に係るエッチング量算出方法を示すフローチャートである。It is a flowchart which shows the etching amount calculation method which concerns on this Embodiment. 本発明の第2の実施の形態に係るエッチング量算出方法を示すフローチャートである。It is a flowchart which shows the etching amount calculation method which concerns on the 2nd Embodiment of this invention. 算出されたエッチング量と実測されたエッチング量との間における誤差を示す図である。It is a figure which shows the difference | error between the calculated etching amount and the measured etching amount. 開口率が変化した場合における重畳干渉波の一部の変化を示す図である。It is a figure which shows a part change of the superposition interference wave when an aperture ratio changes. 開口率が変化した場合における窓の波形の周波数分布の変化を示す図である。It is a figure which shows the change of the frequency distribution of the waveform of the window when an aperture ratio changes. 開口率が0.5%の場合における重畳干渉波を示す図である。It is a figure which shows the superimposed interference wave in case an aperture ratio is 0.5%. 本発明の第3の実施の形態に係るエッチング量算出方法における窓の波形の抽出を説明するための図である。It is a figure for demonstrating extraction of the waveform of the window in the etching amount calculation method which concerns on the 3rd Embodiment of this invention. 重畳干渉波からマスク膜干渉波が占める部分の殆どが除去された波形を示す図である。It is a figure which shows the waveform from which most of the part which a mask film | membrane interference wave occupies was removed from the superimposition interference wave. 重畳干渉波からマスク膜干渉波が占める部分の殆どが除去された波形から周波数解析によって得られた周波数分布を示す図である。It is a figure which shows the frequency distribution obtained by the frequency analysis from the waveform from which most of the part which a mask film | membrane interference wave occupies was removed from the superimposition interference wave. 本実施の形態に係るエッチング量算出方法を示すフローチャートである。It is a flowchart which shows the etching amount calculation method which concerns on this Embodiment. 図6のエッチング量算出方法によって算出されたエッチングレートと干渉波における各極値間の時間から求められたエッチングレートとの比較図である。FIG. 7 is a comparison diagram of the etching rate calculated by the etching amount calculation method of FIG. 6 and the etching rate obtained from the time between each extreme value in the interference wave. 図6のエッチング量算出方法によって算出されたエッチング量及び実際のエッチング量の誤差と、干渉波における各極値間の時間から求められたエッチング量及び実際のエッチング量の誤差との比較図である。FIG. 7 is a comparison diagram of an error between an etching amount and an actual etching amount calculated by the etching amount calculating method in FIG. 6, and an etching amount and an actual etching amount error obtained from the time between each extreme value in an interference wave . 図6のエッチング量算出方法によって算出されたエッチング量と、エッチング開始時から各タイミングまでの重畳干渉波の周波数解析によって求められたエッチング量との比較図である。FIG. 7 is a comparison diagram between an etching amount calculated by the etching amount calculation method of FIG. 6 and an etching amount obtained by frequency analysis of a superimposed interference wave from the start of etching to each timing. 最大エントロピー法を用いて算出されたエッチングレートと高速フーリエ変換法を用いて算出されたエッチングレートとの比較図である。It is a comparison figure of the etching rate computed using the maximum entropy method, and the etching rate computed using the fast Fourier transform method. 図15のエッチング量算出方法を用いて得られた開口率が異なる各ウエハのエッチングレートを示す図である。It is a figure which shows the etching rate of each wafer from which the aperture ratios obtained using the etching amount calculation method of FIG. 15 differ. 図6のエッチング量算出方法を用いて得られた開口率が異なる各ウエハのエッチングレートを示す図である。It is a figure which shows the etching rate of each wafer from which the aperture ratios obtained using the etching amount calculation method of FIG. 6 differ. エッチング中における光の干渉を説明するための図である。It is a figure for demonstrating the interference of the light during an etching. 重畳干渉波を示す図である。It is a figure which shows a superimposed interference wave. 干渉波における各極値間の時間から求められたエッチングレートを示す図である。It is a figure which shows the etching rate calculated | required from the time between each extreme value in an interference wave. 外乱が加わった重畳干渉波を示す図である。It is a figure which shows the superimposed interference wave to which the disturbance was added. 重畳干渉波の周波数解析によって得られる周波数分布であり、図26(A)は重畳干渉波に外乱が加わっていない場合であり、図26(B)は重畳干渉波に外乱が加わっている場合である。FIG. 26A shows a frequency distribution obtained by frequency analysis of the superimposed interference wave. FIG. 26A shows a case where no disturbance is added to the superimposed interference wave, and FIG. 26B shows a case where disturbance is added to the superimposed interference wave. is there.

符号の説明Explanation of symbols

L1 レーザ光
、L、L 反射光
W ウエハ
10 基板処理装置
25 エッチング量算出装置
26 レーザ光源
27 検出器
28 演算部
30、50 重畳干渉波
31 窓
130 被エッチング層
131 マスク膜
132 トレンチ
L1 Laser light L 2 , L 3 , L 4 reflected light W Wafer 10 Substrate processing device 25 Etching amount calculation device 26 Laser light source 27 Detector 28 Calculation unit 30, 50 Superimposed interference wave 31 Window 130 Etched layer 131 Mask film 132 Trench

Claims (13)

マスク膜を用いて凹部を形成する基板のエッチングにおいて前記凹部のエッチング量を算出するエッチング量算出方法であって、
前記基板に光を照射する照射ステップと、
少なくとも前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉光が他の干渉光に重畳された重畳干渉光を受光する受光ステップと、
前記受光された重畳干渉光から重畳干渉波を算出する干渉波算出ステップと、
前記重畳干渉波から所定期間の波形を抽出する波形抽出ステップと、
前記抽出された波形に周波数解析を施す周波数解析ステップと、
前記周波数解析によって得られた周波数分布から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出する干渉周期検出ステップと、
前記所定期間を所定時間だけずらしつつ前記干渉波算出ステップ、前記波形抽出ステップ、前記周波数解析ステップ及び前記干渉周期検出ステップを繰り返し、繰り返しの度に前記検出された干渉波の周期を積算平均する積算平均ステップと、
前記積算平均された干渉波の周期に基づいて前記凹部のエッチング量を算出するエッチング量算出ステップとを有することを特徴とするエッチング量算出方法。
An etching amount calculation method for calculating an etching amount of the recess in etching of a substrate that forms a recess using a mask film,
An irradiation step of irradiating the substrate with light;
A light receiving step of receiving a superimposed interference light in which interference light of at least reflected light from the mask film and reflected light from the bottom of the recess is superimposed on other interference light;
An interference wave calculating step of calculating a superimposed interference wave from the received superimposed interference light;
A waveform extraction step of extracting a waveform of a predetermined period from the superimposed interference wave;
A frequency analysis step of performing frequency analysis on the extracted waveform;
An interference period detection step of detecting a period of interference waves of reflected light from the mask film and reflected light from the bottom of the recess from the frequency distribution obtained by the frequency analysis;
An integration in which the interference wave calculation step, the waveform extraction step, the frequency analysis step, and the interference cycle detection step are repeated while shifting the predetermined period by a predetermined time, and the cycle of the detected interference wave is integrated and averaged each time it is repeated. Average step and
An etching amount calculating method, comprising: an etching amount calculating step of calculating an etching amount of the concave portion based on a cycle of the integrated and averaged interference wave.
前記所定期間は、前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波よりも周期が長い前記他の干渉光の波形の1周期よりも大きいことを特徴とする請求項1記載のエッチング量算出方法。   2. The predetermined period is longer than one period of a waveform of the other interference light having a longer period than an interference wave of reflected light from the mask film and reflected light from the bottom of the recess. The etching amount calculation method described. 前記他の干渉光の波形の周期が前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期よりも長い場合に、前記周波数解析ステップに先立って、前記重畳干渉波から抽出された所定期間の波形から前記他の干渉光の波形が占める部分の殆どを除去する解析前処理ステップをさらに有し、
前記周波数解析ステップでは、前記他の干渉光の波形が占める部分の殆どが除去された波形に周波数解析を施すことを特徴とする請求項1又は2記載のエッチング量算出方法。
When the period of the waveform of the other interference light is longer than the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess, prior to the frequency analysis step, from the superimposed interference wave An analysis pre-processing step for removing most of the portion occupied by the waveform of the other interference light from the extracted waveform of a predetermined period;
3. The etching amount calculation method according to claim 1, wherein in the frequency analysis step, frequency analysis is performed on a waveform from which most of the portion occupied by the waveform of the other interference light is removed.
前記解析前処理ステップでは、前記抽出された波形から、該抽出された波形を二次多項式で近似した波形を除去することを特徴とする請求項3記載のエッチング量算出方法。   4. The etching amount calculation method according to claim 3, wherein, in the pre-analysis processing step, a waveform obtained by approximating the extracted waveform with a second-order polynomial is removed from the extracted waveform. 前記所定期間は前記他の干渉光の波形の1/4周期以下であることを特徴とする請求項3又は4記載のエッチング量算出方法。   The etching amount calculation method according to claim 3 or 4, wherein the predetermined period is equal to or less than ¼ period of the waveform of the other interference light. 前記基板の表面における前記凹部の開口率が0.5%以下であるか、又は前記凹部がディープトレンチであることを特徴とする請求項3乃至5のいずれか1項に記載のエッチング量算出方法。   6. The etching amount calculation method according to claim 3, wherein an opening ratio of the recesses on the surface of the substrate is 0.5% or less, or the recesses are deep trenches. . 前記周波数解析では最大エントロピー法を用いることを特徴とする請求項1乃至6のいずれか1項に記載のエッチング量算出方法。   The etching amount calculation method according to claim 1, wherein a maximum entropy method is used in the frequency analysis. 前記周波数分布から検出された前記干渉波の周期が異常値に該当する場合には該干渉波の周期を除去する干渉周期修正ステップをさらに有することを特徴とする請求項1乃至7のいずれか1項に記載のエッチング量算出方法。   8. The method according to claim 1, further comprising an interference period correcting step of removing the period of the interference wave when the period of the interference wave detected from the frequency distribution corresponds to an abnormal value. The etching amount calculation method according to item. 前記干渉周期修正ステップでは、前記異常値に該当する前記干渉波の周期が求められた前記所定期間の前の前記所定期間又は後の前記所定期間から求められる前記干渉波の周期を、前記異常値に該当する前記干渉波の周期が求められた前記所定期間の干渉波の周期とみなすことを特徴とする請求項8記載のエッチング量算出方法。   In the interference period correcting step, the period of the interference wave obtained from the predetermined period before or after the predetermined period from which the period of the interference wave corresponding to the abnormal value is obtained is set to the abnormal value. The etching amount calculation method according to claim 8, wherein the period of the interference wave corresponding to is considered as the period of the interference wave in the predetermined period obtained. 前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を予め予測し、
前記干渉周期検出ステップでは、前記周波数解析によって得られた周波数分布において、前記予測された周期近傍から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出することを特徴とする請求項1乃至9のいずれか1項に記載のエッチング量算出方法。
Predicting the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess,
In the interference period detection step, in the frequency distribution obtained by the frequency analysis, a period of interference waves of reflected light from the mask film and reflected light from the bottom of the recess is detected from the vicinity of the predicted period. The etching amount calculation method according to claim 1, wherein:
前記他の干渉光は、前記マスク膜の表面からの反射光、並びに前記マスク膜及び前記基板の表面の境界面からの反射光の干渉光であることを特徴とする請求項1乃至10のいずれか1項に記載のエッチング量算出方法。   11. The other interference light is reflected light from a surface of the mask film and interference light of reflected light from a boundary surface between the mask film and the surface of the substrate. The etching amount calculation method according to claim 1. マスク膜を用いて凹部を形成する基板のエッチングにおいて前記凹部のエッチング量を算出するエッチング量算出方法をコンピュータに実行させるプログラムを格納するコンピュータで読み取り可能な記憶媒体であって、前記エッチング量算出方法は、
前記基板に光を照射する照射ステップと、
少なくとも前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉光が他の干渉光に重畳された重畳干渉光を受光する受光ステップと、
前記受光された重畳干渉光から重畳干渉波を算出する干渉波算出ステップと、
前記重畳干渉波から所定期間の波形を抽出する波形抽出ステップと、
前記抽出された波形に周波数解析を施す周波数解析ステップと、
前記周波数解析によって得られた周波数分布から前記マスク膜からの反射光及び前記凹
部の底からの反射光の干渉波の周期を検出する干渉周期検出ステップと、
前記所定期間を所定時間だけずらしつつ前記干渉波算出ステップ、前記波形抽出ステップ、前記周波数解析ステップ及び前記干渉周期検出ステップを繰り返し、繰り返しの度に前記検出された干渉波の周期を積算平均する積算平均ステップと、
前記積算平均された干渉波の周期に基づいて前記凹部のエッチング量を算出するエッチング量算出ステップとを有することを特徴とする記憶媒体。
A computer-readable storage medium storing a program for causing a computer to execute an etching amount calculation method for calculating an etching amount of the recess in etching a substrate on which a recess is formed using a mask film, the etching amount calculating method Is
An irradiation step of irradiating the substrate with light;
A light receiving step of receiving a superimposed interference light in which interference light of at least reflected light from the mask film and reflected light from the bottom of the recess is superimposed on other interference light;
An interference wave calculating step of calculating a superimposed interference wave from the received superimposed interference light;
A waveform extraction step of extracting a waveform of a predetermined period from the superimposed interference wave;
A frequency analysis step of performing frequency analysis on the extracted waveform;
An interference period detection step of detecting a period of interference waves of reflected light from the mask film and reflected light from the bottom of the recess from the frequency distribution obtained by the frequency analysis;
An integration in which the interference wave calculation step, the waveform extraction step, the frequency analysis step, and the interference cycle detection step are repeated while shifting the predetermined period by a predetermined time, and the cycle of the detected interference wave is integrated and averaged each time it is repeated. Average step and
An etching amount calculating step of calculating an etching amount of the concave portion based on the cycle of the integrated averaged interference wave.
マスク膜を用いて凹部を形成する基板のエッチングにおいて前記凹部のエッチング量を算出するエッチング量算出装置において、
前記基板に光を照射する照射部と、
少なくとも前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉光が他の干渉光に重畳された重畳干渉光を受光する受光部と、
前記受光された重畳干渉光から重畳干渉波を算出する干渉波算出部と、
前記重畳干渉波から所定期間の波形を抽出する波形抽出部と、
前記抽出された波形に周波数解析を施す周波数解析部と、
前記周波数解析によって得られた周波数分布から前記マスク膜からの反射光及び前記凹部の底からの反射光の干渉波の周期を検出する干渉周期検出部と、
前記所定期間を所定時間だけずらしつつ前記重畳干渉波の算出、前記所定期間の波形の抽出、前記周波数解析及び前記干渉波の周期の検出を繰り返し、繰り返しの度に前記検出された干渉波の周期を積算平均する積算平均部と、
前記積算平均された干渉波の周期に基づいて前記凹部のエッチング量を算出するエッチング量算出部とを備えることを特徴とするエッチング量算出装置。
In an etching amount calculation apparatus for calculating an etching amount of the recess in etching of a substrate that forms a recess using a mask film,
An irradiation unit for irradiating the substrate with light;
A light receiving unit that receives at least the interference light of the reflected light from the mask film and the reflected light of the reflected light from the bottom of the concave portion superimposed on the other interference light;
An interference wave calculator for calculating a superimposed interference wave from the received superimposed interference light;
A waveform extraction unit that extracts a waveform of a predetermined period from the superimposed interference wave;
A frequency analysis unit for performing frequency analysis on the extracted waveform;
An interference period detection unit that detects the period of the interference wave of the reflected light from the mask film and the reflected light from the bottom of the recess from the frequency distribution obtained by the frequency analysis;
The calculation of the superimposed interference wave, the extraction of the waveform during the predetermined period, the frequency analysis, and the detection of the period of the interference wave are repeated while shifting the predetermined period by a predetermined time, and the period of the detected interference wave is repeated each time it is repeated. A cumulative average part for cumulative average;
An etching amount calculation apparatus comprising: an etching amount calculation unit that calculates an etching amount of the concave portion based on a cycle of the integrated and averaged interference wave.
JP2008192311A 2008-01-17 2008-07-25 Etching amount calculation method, storage medium, and etching amount calculation device Active JP5026363B2 (en)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP2008192311A JP5026363B2 (en) 2008-01-17 2008-07-25 Etching amount calculation method, storage medium, and etching amount calculation device
CN2008101888359A CN101488456B (en) 2008-01-17 2008-12-26 Etching amount calculating method and etching amount calculating apparatus
KR1020090001096A KR101073229B1 (en) 2008-01-17 2009-01-07 Etching amount calculation method, storage medium, and etching amount calculation apparatus
US12/353,457 US8014891B2 (en) 2008-01-17 2009-01-14 Etching amount calculating method, storage medium, and etching amount calculating apparatus
TW098101548A TWI508161B (en) 2008-01-17 2009-01-16 An etching amount calculating method, a memory medium, and an etching amount calculating means

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2008008312 2008-01-17
JP2008008312 2008-01-17
JP2008192311A JP5026363B2 (en) 2008-01-17 2008-07-25 Etching amount calculation method, storage medium, and etching amount calculation device

Publications (2)

Publication Number Publication Date
JP2009194356A true JP2009194356A (en) 2009-08-27
JP5026363B2 JP5026363B2 (en) 2012-09-12

Family

ID=40891270

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008192311A Active JP5026363B2 (en) 2008-01-17 2008-07-25 Etching amount calculation method, storage medium, and etching amount calculation device

Country Status (3)

Country Link
JP (1) JP5026363B2 (en)
CN (1) CN101488456B (en)
TW (1) TWI508161B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2011016525A1 (en) * 2009-08-06 2013-01-17 芝浦メカトロニクス株式会社 Plasma etching apparatus and plasma etching method
US8872059B2 (en) 2010-09-06 2014-10-28 Samsung Electronics Co., Ltd. Etching system and method of controlling etching process condition
JP2018088455A (en) * 2016-11-28 2018-06-07 住友電気工業株式会社 Method of manufacturing semiconductor device

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI582837B (en) 2012-06-11 2017-05-11 應用材料股份有限公司 Melt depth determination using infrared interferometric technique in pulsed laser annealing
CN103745904B (en) * 2013-12-31 2016-08-17 深圳市华星光电技术有限公司 A kind of dry etching machine and lithographic method thereof
CN113496887B (en) * 2020-04-03 2023-06-02 重庆超硅半导体有限公司 Uniform etching method of silicon wafer for integrated circuit

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61124997A (en) * 1984-11-22 1986-06-12 富士通株式会社 Signal analyzer
JPH0271517A (en) * 1988-09-06 1990-03-12 Fujitsu Ltd Etching apparatus
JPH09129619A (en) * 1995-08-31 1997-05-16 Toshiba Corp Etching depth measuring device
JP2000329525A (en) * 1999-05-18 2000-11-30 Toshiba Corp Measurement method of difference in level, and measuring method of etching depth, and devices therefor
JP2001343219A (en) * 2000-05-31 2001-12-14 Toshiba Corp Apparatus for measuring etching depth and method therefor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10204943B4 (en) * 2002-02-07 2005-04-21 Leica Microsystems Jena Gmbh Method for determining layer thicknesses

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61124997A (en) * 1984-11-22 1986-06-12 富士通株式会社 Signal analyzer
JPH0271517A (en) * 1988-09-06 1990-03-12 Fujitsu Ltd Etching apparatus
JPH09129619A (en) * 1995-08-31 1997-05-16 Toshiba Corp Etching depth measuring device
JP2000329525A (en) * 1999-05-18 2000-11-30 Toshiba Corp Measurement method of difference in level, and measuring method of etching depth, and devices therefor
JP2001343219A (en) * 2000-05-31 2001-12-14 Toshiba Corp Apparatus for measuring etching depth and method therefor

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPWO2011016525A1 (en) * 2009-08-06 2013-01-17 芝浦メカトロニクス株式会社 Plasma etching apparatus and plasma etching method
JP5665746B2 (en) * 2009-08-06 2015-02-04 芝浦メカトロニクス株式会社 Plasma etching apparatus and plasma etching method
US8872059B2 (en) 2010-09-06 2014-10-28 Samsung Electronics Co., Ltd. Etching system and method of controlling etching process condition
JP2018088455A (en) * 2016-11-28 2018-06-07 住友電気工業株式会社 Method of manufacturing semiconductor device

Also Published As

Publication number Publication date
JP5026363B2 (en) 2012-09-12
CN101488456A (en) 2009-07-22
TW200949932A (en) 2009-12-01
CN101488456B (en) 2010-12-15
TWI508161B (en) 2015-11-11

Similar Documents

Publication Publication Date Title
JP5026363B2 (en) Etching amount calculation method, storage medium, and etching amount calculation device
KR101073229B1 (en) Etching amount calculation method, storage medium, and etching amount calculation apparatus
US10453695B2 (en) Plasma processing apparatus and plasma processing method
JP4925507B2 (en) Film thickness control using spectral interferometry
KR101460375B1 (en) Substrate processing control method and storage medium
US11424110B2 (en) Plasma processing apparatus and operational method thereof
US20080078948A1 (en) Processing termination detection method and apparatus
TWI615901B (en) Plasma processing apparatus and plasma processing method
TW201724161A (en) Plasma processing apparatus and operating method of plasma processing apparatus
JP5080775B2 (en) Processing end point detection method and processing end point detection device
JP5199981B2 (en) Etching depth detection method, etching monitor apparatus and etching apparatus
US11404253B2 (en) Plasma processing apparatus and analysis method for analyzing plasma processing data
WO2018105462A1 (en) Signal processing method and program
KR101529827B1 (en) Method for detecting endpoint of plasma etching
JP2742446B2 (en) Etching method and apparatus
JP4909372B2 (en) Etching depth detection method, etching monitor apparatus and etching apparatus
CN1380541A (en) Method for measuring temp. and monitoring etching rate by using optical method
JP7201828B2 (en) Plasma processing apparatus and plasma processing method
JP2016009720A (en) Estimation method, and plasma processing apparatus
JP2000111313A (en) Device and method for detecting etching starting point, and plasma etching processor
KR101358538B1 (en) Apparatus for detecting an etching end point and method of detecting an etching end point

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110630

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

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120524

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120620

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

Free format text: PAYMENT UNTIL: 20150629

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 5026363

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