WO2005071144A1 - シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶媒体 - Google Patents

シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶媒体 Download PDF

Info

Publication number
WO2005071144A1
WO2005071144A1 PCT/JP2005/001114 JP2005001114W WO2005071144A1 WO 2005071144 A1 WO2005071144 A1 WO 2005071144A1 JP 2005001114 W JP2005001114 W JP 2005001114W WO 2005071144 A1 WO2005071144 A1 WO 2005071144A1
Authority
WO
WIPO (PCT)
Prior art keywords
silicon single
single crystal
oxygen
heat treatment
concentration
Prior art date
Application number
PCT/JP2005/001114
Other languages
English (en)
French (fr)
Inventor
Kozo Nakamura
Junsuke Tomioka
Tetsuro Akagi
Shiro Yoshino
Original Assignee
Komatsu Denshi Kinzoku Kabushiki Kaisha
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 Komatsu Denshi Kinzoku Kabushiki Kaisha filed Critical Komatsu Denshi Kinzoku Kabushiki Kaisha
Priority to JP2005517323A priority Critical patent/JP4901217B2/ja
Priority to US10/586,445 priority patent/US8246744B2/en
Publication of WO2005071144A1 publication Critical patent/WO2005071144A1/ja

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C30CRYSTAL GROWTH
    • C30BSINGLE-CRYSTAL GROWTH; UNIDIRECTIONAL SOLIDIFICATION OF EUTECTIC MATERIAL OR UNIDIRECTIONAL DEMIXING OF EUTECTOID MATERIAL; REFINING BY ZONE-MELTING OF MATERIAL; PRODUCTION OF A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; SINGLE CRYSTALS OR HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; AFTER-TREATMENT OF SINGLE CRYSTALS OR A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; APPARATUS THEREFOR
    • C30B15/00Single-crystal growth by pulling from a melt, e.g. Czochralski method
    • C30B15/20Controlling or regulating
    • C30B15/206Controlling or regulating the thermal history of growing the ingot
    • CCHEMISTRY; METALLURGY
    • C30CRYSTAL GROWTH
    • C30BSINGLE-CRYSTAL GROWTH; UNIDIRECTIONAL SOLIDIFICATION OF EUTECTIC MATERIAL OR UNIDIRECTIONAL DEMIXING OF EUTECTOID MATERIAL; REFINING BY ZONE-MELTING OF MATERIAL; PRODUCTION OF A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; SINGLE CRYSTALS OR HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; AFTER-TREATMENT OF SINGLE CRYSTALS OR A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; APPARATUS THEREFOR
    • C30B29/00Single crystals or homogeneous polycrystalline material with defined structure characterised by the material or by their shape
    • C30B29/02Elements
    • C30B29/06Silicon
    • CCHEMISTRY; METALLURGY
    • C30CRYSTAL GROWTH
    • C30BSINGLE-CRYSTAL GROWTH; UNIDIRECTIONAL SOLIDIFICATION OF EUTECTIC MATERIAL OR UNIDIRECTIONAL DEMIXING OF EUTECTOID MATERIAL; REFINING BY ZONE-MELTING OF MATERIAL; PRODUCTION OF A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; SINGLE CRYSTALS OR HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; AFTER-TREATMENT OF SINGLE CRYSTALS OR A HOMOGENEOUS POLYCRYSTALLINE MATERIAL WITH DEFINED STRUCTURE; APPARATUS THEREFOR
    • C30B33/00After-treatment of single crystals or homogeneous polycrystalline material with defined structure
    • C30B33/02Heat treatment
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01LSEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
    • H01L21/00Processes or apparatus adapted for the manufacture or treatment of semiconductor or solid state devices or of parts thereof
    • H01L21/02Manufacture or treatment of semiconductor devices or of parts thereof
    • H01L21/04Manufacture or treatment of semiconductor devices or of parts thereof the devices having potential barriers, e.g. a PN junction, depletion layer or carrier concentration layer
    • H01L21/18Manufacture or treatment of semiconductor devices or of parts thereof the devices having potential barriers, e.g. a PN junction, depletion layer or carrier concentration layer the devices having semiconductor bodies comprising elements of Group IV of the Periodic Table or AIIIBV compounds with or without impurities, e.g. doping materials
    • H01L21/30Treatment of semiconductor bodies using processes or apparatus not provided for in groups H01L21/20 - H01L21/26
    • H01L21/322Treatment of semiconductor bodies using processes or apparatus not provided for in groups H01L21/20 - H01L21/26 to modify their internal properties, e.g. to produce internal imperfections
    • H01L21/3221Treatment of semiconductor bodies using processes or apparatus not provided for in groups H01L21/20 - H01L21/26 to modify their internal properties, e.g. to produce internal imperfections of silicon bodies, e.g. for gettering
    • H01L21/3225Thermally inducing defects using oxygen present in the silicon body for intrinsic gettering

Definitions

  • Storage medium for storing a method for predicting oxygen precipitation behavior in silicon single crystal, a method for determining production parameters of silicon single crystal, and a program for predicting oxygen precipitation behavior in silicon single crystal
  • the present invention relates to a method for predicting the behavior of oxygen precipitates in a silicon single crystal for predicting the behavior of oxygen precipitates generated in a silicon single crystal in response to a heat treatment, and determining a manufacturing parameter of the silicon single crystal using the method.
  • the present invention relates to a method and a storage medium for storing a program for predicting oxygen precipitation behavior in a silicon single crystal.
  • a raw material that is, a solid silicon raw material is charged into a quartz crucible, and a raw material melt is generated by heating.
  • oxygen existing on the surface layer of the inner wall of the crucible and on the surface of the material is mixed into the raw material melt.
  • 10 X 10 17 The silicon is produced single crystal pulling method - m 3 oxygen (old ASTM Display) extent is taken N 20 X 10 17 atoms.
  • the incorporated oxygen is supersaturated in the heat treatment process performed in the device manufacturing process, and precipitates in accordance with the heat treatment.
  • the oxygen precipitate thus formed is a small defect, but is effective as a gettering site for impurities. Gettering using oxygen precipitates is particularly called intrinsic gettering (hereinafter referred to as "IG") and is widely adopted as a method for gettering harmful heavy metals.
  • IG intrinsic gettering
  • the IG ability for heavy metals is related to the density and size of oxygen precipitates, that is, the density of oxygen precipitates and the amount of oxygen precipitates.
  • Patent Literature 1 discloses Fe
  • Patent Literatures 2 and 3 disclose Ni and Cu.
  • the oxygen concentration and the silicon concentration of the silicon single crystal are adjusted so that an appropriate IG capability can be obtained by the density and the amount of oxygen precipitate generated in the silicon single crystal in accordance with the heat treatment performed on the wafer in the device manufacturing process. The process was chosen.
  • the density of oxygen precipitates is determined by a selective etching method.
  • the amount of precipitated oxygen can be determined by taking the difference in the amount of infrared absorption of dissolved oxygen before and after the heat treatment.
  • recent elements Manufacturing processes have shifted from conventional high-temperature processes to low-temperature processes, and it has become difficult to evaluate the density of oxygen precipitates and the amount of oxygen precipitates using conventional methods.
  • Patent Document 4 discloses a method of calculating the density of oxygen precipitates and the amount of oxygen precipitates by calculation. This technique is based on three parameters: initial oxygen concentration in silicon single crystal, dopant concentration or resistivity in silicon single crystal, and heat treatment conditions applied to silicon single crystal. This is a method for determining the density and amount of oxygen precipitates generated in a silicon single crystal. A similar method is disclosed in Patent Documents 2 and 3 mentioned above. According to these methods, it is possible to evaluate the density of oxygen precipitates and the amount of oxygen precipitates in a low-temperature process.
  • the precipitation of oxygen strongly depends on the concentration of the thermal donor generated in accordance with the thermal history between 400 ° C and 550 ° C that the silicon single crystal receives during crystal growth.
  • Thermal donors are oxygen clusters composed of several tens of tens of oxygen atoms, and are electrically measured as donors. Thermal donors occur at temperatures between 400 ° C and 550 ° C, and have a rapid rate of formation, especially at temperatures between 450 ° C and 500 ° C. The time (thermal history) of staying at a temperature between 450 ° C and 500 ° C is reflected in the thermal donor concentration.
  • Patent Document 1 JP 2003-257983 A
  • Patent Document 2 JP-A-2000-68280
  • Patent Document 3 Japanese Patent Application Laid-Open No. 2003-318181
  • Patent Document 4 JP-A-11-147789
  • Patent Document 5 JP-A-2-263792
  • Patent Document 6 JP-A-4-130732
  • Patent Document 7 Japanese Patent Application Laid-Open No. H4-298042
  • Patent Document 8 Japanese Patent Application Laid-Open No. 4-175300
  • Patent Document 9 JP-A-5-102167
  • the present invention has been made in view of such circumstances, and reflects the thermal donor concentration, which is the most important factor for oxygen precipitation, in the calculation method, thereby improving the effect of thermal donors on oxygen precipitation.
  • the goal is to improve the reliability of prediction of oxygen precipitates.
  • the first invention is In the method for predicting the behavior of oxygen precipitates in a silicon single crystal, which predicts the behavior of oxygen precipitates generated in a silicon single crystal in response to heat treatment,
  • the initial oxygen concentration in the silicon single crystal the concentration of the thermal donor generated according to the thermal history between 400 ° C and 550 ° C that the silicon single crystal receives during crystal growth, and the heat treatment applied to the silicon single crystal
  • the method is characterized in that an operation is performed using the conditions and parameters as parameters, and the density and the amount of oxygen precipitates generated in the silicon single crystal when the silicon single crystal is heat-treated are obtained.
  • the second invention is a first invention
  • the calculation is performed using the initial oxygen concentration in the silicon single crystal and the concentration of the thermal donor generated according to the thermal history between 400 ° C and 550 ° C received by the silicon single crystal during crystal growth as parameters. Obtain the nucleation rate of oxygen precipitates generated in the silicon single crystal during the heat treatment process when the silicon single crystal is heat-treated, and use the obtained nucleation rate to determine the density and amount of oxygen precipitates.
  • the third invention is
  • the method for predicting the oxygen precipitation behavior in a silicon single crystal according to the first or second invention is used so that the density and amount of oxygen precipitates in the silicon single crystal are set to desired values.
  • the thermal history of the silicon single crystal ingot during the crystal growth between 400 ° C and 550 ° C, and the heat treatment conditions applied to the silicon single crystal is determined.
  • the fourth invention is a
  • a storage medium storing a program for predicting the behavior of oxygen precipitates generated in a silicon single crystal in response to heat treatment by a computer
  • the fifth invention provides:
  • a storage medium storing a program for predicting the behavior of oxygen precipitates generated in a silicon single crystal in response to heat treatment by a computer
  • the nucleation rate of oxygen precipitates generated in the silicon single crystal during the heat treatment is determined, and the density and amount of oxygen precipitates are determined using the determined nucleation rate.
  • Kt Nucleation rate (cm- 3 s-, Ding: temperature, a (T): coefficient determined by the temperature, C: oxygen concentration (X 10 17 cm- 3), TD: Thermal donor concentration (X 10 15 cm “ 3 ) It is.
  • the amount of precipitated oxygen can be calculated by applying a general growth equation representing the growth and extinction of a nucleus generated at a critical nucleus radius at each temperature under diffusion control.
  • the thermal donor concentration can be determined from the change in resistivity before and after the donor erasure heat treatment, or from a thermal history between 400 ° C and 550 ° C and the oxygen concentration using a prediction formula.
  • the thermal donor generated according to the initial oxygen concentration in the silicon single crystal and the thermal history between 400 ° C and 550 ° C that the silicon single crystal receives during crystal growth By determining the concentration and the nucleation rate of oxygen precipitates generated in the silicon single crystal during the heat treatment process when the silicon single crystal is heat-treated, By specifying the heat treatment conditions, it can be said that the density of oxygen precipitates and the amount of oxygen precipitates under arbitrary heat treatment conditions can be predicted by calculation.
  • the above calculation method is programmed and stored in a storage medium.
  • the above operation is performed by the arithmetic unit of the computer. In this way, the calculation efficiency can be improved and an accurate calculation result can be obtained.
  • the invention's effect is programmed and stored in a storage medium.
  • the present invention it is possible to predict the density of oxygen precipitates and the amount of precipitated oxygen under arbitrary heat treatment conditions by calculating the oxygen concentration, the thermal donor concentration, and the heat treatment conditions as parameters.
  • the parameters of the silicon single crystal wafer manufacturing process are determined so that the heat treatment process performed in the device manufacturing process can obtain an oxygen precipitation state exhibiting appropriate IG capability.
  • the reliability of the calculation results is improved. This eliminates the need for labor-intensive work such as determining the conditions for obtaining an appropriate oxygen precipitation state for each heat treatment condition as in the related art. Therefore, it is easy to design an appropriate silicon single crystal wafer for any heat treatment process without reducing work efficiency. Is made.
  • the calculation consists of the following steps 1) -1).
  • the heat treatment process is divided into a plurality of time intervals, and the nucleation rate for each time interval is determined by the nucleation rate equation using the initial oxygen concentration, thermal donor concentration, and temperature, and the density of nuclei generated in each time interval.
  • the oxygen precipitate is a silicon oxygen compound of Si ⁇ and is spherical.
  • the thermal donor concentration of a silicon single crystal has a strong correlation with the time required to pass between 450 ° C and 500 ° C and the oxygen concentration during the cooling process during crystal growth. It is possible to ask. However, it is preferable to obtain it by the following calculation method.
  • TD thermal donor concentration
  • TD eq thermal equilibrium concentration of the thermal donor concentration
  • a coefficient ( 9. 2 X 10- 5 ° )
  • k Boltzmann constant
  • T absolute temperature
  • D diffusion coefficient of oxygen
  • C Oxygen concentration
  • t time
  • the time t is represented by the following equation (3).
  • TD TD (T 2 ) eq ⁇ l -expC- a DC (t 1 2 + t)) ⁇ ... (4)
  • FIG. 1 is a diagram showing a cooling curve at a straight body position in a crystal growth process of a silicon single crystal.
  • Fig. 1 shows the cooling curve for a crystal with a diameter of 150 mm, and the straight body starting position, the straight-line distance separated by a predetermined ⁇ giant separation (50, 150, 250, 350, 450, 550, 650, 750 mm) The cooling curve during crystal growth at the barrel position is shown.
  • the temperature of the silicon single crystal constantly decreases during crystal growth. Therefore, it is necessary to repeat the above-mentioned calculation. That is, the equivalent amount of thermal donor generation (concentration) TD at the temperature T up to the temperature T where the temperature has dropped is determined, and the equivalent time t is added with the time increase of the time step ⁇ t to add the amount of thermal donor generation. (Concentration) The calculation of TD is repeated.
  • the generation of the thermal donor is the most active force between 450 ° C and 500 ° C.
  • the temperature range to be calculated is 400 ° C power, which is wider than 550 ° C.
  • FIG. 2 is a diagram showing the thermal donor concentration obtained by calculation and the thermal donor concentration obtained by actual measurement.
  • the calculated value shown in FIG. 2 is the thermal donor concentration at each straight body position when the above-described calculation method is applied to the temperature history in FIG.
  • the heat treatment process of 2) is divided into a plurality of time intervals at predetermined intervals, for example, 5 second intervals, and the nucleation rate in each time interval is calculated by the nucleation rate equation using the initial oxygen concentration, the thermal donor concentration, and the temperature.
  • the process of obtaining the density of nuclei generated in each time interval by the calculation is described below.
  • the present inventors have found that the nucleation rate at each temperature is shown as a function of oxygen concentration and thermal donor concentration. The relationship is shown as the following equation (5). In the following equation (5),
  • T, C, TD nucleation rate (cm— 3 s—, T: temperature, a (T): coefficient determined by temperature, C: oxygen concentration (X 10 17 cm 3 ) TD: thermal donor concentration ( X 10 15 cm 1-3 )
  • FIG. 3 is a diagram showing the correspondence between the coefficients used in the equation for determining the nucleation rate of oxygen precipitates and the temperature.
  • the present inventors have found that the coefficient a (T) determined by the temperature has the value shown in FIG.
  • the thermal donor concentration TD used in the above equation (5) is an as grown concentration generated during the crystal growth process.
  • As the electrical properties there is a substance that disappears during heating in heat treatment.
  • the thermal donor generated during the heat treatment should also be added.
  • a feature of the present invention is that it has been found that the nucleation rate of oxygen precipitates is expressed by the above equation (5), and for the first time the effect of the thermal donor is applied to the calculation prediction.
  • N (t ') I (T, C, T D) A t... (6)
  • the radius of the generated nucleus is assumed to be a critical nucleus radius Rcri represented by the following equation (7).
  • R W Time radius of the nucleus generated at the time t before the force A t elapses
  • Ci Equilibrium oxygen concentration at the interface of the spherical particles of radius
  • R (, t) can be obtained by the calculation using the above equation (8). Since the change in the oxygen concentration in the silicon single crystal is represented by the following equation (9), the amount of precipitated oxygen can be calculated. Since the oxygen concentration decreases by the amount of precipitated oxygen, the residual oxygen concentration can be obtained as the difference between the initial oxygen concentration and the amount of oxygen precipitation obtained by using the following equation (9). The arithmetic processing of obtaining the residual oxygen concentration for each time interval, obtaining the nucleus generation rate of the oxygen precipitate, and obtaining the growth of the oxygen precipitate is repeated.
  • the above calculation method is programmed and stored in a storage medium.
  • the program is executed by inputting the parameter, the above calculation is performed by the calculation device of the computer. In this way, the calculation efficiency can be improved and an accurate calculation result can be obtained.
  • Example 1 As Example 1, verification of the prediction of oxygen precipitation behavior according to the present invention is performed for a typical heat treatment process.
  • the heat treatment performed is to put the wafer into a furnace at the set temperature and hold it.
  • Silicon single crystals used was an oxygen concentration of 11 X 10 17 - 17 X 10 17 atoms N m 3 (formerly
  • the process that occurs at 800 ° C or lower is called nucleation heat treatment, and the process where nucleus growth mainly occurs at 800 ° C or higher is called growth heat treatment.
  • a growth heat treatment was applied to measure the density of oxygen precipitates and the amount of oxygen precipitates.
  • the measured value and the calculated value of the amount of oxygen precipitate of the oxygen precipitate are compared.
  • the actually measured value was determined by taking the difference in the amount of infrared absorption of dissolved oxygen before and after the heat treatment. Also
  • FIG. 4 is a diagram showing measured values and calculated values of the amount of oxygen deposition when the growth heat treatment 1 is added after the nucleation heat treatment.
  • FIG. 5 shows measured and calculated values of the amount of precipitated oxygen when the growth heat treatment 2 was added after the nucleation heat treatment.
  • Fig. 6 shows the actual amount of oxygen precipitates when the growth heat treatment 3 was added after the nucleation heat treatment. It is a figure which shows a measured value and a calculated value. In each of the figures, the abscissa is the measured value and the ordinate is the calculated value. As the distribution closer to the straight line A is obtained in each figure, the difference between the measured value and the calculated value is smaller. In FIGS. 4, 5, and 6, the calculated value of the amount of precipitated oxygen and the measured value are almost the same, and a distribution close to the straight line A is obtained. From this result, it can be seen that the accuracy of the calculation according to the present invention is high.
  • the measured value was obtained by a selective etching method. Further, the calculated value was obtained by the calculation of the present invention. The results are shown in FIGS. 7, 8, and 9.
  • FIG. 7 is a diagram showing measured and calculated values of the density of oxygen precipitates in the case where the growth heat treatment 1 is added after the nucleation heat treatment.
  • FIG. 8 is a diagram showing measured values and calculated values of the density of oxygen precipitates when the growth heat treatment 2 is added after the nucleation heat treatment.
  • FIG. 9 is a diagram showing measured and calculated values of the density of oxygen precipitates when the growth heat treatment 3 is added after the nucleation heat treatment.
  • the horizontal axis represents measured values and the vertical axis represents computed values. In each figure, the closer the distribution closer to the straight line A is obtained, the smaller the difference between the measured values and computed values.
  • FIGS. 7, 8, and 9 the calculated values of the density of oxygen precipitates and the measured values are almost the same, and a distribution close to the straight line A is obtained. From this result, it can be seen that the accuracy of the calculation according to the present invention is high.
  • Patent Document 6 described above by the present inventors discloses a method for predicting the amount of precipitated oxygen.
  • the method disclosed in Patent Document 6 is a method in which the growth heat treatment under the three conditions described in the present embodiment is limited, and the amount of precipitated oxygen is determined by multivariate analysis. Therefore, this method was not applicable to a wide range of heat treatment conditions, and was not versatile. In addition, no information was obtained on the density and average size of the oxygen precipitates.
  • the density of oxygen precipitates can be predicted with excellent accuracy together with the amount of oxygen precipitates. This means that the average number of oxygen atoms constituting the oxygen precipitate, that is, the average size of the oxygen precipitate can be predicted.
  • Example 2 the prediction of the oxygen precipitation behavior according to the present invention was verified for a heat treatment process different from that of Example 1, and the versatility of the present invention was confirmed.
  • the heat treatment was performed in a low-temperature furnace. ⁇ It is to throw the e-ha and gradually heat it to the set temperature.
  • the wafer was charged into a furnace at 500 ° C., and after the charging, the temperature was raised at a rate of 1 ° C.Z. Then, four temperatures of 650, 700, 750, and 800 ° C were set as the reached temperatures, and after the furnace temperature reached the set temperature, each temperature was maintained for 4 hours. After performing such a nucleation heat treatment on Aha Co., a growth heat treatment was performed at 780 ° C for 3 hours + 1000 ° C for 16 hours, and the density of oxygen precipitates and the amount of oxygen precipitates were measured.
  • the measured value and the calculated value of the amount of oxygen precipitate of the oxygen precipitate are compared.
  • FIG. 10 shows the results.
  • FIG. 10 shows measured and calculated values of the amount of precipitated oxygen when the growth heat treatment was applied after the nucleation heat treatment.
  • Fig. 10 shows the measured values on the horizontal axis and the calculated values on the vertical axis, and the closer the distribution closer to the straight line A is obtained in Fig. 10, the more the difference between the measured value and the calculated value is It can be said that there are few.
  • the calculated value of the amount of precipitated oxygen and the measured value are almost the same, and a distribution close to the straight line A is obtained. From this result, it can be seen that the accuracy of the calculation according to the present invention is high.
  • FIG. 11 is a diagram showing measured and calculated values of the density of oxygen precipitates when a growth heat treatment is applied after a nucleation heat treatment.
  • Fig. 11 shows the measured values on the horizontal axis and the calculated values on the vertical axis.
  • the calculated value of the density of oxygen precipitates and the measured value are almost the same, and a distribution close to the straight line A is obtained. From this result, it can be seen that the accuracy of the operation according to the present invention is high.
  • the density of oxygen precipitates can be predicted with excellent accuracy together with the amount of oxygen precipitates. This means that the average number of oxygen atoms constituting the oxygen precipitate, that is, the average size of the oxygen precipitate can be predicted.
  • FIG. 1 is a diagram showing a cooling curve at a straight body position in a crystal growth process of a silicon single crystal.
  • FIG. 2 is a diagram showing a thermal donor concentration obtained by calculation and a thermal donor concentration obtained by actual measurement.
  • FIG. 3 is a diagram showing a correspondence relationship between a coefficient used in an equation for calculating a nucleation rate of oxygen precipitates and temperature.
  • FIG. 4 is a diagram showing measured and calculated values of the amount of precipitated oxygen when a nucleation heat treatment is performed and then a growth heat treatment is performed at 800 ° C for 4 hours + 1000 ° C for 16 hours.
  • Fig. 4 is a diagram showing measured and calculated values of the amount of precipitated oxygen when a growth heat treatment is performed at 900 ° C for 4 hours + 1000 ° C for 16 hours after the nucleation heat treatment.
  • Fig. 6 is a diagram showing measured and calculated values of the amount of precipitated oxygen when a nucleation heat treatment is performed and then a growth heat treatment is performed at 1000 ° C for 16 hours.
  • FIG. 7 is a graph showing the density and calculated values of oxygen precipitates when a growth heat treatment was applied at 800 ° C for 4 hours + 1000 ° C for 16 hours after a nucleation heat treatment.
  • Fig. 8 shows the measured and calculated values of the density of oxygen precipitates when the growth heat treatment was applied at 900 ° C for 4 hours + 1000 ° C for 16 hours after the nucleation heat treatment. is there.
  • FIG. 9 is a diagram showing measured and calculated values of the density of oxygen precipitates when a growth heat treatment at 1000 ° C. for 16 hours is performed after a nucleation heat treatment.
  • FIG. 10 is a diagram showing measured and calculated values of oxygen precipitation when a nucleation heat treatment is performed and then a growth heat treatment at 780 ° C for 3 hours + 1000 ° C for 16 hours is applied. .
  • Fig. 11 shows the measured and calculated values of the density of oxygen precipitates when the nucleation heat treatment was performed and the growth heat treatment was performed at 780 ° C for 3 hours + 1000 ° C for 16 hours. It is.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Metallurgy (AREA)
  • Materials Engineering (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Power Engineering (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computer Hardware Design (AREA)
  • Manufacturing & Machinery (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Thermal Sciences (AREA)
  • Crystals, And After-Treatments Of Crystals (AREA)

Abstract

 シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける400°Cから550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、を特定することによって、シリコン単結晶を熱処理した場合にその熱処理過程でシリコン単結晶中に発生する酸素析出物の核発生速度を求める。さらにシリコン単結晶に加える熱処理条件を特定することによって、任意の熱処理条件における酸素析出物の密度および酸素析出量を演算によって予測する。

Description

明 細 書
シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメ ータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する 記憶媒体
技術分野
[0001] 本発明は、熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動を予測 するシリコン単結晶中の酸素析出挙動予測方法、この方法を利用したシリコン単結 晶の製造パラメータ決定方法およびシリコン単結晶中の酸素析出挙動予測用プログ ラムを記憶する記憶媒体に関する。
背景技術
[0002] 引き上げ法においては、まず石英製のルツボに素材すなわち固体状のシリコン原 料が投入され、加熱によって原料融液が生成される。この際、ルツボ内壁の表層や 素材表面等に存在していた酸素が原料融液中に混入する。このため、引き上げ法で 製造されるシリコン単結晶には 10 X 1017— 20 X 1017atomsん m3 (旧 ASTM表示)程 度の酸素が取り込まれる。取り込まれた酸素は素子製造工程で行われる熱処理過程 において過飽和であり、熱処理に応じて析出する。こうして形成される酸素析出物は 微小欠陥ではあるが、不純物のゲッタリングサイトとして有効である。酸素析出物を利 用したゲッタリングは特にイントリンシックゲッタリング(以下、 "IG"という)と呼ばれて おり、有害な重金属のゲッタリング方法として広く採用されている。
[0003] 重金属に対する IG能は酸素析出物の密度およびそのサイズ、つまり酸素析出物の 密度および酸素析出量と関係がある。この関係に関して、例えば下記特許文献 1に は Feについて開示されており、下記特許文献 2、 3には Ni、 Cuについて開示されて いる。従来は、素子製造工程でゥヱーハに施される熱処理に応じてシリコン単結晶中 に発生する酸素析出物の密度および酸素析出量によって適切な IG能が得られるよう に、シリコン単結晶の酸素濃度およびプロセスが選ばれていた。
[0004] 酸素析出物の密度は選択エッチング法で求められる。また、酸素析出量は熱処理 前後での固溶酸素の赤外線吸収量の差をとることで求められる。しかし、近年の素子 製造工程は従来の高温プロセスから低温プロセスへと移行しており、従来の方法で は酸素析出物の密度および酸素析出量を評価することが困難となってきている。
[0005] 下記特許文献 4には演算によって酸素析出物の密度および酸素析出量を割り出す 方法が開示されている。この技術は、シリコン単結晶中の初期酸素濃度と、シリコン 単結晶中のドーパント濃度または抵抗率と、シリコン単結晶に加える熱処理条件とい う三つのパラメータを用いて、シリコン単結晶を熱処理した場合にシリコン単結晶中に 発生する酸素析出物の密度および析出量を求める方法である。同じ様な方法は前 述した特許文献 2、 3でも開示されている。これらの方法によれば低温プロセスにおけ る酸素析出物の密度および酸素析出量を評価することが可能ではある。
[0006] 特許文献 2— 4で開示された酸素析出物の密度および析出量を予測する方法は、 次のような仮定 a)、 b)に基づくものである。
a)酸素析出物の核発生は、過飽和酸素の自由エネルギーを駆動力とした均一核 発生プロセスである。
b)均一核発生プロセスによって発生した核は、酸素の拡散律速の過程によって成 長する。
[0007] ところで、酸素の析出はシリコン単結晶が結晶成長の際に受ける 400°Cから 550°C の間の熱履歴に応じて発生するサーマルドナーの濃度に強く依存する、とレ、うこと力 S 従来からよく知られている。サーマルドナーは数個力 数十個の酸素原子からなる酸 素クラスターであり、電気的にドナーとして計測される。サーマルドナーは 400°Cから 550°Cの間の温度において発生し、特に 450°Cから 500°Cの間の温度においてそ の形成速度が速い。この 450°Cから 500°Cの間の温度に留まる時間(熱履歴)は、サ 一マルドナー濃度に反映される。
[0008] 400°Cから 550°Cの間の熱履歴、言い換えればサーマルドナー濃度がその後の熱 処理過程における酸素析出を決定することは、例えば下記特許文献 5 9に開示さ れている。これらの技術は、好ましい酸素析出を得るために 450°Cから 500°Cの間の 熱履歴つまりサーマルドナー濃度を調整すること、に関するものである。
特許文献 1 :特開 2003 - 257983号公報
特許文献 2:特開 2000 - 68280号公報 特許文献 3:特開 2003 - 318181号公報
特許文献 4:特開平 11 - 147789号公報
特許文献 5:特開平 2 - 263792号公報
特許文献 6 :特開平 4 - 130732号公報
特許文献 7:特開平 4 - 298042号公報
特許文献 8:特開平 4 - 175300号公報
特許文献 9:特開平 5 - 102167号公報
発明の開示
発明が解決しょうとする課題
[0009] ここで問題なのは、上記特許文献 4などのように演算によって酸素析出物の密度お よび酸素析出量を求める方法ではサーマルドナーの酸素析出への効果を現すこと ができないことである。つまり、前述したように酸素析出がサーマルドナーに依存する にも関わらず、特許文献 4などではサーマルドナーについて全く考慮されていないこ とが問題である。
[0010] 現在はサーマルドナーと酸素析出との関係について理論的に説明するモデルがな く、当てはめる理論式が存在しないため、その効果を考慮した演算がなされていない 。このように酸素析出物の密度および酸素析出量を求める従来の演算方法には最も 重要な要素が抜けていることから、その演算結果は信頼性が低レ、ものと考えられる。 信頼性の低さを補うためには、熱処理条件毎に適切な酸素析出状態になる条件を 決定する、といった多大な労力を要する作業が必要となり、さらに熱処理プロセス毎 にフィッティングパラメータの調整が必要となる。すると、作業効率の低下が避けられ ない。したがって、従来の酸素析出物の演算方法は汎用性に乏しかった。
[0011] 本発明はこうした実状に鑑みてなされたものであり、酸素析出に対して最も重要な 要素とされているサーマルドナー濃度を演算方法に反映することによってサーマルド ナ一の酸素析出への効果を表現し、酸素析出物の予測に関する信頼性を向上させ ることを角军決課題とするものである。
課題を解決するための手段
[0012] 第 1発明は、 熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動を予測するシリコ ン単結晶中の酸素析出挙動予測方法において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、シリコン単結 晶に加える熱処理条件と、をパラメータとして演算を行い、シリコン単結晶を熱処理し た場合にシリコン単結晶中に発生する酸素析出物の密度および析出量を求めること を特徴とする。
[0013] 第 2発明は、
熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動を予測するシリコ ン単結晶中の酸素析出挙動予測方法において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、をパラメータ として演算を行い、シリコン単結晶を熱処理した場合にその熱処理過程でシリコン単 結晶中に発生する酸素析出物の核発生速度を求め、求めた核発生速度を用いて酸 素析出物の密度および析出量を求めること
を特徴とする。
[0014] 第 3発明は、
シリコン単結晶の製造パラメータ決定方法において、
シリコン単結晶中の酸素析出物の密度および析出量を所望の値にするように、第 1 発明又は第 2発明のシリコン単結晶中の酸素析出挙動予測方法を利用して、シリコ ン単結晶中の初期酸素濃度と、シリコン単結晶インゴットが結晶成長の際に受ける 4 00°Cから 550°Cの間の熱履歴と、シリコン単結晶に加える熱処理条件と、を決定する こと
を特徴とする。
[0015] 第 4発明は、
熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動をコンピュータで 予測するためのプログラムを記憶する記憶媒体において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、シリコン単結 晶に加える熱処理条件と、をパラメータとして演算を行う処理と、
シリコン単結晶を熱処理した場合にシリコン単結晶中に発生する酸素析出物の密 度および析出量を求める処理と、
をプログラムとして記憶すること
を特徴とする。
[0016] 第 5発明は、
熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動をコンピュータで 予測するためのプログラムを記憶する記憶媒体において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、をパラメータ として演算を行う処理と、
シリコン単結晶を熱処理した場合にその熱処理過程でシリコン単結晶中に発生す る酸素析出物の核発生速度を求め、求めた核発生速度を用いて酸素析出物の密度 および析出量を求める処理と、
をプログラムとして記憶すること
を特徴とする。
[0017] 前述したようにサーマルドナーと酸素析出との関係について理論的に説明するモ デルがないため、本発明者らは系統的な実験を行い、その結果を詳細に解析するこ とによって酸素析出挙動を明かにした。本発明者らは上記特許文献 6、 7において種 々の温度の熱処理における酸素析出量へのサーマルドナー濃度の影響について詳 細に示した力 S、本発明を得るためにこれらの実験サンプルをさらに詳細に調査した。 また、上記特許文献 6、 7では酸素析出量のみの評価であつたが、酸素析出物の密 度の評価を加えて行った。その結果以下のような結論 A)、 B)が得られた。
A)各温度での核発生速度は下記(1)式のような酸素濃度とサーマルドナー濃度の 関数として示される。下記(1)式において、
Kt):核発生速度(cm— 3s— 、丁:温度、 a (T):温度によって定まる係数、 C :酸素濃度 ( X 1017cm— 3)、 TD:サーマルドナー濃度( X 1015cm"3) である。
[数 1]
Figure imgf000008_0001
[0018] B)酸素析出量は、各温度において臨界核半径で発生した核の拡散律速での成長 および消滅を表す一般的な成長式を適用することによって演算することができる。サ 一マルドナー濃度はドナー消去熱処理の前後の抵抗率の変化から求められ、また 4 00°Cから 550°Cの間の熱履歴と酸素濃度から予測式を用いても求められる。
[0019] 上記 A)、 B)によって、シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶 成長の際に受ける 400°Cから 550°Cの間の熱履歴に応じて発生するサーマルドナー の濃度と、を特定することによって、シリコン単結晶を熱処理した場合にその熱処理 過程でシリコン単結晶中に発生する酸素析出物の核発生速度を求めることができ、さ らにシリコン単結晶に加える熱処理条件を特定することによって、任意の熱処理条件 における酸素析出物の密度および酸素析出量を演算によって予測することができる 、といえる。
[0020] 上記演算方法はプログラム化され記憶媒体に記憶される。パラメータが入力されプ ログラムが実行されるとコンピュータの演算装置で上記演算が行われる。このようにす れば、演算効率を向上させることができ且つ正確な演算結果を得ることができる。 発明の効果
[0021] 本発明によれば、酸素濃度とサーマルドナー濃度と熱処理条件をパラメータとして 演算することによって、任意の熱処理条件における酸素析出物の密度および酸素析 出量を予測することができる。さらに、素子製造工程で行われる熱処理過程によって 適切な IG能を示す酸素析出状態が得られるように、シリコン単結晶ゥエーハの製造 工程のパラメータが決定される。酸素析出に対して最も重要な要素、すなわちサーマ ノレドナー濃度という要素を用いて演算することによって演算結果の信頼性は向上す る。このため、従来のように熱処理条件毎に適切な酸素析出状態になる条件を決定 する、といった労力を要する作業が必要なくなる。したがって、作業効率を低下させる ことなく任意の熱処理プロセスに対する適切なシリコン単結晶ゥエーハの設計が容易 になされる。
発明を実施するための最良の形態
[0022] 以下、本発明の実施の形態について図面を参照して説明する。
[0023] [1.演算について]
まず、酸素析出物の密度および酸素析出量を求める演算について説明する。演算 は次の 1)一 3)の過程からなる。
1)シリコン単結晶のインゴットが結晶成長中の冷却過程の際に受ける 400°Cから 5 50°Cの間の熱履歴と、初期酸素濃度と、力 シリコン単結晶中に発生するサーマル ドナーの濃度を求める過程
2)熱処理プロセスを複数の時間区間に区切り、各時間区間の核発生速度を初期 酸素濃度、サーマルドナー濃度、温度、を用いた核発生速度式によって求め、各時 間区間で発生する核の密度を求める過程
3)各時間区間で発生する核の成長および消滅を求める過程
なお酸素析出物は Si〇なるシリコン酸素化合物であり、球体であると想定する。
2
[0024] 〔過程 1)について〕
まず、 1)のサーマルドナーの濃度を求める過程について説明する。シリコン単結晶 のサーマルドナー濃度は、結晶成長中の冷却過程において 450°Cと 500°Cとの間を 通過するのに要した時間および酸素濃度に対して強い相関があり、経験式を用いて 求めることが可能である。しかし、以下の計算方法によって求めることが好ましい。
[0025] 基本式として [K.Wada, N.Inoue, Semiconductor Silicon, ed. H.R. Huff et
al.,(Electrochem. So , Pennington N.J" 1986)p778]によって提案された式を用いた 。この式は下記(2)式のように示される。下記(2)式は、温度一定の条件での熱処理 時間とサーマルドナー濃度との関係を示し、サーマルドナー濃度 TDは酸素濃度 Cと 温度(絶対温度) Tとの関数で表されることを示している。下記(2)式において、
TD :サーマルドナー濃度、 TDeq:サーマルドナー濃度の熱平衡濃度、 & :係数(= 9. 2 X 10— 5°)、 k:ボルツマン定数、 T:絶対温度、 D:酸素の拡散係数、 C:酸素濃度、 t :時間
である。 T D = T De q {l -exp (- a D C t)}
Ί. 535 e V、
但し T De ¾= a C3 exp
, k T
[0026] 式の形から、上記(2)式は時間 t =∞で TDeqに収束する関数である、ということが分 かる。なお、上記(2)式は温度一定という条件で成立するため、この式を結晶成長中 の冷却過程という温度が変化する条件で適用するために、次のような工夫をした。
[0027] 温度変化と共にサーマルドナーが発生する場合は、上記(2)式の時間 tを温度毎 の実効時間に換算しなくてはならない。ここで、温度 Tで時間 t (sec)の熱処理を施
1 1
し、サーマルドナー濃度が TDとなった後に、温度 Tで時間 t (sec)の熱処理を施す
1 2 2
場合を考える。このような場合は、温度 τで発生したサーマルドナーの生成量 (濃度
) TDを温度 Tで発生させる場合を考え、その時間 t を演算する。上記(2)式から、
1 2 12
時間 t は下記(3)式のように示される。
12
[数 3]
Figure imgf000010_0001
[0028] すると、温度 Tで時間 t (sec)の熱処理を追加する場合は、サーマルドナー生成量 (
2
濃度) TDは下記 (4)式のように示される、と考えられる。
[数 4]
T D = T D(T 2)e q {l -expC- a D C(t 1 2 + t))} … (4)
図 1はシリコン単結晶の結晶成長過程における直胴位置での冷却曲線を示す図で ある。図 1では直径 150mmの結晶についての冷却曲線が示されており、直胴開始位 置力、ら所定 β巨離(50、 150、 250、 350、 450、 550、 650、 750mm)だけ離隔した直 胴位置における結晶成長中の冷却曲線が示されている。
図 1に示されるように、結晶成長中にシリコン単結晶の温度は常に降下し続ける。こ のため、前述した演算を繰り返し行う必要がある。すなわち、降下した温度 Tまでの サーマルドナー生成量 (濃度) TDを温度 Tで生成する場合の等価時間 tを求め、 等価時間 tへそのタイムステップ Δ t分の時間増を加えてサーマルドナー生成量 (濃 度) TDを求める、といった演算を繰り返すことになる。サーマルドナーの発生は 450 °Cから 500°Cの間が最も活発である力 演算する温度範囲としてはそれよりも広い 40 0°C力、ら 550°Cとすること力 S好ましレ、。
[0030] 図 2は演算によって求めたサーマルドナー濃度と実測によって求めたサーマルドナ 一濃度とを示す図である。図 2で示される演算値は、前述した演算方法を図 1の温度 履歴に適用した場合の各直胴位置におけるサーマルドナー濃度である。
図 2によれば、サーマルドナー濃度の演算値と実測値との誤差が微小であることが 分かる。このこと力 、前述した演算方法はサーマルドナー濃度を精度よく予測できる 、といえる。
[0031] 〔過程 2)について〕
次に、 2)の熱処理プロセスを所定間隔、例えば 5秒間隔の複数の時間区間に区切 り、各時間区間の核発生速度を初期酸素濃度、サーマルドナー濃度、温度、を用い た核発生速度式によって求め、各時間区間で発生する核の密度を求める過程につ いて説明する。本発明者らは、各温度での核発生速度が酸素濃度とサーマルドナー 濃度の関数として示されることを見出した。その関係は下記(5)式のように示される。 下記(5)式において、
I (T, C, TD):核発生速度(cm— 3s— 、 T:温度、 a (T):温度によって定まる係数、 C : 酸素濃度( X 1017cm 3) TD:サーマルドナー濃度( X 1015cm一3)
である。
[数 5]
1
I , C, T D)= a(T)C9T D3 (5)
[0032] 上記(5)式にっレ、て補足説明をする。 図 3は酸素析出物の核発生速度を求める式に用いられる係数と温度との対応関係 を示す図である。本発明者らは、温度によって定まる係数 a (T)は図 3に示す値となる ことを見出した。また、上記(5)式で用いるサーマルドナー濃度 TDは結晶成長過程 において生じる as grown濃度である。それらは電気的性質としては熱処理における昇 温中に消滅する力 実体作用は残るため、 as grownでの濃度を用いる。さらに、熱処 理においても 400°Cから 550°Cのプロセスを加える場合には、熱処理中に発生する サーマルドナーもこれに加える。本発明の特徴は、酸素析出物の核発生速度が上記 (5)式で示されることを見出し、初めてサーマルドナーの効果を演算予測に適用した ことである。
[0033] 時間 力 A t経過するまでの間に核発生する核の密度 N ( )は下記(6)式の ように示される。
N(t ')= I(T, C,T D)A t … (6)
[0034] そして、発生した核の半径は下記(7)式で示される臨界核半径 Rcriであるとする。
下記(7)式において、
1^^ :臨界核半径、¥: 31〇分子体積、 σ: Si-SiO間の界面エネルギー、 C:酸素
2 2
濃度、 Ceq:酸素の熱平衡濃度
である。
[数 7]
Figure imgf000012_0001
[0035] 上記(6)式、(7)式で示される密度およびサイズにて発生した時間区間毎の核のデ ータをメモリーに記録して、それぞれの核について成長および消滅を計算していく。
[0036] 〔過程 3)について〕
次に、 3)の各時間区間に発生する核の成長および消滅を求める過程について説 明する。時間 力 A t経過するまでに核発生する核の時間 tにおける成長速度は 下記(8)式のように示される。下記(8)式にぉレ、て、
R W , :時間 力 A t経過するまでの間に核発生する核の時間 tにおける半径 、 Ci :半径 Rの球状粒子の界面における平衡酸素濃度
である。
[数 8]
Figure imgf000013_0001
[0037] 上記(8)式を用いた演算によって、 R ( , t)を求めることができる。シリコン単結晶 中の酸素濃度の変化は下記(9)式のように示されるため、演算によって酸素析出量 を求めることができる。酸素濃度は析出した酸素の量だけ減少するため、残存酸素 濃度は、初期酸素濃度から下記(9)式を用いて求めた酸素析出量の差として求める こと力 Sできる。時間区間毎に残存酸素濃度を求め、酸素析出物の核発生速度を求め 、酸素析出物の成長を求めるという演算処理を繰り返す。
dC
4 πΌ t)(C— Ci)d t ' ■■■ (9) 8 t J N(t ,] (t ',
[0038] 以上の演算を熱処理プロセスに適用することによって、酸素析出物の密度および サイズ分布そして酸素析出量を求めることができる。具体的には、上記演算方法はプ ログラム化され記憶媒体に記憶される。ノ メータが入力されプログラムが実行される とコンピュータの演算装置で上記演算が行われる。このようにすれば、演算効率を向 上させることができ且つ正確な演算結果を得ることができる。
[0039] 逆に、酸素析出物の密度や酸素析出量を所望値にするためには、演算で用いられ る熱処理プロセスの各パラメータを調整すればよい。 [0040] [2.実施例]
次に、本発明の具体的な実施例を示す。なお、本発明はこれらに限定されるもので はない。
実施例 1
[0041] 実施例 1として、代表的な熱処理プロセスに対して本発明による酸素析出挙動予測 の検証を行う。行った熱処理は、設定温度の炉にゥエーハを投入して保持するという ものである。
[0042] 使用したシリコン単結晶は、酸素濃度が 11 X 1017— 17 X 1017atomsん m3 (旧
ASTM表示)であり、 as groun状態でのサーマルドナー濃度が 1 X 10"— 2 X 1015cm— 3である直径 150mmの P型結晶である。以下では、熱処理プロセスのうち、核発生が 活発に起こる 800°C以下のプロセスを核発生熱処理といい、核の成長が主に起こる 8 00°C以上のプロセスを成長熱処理という。本発明者らは以下に示す種々の条件で 核発生熱処理を施した後に成長熱処理を加えて、酸素析出物の密度と酸素析出量 を測定した。
[0043] 最初に、下記各温度と各時間とを組合せた条件において核発生熱処理を施した。 · 温度: 575、 600、 625、 650、 700、 750°C
•時間: 0、 0. 5、 2、 4時間
次に、下記三条件の成長熱処理を加えた。
•成長熱処理 1 : 800°Cで 4時間 + 1000°Cで 16時間
•成長熱処理 2: 900°Cで 4時間 + 1000°Cで 16時間
•成長熱処理 3 : 1000°Cで 16時間
まず、酸素析出物の酸素析出量について実測値と演算値とを比較する。
[0044] 実測値を熱処理前後での固溶酸素の赤外線吸収量の差をとることで求めた。また
、演算値を本発明の演算によって求めた。その結果を図 4、図 5、図 6に示す。
[0045] 図 4は上記核発生熱処理を施した後に上記成長熱処理 1を加えた場合の酸素析 出量の実測値と演算値を示す図である。図 5は上記核発生熱処理を施した後に上記 成長熱処理 2を加えた場合の酸素析出量の実測値と演算値を示す図である。図 6は 上記核発生熱処理を施した後に上記成長熱処理 3を加えた場合の酸素析出量の実 測値と演算値を示す図である。何れの図も横軸が実測値、縦軸が演算値とされてお り、各図において直線 Aに近い分布が得られるほど、実測値と演算値の差が少ないと レ、える。図 4、図 5、図 6においては酸素析出量の演算値と実測値とがほぼ一致し、直 線 Aに近い分布が得られている。この結果から本発明による演算の精度が高いことが 分かる。
[0046] 次に酸素析出物の密度について実測値と演算値とを比較する。
実測値を選択エッチング法によって求めた。また、演算値を本発明の演算によって 求めた。その結果を図 7、図 8、図 9に示す。
図 7は上記核発生熱処理を施した後に上記成長熱処理 1を加えた場合の酸素析 出物の密度の実測値と演算値を示す図である。図 8は上記核発生熱処理を施した後 に上記成長熱処理 2を加えた場合の酸素析出物の密度の実測値と演算値を示す図 である。図 9は上記核発生熱処理を施した後に上記成長熱処理 3を加えた場合の酸 素析出物の密度の実測値と演算値を示す図である。何れの図も横軸が実測値、縦 軸が演算値とされており、各図において直線 Aに近い分布が得られるほど、実測値と 演算値の差が少ないといえる。図 7、図 8、図 9においては酸素析出物の密度の演算 値と実測値とがほぼ一致し、直線 Aに近い分布が得られている。この結果から本発明 による演算の精度が高いことが分かる。
[0047] 本発明者らによる上記特許文献 6では、酸素析出量の予測方法について開示され ている。しかし、上記特許文献 6で開示されているのは、本実施例で示した上記三条 件の成長熱処理に限定し、酸素析出量を多変量解析によって求める方法である。よ つて、この方法は広範囲の熱処理条件に適用できず、汎用性がなかった。さらに、酸 素析出物の密度および平均サイズについての情報も得られなかった。
[0048] 本発明によれば、酸素析出物の密度を酸素析出量ともに優れた精度で予測するこ とができる。それは酸素析出物を構成する平均酸素原子数、すなわち酸素析出物の 平均サイズをも予測できることを意味する。
実施例 2
[0049] 実施例 2として、実施例 1と異なる熱処理プロセスに対して本発明による酸素析出 挙動予測の検証を行い、本発明の汎用性を確認した。行った熱処理は、低温の炉に ゥエーハを投入して設定温度まで徐熱するというものである。
[0050] 使用したシリコン結晶は、酸素濃度が 11. 2 X 1017— 13. 5 X 1017aton^cm3 (旧
ASTM表示)である直径 200mmの P型結晶である。本実施例で行った熱処理では、ゥ エーハを 500°Cの炉に投入し、投入後 1°CZ分の割合で昇温した。そして、到達温度 として 4つの温度 650、 700、 750、 800°Cを設定し、炉内温度が設定温度に到達し た後、それぞれの温度を 4時間維持した。ゥエー八にこのような核発生熱処理を施し た後に、 780°Cで 3時間 + 1000°Cで 16時間の成長熱処理をカ卩えて、酸素析出物の 密度と酸素析出量を測定した。
[0051] まず、酸素析出物の酸素析出量について実測値と演算値とを比較する。
測定方法及び演算方法は実施例 1と同じである。その結果を図 10に示す。 図 10は核発生熱処理を施した後に成長熱処理を加えた場合の酸素析出量の実測 値と演算値を示す図である。図 4一図 6と同様に、図 10は横軸が実測値、縦軸が演 算値とされており、同図において直線 Aに近い分布が得られるほど、実測値と演算値 の差が少ないといえる。図 10においては酸素析出量の演算値と実測値とがほぼ一 致し、直線 Aに近い分布が得られている。この結果から本発明による演算の精度が 高いことが分かる。
[0052] 次に、酸素析出物の密度について実測値と演算値とを比較する。
測定方法及び演算方法は実施例 1と同じである。その結果を図 11に示す。 図 11は核発生熱処理を施した後に成長熱処理を加えた場合の酸素析出物の密度 の実測値と演算値を示す図である。図 7—図 9と同様に、図 11は横軸が実測値、縦 軸が演算値とされており、同図において直線 Aに近い結果が得られるほど、実測値と 演算値の差が少なレ、といえる。図 11におレ、ては酸素析出物の密度の演算値と実測 値とがほぼ一致し、直線 Aに近い分布が得られている。この結果から本発明による演 算の精度が高レ、ことが分かる。
[0053] 本発明によれば、酸素析出物の密度を酸素析出量ともに優れた精度で予測するこ とができる。それは酸素析出物を構成する平均酸素原子数、すなわち酸素析出物の 平均サイズをも予測できることを意味する。
図面の簡単な説明 [図 1]図 1はシリコン単結晶の結晶成長過程における直胴位置での冷却曲線を示す 図である。
[図 2]図 2は演算によって求めたサーマルドナー濃度と実測によって求めたサーマル ドナー濃度とを示す図である。
[図 3]図 3は酸素析出物の核発生速度を求める式に用いられる係数と温度との対応 関係を示す図である。
[図 4]図 4は核発生熱処理を施した後に 800°Cで 4時間 + 1000°Cで 16時間の成長 熱処理を加えた場合の酸素析出量の実測値と演算値を示す図である。
[図 5]図 4は核発生熱処理を施した後に 900°Cで 4時間 + 1000°Cで 16時間の成長 熱処理を加えた場合の酸素析出量の実測値と演算値を示す図である。
[図 6]図 6は核発生熱処理を施した後に 1000°Cで 16時間の成長熱処理を加えた場 合の酸素析出量の実測値と演算値を示す図である。
[図 7]図 7は核発生熱処理を施した後に 800°Cで 4時間 + 1000°Cで 16時間の成長 熱処理を加えた場合の酸素析出物の密度と演算値を示す図である。
[図 8]図 8は核発生熱処理を施した後に 900°Cで 4時間 + 1000°Cで 16時間の成長 熱処理を加えた場合の酸素析出物の密度の実測値と演算値を示す図である。
[図 9]図 9は核発生熱処理を施した後に 1000°Cで 16時間の成長熱処理を加えた場 合の酸素析出物の密度の実測値と演算値を示す図である。
[図 10]図 10は核発生熱処理を施した後に 780°Cで 3時間 + 1000°Cで 16時間の成 長熱処理を加えた場合の酸素析出量の実測値と演算値を示す図である。
[図 11]図 11は核発生熱処理を施した後に 780°Cで 3時間 + 1000°Cで 16時間の成 長熱処理を加えた場合の酸素析出物の密度の実測値と演算値を示す図である。

Claims

請求の範囲
[1] 熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動を予測するシリコ ン単結晶中の酸素析出挙動予測方法において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、シリコン単結 晶に加える熱処理条件と、をパラメータとして演算を行い、シリコン単結晶を熱処理し た場合にシリコン単結晶中に発生する酸素析出物の密度および析出量を求めること を特徴とするシリコン単結晶中の酸素析出挙動予測方法。
[2] 熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動を予測するシリコ ン単結晶中の酸素析出挙動予測方法において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、をパラメータ として演算を行い、シリコン単結晶を熱処理した場合にその熱処理過程でシリコン単 結晶中に発生する酸素析出物の核発生速度を求め、求めた核発生速度を用いて酸 素析出物の密度および析出量を求めること
を特徴とするシリコン単結晶中の酸素析出挙動予測方法。
[3] シリコン単結晶中の酸素析出物の密度および析出量を所望の値にするように、請 求項 1又は請求項 2記載のシリコン単結晶中の酸素析出挙動予測方法を利用して、 シリコン単結晶中の初期酸素濃度と、シリコン単結晶インゴットが結晶成長の際に受 ける 400°Cから 550°Cの間の熱履歴と、シリコン単結晶に加える熱処理条件と、を決 定すること
を特徴とするシリコン単結晶の製造パラメータ決定方法。
[4] 熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動をコンピュータで 予測するためのプログラムを記憶する記憶媒体において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、シリコン単結 晶に加える熱処理条件と、をパラメータとして演算を行う処理と、
シリコン単結晶を熱処理した場合にシリコン単結晶中に発生する酸素析出物の密 度および析出量を求める処理と、
をプログラムとして記憶すること
を特徴とするシリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶 媒体。
熱処理に応じてシリコン単結晶中に発生する酸素析出物の挙動をコンピュータで 予測するためのプログラムを記憶する記憶媒体において、
シリコン単結晶中の初期酸素濃度と、シリコン単結晶が結晶成長の際に受ける 400 °Cから 550°Cの間の熱履歴に応じて発生するサーマルドナーの濃度と、をパラメータ として演算を行う処理と、
シリコン単結晶を熱処理した場合にその熱処理過程でシリコン単結晶中に発生す る酸素析出物の核発生速度を求め、求めた核発生速度を用いて酸素析出物の密度 および析出量を求める処理と、
をプログラムとして記憶すること
を特徴とするシリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶 媒体。
PCT/JP2005/001114 2004-01-27 2005-01-27 シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶媒体 WO2005071144A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2005517323A JP4901217B2 (ja) 2004-01-27 2005-01-27 シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶媒体
US10/586,445 US8246744B2 (en) 2004-01-27 2005-01-27 Method for predicting precipitation behavior of oxygen in silicon single crystal, method for determining production parameter of silicon single crystal, and storage medium for storing program for predicting precipitation behavior of oxygen in silicon single crystal

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2004-018313 2004-01-27
JP2004018313 2004-01-27

Publications (1)

Publication Number Publication Date
WO2005071144A1 true WO2005071144A1 (ja) 2005-08-04

Family

ID=34805560

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/001114 WO2005071144A1 (ja) 2004-01-27 2005-01-27 シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶媒体

Country Status (4)

Country Link
US (1) US8246744B2 (ja)
JP (1) JP4901217B2 (ja)
TW (1) TWI290182B (ja)
WO (1) WO2005071144A1 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008038786A1 (fr) 2006-09-29 2008-04-03 Sumco Techxiv Corporation Procédé de traitement thermique de plaquettes en silicium
WO2013084410A1 (ja) * 2011-12-06 2013-06-13 信越半導体株式会社 シリコン単結晶の評価方法およびシリコン単結晶の製造方法
WO2018012019A1 (ja) * 2016-07-11 2018-01-18 株式会社Sumco シリコンウェーハの評価方法及び製造方法
CN108133072A (zh) * 2017-11-17 2018-06-08 上海同继地质科技有限公司 云母Ar-Ar年代学热史模拟方法及系统
WO2018235678A1 (ja) * 2017-06-23 2018-12-27 株式会社Sumco シリコンウェーハのサーマルドナー生成挙動予測方法、シリコンウェーハの評価方法およびシリコンウェーハの製造方法
JP2019019030A (ja) * 2017-07-18 2019-02-07 信越半導体株式会社 シリコン単結晶の評価方法およびシリコン単結晶の製造方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2997096B1 (fr) * 2012-10-23 2014-11-28 Commissariat Energie Atomique Procede de formation d'un lingot en silicium de resistivite uniforme
US20150294868A1 (en) * 2014-04-15 2015-10-15 Infineon Technologies Ag Method of Manufacturing Semiconductor Devices Containing Chalcogen Atoms
KR101870702B1 (ko) 2017-01-16 2018-06-25 에스케이실트론 주식회사 실리콘 단결정 내의 산소 석출물 예측 방법
FR3075379B1 (fr) * 2017-12-15 2019-11-22 Commissariat A L'energie Atomique Et Aux Energies Alternatives Methode de validation de l'histoire thermique d'un lingot semi-conducteur
EP3929335A1 (en) * 2020-06-25 2021-12-29 Siltronic AG Semiconductor wafer made of single-crystal silicon and process for the production thereof

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04130732A (ja) * 1990-09-21 1992-05-01 Komatsu Electron Metals Co Ltd 半導体ウェーハの熱処理方法
JPH04298042A (ja) * 1991-03-27 1992-10-21 Komatsu Electron Metals Co Ltd 半導体の熱処理方法
JPH11147789A (ja) * 1997-11-11 1999-06-02 Shin Etsu Handotai Co Ltd シリコン単結晶中の酸素析出挙動を割り出す方法、およびシリコン単結晶ウエーハ製造工程の決定方法、並びにシリコン単結晶中の酸素析出挙動を割り出すためのプログラムを記録した記録媒体

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02263792A (ja) 1989-03-31 1990-10-26 Shin Etsu Handotai Co Ltd シリコンの熱処理方法
JPH04174300A (ja) 1990-11-06 1992-06-22 Mitsubishi Electric Corp 誘導飛しょう体
JPH05102167A (ja) 1991-10-07 1993-04-23 Shin Etsu Handotai Co Ltd シリコンの熱処理方法
JP3279527B2 (ja) 1998-08-24 2002-04-30 住友金属工業株式会社 半導体シリコン基板におけるig能の評価方法、及び半導体シリコン基板の製造方法
JP3874255B2 (ja) 2002-02-28 2007-01-31 信越半導体株式会社 シリコンウェーハ中のbmdサイズの評価方法
JP2003318181A (ja) 2002-04-25 2003-11-07 Sumitomo Mitsubishi Silicon Corp 半導体シリコン基板におけるig能の評価方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04130732A (ja) * 1990-09-21 1992-05-01 Komatsu Electron Metals Co Ltd 半導体ウェーハの熱処理方法
JPH04298042A (ja) * 1991-03-27 1992-10-21 Komatsu Electron Metals Co Ltd 半導体の熱処理方法
JPH11147789A (ja) * 1997-11-11 1999-06-02 Shin Etsu Handotai Co Ltd シリコン単結晶中の酸素析出挙動を割り出す方法、およびシリコン単結晶ウエーハ製造工程の決定方法、並びにシリコン単結晶中の酸素析出挙動を割り出すためのプログラムを記録した記録媒体

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
TAKENO H. ET AL: "Practical Computer Simulation Technique to Predict Oxygen Precipitation Behavior in Czochralski Silicon Wafers for Various Thermal Processes", JOURNAL OF ELECTROCHEMICAL SOCIETY, vol. 144, no. 12, 1997, pages 4340 - 4345, XP002094508 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008038786A1 (fr) 2006-09-29 2008-04-03 Sumco Techxiv Corporation Procédé de traitement thermique de plaquettes en silicium
US8573969B2 (en) 2006-09-29 2013-11-05 Sumco Techxiv Corporation Silicon wafer heat treatment method
WO2013084410A1 (ja) * 2011-12-06 2013-06-13 信越半導体株式会社 シリコン単結晶の評価方法およびシリコン単結晶の製造方法
JP2013119486A (ja) * 2011-12-06 2013-06-17 Shin Etsu Handotai Co Ltd シリコン単結晶の評価方法およびシリコン単結晶の製造方法
US9111883B2 (en) 2011-12-06 2015-08-18 Shin-Etsu Handotai Co., Ltd. Method for evaluating silicon single crystal and method for manufacturing silicon single crystal
WO2018012019A1 (ja) * 2016-07-11 2018-01-18 株式会社Sumco シリコンウェーハの評価方法及び製造方法
JP2018008832A (ja) * 2016-07-11 2018-01-18 株式会社Sumco シリコンウェーハの評価方法及び製造方法
WO2018235678A1 (ja) * 2017-06-23 2018-12-27 株式会社Sumco シリコンウェーハのサーマルドナー生成挙動予測方法、シリコンウェーハの評価方法およびシリコンウェーハの製造方法
US11121003B2 (en) 2017-06-23 2021-09-14 Sumco Corporation Method of predicting thermal donor formation behavior in silicon wafer, method of evaluating silicon wafer, and method of producing silicon wafer
JP2019019030A (ja) * 2017-07-18 2019-02-07 信越半導体株式会社 シリコン単結晶の評価方法およびシリコン単結晶の製造方法
CN108133072A (zh) * 2017-11-17 2018-06-08 上海同继地质科技有限公司 云母Ar-Ar年代学热史模拟方法及系统
CN108133072B (zh) * 2017-11-17 2021-08-31 上海同继地质科技有限公司 云母Ar-Ar年代学热史模拟方法及系统

Also Published As

Publication number Publication date
US20090210166A1 (en) 2009-08-20
JPWO2005071144A1 (ja) 2007-09-06
US8246744B2 (en) 2012-08-21
JP4901217B2 (ja) 2012-03-21
TW200526822A (en) 2005-08-16
TWI290182B (en) 2007-11-21

Similar Documents

Publication Publication Date Title
WO2005071144A1 (ja) シリコン単結晶中の酸素析出挙動予測方法、シリコン単結晶の製造パラメータ決定方法、シリコン単結晶中の酸素析出挙動予測用プログラムを記憶する記憶媒体
Vinci et al. Thermal strain and stress in copper thin films
Kulkarni et al. Quantification of defect dynamics in unsteady-state and steady-state Czochralski growth of monocrystalline silicon
JP2023502716A (ja) インサイチュで収集された情報に基づく熱処理の調整制御方法及びその応用
JP3279527B2 (ja) 半導体シリコン基板におけるig能の評価方法、及び半導体シリコン基板の製造方法
CN111479957B (zh) 用于验证半导体晶锭热史的方法
Brambilla et al. Solid–liquid interface velocity and diffusivity in laser-melt amorphous silicon
KR100912345B1 (ko) 단결정 성장 공정 파라미터를 이용한 산소농도 예측방법 및그 프로그램이 기록된 기록매체
TWI307924B (ja)
US6803242B2 (en) Evaluation method of IG effectivity in semiconductor silicon substrates
Huang et al. Evaluating cooling tube effects on defect formation and distribution in Czochralski crystal growth
Riley et al. Influence of microstructure and phase morphology on the stability of high temperature irradiation resistant thermocouples
Perry et al. An In Situ Resistance-Based Method for Tracking the Temporal Evolution of Recovery and Recrystallization in Ni-Base Single-Crystal Superalloy at Super-Solvus Temperatures
Yin et al. Improved Mechanical and Electrical Properties of Phase-Change-Memory Ga–Sb Thin Films
JP2010083712A (ja) 結晶欠陥状態予測方法、シリコンウェーハの製造方法
JPS61219795A (ja) 析出核の形成速度が速いシリコン単結晶ウエハおよびその製造法
JP7437165B2 (ja) エピタキシャルシリコンウェーハにおける欠陥の発生予測方法およびエピタキシャルシリコンウェーハの製造方法
WO2022172367A1 (ja) エピタキシャルシリコンウェーハにおける欠陥の発生予測方法およびエピタキシャルシリコンウェーハの製造方法
JP3073537B2 (ja) 金属薄膜形成条件の制御方法
Jakubczyk Phase transitions in Co78Si9B13 and Fe78Si9B13 metallic glasses induced by isochronal annealing
Mikhin et al. Effect of Crystal Quality on HCP-BCC Phase Transition in Solid 4 He
JP3508541B2 (ja) シリコン基板の熱処理条件を設定する方法、およびシリコン基板を熱処理する方法、並びにシリコン基板の製造方法
Molodov Motion of grain boundaries: experiments on bicrystals
JP5072185B2 (ja) 熱処理ボートの評価方法
Lang et al. Vacancy Type Defects and Related Phenomena in Face-Centered Cubic Copper-Zinc Alloys. EUR 4622.

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2005517323

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 10586445

Country of ref document: US

122 Ep: pct application non-entry in european phase