Semiclassical Simulation of Homogeneous Emitter Ensembles with Local Dissipation

Probleemstelling

Homogene ensembles van emitters (zoals atomen of vaste-stof defecten) zijn fundamenteel voor kwantumoptische technologieën. In moderne experimenten coëxisteren collectieve interacties vaak met lokale dissipatie (bijvoorbeeld lokale dephasing of decay). Hoewel exacte simulaties mogelijk zijn voor kleine systemen, worden deze beperkt tot bescheiden systeemgroottes door de exponentiële schaal van de Hilbertruimte. Bestaande benaderingen zoals cumulant-ontwikkelingen, tensornetwerken en fase-ruimtemethoden hebben moeite om kwantitatief gecontroleerde voorspellingen te doen voor grotere systemen, vooral wanneer lokale dissipatie een rol speelt in permutatie-invariante modellen. Er is een dringende behoefte aan schaalbare methoden die de overgang van microscopische modellen naar emergent gedrag kunnen overbruggen.

Methodologie: Truncated Wigner Approximation (TWA)

De auteurs ontwikkelen een Truncated Wigner Benadering (TWA) specifiek voor permutatie-invariante emitter-ensembles die onderhevig zijn aan lokale dissipatie. De kern van de methode is een kwantum-naar-klassiek mapping in een uitgebreide fase-ruimte.

1. Uitgebreide Fase-ruimte:

   • Traditionele TWA voor collectieve spins gebruikt de Bloch-sfeer (variabelen $\phi, \theta$).
   • Om lokale dissipatie te modelleren, introduceren de auteurs twee extra variabelen: een extra hoek $\psi$ en een totale spin-grootte $J$.
   • De volledige fase-ruimte wordt beschreven door $Z = (\phi, \theta, \psi, J)$.
   • De variabelen $\phi$ en $\theta$ beschrijven de collectieve dynamiek (op de oppervlakte van de Bloch-sfeer), terwijl $\psi$ en $J$ nodig zijn om de veranderingen in de totale spin te vangen die worden veroorzaakt door lokale dissipatieprocessen (die de spin naar het binnenste van de sfeer drijven).

2. Wiskundige Formulering:

   • De dichtheidsmatrix wordt uitgedrukt in de hoekmomentum-representatie ($|JM\rangle$).
   • Er wordt een Stratonovich-Weyl (SW) correspondentie toegepast voor systemen met variabele spinlengte.
   • De lokale dissipatie (beschreven door jump-operatoren $\hat{l}_{jq}$) koppelt de spin-inversie $M$ en het totale spin-getal $J$.
   • Door een semi-klassieke expansie toe te passen (waarbij $J \gg 1$ als continu wordt behandeld), wordt de Master-vergelijking omgezet in een Fokker-Planck vergelijking met niet-negatieve diffusie.
   • Dit wordt ontvouwd naar een systeem van stochastische differentiaalvergelijkingen (SDE's):
     $$\frac{dZ}{dt} = \{Z, H\} + \text{Re}\left[\sum_k \{Z, \sqrt{r_k}L_k\} \circ (\eta_k + i\sqrt{r_k}L_k^*)\right]$$
     Waar $\{,\}$ de Poisson-haakjes zijn, $H$ de Hamiltoniaan, $L_k$ de Weyl-symbolen van de jump-operatoren, $r_k$ de rates, en $\eta_k$ complexe witte ruis.

3. Schaalbaarheid:

   • De dimensie van de fase-ruimte schaalt lineair met het aantal verschillende ensembles in een samengesteld systeem.
   • Voor een systeem met $M$ ensembles worden $4M$ variabelen gesimuleerd, wat mogelijk maakt om systemen met tienduizenden emitters te modelleren.

Belangrijkste Bijdragen

• Generalisatie van TWA: Voor het eerst wordt een TWA geformuleerd die zowel collectieve als lokale dissipatie kan behandelen in permutatie-invariante systemen door de fase-ruimte uit te breiden met de variabele $J$ en de hoek $\psi$.
• Schaalbare Simulatie: De methode maakt simulaties mogelijk voor systemen met tot $10^6$ emitters, een schaal die ontoegankelijk is voor exacte diagonalisatie of Monte Carlo-simulaties van quantum-trajecten.
• Kwantum-Klassieke Overgang: Het biedt een brug tussen microscopische Master-vergelijkingen en emergent collectief gedrag in uitgebreide licht-materie systemen.

Resultaten

De auteurs valideren de methode via benchmarks en toepassing op complexe scenario's:

1. Nauwkeurigheid en Schaal:

   • Vergelijkingen met exacte resultaten voor kleine systemen ($N=10$) tonen een afwijking van minder dan 10%.
   • De nauwkeurigheid verbetert systematisch naarmate het aantal emitters $N$ toeneemt (fouten schalen als $1/N$).
   • Voor $N \approx 100$ komt de TWA zeer dicht bij de exacte dynamiek.

2. Niet-klassieke Signatures:

   • De methode slaagt erin kwantum-squeezing en subradiantie (onderdemping) correct te vangen.
   • Subradiantie, waarbij destructieve interferentie de emissie onderdrukt, is vaak afwezig in gemiddelde-veldbenaderingen, maar wordt hier correct gereproduceerd.
   • De methode faant alleen in het extreem subradiante regime waar de totale spin $J$ zeer klein wordt ($O(1)$), wat de grens van de semi-klassieke expansie markeert.

3. Samengestelde Systemen (1D Keten):

   • De auteurs simuleren een oneindige 1D-rij van gekoppelde resonatoren met $M$ ensembles.
   • Ze observeren ruimtelijke coherentie en een overgang naar selectieve directionaliteit van de emissie.
   • Bij het toepassen van een zwakke resonante probe op een gepompt systeem, vertoont het systeem een sterke, systeemgrootte-afhankelijke respons waarbij de emissie sterk in één richting wordt gericht (unidirectionele emissie) voor grote $M$ (tot 500 ensembles).

4. Boundary Time Crystals:

   • De methode wordt toegepast op boundary time crystals onder lokale dephasing en decay, waarbij het de demping van oscillaties correct voorspelt voor systemen met $N=10^4$.

Significantie

Deze studie opent de deur tot het bestuderen van uitgebreide licht-materie systemen in het mesoscopische regime, waar exacte methoden falen.

• Technologische Impact: Het biedt een krachtig numeriek gereedschap voor het ontwerpen van kwantumtechnologieën die gebruikmaken van collectieve fenomenen, zoals kwantumbatterijen, superradiante lasers en kwantummetrologie.
• Fundamenteel Inzicht: Het verduidelijkt hoe realistische dissipatie emergent gedrag zoals synchronisatie en tijd-kristal-dynamica beïnvloedt.
• Toekomstige Richtingen: De framework is uitbreidbaar naar (ultra)sterke koppeling, interne niveaus en niet-Markovse dynamica, wat het een veelzijdig platform maakt voor toekomstig onderzoek in open kwantum-systemen.

Kortom, de paper presenteert een robuuste, schaalbare semi-klassieke methode die de kloof overbrugt tussen kleine, exact oplosbare kwantumsystemen en grote, experimenteel relevante ensembles met complexe dissipatieve interacties.