Orbital-optimized spin-adapted multistate contracted VQE for excited states and properties on quantum hardware

Diese Arbeit stellt die orbitale optimierte multistatige kontrahierte VQE-Methode (oo-MC-VQE) vor, welche Spin-adaptierte Operatoren nutzt, um Grund- und angeregte Zustände sowie deren Eigenschaften auf Quantenhardware effizient zu berechnen, während sie durch eine lineare Parameterskalierung mit der Anzahl der Zustände ein Gleichgewicht zwischen Genauigkeit und Schaltkreiskomplexität herstellt.

Das Wesentliche

Stellen Sie sich vor, Sie versuchen, ein riesiges, komplexes Orchester so zu stimmen, dass es eine perfekte Sinfonie spielt. In der Welt der Chemie ist das „Orchester“ ein Molekül, und die „Musik“ ist die Art und Weise, wie sich seine Elektronen bewegen und interagieren. Um zu verstehen, wie ein Molekül Licht absorbiert (was Farben erzeugt und Prozesse wie die Photosynthese antreibt), müssen Wissenschaftler die exakten Noten berechnen, die diese Elektronen spielen.

Lange Zeit war es, die Musik für Moleküle mit vielen Elektronen zu berechnen, so, als würde man versuchen, ein Puzzle zu lösen, das mit jedem zusätzlichen Teil exponentiell schwieriger wird. Klassische Computer (die, die wir heute verwenden) stoßen schließlich an eine Wand und können diese Rätsel für komplexe Moleküle nicht mehr lösen.

Dieses Paper stellt einen neuen Weg vor, um diese Rätsel mithilfe von Quantencomputern zu lösen – speziellen Maschinen, die darauf ausgelegt sind, mit dieser Art von Komplexität auf natürliche Weise umzugehen. Hier ist eine einfache Aufschlüsselung dessen, was die Autoren getan und herausgefunden haben:

1. Das Problem: Viele Noten gleichzeitig stimmen

Normalerweise versuchen Wissenschaftler, das Orchester so zu stimmen, dass es nur eine einzige Note (den Grundzustand) perfekt spielt. Aber um zu verstehen, wie ein Molekül auf Licht reagiert, müssen sie über viele verschiedene Noten (angeregte Zustände) gleichzeitig Bescheid wissen.

• Die Herausforderung: Wenn man versucht, 10 verschiedene Stücke gleichzeitig zu stimmen, werden die Anweisungen (der Computerkreislauf) unglaublich lang und kompliziert. Wenn die Anweisungen zu lang sind, wird der Quantencomputer durch „Rauschen“ (Störungen oder Fehler) verwirrt, und die Musik bricht zusammen.
• Der Kompromiss: Man benötigt einen komplexen Kreislauf, um eine genaue Antwort zu erhalten, aber ein komplexer Kreislauf neigt auf aktuellen, verrauschten Maschinen eher dazu, fehlerhaft zu sein.

2. Die Lösung: Ein intelligenter, symmetrischer Dirigent

Die Autoren entwickelten eine neue Methode namens oo-MC-VQE. Betrachten Sie dies als einen „intelligenten Dirigenten“ für das Quantenorchester.

• Spin-angepasst: In der Quantenchemie besitzen Elektronen eine Eigenschaft namens „Spin“ (wie kreiselnde Spielzeuge). Die Autoren bauten ihre Methode so, dass der Dirigent immer sicherstellt, dass die Kreisel auf die korrekte, symmetrische Weise rotieren. Dies verhindert, dass die Musik aufgrund von Symmetriefehlern „verstimmt“ wird.
• Orbital-optimiert: Sie ließen den Dirigenten zudem den Sitzplan der Musiker (die Orbitale) neu arrangieren, um die Musik schon vor dem Beginn der komplexen Abstimmung besser klingen zu lassen.
• Multistate-kontrahiert: Anstatt zu versuchen, 10 Stücke mit 10 separaten, massiven Handbüchern abzustimmen, fanden sie einen Weg, ein gemeinsames, effizientes Set an Anweisungen zu nutzen, das für alle Stücke gleichzeitig funktioniert.

3. Die Entdeckung: Lineares Wachstum

Eine der größten Fragen war: Wenn ich 10 Zustände statt 1 berechnen möchte, brauche ich dann 10-mal mehr Computerleistung?

• Das Ergebnis: Die Autoren fanden heraus, dass die Antwort überraschend einfach ist. Der Menge an „Aufwand“ (Kreislaufparameter), die für den Computer benötigt wird, wächst linear. Wenn Sie die Anzahl der zu berechnenden Zustände verdoppeln, verdoppelt sich im Wesentlichen auch die Länge des Handbuchs. Das ist eine gute Nachricht, denn es bedeutet, dass die Methode skalierbar ist.

4. Der Realitätstest: Spielen auf einer verrauschten Bühne

Die Autoren haben ihre Methode nicht nur auf einem perfekten Computer simuliert; sie haben sie tatsächlich auf echter Quantenhardware (IBM-Quantencomputer) ausgeführt.

• Der Aufbau: Sie testeten zwei kleine Moleküle: Formaldehyd (ein häufig vorkommendes chemisches Molekül) und ein Trihydrogen-Kation ($H_3^+$).
• Das Rausch-Problem: Echte Quantencomputer sind wie eine Bühne mit einer lauten Menge und flackernden Lichtern. Ohne Hilfe waren die Ergebnisse ungeordnet.
• Die Lösung: Sie nutzten „Error Mitigation“-Techniken (Fehlerminimierung). Stellen Sie sich das wie einen Tontechniker vor, der eine Software verwendet, um das Publikumsrauschen und das Flackern der Lichter nach der Aufführung herauszufiltert.
• Das Resultat:
  • Bei Formaldehyd funktionierte die Methode recht gut. Selbst mit dem Rauschen konnten sie die „Peaks“ (Spitzen) im Absorptionsspektrum (den Farben, die das Molekül absorbiert) deutlich erkennen.
  • Bei $H_3^+$ war das Rauschen ein größeres Problem und verschob die Ergebnisse erheblich. Die Autoren merkten an, dass die Mathematik dieses spezifischen Moleküls empfindlicher auf Rauschen reagiert (wie ein empfindliches Instrument, das leicht verstimmt).
  • Wichtigste Erkenntnis: Obwohl die Zahlen auf den echten Maschinen nicht perfekt waren, war die Form der Ergebnisse korrekt. Sie konnten immer noch die Hauptmerkmale des Verhaltens des Moleküls erkennen.

Zusammenfassung

Das Paper zeigt, dass Wissenschaftler durch einen intelligenten, symmetrischen Ansatz in der Lage sind, das Verhalten angeregter Elektronen in Molekülen mithilfe aktueller, unvollkommener Quantencomputer zu berechnen. Sie haben bewiesen, dass die Berechnung mehrerer Zustände keine unmögliche Menge an Ressourcen erfordert, und dass man mit einigen „Rauschunterdrückungs-Tricks“ bereits heute nützliche chemische Erkenntnisse aus echten Quantengeräten gewinnen kann.

Was sie NICHT behauptet haben:
Das Paper behauptet nicht, dass diese Methode sofort neue Solarzellen entwerfen, Krankheiten heilen oder neue Materialien erschaffen kann. Es konzentriert sich strikt darauf, zu beweisen, dass die Methode funktioniert, um Spektren auf Quantenhardware zu berechnen. Jegliche zukünftigen Anwendungen sind durch das allgemeine Fachgebiet impliziert, stellen aber keine spezifischen Behauptungen dieser Studie dar.

Technische Zusammenfassung

Technischer Überblick: Orbital-optimiertes, spin-adaptiertes Multistate-kontrahiertes VQE

Problemstellung
Die robuste und ausgewogene Berechnung elektronisch angeregter Zustände und ihrer Übergangseigenschaften bleibt eine erhebliche Herausforderung, insbesondere in Regimen mit starker elektronischer Reorganisation, Nahezu-Entartungen oder multikonfigurativem Charakter. Lineare-Antwort-Methoden (LR) sind in vielen Kontexten erfolgreich, behandeln angeregte Zustände jedoch als Störungen eines Grundzustands-Referenzzustands, was zu konzeptionellen Einschränkungen in der Nähe von konischen Durchschneidungen und in stark korrelierten Systemen führt. Darüber hinaus stehen bestehende Variational Quantum Eigensolver (VQE)-Ansätze für angeregte Zustände, wie etwa Subspace-Search VQE (SS-VQE) und Multistate-kontrahiertes VQE (MC-VQE), oft vor einem Kompromiss zwischen der Genauigkeit bei der Beschreibung einer steigenden Anzahl elektronischer Zustände und der Expressivität (Schaltkreiskomplexität) des zugrunde liegenden Ansatzes. Zudem ist die aktuelle Quantenhardware durch Rauschen und Dekohärenz begrenzt, was die Schaltungstiefe einschränkt und die Entwicklung von Methoden erforderlich macht, die sowohl effizient als auch mit Fehlerminderungsverfahren kompatibel sind.

Methodik
Die Autoren führen die Methode des orbital-optimierten Multistate-kontrahierten variativen Quanten-Eigenwert-Lösers (oo-MC-VQE) ein, die spin-adaptierte Operatoren zur Berechnung von Grund- und angeregten Zuständen sowie zustands- und übergangsspezifischen Eigenschaften nutzt.

• Spin-Adaption: Um die Spin-Symmetrie der Referenzzustände während der VQE-Optimierung zu gewährleisten, verwendet die Methode einen spin-adaptierten Operator-Pool. Konkret nutzen die Autoren spin-adaptierte Ein-Anregungs-Operatoren für Singletts und Paar-Zwei-Anregungs-Operatoren. Dies vermeidet die Notwendigkeit von Spin-Straftermen und stellt die korrekte Spin-Symmetrie aller inkludierten Zustände sicher, ohne redundante Parameter einzuführen.
• Multistate-kontrahierter Rahmen: Die Methode optimiert mehrere elektronische Zustände gleichzeitig unter Verwendung eines zustandsgemittelten Energie-Funktionals. Jeder Zustand wird unter Verwendung eines festen, parameterfreien Zustandspräparations-Schaltkreises (Initialisierung orthogonaler Basiszustände) gefolgt von einem parametrisierten Korrelations-Schaltkreis (unter Verwendung eines Unitary Product State, oder UPS, Ansatzes) konstruiert, der Elektronenkorrelation einführt. Die Optimierung zielt auf die zustandsgemittelte Energie ab, gefolgt von einer klassischen Diagonalisierung der Multistate-Hamilton-Matrix, um die einzelnen Zustandsenergien und Wellenfunktionen aufzulösen.
• Orbital-Optimierung: Der Rahmen wird durch die Einbeziehung der Orbital-Optimierung (oo-MC-VQE) erweitert. Orbital-Rotationen werden als klassische Kostenfunktion behandelt, was die gleichzeitige Optimierung der Schaltkreisparameter ($\theta$) und der Orbital-Rotationsparameter ($\kappa$) ermöglicht. Die Autoren zeigen, dass für den spezifischen verwendeten Ansatz benachbarte Ein-Anregungs-Schaltkreisparameter redundant sind und entfernt werden können, was die Komplexität des Schaltkreises reduziert.
• Bewertung von Matrixelementen: Um Matrixelemente zwischen verschiedenen Zuständen zu bewerten (die für die Multistate-Hamilton-Matrix und Übergangseigenschaften erforderlich sind), verwenden die Autoren eine Superpositions-Dekompositionstechnik. Dieser Ansatz vermeidet tiefe Zustandspräparations-Schaltkreise für Multi-Determinanten-Zustände, indem er einfachere Superpositionen von zwei Determinanten auswertet, wodurch die Anzahl der CNOT-Gatter reduziert wird, allerdings auf Kosten erhöhter Erwartungswert-Messungen.
• Fehlerminderung: Die Implementierung integriert Mansatz0, eine ansatzbasierte Readout- und Gate-Fehlerminderungstechnik. Für Multistate-Berechnungen unter Verwendung von Superpositionen von Referenzzuständen schlagen die Autoren eine Erweiterung (M0+) vor, die Fehler für alle eindeutigen Multi-Determinanten-Referenzzustände charakterisiert, weisen jedoch auf die exponentielle Skalierung dieses Ansatzes hin.

Wesentliche Beiträge

1. Spin-adaptiertes oo-MC-VQE: Die Entwicklung einer spin-adaptierten Formulierung für den orbital-optimierten zustandsgemittelten VQE, die eine korrekte Erhaltung der Spin-Symmetrie über mehrere Zustände hinweg gewährleistet.
2. Skalierungsanalyse: Eine explizite Untersuchung des Kompromisses zwischen Genauigkeit und Schaltkreiskomplexität. Die Autoren stellen fest, dass die Anzahl der benötigten Schaltkreisparameter, um genaue Ergebnisse zu erzielen, annähernd linear mit der Anzahl der in die Zustandsmittelung einbezogenen Zustände skaliert.
3. Hardware-Implementierung: Eine explizite Quantenschaltkreis-Implementierung der Methode, einschließlich der Reduktion von CNOT-Gattern durch Determinanten-Superposition und der Integration von Fehlerminderungstechniken.
4. Demonstration auf realer Hardware: Durchführung der Methode auf realer Quantenhardware (IBM ibm marrakesh und ibm aachen), um Absorptionsspektren für Benchmark-Molekülsysteme (Formaldehyd und $H_3^+$) zu berechnen.

Ergebnisse

• Schaltkreiskomplexität: Simulationen an LiH, $H_2O$ und $NH_3$ unter Verwendung der ADAPT- und oo-ADAPT-Ansätze zeigen, dass die Anzahl der benötigten Parameter in etwa linear mit der Anzahl der Zustände ansteigt. Beispielsweise erforderte das Erreichen eines Fehlers unter $10^{-4}$ Hartree für $H_2O$ jeweils 27, 32, 24 und 30 Parameter pro Zustand für 1, 2, 3 bzw. 4 Zustände.
• Einfluss der Orbital-Optimierung: Während die Orbital-Optimierung die Konvergenz in der frühen Phase des ADAPT-Verfahrens verbessert, hat sie einen vernachlässigbaren Effekt auf die endgültige Anzahl der benötigten Schaltkreisparameter für eine engere Konvergenz. In einigen Fällen (z. B. $NH_3$ mit 3 oder 4 Zuständen) hatte das Einbetten von Ein-Anregungen in die Orbital-Optimierung einen negativen Effekt auf die Parameteranzahl.
• Leistung der Quantenhardware:
  • Simuliertes Rauschen: Auf emulierter Hardware unter Verwendung von IBM-Rauschmodellen erzeugten Fehlerminderungstechniken (M0 und M0+) Ergebnisse, die nahezu identisch mit rauschfreien Simulationen waren. Das $H_3^+$-System wurde durch Rauschen stärker beeinträchtigt als Formaldehyd, was auf eine höhere Konditionszahl der Multistate-Hamilton-Matrix zurückzuführen ist.
  • Reale Hardware: Berechnungen auf realen Geräten zeigten, dass die fehlermindernden Ergebnisse zwar die qualitativen Trends und primären Merkmale der Absorptionsspektren bewahrten, die quantitative Genauigkeit jedoch signifikant niedriger war als auf simulierten Geräten. Diese Diskrepanz verdeutlicht die Einschränkungen aktueller Rauschmodelle hinsichtlich zeitabhängigen Rauschens und Übersprechen (Crosstalk). Dennoch konnte die Methode trotz Hardware-Rauschens aussagekräftige chemische Erkenntnisse gewinnen.

Bedeutung und Ansprüche
Das Paper behauptet, dass die oo-MC-VQE-Methode eine ausgewogene Beschreibung der Elektronenstruktur bietet, die für anspruchsvolle Szenarien wie konische Durchschneidungen geeignet ist, und damit die Einschränkungen herkömmlicher Linear-Response-Ansätze adressiert. Die Autoren betonen, dass die lineare Skalierung der Schaltkreisparameter mit der Anzahl der Zustände darauf hindeutet, dass die Methode für das Anvisieren mehrerer angeregter Zustände auf Near-Term-Quantenhardware praktikabel ist. Während sie einräumen, dass die quantitative Genauigkeit auf heutigen Geräten eine Herausforderung bleibt, zeigt die Arbeit, dass durch Kombination der Methode mit geeigneten Fehlerminderungstechniken aussagekräftige chemische Erkenntnisse, wie etwa Absorptionsspektren, gewonnen werden können. Die Studie dient als Proof-of-Concept für die Durchführung von Multistate-, spin-adaptierten und orbital-optimierten Berechnungen auf realen Quantengeräten.