Brevet FR2849505A1 — publié le 20040702

PROCEDE DE CARACTERISATION OPTIQUE DE MATERIAUX SANS UTILISATION DE MODELE PHYSIQUE DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne un procédé de caractérisation optique de matériaux. Ce procédé permet de caractériser des couches minces ou épaisses de ces matériaux, qui sont formées sur des substrats. Les grandeurs physiques, que ce procédé permet de déterminer, sont: - l'épaisseur d'une couche d'un matériau, - l'indice de réfraction de ce matériau, et - le coefficient d'absorption de ce matériau. La caractérisation optique des matériaux est utile pour l'analyse chimique de ces matériaux (notamment étude des bandes d'absorption, des propriétés de densification et des propriétés d'oxydation), dans les domaines de la micro- électronique, des capteurs, de la biologie, de la médecine), ou bien pour l'analyse des épaisseurs de dépôts de ces matériaux. Pour des exemples d'application, on se reportera au document [1] qui, comme les autres documents cités par la suite, est mentionné à la fin de la présente description. La caractérisation des propriétés optiques d'un matériau est également utile lorsque le matériau est par la suite structuré ( pour y former par exemple des gravures ou des rugosités) et que les propriétés de diffraction optique de la structure obtenue doivent être calculées (voir le document [2]). Indiquons dès à présent que l'invention est particulièrement utile lorsque la loi physique suivie par l'indice de réfraction complexe du matériau que l'on veut caractériser n'est, a priori, pas connue. ETAT DE LA TECHNIQUE ANTERIEURE On rappelle que les mesures optiques peuvent être de diverses natures : Il peut s'agir de mesures réflectométriques. Dans ce cas, le coefficient de réflexion en intensité d'une structure est mesuré sur un spectre (c'est-à-dire un intervalle) de longueurs d' ondes [#m, #M]. L'angle d'incidence de la lumière d'éclairement peut être non nul. Le coefficient de réflexion peut être mesuré pour plusieurs angles d'incidence #. Nous noterons R (#, À, p) le spectre réflectométrique, où p est la polarisation du faisceau incident et # la longueur d'onde de ce dernier. Généralement, l'angle # est nul et la polarisation p indéterminée. Dans le cas où # n'est pas nul, il faut connaître cette polarisation p. En général, cette dernière est de type (S) où (P). Il peut s'agir aussi de mesures ellipsométriques. Les grandeurs mesurées sont alors les parties réelle et imaginaire du rapport du coefficient de réflexion en polarisation (P) au coefficient de réflexion en polarisation (S). On note généralement p=pexp(j) ce rapport complexe (avec j2=-1). Les grandeurs généralement traitées sont 1 pl, que l'on note tan (\|/) , et cos(), ou des combinaisons des deux. Par exemple, les variables issues d'un ellipsomètre à modulation de phase sont les suivantes : Is=sin (2yl} sin (A) et Ic=cos (2yl) . Un ellipsomètre classique fournit, quant à lui, les variables suivantes : a=(tan2\|/-l) /(tan 2V+ 1) et (3=cosÀ (1-a2) 1/2 . Par souci de généralité, nous noterons Si et S2 les grandeurs traitées. Les spectres Si, i # [1,2], sont mesurés sur une plage de longueurs d' ondes [An,, zon]. L'angle d'incidence peut être quelconque. Plusieurs spectres peuvent être mesurés à différents angles d'incidence afin d'obtenir un spectre plus riche. Nous noterons s (6, 7) =S1 (6, ,) , SZ (A,,) } le spectre ellipsométrique . De façon complémentaire, des mesures goniométriques (coefficient de réflexion en fonction de l'angle d'incidence) peuvent être ajoutées aux mesures servant à la caractérisation, afin de déterminer l'épaisseur des diverses couches, pour une ou plusieurs longueurs d'ondes. Ces mesures ne sont pas suffisantes en elles-mêmes puisque l'on veut déterminer l'indice de réfraction complexe sur une gamme spectrale allant de 7v,I, à 7r,, . Afin de simplifier l'exposé, nous noterons # un ensemble de spectres réflectométrique ou/et ellipsométrique(s) (et éventuellement goniométriques pour quelques longueurs d'ondes). Sans perdre aucune généralité, nous n'exposerons, dans la présente description, le mode d'utilisation des procédés de l'art antérieur et de la présente invention que dans le cas d'une seule couche mince d'un matériau, formée sur un substrat connu. L'épaisseur de cette couche est notée # et l'indice de réfraction complexe du matériau à la longueur d' onde , est noté n* (,) . On rappelle à ce propos que la partie réelle (respectivement imaginaire) de cet indice de réfraction complexe est notée n(#) (respectivement k(#)) et appelée "indice de réfraction" (respectivement "coefficient d'extinction"). En outre on note Er ((1), T" une fonction d'erreur (par exemple l'écart quadratique moyen) entre deux spectres (1) et 'Y2 . Par exemple, on pourra prendre, lorsque l'on dispose de spectres ellipsométriques sur plusieurs angles #i, ie {l..n}, et d'un spectre réflectométrique : Er{' (1), ' (2))= avec (1) ) -S1 m (8i, À) , S2 (1) (Si, À) , R (1) (À)} Y(2> (2) (;) = S 1 (21 ),, ;) , s22 (0., ,) , R (2) (2,) et i # {1...n} Des facteurs de pondération peuvent être apportés à l'intégrale de façon à ce que la fonction d'erreur puisse tenir compte de variations sur la précision de mesure des spectres. La caractérisation optique de couches de matériaux s'articule généralement autour de deux applications : La première application est le contrôle dimensionnel du dépôt de couches minces que l'on utilise en microélectronique. Généralement on connaît bien le matériau déposé, c'est-à-dire que l'on connaît bien l'indice de réfraction complexe de ce matériau aux longueurs d'ondes de la lumière utilisée pour la caractérisation. Les lois suivies par l'indice de réfraction complexe sont soit tabulées, soit approchées par des lois physiques connues telles que, par exemple, le modèle de Cauchy, le modèle de Sellmeier (voir le document [3]), les lois de Forouhi (voir le document [4]), et les lois lois d'oscillateurs harmoniques (voir le document [5]). Ces lois sont définies par un nombre fini de paramètres. Par exemple, une loi de type Cauchy sans absorption, à deux paramètres, est définie de la manière suivante : Re[*(A)] =n(À)=ao+a12 Imn*(1.)]=k (7,)=0 Lorsque l'on est sûr de la valeur des coefficients a~j~(iE{0,l}) mais que l'on ne connaît pas l'épaisseur, un algorithme de recherche est utilisé afin de trouver l'épaisseur qui minimise l'erreur entre la mesure Y et la réponse théorique ' compte tenu de l'indice modélisé. L'algorithme de recherche peut être, par exemple, la méthode du Simplex, la recherche Tabou, la méthode de Levendt-Marquart ou la méthode du recuit simulé (voir le chapitre 10 du document [6]). Lorsque l'indice de réfraction est approximatif, les coefficients ai sont intégrés dans la procédure d'ajustement de # et #. La recherche des coefficients ai constitue un procédé de caractérisation de l'indice de réfraction. Cependant, lorsque la loi suivie par cet indice de réfraction est inconnue (il arrive que le matériau soit inconnu ou qu'il ne soit pas bien décrit par une loi physique connue), ce procédé reste approximatif et l'épaisseur risque d'être fausse. La deuxième application est la caractérisation des matériaux. Le procédé utilisé reste le même, mis à part que le matériau n'est pas bien connu. C'est précisément la fonction d'indice de réfraction complexe la plus proche de la réalité qui est visée. Le type de loi peut être choisi par analogie avec d'autres matériaux. Cependant, la loi suivie par l'indice de réfraction complexe peut être compliquée, ce qui est par exemple le cas d'une loi d'oscillateurs harmoniques : Dans l'expression ci-dessus, j2=-1 et l'indice de réfraction et le coefficient d'extinction sont exprimés non pas en fonction de # mais de E, avec E= 1240/À (# en nm) . Dans ce cas, les coefficients des oscillateurs sont difficiles à trouver si l'on n' a pas leur ordre de grandeur. La recherche est difficilement automatisable, les algorithmes de recherche pouvant donner des réponses erronées et le temps perdu peut être considérable. Il existe une alternative à la recherche des coefficients : la méthode point-à-point (PAP). Cette méthode PAP propose de ne pas choisir de loi physique et de rechercher l'indice de réfraction complexe du matériau pour chaque longueur d'onde #i, où ie [l...n], avec À1=Àm et /.n=%1,M . Pour chaque #i, un algorithme de recherche tente de trouver l'épaisseur, l'indice n(#i) et le coefficient d'extinction k(#i) qui minimisent l'erreur entre la mesure Y (#i) et la réponse théorique V(Âi,n(Âi),k(Â.),e) . Un tel procédé pose un problème parce que les divers points (,i, E, n (,i) , k (Ài) ) ne sont pas forcément physiquement compatibles entre eux : par exemple, l'épaisseur trouvée peut varier en fonction de la longueur d'onde et la loi suivie par l'indice de réfraction complexe, plus simplement appelée loi d'indice, peut présenter des discontinuités. Ce procédé est généralement valable seulement lorsque l'épaisseur est très bien connue et que les mesures sont de très bonne qualité. EXPOSÉ DE L'INVENTION La présente invention a pour but de remédier aux inconvénients précédents. Le procédé objet de l'invention permet de caractériser un matériau sans utiliser un modèle physique, c'est-à-dire sans utiliser une loi physique suivie par l' indice de réfraction complexe du matériau étudié. Il est donc tout particulièrement utile lorsqu'une telle loi n'est pas connue. Ce procédé constitue une alternative aux procédés de caractérisation connus, mentionnés plus haut. Il peut être appelé "méthode des n#uds" car il utilise des "n#uds" c'est-à-dire des points de coordonnées (Xi n j ) , où n est la valeur prise par l'indice de réfraction complexe à la longueur d'onde #i et i prend un nombre limité de valeurs (entières). De façon précise, la présente invention a pour objet un procédé de caractérisation optique d'au moins une couche d'un matériau dans un intervalle A de valeurs prises par une fonction a d'une longueur d'onde optique #, lorsque À varie dans un intervalle de longueurs d'ondes, cette couche étant formée sur un substrat, ce procédé étant caractérisé en ce qu'il comprend les étapes suivantes : 1) on effectue un ensemble de mesures de réflectométrie et/ou d'ellipsométrie sur l'intervalle A, cet ensemble de mesures conduisant à un spectre mesuré, noté #, et l'on choisit les méthodes de calcul associées à la nature des mesures et au type de couche à caractériser, 2) on choisit m valeurs initiales [alpha]1 ...[alpha]m de la fonction a, appartenant à cet intervalle A, m étant un nombre entier au moins égal à 1, et l'on définit un intervalle B comme étant l'ensemble des points a de l'intervalle allant du plus petit au plus grand des nombres [alpha]1...[alpha]m, lorsque m est supérieur à 1, et comme étant l'intervalle A lorsque m vaut 1, 3) on choisit m valeurs initiales complexes d'un indice de réfraction complexe n*=n+jk aux m points [alpha]i, i allant de 1 à m, 4) lorsque m est différent de 1, on choisit une loi d'interpolation qui permet de calculer l'indice de réfraction n([alpha]) du matériau sur l'intervalle B, à partir des points ([alpha]i, ni), avec ni=n (ai) , i allant de 1 à m, et lorsque m vaut 1, n([alpha]) est pris égal au nombre ni (ai) sur tout l' intervalle B, 5) on choisit M paramètres variables, M étant inférieur ou égal à 2m+l, 6) on choisit une fonction d'erreur Er (T, # ) qui caractérise l'écart entre un spectre mesuré T et un spectre expérimental F, 7) à l'aide d'une fonction de minimisation de Er (Y, F) à M paramètres, on effectue la série d'étapes suivante : a) à l'aide de la loi d'interpolation des ([alpha]i,ni) sur l' intervalle B, on déduit n(a), a appartenant à B, b) à l' aide de n(a) et de l' épaisseur # de la couche, et de méthodes de calcul de spectres, on calcule un spectre théorique 'Y(n(a),s), c) on compare T et Y à l' aide de Er (T, 'F ) et, si Er(#, #) est suffisamment petite, c'est-à-dire inférieure à une valeur prédéfinie e, ou est minimale, on passe à l'étape e), sinon on passe à l'étape d), d) on fait varier les M paramètres variables de façon à tendre vers le minimum de Er(#, #) et l'on retourne à l'étape a), e) si Er (T, 'I ) est inférieure à e, on obtient donc un ensemble de M paramètres variables, pour lequel Er (, (n (a, M) , E)) est minimum et l'indice de réfraction est alors pris égal à celui qui a été obtenu en dernier, et si Er(#, #) est supérieure ou égale à e on passe à l'étape 8), 8) on augmente le nombre m de valeurs initiales de la fonction a et l'on retourne à l'étape 2) . Il est donc par exemple possible d'effectuer la caractérisation optique : - sur un intervalle de longueurs d'ondes #, en l' occurrence l' intervalle [#min, Àmax] , - ou sur un intervalle de longueurs d'ondes inverses 1/#, en l'occurrence sur un intervalle [ (1/.}min, (1/7)max], où (1/7}min est égal à 1/(,max) et (1/7)max à 1/(,min), - ou sur un intervalle d'énergies E (avec E = hv = hc/# où h est la constante de Planck, c la vitesse de la lumière dans le vide et v la fréquence correspondant à #), en l'occurrence sur un intervalle [Emin, Emax], où Emin est égal à hc/ (Àmax) et Emax à hc/ (?min) , - ou, plus généralement, sur un intervalle [amin, amax] de valeurs prises par une fonction a de la variable #. Il convient en outre de noter que l'invention est utilisable pour caractériser un spectre ou une partie de spectre. Chaque loi d'interpolation peut être choisie parmi les lois d'interpolation linéaires, les lois d'interpolation cubiques, les lois d'interpolation polynômiales et les lois d'interpolation par exemple de type fonction spline . Selon un mode de mise en #uvre préféré du procédé objet de l'invention, le spectre est échantillonné uniformément en a (#), c'est-à-dire que les valeurs initiales de la fonction a (voir l'étape 2) mentionnée plus haut) sont réparties de manière uniforme dans l' intervalle A, la répartition des n#uds étant ainsi homogène. Comme on l'a vu, [alpha](#) peut être choisi parmi 7, 1/Â et hc/Â ou toute autre fonction de 7, où h est la constante de Planck et c est la vitesse de la lumière dans le vide. De préférence, à l'étape 6) mentionnée plus haut, on mesure l'erreur sur un intervalle d'intérêt C qui est inclus dans l'intervalle B ou égal à cet intervalle B. Les M paramètres variables peuvent être les parties réelles des indices de réfraction aux points ai, i allant de 1 à m, ou les parties imaginaires de ces indices de réfraction, ou ces M paramètres variables peuvent être constitués par l'épaisseur du matériau dont on cherche l'indice de réfraction. BRÈVE DESCRIPTION DES DESSINS La présente invention sera mieux comprise à la lecture de la description d'exemples de réalisation donnés ci-après, à titre purement indicatif et nullement limitatif, en faisant référence aux dessins annexés, sur lesquels : - la figure 1 est une vue schématique de dispositifs permettant de caractériser une couche conformément à l'invention, - la figure 2 montre les variations de l'indice de réfraction en fonction de la longueur d'onde, pour un matériau suivant une loi de Cauchy (courbe I) et pour un matériau caractérisé conformément à l'invention (courbe II), - la figure 3A (respectivement 3B) montre les variations de l'indice de réfraction (respectivement du coefficient d'extinction) en fonction de la longueur d'onde, pour un matériau suivant une loi à deux oscillateurs harmoniques (courbe I) et pour un matériau caractérisé conformément à l'invention, et - la figure 4 illustre schématiquement les paramètres utilisés dans une généralisation d'exemples de l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS L'invention propose une alternative aux procédés classiques, mentionnés plus haut. Elle permet d'allier la cohérence d'un modèle de couche (correspondant à une loi d'indice continue et à des épaisseurs physiques constantes), à la généralité concernant la loi d'indice à trouver (comme dans la méthode PAP). En outre, la résolution n'est limitée que par la résolution du spectre mesuré. Dans le procédé objet de l'invention, le spectre d'indice n* (À) est caractérisé par : - un nombre réduit de "n#uds", qui sont des points de coordonnées (,i,ni,kl) ou (rai, n*i) , avec nl = n (Ai), ki = k (Ai) et n*i=ni+jki, où j2=-1, et - une loi d'interpolation entre les n#uds, qui peut être, par exemple, linéaire, cubique, de type spline ou polynômiale (de degré quelconque). Cette loi d'interpolation permet de calculer, à partir des n#uds, les indices de réfraction et les coefficients d'extinction pour les longueurs d'ondes situées entre les n#uds. Par exemple, lorsque l'indice de réfraction est caractérisé par un ensemble de valeurs pour des longueurs d' ondes #1 ... #m, on peut utiliser une interpolation linéaire entre deux longueurs d'ondes #i et #i+1 pour calculer l'indice n à la longueur d'onde # (voir le document [6] chapitre 3) : On peut faire de même pour le coefficient d'extinction. Lorsque le nombre de n#uds le permet, des formules d'interpolation plus complexes, faisant intervenir les n#uds voisins, peuvent être utilisées (voir le document [6] chapitre 3). Un modèle de couche est donc caractérisé par une épaisseur E et une famille de n#uds. On donne ci-après un exemple du procédé objet de l'invention. Dans cet exemple, les mesures # sont constituées d'une mesure réflectométrique R(#) et d' une mesure ellipsométrique Si,2 (#,#) où # est l'angle d'incidence du faisceau lumineux que l'on envoie sur la couche à étudier lors de la mesure d'ellipsométrie. Cette couche est une couche mince de sorte que l'épaisseur de cette couche est aussi une variable du problème. En outre, on suppose qu'une seule couche est inconnue, cette couche étant formée sur un substrat connu. Expliquons d'abord brièvement cet exemple, qui utilise un algorithme (algorithme de la méthode des noeuds conformément à l'invention). A partir d'informations supposées sur l'épaisseur # de la couche étudiée et sur l'indice de réfraction n (,) et le coefficient d'extinction k(Â) du matériau de cette couche, on construit des n#uds de départ (en faible nombre) et une épaisseur de départ #. On dispose ainsi de m n#uds et, par interpolation, on peut connaître n (,) et ka) en dehors des valeurs des longueurs d'ondes associées aux n#uds. A partir de l'épaisseur de départ E et de ces valeurs de départ n(.) et k (,) , on détermine le spectre théorique # en utilisant des calculs ellipsométriques et réflectrométriques. Par ailleurs, au moyen de dispositifs d'ellipsométrie et de réflectrométrie et d'un spectromètre, on obtient Si,2 (8, À) et R (À) et l'on en déduit les mesures notées (pour des conditions de mesure # et #). On compare ensuite T et F en utilisant une fonction d'erreur Er et l'on optimise la valeur de l'indice de réfraction et la valeur du coefficient d'extinction aux différents n#uds, ainsi que la valeur de l'épaisseur, en cherchant à minimiser Er (#,#). Lorsque ces valeurs sont optimisées et si la précision sur le spectre n(Â), le spectre k(X) et l'épaisseur # n'est pas suffisante, on ajoute de nouveaux n#uds, on fait varier l'épaisseur #, et l'on recommence la détermination de #, la comparaison de # et # et l'optimisation que l'on a mentionnées plus haut, etc. On arrête la boucle ainsi définie lorsque la précision sur chacun des spectres n (À) et k(#) et sur l'épaisseur E est jugée suffisante (ajustement satisfaisant de T et # ) . Les spectres n(À), k (,) et l' épaisseur e sont ainsi caractérisés. Dans un exemple donné à titre purement indicatif et nullement limitatif, on trouve une épaisseur E égale à 212,3nm pour la couche. La figure 1 montre de façon schématique la couche étudiée 2, formée sur un substrat 4. On voit le dispositif d'ellipsométrie 5,6, le dispositif de réflectrométrie 8 et le spectromètre 10. On voit en outre des moyens électroniques de traitement 12, comprenant un ordinateur et permettant de caractériser n (,) , k (À) et s en fonction des informations fournies par le spectromètre 10 et conformément au procédé de l'invention. Ces moyens 12 sont munis de moyens d'affichage 14 qui permettent, en particulier, d'afficher la courbe des variations de n en fonction de # et la courbe des variations de k en fonction de #. Revenons de façon plus détaillée sur l'exemple donné. Phase 1 Le procédé de cet exemple comporte d'abord une étape d'initialisation. L'algorithme débute avec un nombre réduit de n#uds, plus précisément au moins un n#ud. On peut donc débuter par un seul n#ud, en imposant un indice de réfraction et un coefficient d'extinction qui restent constants lorsque la longueur d'onde varie. On choisit des positions de n#uds de façon à pouvoir, à partir de cette famille de n#uds, déduire tout le spectre par interpolation. Le modèle de couche est donc à 3 paramètres ou plus, puisque l'épaisseur est aussi une variable à déterminer. La table d'indice sur tout le spectre est donc déduite des n#uds par interpolation. C'est le cas lorsque l'épaisseur des couches est de l'ordre de grandeur de la longueur d'onde (couches minces). Mais lorsque l'épaisseur de la couche est supérieure à la longueur de cohérence de la source de lumière, l'épaisseur n'intervient quasiment plus dans le calcul de la réponse de la couche, et n'est donc plus une variable du problème. Par exemple, dans le cas de disques optiques (CDROM), sur lesquels une couche très importante est déposée (son épaisseur est de l'ordre de 1 millimètre), le coefficient de réflexion d'une telle couche n'est plus fonction de l'épaisseur de la couche mais seulement de l'indice de réfraction de celle-ci, la longueur de cohérence du faisceau lumineux incident étant inférieure à l'épaisseur de la couche. Dans ce cas précis, la longueur de cohérence du faisceau incident est déterminée par la rugosité des couches. On choisit, par exemple, de placer les deux premiers n#uds aux extrémités À min et # max du spectre. Les valeurs de l'indice complexe en ces extrémités sont choisies en fonction du type du matériau étudié. Par exemple, sur un spectre ellipsométrique entre 300nm et 800nm d'une couche mince de résine photosensible, on prend n (300nm) =n (800nm) =1, 5 et k(300nm)=k(800nm)=0. Lorsque le spectre n'est caractérisé que par deux n#uds, l'indice entre les extrêmes est déterminé par interpolation linéaire. Dans le cas considéré, on a donc n((,) =1,5 et k (,) =0 pour À appartenant à [300nm, 800nm]. A partir de trois n#uds, on choisit plutôt une interpolation cubique afin d'obtenir des formes de loi d'indice plus douces que celles obtenues par interpolation linéaire. L'épaisseur de départ est, quant à elle, choisie aussi proche que possible de l'épaisseur réelle. Phase 2 On procède ensuite à une détermination optimale des valeurs de l'indice de réfraction et du coefficient d'extinction sur les n#uds et de la valeur de l'épaisseur. Pour ce faire, les spectres (À) sont calculés à l'aide du modèle de couche utilisé, résultant du choix des n#uds, de la loi d'interpolation et de l'épaisseur de la couche. Le modèle physique utilisé pour le calcul de F est bien sûr fonction de la méthode de mesure utilisée, c'est-à-dire notamment de l'angle d'incidence de la lumière, du spectre utilisé et du modèle de couches minces ou de couches épaisses le cas échéant (voir par exemple les modèles de couches empilées dans le document [3]). Les spectres T étant constitués par un ensemble de mesures de natures diverses (par exemple des mesures ellipsométriques et réflectométriques), on utilise pour les mesures de réflectométrie (respectivement d'ellipsométrie) une méthode de calcul de réflectométrie (respectivement d'ellipsométrie). Les mesures réflectométriques et ellipsométriques sont combinées par l'intermédiaire d'une fonction d'erreur Er(#,#) qui est par exemple du genre de celle qui est définie par l'équation (1). Grâce à une fonction de recherche, on minimise l'écart entre T (#) et #(#), en faisant varier la valeur de l'indice de réfraction et la valeur du coefficient d'extinction à la position de chacun des n#uds ainsi que l'épaisseur de couche (si cette épaisseur est un facteur important dans le calcul de #). Lorsque l'écart est minimum c'est-à-dire lorsque Er (', I' ) est minimum, cela signifie que les mesures réflectométriques et ellipsométriques coïncident au mieux (pour un nombre de n#uds donné). A ce stade, on obtient, pour un nombre connu de n#uds et une position spectrale connue pour chacun de ces des n#uds, le modèle de couche (indice de réfraction, coefficient d'extinction et épaisseur) qui correspond au mieux à la couche réelle. La validité du modèle trouvé est d'autant plus assurée que le nombre de mesures est grand. Pour avoir un grand nombre de mesures, on peut par exemple utiliser plusieurs angles d'incidence de la lumière #i, 1 # i #, pour les spectres ellipsométriques, faire une mesure réflectométrique et faire des mesures goniométriques complémentaires. Phase 3 Ensuite, on augmente le nombre de n#uds. On ajoute un nombre fini de n#uds. Dans un premier mode de réalisation, les n#uds rajoutés sont positionnés de façon à représenter au mieux les spectres n(,) et k (,) . A titre d'exemple, on place ces n#uds supplémentaires aux endroits où l'écart entre # et # est maximum ou aux endroits en lesquels les n#uds sont les plus espacés. Et l'on retourne à la phase 2 tant que la précision sur chacun des spectres n(#) et k (À) et sur l'épaisseur E n'est pas suffisante, c'est-à-dire n'est pas égale à une précision prédéfinie. Dans un deuxième mode de réalisation, on intercalle de nouveaux n#uds entre deux n#uds de l'ensemble de n#uds préalablement choisis, ces nouveaux noeuds se répartissant de façon homogène sur le spectre. Il convient toutefois de noter ce qui suit. Lorsque l'on augmente le nombre de n#uds, il n'est pas nécessaire de conserver la position des anciens n#uds. Par exemple, admettons que l'on échantillonne uniformément un spectre de 400nm à 800nm avec 3 n#uds. Ces n#uds se situent donc respectivement à 400nm, 600nm et 800nm. Lorsque l'on passe à 6 n#uds, les trois n#uds additionnels ne peuvent être placés de façon à ce que le spectre soit uniformément échantillonné si l'on conserve la position des anciens n#uds. La position des 6 n#uds peut être définie, si l'on souhaite un échantillonnage uniforme par les valeurs 400,480, 560,640, 720, 800nm. L'ancien n#ud central à 600nm disparaît donc. Pour calculer la valeur de l'indice à ces positions à partir des anciens n#uds, on procède par interpolation. Dans ce qui suit, nous donnons deux exemples courants d'application de l'invention. Ces exemples mettent en jeu deux types de matériaux qui suivent des lois différentes. A partir de mesures ellipsométriques et réflectométriques, nous proposons de retrouver les lois physiques suivies par ces matériaux. Nous procédons de la façon suivante : Un matériau fictif est créé, ce matériau suivant une loi théorique connue (une loi de Cauchy ou une loi d'oscillateurs harmoniques), avec des paramètres que nous fixons arbitrairement. Les variations de l'indice de réfraction complexe en fonction de la longueur d'onde sont ainsi parfaitement connues. Nous imposons en outre au matériau une épaisseur de 200,00nm, sur un substrat de silicium, ce dernier étant également très bien connu. Des mesures fictives (mesures ellipsométriques, mesures réflectométriques) sont calculées, puis bruitées de façon à introduire un défaut d'appareillage. Tout se passe comme si nous avions des mesures réelles, effectuées sur le matériau. Mais, contrairement à la réalité, nous connaissons parfaitement l'indice de réfraction complexe puisque nous l'avons fixé, de même que nous avons fixé l'épaisseur de la couche de matériau. Nous testons ici la méthode "en aveugle", c'est-à-dire que nous partons d'une épaisseur fausse (220nm) et d'indices de réfraction complexes faux, puisqu'ils sont sensés être inconnus. Nous appliquons le procédé objet de l'invention puis nous comparons l'indice de réfraction complexe trouvé à l'indice de réfraction complexe théorique. Nous retrouvons bien les mêmes lois, très précisément, ainsi que la même épaisseur de couche. Comme premier exemple, prenons un matériau dont l'indice de réfraction complexe suit une loi de Cauchy telle que : Cette loi d'indice est typique des résines photosensibles (dans la plage spectrale allant de 300nm à 800nm) . Afin de retrouver à l'aide de la méthode des n#uds (c'est-à-dire du procédé objet de l'invention) la loi d'indice mentionnée ci-dessus, nous effectuons deux mesures, à savoir une mesure ellipsométrique à un angle de 70 et une mesure réflectométrique. Les conditions de traitement sont les suivantes : - le spectre traité est compris entre 300nm et 800nm, - au départ, les n#uds sont aux positions (300nm, 1,6) et (800nm, 1,6), c'est-à-dire que l'indice est considéré comme variant linéairement entre 300nm et 800nm, et que sa valeur est constante (égale à 1,6), - le nombre de n#uds est augmenté itérativement selon la séquence 2#4#6, - lors de la procédure d'augmentation des n#uds, la position en longueur d' onde de chaque n#ud est calculée de façon que l' échantillonage en 1/# soit uniforme (À : longueur d'onde), la densité de points étant donc d'augmentée vers les faibles longueurs d'ondes, - la loi d'interpolation est une loi cubique, lorsque le nombre de n#uds est supérieur à 2, sinon elle est linéaire, et - l'algorithme de minimisation utilisé est un algorithme de type Simplex. La figure 2 permet de comparer l'indice de réfraction correspondant au matériau fictif qui suit parfaitement une loi de Cauchy (courbe I) à l'indice de réfraction que nous trouvons par la méthode des n#uds (courbe II), à l'aide de 6 n#uds (représentés par des cercles sur la figure 2). Nous avons utilisé un # composé d'une mesure ellipsométrique {S1 (7) , SZ (7) à 70 et d'une mesure réflectométrique R(#). Nous aboutissons à une épaisseur de 199,8384nm. Comme deuxième exemple, prenons un matériau dont l'indice de réfraction complexe suit une loi à deux oscillateurs harmoniques, telle que : avec j2 = -1 et E = 1240/À (, en nm) Al=0, 25x1, 52 AZ=0, 25x1, 52 Ei=1240/400 E2=1240/300 Gi=0,3 G2=0,3 Dans ce deuxième exemple, la méthode des n#uds est appliquée à un ensemble de mesures ellipsométriques effectuées entre 250nm et 800nm, à 75 , 70 , 60 et 45 , avec en plus une mesure réflectométrique. L'épaisseur réelle du matériau étant de 200nm, on trouve une épaisseur de 200,25 nm avec la méthode des n#uds. L'ajustement sur la loi d'indice considérée dans ce deuxième exemple est très bon, comme le montrent les figures 3A et 3B. Ces figures 3A et 3B illustrent respectivement les reconstructions des courbes n(,) et k(#) du matériau à l'aide de la méthode des n#uds. La reconstruction est effectuée à l'aide de quatre spectres ellipsométriques et d'un spectre réflectométrique. Les pics d'absorption réels sont très bien représentés par la courbe obtenue par interpolation cubique entre les n#uds (représentés par des cercles sur les figures 3A et 3B). Sur la figure 3A, la courbe I (respectivement II) correspond à un indice de réfraction n qui suit parfaitement la loi choisie (respectivement à un indice de réfraction n trouvé par la méthode des n#uds). Sur la figure 3B, la courbe I (respectivement II) correspond à un coefficient d'extinction k qui suit parfaitement la loi choisie (respectivement à un coefficient d'extinction k trouvé par la méthode des n#uds) . On vient de décrire des exemples de l'invention. On remarquera d'une manière plus générale que, dans cette dernière, on considère un ensemble X de valeurs, avec X= {nI, n2, ..., ni, ...nm, kir k2, rkir...,km,F-1, où ni est la valeur de l'indice de réfraction (réel) au noeud correspondant à #i, ie{l...m}, m étant le nombre de noeuds ki est la valeur du coefficient d'absorption au n#ud correspondant à 7i, iE {l...m, # est l'épaisseur de la couche étudiée. Dans ce cas, l'opération de minimisation d'erreur Er (', I' ) revient à trouver l' ensemble ou vecteur X tel que Er soit minimum. Lorsque l'on n'impose pas de contrainte particulière, la minimisation est une minimisation à 2xm+l paramètres. On peut bien sûr ajouter des contraintes afin de diminuer le nombre de variables. En particulier, si l'on sait que le matériau est non absorbant, on impose ki=0 pour tout i de {l...m} et X devient : X={ni,n2...,ni,...nm, . Si, par une mesure complémentaire (par exemple une mesure de goniométrie ou une mesure directe non optique), on connaît l'épaisseur de la couche considérée avec une précision suffisante, l'épaisseur # n'est plus une variable et l'on a : X={ni,n2,...,ni,...nm,ki, r k2, ..., ki, ..., kmo . Bien entendu, les deux options précédentes peuvent être combinées. On explique dans ce qui suit un mode de mise en #uvre de l'invention sous une forme plus générale que les exemples donnés précédemment. Soit a (#) une fonction de la longueur d'onde # de la lumière utilisée pour les mesures. On peut par exemple choisir : [alpha](#)=# (voir les figures 3A et 3B où le spectre est uniformément échantillonné en #) oc(,)=1/, (voir la figure 2 où le spectre est uniformément échantillonné en 1/#) a (7) =hc/, où h est la constante de Planck et c la vitesse de la lumière dans le vide, [alpha](#) étant alors homogène à une énergie. Soit A l'intervalle du spectre de mesure, B l'intervalle spectral décrit par les noeuds et C l'intervalle d'intérêt. L'intervalle C est inclus dans l'intervalle B ou égal à cet intervalle B. De même, l'intervalle B est inclus dans l'intervalle A ou égal à cet intervalle A. On précise que chacun des intervalles A, B et C est du type [0,, ocM] où #m est inférieur à ocM et il existe deux longueurs d' ondes Àk et Ai telles que #m=a(Àk) et aM=a(À1). A titre purement indicatif et nullement limitatif, la figure 4 montre un exemple des intervalles A, B et C et des courbes des variations de T et de T en fonction de oc (.} , I' étant en fait une fonction de n*(a(À)). Les cercles N représentent les n#uds. n* est l'indice complexe, qui est exprimé ici en fonction de a (À) et dont les parties réelle et imaginaire sont respectivement notées n(a(À)) et k ([alpha](#)). On voit également, à titre purement indicatif et nullement limitatif, un exemple de la courbe des variations de n (respectivement k) en fonction de [alpha](#), passant par des points de coordonnées (ni, ai) (respectivement (ki, ai) , où [alpha]i=[alpha](#i), ni=n (ai) , k@=k([alpha]i), lim (m entier naturel non nul). Sur la figure 4, on notera la correspondance entre les n#uds N, les points (ni, ai) et les points (ki, ai) . Dans le mode de mise en #uvre considéré, on utilise un algorithme comprenant les étapes suivantes : 1. on fait les mesures # sur l'intervalle A et l'on choisit les méthodes de calcul associées aux mesures (calculs ellipsométriques ou réflectométriques) ; 2. on choisit m nombres ai (constituant m valeurs initiales de la fonction a), i appartenant à {l,...,m}, avec m#1, et {[alpha]i} C A (les ai correspondant aux n#uds ) ; lorsque m>l, B est défini comme l'ensemble des points a tels que min (ai) 3. on choisit m valeurs initiales d'indice complexe n* aux m points ai, i appartenant à {l,...,m} ; 4. si m#1 on choisit une loi d'interpolation qui permet de calculer l'indice de réfraction n([alpha]) sur l' intervalle B à partir des points (ai, ni) , i appartenant à {l,...,m}; si m=l, alors n (oc) =nl (al) sur tout l' intervalle B ; 5. on choisit M paramètres variables avec M#2m+1 ; ces paramètres peuvent être par exemple : - les parties réelles des indices de réfraction aux points [alpha]i,i appartenant à {l,...,m}, ou - les parties imaginaires des indices de réfraction en ces mêmes points, ou - l'épaisseur du matériau dont on cherche l'indice de réfraction ; 6. on choisit une fonction d'erreur Er(#, #) qui caractérise l'écart entre un spectre expérimental et un spectre mesuré ; en général, l'erreur est mesurée sur l'intervalle C ; 7. à l'aide d'une fonction de minimisation de Er(#,#) à M paramètres, on effectue la série d'étapes suivantes : a) à l'aide de la fonction d'interpolation des ([alpha]i, ni) sur B, on déduit n([alpha]) avec a appartenant à B ; b) à l' aide de n([alpha]) de l'épaisseur E et des méthodes de calcul de spectre, on calcule le spectre théorique T (n(a),e) ; c) on compare et # à l' aide de Er=Er(#,#) ; si Er est suffisamment petite (c'est-àdire si Er est inférieure à une valeur prédéfinie e), ou bien si Er est minimale, on passe à l'étape e), sinon on passe à l'étape d) ; d) on fait varier les M paramètres de façon à tendre vers le minimum de Er et l'on revient à l'étape a) ; e) si Er est inférieure à e on obtient donc un ensemble de M paramètres tel que ER(T,T (n(a,M) , ) ) est minimum et le calcul de l' indice de réfraction est achevé : cet indice est pris égal à celui qui a été obtenu en dernier ; et si Er est supérieure ou égale à e, on passe à l'étape 8) ; 8. on augmente m et l'on revient à l'étape 2) . Il convient de noter que la présente invention peut utiliser non seulement des échantillonnages en longueur d'onde (#) mais encore des échantillonnages en fréquence (c/# ou, plus simplement 1/#), en énergie (hc/#) et, plus généralement, en paramètre fonction de la longueur d'onde. Il convient en outre de noter qu'une étape essentielle de l'invention (étape 8) mentionnée plus haut) n'est pas limitée à l'ajout d'un ensemble de n#uds aux n#uds déjà existants : elle englobe le cas plus général où le nombre de n#uds augmente . Cela veut dire, dans un mode de réalisation particulier de l'invention, qu'après minimisation avec 3 n#uds, si l'on veut passer à 6 n#uds au total, la position des 3 anciens n#uds est effacée de façon, par exemple, à avoir une densité de n#uds constante sur le spectre. En pratique, c'est la meilleure option. L'information sur l'ancienne position des anciens n#uds n' est pas perdue car la valeur des indices aux anciens n#uds est utilisée pour calculer la valeurs des 6 nouveaux n#uds (en fait 3 nouveaux n#uds plus trois anciens n#uds). Ainsi, selon un mode de mise en #uvre particulier du procédé objet de l'invention, on peut augmenter le nombre de valeurs initiales de la fonction a en ajoutant une ou plusieurs valeurs aux valeurs initiales existantes ; mais, selon un mode de mise en #uvre préféré, on peut augmenter le nombre de valeurs initiales de la fonction a en remplaçant les valeurs initiales existantes par de nouvelles valeurs initiales dont le nombre est supérieur au nombre des valeurs initiales existantes. La présente invention n'est pas limitée à la caractérisation de couches minces. Elle s'applique aussi à la caractérisation de couches épaisses. En outre, la présente invention n'est pas limitée à la caractérisation d'une seule couche, formée sur un substrat. Elle s'applique aussi à la caractérisation de deux, ou plus de deux, couches formées sur un substrat. Les documents cités dans la présente description sont les suivants : [1] R. M.A. Azzam and N. M. Bashara, Ellipsometry and Polarized Light, North-Holland Physics Publishing, 1997, chapitre 6. [2] B. K. Minhas, S.A. Coulombe, S. Sohail, H. Naqvi and J.R. McNeil, Ellipsometric scatterometry for the metrology of sub-0.1-m-linewidth structures, Applied Optics, 37(22) : 5112-5115,1998. [3] M. Born and E. Wolf, Principle of Optics, Cambridge University Press edition. [4] A. R. Forouhi and I. Bloomer, Optical dispersion relations for amorphous semiconductors and amorphous dielectrics, Physical Review B, 34 (10) 7018-7026, November 1986. [5] F. L. Terry, Jr., A modified harmonic oscillator approximation scheme for the dielectric constants of AlxGa1-xAs, Journal of Applied Physics, 70(1), 1991, pages 409-417. [6] W. H. Press, S. A. Teukolsky, W. T. Vetterling and B. P. Flannery, Numerical Recipes in C, Cambridge University Press, 1992, chapitres 3 et 10. REVENDICATIONS 1. Procédé de caractérisation optique d'au moins une couche d'un matériau dans un intervalle A de valeurs prises par une fonction a d'une longueur d'onde optique #, lorsque À varie dans un intervalle de longueurs d'ondes , cette couche étant formée sur un substrat, ce procédé étant caractérisé en ce qu'il comprend les étapes suivantes : 1) on effectue un ensemble de mesures de réflectométrie et/ou d'ellipsométrie sur l'intervalle A, cet ensemble de mesures conduisant à un spectre mesuré, noté #, et l'on choisit les méthodes de calcul associées à la nature des mesures et au type de couche à caractériser, 2) on choisit m valeurs initiales [alpha]1 ...On, de la fonction a, appartenant à cet intervalle A, m étant un nombre entier au moins égal à 1, et l'on définit un intervalle B comme étant l'ensemble des points a de l'intervalle allant du plus petit au plus grand des nombres [alpha]1...[alpha]m, lorsque m est supérieur à 1, et comme étant l'intervalle A lorsque m vaut 1, 3) on choisit m valeurs initiales complexes d'un indice de réfraction complexe n*=n+jk aux m points [alpha]i, i allant de 1 à m, 4) lorsque m est différent de 1, on choisit une loi d'interpolation qui permet de calculer l'indice de réfraction n(a) du matériau sur l'intervalle B, à partir des points (ai, ni), avec ni=n(ai), i allant de 1 à m, et lorsque m vaut 1, n([alpha]) est pris égal au nombre ni (ai) sur tout l'intervalle B, 5) on choisit M paramètres variables, M étant inférieur ou égal à 2m+l, 6) on choisit une fonction d'erreur Er(#, #) qui caractérise l'écart entre un spectre mesuré # et un spectre expérimental #, 7) à l'aide d'une fonction de minimisation de Er(#, #) à M paramètres, on effectue la série d'étapes suivante : a) à l' aide de la loi d'interpolation des ([alpha]i,ni) sur l'intervalle B, on déduit n(a), a appartenant à B, b) à l' aide de n(a) et de l' épaisseur # de la couche, et de méthodes de calcul de spectres, on calcule un spectre théorique (n (oc) ,e , c) on comparez et Y à l' aide de Er (Y, P ) et, si Er(#, #) est suffisamment petite, c'est-à-dire inférieure à une valeur prédéfinie e, ou est minimale, on passe à l'étape e), sinon on passe à l'étape d), d) on fait varier les M paramètres variables de façon à tendre vers le minimum de Er(#, #) et l'on retourne à l'étape a), e) si Er(#, #) est inférieure à e, on obtient donc un ensemble de M paramètres variables, pour lequel Er(#, # (n([alpha],M),#)) est minimum et l' indice de réfraction est alors pris égal à celui qui a été obtenu en dernier, et si Er (#, #) est supérieure ou égale à e on passe à l'étape 8), 8) on augmente le nombre m de valeurs initiales de la fonction a et l'on retourne à l'étape 2). 2. Procédé selon la revendication 1, dans lequel on augmente le nombre de valeurs initiales de la fonction a en ajoutant une ou plusieurs valeurs aux valeurs initiales existantes. 3. Procédé selon la revendication 1, dans lequel on augmente le nombre de valeurs initiales de la fonction a en remplaçant les valeurs initiales existantes par de nouvelles valeurs initiales dont le nombre est supérieur au nombre des valeurs initiales existantes. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel chaque loi d'interpolation est choisie parmi les lois d'interpolation linéaires, les lois d'interpolation cubiques, les lois d'interpolation polynômiales et les lois d'interpolation par exemple de type fonction spline. 5. Procédé selon l'une quelconque des revendications 1 à 4 , dans lequel les valeurs initiales de la fonction a sont réparties de manière uniforme dans l'intervalle A, la répartition des n#uds étant ainsi homogène. 6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel [alpha](#) est choisi parmi À, 1/# et hc/À, où h est la constante de Planck et c est la vitesse de la lumière dans le vide. 7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel on mesure l'erreur, à l'étape 6), sur un intervalle d'intérêt C qui est inclus dans l'intervalle B ou égal à cet intervalle B. 8. Procédé selon l'une quelconque des revendications 1 à 7, dans lequel les M paramètres variables sont les parties réelles des indices de réfraction aux points ai, i allant de 1 à m. 9. Procédé selon l'une quelconque des revendications 1 à 7, dans lequel les M paramètres variables sont les parties imaginaires des indices de réfraction aux points ai, i allant de 1 à m. 10. Procédé selon l'une quelconque des revendications 1 à 7, dans lequel les M paramètres variables sont constitués par l'épaisseur du matériau dont on cherche l'indice de réfraction.