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

Sintesi Tecnica di WGQED1D

Definizione del Problema
La propagazione di luce debole attraverso un insieme atomico è un problema fondamentale nell'ottica quantistica, rilevante per fenomeni quali la trasparenza indotta elettromagneticamente, la propagazione della luce lenta e i tempi di permanenza negativi post-selezionati. Sebbene la teoria della risposta lineare macroscopica (ad esempio, l'assorbimento di Beer–Lambert e la dispersione di Kramers–Kronig) sia ben consolidata, colmare il divario tra la dinamica microscopica di $N$ singoli atomi e la risposta ottica integrata di una nuvola richiede un attento trattamento di input–output. Gli autori osservano che le implementazioni dirette delle formule standard spesso producono disaccordi quantitativi con il limite di Beer–Lambert con fattori di ordine unitario. Tali discrepanze sono attribuite a sottili scelte di convenzione riguardanti la struttura bidirezionale dell'accoppiamento in guida d'onda e alla calibrazione tra l'ottica di profondità e l'efficienza di accoppiamento per atomo.

Metodologia
Il documento introduce WGQED1D, un toolkit Python leggero e validato, progettato per la simulazione microscopica dello scattering di un singolo fotone da parte di $N$ atomi a due livelli in una guida d'onda unidimensionale bidirezionale. La metodologia si basa sui seguenti componenti principali:

1. Modello Fisico: Il sistema è modellato nel sottospazio a singola eccitazione, giustificato da input a singolo fotone e bassa saturazione. La dinamica è governata da un'Hamiltoniana efficace non hermitiana standard ($H_{\text{eff}}$) che agisce su $N$ ampiezze complesse ($b_j$). Questa Hamiltoniana include:
   • Interazioni collettive dipolo–dipolo mediate dallo scambio di fotoni ritardato ($e^{ik_0|z_j-z_k|}$).
   • Perdita diagonale verso lo spazio libero.
   • Un termine di guida che rappresenta il pacchetto d'onda gaussiano in ingresso.
2. Relazione Input–Output: Il toolkit ricostruisce il campo in uscita ($\phi_{\text{out}}$) partendo dalle ampiezze atomiche utilizzando una specifica relazione input–output che tiene conto dell'accumulo di fase degli atomi rispetto al confine della nuvola.
3. Calibrazione Critica: Un contributo metodologico centrale è l'identificazione e la verifica di una specifica calibrazione tra l'efficienza di accoppiamento in guida d'onda per atomo ($\beta = \Gamma_{1D}/\Gamma$) e la target ottica profondità ($OD$). Gli autori stabiliscono che l'accordo quantitativo con il limite di Beer–Lambert in spazio libero richiede $\beta = OD/(2N)$. Questo fattore due deriva dalla natura bidirezionale della guida d'onda, dove gli atomi emettono in avanti e all'indietro con uguale probabilità, contrastando con le idealizzazioni unidirezionali ("chirali").
4. Implementazione Numerica: Il toolkit utilizza l'evoluzione deterministica tramite un metodo esplicito Runge–Kutta 4(5). Opera con una complessità di $O(N^2)$ in memoria e operazioni per timestep, utilizzando l'algebra lineare densa. Include strumenti di analisi spettrale per estrarre il tempo di eccitazione atomica post-selezionato ($\tau_T$) utilizzando la formula del weak-value di Wiseman–Steinberg.

Contributi Chiave

• Toolkit WGQED1D: Un modulo Python riproducibile e open-source (~400 righe) che implementa il modello microscopico, la ricostruzione input–output e l'analisi del weak-value spettrale. Dipende solo da librerie scientifiche standard (NumPy, SciPy, Matplotlib).
• Calibrazione Microscopica: La derivazione esplicita e la verifica numerica della relazione $\beta = OD/(2N)$, che risolve i sistematici disaccordi quantitativi tra le simulazioni microscopiche e le predizioni macroscopiche di Beer–Lambert.
• Infrastruttura di Validazione: Il documento documenta tre distinte validazioni:
  1. Accordo con la trasmissione analitica di Beer–Lambert mediata sull'impulso attraverso ottiche profondità $OD \in [0.25, 3]$ con errori $<0.1\%$.
  2. Verifica della formula del weak-value di Wiseman–Steinberg per i tempi di eccitazione post-selezionati.
  3. Confronto diretto con i dati sperimentali di Angulo et al. (2026) riguardanti i tempi di eccitazione atomica negativi.

Risultati

• Accordo con Beer–Lambert: Le simulazioni con $N=500$ atomi riproducono la trasmissione analitica di Beer–Lambert con una precisione superiore allo 0.1% attraverso l'intervallo di ottica profondità testato, confermando la validità della calibrazione $\beta = OD/(2N)$ nel limite diluito.
• Estrazione del Weak-Value: Il toolkit estrae con successo i tempi di eccitazione post-selezionati ($\tau_T$) coerentemente con la formula analitica del weak-value. Le discrepanze sono entro il 10% e sono attribuite alla risoluzione spettrale finita e alle correzioni collettive di finite-$N$ ad alte ottiche profondità.
• Confronto Sperimentale: Quando confrontato con i dati di Angulo et al. sui tempi di permanenza post-selezionati negativi in nubi di rubidio freddo, il toolkit riproduce tre dei sei punti dati entro i loro errori $1\sigma$. Per i punti rimanenti (principalmente a $OD=4$), la simulazione segue il trend qualitativo e, in alcuni casi, si allinea più strettamente alle misurazioni sperimentali rispetto alla formula della risposta lineare.
• Prestazioni: Il toolkit gira efficientemente su hardware standard, producendo risultati convergenti per ensemble realistici ($N \sim 200–1000$) in secondi o minuti.

Significatività e Rivendicazioni
Il documento posiziona WGQED1D come una "piattaforma numerica riproducibile" per studiare la misura debole, la dispersione anomala, il ritardo di gruppo e gli osservabili post-selezionati nella trasmissione di un singolo fotone attraverso nubi atomiche con alta ottica profondità. Gli autori sottolineano che il toolkit privilegia l'affidabilità numerica e la verificabilità rispetto alle prestazioni pure, offrendo una corrispondenza riga per riga con le equazioni governanti.

La significatività del lavoro risiede nel fornire uno strumento validato che gestisce correttamente le sottigliezze dell'accoppiamento in guida d'onda bidirezionale, consentendo così ai ricercatori di eseguire simulazioni microscopiche che corrispondano quantitativamente ai limiti macroscopici e alle recenti osservazioni sperimentali. Gli autori dichiarano esplicitamente che per utenti che richiedono dimensioni di sistema maggiori ($N \gtrsim 10^4$), la struttura algoritmica può essere estesa all'algebra lineare sparsa o a backend GPU, sebbene questi non siano attualmente implementati. Il toolkit è rilasciato sotto licenza MIT per facilitare l'indagine scientifica aperta.