Nuclear scattering phase shifts with factorized geometric-time RODEO on a quantum processor

Résumé technique : Déphasages de diffusion nucléaire par la méthode RODEO à temps géométrique factorisé sur un processeur quantique

Énoncé du problème
La reconstruction des déphasages de diffusion nucléaire à partir de spectres calculés par voie quantique représente un défi important. Bien que les solveurs quantiques d'états propres puissent estimer des niveaux d'énergie discrets dans un système confiné (un « piège »), relier ces états discrets aux observables du continuum (diffusion en espace libre) nécessite un post-traitement classique, plus précisément l'expansion de portée effective modifiée (MERE). La performance d'un circuit quantique dans ce contexte ne peut être jugée uniquement par l'erreur quadratique moyenne (RMS) des niveaux d'énergie individuels. Les erreurs spectrales peuvent être amplifiées, supprimées ou réorganisées lors de l'interpolation d'énergie et de l'extrapolation à confinement nul subséquentes. Par conséquent, un circuit quantique à ressources réduites doit être évalué par sa capacité à préserver l'information nécessaire à la reconstruction finale de la diffusion, plutôt que par la seule fidélité spectrale.

Méthodologie
Les auteurs reconstruisent les déphasages de diffusion élastique s-onde neutron-proton (np) en utilisant des données du processeur quantique IBM Aachen. L'étude compare deux implémentations de l'algorithme RODEO (Repetitive Operation for Discrete Energy Optimization) :

1. Direct R10 Dynamique : Une implémentation de référence utilisant dix cycles avec des pas de temps à quantiles gaussiens fixes. Pour chaque énergie de requête, l'évolution contrôlée est resynthétisée et compilée séparément.
2. RODEO à temps géométrique factorisé (FG–R6) : Une implémentation compressée utilisant six cycles avec des pas de temps échantillonnés géométriquement. Cette méthode emploie une séparation exacte de la phase de requête, où l'évolution hamiltonienne est fixée à une énergie de référence ($E_{ref}$) et compilée une seule fois. La phase dépendante de la requête ($\Delta E = E_q - E_{ref}$) est liée au qubit ancilla après la compilation.

Les deux méthodes utilisent un registre de trois qubits (deux qubits de système et un ancilla recyclé) sur le backend IBM Aachen. Le modèle physique est une interaction de type puits carré s-onde np schématique. Des hamiltoniens effectifs d'Okubo–Lee–Suzuki (OLS) en quatre dimensions, construits classiquement, sont utilisés pour représenter le système. L'expérience mesure quatre niveaux d'énergie positifs à cinq intensités de piège différentes ($\omega = 3,6$ à $4,0$ MeV), soit vingt entrées spectrales. Ces entrées sont injectées dans une analyse MERE classique pour l'extrapolation vers l'espace libre et le calcul des déphasages sur une grille d'énergie de 0,1 à 30,0 MeV.

Contributions clés

• Compression de circuit via la factorisation : L'article introduit le circuit FG–R6, qui combine la réutilisation de l'ancilla avec une stratégie d'évolution temporelle factorisée. En séparant la phase de requête de l'évolution hamiltonienne, le circuit réduit le nombre de cycles d'évolution contrôlée de dix à six et élimine la nécessité de recompiler l'opérateur d'évolution pour chaque énergie de requête.
• Évaluation matérielle (Benchmarking) : L'étude fournit un benchmark complet dérivé du matériel pour la diffusion nucléaire, allant au-delà des tests de niveau unique pour englober un ensemble complet de vingt entrées spectrales nécessaires à la reconstruction en espace libre.
• Analyse du compromis Ressources vs Précision : Le travail démontre que la réduction de la profondeur du circuit et du nombre de portes à deux qubits ne dégrade pas nécessairement la précision de l'observable de diffusion finale, même si l'erreur RMS des niveaux d'énergie intermédiaires augmente.

Résultats

• Réduction des ressources : Par rapport à l'implémentation Direct R10, l'implémentation FG–R6 réduit la profondeur médiane du circuit compilé de 33,1 % (de 459 à 307) et le nombre de portes à deux qubits de 37,2 % (de 113 à 71).
• Incertitude d'énergie : FG–R6 atteint une incertitude d'énergie médiane à tirage fini plus faible (0,534 keV) que R10 (0,610 keV), malgré une amplitude ajustée médiane et un décalage de centre RMS ($\Delta E_c$) plus importants par rapport à l'énergie de référence (1,648 keV contre 0,793 keV).
• Précision du déphasage de diffusion : Sur l'ensemble de la grille 0,1–30,0 MeV, la déviation maximale du déphasage de diffusion reconstruit par rapport à la référence MERE d'énergie exacte diminue de manière significative, passant de 0,842° (R10) à 0,163° (FG–R6). La déviation maximale par rapport à la solution analytique du puits carré chute également de 0,981° à 0,156°.
• Performance dépendante de l'énergie : L'amélioration est principalement pilotée par la région d'extrapolation à basse énergie (inférieure à 10 MeV). Dans le sous-ensemble d'interpolation commun 10,0–30,0 MeV, l'implémentation Direct R10 produit en fait des résidus légèrement plus faibles que FG–R6.
• Sensibilité de la forme de l'ajustement (Fit-Form) : L'étude note que la sensibilité à la forme de l'ajustement reste appréciable. L'écart des résultats entre les différents ordres polynomiaux pour les ajustements d'énergie et de piège est substantiel (par exemple, une demi-plage de 20,12° sur la grille complète pour R10 contre 6,75° pour FG–R6), indiquant que le choix du modèle d'extrapolation impacte significativement le résultat final.

Signification et affirmations
L'article affirme que ce benchmark à espace réduit réussit à connecter la compression de circuit RODEO aux observables de la diffusion nucléaire. La signification primaire réside dans la démonstration que la performance d'un circuit doit être évaluée par la reconstruction de l'observable de diffusion et de sa plage d'énergie spécifique, plutôt que par les seules erreurs RMS spectrales. L'implémentation FG–R6 est présentée comme un compromis précision-ressource validé pour le spectre et la fenêtre testés, plutôt que comme une solution universelle à nombre de cycles minimal.

Les auteurs déclarent explicitement que ces résultats ne démontrent pas d'avantage quantique ni de scalabilité multi-corps. Les données matérielles ne déterminent pas séparément les effets de la compilation de phase par rapport à l'échantillonnage temporel, et n'établissent pas non plus la reproductibilité des exécutions indépendantes. Ce travail sert de point de départ pour étendre ces méthodes à des espaces de modèles plus larges, des interactions réalistes et la dynamique de réaction de quelques corps, telles que la diffusion résonante à ondes partielles supérieures (par exemple, la diffusion neutron-alpha) et les systèmes impliquant des interactions Coulombiennes à longue portée.