WGQED1D: a microscopic N-atom waveguide-QED toolkit for single-photon scattering and weak-measurement experiments

Technische Samenvatting van WGQED1D

Probleemstelling
De propagatie van zwak licht door een atomaire ensemble is een fundamenteel probleem in de kwantumoptica, relevant voor fenomenen zoals elektromagnetisch geïnduceerde transparantie, slow-light propagatie en negatieve post-geselecteerde verblijftijden. Hoewel macroscopische lineaire-respons-theorie (bijv. Beer–Lambert absorptie en Kramers–Kronig dispersie) goed gevestigd is, vereist het overbruggen van de kloof tussen de microscopische dynamica van $N$ individuele atomen en de geïntegreerde optische respons van een wolk een zorgvuldige input–output-behandeling. De auteurs merken op dat directe implementaties van standaardformules vaak kwantitatieve afwijkingen vertonen met de Beer–Lambert-limiet van de orde grootte van eenheden. Deze discrepanties worden toegeschreven aan subtiele conventieke keuzes met betrekking tot de bidirectionele structuur van golfgeleiderkoppeling en de kalibratie tussen optische diepte en koppelingsefficiëntie per atoom.

Methodologie
Het artikel introduceert WGQED1D, een lichtgewicht, gevalideerde Python-toolkit ontworpen voor de microscopische simulatie van enkelvoudige foton-verstrooiing door $N$ twee-niveau atomen in een eendimensionale bidirectionele golfgeleider. De methodologie steunt op de volgende kerncomponenten:

1. Fysisch Model: Het systeem wordt gemodelleerd in de enkelvoudige excitatie-subruimte, gerechtvaardigd door enkelvoudige foton-inputs en lage verzadiging. De dynamica wordt beheerst door een standaard niet-Hermitische effectieve Hamiltoniaan ($H_{\text{eff}}$) werkend op $N$ complexe amplitudes ($b_j$). Deze Hamiltoniaan omvat:
   • Collectieve dipool–dipool interacties gemedieerd door geretardeerde foton-uitwisseling ($e^{ik_0|z_j-z_k|}$).
   • Diagonale verlies naar de vrije ruimte.
   • Een drifterm die de input Gaussische golfpakket representeert.
2. Input–Output Relatie: De toolkit reconstrueert het outputveld ($\phi_{\text{out}}$) uit atomaire amplitudes met behulp van een specifieke input–output relatie die rekening houdt met de fase-accumulatie van atomen ten opzichte van de wolk-grens.
3. Kritische Kalibratie: Een centrale methodologische bijdrage is de identificatie en verificatie van een specifieke kalibratie tussen de koppelingsefficiëntie per atoom in de golfgeleider ($\beta = \Gamma_{1D}/\Gamma$) en de doel-optische diepte ($OD$). De auteurs stellen vast dat kwantitatieve overeenstemming met de vrije-ruimte Beer–Lambert-limiet vereist dat $\beta = OD/(2N)$. Deze factor twee komt voort uit de bidirectionele aard van de golfgeleider, waarbij atomen met gelijke waarschijnlijkheid voorwaarts en achterwaarts emitteren, in tegenstelling tot unidirectionele ("chirale") idealisaties.
4. Numerieke Implementatie: De toolkit gebruikt deterministische evolutie via een expliciete Runge–Kutta 4(5) methode. Het opereert in $O(N^2)$ geheugen en operaties per tijdstap, gebruikmakend van dense lineaire algebra. Het bevat spectrale analyse-instrumenten om de post-geselecteerde atomaire excitatietijd ($\tau_T$) te extraheren met behulp van de Wiseman–Steinberg weak-value formule.

Belangrijkste Bijdragen

• WGQED1D Toolkit: Een reproduceerbare, open-source Python-module (~400 regels) die het microscopische model, de input–output reconstructie en de spectrale weak-value analyse implementeert. Het is afhankelijk van slechts standaard wetenschappelijke bibliotheken (NumPy, SciSciPy, Matplotlib).
• Microscopische Kalibratie: De expliciete afleiding en numerieke verificatie van de relatie $\beta = OD/(2N)$, die systematische kwantitatieve onenigheid tussen microscopische simulaties en macroscopische Beer–Lambert voorspellingen oplost.
• Validatie-infrastructuur: Het artikel documenteert drie verschillende validaties:
  1. Overeenstemming met de analytische pulse-gemiddelde Beer–Lambert transmissie over optische dieptes $OD \in [0.25, 3]$ met fouten $<0.1\%$.
  2. Verificatie van de Wiseman–Steinberg weak-value formule voor post-geselecteerde excitatietijden.
  3. Directe vergelijking met experimentele data van Angulo et al. (2026) betreffende negatieve atomaire excitatietijden.

Resultaten

• Beer–Lambert Overeenstemming: Simulaties met $N=500$ atomen reproduceren de analytische Beer–Lambert transmissie tot beter dan 0,1% over de geteste range van optische dieptes, wat de geldigheid van de $\beta = OD/(2N)$ kalibratie in de ijle limiet bevestigt.
• Weak-Value Extractie: De toolkit extraheert succesvol post-geselecteerde excitatietijden ($\tau_T$) die consistent zijn met de analytische weak-value formule. Discrepanties liggen binnen 10% en worden toegeschreven aan de eindige spectrale resolutie en collectieve finite-$N$ correcties bij hoge optische dieptes.
• Experimentele Vergelijking: Bij vergelijking met de data van Angulo et al. over negatieve post-geselecteerde verblijftijden in koude rubidiumwolken, reproduceert de toolkit drie van de zes datapunten binnen hun $1\sigma$ foutmarges. Voor de resterende punten (voornamelijk bij $OD=4$) volgt de simulatie de kwalitatieve trend en komt deze in sommige gevallen nauwer overeen met experimentele metingen dan de lineaire-respons analytische formule.
• Prestaties: De toolkit draait efficiënt op standaard hardware en produceert geconvergeerde resultaten voor realistische ensembles ($N \sim 200–1000$) binnen seconden tot minuten.

Betekenis en Claims
Het artikel positioneert WGQED1D als een "reproduceerbaar numeriek platform" voor het bestuderen van weak measurement, anomalieën in dispersie, group delay, en post-geselecteerde observabelen in enkelvoudige foton transmissie door optisch dikke atomaire wolken. De auteurs benadrukken dat de toolkit prioriteit geeft aan numerieke betrouwbaarheid en verifieerbaarheid boven ruwe prestatie, en biedt een directe regel-voor-regel match met de leidende vergelijkingen.

De betekenis van het werk ligt in het bieden van een gevalideerde tool die de subtiliteiten van bidirectionele golfgeleiderkoppeling correct afhandelt, waardoor onderzoekers in staat zijn om microscopische simulaties uit te voeren die kwantitatief overeenkomen met macroscopische limieten en recente experimentele observaties. De auteurs geven expliciet aan dat voor gebruikers die grotere systeemgroottes vereisen ($N \gtrsim 10^4$), de algoritmische structuur kan worden uitgebreid naar sparse lineaire algebra of GPU-backends, hoewel deze momenteel niet zijn geïmplementeerd. De toolkit wordt uitgebracht onder de MIT-licentie om open wetenschappelijk onderzoek te faciliteren.