Real-time Krylov Diagonalisation for Open Quantum Systems

Cet article présente une modification des méthodes de sous-espace quantique en temps réel pour simuler des systèmes quantiques ouverts sous forme de Lindblad, démontrant leur application à un résonateur Kerr supraconducteur piloté par deux photons afin d'estimer le gap de Liouvillien dans le régime du qubit Chat de Kerr.

L'essentiel

La vue d'ensemble : Écouter une pièce bruyante

Imaginez que vous essayez de comprendre la musique qui joue dans une pièce très bruyante et chaotique. Dans le monde de la physique quantique, cette « pièce » est un système quantique ouvert — une machine qui perd constamment de l'énergie ou qui est perturbée par son environnement (comme la chaleur ou le frottement).

Habituellement, les scientifiques utilisent des ordinateurs puissants pour simuler ces systèmes. Mais à mesure que ces systèmes grandissent, les ordinateurs classiques se retrouvent bloqués. Ce document propose une nouvelle façon d'utiliser les ordinateurs quantiques pour écouter ces systèmes en temps réel, spécifiquement pour déterminer à quelle vitesse ils se stabilisent ou combien de temps ils peuvent conserver l'information.

Le problème : Le raccourci « Krylov »

Pour comprendre le document, vous devez d'abord connaître un outil appelé diagonalisation du sous-espace de Krylov.

• L'analogie : Imaginez que vous voulez connaître la hauteur exacte d'une note de guitare, mais que vous ne pouvez pas mesurer la corde directement. Au lieu de cela, vous pincez la corde et vous écoutez le son qu'elle produit au fil du temps.
• L'ancienne méthode : Vous écoutez pendant un court instant, prenez une photo instantanée et essayez de deviner la note.
• La méthode du document : Vous écoutez l'évolution du son sur une longue période. Vous enregistrez une série de photos instantanées : le son à 1 seconde, 2 secondes, 3 secondes, etc. En examinant comment ces instantanés se rapportent les uns aux autres, vous pouvez reconstruire mathématiquement la « vraie » note de la corde sans avoir besoin de la mesurer directement.

En termes quantiques, cette méthode construit une « bibliothèque » des états passés du système (le sous-espace de Krylov) pour déterminer ses propriétés cachées.

La surprise : Gérer les systèmes « ouverts »

La plupart des ordinateurs quantiques sont conçus pour des systèmes « fermés » (parfaitement isolés, comme dans le vide). Mais les dispositifs quantiques réels sont « ouverts » — ils fuient de l'énergie et deviennent désordonnés.

L'auteur, D. A. Herrera-Martí, explique comment modifier la méthode d'« écoute » pour qu'elle fonctionne dans ces environnements désordonnés et ouverts.

• Le défi : Dans un système fermé, les ondes sonores rebondissent simplement d'avant en arrière. Dans un système ouvert, le son s'estompe (se dissipe).
• L'intuition : L'auteur a réalisé que, parce que le système s'estompe, vous pouvez en fait écouter plus longtemps que d'habitude.
  • Analogie : Imaginez une toupie qui tourne. Dans une pièce sans frottement, elle tourne éternellement. Dans une pièce avec frottement, elle ralentit et finit par s'arrêter. Si vous voulez étudier la partie la « plus lente » de la rotation (le moment juste avant qu'elle ne s'arrête), vous n'avez pas besoin de l'observer pendant un court éclat ; vous devez l'observer jusqu'à ce que les parties rapides et vacillantes disparaissent, ne laissant que le balancement lent et régulier.
  • Le document montre qu'en laissant le système quantique évoluer pendant plus de temps, le « bruit » rapide disparaît, et les signaux « lents » et importants deviennent clairs.

Le cas de test : Le qubit « Chat de Kerr »

Pour prouver que cela fonctionne, l'auteur l'a testé sur un type spécifique de bit quantique appelé qubit Chat de Kerr.

• Qu'est-ce que c'est ? Imaginez ce qubit comme un pendule qui peut osciller dans deux directions à la fois (gauche ou droite). Il est conçu pour être très stable face aux erreurs.
• Le « gap » : En physique, il existe un concept appelé « gap » (écart). Imaginez une vallée entre deux collines. Le « gap » est la hauteur de la colline que vous devez gravir pour passer d'un côté à l'autre.
  • Si le gap est large, le système est stable et change lentement.
  • Si le gap est étroit (ou en train de se refermer), le système est au bord d'une transition de phase et devient très sensible.
• Le résultat : L'auteur a utilisé sa nouvelle méthode d'« écoute prolongée » pour mesurer cet écart. Il a constaté qu'en augmentant la puissance de l'entraînement (la « poussée » sur le pendule), l'écart devenait de plus en plus petit. Cela a confirmé que le système entrait dans un état spécial « protégé » où l'information est difficile à détruire.

Pourquoi cela compte (selon le document)

Le document ne prétend pas que cela construira immédiatement une meilleure intelligence artificielle ou guérira des maladies. Au contraire, il affirme :

1. Meilleurs outils : Nous avons maintenant un moyen d'utiliser les ordinateurs quantiques pour étudier les systèmes « bruyants » avec plus de précision qu'auparavant.
2. Compréhension de la stabilité : Nous pouvons mieux comprendre comment construire des ordinateurs quantiques qui ne se cassent pas facilement (comme le qubit Chat).
3. Efficacité : En écoutant plus longtemps, nous pouvons obtenir de meilleures réponses avec moins d'étapes mathématiques, ce qui est crucial car les ordinateurs quantiques sont actuellement très fragiles et sujets aux erreurs.

Résumé

Le document est comme un nouveau mode d'emploi pour un détective. Au lieu d'essayer d'attraper un voleur (l'état quantique) en prenant une photo rapide, le détective sait maintenant installer une caméra de surveillance à long terme. En observant comment les mouvements du voleur ralentissent et s'estompent au fil du temps, le détective peut identifier la véritable identité du voleur (les propriétés du système) beaucoup plus clairement, même dans une ville bondée et bruyante.

L'auteur a appliqué avec succès cette technique de « surveillance à long terme » à un dispositif quantique spécifique (le Chat de Kerr) et a prouvé qu'il pouvait mesurer les limites de stabilité de l'appareil, une étape cruciale pour la construction d'ordinateurs quantiques futurs.

Résumé technique

Résumé technique : Diagonalisation de Krylov en temps réel pour les systèmes quantiques ouverts

Énoncé du problème
L'article aborde le défi de la simulation des systèmes quantiques ouverts décrits par des équations maîtresses de Lindblad, en se concentrant spécifiquement sur l'estimation de la lacune de Liouvillien. Alors que la diagonalisation de sous-espace quantique (QSD) est devenue un algorithme de premier plan, non tolérant aux pannes, pour les systèmes fermés (hamiltoniens), son application aux systèmes dissipatifs n'est pas triviale. Les systèmes ouverts sont régis par des super-opérateurs de Liouvillien non hermitiens et non normaux, qui introduisent des différences structurelles telles que des vecteurs propres bi-orthogonaux, des valeurs propres complexes et une dynamique non unitaire. Ces caractéristiques compliquent les méthodes standard de sous-espace de Krylov, notamment en ce qui concerne la stabilité numérique (mauvais conditionnement) et l'interprétation des filtres spectraux. Les auteurs visent à adapter la QSD en temps réel pour estimer la lacune de Liouvillien, un paramètre critique définissant le temps de relaxation et la présence de modes lents dans les transitions de phase dissipatives.

Méthodologie
Les auteurs proposent une modification de la diagonalisation de sous-espace quantique (QSD) en temps réel pour gérer la dynamique de Liouvillien. La méthodologie centrale comprend :

1. Construction du sous-espace de Krylov : Au lieu d'évoluer un état sous un hamiltonien unitaire $H$, le système évolue sous le super-opérateur de Liouvillien $\mathcal{L}$. Un sous-espace de Krylov est généré en appliquant des puissances de l'opérateur d'évolution temporelle $e^{\mathcal{L}\tau}$ à un état initial $|\rho_0)$ :
   $$\text{span}\{|\rho_0), e^{\mathcal{L}\tau}|\rho_0), \dots, e^{\mathcal{L}(D-1)\tau}|\rho_0)\}$$
2. Problème de valeurs propres généralisé (GEVP) : Étant donné que les circuits quantiques ne peuvent pas effectuer de réorthogonalisation en temps réel, la méthode résout un GEVP $\bar{H}v = \lambda Sv$. Pour les Liouvilliens, la matrice de recouvrement $S$ et la matrice hamiltonienne $\bar{H}$ sont construites en utilisant des vecteurs propres gauches et droits bi-orthogonaux.
3. Interprétation du filtre spectral : La méthode est présentée comme l'optimisation d'un filtre spectral. Contrairement aux systèmes hamiltoniens où les états propres sont réels et l'évolution temporelle est oscillatoire (nécessitant un échantillonnage de Nyquist strict), la dynamique de Liouvillien est dissipative. Les modes à décroissance rapide spiralent vers l'intérieur dans le plan complexe, permettant des temps d'évolution $\tau$ plus longs sans artefacts de repliement. Cela permet d'isoler les modes lents (valeurs propres avec de petites parties réelles) même avec une dimension de Krylov $D$ réduite.
4. Gestion de la non-normalité : Les auteurs reconnaissent que les opérateurs non normaux entraînent une forte sensibilité aux perturbations (théorème de Bauer-Ficke). La matrice de recouvrement $S$ devient exponentiellement mal conditionnée à mesure que $D$ augmente. Pour atténuer cela, la méthode s'appuie sur la nature dissipative du système pour supprimer les modes rapides, permettant des $\tau$ plus grands et des $D$ plus petits, réduisant ainsi l'instabilité numérique associée aux grandes dimensions de Krylov.
5. Sélection de l'état initial : De manière cruciale, l'état initial doit avoir une projection non nulle sur le secteur de trace nulle (cohérences ou déséquilibres de population) pour exciter les modes lents. Le secteur de trace 1 (état stationnaire) ne contribue pas au signal pour l'estimation de la lacune dans ce cadre.

Contributions clés

• Extension de la QSD aux Liouvilliens : L'article fournit un cadre théorique pour appliquer la QSD en temps réel aux super-opérateurs de Liouvillien non hermitiens et non normaux sous forme de Lindblad.
• Stratégie d'estimation de la lacune : Il démontre comment estimer la lacune de Liouvillien ($\Delta = -\max \Re(\lambda_k)$) en exploitant la nature dissipative du système pour filtrer les modes à décroissance rapide, contournant efficacement certaines contraintes de Nyquist qui s'appliquent à l'évolution unitaire.
• Analyse de la stabilité numérique : Le travail met en évidence le compromis entre le temps d'évolution $\tau$ et la dimension de Krylov $D$. Il soutient que l'augmentation de $\tau$ (plutôt que de $D$) est une stratégie viable pour les systèmes ouverts afin d'améliorer la résolution de la lacune tout en maintenant la stabilité numérique, car les modes rapides décroissent naturellement.
• Application aux qubits de chat de Kerr : La méthode est appliquée à un résonateur supraconducteur de Kerr piloté à deux photons (qubit de chat de Kerr), un système pertinent pour les architectures de qubits protégés.

Résultats
Les auteurs ont appliqué leur méthode à un modèle simplifié d'un qubit de chat de Kerr avec perte de photon unique.

• Paramètres de simulation : Les simulations ont été réalisées avec un maximum de 30 photons et une dimension de Krylov $D=20$, résultant en une matrice de Krylov de $900 \times 20$.
• Reconstruction de la lacune : La méthode a permis de reconstruire avec succès la lacune de Liouvillien à travers différentes forces de pilotage. Les résultats ont montré que l'augmentation du temps d'évolution $\tau$ améliorait la précision de l'estimation de la lacune, en particulier dans les régimes où le commutateur de Liouvillien (une mesure de la non-normalité) était élevé.
• Identification des modes lents : Les simulations ont confirmé l'émergence d'un mode lent correspondant à la décroissance des cohérences et des déséquilibres de population. Dans le régime de « bruit biaisé », la lacune rétrécit exponentiellement avec l'amplitude de l'état cohérent ($\propto e^{-\alpha^2}$), indiquant la formation d'un sous-espace de calcul protégé.
• Représentation de Wigner : Les états propres reconstruits à partir de la matrice de Krylov ont été utilisés pour visualiser la représentation de Wigner du mode lent, montrant la transition d'un régime amorti à un régime avec deux lobes distincts (brisure de symétrie), puis à nouveau à un régime comprimé à fort pilotage.

Signification et revendications
L'article revendique démontrer la faisabilité de l'utilisation de la QSD en temps réel pour étudier les systèmes quantiques ouverts, spécifiquement pour estimer la lacune de Liouvillien dans le contexte des transitions de phase dissipatives. Les auteurs positionnent ce travail comme une adaptation nécessaire des algorithmes de l'ère NISQ pour la conception et la simulation de matériel quantique, qui implique intrinsèquement de la dissipation.

La signification est présentée autour du potentiel d'utiliser les ordinateurs quantiques pour simuler les modes lents des systèmes quantiques qui sont coûteux en calcul à simuler classiquement, en particulier dans le contexte des qubits protégés comme le qubit de chat de Kerr. Les auteurs notent modestement que si les qubits de chat de Kerr individuels sont simulables classiquement, les réseaux de tels qubits pourraient ne pas l'être, suggérant l'utilité de la méthode pour des systèmes plus grands et plus complexes. Le travail est présenté comme une étape vers la « découverte automatisée de matériel quantique » en fournissant des outils pour caractériser les propriétés de bruit et la stabilité des architectures de qubits proposées. L'article ne revendique pas avoir résolu le problème général de la simulation d'opérateurs non normaux, mais offre plutôt une approche spécifique et modifiée pour l'estimation de la lacune dans les systèmes dissipatifs.