Quantencomputer simuliert Teilchenkollision: IBM-Prozessor liefert Hinweise auf Teilchenproduktion

Zwei Wellenpakete, eine digitale Kollision

Roland C. Farrell, Nikita A. Zemlevskiy, Marc Illa und John Preskill führten auf IBMs 104-Qubit-Prozessor ibm_marrakesh eine digitale Echtzeit-Streusimulation durch. Modelliert wurde eine eindimensionale Ising-Feldtheorie: Zwei räumlich getrennte, leichte Teilchenwellenpakete bewegten sich aufeinander zu, kollidierten und entwickelten sich anschließend weiter auseinander. Die Berechnung erreichte bis zu 45 Trotter-Schritte und 5.589 Zwei-Qubit-Gatter. 7

Dabei zerlegten die Forschenden die kontinuierliche Zeitentwicklung des Hamiltonoperators in kurze digitale Entwicklungsschritte – die sogenannte Trotterisierung. Während der gesamten Entwicklung maßen sie die lokale Energiedichte im System. So ließ sich verfolgen, wie sich die Energie vor, während und nach dem Zusammenstoß räumlich verteilte. 7

Wie die beiden einlaufenden Teilchen vorbereitet wurden

Jedes einzelne Wellenpaket entstand in zwei Schritten:

1. Zunächst wurde ein maßgeschneiderter, W-Zustands-ähnlicher Zustand ( |W(k_0)\rangle ) erzeugt. Dieser legt die räumliche Form des Pakets, seinen Impulsinhalt und seine Quantenzahlen fest. Zu diesem Zeitpunkt kann der Zustand allerdings noch Mehrteilchenanteile enthalten.
2. Anschließend brachten symmetrieerhaltende variationale Schaltungen den Zustand durch Energieminimierung näher an den gewünschten Einteilchen-Eigenzustand. 7

Die Schaltungen für die Wellenpakete auf dem großen Ising-Gitter wurden zunächst klassisch mithilfe von Matrixproduktzuständen entworfen und danach als Anfangszustand auf der 104-Qubit-Hardware ausgeführt. 7

Woran das Team eine inelastische Kollision erkannte

Bei einer elastischen Streuung würden zwei leichte Teilchen nach dem Zusammenstoß wieder auseinanderlaufen. Ihre Impulse können sich ändern, aber es entsteht keine neue Teilchenart; die Energiedichteprofile wären vergleichsweise symmetrisch. 7

Das Experiment zeigte dagegen eine Asymmetrie im Energiedichteprofil nach der Kollision. Quantifiziert wurde sie über die Schiefe der gemessenen Energiedichte. Ein von null verschiedener Wert, der zur vorhergesagten Bahn eines langsamen schweren Teilchens passte, unterschied den beobachteten Kanal vom elastischen Vergleichsfall. 7

Die Interpretation lautet daher: Ein Teil der Bewegungsenergie der beiden leichten einlaufenden Teilchen wurde in einen Endzustand umgewandelt, der ein leichtes und ein schwereres Teilchen enthält. Das ist ein Hinweis auf inelastische Teilchenproduktion – nicht bloß auf eine Umverteilung von Energie. 7

Warum W-Zustände den Ansatz praktikabler machen

Breitet sich ein Wellenpaket über viele Gitterplätze aus, werden herkömmliche Präparationsmethoden zunehmend tief und damit auf fehleranfälliger Hardware schwer ausführbar. Das Team übertrug deshalb Verfahren zur Erzeugung von W-Zuständen auf die Präparation von Teilchenwellenpaketen. 27

Der zentrale Kniff sind Mittelschaltmessungen: Messungen finden bereits während der Quantenrechnung statt. Ihre klassischen Ergebnisse steuern anschließend bedingt weitere Operationen – ein Verfahren, das als Feed-forward bezeichnet wird. Damit gelang ein ancillafreies W-Zustands-Protokoll mit konstanter Tiefe. Die Schaltungstiefe hängt, abgesehen von der physikalischen Korrelationslänge, nicht von Größe oder räumlicher Dimension des Wellenpakets ab. 7

Die Autoren bezeichnen dies als superexponentielle Verbesserung gegenüber früheren Ansätzen, die entweder exponentiell aufwendige klassische Vorverarbeitung benötigten oder deren Schaltungstiefe polynomial mit der Paketgröße anwuchs. 7

Ein Schritt Richtung komplexerer Quantenfeldtheorien – aber keine QCD-Simulation

Das Präparationsverfahren wurde nicht nur für die eindimensionale Ising-Feldtheorie entwickelt beziehungsweise demonstriert, sondern auch für eine eindimensionale skalare Feldtheorie, das Schwinger-Modell und die zweidimensionale Ising-Feldtheorie. 7

Langfristig geht es um stark gekoppelte Quantenfelddynamik in Echtzeit. Bei hochenergetischen Kollisionen kann Bewegungsenergie neue Teilchen erzeugen; die Vorhersage ihrer Häufigkeiten und ihres kollektiven Verhaltens ist auch für die Hadronenproduktion und letztlich für die Quantenchromodynamik (QCD) relevant. Gerade allgemeine Echtzeitdynamik ist für klassische Verfahren schwierig – hier könnten Quantencomputer künftig Vorteile bieten. 7

Die Grenzen des Ergebnisses sind jedoch klar: Simuliert wurde ein vereinfachtes eindimensionales Ising-Modell, nicht QCD, keine dreidimensionalen Hadronjets und kein realistisches Ereignis eines Teilchenbeschleunigers. Die Hardware war zudem verrauscht und arbeitete mit Fehlerunterdrückung statt vollständiger, fehlertoleranter Quantenfehlerkorrektur. Das Resultat ist somit ein experimenteller und algorithmischer Fortschritt mit Evidenz für den modellierten Prozess, aber noch keine systematisch präzise Streurechnung auf Collider-Niveau. 7

Für höherdimensionale Eichtheorien, insbesondere QCD, wären deutlich größere und genauere, wahrscheinlich fehlertolerante Quantencomputer nötig – ebenso skalierbare Zustandspräparation, Kodierungen von Eichfeldern, lange Echtzeitentwicklungen und präzise Verfahren zur Auswertung von Vielteilchenobservablen. 7