FR 3114672 A1 20220401 FR 2009795 A 20200925 Système de stabilisation autonome d’états quantiques comprenant : - un oscillateur harmonique quantique, - un dispositif auxiliaire configuré pour être dans un état fondamental |g> ou dans un état excité |f>, le dispositif étant couplé à l’oscillateur de manière dispersive, de sorte qu’une différence de fréquence entre l’état fondamental |g> et l’état excité |f> dépend de manière linéaire du nombre de bosons, - un dissipateur, - un excitateur configuré pour générer des premier (405) et deuxième (406) peignes de fréquences, destinés à être envoyés vers au moins l’un parmi l’oscillateur harmonique quantique, le dispositif auxiliaire et le dissipateur caractérisé en ce qu’un premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences (405) est égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne (406). Figure pour l’abrégé : Fig. 3 Système de stabilisation autonome d’états quantiques ayant une parité prédéterminée pour la correction d’erreur DOMAINE DE L'INVENTION L’invention concerne la stabilisation autonome d’états quantiques présentant une parité prédéterminée et la correction d’erreurs dans le codage d’information quantique, en particulier dans le codage qui encode de l’information dans un système bosonique ayant une dimension, au sens de l’espace de Hilbert, plus grande que les deux états nécessaires pour encoder un bit d’information. ETAT DE LA TECHNIQUE La mise en œuvre d’un traitement d’information quantique nécessite de mettre en place des protocoles de correction d’erreur permettant de corriger les changements indésirables de l’état des qubits qui se produisent inévitablement en raison de leur couplage au monde extérieur. Dans les protocoles d’information quantique de type bosonique, l’information est encodée dans des états d’un oscillateur harmonique quantique. Dans l’ensemble de ce texte, on entend par oscillateur harmonique tout oscillateur dont la fréquence caractéristique ne change pas plus de 20% par quantum d’excitation. Les états de l’oscillateur peuvent être décomposés sur la base des états de Fock, chaque état de Fock correspondant à un nombre particulier de bosons compris dans l’oscillateur. La parité du nombre de bosons peut permettre de définir une parité pour les états de l’oscillateur harmonique quantique. Un état pair de l’oscillateur peut être défini lorsque dans la décomposition de cet état sur la base des états de Fock, les seuls coefficients non nuls sont les coefficients des états de Fock correspondant à un nombre pair de bosons. Des procédés de correction d’erreur quantique sont connus pour corriger certaines erreurs, notamment dans le cas de protocole où le codage d’information utilise uniquement des états de l’oscillateur harmonique quantique de parité particulière. Ces procédés de correction d’erreur sont fondés sur un maintien de la parité des états qui codent l’information, voir par exemple les documents OFEK Nissim et al. « Extending the lifetime of a quantum bit with error correction in superconducting circuits », Nature 536 pages 441 à 445; MICHAEL Marios H.et al. « New Class of Quantum Error-Correcting Codes for a Bosonic Mode », Physical Review X Vol.6 pages 031006; GERTLER Jeffrey M. et al. « Protecting a Bosonic Qubit with Autonomous Quantum Error Correction » .Par exemple si une erreur se produit et que la parité est changée, l’erreur peut être corrigée en ajoutant un boson dans l’oscillateur harmonique quantique. Ces procédés de correction d’erreur peuvent cependant être extrêmement difficiles à réaliser expérimentalement. C’est le cas quand ils nécessitent de mesurer continuellement la parité afin de détecter en temps réel les changements de parité lorsqu’ils se produisent et appliquer immédiatement les corrections, comme par exemple dans les documents précédemment cités de OFEK Nissim et GERTLER Jeffrey M.. Des procédés de correction autonome, plus simples à mettre en œuvre, sont connus, voir par exemple le document précédemment cité de GERTLER Jeffrey M.. Ces procédés reposent sur une ingénierie du système permettant que la détection et la correction des erreurs soient réalisées de façon autonome, sans intervention de l’expérimentateur, à l’intérieur même du système quantique. Les procédés de correction autonome de l’art antérieur ne donnent cependant pas entière satisfaction car dans les cas où aucune erreur ne s’est produite, une correction peut tout de même être apportée et un boson peut être ajouté dans l’oscillateur harmonique. Ce boson ajouté induit un changement de la parité alors qu’elle avait la valeur prédéterminée, ce qui provoque une erreur. Il existe donc un besoin d’un protocole de correction d’erreur plus sélectif sur la parité, c’est-à-dire un protocole de correction d’erreur qui n’apporte une correction que lorsque la parité n’a pas une valeur prédéterminée. Un but de l’invention est de proposer un protocole de correction d’erreur plus sélectif sur la parité que dans l’art antérieur. Le but est atteint dans le cadre de la présente invention grâce à un système de stabilisation autonome d’états quantiques comprenant : - un oscillateur harmonique quantique configuré pour porter une information encodée dans des états de l’oscillateur ayant une parité prédéterminée d’un nombre de bosons compris dans l’oscillateur harmonique, - un dispositif auxiliaire configuré pour être dans un état fondamental |g> ou dans un état excité |f> parmi une pluralité d’états excités, le dispositif étant couplé à l’oscillateur de manière dispersive, de sorte qu’une différence de fréquence entre l’état fondamental |g> et l’état excité |f> dépend de manière linéaire d’un décalage de fréquence χ et du nombre de bosons, - un dissipateur configuré pour passer d’un état excité |1> dans un état fondamental |0>, - un excitateur configuré pour générer des premier et deuxième peignes de fréquences, destinés à être envoyés vers au moins l’un parmi l’oscillateur harmonique quantique, le dispositif auxiliaire et le dissipateur de sorte à ajouter un boson dans l’oscillateur harmonique si le nombre de bosons n’a pas la parité prédéterminée, un premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences étant égal au double du décalage de fréquence χ et étant égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne. Un tel système est avantageusement complété par les différentes caractéristiques suivantes prises seules ou en combinaison : - les premier et deuxième peignes présentent des pics dans le domaine temporel, les pics du premier peigne étant décalés temporellement par rapport aux pics du deuxième peigne ; - le système est configuré pour être, pour chaque entier k d’un ensemble d’entiers, dans : - un premier niveau défini par un premier nombre de bosons égal à l’entier k, l’état fondamental |g> et l’état fondamental |0>, - un deuxième niveau défini par un deuxième nombre de bosons égal à l’entier k+1, l’état excité |g> et l’état fondamental |0>, - un troisième niveau dans lequel le dispositif auxiliaire est dans l’état excité |f> et le dissipateur est dans l’état stable |0>, - un quatrième niveau dans lequel le dispositif auxiliaire est dans l’état fondamental |g> et le dissipateur est dans l’état excité |1>, l’excitateur étant configuré pour coupler, pour chaque entier k, les premier, troisième et quatrième niveaux si l’entier k n’a pas la parité prédéterminée, de sorte à déplacer le système du premier niveau dans le deuxième niveau, l’excitateur étant configuré pour générer les premier et deuxième peignes de fréquences de sorte que, pour chaque entier k : - le premier peigne comprend une première raie si l’entier k n’a pas la parité prédéterminée, la première raie ayant une fréquence égale à une fréquence de transition entre les premier et troisième niveaux, et - le deuxième peigne comprend une deuxième raie ayant une fréquence égale à une fréquence de transition entre les troisième et quatrième niveaux ; - l’oscillateur harmonique comprend le deuxième nombre de bosons égal à l’entier k+1 dans le quatrième niveau ; - l’oscillateur harmonique quantique est un premier résonateur micro-onde supraconducteur, le système auxiliaire est un transmon, le dissipateur est un deuxième résonateur micro-onde supraconducteur ayant un facteur de qualité inférieur à un facteur de qualité du premier résonateur micro-onde supraconducteur ; - une amplitude des raies du premier peigne de fréquences et du deuxième peigne de fréquences est caractérisée par une fréquence de couplage, la fréquence de couplage étant inférieure au quart du deuxième écart fréquentiel ; - l’information est encodée dans des états de l’oscillateur correspondant à un nombre de bosons compris entre un entier plancher ayant la parité prédéterminée et un entier plafond ayant la parité prédéterminée, l’ensemble d’entiers étant égal aux entiers compris entre les entiers et . L’invention porte également sur un procédé de stabilisation autonome d’états quantiques comprenant les étapes suivantes : - encoder une information dans des états d’un oscillateur harmonique quantique, les états ayant une parité prédéterminée d’un nombre de bosons compris dans l’oscillateur harmonique, - coupler de manière dispersive l’oscillateur harmonique à un dispositif auxiliaire, de sorte qu’une différence de fréquence entre un état fondamental |g> du dispositif auxiliaire et un état excité |f> du dispositif auxiliaire dépend de manière linéaire d’un décalage de fréquence χ et du nombre de bosons, - coupler le dispositif auxiliaire à un dissipateur configuré pour passer d’un état excité |1> dans un état fondamental |0> - générer des premier et deuxième peignes de fréquences, un premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences étant égal au double du décalage de fréquence χ et égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne et - envoyer les premier et deuxième peignes de fréquences vers au moins l’un parmi l’oscillateur harmonique quantique, le dispositif auxiliaire et le dissipateur de sorte à ajouter un boson dans l’oscillateur harmonique si le nombre de bosons n’a pas la parité prédéterminée. Un tel procédé est avantageusement complété par les différentes caractéristiques ou étapes suivantes prises seules ou en combinaison : - une optimisation de l’amplitude des raies des premier et deuxième peignes de fréquences, l’amplitude des raies étant caractérisée par une fréquence de couplage, comprenant les étapes suivantes : - (E1) choisir une valeur de fréquence de couplage inférieure au quart du deuxième écart fréquentiel, - (E2) mesurer un temps de vie des états ayant un nombre de bosons de la parité prédéterminée et - (E3) modifier la valeur de la fréquence de couplage en fonction du temps de vie mesuré selon un algorithme d’optimisation du temps de vie de sorte à identifier une valeur optimale de la fréquence de couplage correspondant à un maximum du temps de vie ; - une optimisation d’un taux de désexcitation égal à la probabilité de passer du niveau vers le niveau par unité de temps, comprenant les étapes suivantes : - (S1) choisir une valeur du taux de désexcitation inférieure au deuxième écart fréquentiel, - (S2) mesurer un temps de vie des états ayant un nombre de bosons de la parité prédéterminée et - (S3) modifier la valeur du taux de désexcitation en fonction du temps de vie mesuré selon un algorithme d’optimisation du temps de vie de sorte à identifier une valeur optimale du taux de désexcitation correspondant à un maximum du temps de vie. DESCRIPTION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit être lue en regard des dessins annexés sur lesquels : la figure 1 est une représentation schématique des états d’un oscillateur harmonique. la figure 2 est une représentation schématique d’un système de stabilisation d’états quantiques. la figure 3 est une représentation schématique d’une architecture de niveaux d’énergie d’un système de stabilisation d’états quantiques. la figure 4 est une représentation schématique de fonctions caractérisant l’évolution temporelle de peignes de fréquences générés par le système de stabilisation d’états quantiques. DESCRIPTION DETAILLEE DE L'INVENTION La présente invention a pour but de stabiliser une parité prédéterminée du nombre de bosons d’un oscillateur harmonique quantique de manière autonome. Par souci de simplicité, il est supposé dans la suite que c’est la parité paire qu’il faut stabiliser. Un nombre impair de bosons correspond donc à une erreur qu’il convient de corriger. Des changements triviaux permettraient de stabiliser la parité impaire. Canaux d’erreur et de correction Le principal canal d’erreur conduisant à un changement de parité est la perte d’un boson de l’oscillateur harmonique quantique. La figure 1 représente schématiquement les états d’un oscillateur harmonique et le canal d’erreur. Un oscillateur harmonique peut être dans un état |k>, où k est un entier et représente le nombre de bosons k contenus dans l’oscillateur. Les états |0>, |1>, etc., de l’oscillateur harmonique quantique sont représentés par des traits horizontaux sur la figure 1. Deux états successifs de l’oscillateur sont séparés d’une différence d’énergie constante où est la fréquence de résonance de l’oscillateur harmonique quantique, l’état |0> étant l’état de plus basse énergie. Le niveau d’énergie d’un état est représenté verticalement, de sorte que les traits horizontaux des états |0>, |1>, etc., sont séparés verticalement de la même distance sur la figure 1. Le processus de relaxation 100 est une transition d’un état |k> vers un état |k-1> et correspond à la perte d’un boson. Le processus de relaxation 100 est caractérisé par un taux de relaxation , de sorte que la probabilité de transition d’un état |k> vers un état |k-1> par unité de temps est égale à . Une manière de stabiliser la parité contre la perte d’un boson de manière autonome consiste à réaliser un canal de dissipation effectif 101 ajoutant un boson dans l’oscillateur harmonique quantique si le nombre de bosons est impair. Dans la figure 1, ce canal est représenté par des flèches continues entre un état |k> vers un état |k+1> lorsque k est impair. Dans le cas de la figure 1 et en accord avec le choix de stabiliser la parité paire, il est choisi de stabiliser des nombres k de bosons pairs, |0>, |2>, |4>, etc. Le canal de dissipation effectif 101 ou processus de correction 101 est caractérisé par un taux de correction qui est égal à la probabilité d’ajout d’un boson par unité de temps. Selon la manière dont le processus de correction 101 est réalisé, le taux de correction peut dépendre du nombre de bosons k. Le processus de stabilisation de la parité est d’autant plus efficace que le taux de correction , est grand devant le taux de relaxation . Cependant, dans l’art antérieur, les réalisations d’un canal de dissipation effectif 101 induisent un canal d’altération 102, comme par exemple dans le document précédemment cité GERTLER Jeffrey M.. Le canal d’altération 102 a comme effet d’ajouter un boson dans l’oscillateur alors que la parité est paire ce qui est indésirable. Une parité paire du nombre de bosons correspond en effet à un état de l’oscillateur qu’il convient de ne pas modifier. Dans la figure 1, le canal d’altération 102 est représenté par des flèches en pointillés depuis un état |k> vers un état |k+1> lorsque k est pair. Le canal d’altération 102 ou processus d’altération 102 est caractérisé par un taux d’altération , qui est égal à la probabilité d’ajout d’un boson par unité de temps. Selon la manière dont le processus de correction 101 est réalisé, le taux d’altération peut dépendre du nombre de bosons k. Il est important de minimiser le taux d’altération pour au moins les trois raisons suivantes. Le canal d’altération 102 engendre des erreurs de parité. Dans les réalisations d’un canal de dissipation effectif 101 de l’art antérieur, le taux d’altération croît lorsque le taux de correction croît. Cela limite la capacité à corriger le processus de relaxation 100, car augmenter le taux de correction augmente également le taux d’altération . Dans certains cas comme dans le code de chat à quatre pattes voir MIRRAHIMI NJP et al. “ Dynamically protected cat-qubits: a new paradigm for universal quantum computation ”, New Journal of Physics, vol. 16 et ci-après, les erreurs de parité engendrées par le processus d’altération 102 ne sont pas corrigées par le processus de correction 101 mais au contraire aggravées. Code de chats à quatre pattes Un chat à quatre pattes est défini comme la superposition de quatre états cohérents de l’oscillateur harmonique quantique dont les phases diffèrent d’une valeur de π/2. Lorsqu’un chat à quatre pattes est décomposé sur les états |k>, les seuls coefficients non nuls correspondent à des nombres p de bosons congrus à n modulo 4, où n est un entier prédéterminé qui peut valoir 0, 1, 2 ou 3. Chaque valeur de l’entier prédéterminé n définit un chat à quatre pattes. Il exste donc quatre chats à quatre pattes. Les chats à quatre pattes correspondant à l’entier prédéterminé n égal à 0 ou à 2 peuvent permettre de définir un premier qubit logique. Les chats à quatre pattes correspondant à l’entier prédéterminé n égal à 1 ou à 3 définissent un second qubit logique. Un changement de parité (causé par exemple par une perte ou un gain d’un boson dans l’oscillateur harmonique) dans un qubit logique est équivalent à un changement du qubit logique. Autrement dit la perte ou le gain d’un boson dans l’oscillateur harmonique fait passer du premier qubit logique au second qubit logique et inversement. Un code de chat à quatre pattes consiste à corriger une telle erreur. La perte d’un boson peut être corrigée en ajoutant immédiatement un boson dans l’oscillateur harmonique. Le gain d’un boson peut être corrigé en retirant immédiatement un boson de l’oscillateur harmonique. La perte ou l’ajout consécutif de deux bosons de l’oscillateur harmonique ne peut pas, en revanche, être corrigé et provoque une erreur logique. Pour qu’une perte successive de deux bosons se produise, il est nécessaire de perdre un premier boson par le processus de relaxation 100, puis de perdre un deuxième boson par le processus de relaxation 100 avant que le gain d’un boson par processus de correction 101 ne se produise. La probabilité que le deuxième boson soit perdu avant que le processus de correction agisse vaut le rapport du taux de relaxation sur le taux de correction , dans la limite où est grand devant . C’est dans cette limite que le procédé de stabilisation de la parité prédéterminée est efficace. Pour qu’un gain successif de deux bosons se produise, il est nécessaire de gagner un premier boson par le processus d’altération 102, puis par exemple de gagner un deuxième boson par le processus de correction 101 avant que la perte d’un boson par le processus de relaxation 100 ne se produise. La probabilité que le deuxième boson soit gagné avant que le processus de relaxation agisse vaut , dans la limite où est grand devant . Le processus de correction 101 ne corrige pas une erreur causée par le processus d’altération 102, mais conduit au contraire avec une forte probabilité à l’ajout d’un deuxième boson, ce qui engendre une erreur logique. Sélectivité La présente invention permet d’augmenter significativement le rapport du taux d’altération sur le taux de correction. Ce rapport peut être nommé sélectivité du protocole de correction d’erreur. La sélectivité donne le rapport de la probabilité d’ajouter un boson lorsque c’est souhaitable c’est-à-dire lorsque le nombre de bosons est impair sur la probabilité d’ajouter un boson lorsque ce n’est pas souhaitable c’est-à-dire lorsque le nombre de bosons est pair. Plus la sélectivité est importante et plus le protocole corrige les erreurs dues au canal de relaxation 100 en ajoutant un boson lorsque le nombre de bosons est impair, sans ajouter un boson alors que le nombre de bosons est pair. Système de stabilisation d’états quantiques Un système 200 de stabilisation d’états quantiques configuré pour mettre en œuvre un protocole de correction d’erreur quantique est représenté schématiquement en figure 2. Le système 200 comprend un oscillateur harmonique quantique 202, un dispositif auxiliaire non-linéaire 201 qui est couplé de manière dispersive à l’oscillateur harmonique quantique 202. L’oscillateur harmonique quantique a une fréquence de résonance . Son hamiltonien peut être représenté comme suit : Où est la constante de Planck réduite et et sont les opérateurs de création et d’annihilation de l’oscillateur harmonique quantique vérifiant L’oscillateur harmonique sert à encoder de l’information quantique. Le dispositif auxiliaire 201 peut être dans un niveau stable , et au moins un niveau excité . La fréquence de transition entre les états et lorsque l’oscillateur harmonique quantique est dans son état fondamental est notée . La différence d’énergie ou, de manière équivalente, la différence de fréquence entre l’état ou le niveau fondamental et l’état ou le niveau excité dépend de manière linéaire du nombre k de bosons dans l’oscillateur harmonique quantique (voir BLAIS Alexandre et al. “ Cavity quantum electrodynamics for superconducting electrical circuits: An architecture for quantum computation ”, Physical Review A vol. 69, page 062320). Une telle différence d’énergie peut être caractérisée par un décalage de fréquence sous la forme . La différence de fréquence entre l’état fondamental |g> et l’état excité |f> dépend de manière linéaire du décalage de fréquence χ et du nombre de bosons dans l’oscillateur. Le décalage de fréquence est égal à la différence entre la différence d’énergie entre le niveau excité et le niveau stable pour un nombre k+1 de bosons et la différence d’énergie entre le niveau excité et le niveau stable pour un nombre k de bosons. Le décalage de fréquence peut être positif ou négatif. Le dispositif auxiliaire 201 peut permettre une action qui dépend de la parité du nombre de bosons. Le hamiltonien dispersif peut être représenté comme suit : Où et sont les opérateurs de création et d’annihilation de l’oscillateur harmonique quantique vérifiant Le système 200 comprend aussi un dissipateur 203, comportant au moins deux niveaux ou états, dont l’un, le niveau ou état , est stable et l’autre, le niveau ou état est un état excité relaxe spontanément vers . La différence d’énergie entre ces deux états correspond à une fréquence de transition entre les états et . La relaxation (ou désexcitation) spontanée du niveau vers le niveau se produit avec un taux de désexcitation , qui est égal à la probabilité de passer du niveau vers le niveau par unité de temps. La relaxation spontanée du niveau vers le niveau peut permettre de rendre, au cours du protocole de correction, des déplacements de population irréversibles. Chaque niveau d’énergie du système 200 est noté et dépend : - du niveau d’énergie de l’oscillateur harmonique quantique 202, désignant le nombre de bosons de l’oscillateur harmonique quantique, - du niveau d’énergie du dispositif auxiliaire 201, désignant le niveau stable ou le niveau excité et - du niveau d’énergie du dissipateur 203, désignant le niveau stable ou le niveau . Le système 200 comprend également un excitateur 204 configuré pour créer deux peignes de fréquences 205 et 206. Il n’est pas nécessaire de prendre en compte l’état de l’excitateur 204 pour caractériser l’état quantique du système 200 qui permet de porter l’information dans le protocole. Les peignes de fréquences 205 et 206 générés par l’excitateur 204 sont envoyés vers au moins l’un de l’oscillateur harmonique quantique 202, du dispositif auxiliaire non-linéaire 201 et du dissipateur 203. Les peignes de fréquence sont en résonance avec certaines transitions des niveaux d’énergie du système 200 de sorte à ajouter un boson dans l’oscillateur harmonique si l’entier k n’a pas la parité prédéterminée. L’écart fréquentiel séparant deux raies successives du premier peigne de fréquences est égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne. De cette manière, il est possible de séparer temporellement l’action du premier peigne et l’action du deuxième peigne. En particulier, cela permet un contrôle plus fin des transferts de population entre les différents niveaux d’énergie du système 200. Cela diminue le canal d’altération 102. Le canal d’altération est même potentiellement complètement supprimé en l’absence d’autres effets parasites comme, par exemple, les excitations thermiques du dispositif auxiliaire. Architecture de niveaux et peignes de fréquences Le système 200 peut être configuré pour être, pour chaque entier k d’un ensemble d’entiers, dans : - un premier niveau défini par un premier nombre de bosons égal à l’entier k, l’état fondamental |g> et l’état fondamental |0>, - un deuxième niveau défini par un deuxième nombre de bosons égal à l’entier k+1, l’état excité |g> et l’état fondamental |0>, - un troisième niveau dans lequel le dispositif auxiliaire est dans l’état excité |f> et le dissipateur est dans l’état stable |0>, - un quatrième niveau dans lequel le dispositif auxiliaire est dans l’état fondamental |g> et le dissipateur est dans l’état excité |1>. Pour chaque entier k de l’ensemble d’entiers, le troisième niveau est défini par un nombre de bosons égal à l’entier k+p, où p est un entier positif ou négatif. Pour chaque entier k de l’ensemble d’entiers, le quatrième niveau est défini par un nombre de bosons égal à l’entier k+q, où q est un entier positif ou négatif. L’excitateur est configuré pour coupler, pour chaque entier k, les premier, troisième et quatrième niveaux si l’entier k n’a pas la parité prédéterminée, de sorte à déplacer le système du premier niveau dans le deuxième niveau. L’excitateur est configuré pour générer des premier et deuxième peignes de fréquences de sorte que, pour chaque entier k : - le premier peigne comprend une première raie si l’entier k n’a pas la parité prédéterminée, la première raie ayant une fréquence égale à une fréquence de transition entre les premier et troisième niveaux, et - le deuxième peigne comprend une deuxième raie ayant une fréquence égale à une fréquence de transition entre les troisième et quatrième niveaux. L’ensemble d’entiers peut notamment être un intervalle d’entiers qui couvre le nombre moyen de bosons initialement généré dans l’oscillateur harmonique. L’ensemble d’entiers peut également ne comprendre qu’un seul entier. La figure 3 représente schématiquement un exemple d’une telle architecture de niveaux d’énergie du système 200 pour un procédé de stabilisation d’états quantiques. Le niveau d’énergie est représenté verticalement et croissant vers le haut. La figure 3 représente - le premier niveau ou état référencé 401, qui est un état que l’on souhaite conserver si k est pair et corriger en lui ajoutant un boson si k est impair, - le deuxième niveau ou état , référencé 402, - une pluralité de troisièmes niveaux ou états référencée 403, avec p l’entier positif ou négatif; l’entier p est supposé positif sur la figure 3 ; les niveaux sont séparés verticalement selon la valeur de k en raison du couplage dispersif entre l’oscillateur harmonique 202 et le dispositif auxiliaire 201 ; le décalage de fréquence est supposé négatif sur cette figure, et sépare verticalement les états correspondant à deux valeurs successives de k, - un quatrième niveau ou état référencé 404 c’est-à-dire un quatrième niveau pour lequel l’entier q a été choisi égal à 1. L’excitateur est configuré pour générer : - un premier peigne de fréquences 405 qui permet des transitions entre le premier niveau 401 et les troisièmes niveaux 403 ; il est configuré pour exciter le dispositif auxiliaire et, éventuellement, décaler de « p » le nombre de bosons de l’oscillateur harmonique quantique ; si p est positif des bosons sont ajoutés par le premier peigne 405, si p est négatif des bosons sont enlevés par ce peigne. - un deuxième peigne 406 qui permet des transitions entre les troisièmes niveaux 403 et le quatrième niveau 404 ; il est configuré pour désexciter le dispositif auxiliaire, exciter le dissipateur et décaler le nombre de bosons de manière à ce que le quatrième niveau ait un boson supplémentaire par rapport au premier niveau. L’état 404 se désexcite spontanément vers l’état 402 grâce à la relaxation spontanée du niveau vers le niveau selon le taux de désexcitation . Cette désexcitation est représentée par une flèche ondulée 407. L’excitateur génère les premier et deuxième peignes de fréquences de sorte qu’un premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences est égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne. Le premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences est égal au double du décalage de fréquence et est égal au double du deuxième écart fréquentiel séparant deux raies successives du deuxième peigne. L’effet de cette différence d’écart fréquentiel est de diminuer le canal d’altération 102 qui correspond à l’ajout de bosons dans l’oscillateur harmonique lorsque celui-ci comprend un nombre de bosons présentant la parité prédéterminée. Le premier peigne de fréquences 405 comprend des raies de fréquence qui sont en résonance avec les transitions d’un niveau 401 à un niveau 403, uniquement lorsque le niveau 401 correspond à un nombre k impair de bosons. Le premier peigne de fréquences 405 est donc configuré pour comprendre des raies de fréquences égales à l’expression , avec k impair. Le premier peigne de fréquences est un peigne de fréquences dont les raies sont séparées d’un intervalle fréquentiel égal à .Le deuxième peigne de fréquences 406 comprend des raies de fréquence qui sont en résonance avec les transitions d’un niveau 403 à un niveau 404, pour des valeurs du nombre k de bosons pair et impair. Le deuxième peigne de fréquences 406 est donc configuré pour comprendre des raies de fréquences égales à l’expression , avec k pair ou impair. Le deuxième peigne de fréquences est un peigne de fréquences dont les raies sont séparées d’un intervalle fréquentiel égal à . Les peignes de fréquences 405 et 406 permettent de réaliser un déplacement 408 correspondant au canal de correction 101 pour déplacer l’état du système depuis le niveau 401 dans le niveau 402, lorsque est impair. Le hamiltonien de couplage correspondant à l’action du premier peigne de fréquences 405 peut être formulé comme suit en représentation interaction : Le terme désigne une somme sur des entiers k. La fonction caractérise l’évolution temporelle en représentation interaction du premier peigne de fréquences 405 généré par l’excitateur 204. Cette fonction dépend du nombre de bosons k car la représentation interaction peut être vue en première approximation comme un changement de référentiel, le référentiel dépendant du nombre de bosons k dans l’oscillateur harmonique quantique. Les coefficients sont des nombres complexes dépendant de la manière dont le premier peigne de fréquences 405 se couple au système 200. La fonction du premier peigne de fréquences 405 peut être formulée sous la forme suivante : dans laquelle est un nombre complexe caractérisant la raie du premier peigne. La fréquence de la raie du premier peigne est égale à . L’amplitude de est une fréquence de couplage caractérisant l’amplitude de la raie du premier peigne. La phase de est une phase caractérisant la phase de la raie du premier peigne. Le terme désigne une somme sur l’entier k’ qui peut être positif, négatif ou nul. Le hamiltonien de couplage correspondant à l’action du deuxième peigne de fréquences 406 peut être formulé comme suit en représentation interaction : La fonction caractérise l’évolution temporelle en représentation interaction du deuxième peigne de fréquences 406 généré par l’excitateur 204. Cette fonction dépend du nombre de bosons k car la représentation interaction peut être vue en première approximation comme un changement de référentiel, le référentiel dépendant du nombre de bosons k dans l’oscillateur harmonique quantique. Les coefficients sont des nombres complexes dépendant de la manière dont le deuxième peigne de fréquences 406 se couple au système 200. La fonction du deuxième peigne de fréquences 406 peut être formulée sous la forme suivante : dans laquelle est un nombre complexe caractérisant la raie du deuxième peigne. La fréquence de la raie du deuxième peigne est égale à . L’amplitude de est une fréquence de couplage caractérisant l’amplitude de la raie du deuxième peigne. La phase de est une phase caractérisant la phase de la raie du deuxième peigne. L’effet de la différence d’écart fréquentiel entre le premier peigne de fréquences 405 et le deuxième peigne de fréquences 406 peut être compris à partir de la dépendance temporelle en représentation interaction du premier peigne et du deuxième peigne dans les cas où k=1 et k=2, en rapport avec la figure 4. La figure 4 comprend quatre graphes 4A, 4B, 4C et 4D, représentant respectivement de manière schématique les fonctions , , et . Pour cet exemple et pour simplifier l’étude, il a été supposé que : La fonction permet de comprendre le couplage entre le niveau 401 et le niveau 403, quand k=1, c’est-à-dire dans une situation où k est impair. La fonction permet de comprendre le couplage entre le niveau 401 et le niveau 403, quand k=2, c’est-à-dire dans une situation où k est pair. La fonction permet de comprendre le couplage entre le niveau 403 et le niveau 404, quand k=1, c’est-à-dire dans une situation où k est impair. La fonction permet de comprendre le couplage entre le niveau 403 et le niveau 404, quand k=2, c’est-à-dire dans une situation où k est pair. Dans ces graphes 4A-4D, les lignes continues représentent la partie réelle et les lignes tiretées la partie imaginaire des fonctions , , et . Si dans le domaine fréquentiel un peigne de fréquences est constitué de raies régulièrement séparées par un écart fréquentiel, dans le domaine temporel ce peigne de fréquences est constitué de pics séparés par un écart temporel inversement proportionnel à l’écart fréquentiel. On utilise le terme de raie pour le domaine fréquentiel et le terme de pic pour le domaine temporel. Dans le graphe 4A, la partie réelle de présente des pics régulièrement espacés par un intervalle de temps égal à . Il s’agit des pics composant le premier peigne de fréquences dans sa représentation temporelle. Les pics successifs sont identiques. L’effet associé à chaque pic de est de transférer une partie de la population du premier niveau 401 dans le troisième niveau 403. Dans le graphe 4B, la partie réelle de présente des pics régulièrement espacés par un intervalle de temps égal à . Il s’agit des pics composant le deuxième peigne de fréquences dans sa représentation temporelle. Les pics successifs sont identiques. L’effet associé à chaque pic de est de transférer une partie de la population du troisième niveau 403 au quatrième niveau 404. L’effet combiné des pics de et est, comme attendu, de réaliser le canal de correction 101, c’est-à-dire ajouter un boson dans l’oscillateur harmonique quand initialement celui-ci contient un nombre impair de bosons, en l’occurrence un nombre impair égal à un. Dans le graphe 4C, la partie réelle de présente des pics régulièrement espacés par un intervalle de temps égal à . Contrairement aux pics de qui sont tous de même signe, les pics de ont des signes alternés. Plus précisément deux pics successifs (501 et 502 sur la figure 5) présentent des signes différents. Dans le graphe 4D, la partie imaginaire de présente des pics 503 régulièrement espacés par un intervalle de temps égal à . L’effet associé à chaque pic positif 501 de est d’effectuer un transfert d’une partie de la population du premier niveau 401 dans le troisième niveau 403, et l’effet associé à chaque pic négatif 502 de est d’effectuer le transfert exactement opposé, au sens d’une rotation d’angle exactement opposé dans la sphère de Bloch, d’une partie de la population du troisième niveau 403 dans le premier niveau 401. L’effet associé à chaque pic 503 de est d’effectuer un transfert d’une partie de la population du troisième niveau 403 au quatrième niveau 404. Dans le cas de la figure 4, présente un écart entre ses pics égal à , alors que présente un écart entre ses pics égal à . Cet écart temporel double est la conséquence du fait que les raies du premier peigne de fréquences 405 sont espacées d’un écart en fréquence égal au double de l’écart en fréquence des raies du deuxième peigne de fréquences 406. De cette manière, un couple pic positif 501 / pic négatif 502 de se trouve toujours entre deux pics successifs 503 de . Ainsi, l’action du premier pic positif 501 du premier peigne est annulée par l’action du deuxième pic négatif 502. Comme il n’y a pas de pic 503 du deuxième peigne sur l’intervalle de temps séparant ces deux pics du premier peigne, le deuxième peigne n’interfère pas avec le premier peigne sur cet intervalle. Cela permet de diminuer le canal d’altération 102. Le canal d’altération est même potentiellement complètement supprimé en l’absence d’autres effets parasites comme, par exemple, les excitations thermiques du dispositif auxiliaire. Décalage temporel entre les peignes de fréquences Il est possible de configurer l’excitateur de sorte qu’il génère le premier et le deuxième peignes de fréquences avec les pics du premier peigne décalés temporellement par rapport aux pics du deuxième peigne. L’écart fréquentiel séparant deux raies successives du premier peigne de fréquences est égal à et l’écart fréquentiel séparant deux raies successives du deuxième peigne est égal à . Le décalage temporel peut être caractérisé par la durée D séparant un pic négatif 502 du premier peigne de fréquences 405 et le premier pic 503 du deuxième peigne de fréquence 406 qui suit temporellement le pic négatif 502. Le décalage temporel D peut être choisi librement dans l’intervalle . Le décalage temporel D peut être avantageusement choisi dans l’intervalle ou l’intervalle où et désignent respectivement la largeur à mi-hauteur des pics du premier peigne et du deuxième peigne. Comme précédemment présenté, si l’écart fréquentiel du premier peigne est le double de l’écart fréquentiel du deuxième peigne, alors il y a un écart entre pics de la fonction double de celui de la fonction . Si les pics du premier peigne de fréquences 405 sont décalés temporellement par rapport aux pics du deuxième peigne de fréquences 406, alors il n’y a pas de coïncidence temporelle entre un pic du premier peigne de fréquences 405 et un pic du deuxième peigne de fréquences 406. De cette manière, il n’y a pas de pic 503 du deuxième peigne de fréquences 406 sur l’intervalle de temps séparant un pic 501 du pic suivant 502 du premier peigne de fréquences 405. L’effet du pic 501 du premier peigne de fréquences 405 est annulé par l’effet du pic suivant 502 du premier peigne de fréquences 405, de sorte qu’à l’arrivée du pic 503 du deuxième peigne de fréquences 406, il n’y a aucune population dans le niveau 343 que ce pic peut transférer vers le niveau 404 lorsque le nombre de bosons k est pair. Sans coïncidence temporelle entre un pic du premier peigne de fréquences 405 et un pic du deuxième peigne de fréquences 406, il n’est pas possible non plus qu’il y ait un transfert de population du niveau 401 vers le niveau 403 par le pic du premier peigne qui soit simultanément transféré vers le niveau 404 par le pic du deuxième peigne de fréquences 406. Le canal d’altération 102 est encore diminué par cette caractéristique de non-coïncidence des pics du premier et du deuxième peignes de fréquences. Nombre de raies dans les peignes de fréquences Comme précédemment décrit, le premier peigne de fréquences 405 comprend des raies de fréquence qui sont en résonance avec les transitions du premier niveau 401 au troisième niveau 403, uniquement lorsque le niveau 401 correspond à un nombre k impair de bosons. Ces raies se trouvent aux fréquences où est un entier numérotant la raie. Elles sont séparées d’un intervalle fréquentiel égal à . Si le code de correction d’erreur requiert que soient stabilisés les états de Fock de parité paire ayant un nombre de bosons compris entre un entier plancher pair et un entier plafond pair , alors l’ensemble d’entiers est égal aux entiers compris entre les entiers et . Dans ce cas, le premier peigne de fréquences 405 doit comprendre au minimum les raies correspondant aux entiers , tels que . C’est notamment le cas lorsque l’information est encodée dans des états de l’oscillateur correspondant à un nombre de bosons compris entre un entier plancher pair et un entier plafond pair . Le premier peigne de fréquences 405 peut comprendre avantageusement davantage de raies et notamment comprendre les raies , tels que où est un entier impair inférieur à et est un entier impair supérieur à . Le nombre peut être négatif. De cette manière, le premier peigne de fréquences 405 ne comprend aucune raie résonante entre le premier niveau 401 et la pluralité de troisièmes niveaux 403 lorsque le nombre de bosons k du premier niveau 401 est pair. Plus le nombre de raies comprises dans le premier peigne de fréquences 405 est important et meilleure est la définition des pics 501 et 502 de sa représentation dans le domaine temporel. Cela permet notamment de faciliter et améliorer la condition de non-coïncidence entre les pics 501 et 502 du premier peigne de fréquences 405 et les pics 503 du deuxième peigne de fréquences 406. Comme précédemment décrit, le deuxième peigne de fréquences 406 comprend des raies de fréquence qui sont en résonance avec les transitions du troisième niveau 403 au quatrième niveau 404, lorsque le nombre k est pair ou impair. Ces raies se trouvent aux fréquences où est un entier numérotant la raie. Elles sont séparées d’un intervalle fréquentiel égal à . Si le code de correction d’erreur requiert que soient stabilisés les états de Fock de parité paire ayant un nombre de bosons compris entre l’entier plancher pair et l’entier plafond pair , alors le deuxième peigne de fréquences 406 doit comprendre au minimum les raies . Le deuxième peigne de fréquences 406 peut comprendre avantageusement davantage de raies et notamment comprendre les raies où est un entier pair ou impair inférieur à et est un entier pair ou impair supérieur à . Le nombre peut être négatif. De cette manière, le deuxième peigne de fréquences 406 comprend toujours des pics 503 temporels espacés d’un écart temporel double de l’écart temporel entre les pics 501 et 502 du premier peigne de fréquences 405. Plus le nombre de raies dans le deuxième peigne de fréquences 406 est important et meilleure est la définition des pics 503 de sa représentation dans le domaine temporel. Cela permet notamment de faciliter et améliorer la condition de non-coïncidence entre les pics 501 et 502 du premier peigne de fréquences 405 et les pics 503 du deuxième peigne de fréquences 406. Amplitude des raies des peignes de fréquences Chaque pic 501, 502 du premier peigne de fréquences 405 ou chaque pic 503 du deuxième peigne de fréquence 406 engendre un transfert de population entre des niveaux de l’architecture telle qu’on a pu les présenter plus haut. Un régime intéressant est obtenu lorsque le transfert engendré par les deux pics successifs 501 et 502 et le transfert engendré par le pic 503 correspondent chacun à une rotation dans la sphère de Bloch dont l’angle est plus petit que π. Ce régime correspond à avoir des peignes dont les fréquences de couplage et caractérisant les amplitudes des raies de fréquence sont plus petites que le quart du deuxième écart fréquentiel séparant deux raies successives du deuxième peigne de fréquences 406. Ce régime étant caractérisé par des raies ayant de petites amplitudes, donc par un excitateur envoyant des peignes ayant une puissance réduite, il est plus aisé à réaliser expérimentalement et est moins susceptible de provoquer des effets indésirables comme un échauffement du système, la saturation de ses composants non-linéaires ou l’apparition d’harmoniques. De plus, il existe dans ce régime un unique maximum du temps de vie des états ayant un nombre pair de bosons en fonction des fréquences de couplage et . En effet, augmenter ces fréquences permet d’augmenter le taux de transfert effectif 408 du premier niveau 401 au deuxième niveau 402 lorsque le nombre de bosons k dans le premier niveau est impair. Cela résulte en une augmentation du taux de correction du processus de correction 101. Cependant, comme cela a été précédemment mentionné, les imperfections expérimentales peuvent rendre non nul le taux de transfert effectif 408 du premier au deuxième niveau lorsque le nombre de bosons k dans le premier niveau est pair. Dans ce cas, augmenter les fréquences de couplage et résulte aussi en une augmentation du taux d’altération du processus d’altération 102. Il a été mentionné que le processus d’altération 102 pouvait ne pas être corrigé par le processus de correction 101, mais au contraire aggravé. Dans le régime où les fréquences de couplage et sont plus petites que , il existe ainsi un réglage optimal de l’amplitude des raies correspondant à un maximum unique du temps de vie des états ayant un nombre pair de bosons. Le fait que ce maximum soit unique permet de trouver à coup sûr le réglage optimal, aussi bien expérimentalement que théoriquement, par un algorithme d’optimisation comme, par exemple, l’algorithme du gradient. Une telle optimisation de l’amplitude des raies des premier et deuxième peignes de fréquences, l’amplitude des raies étant caractérisée par une fréquence de couplage, peut notamment être réalisée par les étapes suivantes : - (E1) choisir une valeur de fréquence de couplage inférieure au quart du deuxième écart fréquentiel, - (E2) mesurer un temps de vie des états ayant un nombre pair de bosons, - (E3) modifier la valeur de la fréquence de couplage en fonction du temps de vie mesuré selon un algorithme d’optimisation du temps de vie de sorte à identifier une valeur optimale de la fréquence de couplage correspondant à un maximum du temps de vie. Un exemple d’algorithme d’optimisation consiste à suivre les étapes suivantes : - (F1) choisir un entier a supérieur à 2 et calculer un écart de fréquence de couplage égal au rapport du quart du deuxième écart fréquentiel sur l’entier a, - (F2) réaliser une boucle de sous-étapes indicées par l’entier i initialement égal à un : - (F21) calculer une valeur FC(i) de fréquence du couplage égale au produit de l’entier i par l’écart de fréquence de couplage, les premier et deuxième peigne de fréquences générés et envoyés vers au moins l’un parmi l’oscillateur harmonique quantique, le dispositif auxiliaire et le dissipateur, présentant une amplitude des raies caractérisée par la valeur de fréquence de couplage calculée, - (F22) mesurer un temps de vie TV(i) des états de l’oscillateur ayant un nombre pair de bosons, - (F23) si i est égal à 1 remplacer l’entier i par l’entier i+1 et répéter les étapes F21 et F22, sinon - (F24) lorsque TV(i-1) est inférieur à TV(i), remplacer l’entier i par l’entier i+1 et répéter les étapes F21 à F24, sinon - (F3) arrêter la boucle de sous-étapes et restituer sur une interface la valeur FC(i-1) de fréquence du couplage correspondant à un temps de vie maximum TV(i-1) des états de l’oscillateur. L’interface peut notamment être une interface d’entrée et de sortie permettant la saisie des données ou de restituer les résultats d’un procédé. Il est à noter que plus l’entier a est important et meilleure est la précision de la fréquence du couplage recherchée correspondant au temps de vie maximum. Cependant, au moins en raison de l’erreur de mesure du temps de vie, il existe un entier a au-delà duquel l’application des étapes précédentes n’améliore la précision sur la fréquence du couplage recherchée. Taux de relaxation du dissipateur Le taux de transfert effectif 308 est ultimement limité par le taux de désexcitation , égal à la probabilité de passer du niveau du dissipateur vers le niveau par unité de temps. Il en résulte que le taux de correction est nécessairement inférieur au taux de désexcitation Augmenter peut donc permettre d’augmenter le taux de correction . Cependant, en vertu de la relation d’incertitude temps.énergie de Heisenberg, augmenter entraîne un élargissement du troisième niveau 403 proportionnel à . Lorsque cet élargissement est de l’ordre de la quantité , les niveaux correspondant à un nombre de bosons pair ou impair se recouvrent. Ce recouvrement peut potentiellement augmenter le taux d’altération et nuire à la sélectivité. Dans le régime mentionné précédemment où les fréquences de couplage et sont plus petites que le quart du deuxième écart fréquentiel séparant deux raies successives du deuxième peigne de fréquences 406, l’augmentation simultanée du taux de correction et du taux d’altération avec résulte en l’existence d’un unique maximum du temps de vie des états ayant un nombre pair de bosons en fonction de . Le fait que ce maximum soit unique permet de trouver à coup sûr la valeur optimale de , aussi bien expérimentalement que théoriquement, par un algorithme d’optimisation comme, par exemple, l’algorithme du gradient. Une telle optimisation du taux de désexcitation peut notamment être réalisée par la réalisation des étapes suivantes : - (S1) choisir une valeur du taux de désexcitation inférieure au deuxième écart fréquentiel, - (S2) mesurer un temps de vie des états ayant un nombre pair de bosons, - (S3) modifier la valeur du taux de désexcitation en fonction du temps de vie mesuré selon un algorithme d’optimisation du temps de vie de sorte à identifier une valeur optimale du taux de désexcitation correspondant à un maximum du temps de vie. Un exemple d’algorithme d’optimisation consiste à suivre les étapes suivantes - (T1) choisir un entier b supérieur à 2 et calculer un écart de taux de désexcitation égal au rapport du deuxième écart fréquentiel sur l’entier b, - (T2) réaliser une boucle de sous-étapes indicées par l’entier i initialement égal à un : - (T21) calculer une valeur TD(i) du taux de désexcitation égale au produit de l’entier i par l’écart de taux de désexcitation, le dispositif auxiliaire étant couplé au dissipateur configuré pour passer d’un état excité |1> dans un état fondamental |0> selon un taux de désexcitation égal à la valeur TD(i), - (T22) mesurer un temps de vie TV(i) des états de l’oscillateur ayant un nombre pair de bosons , - (T23) si i est égal à 1 remplacer l’entier i par l’entier i+1 et répéter les étapes T21 et T22, sinon - (T24) lorsque TV(i-1) est inférieur à TV(i), remplacer l’entier i par l’entier i+1 et répéter les étapes T21 à T24, sinon - (T3) arrêter la boucle de sous-étapes et restituer sur une interface la valeur TD(i-1) du taux de désexcitation correspondant à un temps de vie maximum TV(i-1) des états de l’oscillateur L’interface peut notamment être l’interface précédemment évoqué pour la recherche de la fréquence de couplage correspondant au temps de vie maximum. Il est à noter que plus l’entier b est important et meilleure est la précision du taux de désexcitation recherché correspondant au temps de vie maximum. Cependant, au moins en raison de l’erreur de mesure du temps de vie, il existe un entier b au-delà duquel l’application des étapes précédentes n’améliore la précision sur le taux de désexcitation recherché. Autres architectures de niveaux possibles Il a été proposé jusqu’à présent d’utiliser le premier niveau ou état référencé 401, une pluralité de troisièmes niveaux ou états référencée 403, un quatrième niveau référencé 404, et un deuxième niveau référencé 402. Plus généralement, il est possible d’utiliser un quatrième niveau comportant un nombre de bosons k+q différent de k+1, avec q un entier positif ou négatif. Lorsque l’entier q est différent de 1, le passage de l’état |1> à l’état |0> du dissipateur avec le taux s’accompagne alors d’un décalage q-1 du nombre de bosons. Dans ce cas, les fréquences des raies du deuxième peigne sont égales à où est un entier. Technologie des circuits supraconducteurs Le procédé de stabilisation décrit ici peut, par exemple, être réalisé en utilisant la technologie des circuits supraconducteurs. L’oscillateur harmonique quantique peut alors être un résonateur micro-onde supraconducteur. Le système auxiliaire peut être un transmon. Le dissipateur peut être un autre résonateur micro-onde supraconducteur ayant, comparativement à l’oscillateur, un mauvais facteur de qualité. Dans ce cas, le hamiltonien du système peut être écrit de la manière suivante : où est la fréquence de résonance du résonateur micro-onde supraconducteur jouant le rôle de l’oscillateur harmonique quantique encodant le qubit logique ; sont ses opérateurs de création et annihilation vérifiant ; est la fréquence de transition fondamentale du transmon ; son anharmonicité ; sont ses opérateurs de création et annihilation vérifiant ; est la fréquence de résonance du résonateur micro-onde supraconducteur jouant le rôle du dissipateur ; sont ses opérateurs de création et annihilation vérifiant ; est un décalage de fréquence du premier niveau excité du transmon par boson ajouté dans l’oscillateur harmonique quantique. De plus, les opérateurs des trois composantes du système tripartite commutent : Le processus de relaxation 100 de l’oscillateur harmonique quantique et le processus de désexcitation 407 du dissipateur peuvent être décrits, respectivement, par les opérateurs sauts suivants : Où est le taux de relaxation et est le taux de désexcitation. L’état de l’oscillateur harmonique quantique correspond à l’état contenant photons dans le résonateur micro-onde supraconducteur. Les états et du dissipateur correspondent aux états contenant respectivement 0 et 1 photon dans le résonateur micro-onde supraconducteur ayant un mauvais facteur de qualité. L’état du dispositif auxiliaire correspond à l’état fondamental du transmon ( ). L’état peut correspondre à un état excité parmi plusieurs états excités du transmon. Il est possible d’obtenir différents modes de réalisation en faisant varier : l’état excité du transmon choisi pour jouer le rôle de l’état [Math. 301] et la valeur de l’entier p qui définit la pluralité des troisièmes niveaux ou états [Math. 302] . La table 1 liste trois modes de réalisation correspondant à un couple particulier d’état excité et d’entier p . Une fois ce couple fixé, il est possible de déduire, pour chaque mode les expressions des termes , , , , , et . Ces termes ont été introduits précédemment à l’exception des termes et qui sont tels que et que , avec et les fonctions définissant respectivement les premier et deuxième peignes de fréquences en représentation de Schrödinger. Mode 1 Mode 2 Mode 3 Niveau excité du transmon 1 1 0 q 1 1 1 Table 1 : Exemples de modes de réalisation de l’architecture de niveaux pour le protocole de stabilisation de parité L’état excité est défini dans la table 1 comme état propre de l’opérateur . Il est possible, en choisissant d’autres couples état excité / entier p , de déterminer d’autres modes de réalisation. Selon l’expression du et du obtenue par le choix du couple état excité / entier p , le mode de réalisation est plus ou moins simple à réaliser. Dans le cas où l’expression d’un opérateur ou fait intervenir le produit d’un ou de trois opérateurs de création et d’annihilation et , cet opérateur peut être réalisé en générant par l’excitateur un signal micro-onde excitant le transmon (voir par exemple le document précédemment cité de GERTLER Jeffrey M. et PECHAL M. et al. “ Microwave-Controlled Generation of Shaped Single Photons in Circuit Quantum Electrodynamics ” Physical Review X, Vol. 4, page 041010) et ayant des raies aux fréquences : pour un opérateur c’est-à-dire pour le premier peigne et pour un opérateur c’est-à-dire pour le deuxième peigne. Dans le cas où l’expression d’un opérateur ou fait intervenir le produit de deux opérateurs de création et d’annihilation et , cet opérateur peut être réalisé en générant par l’excitateur deux signaux micro-onde excitant le transmon, comme présenté dans CAMPAGNE-IBARCQ P. et al. “ Deterministic Remote Entanglement of Superconducting Circuits through Microwave Two-Photon Transitions ”, Physical Review Letters Vol. 120 page 200501. Pour un opérateur faisant intervenir le produit de deux opérateurs création et annihilation, les deux signaux micro-onde doivent avoir un produit qui est un signal présentant des raies aux fréquences . Par exemple les deux signaux micro-onde peuvent être : - d’une part un peigne de fréquences comprenant les raies (où est une fréquence constante) : , et - d’autre part un signal monochromatique à la fréquence . Pour un opérateur faisant intervenir le produit de deux opérateurs création et annihilation, les deux signaux micro-onde doivent avoir un produit qui est un signal présentant des raies aux fréquences . Par exemple les deux signaux micro-onde peuvent être : - d’une part un peigne de fréquences comprenant les raies (où est une fréquence constante) : , et - d’autre part un signal monochromatique à la fréquence . Système de stabilisation (200) autonome d’états quantiques comprenant : - un oscillateur harmonique quantique (202) configuré pour porter une information encodée dans des états de l’oscillateur ayant une parité prédéterminée d’un nombre de bosons compris dans l’oscillateur harmonique (202), - un dispositif auxiliaire (201) configuré pour être dans un état fondamental |g> ou dans un état excité |f> parmi une pluralité d’états excités, le dispositif (201) étant couplé à l’oscillateur (202) de manière dispersive, de sorte qu’une différence de fréquence entre l’état fondamental |g> et l’état excité |f> dépend de manière linéaire d’un décalage de fréquence χ et du nombre de bosons, - un dissipateur (203) configuré pour passer d’un état excité |1> dans un état fondamental |0>, - un excitateur (204) configuré pour générer des premier (205, 405) et deuxième (206, 406) peignes de fréquences, destinés à être envoyés vers au moins l’un parmi l’oscillateur harmonique quantique (202), le dispositif auxiliaire (201) et le dissipateur (203) de sorte à ajouter un boson dans l’oscillateur harmonique (202) si le nombre de bosons n’a pas la parité prédéterminée, caractérisé en ce qu’un premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences (205, 405) est égal au double du décalage de fréquence χ et est égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne (206, 406). Système selon la revendication 1 dans lequel les premier et deuxième peignes présentent des pics dans le domaine temporel, les pics du premier peigne étant décalés temporellement par rapport aux pics du deuxième peigne. Système selon la revendication 1 ou 2 configuré en outre pour être, pour chaque entier k d’un ensemble d’entiers, dans : - un premier niveau (401) défini par un premier nombre de bosons égal à l’entier k, l’état fondamental |g> et l’état fondamental |0>, - un deuxième niveau (402) défini par un deuxième nombre de bosons égal à l’entier k+1, l’état excité |g> et l’état fondamental |0>, - un troisième niveau (403) dans lequel le dispositif auxiliaire est dans l’état excité |f> et le dissipateur est dans l’état stable |0>, - un quatrième niveau (404) dans lequel le dispositif auxiliaire est dans l’état fondamental |g> et le dissipateur est dans l’état excité |1>, l’excitateur (204) étant configuré pour coupler, pour chaque entier k, les premier, troisième et quatrième niveaux si l’entier k n’a pas la parité prédéterminée, de sorte à déplacer le système du premier niveau (401) dans le deuxième niveau (402), l’excitateur (204) étant configuré pour générer les premier et deuxième peignes de fréquences de sorte que, pour chaque entier k : - le premier peigne (405) comprend une première raie si l’entier k n’a pas la parité prédéterminée, la première raie ayant une fréquence égale à une fréquence de transition entre les premier et troisième niveaux, et - le deuxième peigne (406) comprend une deuxième raie ayant une fréquence égale à une fréquence de transition entre les troisième et quatrième niveaux. Système selon la revendication 3 dans lequel l’oscillateur harmonique (202) comprend le deuxième nombre de bosons égal à l’entier k+1 dans le quatrième niveau (402). Système selon l’une des revendications 1 à 4 dans lequel l’oscillateur harmonique quantique (202) est un premier résonateur micro-onde supraconducteur, le système auxiliaire (201) est un transmon, le dissipateur (203) est un deuxième résonateur micro-onde supraconducteur ayant un facteur de qualité inférieur à un facteur de qualité du premier résonateur micro-onde supraconducteur. Système selon l’une des revendications 1 à 5 dans lequel une amplitude des raies du premier peigne de fréquences et du deuxième peigne de fréquences est caractérisée par une fréquence de couplage, la fréquence de couplage étant inférieure au quart du deuxième écart fréquentiel. Système selon la revendication 3 ou 4 dans lequel l’information est encodée dans des états de l’oscillateur correspondant à un nombre de bosons compris entre un entier plancher ayant la parité prédéterminée et un entier plafond ayant la parité prédéterminée, l’ensemble d’entiers étant égal aux entiers compris entre les entiers et . Procédé de stabilisation autonome d’états quantiques comprenant les étapes suivantes : - encoder une information dans des états d’un oscillateur harmonique quantique, les états ayant une parité prédéterminée d’un nombre de bosons compris dans l’oscillateur harmonique, - coupler de manière dispersive l’oscillateur harmonique à un dispositif auxiliaire, de sorte qu’une différence de fréquence entre un état fondamental |g> du dispositif auxiliaire et un état excité |f> du dispositif auxiliaire dépend de manière linéaire d’un décalage de fréquence χ et du nombre de bosons, - coupler le dispositif auxiliaire à un dissipateur configuré pour passer d’un état excité |1> dans un état fondamental |0> - générer des premier et deuxième peignes de fréquences, un premier écart fréquentiel séparant deux raies successives du premier peigne de fréquences étant égal au double du décalage de fréquence χ et égal au double d’un deuxième écart fréquentiel séparant deux raies successives du deuxième peigne et - envoyer les premier et deuxième peignes de fréquences vers au moins l’un parmi l’oscillateur harmonique quantique, le dispositif auxiliaire et le dissipateur de sorte à ajouter un boson dans l’oscillateur harmonique si le nombre de bosons n’a pas la parité prédéterminée. Procédé selon la revendication 8 comprenant en outre une optimisation de l’amplitude des raies des premier et deuxième peignes de fréquences, l’amplitude des raies étant caractérisée par une fréquence de couplage, comprenant les étapes suivantes : - (E1) choisir une valeur de fréquence de couplage inférieure au quart du deuxième écart fréquentiel, - (E2) mesurer un temps de vie des états ayant un nombre de bosons de la parité prédéterminée et - (E3) modifier la valeur de la fréquence de couplage en fonction du temps de vie mesuré selon un algorithme d’optimisation du temps de vie de sorte à identifier une valeur optimale de la fréquence de couplage correspondant à un maximum du temps de vie. Procédé selon la revendication 8 comprenant en outre une optimisation d’un taux de désexcitation égal à la probabilité de passer du niveau vers le niveau par unité de temps, comprenant les étapes suivantes : - (S1) choisir une valeur du taux de désexcitation inférieure au deuxième écart fréquentiel, - (S2) mesurer un temps de vie des états ayant un nombre de bosons de la parité prédéterminée et - (S3) modifier la valeur du taux de désexcitation en fonction du temps de vie mesuré selon un algorithme d’optimisation du temps de vie de sorte à identifier une valeur optimale du taux de désexcitation correspondant à un maximum du temps de vie.