Quantum instanton approach to metastable collective spins

Titel: Quanten-Instanton-Ansatz für metastabile kollektive Spins

Autoren: Krzysztof Ptaszyński, Maciej Chudak und Massimiliano Esposito
Institutionen: Institut für Molekularphysik, Polnische Akademie der Wissenschaften; Universität Luxemburg

---

1. Problemstellung

Kollektive Spinsysteme, bei denen Spin-Ensembles an ein gemeinsames Reservoir gekoppelt sind und effektiv durch einen einzigen Makrospin beschrieben werden, spielen eine zentrale Rolle in der Atom- und Festkörperphysik. Aufgrund ihrer intrinsischen Nichtlinearität können diese Systeme mehrere langlebige metastabile Zustände aufweisen, die schließlich in einen einzigen wahrscheinlichsten Zustand relaxieren.

• Herausforderung: Bei großen Systemgrößen (großer Spin $J$, analog zum Extensivitätsparameter $V$) wird die Analyse der Metastabilität mittels der quantenmechanischen Master-Gleichung (QME) rechnerisch extrem aufwendig.
• Phänomen: In diesem Regime treten oft erste Ordnungs-Phasenübergänge auf, bei denen sich der wahrscheinlichste Zustand abrupt ändert.
• Limitierung bestehender Methoden: Bisherige Ansätze basierten häufig auf semiklassischen Näherungen, wie der Wigner-Verteilung (SW-Ansatz), die auf Fokker-Planck-Gleichungen beruhen. Diese Näherungen vernachlässigen nicht-gaußsche Fluktuationen (Terme höherer Ordnung in der Ableitung), was zu ungenauen Vorhersagen der Aktivierungsbarrieren und damit der Relaxationsraten führt.

2. Methodik

Die Autoren entwickeln einen Echtzeit-Instanton-Ansatz, der auf der exakten Dynamik quantenmechanischer Quasi-Wahrscheinlichkeitsverteilungen basiert, ohne diese zu truncieren.

• Quasi-Wahrscheinlichkeits-Darstellung: Das Dichtematrix-Problem wird in die Husimi- ($p_H$) und P-Darstellung ($p_P$) überführt. Die Dynamik dieser Verteilungen wird durch Differentialoperatoren $L_\alpha$ beschrieben.
• WKB-Ansatz und Hamilton-Jacobi-Gleichung: Im Limes großer Spins ($J \to \infty$) wird der Propagator der Dynamik mittels eines WKB-Ansatzes ($K \sim e^{-J S}$) analysiert. Dies führt zu einer Hamilton-Jacobi-Gleichung für die Wirkung $S$, wobei ein Hilfs-Hamiltonian $H_\alpha$ definiert wird.
• Instanton-Trajektorien: Die Aktivierungsbarrieren $A_{i \to j}$ zwischen metastabilen Zuständen werden als Wirkung entlang spezieller Trajektorien (Instantonen) im Phasenraum bestimmt. Diese Trajektorien verbinden die stabilen Fixpunkte der Mean-Field-Dynamik mit instabilen Sattelpunkten.
• Auswahlkriterien: Da die Hilfs-Hamiltoniane in diesem System nicht-konvex sind (im Gegensatz zu klassischen stochastischen Systemen), ist die Auswahl der physikalisch relevanten Lösung nicht trivial. Die Autoren leiten Randbedingungen her ($\pi=0$ an den Endpunkten) und nutzen die Symmetrie des Modells, um die Trajektorien effizient zu berechnen (Continuation-Methode).

3. Schlüsselbeiträge

1. Verallgemeinerung des Instanton-Ansatzes: Erstmals wird der Instanton-Ansatz explizit auf kollektive Spinsysteme angewendet, die durch Lindblad-Operatoren beschrieben werden, und nicht nur auf bosonische Resonatoren.
2. Überwindung semiklassischer Limitierungen: Der Ansatz verzichtet auf die Truncierung höherer Ableitungsterme (die in der Wigner-Näherung üblich ist). Dadurch werden nicht-gaußsche Fluktuationen korrekt erfasst.
3. Analytische und numerische Validierung: Die Methode wird an einem konkreten Modell (ein Makrospin mit kohärentem Antrieb, linearer Pumpung und nichtlinearer Dissipation) demonstriert. Die Ergebnisse werden direkt mit exakten QME-Simulationen für endliche $J$ verglichen.

4. Ergebnisse

• Aktivierungsbarrieren: Die berechneten Aktivierungsbarrieren $A_{i \to j}$ stimmen hervorragend mit den asymptotischen Skalierungen der Übergangsraten überein, die aus den QME-Simulationen extrahiert wurden.
• Phasenübergang: Der Ansatz sagt den Punkt des ersten Ordnungs-Phasenübergangs korrekt voraus (dort, wo sich die Barrieren $A_{\ell \to u}$ und $A_{u \to \ell}$ schneiden).
  • Vergleich: Der semiklassische Wigner-Ansatz (SW) sagt den Übergang bei einem falschen Parameterwert voraus ($\Gamma \approx 7.9\gamma$ statt $\approx 8.9\gamma$) und unterschätzt die Barrieren signifikant.
• Liouvillian-Lücke: Die Relaxationszeit des Systems wird durch die Liouvillian-Lücke $\lambda$ charakterisiert. Die Skalierung $\lambda \sim \exp(-J A_{\min})$ wird durch den neuen Ansatz korrekt wiedergegeben, während der SW-Ansatz die Relaxationszeit stark unterschätzt.
• Nicht-Monotonie der Wirkung: Im Gegensatz zu klassischen Systemen kann die Wirkung $S$ entlang der Instanton-Trajektorien in diesem quantenmechanischen Kontext nicht-monoton verlaufen und sogar negative Werte annehmen, bevor sie im Limes $t \to \infty$ gegen den korrekten positiven Aktivierungsbarrieren-Wert konvergiert.

5. Bedeutung und Ausblick

• Theoretische Durchbrüche: Die Arbeit liefert den Beweis, dass Instanton-Methoden auf exakten Quanten-Quasi-Wahrscheinlichkeiten basierend, eine präzise Beschreibung von Metastabilität in offenen Quantensystemen jenseits der semiklassischen Näherung ermöglichen.
• Allgemeine Anwendbarkeit: Obwohl auf Spinsysteme fokussiert, ist das theoretische Gerüst allgemein gültig für Systeme mit Quasi-Wahrscheinlichkeits-Dynamik. Es bietet eine konzeptionell einfachere Alternative zu Keldysh-Pfadintegralen für bosonische Systeme und ist auf Spin-Boson-Komplexe sowie Systeme mit lokaler Dissipation übertragbar.
• Experimentelle Relevanz: Die Ergebnisse sind direkt relevant für Experimente in der Cavity-QED, Circuit-QED und bei superradianten Lasern, wo dissipative Phasenübergänge und metastabile Zustände beobachtet werden.

Zusammenfassend etabliert diese Arbeit einen robusten theoretischen Rahmen, um Quanteneffekte in der Dynamik metastabiler makroskopischer Systeme präzise zu quantifizieren, wo herkömmliche semiklassische Methoden versagen.