Polynomial-time thermalization and Gibbs sampling from system-bath couplings

Résumé Technique : Thermalisation en temps polynomial et échantillonnage de Gibbs à partir de couplages système-réservoir

Énoncé du problème
L'article traite de la vitesse de convergence de la dynamique de Lindblad, qui modélise les systèmes quantiques ouverts faiblement couplés à un réservoir thermique, vers leurs états de Gibbs stationnaires. Comprendre ces taux de convergence (temps de mélange) est crucial pour borner les temps d'exécution des algorithmes d'échantillonnage de Gibbs quantique et les échelles de temps de la thermalisation physique.

Deux familles spécifiques de processus sont étudiées :

1. Échantillonnage de Gibbs par interactions répétées : Un algorithme où un système subit des interactions faibles répétées avec des qubits thermiques ancillaires. Bien que conceptuellement simple et adapté à une mise en œuvre à court terme, les garanties de performance précédentes étaient limitées aux Hamiltoniens commutants.
2. Thermalisation de réservoir macroscopique : Un modèle dérivé des interactions faibles système-réservoir où le réservoir est un système large et non contrôlé. Les générateurs de Lindblad résultants satisfont la condition de l'équilibre détaillé de Kubo-Martin-Schwinger (KMS). Cependant, pour que l'approximation soit précise (couplage $\alpha$ faible), les opérateurs de saut deviennent hautement non locaux, rendant les preuves existantes sur les temps de mélange inapplicables.

Le défi technique central est que les preuves existantes pour les gaps spectraux (qui régissent la vitesse de convergence) de Lindbladiens symétriques KMS reposent sur le fait que les opérateurs de saut soient quasi-locaux. Dans les deux modèles ci-dessus, les régimes où l'approximation physique est précise (grand $\kappa$ pour les interactions répétées, faible $\alpha$ pour les réservoirs macroscopiques) produisent des générateurs non locaux, créant un écart entre la précision physique et la convergence prouvable.

Méthodologie
Les auteurs introduisent un nouveau lemme technique (Lemme 1) qui établit la monotonie des gaps spectraux pour des familles spécifiques de Lindbladiens.

• Pour la famille des interactions répétées $\{LG_\kappa\}$, le gap spectral $\lambda_{Gap}$ est non décroissant avec le paramètre $\kappa$ (la largeur de la fonction d'interaction).
• Pour la famille du réservoir macroscopique $\{L^{MB}_\alpha\}$, le gap spectral est non croissant avec la force de couplage $\alpha$.

Cette monotonie permet aux auteurs de « bootstrarapper » (amorcer) les bornes inférieures du gap spectral. Ils peuvent prouver l'existence d'un gap dans un régime où les opérateurs sont quasi-locaux (où les techniques de preuve standard s'appliquent) et extrapoler cette borne vers le régime non local où les modèles physiques sont précis. Ceci est réalisé en exprimant les formes de Dirichlet des Lindbladiens comme des convolutions.

Contributions clés et résultats

1. Préparation efficace via des interactions répétées (Théorème 1) :
   Les auteurs prouvent que l'algorithme d'échantillonnage de Gibbs par interactions répétées prépare des états de Gibbs en temps polynomial pour plusieurs systèmes non commutants, étendant ainsi les cas précédemment connus de systèmes commutants. Le temps de simulation de l'Hamiltonien total passe à $\tilde{O}(N^{10}/\epsilon^4)$ pour :

   • Des modèles de qubits locaux à haute température.
   • Des modèles fermioniques faiblement interagissants à n'importe quelle température constante.
   • Des chaînes de spins 1D à n'importe quelle température constante.
     Ceci démontre que des schémas simples, adaptés à l'ère de la tolérance aux fautes précoce, peuvent préparer efficacement des états de Gibbs complexes.

2. Thermalisation physique efficace (Théorème 2) :
   L'article prouve que les systèmes à corps nombreux ouverts faiblement couplés à des réservoirs à haute température atteignent leur état stationnaire de Gibbs en temps polynomial. Spécifiquement, pour des Hamiltoniens locaux à hautes températures ($\beta \leq \beta^*$), l'état du système devient $\epsilon$-proche de l'état de Gibbs en norme 1 au cours d'un temps $t = \tilde{O}(N^7 \epsilon^{-2})$ avec une force de couplage $\alpha = \tilde{O}(N^{-3}\epsilon)$. Cela fournit une justification de premier principe selon laquelle les Lindbladiens KMS capturent avec précision la relaxation thermique naturelle jusqu'à l'état stationnaire.

3. Thermalisation plus rapide pour les modèles commutants (Théorème 3) :
   Pour les Hamiltoniens commutants (spécifiquement les codes Toriques 2D et les codes CSS apparentés), les auteurs combinent leur technique d'extrapolation avec des résultats récents sur les Inégalités de Log-Sobolev Modifiées (MLSI). Ils dérivent des bornes plus serrées montrant que ces systèmes se thermalisent en un temps $t = \tilde{O}(\gamma_{max}^4 N^4 \tau^2 / \epsilon^2)$. Cela offre une justification dynamique à l'affirmation selon laquelle les codes CSS de basse dimension ne peuvent pas servir de bonnes mémoires quantiques.

Signification
L'article affirme fournir le fondement théorique pour la préparation dissipative de l'état de Gibbs de systèmes complexes dans l'ère de la tolérance aux fautes précoce, une direction qui a connu des progrès expérimentaux récents. En surmontant la barrière de la localité dans les preuves de gap spectral, les auteurs démontrent que :

• Des algorithmes quantiques dissipatifs simples sont capables de préparer efficacement des états de Gibbs complexes.
• La dynamique de Lindblad dérivée de couplages faibles système-réservoir capture avec précision les échelles de temps de la relaxation thermique pour des systèmes à corps nombreux non commutants.
• La monotonie des gaps spectraux est un outil puissant qui peut probablement être appliqué à d'autres Lindbladiens actuellement hors de portée des techniques de preuve standard, suggérant une similitude qualitative des temps de mélange à travers divers échantillonneurs de Gibbs proposés.

Le travail ne propose pas de nouveaux dispositifs expérimentaux mais valide plutôt l'efficacité théorique de modèles algorithmiques et physiques existants, comblant le fossé entre les approximations physiques et les garanties de convergence rigoureuses.