FR 3060752 A1 20180622 FR 1662509 A 20161215 Procédé de mise en oeuvre d'une technique de caractérisation CD-SEM DOMAINE La présente invention concerne le domaine de la métrologie dimensionnelle et a pour objet un procédé de mise en oeuvre d'une technique de caractérisation par microscopie à balayage d'électrons pour la détermination d'une dimension critique, dite technique CD-SEM. Les techniques CD-SEM sont aujourd'hui largement utilisées pour la caractérisation de structures micrométriques ou nanométriques comme, par exemple, les motifs utilisés dans les circuits imprimés en microélectronique. ETAT DE L'ART La caractérisation de structures de taille micrométrique ou nanométrique est aujourd'hui d'importance primordiale dans l'industrie des semiconducteurs, notamment en microélectronique, et plus généralement dans le domaine des nanotechnologies. Si l'on considère le cas de la microélectronique, la mesure précise et fiable des dimensions caractéristiques des nanostructures se rend nécessaire tout au long des différentes étapes du processus de production, à partir de la phase de recherche et développement jusqu'à la phase de production. De plus, la miniaturisation des composantes microélectroniques demande des outils de caractérisation toujours plus performants. L'industrie des semiconducteurs définit la taille de ses produits en termes d'une taille caractéristique appelée dimension critique ou « critical dimension » (CD) en anglais. On entend par dimension critique la ou les dimensions critiques d'un motif représentant une ou plusieurs longueurs caractéristiques du motif qui s'avèrent critiques soit pour la maitrise du procédé de fabrication soit pour garantir les performances électriques du dispositif électronique final constitué desdits motifs. La mesure fiable de cette dimension critique est aujourd'hui l'un des enjeux principaux pour la microscopie électronique à balayage. Lors de l'acquisition d'une image par microscopie électronique à balayage, un faisceau d'électrons, dit faisceau d'électrons primaires, balaye la surface d'un échantillon d'intérêt. Chaque point de la surface est identifiée par un couple de cordonnées spatiales par exemple les cordonnées cartésiennes (x, y). Les collisions 1 3060752 entre électrons primaires et les atomes de l'échantillon produisent plusieurs phénomènes physiques qui peuvent être très complexes. Parmi ces phénomènes il y a l'émission d'électrons secondaires et l'émission d'électrons rétrodiffusés. Ces électrons sont captés par des détecteurs sélectifs à proximité de l'échantillon. En 5 balayant point par point l'échantillon on parvient à la formation d'une image électronique reflétant les différentes formes et différents matériaux de l'échantillon. Dans ce qui suit on se réfèrera aux différentes formes de l'échantillon et leurs différents matériaux comme à la structure de l'échantillon. L'image SEM ainsi obtenue est une image en niveau de gris et les mesures 10 géométriques des objets sont déterminées par une analyse des contours des objets obtenus après le choix d'un ou plusieurs seuils. On peut par exemple reconstruire un graphique comme illustré sur la figure 1 représentant un profil d'intensité d'électrons secondaires PI comportant en ordonnées un pourcentage d'électrons secondaires recueillis et en abscisse une dimension en nm représentative de la position de 15 balayage. Mais un des problèmes majeurs est l'interprétation des images par des algorithmes à seuillage : c'est en effet le choix de ce seuil qui va déterminer la valeur de la taille caractéristique, ou dimension critique, d'intérêt. Or le choix du seuil est souvent arbitraire et le choix optimal est différent pour chaque objet, notamment en fonction de la hauteur des motifs, des matériaux en jeu, des dimensions des objets, 20 de la densité de motifs, etc. D'autres types d'algorithmes sont basés non pas directement sur le niveau de gris mais sur l'analyse de la forme des pics d'intensité et le choix d'un seuil relatif à l'amplitude des pics. Là encore, le choix du seuil relatif à l'analyse du pic est empirique. Actuellement il est connu que les techniques CD-SEM donnent des mesures peu fiables pour des objets de taille inférieure à 100 nm 25 (écart à la réalité du nanomètre à plusieurs nanomètres). Un seuil fixe empirique appliqué à tout type de motif est d'autant moins satisfaisant qu'il n'y a aucun lien physique direct entre le pourcentage d'électrons secondaires recueillis et la hauteur réelle du motif à laquelle on cherche à déterminer la dimension critique. En d'autres termes, si une mesure à 80% d'électrons secondaires est réalisée, ce n'est pas pour 30 autant que l'on mesure une dimension critique à 80% de la hauteur du motif. Une analyse des images SEM avec un algorithme de seuillage peut donc très facilement fausser les mesures des dimensions critiques. 2 3060752 Pour remédier à ce problème, deux approches ont été envisagées. La première est basée sur des simulations physiques de type Monte Carlo. La mise en oeuvre de cette méthode requiert la simulation des processus physiques d'interaction entre les électrons primaires du faisceau incident sur l'échantillon et l'échantillon lui- 5 même. La réalisation de ces simulations très complexes demande une grande puissance de calcul et des temps d'exécution très longs. Ces caractéristiques rendent les simulations physiques de type Monte Carlo généralement peu accessibles et très peu adaptées au traitement d'un grand nombre d'images CD-SEM. La seconde solution consiste à réaliser une modélisation mathématique 10 paramétrique des images ; une telle solution est notamment décrite dans l'article « CD characterization of nanostructures in SEM metrology » (C. G. Frase, E. Buhr, et K. Dirscherl, Meas. Sci. Technol., vol. 18, n° 2, p. 510, février. 2007). Un exemple de telles fonctions est illustré en figure 2b pour un motif ayant un profil selon la direction x trapézoïdale, représenté en figure 2a. Ces fonctions peuvent être choisies selon un 15 modèle phénoménologique simple ou bien par l'analyse de modélisations Monte- Carlo préliminaires décrivant la formation des images SEM (plus complexes et précises mais nécessitant des temps de calcul longs). Selon l'exemple de la figure 2, le motif est décomposé suivant des sections selon la direction x (ici six sections), chaque section correspondant à une fonction de forme particulière. La fonction M(x) 20 peut ainsi se décomposer comme suit (en six fonctions selon chacune des sections) : Chacune des six fonctions paramétriques reproduit la hauteur du profil de l'échantillon en fonction de la cordonnée x dans une certaine gamme de valeurs de 25 x. Le modèle présenté ci-dessus comporte donc une pluralité de paramètres : x1, x2, 3 3060752 X3, x4, x5, So, S1, S2, S3, tO, t1, t2, t3, t4, t5, A, B, C, D. L'existence de ce jeu de paramètres justifie l'appellation de fonction mathématique paramétrique. Souvent, pour des applications en métrologie dimensionnelle, les paramètres les plus intéressants sont ceux relatifs aux dimensions latérales du motif à savoir : x1, x2, x4 et 5 X5. Une autre méthode de modélisation paramétrique des images CD-SEM est décrite dans l'article « Analytical Linescan Model for SEM Metrology » (C. A. Mack et B.D. Bunday, Metrology, Inspection and Process Control for Microlitography XXIX, Proc. SPIE Vol. 9424 2015). Selon ce modèle la fonction M(x) peut se décomposer 10 comme suit : Dans ce modèle les paramètres sont So, S1, S2, af, , ab, ae,av, Gf, Gb, Ge, Gy, d1, d2,X2. De plus, pour tenir compte de la taille du faisceau d'électrons primaires, qui 15 n'est pas nulle, il convient d'associer au profil décrit par M(x) un produit de convolution incluant dans le modèle les paramètres décrivant la taille et la forme du faisceau. Le plus souvent le profil du faisceau d'électrons primaires est caractérisé à l'aide d'un profil gaussien (« point spread fonction » ou PSF en anglais), le paramètre apsF étant proportionnel à la largeur à mi-hauteur de la 20 gaussienne. Le modèle paramétrique ainsi modifié est calculé à l'aide de la formule suivante : (x-x')2 M(x) = f e 2upsF2 * M(x')dx' Dans la suite l'expression « jeu de paramètres du modèle théorique » se 25 réfèrera à l'ensemble des paramètres utilisés dans la formulation mathématique du 4 3060752 modèle théorique. Parmi ces paramètres, certains sont utilisés dans la description géométrique de la structure de l'échantillon, d'autres décrivent la réponse instrumentale. Une fonction paramétrique telle que la fonction M(x) sera appelée dans la suite profil modèle ou profil paramétrique. 5 Un autre type de modélisation consiste à modéliser la probabilité d'émission d'un électron secondaire en fonction du matériau et la propagation de ces électrons vers les détecteurs en fonction de l'angle des détecteurs (voir par exemple les modèles de type AMAT, « SEM Simulation for 2D and 3D Inspection Metrology and Defect Review » de Shimon Levi, Ishai Schwarzband, Sergey Khristo, Yan Ivanchenko, Ofer 10 Adan, Advances in Patterning Materials and Processes XXXI, 2014). Des modèles plus simples ou plus complexes peuvent être élaborés en fonction du type d'échantillon, de la qualité des images, ou bien des informations a priori sur les matériaux ou les profils. Une fois la modélisation terminée, il faut résoudre un problème inverse, à 15 savoir l'ajustement des données expérimentales au modèle en faisant varier les paramètres. Une des limitations majeures de ces techniques est dans leur incapacité à prendre en compte et à reproduire certains artefacts observés sur les images CD-SEM. C'est par exemple le cas des images acquises en mode de balayage dit TV (ou « raster scan » en anglais). Ces images montrent des marques sombres dans les 20 régions de la surface de l'échantillon comprises entre deux ou plusieurs structures voisines, comme les plots carrés ou rectangulaires de la figure 3a. Sur cette figure des directions x et y sont indiquées, la direction x étant celle de balayage et la direction y étant perpendiculaire à la direction x. Plus précisément ces marques sombres peuvent être classées suivant deux types : 25 - Des marques plutôt linéaires en gris foncé ou en noir dans le prolongement des arêtes horizontales alignées avec la direction de balayage et indiquées par le repère 1 de la figure 3a ; - Des marques surfaciques en gris moins foncé le long des motifs, toujours dans la direction de balayage et indiquées par le repère 2 de la figure 3a. 30 Ces marques sombres sont des artefacts créés pendant le processus de formation de l'image. 5 3060752 On peut remarquer que la méthode des profils paramétriques ne rend compte que d'une modélisation des profils suivant la direction perpendiculaire aux arêtes et n'est donc pas en mesure de reproduire les effets de traînées noires présentes dans le prolongement des arêtes mêmes. Il est intéressant de noter que, comme illustré 5 sur les figures 3b, 3c et 3d, la présence des traînées noires augmente quand plusieurs structures rapprochées sont présentes. Plus précisément, la figure 3b montre une structure carrée isolée. Des traces sombres sont visibles à gauche et à droite de la structure, indiquées par le repère li, avec des traînées noires dans le prolongement des arêtes horizontales. La figure 3c montre une image SEM de trois 10 structures carrées rapprochées et alignées suivant la direction horizontale. Les marques sombres, indiquées par 12, sont plus intenses que sur la figure 3b. La figure 3d montre une image SEM d'une série de plusieurs structures carrées. Les artéfacts de type marques sombres selon la direction de balayage, indiquées par 13, sont cette fois plus intenses que sur les images 3b et 3c. Ces altérations ne peuvent pas être 15 décrites suivant une vision purement locale comme celle des profils paramétriques. PROBLEME TECHNIQUE Dans ce contexte, la présente invention a pour but d'améliorer l'exploitation des images CD-SEM obtenues expérimentalement par une technique CDSEM paramétrique. zo RESUME DE L'INVENTION A cette fin, l'invention propose un procédé de mise en oeuvre d'une technique de caractérisation par microscopie électronique à balayage pour la détermination d'au moins une dimension critique de la structure d'un échantillon dans le domaine de la métrologie dimensionnelle, dite technique CD-SEM, ledit procédé comportant 25 les étapes suivantes : - réalisation d'une image expérimentale représentative de la structure de l'échantillon et issue d'un microscope électronique à balayage ; - à partir d'un premier modèle théorique basé sur des fonctions mathématiques paramétriques, calcul d'un second modèle théorique obtenu par somme 30 algébrique d'un terme correctif, ledit terme correctif étant le produit de convolution entre un noyau de convolution donné et le premier modèle 6 3060752 théorique, ledit second modèle théorique comportant un jeu de paramètres à déterminer ; - détermination du jeu de paramètres présents dans le second modèle théorique au moyen d'un ajustement entre ledit second modèle théorique et 5 ladite image expérimentale. On entend par dimension critique (« critical dimension » en anglais) la ou les dimensions critiques d'un motif représentant une ou plusieurs longueurs caractéristiques du motif qui s'avèrent critiques soit pour la maitrise du procédé de fabrication soit pour garantir les performances électriques du dispositif électronique 113 final constitué desdits motifs. On entend par modèle théorique basé sur des fonctions paramétriques mathématiques une représentation mathématique de la structure de l'échantillon, par exemple par une des fonctions M(x) définies précédemment et décrites en détail dans les articles « CD characterization of nanostructures in SEM metrology » (C. G. 15 Frase, E. Buhr, et K. Dirscherl, Meas. Sci. Technol., vol. 18, n° 2, p. 510, février. 2007) ou « Analytical Linescan Model for SEM Metrology » (C. A. Mack et B.D. Bunday, Metrology, Inspection and Process Control for Microlitography XXIX, Proc. SPIE Vol. 9424 2015) ou AMAT par exemple. On entend par jeu de paramètres présents dans le second modèle théorique 20 l'ensemble des paramètres utilisés dans la formulation mathématique du second modèle théorique. Ce jeu de paramètres inclut à la fois les paramètres présents dans le premier modèle théorique paramétrique et les paramètres présents dans le terme correctif ajouté lors du calcul du second modèle théorique. Dans ce jeu de paramètres on aura certains paramètres décrivant la géométrie de la structure de 25 l'échantillon, d'autres décrivant la réponse instrumentale. On entend par ajustement entre le second modèle théorique et l'image expérimentale un procédé de calcul itératif visant à trouver le jeu de paramètres qui minimise l'écart entre ledit modèle et ladite image expérimentale. Dans la pratique, ce procédé de calcul est par exemple mis en oeuvre par calculateur en choisissant 30 un des algorithmes connus pour minimiser l'écart entre l'image expérimentale et le modèle théorique, suivant par exemple la méthode dite des moindres carrées. 7 3060752 De façon générale, l'invention consiste à obtenir les dimensions critiques de la structure de l'échantillon au moyen d'un ajustement entre une image expérimentale représentative de la structure de l'échantillon et un modèle théorique, ledit modèle théorique étant calculé de façon à prendre en compte les artefacts présents sur les 5 images expérimentales. Le but de l'invention est l'obtention plus efficace et plus précise des dimensions critiques. Le procédé selon l'invention est un procédé de mise en oeuvre d'une technique de caractérisation CD-SEM qui utilise une modélisation des images 10 expérimentales, ladite modélisation étant capable de décrire, par exemple, les altérations de type variations d'intensité dans le prolongement des arêtes qui apparaissent sur les images. Ces altérations d'intensité sont le résultat de phénomènes qui ne peuvent pas être décrits suivant une vision purement locale comme celle des profils paramétriques. Parmi les phénomènes capables de produire 15 de telles altérations on peut par exemple citer les effets d'accumulation de charge (« charging » en anglais). Ces charges électriques en excès produisent des champs électriques qui peuvent dévier les électrons (primaires et secondaires) et donc produire des baisses du flux d'électrons secondaires correspondant aux traînées noires visibles sur les figures 3a à 3d. Une modélisation plus complexe des images 20 CD-SEM est donc requise. Pour prendre en compte ces effets, le procédé selon l'invention modifie les modèles paramétriques en rajoutant par somme algébrique un terme correctif, ledit terme correctif étant un produit de convolution entre le modèle paramétrique même et un noyau de convolution donné. Grâce à la présence du produit de convolution, ce 25 nouveau modèle permet de prendre en compte la densité d'arêtes ainsi que leur positionnement par rapport, par exemple, à la direction de balayage du faisceau d'électrons primaires. De plus, le choix du noyau de convolution permet d'adapter le procédé selon l'invention à une grande variété de motifs et de méthodes de balayage. 30 Le procédé selon l'invention peut également présenter une ou plusieurs des caractéristiques ci-dessous, considérées individuellement ou selon toutes les combinaisons techniquement possibles : 8 3060752 - le second modèle théorique est calculé selon la formule suivante : +G° Img(x, y) = G(x, y) - C ifF (x', y') * [G(x' , y') - Go]clx'dy' -oe dans laquelle : (x,y) sont les cordonnées spatiales associées à chaque point de 5 l'image suivant les deux directions orthogonales x et y ; G(x,y) est le premier modèle théorique basé sur des fonctions mathématiques paramétriques ; Img(x, y) est le second modèle théorique obtenu par l'application du terme correctif; 10 F (x, y) est le noyau de convolution donné ; Go est la moyenne (« background » en anglais) des intensités présentes sur l'image expérimentale ; C est une constante d'ajustement ; - le faisceau d'électrons primaires balaye la surface de l'échantillon suivant une 15 ligne droite désignée axe x ou axe horizontale selon la méthode de balayage dite TV ou « raster scan » en anglais, la direction de balayage étant la direction horizontale ou axe x ; - le noyau de convolution peut revêtir une forme mono-dimensionnelle symétrique, par l'intermédiaire par exemple d'une gaussienne portée selon 20 l'axe x, ladite gaussienne ayant une largeur à mi-hauteur proportionnelle à sigma selon la formule suivante : +oe (x-x')2 Img(x, y) = G (x, y) - C fe 2,12 * [G (x', y) - Go]dx' -oe - le noyau de convolution peut revêtir une forme bidimensionnelle symétrique, 25 par l'intermédiaire par exemple de deux gaussiennes portées respectivement par les axes x' et y' respectivement parallèle et normale à la direction de balayage du faisceau d'électrons primaires selon la formule suivante : 9 3060752 +co (x-x')2 (Y-37Y Img(x,y) = G(x,y) - C if e 2ex2 e 2°-Y2 * [G(x',y') - Go]clx'dy' - le noyau de convolution peut revêtir une forme monodimensionnelle asymétrique suivant la direction de balayage du faisceau d'électrons primaires. 5 - le noyau de convolution peut revêtir une forme bidimensionnelle asymétrique suivant la direction de balayage du faisceau d'électrons primaires. - le procédé est mis en oeuvre pour calibrer la technique de caractérisation CD-SEM, ladite calibration comportant les étapes suivantes : - réalisation d'une image expérimentale représentative de la structure d'un échantillon de référence dont les dimensions géométriques sont connues, ladite image issue d'un microscope électronique à balayage ; - à partir d'un premier modèle théorique basé sur des fonctions mathématiques paramétriques, calcul d'un second modèle théorique obtenu par somme algébrique d'un terme correctif, ledit terme correctif étant le produit de convolution entre un noyau de convolution donné et le premier modèle théorique, ledit second modèle théorique comportant un jeu de paramètres, ledit jeu de paramètres comprenant à la fois les paramètres connus qui décrivent la structure géométrique de l'échantillon de référence et les paramètres à déterminer qui décrivent la réponse instrumentale ; - détermination des paramètres présents dans le second modèle théorique et décrivant la réponse instrumentale au moyen d'un ajustement entre ledit second modèle théorique et ladite image expérimentale représentative de la structure de l'échantillon de référence. Le procédé de mise en oeuvre d'une technique de caractérisation CD-SEM comporte avantageusement les étapes suivantes : - réalisation d'une image expérimentale représentative de la structure de l'échantillon et issue d'un microscope électronique à balayage 10 15 20 25 10 3060752 - à partir d'un premier modèle théorique basé sur des fonctions mathématiques paramétriques, calcul d'un second modèle théorique obtenu par somme algébrique d'un terme correctif, ledit terme correctif étant le produit de convolution entre un noyau de convolution donné et le premier modèle 5 théorique, ledit second modèle théorique comportant un jeu de paramètres, ledit jeu de paramètres comprenant à la fois les paramètres à déterminer qui décrivent la structure géométrique de l'échantillon d'intérêt et les paramètres déterminés par la calibration décrite ci-dessus qui décrivent la réponse instrumentale ; 10 - détermination des paramètres présents dans le second modèle théorique et décrivant la structure de l'échantillon d'intérêt au moyen d'un ajustement entre ledit second modèle théorique et ladite image expérimentale. LISTE DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront clairement de la 15 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 illustre schématiquement un profil d'intensité d'électrons secondaires en fonction du profil d'un motif obtenu via une instrumentation de type CD-SEM ; 20 - les figures 2a et 2b illustrent un exemple de modélisation mathématique paramétrique d'une image CD-SEM ; - les figures 3a, 3b, 3c et 3d illustrent des exemples d'artefacts de type marques sombres observées sur des images CD-SEM expérimentales ; - la figure 4 représente les différentes étapes du procédé selon l'invention ; 25 - la figure 5a représente une image CD-SEM ; - la figure 5b représente un modèle théorique paramétrique de l'image 5a ; - la figure 5c représente un modèle théorique de l'image 5a obtenu par application du procédé selon l'invention ; - la figure 6a illustre la méthode de balayage dit TV (« raster scan » en 30 anglais) ; - la figure 6b illustre la méthode de balayage en dent de scie ; - la figure 6c illustre la méthode de balayage bidirectionnelle ; 11 3060752 - la figure 6d illustre la méthode de balayage adaptée aux objets mesurés ; - la figure 7a représente une image expérimentale CD-SEM - la figure 7b montre un agrandissement d'une portion de la figure 7a illustrant un effet local particulier en correspondance d'un coin dans la structure de 5 l'échantillon ; - La figure 8 illustre un exemple de profil symétrique du noyau de convolution ; - La figure 9 illustre un exemple de profil asymétrique du noyau de convolution. DESCRIPTION DETAILLEE 10 La figure 4 illustre schématiquement les différentes étapes du procédé 100 selon l'invention. Le procédé 100 selon l'invention vise la mise en oeuvre d'une technique de caractérisation par microscopie électronique à balayage pour la détermination 15 d'au moins une dimension critique de la structure d'un échantillon dans le domaine de la métrologie dimensionnelle. Une application peut être par exemple la mesure d'une dimension critique d'un motif formant un circuit imprimé en microélectronique. La forme du motif peut être quelconque. Le matériau du motif peut être également quelconque. Ce motif peut être par exemple un motif isolé ou appartenant à un 20 réseau de motifs répétés périodiquement. Il peut s'agir d'un motif obtenu après n'importe quelle étape (lithographie, gravure,...) d'un procédé de fabrication. Une étape 101 du procédé 100 selon l'invention consiste à réaliser une image représentative de la structure de l'échantillon et issue d'un microscope électronique à balayage. 25 Selon un mode de réalisation préférentiel de l'étape 101 selon l'invention, le faisceau d'électrons primaires balaye l'échantillon suivant la méthode dite mode TV (ou « raster scan » en anglais). Ce mode de balayage est représenté par la figure 6a. Selon ce mode de balayage on peut identifier une direction de balayage rapide selon la direction x indiquée sur la figure 6a et une direction de balayage lente selon la 30 direction y indiquée sur la figure 6a. Selon cette méthode de balayage, l'acquisition des données a lieu seulement en correspondance des traits horizontaux, indiqués (1) sur la figure 6a. Dans ce qui suit on se référera à la direction de balayage rapide, x sur la figure 6a, comme à la direction de balayage du faisceau d'électrons primaires. 12 3060752 Selon un autre mode de réalisation de l'étape 101 selon l'invention, le faisceau d'électrons primaires balaye l'échantillon suivant la méthode dite en dent de scie. Ce mode de balayage est représenté par la figure 6b. Sur cette figure, le trait plein représente la phase d'acquisition des données (donc de formation de l'image), le trait 5 en pointillé désignant la phase de simple déplacement du faisceau d'électrons primaires sans acquisition des données. Les flèches en gris représentent le sens de balayage du faisceau d'électrons primaires. Selon un autre mode de réalisation de l'étape 101 selon l'invention, le faisceau d'électrons primaires balaye l'échantillon suivant la méthode dite bidirectionnelle. Ce 10 mode de balayage est représenté par la figure 6c. Sur cette figure, le trait plein représente la phase d'acquisition des données (donc de formation de l'image). A la différence des figures 6a et 6b, dans ce cas l'acquisition des données est faite en balayant la surface à la fois dans le sens des valeurs croissantes de la cordonnée x et dans le sens des valeurs décroissantes de la cordonnée x. Les flèches en gris 15 représentent le sens de balayage du faisceau d'électrons primaires. Selon un autre mode de réalisation de l'étape 101 selon l'invention, le faisceau d'électrons primaires balaye l'échantillon suivant une méthode adaptée à la forme des objets à caractériser. Ce mode de balayage est représenté par la figure 6d. Selon une étape 102 du procédé 100 selon l'invention, à partir d'un modèle 20 théorique paramétrique de la surface de l'échantillon on calcule un second modèle théorique capable de prendre en compte les altérations d'intensité de type marques sombres, comme celles présentes dans le prolongement des arêtes et ayant souvent la forme de traînées noires. Plus précisément, l'étape 102 du procédé 100 selon l'invention comporte la 25 modification d'un premier modèle théorique appelé G(x,y) pour obtenir un second modèle théorique appelé Img(x,y). Le modèle G(x,y) est un modèle paramétrique de la réponse du microscope à la structure de l'échantillon obtenu par exemple par une des méthodes exposées dans les articles « CD characterization of nanostructures in SEM metrology » (C. G. Frase, E. Buhr, et K. Dirscherl, Meas. Sci. 30 Technol., vol. 18, n° 2, p. 510, février. 2007) ou « Analytical Linescan Model for SEM Metrology » (C. A. Mack et B.D. Bunday, Metrology, Inspection and Process Control for Microlitography XXIX, Proc. SPIE Vol. 9424 2015) ou AMAT ou autre. La fonction 13 3060752 G(x,y) représente donc une image théorique de la structure de l'échantillon sans effets liés au sens de balayage. Par exemple, l'image théorique sans effet de balayage G(x,y) est obtenue par l'application du modèle M(ô(x,y)) où ô(x,y) est la distance signée aux arrêtes de l'objet imagé. A l'instar de la fonction M(x), ce modèle 5 théorique contient un jeu de paramètres à déterminer et représentant à la fois les informations sur la géométrie de l'échantillon et sur la réponse instrumentale sans la prise en compte du sens de balayage. Ensuite le second modèle théorique, Img(x,y), est calculé. Ledit second 10 modèle théorique est obtenu par somme algébrique d'un terme correctif au premier modèle théorique, ledit terme correctif étant un produit de convolution entre un noyau de convolution donné et le premier modèle théorique. Selon un mode de réalisation de l'étape 102 du procédé 100 selon l'invention, le second modèle théorique est calculé selon la formule suivante (1) : +G° Img(x,y)= G(x,y)-C ifF(x',y')*[G(x',y')-Go]clx'dy' -oe 15 Dans laquelle : - G(x,y) est le premier modèle théorique basé sur des fonctions mathématiques paramétriques ; - (x,y) sont les cordonnées spatiales associées à chaque point de l'image suivant les deux directions orthogonales x et y; 20 - Img(x,y) est le second modèle théorique obtenu par l'application du terme correctif ; - F(x,y) est le noyau de convolution donné ; - Go est la moyenne (« background » en anglais) des intensités présentes sur l'image expérimentale ; 25 - C est une constante d'ajustement. Un avantage de ce mode de réalisation provient dans la formule (1) du produit de convolution calculé entre le noyau F(x,y) et la différence G(x',y')- 60. La soustraction de l'intensité moyenne do permet en effet d'annuler le produit de convolution là où la réponse est constante, et donc le terme correctif, quand 14 3060752 G(x' ,y') = 60. Cette condition est vérifiée pour les points de la surface de l'échantillon qui sont éloignés des motifs, donc pour des portions de l'image sur lesquelles les traînées noires ne sont pas présentes et pour lesquelles la correction n'est pas nécessaire. 5 De façon générale, le noyau de convolution F(x,y) est une fonction des cordonnées spatiales (x,y) ou bien des cordonnées spatiales choisies selon un quelconque système de cordonnées. La présente invention peut être implémentée avec une fonction F(x,y) quelconque, à condition que le résultat de la formule (1) présentée plus haut soit Io bien défini. A titre d'exemple, selon les différents modes de réalisation de la présente invention, des noyaux de convolution symétriques, asymétriques ou anisotropes peuvent être utilisés. Selon un mode de réalisation de l'étape 102 du procédé 100 selon l'invention, le faisceau d'électrons balaye la surface de l'échantillon selon la méthode dite 15 balayage TV (« raster scan » en anglais) illustrée par la figure 6a. Selon un mode de réalisation de l'étape 102, le noyau de convolution est une fonction monodimensionnelle symétrique selon une direction de l'espace. La figure 8 montre un exemple de noyau symétrique 300 au sens de l'invention. Selon un mode de réalisation de l'étape 102, le noyau de convolution est une 20 gaussienne portée selon l'axe x, ledit axe x correspondant à la direction de balayage du faisceau d'électrons primaires, ladite gaussienne ayant une largeur à mi-hauteur proportionnelle à sigma, selon la formule suivante : +oe (x-x')2 Img(x, y) = G (x, y) - C fe 2,12 * [G(x',y) - Go]dx' -oe On remarque que dans ce cas le noyau de convolution F(x',y') est une fonction constante selon la direction verticale y', ce qui explique pourquoi la 25 variable y' n'apparaît pas dans le produit de convolution. Un avantage de ce mode de réalisation est d'être particulièrement adapté à la prise en compte des baisses d'intensité dans le prolongement des arêtes dans la direction horizontale dans le cas de balayage en mode TV. 15 3060752 Selon un mode de réalisation de l'étape 102 du procédé 100 selon l'invention le noyau de convolution est une fonction bidimensionnelle symétrique. Selon un mode de réalisation de l'étape 102 du procédé 100 selon l'invention, le noyau de convolution est symétrique dans les directions respectivement 5 parallèle et normale à la direction de balayage du faisceau d'électrons primaires selon la formule suivante : +,,,, (x-x')2 CY-A2 Img(x,y) = G (x, y) - C If e 20-x2 e 2e2 * [G(x',Y) - Go]clx'dy' -Go Un avantage de ce mode de réalisation est de pouvoir choisir un noyau de convolution capable de rendre compte des effets locaux particuliers, par exemple ceux en correspondance des coins des motifs formant la structure de l'échantillon. 10 Les figures 7a et 7b illustrent un exemple de ces effets. La figure 7a est une image CD-SEM d'une structure de forme carrée. La figure 7b représente un agrandissement de la portion encadrée 10 de la figure 7a. Sur cette figure on peut voir un exemple particulier de marque sombre indiqué par la ligne en pointillé 20. Plus particulièrement on peut reconnaître une dépendance en y, qui s'ajoute à la 15 dépendance en x. Ces observations motivent la nécessité d'un noyau de convolution bidimensionnel et éventuellement anisotrope suivant les deux directions orthogonales x et y. On entend par noyau de convolution anisotrope une fonction bidimensionnelle ayant un profil en x différent du profil en y (par exemple une gaussienne bidimensionnelle avec largeur o-x différente de la largeur ay) 20 Un autre avantage de ce mode de réalisation est de pouvoir être adapté à une méthode de balayage différente du mode TV. Selon un mode de réalisation de l'étape 102 du procédé 100 selon l'invention, le noyau de convolution a un profil asymétrique selon au moins une des directions x ou y. Un exemple est illustré par la figure 9 où l'on voit un exemple de noyau 200 25 de convolution asymétrique selon l'axe x. Dans ce cas on a, par exemple, un profil qui est en partie gaussien 201 et en partie décroissant suivant une dépendance en x non gaussienne 202. Un avantage de ce mode de réalisation est de pouvoir adapter le noyau de convolution au mode de balayage employé et au type de motifs présents sur 30 l'échantillon. 16 3060752 Selon une étape 103 du procédé 100 selon l'invention l'ajustement entre la formule (1) et l'image expérimentale est réalisé. Cet ajustement permet de trouver le jeu de paramètres présents en (1) qui minimise les écarts entre le modèle théorique et l'image expérimentale, par exemple en utilisant un des algorithmes 5 connus qui appliquent la méthode des moindres carrés. A l'instar du modèle G(x,y), le modèle Img(x,y) contient plusieurs paramètres, y-compris les paramètres décrivant la géométrie de l'échantillon et utilisés pour dériver la dimension critique d'intérêt. Selon un mode de réalisation de l'étape 103 du procédé selon l'invention on 10 réalise une étape de calibration pour déterminer les paramètres présents dans le modèle et liés à la réponse instrumentale. Ces paramètres font partie du jeu de paramètres déjà présents dans le premier modèle théorique et ils prennent en compte, par exemple, le fait que le faisceau d'électrons primaires a une taille non nulle. Souvent la forme du faisceau d'électrons primaires est décrite comme un 15 profil gaussien. Cette fonction gaussienne prenant en compte les caractéristiques du faisceau d'électrons primaires est appelée « point spread function » en anglais ou PSF. Les paramètres décrivant la réponse instrumentale, par exemple les paramètres de la fonction PSF, peuvent avantageusement être déterminés lors d'une étape préliminaire de calibration. La calibration est obtenue en réalisant un 20 ajustement entre la formule (1) et l'image expérimentale d'un échantillon dont la structure est connue. Cela permet de fixer les paramètres géométriques lors de l'étape de calibration et d'obtenir de façon plus fiable les paramètres décrivant la réponse instrumentale. Les valeurs des paramètres décrivant la réponse instrumentale seront ensuite fixées lors de la mise en oeuvre de la technique CD- 25 SEM pour la caractérisation d'une image expérimentale d'intérêt. Un avantage de la réalisation de l'étape de calibration est de pouvoir déterminer de façon plus précise et fiable les paramètres décrivant la réponse instrumentale. Ensuite ces paramètres décrivant la réponse instrumentale seront fixés lors de la mise en oeuvre de la technique CD-SEM pour la caractérisation 30 d'une image expérimentale d'intérêt, ce qui permet d'atteindre une détermination plus précise et fiable des paramètres décrivant la structure de l'échantillon et donc des dimensions critiques d'intérêt 17 3060752 Les figures 5a à 5c illustrent les effets obtenus grâce au procédé selon l'invention. La figure 5a est une image CD-SEM d'un échantillon que l'on veut caractériser. Dans cet exemple on considère le cas d'un balayage en mode TV (« raster scan » en anglais). Selon ce mode de balayage, illustré à la figure 6a, le 5 faisceau d'électrons primaires balaye rapidement la surface de l'échantillon suivant la direction x et plus lentement suivant la direction y. Dans ce cas on désigne la direction x simplement comme la direction de balayage du faisceau d'électrons primaires. Le repère Ta indique les traînées noires dans le prolongement des arêtes horizontales, suivant la direction de balayage. Ces 10 artéfacts sont très fréquents sur les images CD-SEM obtenues avec cette méthode de balayage. La figure 5b représente le modèle théorique de l'image CD-SEM basé sur des fonctions paramétriques, ce qu'on a désigné G (x, y) plus haut. Le repère Tb indique les marques sombres qui sont prédites par ce type de modèle. On voit sur cette figure comment l'intensité de ces marques sombres est 15 dépendante uniquement de la distance aux arêtes. De plus, aucune traînée noire du type indiqué par les flèches Ta n'est présente. La figure 5c représente le modèle théorique de l'image CD-SEM calculé au moyen du procédé selon l'invention, à l'aide de la formule suivante : +oe (x-x')2 Img (x, y) = G (x, y) - C fe 2,12 * [G (x', y) - G0] dx' -oe Le noyau de convolution choisi dans ce cas est une gaussienne 20 unidimensionnelle portée selon l'axe x. L'image 5c issue du procédé selon l'invention reproduit les caractéristiques de l'image expérimentale 5a de façon plus précise par rapport au modèle théorique issu des fonctions paramétriques (image 5b). Plus précisément, sur l'image 5c, on peut voir des trainées noires, indiquées par les flèches Tc, qui reproduisent les traînées noires observées sur l'image CD- 25 SEM 5a. Ces images montrent comment le procédé selon l'invention reproduit précisément les artéfacts de type marques sombres ou traînées noires dans le prolongement des arêtes suivant la direction de balayage du faisceau d'électrons primaires. La mise en oeuvre du procédé selon l'invention montrée sur les figures 5a, 5b 30 et 5c peut également être accompagnée par une étape de calibration permettant d'obtenir de façon précise les paramètres décrivant la réponse instrumentale du 18 3060752 microscope électronique à balayage. Pour ce faire il faut préalablement acquérir une image CD-SEM expérimentale d'un échantillon de référence dont les dimensions géométriques sont connues. Cette image est l'image de référence pour l'étape de calibration. Un ajustement entre le second modèle théorique et 5 l'image de référence est ensuite réalisé, en fixant les paramètres géométriques, qui sont connus dans le cas de l'échantillon de référence. Cet ajustement permet d'obtenir les valeurs des paramètres décrivant la réponse instrumentale, par exemple les paramètres de la fonction PSF, ou l'amplitude des signaux. Le fait d'avoir fixé les paramètres géométriques lors de cette étape de calibration permet 10 d'obtenir de façon plus fiable et précise les paramètres décrivant la réponse instrumentale. Les valeurs desdits paramètres décrivant la réponse instrumentale seront ensuite fixées lors de la mise en oeuvre du procédé selon l'invention sur un échantillon d'intérêt. 15 19 REVENDICATIONS1. Procédé (100) de mise en oeuvre d'une technique de caractérisation par microscopie électronique à balayage pour la détermination d'au moins une dimension critique de la structure d'un échantillon dans le domaine de la métrologie dimensionnelle, dite technique CD-SEM, ledit procédé comportant les étapes suivantes : - réalisation (101) d'une image expérimentale représentative de la structure de l'échantillon et issue d'un microscope électronique à balayage ; - à partir d'un premier modèle théorique basé sur des fonctions mathématiques paramétriques, calcul (102) d'un second modèle théorique obtenu par somme algébrique d'un terme correctif, ledit terme correctif étant le produit de convolution entre un noyau de convolution donné et le premier modèle théorique, ledit second modèle théorique comportant un jeu de paramètres à déterminer ; - détermination (103) du jeu de paramètres présents dans le second modèle théorique au moyen d'un ajustement entre ledit second modèle théorique et ladite image expérimentale. 2. Procédé (100) selon la revendication précédente caractérisé en ce que le second modèle théorique est calculé selon la formule suivante : +G° Img(x, y) = G(x, y) - C ifF(x', y') * [G(x' , y') - Go]clx'dy' -oe dans laquelle : (x,y) sont les cordonnées spatiales associées à chaque point de l'image suivant les deux directions orthogonales x et y ; G(x,y) est le premier modèle théorique basé sur des fonctions mathématiques paramétriques ; Img(x,y) est le second modèle théorique obtenu par l'application du terme correctif; F(x, y) est le noyau de convolution donné ; 20 3060 752 Go est la moyenne (« background » en anglais) des intensités présentes sur l'image expérimentale ; C est une constante d'ajustement ; 5 3. Procédé (100) selon l'une des revendications précédentes caractérisé en ce que le faisceau d'électrons primaires balaye la surface de l'échantillon suivant la méthode de balayage dite TV ou « raster scan » en anglais, la direction de balayage étant la direction horizontale ou axe x. 10 4. Procédé (100) selon l'une des revendications précédentes caractérisé en ce que le noyau de convolution est une fonction monodimensionnelle symétrique. 5. Procédé (100) selon la revendication 4 caractérisé en ce que le noyau de convolution est une gaussienne portée selon l'axe x, ladite gaussienne ayant 15 une largeur à mi-hauteur proportionnelle à sigma selon la formule suivante : +oe (x-x')2 Img(x, y) = G (x, y) - C fe 2,12 * [G(x' , y) - Go]dx' -oe 6. Procédé (100) selon l'une des revendications 1 à 3 caractérisé en ce que le noyau de convolution est une fonction bidimensionnelle symétrique. 7. Procédé (100) selon la revendication 6 caractérisé en ce que le noyau de 20 convolution est bidimensionnel symétrique dans les directions respectivement parallèle et normale à la direction de balayage du faisceau d'électrons primaires selon la formule suivante : +co (x-x')2 (Y-37Y Img(x, y) = G (x, y) - C if e 2x2 e 2°-Y2 * [G(x',y') - Go]clx'dy' 8. Procédé (100) selon l'une des revendications 1 à 3 caractérisé en ce que le 25 noyau de convolution a un profil monodimensionnel asymétrique suivant la direction de balayage du faisceau d'électrons primaires ou un profil bidimensionnel asymétrique suivant la direction de balayage du faisceau d'électrons primaires. 21 3060752 9. Procédé selon l'une des revendications 1 à 8, ledit procédé étant mis en oeuvre pour une calibration de la technique de caractérisation CD-SEM, ladite calibration comportant les étapes suivantes : 5 Réalisation d'une image expérimentale représentative de la structure d'un échantillon de référence dont les dimensions géométriques sont connues, ladite image issue d'un microscope électronique à balayage ; - à partir d'un premier modèle théorique basé sur des fonctions mathématiques paramétriques, calcul d'un second modèle théorique obtenu par somme 10 algébrique d'un terme correctif, ledit terme correctif étant le produit de convolution entre un noyau de convolution donné et le premier modèle théorique, ledit second modèle théorique comportant un jeu de paramètres, ledit jeu de paramètres comprenant à la fois les paramètres connus qui décrivent la structure géométrique de l'échantillon de référence et les 15 paramètres à déterminer qui décrivent la réponse instrumentale ; - détermination des paramètres présents dans le second modèle théorique et décrivant la réponse instrumentale au moyen d'un ajustement entre ledit second modèle théorique et ladite image expérimentale représentative de la structure de l'échantillon de référence. 20 10. Procédé de mise en oeuvre d'une technique de caractérisation CD-SEM selon l'une des revendications 1 à 8 et la revendication précédente, ledit procédé comportant les étapes suivantes : - réalisation d'une image expérimentale représentative de la structure de 25 l'échantillon et issue d'un microscope électronique à balayage ; - à partir d'un premier modèle théorique basé sur des fonctions mathématiques paramétriques, calcul d'un second modèle théorique obtenu par somme algébrique d'un terme correctif, ledit terme correctif étant le produit de convolution entre un noyau de convolution donné et le premier modèle 30 théorique, ledit second modèle théorique comportant un jeu de paramètres, ledit jeu de paramètres comprenant à la fois les paramètres à déterminer qui décrivent la structure géométrique de l'échantillon d'intérêt et les paramètres déterminés selon la calibration selon la revendication précédente qui décrivent la réponse instrumentale ; 22 3060752 - détermination des paramètres présents dans le second modèle théorique et décrivant la structure de l'échantillon d'intérêt au moyen d'un ajustement entre ledit second modèle théorique et ladite image expérimentale. 23