Solving Einstein Field Equations on a Digital Quantum Computer

Énoncé du Problème
La relativité numérique (NR) implique la résolution des équations du champ d'Einstein (EFE), un système d'équations aux dérivées partielles (EDP) hautement non linéaires, afin de simuler des systèmes gravitationnels relativistes. Actuellement, ces simulations sont extrêmement coûteuses en termes de calcul, nécessitant des supercalculateurs et des clusters pour discrétiser et résoudre les équations de manière récursive. Bien que les algorithmes quantiques aient montré des promesses pour la résolution d'EDP dans d'autres domaines (par exemple, la dynamique des fluides), ils n'ont pas encore été appliqués au problème spécifique de la résolution numérique des EFE dans le cadre de la relativité numérique. Le défi réside dans la transposition des opérations classiques de la NR, souvent complexes, couplées et non unitaires (telles que la copie d'états ou l'application de fonctions non linéaires), vers les opérations unitaires et réversibles requises par les ordinateurs quantiques numériques.

Méthodologie
Les auteurs proposent un algorithme quantique de preuve de concept pour résoudre les EFE en utilisant le formalisme de la tétrade hyperbolique de Wahlquist-Estabrook-Buchman-Bardeen (WEBB). Ce formalisme a été choisi de préférence aux approches basées sur la métrique (comme ADM ou BSSN) car il fait évoluer les coefficients de connexion (coefficients de rotation de Ricci) qui se transforment comme des champs scalaires, permettant un codage d'amplitude compact dans une structure vectorielle plutôt qu'en tenseur.

La stratégie centrale consiste à mapper les équations d'évolution de WEBB discrétisées vers un cadre de marche quantique (Quantum Walk) via la simulation hamiltonienne. Les étapes spécifiques incluent :

1. Sélection du formalisme : Utilisation du formalisme WEBB, qui produit un système hyperbolique symétrisable de premier ordre. Les équations d'évolution prennent la forme $D_0 q + M^a D_a q = S$, où $q$ représente le vecteur d'état des coefficients de connexion, $M^a$ sont des matrices unitaires creuses, et $S$ contient des termes sources non linéaires.
2. Mapping par opérateur scindé (Split-Operator Mapping) : Les auteurs mappent les EDP discrétisées vers un algorithme de marche quantique, de manière analogue à la méthode de l'opérateur scindé utilisée pour l'équation de Dirac. Cela permet de décomposer l'évolution en une étape d'advection (propagation linéaire) et une étape de collision (source non linéaire).
3. Composants du circuit quantique :
   • Préparation de l'état : Les champs sont encodés en amplitude dans un état spinorique $|\Psi\rangle = \sum \psi_c(x) |c\rangle|x\rangle$, où $|c\rangle$ encode les composantes du champ (registre de la pièce/coin) et $|x\rangle$ encode les points de la grille spatiale (registre de la grille).
   • Étape d'advection : Implémentée via des opérateurs de décalage contrôlés sur le registre de la grille, conditionnés par le registre de la pièce, simulant efficacement la propagation linéaire de la fonction d'onde.
   • Étape de collision : Le terme source non linéaire $S$ est traité à l'aide d'un opérateur de collision heuristique ou d'une linéarisation basée sur le Jacobien. Pour garantir l'unitariété, les opérateurs non unitaires sont intégrés dans une matrice hermitienne plus large agissant sur un qubit ancilla, puis exponentiés sous la forme d'une porte hamiltonienne ($e^{-iHt}$).
   • Lapse et Shift (Relâchement et Décalage) : Les fonctions de lapse et de shift statiques sont implémentées en les encodant en amplitude sur des qubits ancillas et en conditionnant les opérations d'advection/collision sur ces ancillas.
4. Conditions aux limites : L'algorithme utilise des limites d'excision aux horizons de trous noirs (où aucune information ne s'échappe) et des conditions absorbantes à la limite extérieure pour éviter les réflexions parasites, adaptées à la nature périodique de la grille par défaut de Qiskit.

Principales Contributions

• Premier algorithme quantique numérique pour la relativité numérique : Le papier présente la première tentative de numérisation de la résolution des équations du champ d'Einstein spécifiquement pour la relativité numérique sur un ordinateur quantique à portes.
• Adaptation du formalisme WEBB : Les auteurs démontrent que le formalisme de la tétrade WEBB est particulièrement bien adapté à une implémentation quantique en raison de sa structure hyperbolique de premier ordre et de ses variables d'évolution de type scalaire.
• Implémentation algorithmique : Le travail détaille les composants spécifiques du circuit quantique requis pour résoudre ces équations, incluant la préparation d'état, l'advection contrôlée et la gestion des termes sources non linéaires via la simulation hamiltonienne assistée par ancilla.
• Validation physique : L'algorithme est testé en faisant évoluer un espace-temps de trou noir de Schwarzschild et en le perturbant pour générer des modes quasi-normaux gravitationnels (QNM). Les auteurs extraient le scalaire de Weyl $\Psi_4$ pour analyser les fréquences de résonance du trou noir perturbé.

Résultats
Les auteurs ont implémenté les composants de l'algorithme en utilisant le framework logiciel Qiskit. Ils ont exécuté des simulations sur :

• Des simulateurs d'état vectoriel classiques.
• Des ordinateurs quantiques physiques IBM (accessibles via le centre national britannique de calcul quantique UKRI).

Les résultats ont démontré avec succès l'évolution de la solution de Schwarzschild et l'extraction des modes quasi-normaux. Le signal extrait de la simulation quantique (spécifiquement la phase de ringdown) a été comparé à des solutions classiques par différences finies, montrant que l'algorithme quantique pouvait reproduire le comportement physique attendu (la perturbation voyageant vers l'horizon et les QNM voyageant vers l'extérieur). Les auteurs ont quantifié les ressources de calcul et le temps d'exécution, notant les contraintes imposées par le matériel actuel de type NISQ (Noisy Intermediate Scale Quantum), telles que la décohérence et le besoin de mitigation d'erreurs.

Signification et Revendications
Le papier cadre explicitement ce travail comme un exercice de preuve de concept. Les auteurs ne prétendent pas que leur méthode actuelle offre un « avantage quantique » (c'est-à-dire une accélération par rapport aux supercalculateurs classiques) en termes de temps d'exécution ou d'efficacité des ressources. Au contraire, la signification réside dans la démonstration que les équations complexes de la relativité générale peuvent être implémentées avec succès dans un cadre de calcul quantique numérique.

Ce travail sert de étape fondamentale pour les communautés de la relativité numérique et de l'informatique quantique, prouvant que :

1. Les EFE peuvent être numérisées en utilisant des techniques de simulation hamiltonienne.
2. Le formalisme WEBB offre une voie viable pour une telle numérisation.
3. Des phénomènes physiques tels que les perturbations de trous noirs et les QNM peuvent être simulés sur le matériel quantique bruyant actuel.

Les auteurs concluent que, bien que les dispositifs NISQ actuels limitent l'échelle et la précision de ces simulations, les composants algorithmiques développés et la cartographie réussie du problème fournissent une boîte à outils utile pour les améliorations futures, le développement de la correction d'erreurs et le potentiel de mise à l'échelle à mesure que le matériel quantique mûrira.