FR 2731074 A1 19960830 FR 9502263 A 19950227 La présente invention conceme un procédé et un appareil d'ellipsométrie rapide, ainsi qu'un dispositif de contrôle d'élaboration de couches les mettant en oeuvre. L'ellipsométrie est une technique de mesure non destructive permettant la caractérisation optique d'un échantillon disposant d'une surface spéculaire ou quasi spéculaire. L'ellipsométrie peut être mise en oeuvre in situ et permet alors l'étude des mécanismes de croissance des couches minces, de formation des interfaces et le contrôle de procédé d'élaboration de ces couches et interfaces. L'ellipsométrie est, par exemple, utilisée pour l'étude et le contrôle de la fabrication des semisonducteurs. Les mesures ellipsométriques peuvent être réalisées à une longueur d'onde fixe, ou à plusieurs longueurs d'onde (ellipsométrie spectroscopique). Selon le domaine de longueur d'onde de la source, proche ultraviolet, visible, proche infrarouge, infrarouge, etc..., il est possible d'accéder à des propriétés différentes des couches, des matériaux ou d'explorer des matériaux différents. Dans le domaine de l'ultraviolet et du visible, la profondeur de pénétration du rayonnement est souvent faible. Cela constitue des conditions favorables pour l'étude des surfaces et des interfaces, et pour les contrôles en temps réel. Cela ne permet généralement pas d'accéder aux propriétés volumiques des couches et des matériaux qui peuvent, au contraire1 être obtenues par des mesures dans le domaine de l'infrarouge. L'infrarouge est bien adapté aux mesures d'absorption vibrationnelles (liaisons chimiques). Pour la réalisation des mesures ellipsométriques, la surface d'un échantillon est éclairée par un faisceau lumineux et l'état de polarisation d'un faisceau incident i est comparé à celui du faisceau réfléchi r ou transmis. Un vecteur de polarisation E est généralement représenté par ses projections Es et Epî respectivement perpendiculaires et parallèles au plan d'incidence. Les projections Ep et Es sont des amplitudes complexes. Dans le domaine de l'ellipsométrie, on représente généralement le rapport p = (Ep/ES)r/ (Ep/ESi, significatif des modifications de l'état de polarisation produites par la surface étudiée, sous la forme: p = tg.exp(iA) = (EdEs)r/ (Ep/ES)i Les deux angles w et A décrivant le changement de polarisation sont ainsi combinés dans la quantité complexe p. Les angles v et A, donc le nombre p, dépendent à la fois de propriétés de l'échantillon, de l'angle d'incidence d'un faisceau et de la longueur d'onde de mesure. L'expression de g et A, ou de p, en fonction de ces paramètres, est donnée par les équations de Fresnel citées, par exemple, par D. CHARLOT et A. MARUANI dans Appl. Opt. 24, 3368, 1985. Dans un ellipsomètre à modulation de phase, un rayon incident a sa polarisation modulée par une différence de phase générée entre deux axes propres d'un modulateur de phase. Le déphasage 5(0 évolue typiquement avec le temps t selon une loi périodique de pulsation w, 5(0 étant proportionnel au premier ordre à sin(ot). Dans un ellipsomètre à modulation de phase, I'intensité d'un flux lumineux réfléchi par un échantillon permet de déduire, de façon connue, les valeurs de w et A. L'ellipsométrie, et plus particulièrement l'ellipsométrie spectros- copique à modulation de phase (ESMP), est une technique performante pour mesurer en temps réel la croissance de couches sur un substrat. Cette technique présente l'avantage de ne pas perturber des réactions en cours. Elle est par ailleurs très sensible à des paramètres physiques de l'échantillon mesuré, tels qu'une épaisseur d de couche et un indice n de réfraction. D'autre part1 elle permet des mesures rapides. Selon un procédé connu, les angles w et A, ou p, sont déduits des mesures d'intensité. Ces quantités w et A dépendent de paramètres physiques de l'échantillon mesuré, tels que l'indice n et l'épaisseur d de la couche supérieure. Ces demiers peuvent ainsi être calculés ensuite à partir de w et A, par une inversion directe des équations de Fresnel. Cette inversion doit être en général réalisée de manière itérative. L'application à un diagnostic in situ et à un contrôle de croissance de l'ellipsométrie spectroscopique à modulation de phase est, par exemple, décrite dans le document "Highspeed spectral ellipsometry for in situ diagnostics and process control", DUNCAN et al., J. Vac. Sci. Technol. B., 12(4), 1994. Malgré son efficacité, cette méthode présente l'inconvénient de générer dans certaines circonstances des incertitudes dans les mesures de paramètres physiques. Ces incertitudes peuvent apparaître, en particulier, lors de la croissance d'un matériau transparent sur un substrat absorbant. Elle détériore alors sensiblement la précision des mesures. L'objet de la présente invention est de pratiquer l'ellipsométrie à modulation de phase en s'affranchissant des problèmes d'oscillations périodiques, et d'améliorer ainsi la précision des mesures. Un autre but de l'invention consiste à contrôler sans difficulté particulière le dépôt de matériau transparent sur un substrat. C'est aussi un objectif de l'invention de mesurer des paramètres physiques lors du dépôt d'un film sur un substrat à une vitesse élevée, par exemple supérieure à 30As-1, jusqu'à une épaisseur importante, par exemple supérieure à 4000 A, avec une bonne précision. Un but supplémentaire de l'invention est de mesurer des paramètres physiques d'un échantillon, tels que n et d, sans avoir à mesurer les angles T et A, ou le nombre p. A cet effet, I'invention conceme un procédé de mesure ellipsométrique de paramètres physiques représentatifs d'un échantillon. Dans ce procédé: - on polarise linéairement un faisceau lumineux incident, ledit faisceau étant défini par un vecteur polarisation, - on module le faisceau lumineux incident de lumière polarisée de façon à introduire, entre des composantes perpendiculaires du vecteur polarisation, un déphasage a(fl dépendant du temps t selon une variation périodique de pulsation o, 6(t) étant proportionnel au premier ordre à sin (w0, -on éclaire l'échantillon avec le faisceau lumineux incident de lumière polarisée modulée, - on analyse le vecteur polarisation du faisceau lumineux renvoyé par l'échantillon, - on mesure par au moins un photodétecteur un flux du faisceau lumineux, on effectue par une unité de traitement électronique reliée au photodétecteur, des calculs sur des mesures du flux. Le flux détecté a une intensité I(t) de la forme: 1(t) = 1om + 1sm sin 8(t) +1cm cos 8(t), Iom, 1sm. et Icm étant des valeurs mesurées dans l'unité de traitement à partir de l'intensité I(t), et dépendant des paramètres physiques. Selon l'invention: - dans une première étape, on produit des valeurs théoriques IstlIot et ICt/Iot initiales à partir d'estimations initiales des paramètres physiques, - les valeurs théoriques servent à déterminer, dans une deuxième étape, les estimations suivantes des paramètres physiques dont on déduit des valeurs théoriques Istilot et IctlIot suivantes, - la deuxième étape est réitérée jusqu'à une Nième estimation des paramètres physiques, de façon à minimiser l'écart entre les valeurs théoriques et celles mesurées. Les paramètres physiques sont alors évalués à partir des valeurs ISt/Iot et Ictilot au cours de la Nième estimation. Contrairement aux procédés connus, le procédé selon l'invention ne nécessite donc pas la connaissance de v et A. La détermination des paramètres physiques est effectuée directement à partir des valeurs mesurées Iom, Ism( et Icm à partir de l'intensité I(t). Au lieu d'inverser les équations de Fresnel, on ajuste ainsi les équations donnant les valeurs précédentes en fonction des paramètres physiques1 le processus étant dans ce cas aussi itératif. Au lieu de prendre en compte un unique point, on effectue une moyenne ou une intégration sur un ensemble de points pour les ajustements. Typiquement, on traite plusieurs dizaines de points à la fois. Ces points sont obtenus pour une série de mesures à des instants successifs, ou à plusieurs longueurs d'onde. L'obtention de points à des instants successifs est de préférence employée pour des mesures cinétiques, et l'obtention de points à différentes longueurs d'onde, pour des mesures statiques. En cas d'un dépôt de plusieurs couches successives sur un substrat, une réactualisation de la série de points moyennée est effectuée lors d'un passage à une nouvelle couche. Cette technique permet d'éviter l'utilisation des valeurs gênantes de "s et A entraînant des imprécisions des mesures. Ceci advient, en particulier, lorsque T est voisin de 45". Or, lors du dépôt d'un matériau transparent sur un substrat absorbant, un rayon incident subit de multiples interférences dans la couche transparente. Une couche, même relativement épaisse, n'amortit pas ces interférences, si bien qu'elle n'atténue pas l'influence du substrat par rapport à celle de la couche sur les propriétés de réflexion. Il en résulte que "s et A ne convergent pas vers un point d'arrivée correspondant aux caractéristiques optiques de la couche, comme c'est le cas pour le dépôt d'un matériau absorbant. Au contraire, w et A subissent des variations importantes et périodiques au cours de la croissance. Ainsi, les méthodes d'inversion des équations de Fresnel divergentelles pour des épaisseurs d'une couche d'indice n telles que: avec k entier, o étant l'angle d'incidence d'un faisceau incident et X la longueur d'onde de mesure. Les variations importantes de v et A lors du dépôt d'un matériau transparent ne permettent pas d'éviter les zones d'instabilité engendrant les imprécisions mentionnées plus haut. En s'affranchissant de calculs de v et A et en recourant à une moyenne, il est possible de réduire considérablement ces inconvénients. D'autre part, il fait gagner une étape dans l'acquisition de données, car v et A sont déduits de Iom, Ism, et Icm dans les procédés connus. Un test d'arrêt, en soi connu, est nécessaire pour arrêter les itérations. Il consiste typiquement à comparer les valeurs théoriques Istilot et Ict/Iot aux valeurs mesurées Ism/lom et ICm/Iom. Lorsque leur écart devient inférieur à un pourcentage d'erreurs fixé, on arrête les itérations. D'autre part, les estimations initiales des paramètres physiques ne nécessitent pas une grande précision, les itérations les faisant rapidement converger. De façon connue, il est avantageux de choisir comme estimations initiales, les paramètres physiques précédemment calculés. Quelques itérations suffisent pour converger vers les paramètres physiques. Leur évaluation est ainsi très rapide et peut être effectuée avec un micro-ordinateur standard. La détermination de l'épaisseur d et de l'indice de réfraction n d'une couche déposée sur un substrat est typiquement effectuée en environ deux secondes avec une précision de l'ordre de 1 % sur le produit n x d, avec un micro-ordinateur tel que ceux connus sous la dénomination "PC 486". Selon un mode de réalisation préféré du procédé de mesure selon l'invention, des premiers moyens de calcul de l'unité de traitement donnent des composantes de Fourier So, S1 et S2 de l'intensité du flux, respectivement continue, à la pulsation w et à la pulsation 2. Des deuxièmes moyens de calcul permettent de déduire linéairement des composantes So, S1 et S2, les valeurs mesurées Ioml ISmg et 1cm. On trouvera une description de ces étapes par B. DREVILLON dans "Progress in Crystal Growth and Characterization of Materials", vol. 27, pp. 1-87, 1993. Selon un mode d'utilisation préféré de l'invention, l'échantillon comporte au moins une couche transparente éclairée par un faisceau incident. En effet le procédé selon l'invention est alors particulièrement avantageux par rapport aux procédés connus, comme on l'a vu précédemment. Le matériau transparent est, par exemple, de la silice, le substrat pouvant être constitué de silicium. Le procédé selon l'invention est également applicable à des matériaux absorbants. Selon un mode de réalisation préféré de l'invention, l'écart entre valeurs théoriques et valeurs mesurées est minimisé par une méthode de moindres carrés. Ce choix des moindres carrés donne l'écart entre les valeurs théoriques et mesurées. Pour minimiser cet écart, on utilise un procédé classique, tel que la méthode de Levenberg-Marquardt. Selon un mode de réalisation particulièrement avantageux du procédé de mesure, celui-ci est spectroscopique. Il est ainsi possible de faire varier la longueur d'onde de mesure. Typiquement, une source lumineuse à large spectre émet un faisceau incident, celui-ci est réfléchi par l'échantillon, puis une longueur d'onde du faisceau réfléchi est sélectionnée grâce à un monochromateur. Un autre moyen de choisir une longueur d'onde donnée consiste à faire varier la longueur d'onde de la source. Dans des appareils spectroscopiques séquentiels, on sélectionne successivement une longueur d'onde après l'autre. Au lieu de ces dispositifs traditionnels, on peut aussi utiliser des appareils spectroscopiques de type multiplexé. On détecte alors des faisceaux dispersés sur plusieurs photodétecteurs et on effectue un multiplexage électronique de signaux captés. Ces appareils spectroscopiques de type multiplexé permettent ainsi d'effectuer des mesures simultanées avec plusieurs longueurs d'onde. L'invention conceme également un ellipsomètre comportant: - une source lumineuse émettant un faisceau lumineux incident, - un polariseur polarisant linérairement le faisceau lumineux, - un modulateur de phase générant un déphasage 6(t) dépendant au premier ordre du temps t selon une variation périodique de pulsation o, 8(t) étant proportionnel au premier ordre à sin(wt), - un analyseur analysant l'état de polarisation du faisceau lumineux renvoyé par un échantillon éclairé par le faisceau lumineux incident de lumière polarisée, - un photodétecteur mesurant un flux de faisceau lumineux, et - une unité de traitement électronique reliée au photodétecteur, effectuant des calculs sur des mesures du flux. Le flux détecté a une intensité I(t) de la forme: I(f) Iom + Ism sin 6(f)+ Icm cos 6(t), Iom, Ism, et Icm étant des valeurs mesurées dans l'unité de traitement à partir de l'intensité I(t), et dépendant des paramètres physiques. Dans l'ellipsomètre selon l'invention: - l'unité de traitement produit des valeurs théoriques Istilot et IctlIat initiales à partir d'estimations initiales des paramètres physiques, - les valeurs théoriques servent à déterminer les estimations suivantes des paramètres physiques dont on déduit des valeurs théoriques Ist/Iot et 1ct11ot suivantes, - cette opération est réitérée jusqu'à une Nième estimation des paramètres physiques, de façon à minimiser l'écart entre les valeurs théoriques et celles mesurées. Les paramètres physiques sont évalués à partir des valeurs IstlIot et Id/Iot au cours de la Nième estimation. Dans un mode de réalisation préféré de l'ellipsomètre selon l'invention, celui-ci comporte au moins une fibre de l'ensemble formé par une première et une deuxième fibres optiques, la première fibre optique étant placée entre la source et le polariseur et la deuxième fibre optique étant placée entre l'analyseur et un système de détection comprenant le photodétecteur. En plus du photodétecteur, le système de détection est susceptible de contenir un monochromateur pour des mesures ellipsométriques à une longueur d'onde, ou un spectrographe pour des mesures simultanées à plusieurs longueurs d'onde. Dans le cas d'une émission d'un faisceau laser, le système de détection peut ne contenir que le photodétecteur. L'invention conceme également d'un dispositif de contrôle d'élaboration de couches sur un substrat comportant des moyens de contrôle de paramètres physiques représentatifs de l'élaboration de couches. Le substrat et les couches constituant un échantillon, ces moyens de contrôle comprennent: - une source lumineuse émettant un faisceau lumineux incident, - un polariseur polarisant linénairement le faisceau lumineux, - un modulateur de phase générant un déphasage 8(t) dépendant du temps t selon une variation périodique de pulsation e3, 5(t) étant proportionnel au premier ordre à sin(wt), - un analyseur analysant l'état de polarisation du faisceau lumineux renvoyé par un échantillon éclairé par le faisceau lumineux incident de lumière polarisée, - au moins un photodétecteur mesurant un flux de faisceau lumineux, et - une unité de traitement électronique reliée au photodétecteur, effectuant des calculs sur des mesures du flux. Le flux détecté a une intensité I(t) de la forme: 1(t) = Iom + Ism sin 8(t) +1cm cos 8(t), Iom, Ils,, et Icm étant des valeurs mesurées dans l'unité de traitement à partir de l'intensité I(t), et dépendant des paramètres physiques. Dans le dispositif de contrôle selon l'invention: - l'unité de traitement produit des valeurs théoriques IStnot et ktJIot initiales à partir d'estimations initiales des paramètres physiques, - les valeurs théoriques servent à déterminer les estimations suivantes des paramètres physiques dont on déduit des valeurs théoriques Ist/Iot et Ictîlot suivantes, - cette opération est réitérée jusqu'à une Nième estimation des paramètres physiques, de façon à minimiser l'écart entre les valeurs théoriques et celles mesurées. Les paramètres physiques sont évalués à partir des valeurs Istflot et Ictilot au cours de la Nième estimation. L'élaboration de couches peut consister en un dépôt sur un substrat. Ce dispositif est ainsi particulièrement bien adapté au contrôle de dépôt de couches effectué par des techniques telles que MOCVD (Metalorganic Chemical Vapour Deposition), PECVD (Plasma Enhanced Chemical Vapour Deposition) ou MBE (Molecular Beam Epitaxy). L'élaboration de couches peut aussi être réalisée par gravure. Des mesures selon l'invention permettent non seulement d'observer des mécanismes de croissance, mais aussi d'effectuer un procédé de contrôle en temps réel, comprenant des contre-réactions. Il est ainsi possible d'agir en temps réel sur les paramètres d'élaboration des couches. Des mesures précises sont ainsi effectuées sans perturber les mécanismes de croissance. Dans le dispositif de contrôle d'élaboration de couches, les couches ayant chacune une épaisseur, un indice de réfraction et une fonction diélectrique, les paramètres physiques comprennent, de préférence, au moins certains d'entre eux. Dans le dispositif de contrôle selon l'invention, l'élaboration étant un dépôt, le dépôt d'une couche a avantageusement une vitesse supérieure à3oAs1. Pour une telle vitesse, des mesures obtenues avec le dispositif suivant l'invention sont particulièrement précises comparées à des mesures effectuées avec un dispositif connu, avec une même capacité de calcul. Un mode de réalisation de l'invention sera décrit en détail, en référence aux dessins annexés, dans lesquels: - la Figure 1 est une représentation schématique d'un appareil de mesure ellipsométrique à une longueur d'onde selon l'invention; - la Figure 2 est une représentation schématique d'un appareil de mesure ellipsométrique à plusieurs longueurs d'onde selon l'invention; - la Figure 3 est une représentation schématique d'un dispositif de contrôle d'élaboration de couches sur un substrat selon l'invention; - la Figure 4 est un organigramme du procédé de mesures ellipsométriques de paramètres physiques selon l'invention. L'ellipsomètre à une longueur d'onde, représenté sur la Figure 1, est destiné à la mesure de paramètres physiques d'un échantillon 1. Il comporte un groupe d'excitation 2, un groupe d'analyse 3 et une unité de traitement électronique 4. Le groupe d'excitation 2 comporte une source lumineuse 20 reliée à un système optique 22 par l'intermédiaire d'une fibre optique 21, le système optique 22 dirigeant un faisceau lumineux émis par la source 20 vers l'échantillon 1. Le groupe d'excitation 2 comporte également un polariseur 23 suivi d'un modulateur de phase 24 entre le système optique 22 et l'échantillon 1. Il introduit un déphasage 6(f). Le groupe d'analyse 3 comporte un analyseur 30 analysant un faisceau réfléchi par l'échantillon 1, suivi d'un système optique 31 adressant le faisceau réfléchi sur un monochromateur 33 par l'intermédiaire d'une fibre optique 32. Le monochromateur 33 est relié à un photodétecteur 34 qui transforme l'intensité d'un flux détecté en un signal électrique. Ce signal est représenté par I(t). Ce signal est foumi à l'unité de traitement électronique 4 par le biais d'une connexion 5. L'unité de traitement électronique 4 reçoit également une référence de fréquence et de phase du modulateur de phase 24 par une ligne 35. En fonctionnement, la source lumineuse 20 émet un faisceau de lumière incidente 10 dans une gamme de longueur d'onde donnée, ce faisceau étant polarisé par le polariseur 23, puis subissant une modulation par le modulateur de phase 24. Le modulateur de phase 24 consiste typiquement en une barre de silice fondu soumise à une contrainte périodique générée par un transducteur piézo-électrique. On crée ainsi un déphasage 8(t) modulé avec le temps t entre deux axes propres de cette barre. On module ainsi la polarisation de rayons lumineux émergents. Le faisceau lumineux incident 10 polarisé et modulé devient, après réflexion sur l'échantillon 1, un faisceau réfléchi 11 ayant une amplitude et une phase résultant de propriétés physiques de l'échantillon 1. Ce faisceau réfléchi 11 est analysé dans l'analyseur 30, puis une longueur d'onde X est sélectionnée par le monochromateur 33. Un flux du faisceau lumineux 11 est mesuré par le photodétecteur 34 qui foumit un signal électrique généré par l'intensité du flux à l'unité de traitement 4. L'ellipsomètre à plusieurs longueurs d'onde, représenté sur la Figure 2, diffère de l'ellipsomètre à une longueur d'onde de la Figure 1 par le groupe d'analyse 3 et les connexions entre ce demier et l'unité 4 de traitement électronique. Au lieu du monochromateur 33 et du photodétecteur 34, le groupe d'analyse 3 comprend un spectrographe 36 et une série de photodétecteurs 37. Chacun des photodétecteurs 37 permet la mesure d'une longueur d'onde et est relié à l'unité 4 par une connexion 6. L'ellipsomètre à plusieurs longueurs d'onde de la Figure 2 permet ainsi de mesurer simultanément plusieurs longueurs d'onde, un multiplexage pouvant être réalisé dans l'unité 4. L'ellipsomètre spectroscopique à modulation de phase, représenté sur les Figures 1 ou 2, est avantageusement utilisé dans un dispositif de contrôle d'élaboration de couches sur un substrat, représenté sur la Figure 3. L'échantillon 1 consiste en un substrat sur lequel on fait croître un dépôt par une technique connue. L'échantillon 1 est positionné sur un suscepteur 41 et contenu dans une chambre 40. L'appareil de mesure ellipsométrique précédemment décrit est utilisé pour contrôler la croissance de couches sur le substrat. Des fenêtres 44 et 45 sont pratiquées sur des chemins optiques conduisant, respectivement, du groupe d'excitation 2 vers la chambre 40, et de la chambre 40 vers le groupe d'analyse 3. Un faisceau lumineux incident 10 est ainsi directement conduit à la chambre 40, puis le faisceau réfléchi 11 par l'échantillon 1 est transmis au groupe d'analyse 3. Par ce moyen, on effectue des mesures in situ qui ne perturbent pas le processus de croissance. Au lieu d'une croissance de couches, l'élaboration peut consister en une gravure. Le signal reçu par l'unité de traitement 4, après détection du faisceau réfléchi 11, est ensuite soumis aux opérations qui suivent, schématisées sur la Figure 4. Tout d'abord, on extrait d'un signal 50 représentatif d'une intensité I(t) mesurée, des valeurs mesurées Iom, ISm et Icm par un procédé en soi connu. Iom, Ism( et Icm, connus sous le nom de coefficients de la matrice de Mueller dans le cas d'une réflexion spéculaire, sont reliés à I(P) par la relation I(0 = IOm + Ism sin 6 Icm cos 6(f) On calcule, dans une première étape 51, des composantes de Fourier So, S1 et S2 de l'intensité t), respectivement continue, à la pulsation w de modulation et à la pulsation 2w. On extrait classiquement ces composantes, par exemple par une transformée de Fourier discrète. On déduit, dans une étape 52, les valeurs mesurées Iom, ISm et Icm à partir des composantes So, S1 et S2 par des relations linéaires connues. Selon l'invention, il n'est alors pas nécessaire de connaître w et A, ou p. On applique directement aux valeurs mesurées Iom, Ism, et Icm, un processus itératif 53 qui conduit à une évaluation 54 de paramètres physiques représentatifs de l'échantillon 1. Ces paramètres physiques consistent typiquement en l'épaisseur d d'un film déposé sur un substrat, son indice n de réfraction et sa fonction diélectrique. Selon le processus itératif 53, on utilise des estimations initiales 56 des paramètres physiques en entrée 55 de paramètres. On en déduit des valeurs théoriques ISt/Iot et Ict/Iot initiales, selon des formules connues. Les valeurs théoriques Ist/Iot et Ict/Iot s'expriment en effet comme des combinaisons de fonctions trigonométriques de v et A. Par ailleurs, les angles v et A sont reliés aux paramètres physiques recherchés par les équations de Fresnel. Les valeurs théoriques Isthot et Ictîlot sont donc directement exprimables en fonction des paramètres physiques. Les valeurs théoriques obtenues dans l'étape 57 sont comparées aux valeurs mesurées IsmlIom et Ia"/Iom déterminées dans l'étape 52. Cette comparaison 58 nécessite de quantifier l'écart entre les valeurs théoriques et mesurées. On utilise pour ce faire une méthode de moindres carrés en soi connue. D'autre part, on effectue une moyenne sur une série de M points. En indexant les points de mesure par j, on obtient classiquement pour le carré de l'écart entre valeurs théoriques et mesurées: On peut utiliser une autre méthode d'ajustement que les moindres carrés, telle que le maximum de vraisemblance ou le remplacement des carrés par d'autres puissances. L'écart calculé est comparé à un seuil prédéterminé en fonction de la précision souhaitée. Si l'écart est inférieur au seuil, l'approximation des valeurs mesurées est satisfaisante et les paramètres physiques utilisés pour les calculer donnent l'évaluation 54. Dans le cas contraire, l'erreur est trop importante et une nouvelle itération est nécessaire. La comparaison 58 conduit alors à définir des estimations suivantes 59 des paramètres physiques. La détermination de ces estimations suivantes relève de méthodes connues permettant de minimiser l'écart, telle que la méthode de Levenberg-Marquardt. Ces estimations suivantes 59 sont elles-mêmes utilisées comme entrée 55 de paramètres pour calculer des valeurs théoriques Istîlot et Ictîlot suivantes en réitérant l'étape 57. Ces valeurs théoriques suivantes servent elles-mêmes de base à la comparaison 58 avec les valeurs mesurées IsmlIom et Icrnnom- boucle 59, 55, 57, 58 jusqu'à ce que l'écart obtenu lors de la comparaison 58 soit inférieur au seuil prédéterminé. De la Nième estimation des paramètres physiques associés à ces valeurs théoriques, on déduit alors l'évaluation 54 des paramètres physiques. Ce procédé de mesure ellipsométrique selon l'invention affranchit de problèmes posés par la méthode connue reposant sur l'inversion directe des équations de Fresnel, tels que pour des épaisseurs et pour W proche de 45" D'autre part, ce procédé selon l'invention ne nécessite pas de calculer w et A, mais permet d'obtenir plus directement les paramètres physiques recherchés. Généralement, trois ou quatre itérations suffisent à obtenir les paramètres, pour chaque mesure. L'utilisation d'un micro-ordinateur "PC 486" permet ainsi de déterminer l'épaisseur d et l'indice de réfraction n d'une couche déposée sur un substrat en moins de deux secondes. Typiquement, le procédé selon l'invention permet des mesures avec une précision de 3 % jusqu'à une épaisseur de couche de 5 à 6000 A, pour une vitesse de dépôt élevée de l'ordre de 30 A s1 La précision de mesure décroît avec l'épaisseur de la couche déposée. D'autre part, le procédé selon l'invention permet d'obtenir une très bonne précision sur le produit n x d. Cette précision est de l'ordre de 1 % dans le cas évoqué précédemment. Le procédé selon l'invention permet de contrôler l'homogénéité d'une couche déposée sur un substrat en épaisseur et en composition. Il est valable aussi bien pour le dépôt de matériau absorbant que pour celui de matériau transparent. Il permet également le contrôle de l'élaboration de systèmes multicouches comprenant des empilements. REVENDICATIONS 1. Procédé de mesure ellipsométrique de paramètres physiques représentatifs d'un échantillon (1) dans lequel: - on polarise linéairement un faisceau lumineux incident, ledit faisceau étant défini par un vecteur polarisation, - on module le faisceau lumineux incident de lumière polarisée de façon à introduire, entre des composantes perpendiculaires du vecteur polarisation, un déphasage 6(0 dépendant du temps t selon une variation périodique de pulsation o, 5(0 étant proportionnel au premier ordre à sin (cut), - on éclaire l'échantillon (1) avec le faisceau lumineux incident de lumière polarisée modulée, - on analyse le vecteur polarisation du faisceau lumineux renvoyé par l'échantillon (1), - on mesure par au moins un photodétecteur (34, 37) un flux du faisceau lumineux, - on effectue par une unité de traitement électronique (4) reliée au photodétecteur (34, 37), des calculs sur des mesures du flux, le nux détecté ayant une intensité I(t) de la forme: I(f) = Iom + Ism sin #(t)+Icm cos #(t), Iom, Ism, et Icm étant des valeurs mesurées dans l'unité de traitement (4) à partir de l'intensité I(fl, et dépendant desdits paramètres physiques, caractérisé en ce que:: - dans une première étape, on produit des valeurs théoriques Istflot et 1ctt initiales à partir d'estimations initiales des paramètres physiques, - lesdites valeurs théoriques servant à déterminer, dans une deuxième étape, des estimations suivantes des paramètres physiques dont on déduit des valeurs théoriques Ist/lot et Ict/Iot suivantes, - la deuxième étape étant réitérée jusqu'à une Nième estimation des paramètres physiques, de façon à minimiser l'écart entre les valeurs théoriques et celles mesurées, les paramètres physiques étant évalués à partir des valeurs Ist/Iot et Ict/Iot au cours de la Nième estimation. 2. Procédé de mesure selon la revendication 1, caractérisé en ce que: des premiers moyens de calcul de l'unité de traitement (4) donnent des composantes de Fourier So, S1 et S2 de ladite intensité du flux, respectivement continue, à la pulsation w et à la pulsation 2w, des deuxièmes moyens de calcul permettant de déduire linéairement des composantes So, S1 et S2 les valeurs mesurées Iom, Ism, et1cm. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que l'échantillon (1) comporte au moins une couche transparente éclairée par un faisceau incident. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que ledit écart entre valeurs théoriques et valeurs mesurées est minimisé par une méthode de moindres carrés. 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le procédé de mesure est spectroscopique. 6. Ellipsomètre comportant: - une source lumineuse (1) émettant un faisceau lumineux incident, - un polariseur (23) polarisant linéairement le faisceau lumineux, - un modulateur de phase (24) générant un déphasage 6(4 dépendant au premier ordre du temps t selon une variation périodique de pulsation o, t) étant proportionnel au premier ordre à sin(wt), -un analyseur (30) analysant l'état de polarisation du faisceau lumineux renvoyé par un échantillon éclairé par le faisceau lumineux incident de lumière polarisée, - au moins un photodétecteur (34, 37) mesurant un flux de faisceau lumineux, et - une unité de traitement électronique (4) reliée au photodétecteur (34, 37), effectuant des calculs sur des mesures du flux, le flux détecté ayant une intensité I(0 de la forme: I(0 = Iom + Ism sin 6(fol +1cm cos t), Iom, Ism, et Icm étant des valeurs mesurées dans l'unité de traitement (4) à partir de l'intensité t), et dépendant des paramètres physiques, caractérisé en ce que:: - L'unité de traitement produit des valeurs théoriques IstlIot et ICt/Iot initiales à partir d'estimations initiales des paramètres physiques, - lesdites valeurs théoriques servant à déterminer des estimations suivantes des paramètres physiques dont on déduit des valeurs théoriques Ist/Iot et Ictilot suivantes, - cette opération étant réitérée jusqu'à une Nième estimation des paramètres physiques, de façon à minimiser l'écart entre les valeurs théoriques et celles mesurées, les paramètres physiques étant évalués à partir des valeurs 1st11ot et Ictîlot au cours de la Nième estimation. 7. Ellipsomètre selon la revendication 6 caractérisé en ce qu'il comporte au moins une fibre de l'ensemble formé par une première et une deuxième fibres optiques, la première fibre optique étant placée entre la source et le polariseur et la deuxième fibre optique étant placée entre l'analyseur et un système de détection comprenant le photodécteur. 8. Dispositif de contrôle d'élaboration de couches comportant des moyens de contrôle de paramètres physiques représentatifs de l'élaboration de couches, le substrat et les couches constituant un échantillon (1), lesdits moyens de contrôle comprenant: - une source lumineuse (1) émettant un faisceau lumineux incident, - un polariseur (23) polarisant linéairement le faisceau lumineux, -un modulateur de phase (24) générant un déphasage 8( dépendant au premier ordre du temps t selon une variation périodique de pulsation o, 8( étant proportionnel au premier ordre à sin(wt), - un analyseur (30) analysant l'état de polarisation du faisceau lumineux renvoyé par l'échantillon éclairé par le faisceau lumineux incident de lumière polarisée, - au moins un photodétecteur (34, 37) mesurant un flux de faisceau lumineux, et - une unité de traitement électronique (4) reliée au photodétecteur (34, 37), effectuant des calculs sur des mesures du flux, le flux détecté ayant une intensité I(t) de la forme: I(f) = Iom + Ism sin t) + Icm cos 8(fl, Iom, Ism, et Icm étant des valeurs mesurées dans l'unité de traitement (4) à partir de l'intensité I(t), et dépendant des paramètres physiques, caractérisé en ce que:: - L'unité de traitement produit des valeurs théoriques ISt/Iot et Ictn(ot initiales à partir d'estimations initiales des paramètres physiques, - lesdites valeurs théoriques servant à déterminer des estimations suivantes des paramètres physiques dont on déduit des valeurs théoriques Ist/Iot et ICt/lot suivantes, - cette opération étant réitérée jusqu'à une Nième estimation des paramètres physiques, de façon à minimiser l'écart entre les valeurs théoriques et celles mesurées, les paramètres physiques étant évalués à partir des valeurs Istilot et Ictîlot au cours de la Nième estimation. 9. Dispositif selon l'une des revendications 7 ou 8, caractérisé en ce que les couches ont chacune une épaisseur, un indice de réfraction et une fonction diélectrique, les paramètres physiques comprenant au moins certains d'entre eux. 10. Dispositif selon l'une quelconque des revendications 7 à 9, caractérisé en ce que, I'élaboration étant un dépôt1 le dépôt d'une couche a une vitesse supérieure à 30 A s-1.