Comment une puce IBM de 104 qubits a simulé une collision de particules

L’équipe formée par Roland C. Farrell, Nikita A. Zemlevskiy, Marc Illa et John Preskill a réalisé une simulation numérique en temps réel de diffusion dans une théorie de champ d’Ising à une dimension, sur le processeur quantique ibm_marrakesh d’IBM, qui compte 104 qubits. Elle a suivi la dynamique avant, pendant et après une collision entre deux particules, avec jusqu’à 45 pas de Trotter et 5 589 portes à deux qubits. Le résultat central est un indice de production inélastique de particules : une part de l’énergie cinétique initiale semble avoir produit une excitation lourde, accompagnée d’une particule légère. 7

Préparer deux paquets d’ondes qui se font face

L’expérience ne consistait pas à envoyer des particules réelles dans un détecteur. Les chercheurs ont encodé dans les qubits un modèle de théorie quantique des champs sur un réseau à une dimension.

Ils ont d’abord préparé deux paquets d’ondes de particules légères, séparés dans l’espace et dotés d’impulsions opposées afin qu’ils se dirigent l’un vers l’autre. Le Hamiltonien de la théorie d’Ising a ensuite été fait évoluer en temps réel au moyen d’une décomposition de Trotter, une méthode qui approxime l’évolution continue par une succession de petites opérations quantiques. La densité d’énergie locale a été mesurée au fil de l’approche, du choc et de l’éloignement des paquets. 7

La préparation de chaque paquet suivait deux étapes :

1. créer un état de type (|W(k_0)\rangle), réglé pour imposer l’enveloppe spatiale, le contenu en impulsion et les nombres quantiques voulus ; cet état initial peut toutefois contenir des composantes à plusieurs particules ;
2. appliquer des circuits variationnels qui préservent les symétries afin de minimiser l’énergie et de rapprocher l’état de l’état propre à une particule recherché. 7

Avant leur exécution sur la machine IBM, les circuits de préparation des paquets pour le grand réseau d’Ising ont été conçus avec des simulations classiques par états produits matriciels, ou MPS. 7

Pourquoi les états W ont été décisifs

Préparer un paquet d’ondes étendu est habituellement coûteux : plus il couvre de sites du réseau, plus le circuit tend à devenir profond. La méthode proposée adapte une préparation d’états W à profondeur constante pour générer l’état initial de type (|W(k_0)\rangle). 27

Elle exploite des mesures effectuées en cours de circuit : certains qubits sont mesurés avant la fin du calcul, puis les résultats classiques déterminent conditionnellement les opérations suivantes — un mécanisme appelé feed-forward. Cette approche permet un protocole de préparation d’état W sans qubit ancillaire, dont la profondeur ne dépend ni de la taille du paquet ni de la dimension spatiale, à l’exception d’une dépendance à la longueur de corrélation physique. 7

Les auteurs décrivent ce gain comme une amélioration superexponentielle par rapport à des méthodes antérieures exigeant soit un prétraitement classique exponentiellement coûteux, soit une profondeur de circuit croissant polynomialement avec la taille du paquet. 7

Comment le caractère inélastique a été repéré

Dans une diffusion élastique, deux particules légères ressortent après le choc avec des impulsions redistribuées, sans apparition d’une nouvelle espèce de particule. Les profils de densité d’énergie en sortie devraient alors être relativement symétriques.

À l’inverse, une particule lourde se déplaçant lentement produit un profil asymétrique après la collision. Les chercheurs ont quantifié cette dissymétrie par le coefficient d’asymétrie, ou skewness, de la densité d’énergie mesurée. Une valeur non nulle, cohérente avec la trajectoire attendue de la particule lourde, a permis de distinguer le canal inélastique du comportement de référence élastique. 7

L’observation ne signifie donc pas seulement que l’énergie s’est déplacée dans le système : elle est compatible avec la conversion de l’énergie cinétique des deux particules légères incidentes en un état final comportant une particule légère et une excitation plus lourde. 7

Ce que cette expérience ne démontre pas encore

Cette démonstration porte sur une théorie de champ d’Ising à une dimension. Ce n’est ni une simulation de la chromodynamique quantique (QCD), ni la reproduction de jets hadroniques en trois dimensions, ni un événement réaliste de collisionneur. 7

La machine utilisée reste un processeur quantique bruité, associé à des techniques d’atténuation des erreurs plutôt qu’à une correction d’erreurs quantiques pleinement tolérante aux fautes. Le résultat constitue donc un indice en faveur du processus modélisé, pas un calcul de diffusion systématiquement précis au niveau requis pour la physique des collisionneurs. 7

La méthode a aussi été construite ou illustrée pour une théorie scalaire à une dimension, le modèle de Schwinger et la théorie de champ d’Ising à deux dimensions. 7 Mais passer à des théories de jauge de dimension supérieure, et plus encore à la QCD, demanderait des machines nettement plus grandes et plus précises, probablement tolérantes aux fautes, ainsi que des méthodes extensibles pour préparer les états, encoder les champs de jauge, faire évoluer le système longtemps et extraire des observables à plusieurs particules. 7