Solving Dicke superradiance analytically: A compendium of methods

Titel: Analytisch oplossen van Dicke-superradianse: Een compendium van methoden

Auteurs: R. Holzinger, N. S. Bassler, J. Lyne, F. G. Jimenez, J. T. Gohsrich, en C. Genes.

1. Het Probleem

Het artikel behandelt de Dicke-superradianse, een fenomeen waarbij een ensemble van $N$ identieke twee-niveausystemen (TLS) die collectief vervallen, een plotselinge, intense stralingsuitbarsting vertoont.

• Kernmechanisme: In tegenstelling tot een ensemble van onafhankelijke emitters (waar de uitstralingskracht exponentieel afneemt met een snelheid $\propto N\Gamma e^{-\Gamma t}$), vertoont een collectief vervallend systeem een niet-lineaire tijdsafhankelijkheid. Door correlaties tussen de deeltjes neemt de emissiesnelheid toe naarmate het systeem evolueert naar een half-bevolkte toestand, met een piek van $\propto N^2$ in de uitgestraalde vermogen.
• Wiskundige Formulering: Het probleem bestaat uit het bepalen van de tijdsontwikkeling van de dichtheidsoperator $\rho(t)$ voor een initieel volledig geïnverteerd ensemble ($|N\rangle$), beschreven door de Lindblad-mastervergelijking in het interactiekader:
  $$\dot{\rho}(t) = \Gamma \left( S\rho S^\dagger - \frac{1}{2}\{S^\dagger S, \rho\} \right)$$
  Hierbij is $S = \sum \sigma_j$ de collectieve instortingsoperator. Door de permutatiesymmetrie is de dynamica beperkt tot de symmetrische deelruimte, opgespannen door $N+1$ collectieve Dicke-toestanden $|m\rangle$ (waar $m$ het aantal excitaties aangeeft).
• Uitdaging: Hoewel er eerdere analytische oplossingen waren (bijv. door Lee in 1977), werden deze vaak als te complex of computationally demanding beschouwd, wat leidde tot een verschuiving naar benaderingen voor grote $N$ (thermodynamische limiet). Dit artikel beoogt een volledig analytische oplossing te bieden die geldig is voor elk tijdstip en elk aantal spins $N$.

2. Methodologie

De auteurs presenteren vijf verschillende analytische benaderingen om tot dezelfde exacte oplossing te komen. Alle methoden maken gebruik van de observatie dat de dichtheidsoperator in de Dicke-basis diagonaal blijft ($\rho(t) = \sum \rho_m(t) |m\rangle\langle m|$), waardoor het probleem reduceert tot het vinden van de populaties $\rho_m(t)$.

1. Recursieve Vergelijkingen:

   • De oplossing wordt ontwikkeld als een oneindige tijdsreeks $\rho(t) = \sum \frac{t^n}{n!} L^n[\rho_0]$.
   • Er wordt een recursieve relatie afgeleid die de $n+1$-de orde component van een toestand koppelt aan de $n$-de orde componenten van dezelfde en de bovenliggende toestand.
   • Door deze recursie op te lossen en te sommeren, wordt een gesloten vorm gevonden.

2. Combinatorische Benadering:

   • De tijdsreeks wordt opgesplitst in een binomiale expansie van de niet-commuterende superoperatoren $B$ (verlaagt excitatie) en $C$ (behoudt excitatie).
   • Door alle mogelijke paden (volgorde van $B$ en $C$ operatoren) te tellen en te sommeren, wordt een exacte uitdrukking voor $\rho_m(t)$ verkregen.

3. Probabilistische Benadering:

   • De tijd wordt gediskretiseerd. Voor elk klein interval $\Delta t$ wordt de kans berekend op een "jump" (overgang naar een lagere toestand) versus het blijven in dezelfde toestand.
   • De populatie $\rho_m(t)$ wordt gevonden door te sommeren over alle mogelijke paden van de initiële toestand naar de doeltoestand $m$, wat leidt tot een genererende functie die geëvalueerd kan worden via contourintegratie.

4. Niet-Hermitische Hamiltoniaan:

   • De recursieve vergelijkingen worden herschreven als een stelsel gekoppelde differentiaalvergelijkingen in matrixvorm: $\dot{\rho} = \Gamma H \rho$.
   • De matrix $H$ is niet-Hermitisch en bezit uitzonderlijke punten (Exceptional Points - EPs) vanwege de degeneratie van de eigenwaarden ($h_m = h_{\bar{m}}$).
   • De oplossing vereist een Jordan-decompositie van $H$. De aanwezigheid van EP2-punten (Jordan-blokken van grootte 2) introduceert termen van de vorm $t e^{-\lambda t}$ in de tijdsontwikkeling, wat afwijkt van puur exponentieel verval.

5. Quantum Jump Unraveling (De centrale methode):

   • De mastervergelijking wordt "ontward" in kwantumtrajectoires. Elke trajectorie beschrijft een systeem dat $N$ jumps ondergaat van $|N\rangle$ naar $|0\rangle$ op willekeurige tijdstippen.
   • De tijdstippen van de jumps worden bepaald door uniforme willekeurige getallen.
   • Door te middelen over oneindig veel trajectoires (equivalent aan integreren over de ruimte van willekeurige getallen), wordt de exacte oplossing verkregen als een gesloten contourintegraal in het complexe vlak.

3. Belangrijkste Resultaten

De kern van het artikel is de afleiding van een compacte, exacte analytische oplossing voor de populatie $\rho_m(t)$ van een toestand $m$ op tijdstip $t$:

$$\rho_m(t) = \frac{1}{2\pi i} \oint_C dz \, f_m(z, t) = \sum_{p \in \text{poles}} \text{Res}\left[ f_m(z, t) \right]_{z=p}$$

Waarbij de integrand $f_m(z, t)$ is:
$$f_m(z, t) = \frac{(-1)^{N-m} h_N \cdots h_{m+1}}{(z-h_N)\cdots(z-h_m)} e^{-z\Gamma t}$$

• Residusom: De oplossing wordt uitgedrukt als een som van residuen van de functie $f_m(z, t)$ over de polen $p$ in het complexe vlak.
• Polen: De polen corresponderen met de waarden $h_k = k(N+1-k)$.
  • Voor toestanden boven de "evenaar" ($m \ge (N+1)/2$) zijn alle polen enkelvoudig.
  • Voor toestanden onder de evenaar treden gedegenereerde paren op ($h_m = h_{\bar{m}}$), wat leidt tot dubbele polen. De berekening van het residu voor deze dubbele polen vereist een afgeleide, wat de niet-exponentiële tijdsafhankelijkheid (zoals $t e^{-\Gamma t}$) verklaart die voortkomt uit de uitzonderlijke punten.
• Generalisatie: De methode is direct generaliseerbaar voor willekeurige initiële Dicke-toestanden $|m_0\rangle$, niet alleen voor de volledig geïnverteerde toestand.

4. Bijdragen en Significantie

• Unificatie van Methoden: Het artikel biedt een uniek overzicht dat vijf verschillende wiskundige paden (van recursie tot probabilistische trajectoires) unifyt naar één elegante oplossing. Dit bevestigt de robuustheid van het resultaat.
• Toegang tot Exacte Oplossingen: Het demonstreert dat complexe open kwantumsysteemproblemen, die vaak als numeriek moeten worden opgelost, volledig analytisch kunnen worden opgelost voor een brede klasse van systemen met symmetrie.
• Rol van Uitzonderlijke Punten (EPs): Het artikel legt een expliciet verband tussen de fysica van Dicke-superradianse en de theorie van niet-Hermitische systemen. De aanwezigheid van EP2-punten in de effectieve Hamiltoniaan wordt geïdentificeerd als de oorsprong van de specifieke tijdsafhankelijkheid (afwijking van puur exponentieel verval) in de dynamica.
• Praktische Toepassingen: De analytische oplossing is relevant voor moderne experimentele platforms zoals:
  • Cavity QED: Waar atomen in een verliesrijke optische holte worden gevangen (cavity-induced superradiance).
  • Superradiante lasers: Systemen waar de coherentie in het versterkingsmedium wordt opgeslagen in plaats van in het veld.
  • Supergeleidende circuits: Realisatie van superradianse met kunstmatige atomen.

Conclusie

Dit werk levert een fundamentele bijdrage aan de theorie van open kwantumsystemen door een volledig analytische, exacte oplossing te bieden voor het klassieke Dicke-superradianse-probleem. Door de oplossing te formuleren als een contourintegraal en residusom, bieden de auteurs een krachtig wiskundig gereedschap dat niet alleen de dynamica van collectief verval exact beschrijft, maar ook een brug slaat naar de moderne fysica van niet-Hermitische systemen en uitzonderlijke punten. De methode opent de deur voor het analytisch oplossen van complexere, aangedreven dissipatieve systemen.