Quantum Simulation of Two-Level $PT$-Symmetric Systems Using Hermitian Hamiltonians

Résumé technique : Simulation quantique de systèmes à symétrie PT à deux niveaux utilisant des Hamiltoniens hermitiens

Énoncé du problème
La simulation de systèmes quantiques ouverts (SQO) sur des ordinateurs quantiques à portes unitaires présente un défi fondamental car les SQO présentent souvent des propagateurs temporels non unitaires, qui ne peuvent pas être directement mappés sur le cadre unitaire du matériel quantique actuel. Les systèmes à symétrie PT (Parité-Temps) constituent une classe spécifique de SQO caractérisés par des Hamiltoniens non hermitiens qui, dans leur phase « non brisée », possèdent des spectres d'énergie réels. Bien que ces systèmes offrent des propriétés uniques pour la détection quantique et la génération d'intrication — telles qu'une évolution plus rapide entre les états par rapport aux systèmes hermitiens — leur dynamique non unitaire complique la simulation directe sur les processeurs quantiques. Les méthodes existantes pour la simulation de systèmes ouverts, incluant l'évolution en temps imaginaire et l'encodage par blocs, nécessitent des approches efficaces spécifiquement adaptées à la dynamique PT-symétrique pour exploiter leurs avantages uniques.

Méthodologie
Les auteurs proposent une méthode pour simuler la dynamique PT-symétrique en exploitant des transformations de similitude pour mapper des Hamiltoniens PT-symétriques non hermitiens sur des Hamiltoniens hermitiens équivalents. Ce mappage est valide dans le régime de symétrie PT non brisé où l'Hamiltonien est pseudo-hermitien. Le cœur de l'approche repose sur une matrice de transformation de similitude, $\tau = \sqrt{\eta}$, qui relie l'Hamiltonien PT-symétrique non hermitien ($\hat{H}_{PT}$) à un équivalent hermitien ($\hat{h}$). L'évolution temporelle est ensuite exprimée par :
$$e^{-i\hat{H}_{PT}t} = \tau^{-1} e^{-i\hat{h}t} \tau$$
Cette formulation permet de décomposer la simulation en trois étapes : l'initialisation via $\tau$, l'évolution unitaire via $e^{-i\hat{h}t}$, et une transformation finale via $\tau^{-1}$.

Le document introduit deux algorithmes distincts pour implémenter ce cadre :

1. Algorithme hybride classique-quantique : Cette approche effectue l'initialisation et l'évolution temporelle unitaire sur un seul qubit. La dernière étape, l'application de $\tau^{-1}$, est gérée par un post-traitement classique. Cela nécessite une tomographie complète de l'état quantique pour reconstruire la matrice densité avant que la multiplication classique par $\tau^{-1}$ ne soit effectuée.
2. Algorithme de dilatation : Pour éviter les limitations de passage à l'échelle de la tomographie complète, cette méthode implémente $\tau^{-1}$ comme un opérateur unitaire dilaté à l'aide d'un qubit ancillaire. En utilisant la décomposition spectrale de la matrice hermitienne $\tau$, l'opérateur non unitaire est intégré dans un espace unitaire plus large. L'évolution est conditionnée par la mesure du qubit ancillaire dans l'état $|0\rangle$. Cette approche élimine le besoin de une tomographie complète de la matrice densité.

Pour étendre ces méthodes à des systèmes plus larges, les auteurs appliquent la théorie des perturbations du premier ordre pour simuler deux systèmes PT-symétriques faiblement interagissants. Ils construisent une nouvelle transformation de similitude basée sur les états propres perturbés du système couplé, permettant de mapper la dynamique sur une évolution hermitienne adaptée à l'algorithme hybride.

Résultats clés
Les algorithmes ont été validés par des simulations et des expériences sur l'unité de traitement quantique (QPU) IBM-Sherbrooke et des simulateurs avec bruit :

• Dynamique d'un qubit unique : Les deux algorithmes ont reproduit avec succès la dynamique de population des états fondamental et excité d'un système PT-symétrique proche d'un point exceptionnel. Les données expérimentales (cercles) et les simulations de modèles de bruit (lignes pointillées) montrent un fort accord avec les solutions classiques exactes (lignes pleines).
• Intrication de deux qubits : En utilisant l'algorithme hybride et la théorie des perturbations, les auteurs ont simulé deux qubits PT-symétriques faiblement couplés. Ils ont démontré avec succès la génération d'intrication, quantifiée par la concurrence, qui atteint un maximum lorsque les populations des états de base deviennent approximativement égales.
• Analyse d'erreur : L'étude a analysé la fidélité des deux algorithmes sous diverses erreurs de mesure.
  • L'algorithme hybride s'est avéré très sensible aux erreurs de mesure, particulièrement dans le régime d'évolution rapide où les erreurs sont amplifiées par l'étape de post-traitement classique. La fidélité chute significativement (jusqu'à ~20 %) à des taux d'erreur élevés.
  • L'algorithme de dilatation, bien que souffrant d'une augmentation du volume de circuit et d'erreurs de portes à deux qubits, a démontré une plus grande robustesse face aux erreurs de mesure, maintenant une fidélité plus élevée (~70 %) à des taux d'erreur élevés par rapport à l'approche hybride. Cependant, il reste sensible aux limitations de cohérence inhérentes aux portes à deux qubits.

Signification et revendications
L'article affirme offrir une approche évolutive pour simuler des systèmes PT-symétriques et des opérateurs pseudo-hermitiens sur des ordinateurs quantiques à portes unitaires. En construisant des Hamiltoniens hermitiens équivalents, les auteurs démontrent que la dynamique PT-symétrique non unitaire peut être simulée en utilisant l'évolution en temps réel plutôt que des portes non unitaires individuelles.

Le travail souligne que :

• Les systèmes PT-symétriques peuvent être simulés sur des dispositifs quantiques actuels de type NISQ (Noisy Intermediate-Scale Quantum) en utilisant soit des approches hybrides, soit des approches entièrement quantiques (dilatation).
• Ces méthodes peuvent être étendues à des systèmes multi-qubits faiblement couplés via la théorie des perturbations, permettant l'étude de la génération d'intrication dans les systèmes non hermitiens.
• Les caractéristiques d'erreur uniques des simulations PT-symétriques — spécifiquement l'amplification des erreurs dans les régimes d'évolution rapide — reflètent le comportement physique unique de ces systèmes près des points exceptionnels, suggérant une utilité potentielle pour les applications de détection quantique.

Les auteurs concluent que ce travail ouvre une voie pour de futures investigations sur les Hamiltoniens non hermitiens sur les ordinateurs quantiques, fournissant une base pour des travaux futurs visant à améliorer l'efficacité algorithmique et à explorer les applications pratiques de ces systèmes non hermitiens.