FR3075379A1 - Methode de validation de l'histoire thermique d'un lingot semi-conducteur - Google Patents

Methode de validation de l'histoire thermique d'un lingot semi-conducteur Download PDF

Info

Publication number
FR3075379A1
FR3075379A1 FR1762282A FR1762282A FR3075379A1 FR 3075379 A1 FR3075379 A1 FR 3075379A1 FR 1762282 A FR1762282 A FR 1762282A FR 1762282 A FR1762282 A FR 1762282A FR 3075379 A1 FR3075379 A1 FR 3075379A1
Authority
FR
France
Prior art keywords
thermal
concentration
thermal donors
ingot
time step
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
FR1762282A
Other languages
English (en)
Other versions
FR3075379B1 (fr
Inventor
Jordi Veirman
Wilfried Favre
Elenore Letty
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.)
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Commissariat a lEnergie Atomique CEA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
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
Priority to FR1762282A priority Critical patent/FR3075379B1/fr
Application filed by Commissariat a lEnergie Atomique CEA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Commissariat a lEnergie Atomique CEA
Priority to PCT/EP2018/084450 priority patent/WO2019115574A1/fr
Priority to CA3085769A priority patent/CA3085769A1/fr
Priority to US16/772,613 priority patent/US11359305B2/en
Priority to EP18814626.0A priority patent/EP3724377A1/fr
Priority to KR1020207020082A priority patent/KR20200098612A/ko
Priority to CN201880080782.0A priority patent/CN111479957B/zh
Publication of FR3075379A1 publication Critical patent/FR3075379A1/fr
Application granted granted Critical
Publication of FR3075379B1 publication Critical patent/FR3075379B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Materials Engineering (AREA)
  • Metallurgy (AREA)
  • Organic Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Crystals, And After-Treatments Of Crystals (AREA)
  • Testing Or Measuring Of Semiconductors Or The Like (AREA)

Abstract

L'invention concerne une méthode de validation expérimentale d'une histoire thermique d'un lingot semi-conducteur obtenue par simulation d'un procédé de cristallisation. Cette méthode comprend les étapes suivantes : a) mesurer la concentration en oxygène interstitiel dans une portion du lingot semi-conducteur ; b) calculer une valeur théorique de la concentration en donneurs thermiques formés lors du procédé de cristallisation, à partir de la mesure de la concentration en oxygène interstitiel et à partir de l'histoire thermique dans la portion du lingot semi-conducteur ; c) mesurer une valeur expérimentale de la concentration en donneurs thermiques dans la portion du lingot semi-conducteur ; et d) comparer les valeurs théorique et expérimentale de la concentration en donneurs thermiques.

Description

MÉTHODE DE VALIDATION DE L’HISTOIRE THERMIQUE D’UN
LINGOT SEMI-CONDUCTEUR
DOMAINE TECHNIQUE
La présente invention est relative aux procédés de cristallisation de lingots semiconducteurs et concerne plus particulièrement une méthode ou procédé pour valider expérimentalement l'histoire thermique d’un lingot semi-conducteur obtenue par simulation d’un procédé de cristallisation.
ÉTAT DE LA TECHNIQUE
Les cellules photovoltaïques à haut rendement sont pour la plupart fabriquées à partir de plaquettes issues d’un lingot de silicium monocristallin Czochralski (Cz). Bien que réputé pour ses performances électroniques élevées, notamment en termes de durée de vie des porteurs de charge, le silicium CZ n’est pas exempt de défauts et d’impuretés. L’oxygène constitue la principale impureté du silicium CZ et se présente sous la forme d’atomes en positions interstitielles dans le réseau cristallin. Parmi les défauts du silicium CZ, on peut notamment citer les précipités d’oxygène, les lacunes d’oxygène et les donneurs thermiques. Les donneurs thermiques sont des agglomérats (à base d’oxygène) qui se forment à des températures comprises entre 350 °C et 550 °C et qui affectent les propriétés éèctriques du matériau en créant des électrons libres.
Afin d’optimiser à la fois la cadence de production et la qualité des lingots de silicium CZ, les fabricants de lingots ont recours à des outils informatiques permettant de simuler le procédé de cristallisation Czochralski. Grâce à ces simulations, il est par exemple possible de connaître la température dans chaque portion du lingot à un instant donné de la cristallisation, ainsi que son évolution au cours de la cristallisation. Cette évolution, appelée communément « histoire thermique », a une grande influence sur la quantité des défauts formés au cours de la cristallisation dans la portion considérée du lingot. L’histoire thermique du lingot varie en fonction de nombreux paramètres, tels que la hauteur relative (appelée aussi « fraction solidifiée >>) à laquelle se situe la portion considérée du lingot, la géométrie et les matériaux des pièces qui constituent le four de cristallisation, la vitesse de tirage du lingot et la puissance délivrée par les résistances du four.
La simulation thermique d’un lingot fait appel à des modèles physiques complexes. Ces modèles peuvent conduire à des valeurs erronées de l'histoire thermique si, par exemple, l’algorithme de calcul ne converge pas suffisamment ou la définition du maillage est trop faible.
Pour s’assurer de la pertinence des résultats de simulation, et le cas échéant mesurer le degré de précision des modèles, il est nécessaire d’effectuer une validation expérimentale en confrontant l’histoire thermique calculée à des mesures de la température du lingot. Il est cependant difficile de mesurer l’histoire thermique d’un lingot en cours de fabrication, du fait de sa rotation, de son déplacement en translation ainsi que des très hautes températures qui régnent dans l’enceinte du four (>1400 °C).
Le document [« Thermal simulation of the Czochralski silicon grow process by three different models and comparison with experimental results », E. Dornberger et al., Journal of Crystal Growth 180, pp.461-467, 1997] décrit une méthode pour valider expérimentalement l’histoire thermique d’un lingot de silicium CZ. Plusieurs lingots de longueurs différentes sont d’abord tirés en suivant une recette de tirage donnée. Puis, les lingots sont équipés de thermocouples et montés successivement dans un four de cristallisation. Chaque lingot est mis en contact avec un bain de silicium en fusion après que le four ait atteint une température de consigne, à laquelle le tirage du lingot a habituellement lieu. La température dans le lingot est ensuite mesurée au moyen des thermocouples après qu’un état d’équilibre ait été atteint.
Cette méthode de validation de l’histoire thermique est fastidieuse car elle nécessite de pré-fabriquer, pour chaque recette de tirage, des lingots de longueurs différentes (ces longueurs correspondant à différents instants du procédé de cristallisation simulé). Elle n’est en outre pas représentative du procédé de cristallisation Czochralski, car le lingot n’est pas entraîné en rotation lors de l’étape de mesure des températures (à cause du câblage des thermocouples). Enfin, comme aucune cristallisation n’a lieu lors de l’étape de mesure, le lingot ne dissipe pas de chaleur latente de solidification, ce qui peut grandement influencer les températures mesurées dans le lingot. L’histoire thermique calculée et les mesures de température sont alors difficilement comparables.
RÉSUMÉ DE L’INVENTION
Il apparaît ainsi le besoin d’une méthode simple et rapide à mettre en œuvre pour valider (ou invalider) avec certitude une histoire thermique d’un lingot semi-conducteur, cette histoire thermique ayant été obtenue par simulation d’un procédé de cristallisation.
Selon l’invention, on tend à satisfaire ce besoin en prévoyant une méthode de validation comprenant les étapes suivantes :
a) mesurer la concentration en oxygène interstitiel dans une portion du lingot semi-conducteur ;
b) calculer une valeur théorique de la concentration en donneurs thermiques formés lors du procédé de cristallisation, à partir de la mesure de la concentration en oxygène interstitiel et à partir de l'histoire thermique dans la portion du lingot semi-conducteur ;
c) mesurer une valeur expérimentale de la concentration en donneurs thermiques dans la portion du lingot semi-conducteur ; et
d) comparer les valeurs théorique et expérimentale de la concentration en donneurs thermiques.
Les donneurs thermiques sont des défauts à base d’oxygène qui se forment entre 350 °C et 500 °C au cours de la cristallisation etdont la concentration constitue un marqueur de l'histoire thermique du lingot. Ainsi, plutôt que de confronter directement des valeurs de la température dans le lingot, la méthode de validation selon l’invention effectue une comparaison entre des valeurs théorique et expérimentale de la concentration en donneurs thermiques, la valeur théorique de la concentration en donneurs thermiques étant dérivée de l’histoire thermique simulée.
Contrairement à la méthode de l’art antérieur, la méthode de validation selon l’invention ne requiert pas de mesures in situ, i.e. à l’intérieur du four de cristallisation. La concentration en oxygène interstitiel et la concentration en donneurs thermiques peuvent en effet être mesurées après que le lingot ait été extrait du four, et donc bien plus simplement que la température au cours de la cristallisation. La fabrication du lingot n’est donc pas impactée par la méthode de validation selon l’invention, ce qui permet d’obtenir des résultats fidèles au procédé de cristallisation employé.
Dans un mode de mise en œuvre préférentiel de l’invention, le calcul de la valeur théorique de la concentration en donneurs thermiques comporte les opérations suivantes :
- établir au moins une expression empirique décrivant l’histoire thermique (ou représentative de l’histoire thermique) dans la portion du lingot semiconducteur ;
- discrétiser en pas de temps successifs l’histoire thermique dans la portion du lingot semi-conducteur, en associant à chaque pas de temps une température calculée à l’aide de ladite au moins une expression empirique ;
- calculer pour chaque pas de temps dont la température associée est comprise entre 350 °C et 550 °C une teneur en donneurs therniques formés pendant ledit pas de temps ; et
- calculer la valeur théorique de la concentration en donneurs thermiques à partir des teneurs en donneurs thermiques formés pendant les pas de temps successifs.
Ce mode de calcul permet de déterminer avec précision la concentration en donneurs thermiques formés lors du procédé de cristallisation (entre 350 °C et 550 °C), en tenant compte du profil exact de la température dans la portion du lingot semi-conducteur.
Selon un développement de ce mode de mise en œuvre préférentiel, le calcul de la valeur théorique de la concentration en donneurs thermiques est réalisé par itérations et comporte, pour chaque pas de temps, les opérations suivantes :
- sommer la teneur en donneurs thermiques formés pendant ledit pas de temps et une teneur en donneurs thermiques présents dans le matériau semi5 conducteur au pas de temps précédent ;
- comparer ladite somme à une concentration maximale de donneurs thermiques à la température associée audit pas de temps ;
- fixer une teneur en donneurs thermiques présents dans le matériau semiconducteur audit pas de temps égale à ladite somme lorsque celle-ci est inférieure à la concentration maximale de donneurs thermiques ; et
- fixer la teneur en donneurs thermiques présents dans le matériau semiconducteur audit pas de temps égale à la concentration maximale de donneurs thermiques lorsque ladite somme est supérieure à la concentration maximale de donneurs thermiques.
De préférence, la teneur A[DT]n en donneurs thermiques formés pendant ledit pas de temps ΔΤη est calculée à l’aide de la relation suivante :
à[DT]n = ΔΤη* a* DiiTj * [0f]4 * m(TnY2 où a est une constante, Tn la température associée audit pas de temps ΔΤη, Di(Tn) le coefficient de diffusion de l’oxygène à la température Tn et m(Tn) la concentration en électrons libres à la température Tn.
L’histoire thermique dans la portion du lingot semi-conducteur est avantageusement décrite au moyen de deux polynômes du second degré.
Les étapes a) à d) de la méthode de validation selon l’invention peuvent être accomplies pour différentes portions réparties le long du lingot semi-conducteur, afin de vérifier la simulation thermique de l’ensemble du lingot, et pas seulement d’une portion. Ces différentes portions incluent de préférence l’extrémité haute, appelée tête, et l’extrémité basse, appelée queue, du lingot semi-conducteur. Elles sont avantageusement en nombre supérieur ou égal à 5.
BRÈVE DESCRIPTION DES FIGURES
D'autres caractéristiques et avantages de l'invention ressortiront clairement de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles :
- la figure 1 représente des étapes S1 à S4 d’une méthode de validation d’une histoire thermique selon l’invention ; et
- la figure 2 montre un exemple d’histoire thermique d’une portion de lingot semiconducteur, obtenu par simulation du procédé de cristallisation du lingot, ainsi que les résultats de l’étape l’ajustement de cette histoire thermique par deux fonctions polynomiales.
Pour plus de clarté, les éléments identiques ou similaires sont repérés par des signes de référence identiques sur l’ensemble des figures.
DESCRIPTION DÉTAILLÉE D’AU MOINS UN MODE DE RÉALISATION
Dans la description qui suit, on appelle « histoire thermique >> l’évolution de la température d’une portion d’un lingot semi-conducteur au cours de la cristallisation du lingot. Une portion correspond de préférence à une tranche du lingot semi-conducteur orientée perpendiculairement à l’axe longitudinal (ou axe de tirage) du lingot et dont l’épaisseur peut être variable suivant sa position dans le lingot. La position d’une portion ou tranche du lingot est appelée « hauteur relative >> et s’exprime généralement en un pourcentage de la hauteur totale du lingot.
Chacune des portions du lingot possède une histoire thermique spécifique, qu’il est possible de calculer par simulation du procédé de cristallisation. L’ensemble de ces histoires thermiques permet de reconstituer l’évolution du champ de températures dans le lingot semi-conducteur.
La méthode de validation décrite ci-dessous permet de savoir si le calcul de l'histoire thermique pour au moins une portion du lingot semi-conducteur est exact et, le cas échéant, de connaître son degré de précision. Grâce à cette information, il est ensuite possible d’optimiser le procédé de cristallisation en se fiant aux simulations thermiques, d’affiner le modèle physique utilisé lors de ces simulations pour améliorer la précision du calcul de l’histoire thermique ou, dans le cas d’une histoire thermique manifestement erronée, de modifier en profondeur le modèle physique voire de changer de modèle.
Le lingot semi-conducteur est par exemple un lingot en silicium monocristallin obtenu par le procédé de cristallisation Czochralski (aussi appelé silicium CZ).
En référence à la figure 1, la méthode de validation comporte tout d'abord une étape S1 de mesure de la concentration en oxygène interstitiel [Oi] dans la portion du lingot dont on souhaite valider l'histoire thermique.
La mesure de la concentration en oxygène interstitiel [Oi] peut être effectuée en un point, par exemple par spectroscopie infrarouge à transformée de Fourier (FTIR) sur une plaquette épaisse (typiquement entre 1 mm et 2 mm d’épaisseur) prélevée dans la portion du lingot et dont la surface a été polie.
Dans un mode de mise en œuvre particulier de l’étape S1, la concentration en oxygène interstitiel [Oi] est mesurée sur le lingot entier, c’est-à-dire sans découpe préalable de plaquettes. La concentration [Oi] peut être mesurée à l’échelle du lingot par une technique de spectroscopie infrarouge appelée communément « Whole-rod FTIR ». Cette technique, dérivée de la spectroscopie infrarouge à transformée de Fourier (FTIR), comporte l’exposition de la portion du lingot à un faisceau infrarouge. L’absorption du faisceau infrarouge par la portion du lingot permet de déterminer une concentration en oxygène interstitiel moyennée sur le diamètre du lingot.
Une autre technique, basée sur la formation de donneurs thermiques, permet de déterminer la concentration en oxygène [Oi] dans le silicium CZ. Cette technique a été décrite en détail dans les brevets FR2964459 et FR3009380 pour le cas d’une plaquette de silicium. Elle peut aussi être appliquée à l’échelle du lingot.
La résistivité électrique initiale est d’abord mesurée dans la portion du lingot afin de déterminer les concentrations en dopants accepteurs et/ou donneurs. Le lingot est ensuite soumis à un recuit de manière à créer de nouveaux donneurs thermiques, différents de ceux formés lors de la cristallisation. La température de ce recuit est de préférence constante et comprise entre 350 °C et 550 °C. Puis, la résistivité électrique après recuit est mesurée dans la même zone du lingot. À partir de cette deuxième valeur de résistivité et des concentrations en dopants, il est possible de calculer la concentration en donneurs thermiques formés par le recuit. Enfin, la concentration en oxygène interstitiel [Oi] dans la portion du lingot est déterminée à partir de la concentration de donneurs thermiques nouvellement créés et de la durée du recuit entre 350 °C et 550 °C, par exemple au moyen d’un toaque.
Cette dernière technique est précise et particulièrement simple à mettre en œuvre. Elle est avantageuse même lorsqu’elle est appliquée à une plaquette prélevée dans le lingot car, contrairement à la technique FTIR, elle ne requiert pas le polissage de la plaquette et n’est pas limitée en termes d’épaisseur.
La méthode de validation de la figure 1 comporte ensuite une étape S2 permettant de déterminer une valeur théorique [DT]th de la concentration en donneurs thermiques formés dans la portion du lingot, au cours de la cristallisation. Cette valeur théorique [DT]th est calculée à partir de l'histoire thermique de la portion du lingot et à partir de la concentration en oxygène interstitiel [Oi] mesurée à l’étape S1, en utilisant un modèle mathématique qui décrit la cinétique de formation des donneurs thermiques.
Dans un mode de mise en œuvre préférentiel de la méthode de validation, le calcul de la concentration en donneurs thermiques [DT]th se compose des opérations suivantes.
Tout d’abord, l'histoire thermique de la portion du lingot, représentée graphiquement par une courbe de la température T en fonction de la durée de cristallisation t, est décrite analytiquement à l’aide d’une ou plusieurs expressions empiriques. Cette première opération, appelée « ajustement >> (ou « fit >> en anglais), consiste à déterminer au moins une fonction T(t) dont la courbe reproduit l’histoire thermique (i.e. le profil de température) simulée.
Les expressions empiriques décrivant ou représentant l’histoire thermique sont avantageusement choisies parmi des polynômes de degré n (n étant un entier naturel non nul), car les fonctions polynomiales permettent un ajustement précis des courbes d’histoire thermique d’un lingot semi-conducteur.
Pour améliorer la précision du calcul de la concentration en donneurs thermiques, il est préférable que la courbe d’histoire thermique comporte un nombre élevé de points, par exemple entre 300 et 600 points répartis sur une large plage de température (typiquement de 1400°C à 200 °C) dont au moins 20 psints situés dans la plage de températures 350°C-550°C.
La figure 2 montre à titre d’exemple l’histoire thermique d’une portion d’un lingot de silicium CZ située à une hauteur relative de 1 %. Cet exemple d’histoire thermique a été obtenu par simulation du procédé de cristallisation Czochralski, à l’aide de l’outil de simulation ANSYS® Fluent®. Il permet de distinguer deux phases successives du procédé Czochralski : la phase de solidification du lingot et la phase de refroidissement du lingot. La fin de la phase de solidification, à t « 35 h sur le graphe de la figure 2, est marquée par le retrait du lingot du bain de silicium en fusion. Le début de la phase de refroidissement est marqué par une diminution plus rapide de la température.
L’histoire thermique est ici décrite au moyen de deux polynômes du second degré (T(t) = A*t2+B*t + C, avec A, B et C des paramètres d’ajustement constants), l’un pour la phase de solidification, l’autre pour la phase de refroidissement. L’utilisation de deux polynômes de second degré permet de limiter le temps et les moyens de calcul nécessaires à l’opération d’ajustement et s’avère particulièrement adaptée à ce type d’histoire thermique en deux parties.
La deuxième opération du calcul de la valeur théorique [DTjth est une discrétisation de l’histoire thermique en pas de temps Atn successifs. Cette discrétisation peut être réalisée, par exemple à l’aide d’un tableur, en associant à chaque pas de temps Atn une valeur de température Tn calculée à l’aide de l’expression empirique déterminée lors de l’opération précédente. La durée des pas de temps Atn est choisie principalement en fonction de la longueur du lingot et de la vitesse de tirage lors du procédé de cristallisation. Elle est de préférence inférieure à 10 min, par exemple égale à 1 min.
Lorsque l’histoire thermique est décrite par plusieurs expressions, chaque calcul de la température Tn est effectué avec l’expression empirique associée au pas de temps Atn considéré. Dans l’exemple de la figure 2, le premier polynôme (correspondant à la phase de solidification) est utilisé pour les 35*60/Atn(min) premiers pas de temps et le deuxième polynôme est utilisé pour tous les pas de temps suivants.
Ensuite, au moins pour chaque pas de temps Atn dont la température associée Tn est comprise entre 350 °C et 550 °C, on calcule une teieur (i.e. concentration) A[DT]n en donneurs thermiques formés pendant le pas de temps Atn. La teneur A[DT]n en donneurs thermiques est de préférence calculée à l’aide de la relation suivante :
A[DT]n = ΔΤη * a * Dt(Tn~) * [OJ4 * m(Tnr2 (1) où a est une constante, Tn la température associée au pas de temps ΔΤη, Di(Tn) le coefficient de diffusion de l’oxygène à la température Tn et m(Tn) la concentration en électrons libres à la température Tn.
Cette relation (1) est tirée d’un article intitulé [« Unified model for formation kinetics of oxygen thermal donors in silicon », K. Wada, Physical Review B, Vol.30, N. 10, pp. 5885-5895,1984], qui décrit un modèle pour calculer la cinétique de formation des donneurs thermiques dans le silicium CZ à une température constante, par exemple de 450 °C. Ce document, dont le contenu est incorpaé ici à titre de référence, fournit également les valeurs des paramètres de la relation (1) ou les formules permettant de les calculer, à l’exception de la masse effective des trous. Cette dernière est d’après la littérature égale à 0,81 pour le silicium.
La discrétisation en pas de temps rend possible l’utilisation du modèle de Wada, car elle fournit une température Tn constante pour chaque pas de temps.
En pratique, la teneur A[DT]n en donneurs thermiques peut être calculée pour tous les pas de temps de l’histoire thermique simulée, par exemple de 1414 °C à la température ambiante (i.e. 25 °C), mais la contribution des températures en dehors de la plage 350 °C-550 °C sur la formation des donneurs thermiqies est négligeable.
Enfin, lors d’une dernière opération, on calcule la valeur théorique [DT]th de la concentration en donneurs thermiques à partir des teneurs A[DT]n en donneurs thermiques formés pendant les différents pas de temps.
Dans un mode de mise en oeuvre simplifié de l’étape S2, la somme de toutes les teneurs A[DT]n calculées constitue la valeur théorique [DT]th de la concentration en donneurs thermiques formés au cours de la cristallisation :
[on« =^A[DT]n
Toutefois, d’après l’article susmentionné, la concentration en donneurs thermiques dans le lingot à une température donnée ne peut excéder une concentration maximale [DT]max de donneurs thermiques. Autrement dit, il existe une limite de « solubilité >> des donneurs thermiques dans le silicium CZ, qui dépend de la température. Cette concentration maximale s’écrit de la façon suivante :
[DT]maM = l * [Oi]3 * m(Tn)-2 (2) où b est une autre constante fournie par l’article de Wada.
Ainsi, dans un mode de mise en oeuvre particulier de l’étape S2, le calcul de valeur théorique [DT]th de la concentration en donneurs thermiques est réalisé par itérations, en tenant compte à chaque itération de la concentration maximale [DT]maxde donneurs thermiques.
Pour chaque pas de temps Atn, la teneur A[DT]n en donneurs thermiques formés pendant le pas de temps Atn est ajoutée à la teneur en donneurs thermiques présents dans le matériau semi-conducteur au pas de temps précédent, notée [DT]n-i. Puis, cette somme est comparée à la concentration maximale [DTjmax de donneurs thermiques à la température Tn (associée au pas de temps Atn considéré). Si la somme [DT]n-i+A[DT]n est supérieure à la concentration maximale [DTjmax, alors la teneur en donneurs thermiques au pas de temps Atn est considérée comme égale à la concentration maximale [DTjmax des donneurs thermiques (à la température Tn) :
[^Πη - lDT]max(Tn)
Si au contraire la somme [DT]n-i+A[DT]n est inférieure à la concentration maximale [DT]max, alors cette somme devient la teneur en donneurs thermiques présents dans le matériau au pas de temps Atn, soit :
[ΏΤ]η = [ΏΤ]η_1+Δ[ΏΤ]η
On procède ainsi pour tous les pas de temps (au moins ceux dont la température associée Tn est comprise entre 350 °C et 550 °C) et la valeurthéorique [DT]th de la concentration en donneurs thermiques est égale à la teneur en donneurs thermiques du matériau au dernier pas de temps. La teneur initiale en donneurs thermiques du matériau (i.e. avant le premier pas de temps) est supposée nulle ([DT]o = 0).
Cette variante de calcul procure davantage de précision sur le calcul de la valeur théorique [DT]th de la concentration en donneurs thermiques.
Le modèle de Wada n’est pas le seul modèle décrivant la cinétique de formation des donneurs thermiques et permettant de calculer la valeur théorique [DT]th de la concentration en donneurs thermiques formés au cours de la cristallisation. On peut notamment citer le modèle de Y. J. Lee et al., décrit dans l’article [« Simulation of the kinetics of oxygen complexes in crystalline silicon », Physical Review B, Vol.66, 165221,2002],
L’équation du modèle de Y. J. Lee et al. qui donne la teneur A[DT]n en donneurs thermiques pour chaque pas de temps est la suivante :
16 ârd[Ot] â[DT]„ = —^dt . at
fc=0 où [Ok] correspond à la concentration en donneurs thermiques de la famille k (0<k< 16) et dont la dérivée temporelle s’écrit :
J=1 + £(1 + Sk0(-kka'1 [Ot][O(] + k)'1 [Ot+I])
1=1 kJak J et kJdk J sont des constantes d’association et de dissociation pour la réaction de [Oj] et [Ok-j] en [Ok], dont les expressions sont fournis par l’article de Lee et al. ôki est le delta de Kronecker et permet d’éviter de compter deux fois une même réaction.
D’après des expériences menées par la demanderesse, il apparaît que le modèle de Wada tend à surestimer la cinétique de formation des donneurs thermiques par rapport à d’autres modèles plus précis. Afin de prendre en compte cette tendance, les teneurs A[DT]n en donneurs thermiques calculées à l’étape S2 sont avantageusement multipliées par un coefficient réducteur, par exemple de 0,75. Cette pondération des teneurs A[DT]n améliore grandement la précision du calcul de la valeur théorique [DT]th.
Pour la même raison, lorsque la concentration en oxygène interstitiel [Oi] est mesurée par FTIR ou sa variante à l’échelle du lingot à l’étape S1 de la méthode, il est préférable de choisir la même norme de mesure que celle utilisé dans le modèle de formation des donneurs thermiques sélectionné. Sinon, un coefficient pondérateur est avantageusement appliqué à la mesure de la concentration en oxygène interstitiel [Oi] avant qu’elle ne soit utilisée en entrée du modèle de Wada. Par exemple, si la technique de mesure obéit à la norme FTIR recommandée par SEMI avec un coefficient de calibration de 6,28 ppma.cm, la valeur mesurée de la concentration en oxygène interstitiel [Oi] est multipliée par 5,5/6,28, car le modèle de Wada utilise un coefficient de calibration de 5,5 ppma.cm environ. Si le modèle de Lee (coefficient de calibration de 6,28 ppma.cm) est employé et que la concentration en oxygène interstitiel [Oi] est mesurée avec la norme FTIR SEMI, aucun coefficient pondérateur n’est utilisé.
L’étape S3 de la méthode de validation (cf. Fig.1) consiste à déterminer une valeur expérimentale [DT]exp de la concentration en donneurs thermiques formés lors de la cristallisation du lingot.
La valeur expérimentale [DT]exp de la concentration en donneurs thermiques peut être obtenue à partir de la variation de résistivité ou de la variation de la concentration en porteurs de charge, provoquée par un recuit à haute température (> 600 °C). Ce recuit à haute température (typiquement 30 minutes à 650 °C) permet de détruire les donneurs thermiques formés lors de la cristallisation du lingot. La résistivité électrique peut être mesurée (avant et après le recuit de destruction) par la méthode des quatre pointes, la méthode de Van der Pauw ou être dérivée de la mesure du courant de Foucault. Cette technique de mesure est décrite en détail dans le brevet FR3009380, dont le contenu est incorporé par référence. La concentration en porteurs de charge peut être mesurée par effet Hall ou déduite de mesures C-V.
Le lingot est de préférence dopé de façon à présenter une résistivité initiale (i.e. après la cristallisation et avant tout traitement thermique) supérieure à 1 Ω.οιτι, de sorte que la variation de résistivité avant-après recuit de destruction des donneurs thermiques soit détectable avec précision.
L’étape S3 peut être mise en œuvre après l’étape S1 même lorsque la technique des brevets FR2964459 et FR3009380 est utilisée pour mesurer la concentration en oxygène interstitiel [Oi]. Il suffit dans ce cas de considérer la résistivité électrique (ou la concentration en porteurs de charge) initiale du lingot. L’étape S3 peut aussi être mise en œuvre avant l’étape S1, auquel cas il n’y a plus de donneurs thermiques au moment de réaliser le recuit entre 350 °C et 550 °C
En d’autres termes, la méthode de validation n’est limitée à aucun ordre des étapes S1 et S3. Le brevet FR3009380 donne plus de détails sur les façons d’articuler l’étape S1 de mesure de la concentration en oxygène interstitiel [Oi] et l’étape S3 de mesure de la concentration en donneurs thermiques [DT]exp.
Les étapes S1 et S3 peuvent être accomplies par un même équipement, par exemple l’équipement « OxyMap >> commercialisé par la société « AET Solar Tech >>.
Enfin, à l’étape S4 de la figure 1, la valeur théorique [DT]th et la valeur expérimentale [DT]exp de la concentration en donneurs thermiques sont comparées. Si la valeur théorique [DT]th est proche de la valeur expérimentale [DT]exp, typiquement comprise entre 0,7*[DT]exp et 1,3*[DT]exp, l’histoire thermique de la portion du lingot est considérée comme valide. Si au contraire la valeur théorique [DT]th est éloignée de la valeur expérimentale [DT]exp, typiquement > 1,3*[DT]exp, ou < 0,7*[DT]exp, alors l’histoire thermique de la portion du lingot n’est pas validée.
Comme les mesures de la concentration en oxygène et de la concentration en donneurs thermiques [DT]exp sont effectuées après la cristallisation du lingot, la méthode de validation décrite ci-dessus est plus rapide et plus simple à mettre en œuvre que la méthode de l’art antérieur. Les résultats obtenus sont en outre fidèles au procédé de cristallisation employé, le procédé Czochralski dans cet exemple, car la méthode de validation n’interfère pas avec la cristallisation du lingot.
De préférence, les étapes S1 à S4 de la figure 2 sont accomplies pour différentes portions réparties le long du lingot (suivant l’axe longitudinal du lingot), afin de vérifier la validité de la simulation thermique dans sa globalité. La simulation thermique est considérée comme juste après que l’histoire thermique de chacune des portions sélectionnées ait été validée à l’étape S4 de la méthode.
Les différentes portions du lingot sont de préférence en nombre supérieur ou égal à 5. Elles incluent avantageusement l’extrémité haute et l’extrémité basse, appelées respectivement « tête >> et « queue », car leurs histoires thermiques respectives sont très différentes. La concentration en oxygène interstitiel [Oi] et la concentration en donneurs thermiques [DTjexp sont de préférence mesurées au même endroit dans les différentes portions, par exemple sur un bord ou au centre des tranches issues du lingot.
Bien que la méthode de validation ait été décrite en relation avec un lingot en silicium monocristallin CZ, elle pourrait être appliquée à d’autres procédés de cristallisation et/ou matériaux semi-conducteurs (monocristallin ou polycristallin), dès lors que celui contient de l’oxygène. Le germanium et l’alliage silicium-germanium sont des candidats potentiels, car des donneurs thermiques à base d’oxygène sont également formés lors de leur cristallisation.

Claims (8)

  1. REVENDICATIONS
    1. Méthode de validation expérimentale d’une histoire thermique d’un lingot semiconducteur obtenue par simulation d’un procédé de cristallisation, comprenant les étapes suivantes :
    a) mesurer (S1) la concentration en oxygène interstitiel ([Oi]) dans une portion du lingot semi-conducteur ;
    b) calculer (S2) une valeur théorique ([DT» de la concentration en donneurs thermiques formés lors du procédé de cristallisation, à partir de la mesure de la concentration en oxygène interstitiel ([Oi]) et à partir de l’histoire thermique dans la portion du lingot semi-conducteur ;
    c) mesurer (S3) une valeur expérimentale ([DT]exp) de la concentration en donneurs thermiques dans la portion du lingot semi-conducteur ; et
    d) comparer (S4) les valeurs théorique ([DT» et expérimentale ([DT]exp) de la concentration en donneurs thermiques.
  2. 2. Méthode selon la revendication 1, dans laquelle le calcul de la valeur théorique ([DT» de la concentration en donneurs thermiques comporte les opérations suivantes :
    - établir au moins une expression empirique (T(t)) décrivant l’histoire thermique dans la portion du lingot semi-conducteur ;
    - discrétiser en pas de temps (Δ» successifs l’histoire thermique dans la portion du lingot semi-conducteur, en associant à chaque pas de temps (Δ» une température (Tn) calculée à l’aide de ladite au moins une expression empirique ;
    - calculer, pour chaque pas de temps dont la température (Tn) associée est comprise entre 350 °C et 550 °C, une teneur4[DT]n) en donneurs thermiques formés pendant ledit pas de temps ; et
    - calculer la valeur théorique ([DT» de la concentration en donneurs thermiques à partir des teneurs en donneurs thermiques formés pendant les pas de temps successifs.
  3. 3. Méthode selon la revendication 2, dans laquelle le calcul de la valeur théorique ([DT» de la concentration en donneurs thermiques est réalisé par itérations et comporte, pour chaque pas de temps (ΔΤη), les opérations suivantes :
    - sommer la teneur (A[DT]n) en donneurs thermiques formés pendant ledit pas de temps et une teneur ([DT]n-i) en donneurs thermiques présents dans le matériau semi-conducteur au pas de temps précédent ;
    - comparer ladite somme à une concentration maximale ([DT]max) de donneurs thermiques à la température (Tn) associée audit pas de temps (ΔΤη) ;
    - fixer une teneur ([DT]n) en donneurs thermiques présents dans le matériau semi-conducteur audit pas de temps égale à ladite somme lorsque celle-ci est inférieure à la concentration maximale de donneurs thermiques ; et
    - fixer la teneur ([DT]n) en donneurs thermiques présents dans le matériau semiconducteur audit pas de temps égale à la concentration maximale ([DT]max) de donneurs thermiques lorsque ladite somme est supérieure à la concentration maximale de donneurs thermiques.
  4. 4. Méthode selon l’une des revendications 2 et 3, dans laquelle la teneur Δ[ϋΤ]η en donneurs thermiques formés pendant ledit pas de temps ΔΤη est calculée à l’aide de la relation suivante :
    à[DT]n = ΔΤη* a* DiiTj * [0f]4 * m(TnY2 où a est une constante, Tn la température associée audit pas de temps ΔΤη, Di(Tn) le coefficient de diffusion de l’oxygène à la température Tn et m(Tn) la concentration en électrons libres à la température Tn.
  5. 5. Méthode selon l’une quelconque des revendications 2 à 4, dans lequel l’histoire thermique dans la portion du lingot semi-conducteur est décrite au moyen de deux polynômes du second degré.
  6. 6. Méthode selon l’une quelconque des revendications 1 à 5, dans laquelle les étapes
    a) à d) sont accomplies pour différentes portions réparties le long du lingot semiconducteur. 7
  7. 7. Méthode selon la revendication 6, dans laquelle les différentes portions incluent l’extrémité haute, appelée tête, et l’extrémité basse, appelée queue, du lingot semiconducteur.
  8. 8. Méthode selon l’une des revendications 6 et 7, dans laquelle les différentes portions sont en nombre supérieur ou égal à 5.
FR1762282A 2017-12-15 2017-12-15 Methode de validation de l'histoire thermique d'un lingot semi-conducteur Active FR3075379B1 (fr)

Priority Applications (7)

Application Number Priority Date Filing Date Title
FR1762282A FR3075379B1 (fr) 2017-12-15 2017-12-15 Methode de validation de l'histoire thermique d'un lingot semi-conducteur
CA3085769A CA3085769A1 (fr) 2017-12-15 2018-12-12 Methode de validation de l'histoire thermique d'un lingot semi-conducteur
US16/772,613 US11359305B2 (en) 2017-12-15 2018-12-12 Method for validating the thermal history of a semiconductor ingot
EP18814626.0A EP3724377A1 (fr) 2017-12-15 2018-12-12 Méthode de validation de l'histoire thermique d'un lingot semi-conducteur
PCT/EP2018/084450 WO2019115574A1 (fr) 2017-12-15 2018-12-12 Méthode de validation de l'histoire thermique d'un lingot semi-conducteur
KR1020207020082A KR20200098612A (ko) 2017-12-15 2018-12-12 반도체 잉곳의 열 이력을 검증하는 방법
CN201880080782.0A CN111479957B (zh) 2017-12-15 2018-12-12 用于验证半导体晶锭热史的方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1762282A FR3075379B1 (fr) 2017-12-15 2017-12-15 Methode de validation de l'histoire thermique d'un lingot semi-conducteur
FR1762282 2017-12-15

Publications (2)

Publication Number Publication Date
FR3075379A1 true FR3075379A1 (fr) 2019-06-21
FR3075379B1 FR3075379B1 (fr) 2019-11-22

Family

ID=62749020

Family Applications (1)

Application Number Title Priority Date Filing Date
FR1762282A Active FR3075379B1 (fr) 2017-12-15 2017-12-15 Methode de validation de l'histoire thermique d'un lingot semi-conducteur

Country Status (7)

Country Link
US (1) US11359305B2 (fr)
EP (1) EP3724377A1 (fr)
KR (1) KR20200098612A (fr)
CN (1) CN111479957B (fr)
CA (1) CA3085769A1 (fr)
FR (1) FR3075379B1 (fr)
WO (1) WO2019115574A1 (fr)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113447159B (zh) * 2021-08-30 2021-11-09 中国空气动力研究与发展中心设备设计与测试技术研究所 基于全局正态分布及局部离群因子的风洞漏冷监测方法
CN114235899B (zh) * 2021-12-16 2023-11-03 安徽光智科技有限公司 一种超高纯锗单晶载流子浓度的检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3045074A1 (fr) * 2015-12-14 2017-06-16 Commissariat Energie Atomique Procede pour ajuster la resistivite d'un lingot semi-conducteur lors de sa fabrication

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1327041C (zh) * 2002-11-12 2007-07-18 Memc电子材料有限公司 用于生长单晶锭的拉晶机和方法
TWI290182B (en) * 2004-01-27 2007-11-21 Sumco Techxiv Corp Method for predicting precipitation behavior of oxygen in silicon single crystal, determining production parameter thereof, and storage medium storing program for predicting precipitation behavior of oxygen in silicon single crystal
FR2964459B1 (fr) 2010-09-02 2012-09-28 Commissariat Energie Atomique Procede de cartographie de la concentration en oxygene
FR2974180B1 (fr) * 2011-04-15 2013-04-26 Commissariat Energie Atomique Procede de determination de la concentration en oxygene interstitiel.
FR3009380B1 (fr) 2013-08-02 2015-07-31 Commissariat Energie Atomique Procede de localisation d'une plaquette dans son lingot
FR3027676B1 (fr) * 2014-10-22 2016-12-09 Commissariat Energie Atomique Procede de caracterisation de la concentration en oxygene interstitiel dans un lingot semi-conducteur

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3045074A1 (fr) * 2015-12-14 2017-06-16 Commissariat Energie Atomique Procede pour ajuster la resistivite d'un lingot semi-conducteur lors de sa fabrication

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
J. VEIRMAN ET AL: "Thermal History Index as a bulk quality indicator for Czochralski solar wafers", SOLAR ENERGY MATERIALS AND SOLAR CELLS., 1 May 2016 (2016-05-01), NL, XP055297084, ISSN: 0927-0248, DOI: 10.1016/j.solmat.2016.05.051 *
JORDI VEIRMAN ET AL: "A Fast and Easily Implemented Method for Interstitial Oxygen Concentration Mapping Through the Activation of Thermal Donors in Silicon", ENERGY PROCEDIA, ELSEVIER, NL, vol. 8, 17 April 2011 (2011-04-17), pages 41 - 46, XP002717502, ISSN: 1876-6102, DOI: 10.1016/J.EGYPRO.2011.06.099 *

Also Published As

Publication number Publication date
WO2019115574A1 (fr) 2019-06-20
US11359305B2 (en) 2022-06-14
EP3724377A1 (fr) 2020-10-21
CA3085769A1 (fr) 2019-06-20
KR20200098612A (ko) 2020-08-20
US20210079555A1 (en) 2021-03-18
CN111479957B (zh) 2022-08-23
CN111479957A (zh) 2020-07-31
FR3075379B1 (fr) 2019-11-22

Similar Documents

Publication Publication Date Title
EP2612136B1 (fr) Procede de cartographie de la concentration en oxygene
EP2697632B1 (fr) Procede de determination de la concentration en oxygene interstitiel
EP3184673B1 (fr) Procédé d&#39;étalonnage d&#39;un four de recuit utilisé pour former des donneurs thermiques
WO2014064347A1 (fr) Procede de formation d&#39;un lingot en silicium dopé de resistivite uniforme.
FR3075379B1 (fr) Methode de validation de l&#39;histoire thermique d&#39;un lingot semi-conducteur
WO2015015065A1 (fr) Procede de localisation d&#39;une plaquette dans son lingot
EP3510640B1 (fr) Procédé de tri de plaquettes en silicium en fonction de leur durée de vie volumique
EP2834625B1 (fr) Determination de la concentration en oxygene interstitiel dans un echantillon semi-conducteur
EP3210005B1 (fr) Procédé de caractérisation de la concentration en oxygène interstitiel dans un lingot semi-conducteur
WO2017102832A1 (fr) Procédé pour ajuster la résistivité d&#39;un lingot semi-conducteur lors de sa fabrication
EP3210006B1 (fr) Procédé de caractérisation de la concentration en oxygène interstitiel dans un lingot semi-conducteur
EP3218533A1 (fr) Procede de fabrication d&#39;un lingot de silicium monocristallin de type n
FR3059820B1 (fr) Procede de mesure de temperature
FR3059821B1 (fr) Procede de mesure de temperature
FR2978600A1 (fr) Procede et dispositif de fabrication de couche de materiau semi-conducteur
EP4016596A1 (fr) Procédé de détermination de la durée de vie volumique des porteurs de charge d&#39;un substrat et dispositif associé
JP2011114119A (ja) エピタキシャルウェーハ及びその製造方法

Legal Events

Date Code Title Description
PLFP Fee payment

Year of fee payment: 2

PLSC Publication of the preliminary search report

Effective date: 20190621

PLFP Fee payment

Year of fee payment: 3

PLFP Fee payment

Year of fee payment: 4

PLFP Fee payment

Year of fee payment: 5

PLFP Fee payment

Year of fee payment: 6

PLFP Fee payment

Year of fee payment: 7