FR 3000234 A1 20140627 FR 1262521 A 20121221 La présente invention concerne un procédé d'estimation de motifs à imprimer sur plaque ou sur masque par lithographie à faisceau d'électrons. Elle concerne également un dispositif d'impression mettant en oeuvre ce procédé. L'invention s'applique plus particulièrement à un procédé de ce type comportant les étapes suivantes : impression par lithographie à faisceau d'électrons, dans une résine disposée sur un substrat, d'un ensemble de motifs de calibrage, mesure de dimensions caractéristiques de l'ensemble de motifs de calibrage après impression et enregistrement de ces dimensions caractéristiques en mémoire, fourniture, par un processeur ayant accès à la mémoire, d'une estimation de la fonction d'étalement du point caractérisant des électrons diffusés lors de l'impression à partir des dimensions caractéristiques mesurées, estimation des motifs à imprimer par convolution de la fonction d'étalement du point fournie avec une valeur initiale de ces motifs à imprimer. La lithographie par faisceau d'électrons ne peut pas s'affranchir d'une phase de préparation d'un motif à imprimer compte tenu des effets de proximité engendrés par la diffusion des électrons, d'autres espèces chimiques dans la résine, ou autres effets susceptibles de déformer le motif à imprimer et qu'il est nécessaire de compenser. Notamment, l'exposition du substrat au faisceau d'électrons engendre des électrons diffusés, incluant des électrons diffusés dans la résine et des électrons rétrodiffusés par le substrat, produisant ce que l'on appelle « l'effet de proximité ». Cet effet provoque un élargissement du motif final imprimé sur la résine par rapport au motif initialement souhaité. L'effet est d'autant plus gênant que les motifs à placer sont fins et se trouvent dans une zone où la densité de motifs à imprimer est grande. D'une façon générale, pour pouvoir compenser cet effet, une méthode connue consiste à modéliser la fonction d'étalement du point, encore appelée réponse percussionnelle, dite PSF (de l'anglais « Point Spread Function »), comportant la PSF des électrons rétrodiffusés par le substrat, dite PSFBE (de l'anglais « Backward scattered Electrons Point Spread Function ») et la PSF des électrons diffusés dans la résine, aussi connue de l'homme du métier comme la « PSF des électrons diffusés vers l'avant » et dite PSFFE (de l'anglais « Forward scattered Electrons Point Spread Function »). Cette fonction PSF = PSFBE + PSFFE permet ensuite de calculer, lors d'une étape de simulation dite d'estimation, la dose finale à appliquer localement sur la résine pour l'impression du motif. Cette dose finale est estimée moyennant la convolution de la fonction PSF obtenue avec une valeur initiale du motif à imprimer. La dose ainsi obtenue tient compte de la dose parasite induite par l'effet de proximité sur le motif. Une première solution consiste plus précisément à modéliser la forme de la PSF par une somme de fonctions de type gaussienne. Afin de déterminer les coefficients caractérisant ces fonctions gaussiennes, différents motifs de calibrage, par exemple des lignes, sont imprimés sur la résine en appliquant différentes doses de charges électriques au faisceau d'électrons. Des mesures de dimensions caractéristiques sont effectuées sur les motifs de calibrage imprimés et utilisées pour déterminer les paramètres de la fonction gaussienne par simulation à l'aide d'un programme informatique dédié et bien connu de l'homme du métier sous le nom de « solver ». Néanmoins, les doses ainsi calculées sont imprécises et empêchent d'imprimer le motif avec une précision satisfaisante puisqu'en réalité la forme de la PSFBE, n'est pas vraiment gaussienne. En effet, elle dépend de la distribution du nombre d'électrons rétrodiffusés et de la distribution de leur énergie qui ne suivent pas une loi gaussienne. Il peut ainsi être souhaité de prévoir un procédé d'estimation de motifs à imprimer du type précité qui permette de s'affranchir d'au moins une partie des problèmes et contraintes précités, notamment qui permette de ne pas imposer de forme particulière à la PSFBE des électrons rétrodiffusés lors de l'impression et de s'affranchir de l'utilisation d'un « solver » pour le calcul de cette PSFBE. Il est donc proposé un procédé d'estimation de motifs à imprimer sur plaque ou sur masque par lithographie à faisceau d'électrons, comportant les étapes suivantes : impression par lithographie à faisceau d'électrons, dans une résine disposée sur un substrat, d'un ensemble de motifs de calibrage, mesure de dimensions caractéristiques de l'ensemble de motifs de calibrage après impression et enregistrement de ces dimensions caractéristiques en mémoire, fourniture, par un processeur ayant accès à la mémoire, d'une estimation de la fonction d'étalement du point caractérisant des électrons diffusés lors de l'impression à partir des dimensions caractéristiques mesurées, estimation des motifs à imprimer par convolution de la fonction d'étalement du point fournie avec une valeur initiale de ces motifs à imprimer, dans lequel : chaque motif de calibrage imprimé dans la résine comporte une zone centrale exposée au faisceau d'électrons et une pluralité de zones concentriques à symétrie de révolution disposées autour de cette zone centrale, les dimensions caractéristiques mesurées après impression sont des dimensions caractéristiques des zones centrales des motifs, et l'estimation de la fonction d'étalement du point est calculée en tant que somme d'une première portion de fonction d'étalement du point caractérisant des électrons rétrodiffusés par le substrat et d'une seconde portion de fonction d'étalement du point caractérisant des électrons diffusés dans la résine, et comporte une inversion d'une modélisation analytique de l'effet, sur lesdites dimensions caractéristiques, de l'application de la première portion de fonction d'étalement du point à l'ensemble de motifs de calibrage. Grâce à l'invention, le calcul de la PSF des électrons diffusés dans la résine et rétrodiffusés par le substrat ne fait pas d'hypothèse sur la forme de la portion PSFBE (si ce n'est sa symétrie de révolution) et les doses ainsi calculées lors de l'étape d'estimation des motifs à imprimer sont suffisamment précises pour parvenir à imprimer les motifs souhaités avec une précision améliorée et en ayant compensé les effets de proximité. Outre la précision du résultat, le calcul de la PSFBE est aussi accéléré, car il résulte directement des mesures des dimensions caractéristiques des motifs de calibrage de sorte qu'aucun programme de simulation numérique (« solver ») n'est nécessaire. De façon optionnelle, l'application de la première portion de fonction d'étalement du point à l'ensemble de motifs de calibrage est modélisée analytiquement par le produit d'une matrice modélisant l'ensemble de motifs de calibrage avec un vecteur modélisant la première portion de fonction d'étalement du point, l'inversion de la modélisation analytique étant réalisée par une inversion de la matrice modélisant l'ensemble de motifs de calibrage. De façon optionnelle également, la matrice modélisant l'ensemble de motifs de calibrage est une matrice carrée d'ordre n modélisant un ensemble de n motifs à n zones concentriques, chaque ligne de la matrice correspondant à l'un des motifs de calibrage et chaque colonne de la matrice correspondant à l'une de ces zones concentriques. De façon optionnelle également, la matrice modélisant l'ensemble de motifs de calibrage est exprimée à partir d'une matrice de Hadamard. De façon optionnelle également, la première portion de fonction d'étalement du point caractérisant les électrons rétrodiffusés par le substrat, notée PSFBE, est définie par l'expression suivante : 1 0 CD[l] PSFBE oc H - 0 0I ['1 0. 0 [1.1\ . Cdo , où le symbole « » désigne une - \ CD[n] ° S[n] 1 relation de proportionnalité, H est la matrice de Hadamard à partir de laquelle CD [1]S[1] l'ensemble de motifs de calibrage est modélisé, : est un vecteur comportant CD[n] S[n] l'ensemble [CD[1], ...,CD[n]l des dimensions caractéristiques mesurées, pondéré par l'ensemble {S[1], ...S[n]l des surfaces respectives des zones concentriques à symétrie de révolution, et cdo est une constante caractéristique d'une dimension commune de la zone centrale de chaque motif de calibrage lorsque cette zone centrale est exposée de manière isolée au faisceau d'électrons en l'absence de zones concentriques à symétrie de révolution autour de cette zone centrale. De façon optionnelle également, la seconde portion de fonction d'étalement du point est modélisée à partir d'une fonction de type gaussienne, les paramètres de cette fonction gaussienne étant déterminés par simulation sur une pluralité de motifs de calibrage. De façon optionnelle également, les zones concentriques à symétrie de révolution des motifs de calibrage sont toutes de même surface. De façon optionnelle également, la première portion de fonction d'étalement du point caractérisant les électrons rétrodiffusés par le substrat est définie par l'expression suivante : 1 0 0 (1-CD[1]] PSFBE oc H [(.) 0 ?I) [C1:) [n]] cd° 1:1] 0 0 De façon optionnelle également, la zone centrale de chaque motif de calibrage comporte l'une des zones concentriques à symétrie de révolution, celle-ci étant alors placée au centre dudit motif. Il est également proposé un dispositif d'impression de motifs sur plaque ou sur masque par lithographie à faisceau d'électrons, comportant : des moyens d'impression par lithographie à faisceau d'électrons, dans une résine disposée sur un substrat, d'un ensemble de motifs de calibrage, des moyens de mesure de dimensions caractéristiques de l'ensemble de motifs de calibrage après impression et des moyens de stockage de ces dimensions caractéristiques en mémoire, un processeur ayant accès à la mémoire et programmé pour fournir une estimation de la fonction d'étalement du point caractérisant des électrons diffusés lors de l'impression à partir des dimensions caractéristiques mesurées, des moyens d'estimation des motifs à imprimer par convolution de la fonction d'étalement du point fournie avec une valeur initiale de ces motifs à imprimer, dans lequel, chaque motif de calibrage destiné à être imprimé dans la résine comportant une zone centrale exposée au faisceau d'électrons et une pluralité de zones concentriques à symétrie de révolution disposées autour de cette zone centrale : les moyens de mesure sont des moyens de mesure des dimensions caractéristiques des zones centrales des motifs après impression, et le processeur est programmé pour calculer l'estimation de la fonction d'étalement du point en tant que somme d'une première portion de fonction d'étalement du point caractérisant des électrons rétrodiffusés par le substrat et d'une seconde portion de fonction d'étalement du point caractérisant des électrons diffusés dans la résine, et par inversion d'une modélisation analytique de l'effet, sur lesdites dimensions caractéristiques, de l'application de la première portion de fonction d'étalement du point à l'ensemble de motifs de calibrage. L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels : la figure 1 représente schématiquement la structure générale d'un dispositif d'impression, selon un mode de réalisation de l'invention, La figure 2 illustre un ensemble de motifs de calibrage utilisé pour estimer la fonction d'étalement du point caractérisant des électrons diffusés lors de l'impression d'un motif par le dispositif d'impression de la figure 1, et - la figure 3 illustre les étapes successives d'un procédé d'estimation mis en oeuvre par le dispositif d'impression de la figure 1 à l'aide de l'ensemble de motifs de la figure 2. Le dispositif d'impression 10 représenté schématiquement sur la figure 1 comporte un module de traitement 12, un module 14 d'impression de motifs sur plaque ou sur masque par lithographie à faisceau d'électrons et un module 15 de mesure de dimensions caractéristiques de motifs imprimés par le module 14, tous les trois associés de façon classique à une mémoire 16 (par exemple une mémoire RAM). Le module de traitement 12 peut par exemple être mis en oeuvre dans un dispositif informatique tel qu'un ordinateur classique comportant un processeur associé à une ou plusieurs mémoires pour le stockage de fichiers de données et de programmes d'ordinateur. Le module de traitement 12 peut alors lui-même être considéré comme formé d'un processeur associé à une mémoire de stockage des instructions qu'il exécute sous forme de programmes d'ordinateur. Le module de traitement 12 tel qu'illustré sur la figure 1 comporte ainsi fonctionnellement quatre programmes d'ordinateur 18, 20, 22 et 24. Le premier programme d'ordinateur 18 est un programme de préparation préalable d'un ensemble de motifs de calibrage Mc à imprimer par lithographie à faisceau d'électrons dont les paramètres sont stockés en mémoire 16. L'ensemble de motifs de calibrage Mc est par exemple plus précisément modélisé par un ensemble de formes géométriques McF affectées à un ensemble de doses prédéterminées DcF. Par « dose », on entend une dose de charges électriques à appliquer au faisceau d'électrons lors de l'impression individuelle de la forme géométrique qui lui est associée. A partir de l'ensemble McF de formes géométriques et de l'ensemble DcF de doses prédéterminées enregistrées en mémoire 16, ces paramètres définissant théoriquement l'ensemble de motifs de calibrage Mc à imprimer, le premier programme 18 est conçu pour fournir un ensemble de formes géométriques élémentaires McFe modélisant les motifs de calibrage Mc et un ensemble associé de doses DcFe choisies parmi un ensemble de valeurs discrètes, par exemple des valeurs discrètes de l'ensemble DcF si celui-ci est lui-même à valeurs discrètes. Ces ensembles McFe et DcFe sont par exemple stockés en mémoire 16. Une première solution bien connue consiste plus précisément à modéliser l'ensemble de motifs de calibrage Mc par un ensemble de formes géométriques élémentaires McFe juxtaposées et à affecter à chaque forme élémentaire une dose et un temps d'exposition indépendants des autres. Mais comme en pratique les doses doivent être choisies dans un ensemble Du discret et limité de doses prédéterminées, il est alors souvent nécessaire de multiplier les formes géométriques élémentaires et de prévoir certaines formes de très petites tailles pour compenser la discrétisation des doses et obtenir une bonne précision d'impression. Une autre solution est par exemple décrite dans l'article de C. Pierrat et al, intitulé « Mask data correction methodology in the context of model-base fracturing and advanced mask models », publié dans Proceedings SPIE of Optical Microlithography n° 24, vol. 7973, 1-3 mars 2011, &In Jose (CA), US. Cette solution autorise la superposition des formes géométriques élémentaires de sorte que l'on peut ainsi très sensiblement réduire leur nombre et limiter les formes géométriques de petites tailles. Mais la superposition des formes rend les autres paramètres, notamment les doses de charges électriques et les temps d'exposition, dépendants les uns des autres. Une troisième solution plus performante et originale consiste à définir les formes géométriques élémentaires McFe sous la forme de pixels identiques juxtaposés en matrice, formant ainsi un pavage bidimensionnel régulier recouvrant l'ensemble de motifs de calibrage Mc, une correction d'erreur de discrétisation pouvant alors être réalisée par tramage sur ces pixels lors de la détermination des doses DcFe à appliquer aux formes géométriques élémentaires. Le module d'impression 14 tel qu'illustré sur la figure 1 permet d'imprimer par exposition au faisceau d'électrons d'une résine disposée sur un substrat, l'ensemble de motifs de calibrage Mc à partir de l'ensemble McFe de formes géométriques élémentaires et de l'ensemble DcFe associé de doses. Le module de mesure 15 tel qu'illustré sur la figure 1 permet alors de mesurer des dimensions caractéristiques CD de l'ensemble de motifs de calibrage Mc après impression par le module 14. Comme cela sera détaillé ultérieurement en référence à la figure 2, chaque motif de calibrage imprimé dans la résine comporte une zone centrale exposée au faisceau d'électrons et une pluralité de zones concentriques à symétrie de révolution disposées autour de cette zone centrale. Les dimensions caractéristiques CD mesurées après impression correspondent alors aux dimensions des zones centrales des motifs. A titre d'exemple, lorsque la zone centrale d'un motif est un disque, la dimension caractéristique de ce motif est la dimension du diamètre de ce disque. Les dimensions caractéristiques CD de l'ensemble de motifs de calibrage Mc sont par exemple enregistrées en mémoire 16. Le deuxième programme d'ordinateur 20 du module de traitement 12 est un calculateur de la fonction d'étalement du point caractérisant des électrons rétrodiffusés PSFBE dont le fonctionnement sera détaillé en référence à la figure 3. A partir des dimensions caractéristiques CD mesurées et enregistrées en mémoire, ce programme d'ordinateur 20 fournit une estimation de la fonction d'étalement du point PSFBE caractérisant les électrons rétrodiffusés par le substrat lors de l'impression. Cette estimation de la fonction PSFBE est calculée à partir d'une modélisation analytique de l'effet, sur les dimensions caractéristiques CD de l'application de la fonction PSFBE à l'ensemble de motifs de calibrage M. Plus précisément, l'effet, sur les dimensions caractéristiques CD, de l'application de la fonction PSFBE à l'ensemble de motifs de calibrage Mc est modélisé analytiquement par le produit d'une matrice M modélisant l'ensemble de motifs de calibrage Mc et d'un vecteur modélisant la fonction PSFBE. A titre d'exemple, la ième composante de ce vecteur correspond à la quantité d'électrons rétrodiffusés par la ième zone concentrique. Par extension, ce vecteur sera également noté PSFBE, et l'effet précité est défini selon l'expression suivante : CD - cdo oc M *PSFBE, où le symbole « » désigne une relation de proportionnalité et cdo est une constante caractéristique d'une dimension commune de la zone centrale de chaque motif de calibrage lorsque cette zone centrale est exposée de manière isolée au faisceau d'électrons en l'absence de zones concentriques à symétrie de révolution autour de cette zone centrale. Cette constante cdo peut être préétablie et stockée en mémoire 16. On peut donc constater que l'effet de la fonction PSFBE sur l'ensemble de motifs de calibrage Mc est proportionnel à l'élargissement des dimensions caractéristiques CD des zones centrales du motif imprimé, cet élargissement étant provoqué par l'effet de proximité dû à la présence des zones concentriques imprimées autour. L'inversion de cette modélisation analytique, réalisée par inversion de la matrice M modélisant l'ensemble de motifs de calibrage Mc, permet ainsi d'obtenir une estimation de la fonction PSFBE selon l'expression suivante : PSFBE oc M * (CD - cdo). On note que cette inversion impose que la matrice M modélisant l'ensemble de motifs de calibrage soit une matrice carrée et inversible, ce qui impose également quelques contraintes à cet ensemble de motifs de calibrage M. La fonction PSFBE ainsi obtenue est par exemple stockée en mémoire 16. Le troisième programme d'ordinateur 22 du module de traitement 12 est un programme d'estimation de la fonction d'étalement du point PSF caractérisant des électrons diffusés lors de l'impression. Cette fonction est estimée en tant que somme de la première portion précitée de fonction d'étalement du point PSFBE caractérisant les électrons rétrodiffusés par le substrat et d'une seconde portion de fonction d'étalement du point PSFFE caractérisant des électrons diffusés dans la résine. Contrairement à la fonction PSFBE, le comportement de la fonction PSFFE peut être modélisé par une fonction de type gaussienne. De ce fait, la seconde portion PSFFE peut être considérée comme une fonction gaussienne, les paramètres de cette fonction gaussienne étant déterminés par simulation sur une pluralité de motifs de calibrage et à l'aide d'un logiciel dédié. Cette méthode bien connue de l'homme du métier permet ainsi d'obtenir une estimation de la fonction PSFFE suffisamment précise. Les fonctions ainsi calculées, PSFFE et PSF = PSFBE + PSFFE, sont par exemple stockées en mémoire 16. Le quatrième programme d'ordinateur 24 du module de traitement 12 est un programme d'estimation de motifs à imprimer par le dispositif d'impression 10, par convolution de la fonction d'étalement du point PSF, fournie par le programme d'ordinateur 22, avec une valeur initiale de ces motifs à imprimer. Un motif à imprimer comporte par exemple un ensemble de formes géométriques à imprimer, de valeur initiale MpF, associé à un ensemble de doses prédéfinies, de valeur initiale DpF, à appliquer au faisceau d'électrons lors de l'impression dudit motif. Le résultat de cette convolution donne une valeur finale (MF,DF) du motif à imprimer sur lequel les effets des doses parasites induites par l'effet de proximité ont été compensés. Le motif à imprimer (MF,DF) ainsi calculé est ensuite utilisé lors de la modélisation du motif en formes géométriques élémentaires MpFe à partir de l'ensemble de formes géométriques à imprimer MF et plus précisément lors du calcul, par un programme similaire au programme de préparation préalable 18 ou par ce programme lui-même, de l'ensemble de doses DpFe associées audites formes géométriques élémentaires MpFe. Enfin, l'impression sur la résine du motif (M'pF, D'pF) est effectuée par le module d'impression 14 à partir de l'ensemble MpFe de formes géométriques élémentaires et de l'ensemble associé de doses DpFe avec une précision améliorée et en ayant compensé les effets de proximité. On notera par ailleurs que les programmes d'ordinateur 18, 20, 22 et 24 sont présentés comme distincts, mais cette distinction est purement fonctionnelle. Ils pourraient tout aussi bien être regroupés en un ou plusieurs logiciels. Leurs fonctions pourraient aussi être au moins en partie micro programmées ou micro câblées dans des circuits intégrés dédiés. Ainsi, en variante, le dispositif informatique mettant en oeuvre le module de traitement 12 pourrait être remplacé par un dispositif électronique composé uniquement de circuits numériques (sans programme d'ordinateur) pour la réalisation des mêmes actions. La figure 2 illustre à titre d'exemple un ensemble de motifs de calibrage Mc compatible avec une modélisation matricielle inversible selon l'invention. Chaque motif de calibrage Mc [d de l'ensemble de n motifs de calibrage Mc = [Aida 0 En outre, le fait que la matrice de Hadamard H soit utilisée pour exprimer la matrice M modélisant l'ensemble Mc de motifs de calibrage garantit que la moitié n/2 des zones concentriques Ac à symétrie de révolution appartenant à chacun des motifs Mc[i] soit exposée au faisceau d'électrons, ceci étant valable pour tous les motifs de calibrage imprimés dans la résine sauf pour le premier motif Mc[1]. En effet, les valeurs des coefficients de la première ligne de la matrice de Hadamard H étant toutes égales à 1, toutes les zones concentriques Ac à symétrie de révolution de ce motif de calibrage Mc,[1] sont en conséquence exposées au faisceau. Le procédé d'estimation illustré sur la figure 3 et mis en oeuvre par le dispositif d'impression 10 de la figure 1 comporte une première étape 100 d'impression par lithographie à faisceau d'électrons d'un ensemble de motifs de calibrage Mc dans une résine disposée sur un substrat. Lors d'une première sous étape 102 de la première étape 100, le programme 18 de préparation préalable d'un ensemble de motifs de calibrage Mc à imprimer, après avoir extrait préalablement de la mémoire les ensembles McF de formes géométriques et Du de doses prédéterminées, fournit l'ensemble de formes géométriques élémentaires McFe et l'ensemble de doses DcFe associées. Ces ensembles McFe et DcFe sont stockés en mémoire 16. Au cours d'une deuxième sous étape 104 de la première étape 100, le dispositif 14 d'impression par lithographie à faisceau d'électrons imprime l'ensemble de motifs de calibrage Mc par exposition d'une résine disposée sur un substrat à partir desdits ensembles McFe et DcFe préalablement extraits de la mémoire 16 par ce dispositif 14. Au cours d'une deuxième étape 120, le module de mesure 15 fournit les dimensions caractéristiques CD des zones centrales 26 de chaque motif de calibrage Mc[i] imprimé. Ce module peut également participer à la fourniture de la valeur cdo correspondant à une constante caractéristique d'une dimension commune de la zone centrale 26 de chaque motif de calibrage Mc[i] lorsque cette zone centrale 26 est exposée de manière isolée au faisceau d'électrons en l'absence de zones concentriques Ac à symétrie de révolution autour de cette zone centrale 26. La constante cdo peut en effet être obtenue expérimentalement, par exemple en fixant la dose d'exposition de la zone centrale 26 d'un motif de calibrage Mc[i] quelconque et en faisant décroître progressivement la dose d'exposition DcF des zones concentriques Ac à symétrie de révolution de ce motif. Ensuite, une solution possible consiste à tracer une courbe à partir de l'acquisition de la dimension caractéristique CD[i] de la zone centrale du motif de calibrage Mc[i] considéré pour chaque dose différente appliquée aux zones concentriques Ac autour. Enfin, la dimension théorique cdo correspondant à une dose d'exposition des zones concentriques Ac égale à zéro peut être déterminée par extrapolation de ladite courbe, en tenant compte de la continuité en zéro de la fonction PSF et de sa fonction dérivée. Au cours d'une troisième étape 140, les programmes d'ordinateur 20 et 22 du module de traitement 12 permettent d'estimer la fonction PSF d'étalement du point caractérisant les électrons diffusés lors de l'impression à partir des dimensions caractéristiques CD mesurées et de la constante cdo. Lors d'une première sous-étape 142 de la troisième étape 140, le programme d'ordinateur 20, après extraction de la mémoire 16 de l'ensemble des dimensions CD = [CD[i],0 En supposant que la distance d'un point P', appartenant à la surface du ième motif de calibrage M[i], à tout point P situé au voisinage de la zone centrale 26 de ce motif de calibrage peut être approchée par la distance entre ce point P' et le centre de la zone centrale 26 du motif de calibrage M[i], et en passant d'un système de coordonnées cartésiennes (x,y) à un système de coordonnées polaires (r, e), la dose au voisinage de la ième zone centrale D(P) peut être approchée par l'expression suivante : Di(P)= Do (p) + ri, fo2 fRRJ 114,12+1 PSFFE(r,O)rdre +rj2_1f0271- fRRJ ih,J2+1 PSFBE(r,O)rdrclO , où Do(P) correspond à la dose de la zone centrale 26 du ième motif de calibrage M[i] lorsque cette zone centrale 26 est exposée de manière isolée au faisceau d'électrons en l'absence de zones concentriques Ac à symétrie de révolution autour de cette zone centrale 26 et h,j représente les coefficients de la matrice de Hadamard H et donc l'exposition de la ième zone concentrique Ac[j] du ième motif de calibrage M[i]. hu-Fi On notera que c'est l'expression qui permet de passer de la valeur -1 d'un coefficient de la matrice de Hadamard H à la valeur 0 d'un coefficient de la matrice M. On notera ainsi que dans l'exemple de réalisation ci-dessus la matrice M carrée modélisant l'ensemble de motifs de calibrage Mc peut être exprimée à l'aide d'une matrice de Hadamard. Mais dans d'autres modes de réalisation de l'invention, une matrice différente pourrait être utilisée, cette matrice étant nécessairement cependant une matrice inversible. La dose D(P) au voisinage de la ième zone centrale ainsi obtenue tient compte de la dose parasite induite par les zones concentriques Ac appartenant au ième motif et néglige l'influence du reste des motifs autour. Elle est donc valide seulement si les différents motifs de calibrage Mc[i] sont suffisamment éloignés entre eux pour pouvoir négliger leurs influences croisées. Afin de simplifier l'expression, D(P) deux approximations sont ensuite utilisées : le support de la fonction PSFFE étant inférieur au rayon Rc[0], de la zone concentrique Ac[1] la plus proche de la zone centrale 26 d'un motif de calibrage M[i], l'influence de la PSFFE des zones concentriques Ac de ce motif Mc[i] sur sa zone centrale 26 peut être négligée, la valeur de la fonction PSFBE est constante sur la surface d'une zone concentrique Ac[j] et son amplitude dépend de la distance de ladite zone concentrique à la zone centrale 26, PSFBE(j). A titre d'exemple, PSFBE(j) pourrait être égale à la valeur moyenne de la PSFBE sur la zone concentrique Ac[j]. En tenant compte des approximations précédentes, la dose au voisinage de la ième zone centrale Di(P) peut être approchée par l'expression suivante : D(P) = Do(P)+ Sih 2+1 PSFBE(j), où S représente la surface d'une zone concentrique Ac[j] sachant que toutes les zones concentriques sont de même surface. La même expression peut être reprise sous forme matricielle : D(P) = Do(P) + S PSFBE. La dose D(P) au voisinage de la zone centrale est ainsi approchée par la somme de la dose Do(P) de la zone centrale 26 exposée isolément et d'une dose parasite constante modulée par la nature des zones concentriques Ac exposées selon l'expression : PSFBE. Un avantage, provenant du fait que la matrice M modélisant l'ensemble de motifs de calibrage Mc puisse être exprimée à l'aide d'une matrice de Hadamard H d'ordre n, est que la valeur totale de ladite dose parasite varie peu d'un motif à l'autre puisque n/2 zones concentriques Ac sont toujours exposées (à l'exception du premier motif). En considérant que la taille des zones centrales CD des motifs de calibrage varie de façon linéaire avec la dose parasite, on peut en déduire l'expression : CD oc cdo +S-PSFBE, où le signe du facteur de proportionnalité dépend de la polarité de la résine utilisée. Le facteur constant S de l'expression précédente est un coefficient de proportionnalité participant à l'effet de la fonction d'étalement du point PSFBE sur l'élargissement des dimensions caractéristiques CD des zones centrales 26 selon l'expression suivante : (CD - cdo) OC (±I2 )PSFBE. En multipliant cette expression par P) 1- de chaque côté et en appliquant 2 les propriétés des matrices de Hadamard bien connues de l'homme du métier, on obtient l'expression suivante pour la fonction PSFBE H+1 -1 (CD - cdo) = n H - 1 15 --- 0 9 0 --- 1 0 0 1) (CD - cd0[1:1) où CD = PSFBE () OC CD [1] : est le vecteur comportant l'ensemble [CD[1],...,CD[n]l des dimensions CD[n] caractéristiques mesurées. En conséquence, l'estimation de la fonction PSFBE est calculée par inversion de la modélisation analytique de l'effet, sur les dimensions caractéristiques CD, de son application à l'ensemble de motifs de calibrage M. On notera que du fait que la matrice de Hadamard H est utilisée comme matrice participant à la modélisation de l'ensemble de motifs de calibrage Mc, chaque valeur de la fonction PSFBE ainsi obtenue est proportionnelle à une moyenne de n/2 mesures de dimensions caractéristiques CD, correspondant aux n/2 zones concentriques Ac exposées au faisceau d'électrons. Ainsi, selon le théorème central - limite, la fonction PSFBE est une variable aléatoire qui tend à être normalement distribuée lorsque n augmente. L'écart type CiP.- crBE de cette variable aléatoire PSFBE est directement proportionnel à l'écart type cioo de l'ensemble CD de dimensions caractéristiques des motifs de calibrage Mc et inversement proportionnel à la racine du nombre de motifs : o-psFBE oc a CD 7,2 . Un nombre n suffisamment élevé de motifs de calibrage Mc garantira donc la précision de la fonction PSFBE ainsi calculée. De plus, cette précision est améliorée par rapport à un autre mode de réalisation possible où l'influence de chaque zone concentrique serait mesurée séparément sans faire de 20 moyenne. Lors d'une deuxième sous-étape 144 de la troisième étape 140, le programme d'ordinateur 22, après extraction de la mémoire 16 de la fonction PSFBE, fournit une estimation de la fonction d'étalement du point PSFFE caractérisant les électrons diffusés dans la résine, cette fonction PSFFE étant déterminée par simulation 25 de façon connue en soi à partir d'un modèle gaussien. Lors d'une troisième sous-étape 146 de la troisième étape 140, le programme d'ordinateur 22 fournit une estimation de la fonction d'étalement du point PSF caractérisant des électrons diffusés lors de l'impression en effectuant la somme de la fonction PSFBE et de la fonction PSFFE. 30 Au cours d'une quatrième étape 160 le programme d'ordinateur 24 du module de traitement 12 réalise l'estimation d'un motif (MPF,DpF) à imprimer. Lors d'une première sous-étape 162 de cette quatrième étape 160, après extraction de la mémoire 16 de valeurs initiales d'un motif à imprimer (MpF,DpF) et de la fonction d'étalement du point PSF, le programme d'ordinateur 24 fournit par un calcul de convolution une estimation du motif (MPF,DPF) à imprimer sur lequel les effets des doses parasites induites par l'effet de proximité ont été compensés. Lors d'une deuxième sous-étape 164 de la quatrième étape 160, le programme d'ordinateur 18 (ou un autre programme similaire), après extraction de la mémoire 16 du motif à imprimer (MPF,DPF), fournit l'ensemble de formes géométriques élémentaires MpFe et l'ensemble de doses associées DpFe modélisant ce motif et tenant compte de la dose parasite induite. Enfin, au cours d'une cinquième étape 180, le module d'impression 14 effectue l'impression sur la résine du motif (MPF,DPF) à partir de l'ensemble de formes géométriques élémentaires MpFe et de l'ensemble associé de doses DpFe en ayant compensé les effets de proximité et donc avec une précision d'impression améliorée. Bien entendu, des variantes de réalisation pour estimer la fonction PSFBE peuvent être envisagées. Ainsi, la zone centrale 26 de chaque motif de calibrage Mc[i] pourrait comporter la première zone concentrique Ac[1] à symétrie de révolution, celle-ci étant alors placée au centre du motif de calibrage. Dans ce cas et à titre d'exemple, les dimensions caractéristiques CD de la zone centrale 26 de chaque motif de calibrage pour cette variante de réalisation comporteraient la taille du plus petit rayon Rc[0] de la zone concentrique Ac[1] à symétrie de révolution placée au centre du motif de calibrage pourvu qu'elle soit écrite sur tous les motifs de calibrage. La constante caractéristique cdo correspondrait alors à la valeur de Rc[1], commune à chaque motif de calibrage M[i], lorsque la zone concentrique Ac[1] est exposée de manière isolée au faisceau d'électrons en l'absence de zones concentriques Ac = tAc[k],1 REVENDICATIONS1. Procédé d'estimation de motifs (M;,.F,D;,.F) à imprimer sur plaque ou sur masque par lithographie à faisceau d'électrons, comportant les étapes suivantes : impression (100) par lithographie à faisceau d'électrons, dans une résine disposée sur un substrat, d'un ensemble de motifs de calibrage (MCF, D 1 CF) I mesure (120) de dimensions caractéristiques (CD) de l'ensemble de motifs de calibrage (McF,DcF) après impression et enregistrement de ces dimensions caractéristiques (CD) en mémoire, fourniture, par un processeur ayant accès à la mémoire (16), d'une estimation (140) de la fonction d'étalement du point (PSF) caractérisant des électrons diffusés lors de l'impression à partir des dimensions caractéristiques (CD) mesurées, estimation (160) des motifs (M;,.F,D;,.F) à imprimer par convolution de la fonction d'étalement du point (PSF) fournie avec une valeur initiale (MpF,DpF) de ces motifs à imprimer, caractérisé en ce que : chaque motif de calibrage (M[i]) imprimé dans la résine comporte une zone centrale (26) exposée au faisceau d'électrons et une pluralité de zones concentriques à symétrie de révolution disposées autour de cette zone centrale (26), les dimensions caractéristiques (CD) mesurées après impression sont des dimensions caractéristiques (CD) des zones centrales 26 des motifs, et l'estimation de la fonction d'étalement du point (PSF) est calculée en tant que somme d'une première portion de fonction d'étalement du point (PSFBE) caractérisant des électrons rétrodiffusés par le substrat et d'une seconde portion de fonction d'étalement du point (PSFFE) caractérisant des électrons diffusés dans la résine, et comporte une inversion d'une modélisation analytique de l'effet, sur lesdites dimensions caractéristiques (CD), de l'application de la première portion de fonction d'étalement du point (PSFBE) à l'ensemble de motifs de calibrage (MCF, DCF) - 2. Procédé d'estimation de motifs (ME,DE) à imprimer selon la revendication 1, dans lequel l'application de la première portion de fonction d'étalement du point (PSFBE) à l'ensemble de motifs de calibrage 014- D 1 est 'cF, _ cp, est _ modélisée analytiquement par le produit d'une matrice modélisant l'ensemble de motifs de calibrage avec un vecteur modélisant la première portion de fonction d'étalement du point (PSFBE), l'inversion de la modélisation analytique étant réalisée par une inversion de la matrice modélisant l'ensemble de motifs de calibrage. 3. Procédé d'estimation de motifs (M;,.E,DE) à imprimer selon la revendication 2, dans lequel la matrice modélisant l'ensemble de motifs de calibrage (MEE,DEE) est une matrice carrée d'ordre n modélisant un ensemble de n motifs à n zones concentriques, chaque ligne de la matrice correspondant à l'un des motifs de calibrage (MEN) et chaque colonne de la matrice correspondant à l'une de ces zones concentriques (Ac[k]). 4. Procédé d'estimation de motifs (ME,DE) à imprimer selon la revendication 3, dans lequel la matrice modélisant l'ensemble de motifs de calibrage (McF,DcF) est exprimée à partir d'une matrice de Hadamard. 5. Procédé d'estimation de motifs (M;,.E,DE) à imprimer selon la revendication 4, dans lequel la première portion de fonction d'étalement du point caractérisant les électrons rétrodiffusés par le substrat, notée PSFBE, est définie par l'expression suivante : 1 CD [1] PSFBE oc H [(.) 0 001) 711 [11 \ 0 . : - cdo , où le symbole « « » désigne une \ CD [n] 0 0 \ S[n] 1 relation de proportionnalité, H est la matrice de Hadamard à partir de laquelle C D [1] S[1] l'ensemble de motifs de calibrage (MEE,DEE) est modélisé, est un vecteur CD [n] [n] comportant l'ensemble [CD[1], ...,CD[n]l des dimensions caractéristiques mesurées, pondéré par l'ensemble {S[1], ...S[n]l des surfaces respectives des zones concentriques à symétrie de révolution, et cdo est une constante caractéristique d'une dimension commune de la zone centrale (26) de chaque motif de calibrage (MEN) lorsque cette zone centrale (26) est exposée de manière isolée au faisceau d'électrons en l'absence de zones concentriques à symétrie de révolution autour de cette zone centrale (26). 6. Procédé d'estimation de motifs (MPF,DpF) à imprimer selon l'une quelconque des revendications 1 à 5, dans lequel la seconde portion de fonction d'étalement du point (PSFFE) est modélisée à partir d'une fonction de type gaussienne, les paramètres de cette fonction gaussienne étant déterminés par simulation sur une pluralité de motifs de calibrage. 7. Procédé d'estimation de motifs (MPF,DPF) à imprimer selon l'une quelconque des revendications 1 à 6, dans lequel les zones concentriques (11[k]) à symétrie de révolution des motifs de calibrage sont toutes de même surface. 8. Procédé d'estimation de motifs (MPF,DPF) à imprimer selon les revendications 5 et 7, dans lequel la première portion de fonction d'étalement du point (PSFBE) caractérisant les électrons rétrodiffusés par le substrat est définie par l'expression suivante : 1 0 0 ([CDW 1 PSFBE oc H - 0 cdo . 9:1)[CD[n]1 1 9. Procédé d'estimation de motifs (M;,.F,D;,.F) à imprimer selon l'une quelconque des revendications 1 à 8, dans lequel la zone centrale (26) de chaque motif de calibrage (MEN) comporte l'une des zones concentriques à symétrie de révolution, celle-ci étant alors placée au centre dudit motif. 10. Dispositif d'impression (10) de motifs (MPF,DPF) sur plaque ou sur masque par lithographie à faisceau d'électrons, comportant : - des moyens d'impression (14) par lithographie à faisceau d'électrons, dans une résine disposée sur un substrat, d'un ensemble de motifs de calibrage .-CF (11/1 D 1 ,- CF. I - des moyens de mesure (15) de dimensions caractéristiques (CD) de l'ensemble de motifs de calibrage (M DCF. 1 après impression et des ,.- -CF,- moyens de stockage de ces dimensions caractéristiques (CD) en mémoire (16), - un processeur (12) ayant accès à la mémoire (16) et programmé pour fournir une estimation de la fonction d'étalement du point (PSF) caractérisant des électrons diffusés lors de l'impression à partir des dimensions caractéristiques (CD) mesurées, - des moyens d'estimation (24) des motifs (MPF,DPF) à imprimer par convolution de la fonction d'étalement du point (PSF) fournie avec une valeur initiale (MpF,DpF) de ces motifs à imprimer,caractérisé en ce que, chaque motif de calibrage (Mc [d) destiné à être imprimé dans la résine comportant une zone centrale (26) exposée au faisceau d'électrons et une pluralité de zones concentriques à symétrie de révolution disposées autour de cette zone centrale (26) : les moyens de mesure (15) sont des moyens de mesure des dimensions caractéristiques (CD) des zones centrales (26) des motifs après impression, et le processeur (12) est programmé pour calculer l'estimation de la fonction d'étalement du point (PSF) en tant que somme d'une première portion de fonction d'étalement du point (PSFBE) caractérisant des électrons rétrodiffusés par le substrat et d'une seconde portion de fonction d'étalement du point (PSFFE) caractérisant des électrons diffusés dans la résine, et par inversion d'une modélisation analytique de l'effet, sur lesdites dimensions caractéristiques (CD), de l'application de la première portion de fonction d'étalement du point (PSF) à l'ensemble de motifs de calibrage (MCF,DCF) -