FR 3040544 A1 20170303 FR 1558027 A 20150828 PROCEDE DE CONTROLE DE DEPLACEMENT D'UN VORTEX D'ABRIKOSOV L'invention a pour objet un procédé de contrôle de déplacement d'un vortex d'Abrikosov d'un élément constitué d'un matériau supraconducteur, dit élément supraconducteur. Il est connu qu'un élément supraconducteur de type II laisse pénétrer un champ magnétique de valeur comprise entre un premier champ critique et un deuxième champ critique sous la forme de vortex d'Abrikosov. o Un vortex d'Abrikosov est un tube de flux magnétique autour duquel circule un courant. Un élément supraconducteur présente des propriétés physiques spécifiques, particulièrement lorsqu'il est compris dans une jonction Josephson. 15 Une jonction Josephson est réalisée par exemple par un tricouches d'empilement de deux couches minces de l'élément supraconducteur et séparées par un lien faible. Le lien faible est constitué notamment par une couche mince d'un matériau isolant ou d'un métal. 20 Dans une telle jonction Josephson, il existe un courant sans qu'aucune différence de potentiel ne soit appliquée, dû à un passage de paires de Cooper de spins des électrons supraconducteurs par effet tunnel à travers la barrière que constitue le lien faible. Le courant maximal qui peut parcourir la jonction est appelé courant 25 critique. Les jonctions Josephson ont permis des avancées technologiques importantes, comme par exemple la mise au point de détecteurs de champ magnétique de type SQUID. 3040544 2 Une autre application importante des jonctions Josephson est le stockage d'information et permet de surcroit de générer un bit quantique, ou qubit. Il a été démontré que le courant critique d'une jonction Josephson est 5 fortement influencé par la présence de vortex dans son voisinage proche, et peut en particulier être réduit fortement selon la distribution des vortex dans les couches supraconductrices de la jonction Josephson. Pour surmonter cette difficulté, il est connu de modifier la distribution globale des vortex dans les couches supraconductrices, par exemple via o une application contrôlée de champ magnétique. Néanmoins, cette approche globale ne permet pas de contrôler un déplacement d'un vortex unique. Récemment, il a été démontré la possibilité de déplacer un vortex unique par microscopie à force magnétique. 15 Toutefois, cette méthode présente l'inconvénient de ne permettre qu'un déplacement lent du vortex et, de surcroit, à très faible échelle seulement. Un autre inconvenient réside dans le fait qu'une force appliquée par le microscope sur le vortex s'avère souvent inférieure à une force de piégeage du vortex dans l'élément supraconducteur, et qu'alors, le déplacement du 20 vortex par cette méthode est irréalisable. Le but de l'invention est de remédier au moins partiellement à ces inconvenients. A cet effet, l'invention a pour objet un procédé de contrôle de déplacement d'un vortex d'Abrikosov d'un élément constitué d'un matériau 25 supraconducteur, dit élément supraconducteur, comprenant une étape de chauffage inhomogène de l'élément supraconducteur muni du vortex d'Abrikosov dans une zone de chauffage de façon à déplacer le vortex d'Abrikosov vers la zone de chauffage. 3040544 3 Grâce à la présente invention, du fait que le chauffage est inhomogène, un gradient de température est généré et le vortex est attiré vers la zone la plus chaude, c'est-à-dire la zone de chauffage, vers laquelle il se déplace. Ainsi, il suffit de générer la zone chaude vers l'endroit où l'on souhaite 5 positionner le vortex, ce qui assure le déplacement contrôlé du vortex. Selon une autre caractéristique de l'invention, au cours de l'étape de chauffage inhomogène de l'élément supraconducteur, on chauffe la zone de chauffage par application d'un faisceau laser. Selon une autre caractéristique de l'invention, le procédé comprend une étape de focalisation du faisceau laser dans la zone de chauffage. Selon une autre caractéristique de l'invention, au cours de l'étape de chauffage, le faisceau laser génère une force d'attraction du vortex d'une valeur : cDo alici fT = - 47r * -aT * IVT I, où (rio est un quantum de flux magnétique, Hm est un premier champ critique de l'élément supraconducteur, et T est une température associée au faisceau laser. Selon une autre caractéristique de l'invention, le vortex d'Abrikosov se déplace vers la zone de chauffage si la force d'attraction générée par le faisceau laser est supérieure à une force de piégeage du vortex d'Abrikosov dans l'élément supraconducteur. 25 Selon une autre caractéristique de l'invention, la force de piégeage est 3040544 4 égale à : Fp = j,(1)0, où (rio est un quantum de flux magnétique et jc. est un courant critique de l'élément supraconducteur. 5 Selon une autre caractéristique de l'invention, une puissance du laser est comprise entre 1 pW et 1000 pW. Selon une autre caractéristique de l'invention, une longueur d'onde du laser est comprise dans le spectre électromagnétique visible. 10 Selon une autre caractéristique de l'invention, le procédé comprend une étape de détermination de la valeur de la force de piégeage par une force minimale à appliquer par le faisceau laser pour que le vortex d'Abrikosov se déplace. 15 L'invention a également pour objet une utilisation du procédé de contrôle tel que décrit précédemment pour former d'une distribution de vortex dans un élément supraconducteur, dans laquelle on applique le procédé tel que décrit précédemment à l'un au moins des vortex d'une distribution initiale de 20 vortex. De préférence, on applique le procédé tel que décrit précédemment successivement à chaque vortex de la distribution de vortex. Avantageusement, on prévoit une étape d'évitement de collision entre le vortex déplacé par le procédé tel que décrit précédemment et les autres vortex de la distribution de vortex. L'invention a également pour objet un procédé de contrôle tel que décrit précédemment à une jonction Josephson munie d'un vortex d'Abrikosov pour contrôler une valeur d'un courant critique de ladite jonction Josephson. 3040544 5 L'invention a également pour objet une utilisation du procédé tel que décrit précédemment pour un stockage d'information dans laquelle une écriture d'information est réalisée par le déplacement du vortex d'Abrikosov. L'invention a également pour objet un programme d'ordinateur comprenant une suite d'instructions mettant en oeuvre les étapes du procédé tel que décrit précédemment lorsqu'elles sont exécutées par un processeur. L'invention a également pour objet un support de stockage non transitoire lisible par ordinateur comprenant un programme d'ordinateur tel que décrit précédemment. 15 D'autres caractéristiques et avantages de l'invention apparaîtront encore à la lecture de la description qui va suivre. Celle-ci est purement illustrative et doit être lue en regard des dessins annexés sur lesquels : - la figure 1 illustre une représentation schématique d'un système de mise en oeuvre d'un procédé de contrôle de déplacement d'un 20 vortex d'Abrikosov selon la présente invention ; - la figure 2 illustre une vue en coupe transversale d'un échantillon d'un élément supraconducteur pour le système de la figure 1 ; - la figure 3 illustre une vue en perspective d'un déplacement d'un vortex d'Abrikosov par un faisceau laser du système de la figure 1 ; 25 - les figures 4a, 4b et 4c illustrent une séquence expérimentale de déplacement de chaque vortex d'un ensemble de vortex par le système de la figure 1 ; 5 10 3040544 6 - la figure 5 illustre un histogramme d'un rapport d'un nombre (Nd) de vortex déplacés par le faisceau laser de la figure 1 sur l'ensemble de neuf vortex (N) des figures 4a, 4b et 4c en fonction d'une puissance du faisceau laser à trois températures différentes, 5 T = 7,5 K, T = 6,5 K et T = 4,7 K ; - la figure 6a illustre une distribution initiale de vortex dans un échantillon ; et - la figure 6b illustre une distribution de vortex après déplacement par le procédé selon la présente invention. o Système pour la mise en oeuvre d'un procédé de contrôle de déplacement d'un vortex d'Abrikosov Comme visible sur la figure 1, un système 1 pour la mise en oeuvre d'un procédé de contrôle de déplacement 2 d'un vortex d'Abrikosov selon la présente invention comprend une source thermique 3 appliquée sur un 15 échantillon 4 d'un élément supraconducteur. La source thermique 3 est configurée pour chauffer de manière inhomogène l'échantillon 4, comme il sera détaillé ultérieurement. Sur le mode de réalisation illustré, la source thermique 3 de chauffage inhomogène est un laser 5 apte à émettre un faisceau laser F. 20 Comme visible sur la figure 1, le système 1 comprend des moyens de mise en forme 6 du faisceau laser comprenant une lentille de collimation 7 du laser 5 et une optique télécentrique 8. Les moyens de mise en forme comprennent également un diviseur de faisceau 9. 25 Les moyens de mise en forme 6 permettent de focaliser le faisceau laser sur une zone déterminée de l'échantillon 4, dite zone de chauffage 10, comme il va être expliqué. Le système 1 comprend également des moyens d'imagerie 11 pour 3040544 7 visualiser un suivi des déplacements du vortex d'Abrikosov. Comme visible sur la figure 1, les moyens d'imagerie 11, de préférence de type magnéto-optique, comprennent une source de lumière incidente 12, par exemple une source de lumière blanche, dirigée sur une face externe 5 de l'échantillon 4. La lumière incidente 13 émise par la source de lumière incidente 12 présente une polarisation rectiligne P selon une direction appelée X. Les moyens d'imagerie 11 comprennent également un élément de visualisation 14, tel qu'une caméra CCD, orientée pour recevoir un faisceau io lumineux 15 émis par la source 12 et réfléchi par l'échantillon 4, comme il va être expliqué. Le système 1 comprend également un dispositif de refroidissement de l'échantillon 4, non illustré. Le système 1 comprend également un dispositif d'application d'un champ 15 magnétique à l'échantillon 4, non illustré. Bien entendu, l'invention ne se limite pas au système de mise en oeuvre 1 illustré et d'autres systèmes peuvent être utilisés pour le procédé 2 selon la présente invention. Par exemple, selon une variante, le système de mise en oeuvre se 20 présente sous la forme d'un fil comprenant un microcontact pour chauffer l'échantillon 4 à distance suffisante du vortex. Echantillon L'échantillon 4 comprend une couche mince, appelée par la suite film, de l'élément supraconducteur S. 25 L'élément supraconducteur S est de type II, c'est-à-dire que, lorsque la température appliquée à l'échantillon 4 par le dispositif de refroidissement est inférieure à une température de transition supraconductrice, notée Tc., un 3040544 8 champ magnétique appliqué par le dispositif d'application d'un champ magnétique pénètre sous la forme de vortex d'Abrikosov si le champ magnétique appliqué est compris entre un premier champ magnétique critique Hm et un second champ magnétique critique Hc.2 . 5 Comme déjà indiqué, un vortex d'Abrikosov est un tube de flux magnétique autour duquel circule un courant. Par la suite, on fait référence au vortex d'Abrikosov sous le terme vortex V. Sur le mode de réalisation illustré à la figure 2, le film supraconducteur o est un film de niobium (Nb) présentant une épaisseur ds. Avantageusement, le film de niobium est déposé sur un substrat, par exemple de SiO2/Si, non illustré. De préférence, le dépôt de niobium sur le substrat est réalisé par une technique de pulvérisation cathodique magnétron, connue de l'homme du 15 métier. Comme visible sur la figure 2, l'échantillon 4 comprend également une couche d'indication visuelle 16. La couche d'indication visuelle 16 est configurée pour modifier la polarisation de la lumière incidente 13, comme il va être expliqué, ce qui 20 permet le suivi visuel du déplacement du vortex V par les moyens d'imagerie 11. Sur la figure 2, la couche d'indicateur 16 est constituée à base d'un cristal de Bi:LuIG avec magnétisation planaire. La couche de Bi:LuIG forme une face externe 17 de l'échantillon 4. 25 La lumière incidente 13 des moyens d'imagerie magnéto-optique 11 est dirigée sur la couche de Bi:LuIG 16 où la polarisation du faisceau est modifiée du fait d'un effet Faraday d'interaction entre la source lumineuse et 3040544 9 un champ magnétique du vortex V. Avantageusement, la polarisation du faisceau tourne de sorte qu'un angle entre la polarisation du faisceau réfléchi 15 et la polarisation du faisceau incident 13 est un angle OF = 219Bzd, où 19 est une constante de 5 Verdet de la couche 16 de Bi:LuIG, c/ est une épaisseur de la couche 16 de Bi:LuIG et Bz est une composante selon un axe Z du champ magnétique du vortex V. Procédé de déplacement contrôlé de chaque vortex Le procédé de contrôle de déplacement du vortex unique V selon la o présente invention comprend une étape de chauffage inhomogène de l'élément supraconducteur S dans la zone de chauffage 10 de façon à déplacer le vortex V vers la zone de chauffage 10. Par unique, on entend que le procédé 2 permet de contrôler le déplacement individuel de chaque vortex d'une distribution de vortex. 15 Comme visible sur la figure 3, le faisceau laser F est dirigé sur l'échantillon, ce qui crée la zone de chauffage 10. Le laser est avantageusement focalisé sur la zone de chauffage 10. Au cours de l'étape de chauffage, le faisceau laser F génère une force d'attraction du vortex V d'une valeur : 20 fT = - cDo -47r * -aT * IVTI, où (rio est un quantum de flux magnétique, Hm est le premier champ critique de l'élément supraconducteur, T est une température associée au faisceau laser et VT est le gradient de température généré depuis le faisceau laser. On note que le gradient VT généré par le faisceau laser F dépend d'une 25 distance entre la zone de chauffage 10 et du vortex V et de la puissance laser P selon les lois de Fourier de la propagation de la chaleur dans l'échantillon 4. 3040544 10 La dérivée du premier champ critique 141 avec la température s'écrit : , (lnK + 0,5), où Ào est une longueur de pénétration = 47r4iTc magnétique dans le film supraconducteur, et K est un paramètre GinzburgLandu du film supraconducteur. 5 La force de piégeage s'exprime encore : 2 EP() ' fT = (ink+u,D) VT _. 4n2.0 Tc Le vortex V se déplace vers la zone de chauffage 10 si la force d'attraction générée par le faisceau laser F est supérieure à une force de piégeage du vortex d'Abrikosov dans l'élément supraconducteur par des o centres de piégeage. La force de piégeage est égale à : Fp = jc(1)0, où jc. est un courant critique de l'élément supraconducteur. La force minimale de déplacement de vortex /T'in est donc 'DT =(4R)2 iCÂÔ Tc min (lnic+0,5) (Do 15 On note que le courant critique dépend de l'épaisseur ds du film supraconducteur S, ce qui implique que la force minimale dépend de l'épaisseur du film supraconducteur. On note également que le courant critique, et par conséquent la force minimale, dépend également de la température appliquée à l'échantillon 4. 20 On note qu'une puissance maximale du faisceau laser est déterminée par une condition selon laquelle une taille du faisceau, en particulier un rayon du faisceau, est inférieure à une distance entre vortex, quand plusieurs vortex V sont générés dans le film supraconducteur. Résultats expérimentaux 25 Expérience 1 Une première expérience a été menée selon le procédé 2, mis en oeuvre par le système 1 sur l'échantillon 4. Hci aT 3040544 11 Selon cette première expérience, le film de Nb de l'échantillon 4 présente une épaisseur ds de l'ordre de 450 nm. La température est abaissée par le dispositif de refroidissement à une valeur de T = 4,7 K inférieure à la température de transition Tc = 8,4 K. 5 Un champ magnétique extérieur Hext est appliqué d'une valeur de Hext = 4,4 Oe. Comme visible à la figure 4, il résulte de ces conditions de température et champ magnétique une distribution initiale 20 de vortex V. De manière connue, les vortex V se positionnent dans l'échantillon 4 en o fonction des centres de piégeage du film supraconducteur S. Sur la figure 4a, neuf vortex V de la distribution initiale 20 sont illustrés. Selon cette expérience, le laser 5 présente une puissance P et un rayon ro de l'ordre de 0,4 pm. Le faisceau laser F est dirigé vers la face externe 12 de l'échantillon 4. 15 Le faisceau laser F s'étend de préférence transversalement à un plan de la face externe 12. Le faisceau laser F est focalisé sur la face externe 12, successivement à une distance r = 4 pm de chaque vortex V de la distribution 21 de vortex V. Le laser est d'une longueur d'onde visible, de l'ordre de 630 nm. 20 Le faisceau laser chauffe localement l'échantillon 4, ce qui crée la zone chaude 11. Le chauffage par laser dure pendant 1 seconde. Comme visible à la figure 4b, pour une puissance laser de 4,9 pW, deux vortex V se sont déplacés. 25 Comme visible sur la figure 4c, tous les vortex V se sont déplacés pour une puissance Pmin de l'ordre 7,3 pW. Expérience 2 Les conditions expérimentales sont inchangées par rapport à la première 30 expérience. En particulier, l'échantillon 4 est identique à celui de la première 3040544 12 expérience. Le laser est identique à celui de la première expérience. Le laser est focalisé à une distance r = 4 pm de chaque vortex V de la distribution 21. 5 L'échantillon 4 est refroidi à trois températures T, T = 7,5 K, T = 6,5 K et T = 4,7 K. Comme visible sur la figure 5, à T = 4,7 K, 100% des neuf vortex V se déplacent pour une puissance du laser supérieure ou égale à 7,2 pW. A T = 6,5 K, 100% des neuf vortex V se déplacent pour une puissance 10 du laser supérieure ou égale à 3,3 pW. A T = 7,5 K, 100% des neuf vortex V se déplacent pour une puissance du laser supérieure ou égale à 1,8 pW. On constate que plus la température est proche de la température de transition supraconductrice, moins la puissance minimale est élevée, ce qui 15 est dû au fait que la force de piégeage augmente avec une diminution de la température appliquée. Expérience 3 L'échantillon 4 est identique à celui de la première expérience. 20 Le laser est identique à celui de la deuxième expérience. Le laser est focalisé à une distance r = 1 pm de chaque vortex V de la distribution 20. Dans ces conditions, à T = 4,7 K, 100% des neuf vortex V se déplacent pour une puissance du laser supérieure ou égale à 6,5 pW. 25 Dans cette expérience, la puissance maximale mesurée est de l'ordre de 15 pW afin de ne pas dépasser la température de transition supraconductrice. Une modélisation basée sur les lois de Fourier de la propagation de la chaleur dans l'échantillon 4 permet de déterminer que le gradient de 30 température VT généré par le faisceau laser F à la puissance Pmin est de l'ordre de 2 K/pm. 3040544 13 Considérant que la longueur de pénétration magnétique dans le film supraconducteur Ào est égale à 90 nm, et le paramètre Ginzburg-Landau K est égal à 6, alors la force de piégeage est évaluée selon l'équation de l'ordre de 9 pN/ pm. 5 Une mesure indépendante du courant critique donne une valeur de jc. = 0,7 . 106 A/cm2, ce qui correspond à une force de piégeage de l'ordre de 14 pN/ pm. On constate que les deux procédés donnent une valeur du même ordre de grandeur pour la force de piégeage, ce qui assure que le procédé 2 o permette la détermination de la valeur du courant critique. Programme d'ordinateur L'invention a également pour objet un programme d'ordinateur comprenant une suite d'instructions mettant en oeuvre les étapes du procédé 2 lorsqu'elles sont exécutées par un processeur. 15 L'invention a également pour objet un support de stockage non transitoire lisible par ordinateur comprenant un programme d'ordinateur tel que décrit précédemment. Expérience 4 20 Sur la figure 6a est illustrée une autre distribution initiale 20 de vortex V, créée dans l'échantillon 4 sous un champ extérieur de l'ordre de 4,3 Oe. Comme visible sur la figure 6b, le programme d'ordinateur basé sur le procédé 2 permet un réarrangement de la distribution initiale 20 en une distribution 20' souhaitée, triangulaire. 25 De préférence, le procédé comprend également une étape d'évitement de collision entre le vortex déplacé par le procédé 2 et les autres vortex de la distribution de vortex 20. 3040544 14 Estimation numérique Cette estimation numérique est basée sur la loi de Fourier par utilisation du logiciel Comsol Multiphysics (marque déposée). Le film de Nb de l'échantillon 4 présente une épaisseur ds de l'ordre de 5 5 nm. Le laser 5 présente une puissance P = 530 pW. Le vortex se déplace sur une distance de 300 nm avec une constante de temps de 10 ps. Cette estimation démontre que le procédé 2 selon la présente invention o permet un déplacement de vortex à une fréquence de l'ordre de 100 GHz, ce qui est une valeur beaucoup plus élevée que les fréquences de déplacement connues de l'état de la technique. Une telle fréquence de déplacement assure l'utilisation du procédé 2 pour un stockage d'information. 15 Utilisations Le procédé 2 s'applique tout particulièrement à une jonction Josephson munie d'un vortex d'Abrikosov pour contrôler une valeur du courant critique de la jonction Josephson. 20 Comme déjà indiqué, une jonction Josephson un multicouches, de préférence tricouches d'empilement de deux couches minces de l'élément supraconducteur séparées par un lien faible. Le lien faible est constitué notamment par une couche mince d'un matériau isolant ou d'un métal. 25 Dans une telle jonction Josephson, il existe un courant sans qu'aucune différence de potentiel ne soit appliquée, dû à un passage de paires de Cooper de spins des électrons supraconducteurs par effet tunnel à travers la barrier que constitue le lien faible. On rappelle que le courant critique est le courant maximal pouvant 3040544 15 parcourir la jonction Josephson. Grâce au procédé 2, il est possible de changer la valeur du courant critique par le déplacement du vortex, du fait que le courant critique de la jonction Josephson est fortement influencé par la présence de vortex dans 5 son voisinage proche. Le procédé selon la présente invention permet de contrôler le déplacement de chaque vortex d'un ensemble de vortex de manière fiable et rapide et s'applique avantageusement aux utilisations précédemment décrites (distribution de vortex, stockage d'information, jonction Josephson). 10 REVENDICATIONS1. Procédé de contrôle de déplacement d'un vortex d'Abrikosov d'un élément constitué d'un matériau supraconducteur, dit élément supraconducteur, comprenant une étape de chauffage inhomogène de l'élément supraconducteur muni du vortex d'Abrikosov dans une zone de chauffage de façon à déplacer le vortex d'Abrikosov vers la zone de chauffage. 2. Procédé de contrôle selon la revendication 1, dans lequel, au cours de l'étape de chauffage inhomogène de l'élément supraconducteur, on chauffe la zone de chauffage par application d'un faisceau laser. 3. Procédé de contrôle selon la revendication 2, comprenant une étape de focalisation du faisceau laser dans la zone de chauffage. 4. Procédé de contrôle selon l'une des revendications 2 ou 3, dans lequel, au cours de l'étape de chauffage, le faisceau laser génère une force d'attraction du vortex d'une valeur : cl)o fT = _-47r * -ÔT * IVTI, où t est un quantum de flux magnétique, Hc1 est un premier champ critique de l'élément supraconducteur, et T est une température associée au faisceau laser. 5. Procédé de contrôle selon la revendication précédente, dans lequel le vortex d'Abrikosov se déplace vers la zone de chauffage si la force d'attraction générée par le faisceau laser est supérieure à une force de piégeage du vortex d'Abrikosov dans l'élément supraconducteur. 6. Procédé de contrôle selon la revendication précédente, dans lequel la force de piégeage est égale à : Fp = j,c1)0, où D est un quantum de flux magnétique et À est un courant critique de l'élément supraconducteur. 7. Procédé de contrôle selon l'une quelconque des revendications 2 à 6, dans lequel une puissance du laser est comprise entre 1 i_AN et 1000 pW. 3040544 17 8. Procédé de contrôle selon l'une quelconque des revendications 2 à 7, dans lequel une longueur d'onde du laser est comprise dans le spectre électromagnétique visible. 9. Procédé de contrôle selon l'une quelconque des revendications 5 à 8, 5 comprenant une étape de détermination de la valeur de la force de piégeage par une force minimale à appliquer par le faisceau laser pour que le vortex d'Abrikosov se déplace. 10. Utilisation du procédé de contrôle selon l'une des revendications précédentes pour former une distribution de vortex dans un élément 10 supraconducteur, dans laquelle on applique le procédé selon l'une quelconque des revendications précédentes à l'un au moins des vortex d'une distribution initiale de vortex. 11. Utilisation selon la revendication précédente, dans laquelle on applique le procédé selon l'une quelconque des revendications 1 à 9 15 successivement à chaque vortex de la distribution de vortex. 12. Utilisation selon l'une des revendications 10 ou 11, comprenant une étape d'évitement de collision entre le vortex déplacé par le procédé selon l'une des revendications 1 à 9 et les autres vortex de la distribution de vortex. 20 13. Utilisation du procédé de contrôle selon l'une des revendications 1 à 9 à une jonction Josephson munie d'un vortex d'Abrikosov pour contrôler une valeur d'un courant eiiitique de ladite jonction Josephson. 14. Utilisation du procédé de contrôle selon l'une des revendications 1 à 9 pour un stockage d'information dans laquelle une écriture d'information 25 est réalisée par le déplacement du vortex d'Abrikosov. 15. Programme d'ordinateur comprenant une suite d'instructions mettant en oeuvre les étagés du procédé selon l'une des revendications 1 à 9 lorsqu'elles sont exécutées par un processeur. 16. Support de stockage non transitoire lisible par ordinateur comprenant un programme d'ordinateur selon la revendication précédente.