Simple and efficient end-to-end quantum thermal and ground state preparation

Énoncé du problème

La préparation des états thermiques et fondamentaux quantiques est une condition préalable fondamentale pour la simulation de la physique des corps multiples, de la chimie et de la science des matériaux sur les ordinateurs quantiques. Les approches existantes sont confrontées à un compromis persistant entre simplicité, efficacité, généralité et rigueur. Les algorithmes variationnels nécessitent un réglage de paramètres ; les méthodes adiabatiques peuvent être lentes ; et les algorithmes dissipatifs d'ordre élevé (par exemple, basés sur la dynamique de Lindblad) nécessitent souvent des structures de circuits complexes, incluant des évolutions hamiltoniennes contrôlées ou inversées dans le temps ainsi qu'un grand nombre de qubits ancillas. Ces complexités rendent leur mise en œuvre difficile sur les premiers dispositifs quantiques tolérants aux fautes. Inversement, les protocoles simplifiés manquent souvent de garanties théoriques rigoureuses concernant les temps de mélange ou la précision du point fixe. L'objectif est de développer un algorithme suffisamment simple pour le matériel précoce, tout en admettant des garanties de performance rigoureuses de bout en bout pour des hamiltoniens physiquement pertinents.

Méthodologie

Les auteurs proposent un schéma d'interaction système-baignade randomisé qui ne nécessite que l'évolution temporelle vers l'avant sous un hamiltonien système-baignade et un seul qubit ancilla réutilisable. Le protocole fonctionne comme suit :

1. Configuration de l'Hamiltonien : L'hamiltonien total est défini par $H_\alpha(t) = H + H_E + \alpha f(t) (A_S \otimes B_E + A_S^\dagger \otimes B_E^\dagger)$, où $H$ est l'hamiltonien du système, $H_E = -\omega Z/2$ est l'hamiltonien de la baignade dépendant d'une fréquence aléatoire $\omega$, et $A_S$ est un opérateur de couplage choisi aléatoirement dans un ensemble $\mathcal{A}$.
2. Initialisation de la Baignade : Le qubit ancilla (la baignade) est initialisé dans un état $\rho_E \propto \exp(-\beta H_E)$. Pour la préparation de l'état fondamental ($\beta = \infty$), la baignade est initialisée dans $|0\rangle\langle 0|$.
3. Interaction et Évolution : À chaque itération, un opérateur de couplage $A_S$ et une fréquence $\omega$ sont échantillonnés à partir de distributions prescrites. Le système et la baignade évoluent vers l'avant sous l'hamiltonien dépendant du temps $H_\alpha(t)$ avec une enveloppe temporelle lisse (choisie comme une fonction gaussienne).
4. Trace et Réinitialisation : Le qubit ancilla est retracé (trace) et réinitialisé à $\rho_E$. Ce processus définit un canal quantique $\Phi$ agissant sur le système.
5. Itération : Le canal $\Phi$ est itéré $\tau_{mix}$ fois, où $\tau_{mix}$ est le temps de mélange requis pour atteindre le point fixe.

L'implémentation utilise une décomposition de Trotter de second ordre pour simuler l'évolution dépendante du temps, évitant ainsi le besoin d'une évolution hamiltonienne inverse ou d'une logique de contrôle complexe.

Contributions Clés

• Architecture Simple : L'algorithme ne nécessite qu'un seul qubit ancilla et une évolution hamiltonienne vers l'avant, ce qui le rend bien adapté aux dispositifs hybrides analogues-numériques et aux dispositifs tolérants aux fautes précoces.
• Garanties de Point Fixe Rigoureuses : Les auteurs prouvent qu'en concevant soigneusement l'état de la baignade, les opérateurs de couplage et les distributions d'échantillonnage, le point fixe du canal quantique $\Phi$ peut être rendu arbitrairement proche de l'état thermique ou fondamental souhaité.
• Analyse du Temps de Mélange : Le document établit des bornes polynomiales rigoureuses pour le temps de mélange pour plusieurs classes d'hamiltoniens physiquement pertinents, résolvant la tension entre la simplicité algorithmique et la rigueur théorique.
• Efficacité de Bout en Bout : En combinant les bornes d'erreur du point fixe avec les résultats du temps de mélange, les auteurs fournissent une analyse de complexité complète montrant que le temps d'exécution total croît polynomialement avec la taille du système et l'inverse de la précision pour des modèles spécifiques.

Résultats

Le document présente trois théorèmes principaux concernant l'efficacité et la correction de l'algorithme proposé :

1. Approximation du Point Fixe (Théorème 1) : Pour toute précision cible $\epsilon$, il existe des choix de paramètres tels que la distance de trace entre le point fixe du canal et l'état cible est au plus $\epsilon$. Le temps d'évolution total de l'hamiltonien est borné par un polynôme du temps de mélange, de l'inverse de la température (ou du gap spectral) et de $1/\epsilon$.
2. Mélange de l'État Thermique (Théorème 2) : Pour la préparation de l'état thermique de systèmes de spins non interactifs, de systèmes fermioniques quadratiques (à $\beta$ constant) et de systèmes de spins commutants locaux (à haute température), le temps de mélange suit une loi $\text{poly}(N, 1/\epsilon)$.
3. Mélange de l'État Fondamental (Théorème 3) : Pour la préparation de l'état fondamental de systèmes de spins non interactifs et de systèmes fermioniques quadratiques, le temps de mélange suit également une loi $\text{poly}(N, 1/\epsilon)$.

L'analyse repose sur la démonstration que le canal quantique discret $\Phi$ approxime une dynamique de Lindblad en temps continu. Une contribution technique critique consiste à démontrer que le temps de mélange reste borné (indépendant du paramètre de largeur du filtre $\sigma$) même lorsque $\sigma$ est augmenté pour réduire l'erreur du point fixe. Ceci est réalisé en concevant l'interaction pour permettre de grandes transitions d'énergie tout en maintenant la perturbation de l'état de point fixe approximatif très faible.

Signification et Revendications

Les auteurs affirment que ce travail fournit la première garantie d'efficacité de bout en bout dans le cadre de l'interaction système-baignade pour la préparation des états thermiques et fondamentaux. Plus précisément :

• C'est la première étude établissant des garanties de temps de mélange rigoureuses pour les états thermiques et fondamentaux dans ce cadre, posant un fondement théorique pour les résultats empiriques observés dans des études connexes.
• L'algorithme évite la complexité des simulations de Lindblad d'ordre élevé (qui nécessitent souvent une évolution inversée dans le temps et de nombreux ancillas) tout en maintenant des garanties de performance rigoureuses.
• Le protocole sert de référence pour l'ingénierie dissipative et le co-design avec la correction d'erreurs, car il teste la capacité du matériel à supporter une évolution de longue durée vers un point fixe stable en utilisant un circuit minimal.
• Ce travail suggère que le cadre d'interaction faible système-baignade peut être considéré comme une approche approximative d'ordre inférieur à la simulation de Lindblad, offrant une voie pratique pour les dispositifs précoces tolérants aux fautes.

Le document reconnaît que bien que l'implémentation actuelle utilise une Trotterisation de second ordre (ce qui peut entraîner des ralentissements asymptotiques par rapport aux méthodes d'ordre supérieur), la simplicité de l'architecture est un avantage primaire pour les dispositifs de l'ère proche. Il note également que des améliorations futures pourraient impliquer une Trotterisation d'ordre supérieur si l'évolution inverse devient efficace sur des matériels spécifiques.