Exact chiral symmetry with quantum signal processing

Résumé technique : Symétrie chirale exacte avec le traitement du signal quantique

Énoncé du problème
La simulation quantique d'observables non perturbatives et hors équilibre dans la QCD sur réseau (lattice QCD) fait face à des défis importants concernant le doublement des fermions et la préservation de la symétrie chirale dans les formulations hamiltoniennes. Bien que la relation de Ginsparg-Wilson (GW) et l'opérateur de recouvrement (overlap) offrent des solutions élégantes dans l'espace-temps euclidien, un analogue hamiltonien canonique pour les fermions de recouvrement faisait défaut. Les approches existantes, telles que les fermions à paroi de domaine (domain-wall fermions), se transposent bien aux cadres hamiltoniens mais nécessitent une dimension supplémentaire explicite, augmentant ainsi le coût en qubits. Inversement, les fermions de recouvrement offrent une formulation dans les dimensions spatiales physiques, mais impliquent des interactions non locales de type « tout-vers-tout » très complexes, spécifiquement la fonction signe de l'hamiltonien de Wilson-Dirac, $\varepsilon(h_W)$, qui est difficile à implémenter efficacement sur le matériel quantique. Le problème central abordé est la manière de construire des algorithmes quantiques efficaces pour les fermions sur réseau à symétrie chirale, en équilibrant les coûts de mémoire (qubits) par rapport à la complexité des portes, tout en maintenant une symétrie chirale exacte (ou contrôlablement brisée).

Méthodologie
Les auteurs proposent un algorithme de traitement du signal quantique (QSP) pour simuler l'hamiltonien de recouvrement. La méthodologie procède par étapes clés :

1. Formulation de l'hamiltonien : Ce travail utilise l'hamiltonien de recouvrement $H_{ov} = \psi^\dagger (\Gamma_0 + \varepsilon(h_W)) \psi$, où $h_W$ est l'hamiltonien de Dirac-Wilson à particule unique. La fonction signe $\varepsilon(X) = X/\sqrt{X^\dagger X}$ est approximée par un polynôme $p_M(x)$ de degré $M$.
2. Encodage par bloc (Block Encoding) : Les auteurs construisent un encodage par bloc de l'hamiltonien de Wilson à particule unique $\hat{h}_W$ (incluant les champs de jauge) en utilisant un opérateur de préparation $P$ et un opérateur de sélection $S$. Cela encode l'hamiltonien dans un opérateur unitaire $U_{\hat{h}_W}$ agissant sur un espace de Hilbert étendu avec des qubits ancillas.
3. Traitement du signal quantique (QSP) : Pour implémenter l'approximation de la fonction signe, les auteurs emploient la QSP (plus précisément la transformation de la valeur singulière quantique). Cela permet d'appliquer un polynôme de degré $M$, $p_M(U_{\hat{h}_W})$, à l'opérateur encodé par bloc. Le polynôme est choisi pour approximer la fonction signe dans un intervalle spectral $[\kappa, 1]$ avec une erreur $\epsilon_e$.
4. Évolution temporelle : Une fois l'hamiltonien de recouvrement encodé par bloc, l'opérateur d'évolution temporelle $e^{-iHt}$ est approximé en utilisant la QSP, nécessitant un nombre de portes proportionnel au coût de l'encodage par bloc et au temps d'évolution.

Contributions clés et résultats

• Construction algorithmique : Le papier fournit un algorithme concret basé sur la QSP pour l'hamiltonien de recouvrement qui préserve la relation de GW jusqu'à une erreur contrôlable $\epsilon_e$. L'opérateur chiral modifié $\hat{\gamma}_5$ est montré comme commutant avec l'hamiltonien approximatif à une erreur de l'ordre de $\epsilon_e$.
• Analyse de complexité :
  • Complexité des portes : Le coût pour encoder par bloc l'hamiltonien de recouvrement évolue en $O(Q \kappa^{-1} \log(1/\epsilon_e))$, où $Q$ est le nombre de sites du réseau (multiplié par les degrés de liberté internes) et $\kappa$ est lié au gap spectral. La complexité totale des portes pour l'évolution temporelle évolue en $O(Q \kappa^{-1} \log(1/\epsilon_e) [Qt + \log(1/\epsilon_t)])$.
  • Coût en qubits : L'algorithme nécessite $O(Q + \log Q)$ qubits. Il s'agit d'une réduction significative par rapport aux fermions à paroi de domaine, qui nécessitent $O(Q L_5)$ qubits, où $L_5$ est l'étendue de la cinquième dimension.
  • Comparaison : Bien que les fermions à paroi de domaine bénéficient d'une localité géométrique (permettant un coût quasi linéaire $O(Q t \text{polylog})$ via Suzuki-Trotter ou des méthodes similaires), la formulation de recouvrement entraîne une profondeur de portes plus élevée en raison de la nature non locale de l'approximation de la fonction signe. Cependant, l'approche de recouvrement offre une meilleure mise à l'échelle asymptotique en termes de mémoire.
• Interprétation physique de la QSP : Les auteurs démontrent que le degré du polynôme $M$ requis pour atteindre une erreur $\epsilon_e$ évolue comme $M = O(\kappa^{-1} \log(1/\epsilon_e))$. Ils identifient cette mise à l'échelle avec la taille de la dimension supplémentaire $L_5$ dans les fermions à paroi de domaine, où l'erreur évolue selon $e^{-c L_5}$. Ainsi, l'implémentation par QSP « construit » effectivement une dimension supplémentaire à travers la profondeur du circuit, reflétant la correspondance physique entre l'opérateur de recouvrement et la théorie de bord des fermions à paroi de domaine.

Signification et affirmations
Le papier affirme que les simulations quantiques de fermions de Dirac avec une symétrie chirale exacte sont « presque gratuites » dans le sens où l'application de l'hamiltonien de recouvrement ne coûte qu'un facteur logarithmique supplémentaire (en termes de tolérance d'erreur) par rapport à l'hamiltonien de Dirac-Wilson.

L'idée centrale est que la QSP fournit une réalisation algorithmique quantique de la correspondance connue entre l'opérateur de recouvrement et la dimension supplémentaire des fermions à paroi de domaine. Le compromis est explicite :

• Fermions à paroi de domaine : Coût en qubits plus élevé ($O(Q L_5)$) mais profondeur de portes plus faible grâce à la localité.
• Fermions de recouvrement (via QSP) : Coût en qubits plus faible ($O(Q)$) mais profondeur de portes plus élevée ($O(Q^2)$ dans le pire cas pour l'évolution temporelle en raison de la non-localité).

Les auteurs concluent que le choix entre les formulations dépend des contraintes matérielles (nombre de qubits vs profondeur de circuit). Ils notent que si la QCD sur réseau classique utilise souvent des approximations rationnelles (par exemple, Zolotarev) pour la fonction signe, ce qui est plus efficace, celles-ci n'ont pas d'analogues directs dans la QSP standard, suggérant une direction pour des travaux futurs. Le papier précise également que les résultats de mise à l'échelle sont robustes face à l'encodage spécifique des champs de jauge, bien que les facteurs constants puissent varier.