Fast design and scaling of multi-qubit gates in large-scale trapped-ion quantum computers

Résumé Technique : Conception Rapide et Mise à l'Échelle de Portes Multi-qubits dans les Ordinateurs Quantiques à Ions Piégés à Grande Échelle

Énoncé du Problème
Les ordinateurs quantiques à ions piégés offrent une connectivité longue distance et des temps de cohérence longs, ce qui en fait une plateforme de premier plan pour le calcul quantique. Cependant, la mise à l'échelle de ces systèmes vers des centaines d'ions présente un défi conceptuel et technologique majeur. Bien que les ions possèdent une connectivité tous-à-tous via leurs modes de mouvement collectifs (bus de phonons), la conception de portes d'intrication multi-qubits programmables et de haute fidélité pour de grands cristaux est informatiquement prohibitive.

La difficulté centrale réside dans le problème de conception de la porte : trouver un ensemble de commandes de commande qui satisfont simultanément $N(N-1)/2$ contraintes quadratiques (phases d'intrication cibles entre toutes les paires de qubits) tout en minimisant la puissance de commande et en garantissant que les modes de mouvement reviennent à l'état fondamental. Cela se formule comme un problème d'optimisation sous contraintes quadratiques qui évolue en $N^2$ avec le nombre d'ions $N$, rendant le problème NP-difficile. Les approches existantes limitent souvent le problème à de petits sous-ensembles d'ions ou reposent sur des algorithmes de minimisation basés sur le gradient qui peinent face à la complexité et aux contraintes des systèmes à grande échelle.

Méthodologie : La Méthode Large-Scale Fast (LSF)
Les auteurs introduisent un cadre méthodologique nommé « Large-Scale Fast » (LSF) pour atténuer la nature NP-difficile du problème de conception de porte, permettant la conception de portes robustes et programmables pour des cristaux allant jusqu'à des centaines d'ions.

1. Formulation dans le Domaine Spectral : La conception de la porte est formulée dans le domaine spectral, représentant les champs de commande comme une somme de composantes de Fourier désaturées autour de la fréquence de transition du qubit. Le problème est réduit à la recherche d'un vecteur $R$ (amplitudes des tons de commande) qui minimise la norme $\|R\|$ sous des contraintes linéaires (déplacements motionnels nuls) et des contraintes quadratiques (phases d'intrication cibles $\phi_{nn'}$).
2. Graine de Solution à Phase Nulle (ZPS) : La méthode LSF repose sur la génération d'une solution « graine » appelée Solution à Phase Nulle (ZPS). Une ZPS est une solution non triviale du problème homogène où toutes les phases d'intrication cibles sont nulles ($\phi^{(0)}_s = 0$). Ces solutions sont accumulées avant de résoudre pour des portes cibles spécifiques.
3. Conversion et Optimisation : L'algorithme procède en deux étapes principales :
   • Conversion : Une ZPS est mise à l'échelle par un grand facteur $\lambda$ et perturbée par un petit vecteur $D$ pour satisfaire les contraintes quadratiques cibles. En linéarisant les contraintes autour de la ZPS mise à l'échelle, les auteurs dérivent un système linéaire pour résoudre $D$, convertissant la ZPS en une solution approximative $R_{conv}$ qui satisfait les phases cibles.
   • Réduction de Norme : Un processus itératif minimise la norme du vecteur solution tout en respectant strictement les contraintes. Cela implique l'alternance entre des étapes de « réduction de norme » (se déplacer le long de la variété de solution pour réduire la puissance) et des étapes de « réduction d'erreur » (corriger les écarts par rapport aux phases cibles). Les tailles de pas sont calculées analytiquement à l'aide de méthodes de recherche de racines polynomiales (méthodes de Cardano et de Ferrari) pour garantir que l'erreur reste dans une tolérance spécifiée (typiquement $10^{-4}$).
4. Agrégation : La méthode agrège plusieurs ZPS pour assurer un espace de recherche diversifié, permettant à l'algorithme de trouver des solutions efficaces pour diverses cartes d'intrication.

Contributions Clés et Résultats

• Évolutivité en Temps Polynomial : La méthode LSF transforme le problème de conception NP-difficile en une tâche de génération de ZPS (qui peut être faite efficacement pour des symétries spécifiques) suivie d'une optimisation efficace en temps polynomial. Les tests de référence montrent que LSF atteint un temps de résolution environ un ordre de grandeur plus rapide que les optimiseurs standards de type région de confiance et de gradient conjugué, particulièrement pour les systèmes avec $N > 20$ ions.
• Estimation de la Puissance de Commande : Les auteurs dérivent un estimateur heuristique de la fréquence de Rabi totale ($\Omega_{nuc}$) requise pour piloter une porte. Cet estimateur dépend de la norme nucléaire de la matrice d'intrication cible $\|\phi\|_{nuc}$ et évolue linéairement avec le nombre d'ions $N$. Remarquablement, les résultats de simulation montrent que la puissance de commande requise est largement indépendante de la structure spécifique des modes ou de l'espacement des fréquences, se superposant sur une courbe universelle lorsqu'elle est normalisée par l'échelle de temps de la densité moyenne des modes $T_\nu$.
• Temps de Porte Minimal ($T_{min}$) : L'étude identifie une limite inférieure fondamentale sur la durée de la porte. En dessous d'un temps critique $T_{min} \approx 2.9 \cdot T_\nu$ (où $T_\nu$ est proportionnel à l'inverse de l'espacement moyen des modes), aucune solution de haute fidélité ne peut être trouvée. Crucialement, $T_{min}$ évolue linéairement avec le nombre d'ions $N$. Cela implique que bien que le nombre d'opérations d'intrication simultanées évolue en $N^2$, la durée de la porte n'évolue que comme $N$, offrant un avantage d'évolutivité par rapport aux implémentations de portes séquentielles.
• Analyse de l'Évolution des Erreurs : L'article quantifie l'impact de trois sources d'erreurs dominantes sur la fidélité de la porte :
  • Dérives de fréquence des modes de mouvement : Les erreurs évoluent de manière quadratique avec la dérive de fréquence et le nombre de couplages ($N^2$).
  • Fluctuations de l'amplitude de commande : Les erreurs évoluent avec la variance du bruit et le nombre de couplages.
  • Chauffage motionnel : L'évolution des erreurs induites par le chauffage dépend de la corrélation spatiale du bruit. Pour un bruit non corrélé, la fidélité décroît de manière quadratique avec $N$ ; pour un bruit fortement corrélé, elle décroît linéairement.
    Le cadre LSF permet l'inclusion de contraintes linéaires pour rendre les portes robustes contre ces erreurs, réduisant considérablement la sensibilité aux dérives et au chauffage.
• Exemple de Code de Surface : La méthode est démontrée sur un cristal de 49 ions pour implémenter une porte de stabilisateur de code de surface (intriquant 33 ions en une seule impulsion). La solution montre une distribution efficace de la puissance, assignant des amplitudes de commande plus élevées aux ions « nœuds » et des amplitudes plus faibles aux ions « bords », et supprime efficacement l'intrication croisée indésirable entre des paires déconnectées.

Signification et Revendications
L'article affirme tracer une voie vers la mise à l'échelle des ordinateurs quantiques à ions piégés vers des centaines de qubits en surmontant la barrière conceptuelle de la conception de portes multi-qubits. En fournissant un algorithme rapide et robuste (LSF), les auteurs permettent la compilation hors ligne de circuits impliquant des dizaines d'ions et des portes multi-qubits complexes, une tâche auparavant considérée comme redoutable en raison de la nature NP-difficile de l'optimisation.

Ce travail établit que les portes multi-qubits de haute fidélité sont physiquement réalisables pour de grands cristaux, à condition que la durée de la porte dépasse le seuil évoluant linéairement avec $T_{min}$. De plus, les relations dérivées entre la puissance de commande, le temps de porte et les sources d'erreur fournissent des directives de conception et des tolérances pratiques pour les expérimentateurs. Les auteurs notent que bien que leur méthode traite les mécanismes d'erreur dominants dans l'approximation de Lamb-Dicke, les interactions avec les porteurs de qubits et les termes d'ordre supérieur dans l'expansion de Magnus pour les solutions à haute amplitude nécessitent une vérification supplémentaire. Le cadre est présenté comme un outil fondamental pour réaliser le potentiel de la connectivité longue distance dans les architectures à ions piégés pour la correction d'erreurs quantiques et les applications d'optimisation.