Reducing quantum resources for ADAPT-VQE via plateau-operator elimination and correlated mean-field downfolding

Résumé Technique : Réduction des ressources quantiques pour l'ADAPT-VQE via l'élimination d'opérateurs de plateau et le downfolding à champ moyen corrélé

Énoncé du problème
Les algorithmes ADAPT-VQE (Adaptive Derivative-Assembled Problem-Tailored Variational Quantum Eigensolvers) constituent une approche de pointe pour la chimie quantique sur les dispositifs quantiques de l'ère NISQ, offrant un moyen de construire des ansatz compacts et adaptés au système. Cependant, la méthode est confrontée à deux limitations majeures :

1. Opérateurs redondants et stagnation de l'optimisation : Le processus d'optimisation peut stagner en raison de « plateaux stériles » (barren plateaux) ou de l'inclusion d'opérateurs redondants. Plus précisément, des opérateurs avec des paramètres quasi nuls (opérateurs de plateau) peuvent être sélectionnés, ajoutant une profondeur de circuit et des portes CNOT inutiles sans contribuer à la réduction de l'énergie. Cela conduit à une convergence lente et à une inefficacité computationnelle.
2. Négligence de la corrélation dynamique externe : L'ADAPT-VQE standard opère généralement dans un espace actif restreint. Bien que cela réduise les besoins en qubits, cela néglige les effets de corrélation dynamique provenant de l'espace orbital externe (inactif), ce qui entraîne souvent des énergies s'écartant significativement des références FCI (Full Configuration Interaction), particulièrement près des géométries d'équilibre.

Méthodologie
Les auteurs proposent une stratégie à deux volets pour répondre à ces problèmes :

1. Élimination des opérateurs de plateau :

   • Les auteurs introduisent un schéma d'élimination directe intégré au flux de travail ADAPT-VQE. Après l'optimisation classique à chaque itération, les opérateurs dont les angles de rotation $|\theta_i|$ sont inférieurs à un seuil prescrit $\delta$ (proche de zéro) sont identifiés comme des « opérateurs de plateau ».
   • Ces opérateurs sont définitivement retirés du pool d'opérateurs. Les auteurs comparent deux variantes : une où le pool est élagué de façon permanente et une autre où le pool est restauré une fois que la fraction d'opérateurs éliminés dépasse un seuil $\phi$ (par exemple, 20–25 %).
   • Cette approche restreint la croissance des opérateurs à un sous-ensemble actif, réduisant ainsi le nombre d'évaluations de gradient par itération et empêant l'accumulation de portes non contributives.

2. Downfolding à un corps (OBDF - One-Body Downfolding) :

   • Pour capturer les corrélations entre l'espace actif et l'environnement, les auteurs combinent l'ADAPT-VQE avec leur cadre de One-Body Downfolding (OBDF) précédemment développé, basé sur la théorie de la perturbation de second ordre de Møller–Plesset à un corps (OBMP2).
   • Dans ce cadre, l'Hamiltonien moléculaire brut dans l'espace actif est remplacé par un Hamiltonien effectif corrélé ($\hat{H}_{OBDF}$). Cet Hamiltonien effectif incorpore les effets de corrélation dynamique de l'espace externe via un potentiel à un corps corrélé ($\hat{v}^{ext}_{OBMP2}$).
   • Crucialement, parce que la correction de downfolding est purement à un corps, l'Hamiltonien effectif résultant conserve la même structure d'opérateurs (termes à deux corps) que l'Hamiltonien de l'espace actif brut. Cela permet à la méthode d'utiliser les mêmes ressources quantiques (qubits et profondeur de circuit) que le problème standard tout en encodant l'information de corrélation externe.

Contributions clés

• Stratégie d'élimination d'opérateur : Une méthode directe et simple pour supprimer les opérateurs non contributifs du pool ADAPT-VQE, accélérant la convergence et réduisant la profondeur du circuit sans nécessiter de mécanismes de notation complexes comme ceux utilisés dans le Pruned-ADAPT-VQE ou le Param-ADAPT-VQE.
• Intégration de l'OBDF avec l'ADAPT-VQE : La première application du cadre OBDF à l'ADAPT-VQE, créant l'« OBDF-ADAPT-VQE ». Cela permet de simuler des systèmes plus larges ou d'obtenir des résultats plus précis dans le même espace actif en intégrant efficacement les corrélations dynamiques des orbitales externes dans l'Hamiltonien actif.
• Évaluation comparative (Benchmarking) : Évaluations exhaustives sur des chaînes $H_6$ linéaires, des réseaux $H_6$ et des anneaux $H_6$, ainsi que sur la molécule $N_2$ à travers diverses géométries en utilisant le simulateur OpenFermion.

Résultats

• Convergence et profondeur de circuit : Le schéma d'élimination d'opérateurs (spécifiquement la variante sans restauration du pool) améliore significativement le comportement de convergence. Dans les systèmes où l'ADAPT-VQE standard stagne (par exemple, la chaîne $H_6$ à certaines longueurs de liaison), l'algorithme modifié converge de manière fluide.
  • Le nombre de portes CNOT est réduit d'un facteur d'environ 2 à 3 dans tous les systèmes testés. Par exemple, les comptes de portes pour les chaînes, anneaux et réseaux $H_6$ ont été réduits respectivement en dessous de 4000, 3000 et 2000, contre 6000–8000 pour l'approche standard.
  • L'élimination des opérateurs de plateau ne compromet pas l'énergie finale convergée ; dans certains cas, elle produit des énergies légèrement plus basses en ouvrant des directions variationnelles auparavant inaccessibles à cause des opérateurs redondants.
• Performance de la base orbitale : L'utilisation des orbitales OBMP2 comme base pour l'ADAPT-VQE (au lieu des orbitales Hartree-Fock) produit un ansatz plus compact, particulièrement dans la région de dissociation, nécessitant moins d'opérateurs et de portes CNOT par rapport à la base HF.
• Précision via le Downfolding : L'approche OBDF-ADAPT-VQE produit systématiquement des énergies plus proches de la référence FCI que l'ADAPT-VQE standard au sein du même espace actif.
  • Pour les systèmes $H_6$, l'OBDF-ADAPT-VQE corrige les écarts importants près des géométries d'équilibre où la corrélation dynamique est critique.
  • Pour la molécule $N_2$, la méthode abaisse l'énergie pour se rapprocher des résultats CCSD(T) pour des distances inférieures à 1,6 Å et évite la divergence observée avec le CCSD(T) aux géométries étirées, tout en maintenant une forme de courbe de dissociation correcte.

Signification et affirmations
L'article affirme que les stratégies proposées offrent une voie pratique pour améliorer l'efficacité et la précision des simulations quantiques variationnelles sur les dispositifs de l'ère NISQ.

• Le schéma d'élimination d'opérateurs est présenté comme une « modification strictement non dégradante » qui accélère l'optimisation et réduit la surcharge de ressources sans pénalité énergétique.
• La combinaison avec l'OBDF est mise en avant comme une méthode permettant de récupérer une fraction substantielle de la contribution énergétique des orbitales exclues à un coût de calcul modeste, transférant efficacement l'information des orbitales externes inactives vers le calcul de l'espace actif.
• Les auteurs concluent que, bien que la méthode d'élimination directe soit moins sophistiquée que certaines techniques d'élagage existantes, elle donne des résultats prometteurs pour une amélioration pratique. L'intégration du downfolding à un corps montre qu'il est possible d'améliorer systématiquement la précision des surfaces d'énergie potentielle sans nécessiter d'expansion explicite de l'espace orbital actif, ce qui est particulièrement bénéfique lorsque les ressources quantiques sont limitées.