Symmetry-adapted qubit encoding with complete active space and Bravyi--Kitaev mapping for quantum chemistry on a quantum computer

Cet article introduit un codage de qubits adapté à la symétrie avec espace actif complet (SAE-CAS) qui intègre des symétries Z approximatives et le mappage de Bravyi-Kitaev afin de réduire significativement le nombre de qubits et la complexité des circuits pour les simulations de chimie quantique, démontrant une convergence et une efficacité de ressources supérieures aux méthodes standards sur des processeurs quantiques tant de l'ère NISQ que tolérants aux fautes.

L'essentiel

Imaginez que vous essayez de résoudre un puzzle géant et incroyablement complexe pour comprendre comment fonctionne une molécule. Dans le monde de la chimie quantique, ce puzzle est la « structure électronique » d'une molécule. Pour le résoudre sur un ordinateur quantique, nous devons généralement assigner une minuscule partie de l'ordinateur (un « qubit ») à chaque emplacement possible où un électron pourrait se trouver.

Le problème ? Même pour de petites molécules, cela nécessite des milliers de pièces de puzzle (des qubits), et les instructions pour les assembler (le circuit) deviennent si longues et emmêlées que les ordinateurs actuels ne peuvent pas les gérer, et même les futurs ordinateurs auront du mal.

Ce document présente une nouvelle méthode ingénieuse pour résoudre le puzzle appelée SAE-CAS. Voici comment elle fonctionne, en utilisant des analogies simples :

1. La stratégie du « Geler et Ignorer » (La partie CAS)

Considérez une molécule comme un immeuble de bureaux très animé.

• Le Noyau Gelé : Le sous-sol et le dernier étage sont toujours occupés par des gens qui ne partent jamais et n'interagissent jamais avec le reste du bâtiment. En termes quantiques, ce sont des électrons de « cœur gelé » (frozen-core). Ils sont ennuyeux et prévisibles.
• Les Virtuels : Le grenier est complètement vide et restera probablement vide. Ce sont des orbitales « virtuelles ».
• L'Espace Actif : Les étages du milieu sont là où l'action se déroule réellement. Les gens circulent, discutent et changent les choses. C'est l'« espace actif ».

Les méthodes traditionnelles tentent d'assigner un qubit à chaque étage, même au sous-sol ennuyeux et au grenier vide. SAE-CAS dit : « Ignorons simplement le sous-sol et le grenier. » Nous n'assignons des qubits qu'aux étages du milieu, là où la chimie intéressante se produit. Cela réduit immédiatement la taille du puzzle que nous devons résoudre.

2. Le « Raccourci de Symétrie » (La partie SAE)

Même au sein des étages du milieu très animés, il existe des règles. Par exemple, dans une molécule d'eau, le côté gauche est l'image miroir du côté droit. Si vous savez ce qui se passe à gauche, vous savez automatiquement ce qui se passe à droite.

Habituellement, les ordinateurs calculent les deux côtés séparément, gaspillant ainsi du temps et des ressources. SAE-CAS utilise un « tour de magie » mathématique (appelé transformation de Clifford affine) pour réaliser que, grâce à ces règles de miroir, nous n'avons pas besoin de qubits distincts pour les deux côtés. Nous pouvons mathématiquement « plier » le puzzle en deux. Cela supprime encore plus de qubits, rendant le puzzle plus petit et plus facile à résoudre.

3. La « Traduction » (La partie Bravyi–Kitaev)

Une fois que nous avons notre puzzle plié et miniature, nous devons le traduire dans un langage que l'ordinateur quantique comprend. Il existe deux traducteurs principaux :

• Jordan-Wigner (JW) : Le traducteur standard. Il est simple mais rend les instructions très longues (comme lire un livre où chaque mot est répété).
• Bravyi–Kitaev (BK) : Un traducteur plus intelligent. Il organise l'information de manière plus efficace, de sorte que les instructions soient plus courtes et moins emmêlées.

Les auteurs montrent que vous pouvez utiliser leur méthode de « Puzzle Plié » (SAE-CAS) avec soit l'un, soit l'autre traducteur. Ils ont créé une version appelée SAE-CAS-BK qui utilise le traducteur le plus intelligent. Cela ne change pas la réponse finale, mais cela rend souvent le chemin pour y parvenir plus fluide et plus rapide.

Qu'ont-ils trouvé ?

Les auteurs ont testé cette méthode sur neuf petites molécules (comme l'eau, l'oxygène et l'azote) en utilisant deux « stratégies de recherche » différentes (algorithmes) pour trouver l'énergie de la molécule :

1. UCCSD : Une recherche chimiquement précise mais complexe.
2. HE-SCA : Une recherche plus simple, mieux adaptée au matériel (hardware-friendly).

Les Résultats :

• Moins de Qubits : En ignorant les parties ennuyeuses et en pliant les parties symétriques, ils ont eu besoin de nettement moins de qubits (réduisant parfois le nombre de moitié ou plus).
• Circuits plus courts : Les instructions pour exécuter la simulation étaient beaucoup plus courtes et moins emmêlées.
• Succès plus rapide : En utilisant la stratégie de recherche plus simple (HE-SCA), leur méthode a trouvé la bonne réponse pour chaque molécule testée. L'ancienne méthode (JW-CAS) s'est bloquée et a échoué à trouver la réponse pour l'oxygène et le monoxyde de carbone dans les limites de temps imparties.
• Aucune perte de précision : Même en ignorant les « électrons ennuyeux » et en pliant le puzzle, les chiffres d'énergie finaux étaient tout aussi précis que les calculs massifs standards.

L'essentiel

Les auteurs ont construit une boîte à outils « efficace en ressources ». Ils ont prouvé que vous pouvez en toute sécurité écarter les parties de la molécule qui ne changent pas et plier les parties qui sont symétriques, sans perdre la réponse correcte. Cela rend possible l'exécution de ces simulations chimiques complexes sur des ordinateurs quantiques beaucoup plus petits et moins puissants que ce que l'on pensait nécessaire auparavant.

Ils ont également rendu le code de ce « tour de magie » disponible gratuitement (via un package appelé QuantumSymmetry) afin que d'autres puissent l'utiliser pour simuler des molécules sur leurs propres ordinateurs quantiques.

Résumé technique

Résumé technique : Codage de qubits adapté aux symétries avec espace actif complet et mappage de Bravyi–Kitaev

Énoncé du problème
Les simulations de chimie quantique sur les processeurs quantiques à court terme (NISQ) et tolérants aux fautes font face à des contraintes de ressources importantes, principalement dictées par le nombre de qubits requis pour encoder les structures électroniques moléculaires et la complexité (poids de Pauli et profondeur de circuit) des Hamiltoniens résultants. Bien que l'approximation de l'espace actif complet (CAS - Complete Active Space) réduise la taille du problème en figeant les orbitales du cœur et en écartant les orbitales virtuelles, les implémentations standards échouent souvent à exploiter pleinement les symétries pour réduire davantage le nombre de qubits. De plus, différents mappages fermions-vers-qubits, tels que Jordan–Wigner (JW) et Bravyi–Kitaev (BK), offrent des compromis entre localité et exigences de qubits qui ne sont pas toujours intégrés de manière optimale avec les techniques de réduction de symétrie.

Méthodologie
Les auteurs proposent un cadre unifié appelé Codage adapté aux symétries avec espace actif complet (SAE-CAS). Cette approche traite l'approximation CAS non pas simplement comme une troncature de l'espace de Hilbert, mais comme l'imposition de symétries $Z$ approximatives sur les orbitales du cœur figées et les orbitales virtuelles.

1. Cadre théorique :

   • Symétries approximatives : La méthode étend les codages adaptés aux symétries (SAE) exacts pour inclure des symétries approximatives correspondant aux occupancies orbitales fixes (cœur figé et virtuel). Celles-ci sont modélisées comme des symétries $Z$ où des qubits spécifiques sont contraints à des valeurs propres de $-1$ (figé) ou $+1$ (virtuel).
   • Transformations de Clifford affines : La réduction est réalisée via des applications de Clifford affines. Les auteurs définissent une matrice de caractère orbitalaire binaire $S$ représentant les symétries. En construisant une application affine invertible $T$ sur $\mathbb{F}_2$, ils transforment la base de calcul de sorte que les contraintes de symétrie correspondent à des valeurs fixes (typiquement 0) sur des qubits spécifiques.
   • Suppression de qubits : Une fois l'Hamiltonien conjugué par l'opérateur de Clifford, les termes contenant des opérateurs $X$ ou $Y$ sur les qubits contraints sont écartés, et les opérateurs $Z$ contraints sont remplacés par des scalaires. Cela supprime de façon permanente les qubits redondants, réduisant le système de $n$ orbitales de spin à $n-k$ qubits, où $k$ est le nombre de symétries indépendantes.
   • Preuve d'équivalence : Les auteurs prouvent que l'Hamiltonien de qubits de l'espace actif résultant est exactement équivalent à l'Hamiltonien CAS canonique dérivé de la projection des orbitales du cœur et des orbitales virtuelles dans la théorie des structures électroniques moléculaires. L'erreur d'approximation provient uniquement du choix de l'espace actif, et non du codage.

2. Intégration avec Bravyi–Kitaev (SAE-CAS-BK) :

   • Le cadre est étendu au mappage Bravyi–Kitaev (BK). Puisque la transformation BK est elle-même un changement de base de Clifford affine agissant sur les chaînes de bits de calcul, elle peut être composée avec l'application de Clifford SAE-CAS.
   • L'encodage SAE-CAS-BK résultant est unitairement équivalent à SAE-CAS. Il préserve le nombre de qubits, le nombre de termes de Pauli et les paramètres variationnels, mais modifie la structure de localité (poids de Pauli par terme) et la profondeur de circuit/le nombre de CNOT qui en résulte.

3. Évaluation numérique :

   • La méthode a été testée sur neuf petites molécules (incluant H$_2$O, C$_2$H$_4$, O$_2$, N$_2$ et C$_4$H$_6$) utilisant une base minimale STO-3G.
   • Deux ansätze variationnels ont été employés : l'Unitary Coupled Cluster avec excitations simples et doubles (UCCSD) et l'ansatz hardware-efficient shifted-circular-alternating (HE-SCA).
   • Des comparaisons ont été effectuées par rapport à JW standard, JW-CAS (CAS avec JW mais sans réduction de symétrie) et JW-CAS avec filtrage de symétrie (SF) au niveau du circuit.

Résultats clés

• Réduction des ressources : SAE-CAS réduit systématiquement le nombre de qubits, le nombre d'opérateurs de Pauli, la profondeur de circuit, le nombre de CNOT et les paramètres variationnels par rapport à JW-CAS. Par exemple, pour le CAS (4,4) de C$_4$H$_4$, SAE-CAS a réduit le nombre de qubits de 8 (JW-CAS) à 4 et le nombre de termes de Pauli de 97 à 46.
• Convergence et précision :
  • UCCSD : SAE-CAS a atteint des énergies convergées comparables à JW-CAS et JW-CAS (SF) mais a nécessité nettement moins d'itérations d'optimisation (par exemple, 21 contre 109 pour N$_2$). Notamment, pour C$_4$H$_4$ et C$_4$H$_6$, les ansätze UCCSD de JW-CAS et JW-CAS (SF) n'ont pas pu s'exécuter en raison de l'épuisement de la mémoire causé par le grand nombre d'excitations, alors que SAE-CAS a réussi en filtrant la liste des excitations lors de l'étape de mappage.
  • HE-SCA : SAE-CAS a démontré des propriétés de convergence supérieures. Pour toutes les molécules testées, SAE-CAS a convergé vers les énergies de référence CAS dans les budgets de couches testés. En revanche, JW-CAS n'a pas réussi à converger pour O$_2$ et CO même après 200 couches. Les auteurs attribuent cela à l'effet de "plateau stérile" (barren plateau) et à la fuite de secteur dans JW-CAS, où la variété variationnelle s'étend sur des secteurs de parité de spin incorrects, tandis que SAE-CAS restreint la recherche au bon secteur de symétrie de manière native.
  • SAE-CAS vs. SAE-CAS-BK : Les deux encodages ont produit des énergies convergées et des nombres d'itérations identiques, confirmant leur équivalence unitaire. Des différences ont été observées uniquement dans la profondeur de circuit et le nombre de CNOT, qui variaient de moins de 6 % pour la plupart des molécules (sauf CO, qui a vu une augmentation de +20 % pour BK). Les auteurs notent que l'avantage de localité $O(\log n)$ asymptotique de BK n'est pas encore réalisé à ces tailles d'espace actif réduites.

Signification et revendications
L'article affirme que SAE-CAS offre une voie systématique et composable vers des simulations moléculaires économes en ressources. En traitant l'approximation CAS comme un ensemble de symétries approximatives, la méthode unifie la réduction de symétrie exacte (groupe ponctuel et parité de spin) avec la troncature de l'espace actif. La principale signification réside dans :

1. Imposition native des symétries : L'imposition des symétries au niveau de l'encodage (plutement qu'a posteriori dans l'ansatz) élimine les directions de recherche redondantes et empêche l'optimiseur de dériver vers des secteurs incorrects de nombre de particules ou de spin.
2. Scalabilité : La méthode permet des simulations sur des espaces actifs plus larges qui seraient autrement impossibles à cause des contraintes de mémoire lors de la construction de l'ansatz (comme vu dans les cas C$_4$H$_4$ et C$_4$H$_6$).
3. Flexibilité : Le cadre est compatible avec tout mappage fermions-vers-qubits exprimable comme un changement de base de Clifford affine, comme démontré ici avec JW et BK.

Les auteurs fournissent une implémentation open-source dans le package Python QuantumSymmetry et concluent que SAE-CAS est une voie viable pour réduire le nombre de qubits et la complexité des opérateurs sur les processeurs quantiques tant à court terme que tolérants aux fautes.