Determining the ensemble N-representability of Reduced Density Matrices

Technische Zusammenfassung: Bestimmung der Ensemble-N-Repräsentierbarkeit reduzierter Dichtematrizen

Problemstellung
Das N-Repräsentierbarkeitsproblem ist eine fundamentale Herausforderung in der Elektronentheorie, die sich mit den notwendigen und hinreichenden Bedingungen befasst, die eine $p$-Körper-reduzierte Dichtematrix ($p$-RDM) erfüllen muss, um aus einer $N$-Elektronen-Wellenfunktion oder einem gemischten (Ensemble-)Zustand ableitbar zu sein. Während erhebliche Fortschritte hinsichtlich der reinen Zustand-N-Repräsentierbarkeit erzielt wurden, bleibt die Bestimmung der Ensemble-N-Repräsentierbarkeit einer $p$-RDM eine kritische Hürde, insbesondere für Anwendungen mit thermischen Mischungen, offenen Quantensystemen und statistischen Ensembles. Die Verifizierung, ob eine gegebene $p$-RDM einem physikalischen $N$-Elektronen-Zustand entspricht, ist im allgemeinen Fall rechnerisch unpraktikabel und gehört zur Komplexitätsklasse der Quantum Merlin-Arthur-vollständigen Probleme. Bestehende Näherungsbedingungen sind oft unzureichend für stark korrelierte Systeme, und der vollständige Satz notwendiger und hinreichender Bedingungen für $p$-RDMs (wobei $p > 1$) ist unbekannt.

Methodik
Diese Arbeit schlägt ein praktisches Framework zur Bestimmung der Ensemble-N-Repräsentierbarkeit einer $p$-Körper-Matrix vor, indem ein zuvor entwickeltes unitäres Evolutionsverfahren (reine ADAPT-VQA) auf den Ensemble-Fall erweitert wird. Die Kernmethodik beruht auf einer Purifizierungsstrategie:

1. Purifizierungs-Einbettung: Folgend dem Schrödinger-Hughston-Jozsa-Wootters-Theorem wird ein Ensemble-Zustand $\rho$, definiert auf einem Hilbert-Raum $\mathcal{H}_s$, in einen reinen Zustand $|\Psi_{sb}\rangle$ eingebettet, der auf einem erweiterten Hilbert-Raum $\mathcal{H}_{sb} = \mathcal{H}_s \otimes \mathcal{H}_b$ definiert ist, wobei $\mathcal{H}_b$ ein Hilfsraum ist. Die Purifizierung wird so konstruiert, dass $\rho = \text{Tr}_b(|\Psi_{sb}\rangle\langle\Psi_{sb}|)$ gilt, wodurch sichergestellt wird, dass die reduzierten Dichtematrizen des purifizierten Zustands jene des ursprünglichen Ensembles reproduzieren.
2. Variationale unitäre Evolution: Der Algorithmus verwendet einen adaptiven, derivativ-zusammengesetzten Pseudo-Trotter (ADAPT) variativen Quantenalgorithmus (VQA). Er wendet iterativ variationale Unitaritäten auf einen initialen purifizierten Zustand an (konstruiert aus einem initialen $N$-Teilchen-Zustand $|\phi_0\rangle$ und dessen Purifizierung).
3. Operator-Pool: Die Evolution wird durch einen antihermiteschen Anregungs-Deexzitierungs-Operator-Pool getrieben, der für Fermionen geeignet ist und auf kanonische Ensemble-Beschränkungen beschränkt ist (keine Änderung der Teilchenzahl im System oder im Hilfsraum). Diese Operatoren werden via Jordan-Wigner-Transformation auf Pauli-Operatoren abgebildet, was den Algorithmus für die Simulation auf Quantencomputern geeignet macht.
4. Optimierungsziel: Der Algorithmus minimiert die Hilbert-Schmidt-Distanz ($D$) zwischen der $p$-RDM des evolvierten purifizierten Zustands und einer spezifizierten Ziel-$p$-RDM ($p\rho_{\text{target}}$).
   • Wenn der Algorithmus zu einer Distanz $D_{\min} \approx 0$ konvergiert, wird das Ziel als Ensemble-N-repräsentierbar betrachtet.
   • Wenn $D_{\min} > 0$, dient der Betrag der Distanz als quantitatives Maß für den Grad der N-Repräsentierbarkeitsverletzung.
   • Der evolvierte Zustand liefert den „nächsten“ physikalischen, Ensemble-N-repräsentierbaren reduzierten Zustand zum defekten Ziel.

Wesentliche Beiträge

• Ensemble-ADAPT-VQA-Algorithmus: Die Einführung eines hybriden Frameworks, das reine Zustand-Repräsentierbarkeitsalgorithmen mittels Purifizierung auf Ensemble-Zustände erweitert. Dies ermöglicht es dem Algorithmus, zu unterscheiden, ob eine Zielmatrix aus einem reinen Quantenzustand oder einer Ensemble-Mischung stammt.
• Fehlerkorrektur und Verfeinerung: Die Methodik bietet einen Weg zur Verfeinerung von Dichtematrizen. Sie kann defekte Ensemble-reduzierte Dichtematrizen (die die N-Repräsentierbarkeitsbedingungen verletzen) korrigieren, indem sie sie in Richtung des nächstgelegenen physikalischen Zustands evolviert.
• Quantenzustand-Rekonstruktion: Der Algorithmus erleichtert die Quantenzustand-Rekonstruktion auf Quantencomputern und ermöglicht die Rückgewinnung einer kompatiblen Wellenfunktion oder eines purifizierten Zustands aus einer Ziel-RDM.
• Numerische Validierung: Das Framework wird durch numerische Simulationen an sowohl Modellsystemen als auch molekularen Systemen bei endlicher Temperatur validiert, wobei die Fähigkeit zur Handhabung von 1-RDMs und 2-RDMs demonstriert wird.

Ergebnisse
Die Autoren validierten den Algorithmus mittels numerischer Simulationen auf einem rauschfreien Quantensimulator:

• Modellsysteme (Nulltemperatur): Tests wurden an 4-Elektronen-Systemen mit 3 und 4 Orbitalen durchgeführt.
  • Für 1-RDMs konnte der Algorithmus erfolgreich zwischen reinen und Ensemble-N-repräsentierbaren Zuständen unterscheiden. In Fällen, in denen die Klyachko-Bedingungen für reine N-Repräsentierbarkeit (verallgemeinerte Pauli-Bedingungen) verletzt wurden (gemischte Zustände), versagte die reine ADAPT-VQA bei der Konvergenz gegen eine Distanz von Null, während die Ensemble-ADAPT-VQA zu $D_{\min} \approx 0$ konvergierte.
  • Für 2-RDMs, für die vollständige reine N-Repräsentierbarkeitsbedingungen unbekannt sind, identifizierte der Algorithmus korrekt die Ensemble-N-Repräsentierbarkeit für alle konstruierten Ziele. Er demonstrierte ferner die Fähigkeit zu detektieren, ob eine Ensemble-2-RDM zu einer reinen 1-RDM kontrahieren kann.
  • Rauschrobustheit: Durch das Einführen von zufälligem Rauschen in die Zielmatrizen zeigten die Studien, dass $D_{\min}$ mit der Rauschstärke ($\epsilon$) skaliert, was die Nützlichkeit des Algorithmus zur Quantifizierung der Schwere von N-Repräsentierbarkeitsverletzungen bestätigt.
• Molekulare Systeme (Endliche Temperatur): Der Algorithmus wurde auf $H_2$- und lineare $H_3$-Moleküle bei endlichen Temperaturen angewendet. Zielmatrizen wurden aus kanonischen Ensemble-Thermizuständen konstruiert und mit Rauschen perturbiert.
  • Die Ergebnisse deuteten darauf hin, dass die konvergierte Distanz $D_{\min}$ empfindlicher auf das hinzugefügte Rauschen ($\epsilon$) reagiert als auf die intrinsischen Korrelationseigenschaften des zugrunde liegenden rauschfreien Zustands.
  • Mit zunehmendem Rauschen stieg $D_{\min}$ signifikant an (z. B. auf $\sim 10^1$ für 2-RDMs bei $\epsilon=0.1$), was die Kapazität des Algorithmus unterstreicht, Repräsentierbarkeitsdefekte in realistischen molekularen Szenarien zu detektieren und zu quantifizieren.

Bedeutung
Das Paper behauptet, dass diese Arbeit ein praktisches und robustes Framework zur Adressierung des Ensemble-N-Repräsentierbarkeitsproblems bietet, einer Herausforderung, die historisch gesehen rechnerisch unpraktikabel war. Durch die Nutzung von Purifizierung und variativen Quantenalgorithmen bietet die vorgeschlagene Methode eine duale Kapazität: Sie kann die Natur einer gegebenen RDM klassifizieren (rein vs. Ensemble) und defekte Matrizen korrigieren, um physikalische Beschränkungen durchzusetzen. Dies ist signifikant für die Verbesserung der Genauigkeit quantenchemischer Berechnungen, insbesondere für stark korrelierte Elektronensysteme und Anwendungen bei endlicher Temperatur. Die Autoren betonen, dass der Algorithmus mit Quantenhardware kompatibel ist und somit einen Weg für die Verfeinerung von Dichtematrizen und die Quantenzustand-Rekonstruktion in zukünftigen Quantensimulationen eröffnet.