Phase and coherence retrieval from near- and far-field intensities

Sintesi Tecnica: Recupero di Fase e Coerenza da Intensità in Campo Vicino e Lontano

Definizione del Problema
La quantificazione delle proprietà di coerenza di campi ottici complessi è fondamentale per applicazioni che spaziano dai array di laser ad alta potenza ai sistemi quantistici coerenti. Sebbene gli algoritmi iterativi di recupero della fase, come il framework Gerchberg–Saxton (GS), siano uno standard per la ricostruzione di campi complessi da misure di intensità, essi tipicamente assumono campi completamente coerenti. Molti sistemi fisici importanti, inclusi i grandi array di laser, i fasci di raggi X/elettroni e gli stati quantistici strutturati, esibiscono solo una coerenza spaziale parziale. In questi casi, il campo non può essere descritto da una semplice ampiezza complessa, ma richiede una funzione di intensità mutua quadridimensionale, $J(r, r') = \langle E^*(r)E(r') \rangle$, che rappresenta la coerenza spaziale del primo ordine. I metodi di ricostruzione esistenti spesso si affidano ad assunzioni restrittive (ad esempio, coerenza omogenea, basi modali note o matrici di densità a basso rango) o non riescono a recuperare l'intera intensità mutua da misure di intensità in soli due piani (campo vicino e campo lontano).

Metodologia
Gli autori propongono un nuovo paradigma per la ricostruzione dell'intera intensità mutua di campi parzialmente coerenti a partire da misure di intensità in due piani coniugati di Fourier. L'approccio adatta l'algoritmo GS per propagare e proiettare iterativamente l'intensità mutua tra i domini del campo vicino e del campo lontano. Sono stati sviluppati due implementazioni complementari:

1. Algoritmo Tensor GS: Questo metodo ricostruisce direttamente l'intensità mutua applicando proiezioni alternate al tensore dell'intensità mutua quadridimensionale. Il processo prevede:

   • La propagazione dell'intensità mutua tra i piani utilizzando trasformate di Fourier quadridimensionali.
   • L'applicazione dei vincoli di intensità target mediante la proiezione degli elementi diagonali dell'intensità mutua sulle intensità misurate nel campo vicino e nel campo lontano.
   • L'aggiornamento degli elementi fuori diagonale per preservare la struttura di coerenza sottostante, garantendo al contempo la simmetria ermitiana e la coerenza fisica (assicurando che $|\tilde{g}_1(k, k')| \leq 1$).
   • L'iterazione fino alla convergenza.

2. Variante GS Monte Carlo (MC): Per ridurre il costo computazionale associato alle trasformate di Fourier quadridimensionali, questo metodo rappresenta il campo parzialmente coerente come un insieme di $P$ realizzazioni complesse casuali indipendenti.

   • I campi vengono propagati tra i piani utilizzando trasformate di Fourier bidimensionali standard.
   • Le ampiezze di campo sono moltiplicate per una funzione comune ad ogni passaggio per garantire che l'intensità media dell'insieme corrisponda all'intensità target.
   • L'intensità mutua viene stimata tramite la media dell'insieme sulle $P$ realizzazioni.
   • L'accuratezza scala come $1/\sqrt{P}$, permettendo un'approssimazione controllata. La complessità computazionale scala come $PM \log(M)$ (dove $M$ è il numero di modi spaziali), rispetto a $M^2 \log(M)$ per il metodo Tensor GS.

Contributi Chiave

• Framework di Ricostruzione Generale: Il documento presenta un metodo generale per ricostruire l'intera intensità mutua spaziale da misure di intensità in soli due piani, senza richiedere basi modali predefinite, parametrizzazioni a basso rango o molteplici piani di rilevamento.
• Approcci Algoritmici Duali: L'introduzione sia di un metodo ad alta accuratezza (Tensor GS) che di una variante computazionalmente efficiente (MC GS) fornisce flessibilità per diverse dimensioni di sistema e requisiti di accuratezza.
  o l'efficienza computazionale.
• Gestione dell'Inomogeneità: I metodi sono progettati specificamente per gestire campi parzialmente coerenti spazialmente inomogenei, un regime in cui i tradizionali approcci di trasformata inversa di Fourier (IFT) falliscono.

Risultati
Gli autori hanno validato entrambi i metodi attraverso simulazioni numeriche e dimostrazioni sperimentali:

• Simulazioni Numeriche:
  • Array Lineari: Le simulazioni su array lineari di fino a $N=600$ fasci Gaussiani con una prescritta decadenza della coerenza gaussiana hanno mostrato che, mentre i metodi IFT recuperano accuratamente la coerenza solo per lunghezze di coerenza brevi ($\xi \ll Nd$), i metodi Tensor GS e MC GS mantengono l'accuratezza anche quando la lunghezza di coerenza supera la dimensione dell'array ($\xi > Nd$), dove l'intensità mutua diventa fortemente inomogenea.
  • Array ad Anello: Le simulazioni su array ad anello ($N=50$) hanno dimostrato che il metodo IFT fallisce qualitativamente a causa del breakdown delle assunzioni di invarianza per traslazione. Al contrario, entrambi i metodi proposti hanno recuperato accuratamente la corretta struttura dell'intensità mutua, con il Tensor GS che fornisce una ricostruzione ad alta fedeltà e il MC GS che offre un'approssimazione computazionalmente efficiente.
• Dimostrazione Sperimentale:
  • Il metodo Tensor GS è stato applicato a un array triangolare di 130 laser accoppiati con accoppiamento controllato e chiralità scalata.
  • Il sistema ha generato stati con decadimento della coerenza mutua non monotonico e fasi non banali (anti-vortice scalato e sovrapposizione uguale di stati vortice/anti-vortice).
  • La ricostruzione ha ottenuto un eccellente accordo con la teoria, riproducendo i decadimenti di fase attesi (ad esempio, $2\pi/3$) e il decadimento gaussiano della magnitudo della coerenza.
  • L'errore di fase RMS è stato misurato a $2\pi/250$ per lo stato anti-vortice scalato e $2\pi/65$ per lo stato di sovrapposizione uguale.

Significatività e Rivendicazioni
Il documento afferma che questi risultati indicano che l'approccio è praticabile per la caratterizzazione di sistemi ottici parzialmente coerenti su larga scala. I metodi completano gli approcci tomografici esistenti eliminando la necessità di una conoscenza pregressa della struttura di coerenza o di molteplici piani di rilevamento. Gli autori osservano che, sebbene la unicità teorica della ricostruzione dell'intensità mutua da due piani coniugati di Fourier rimanga una questione aperta, i loro risultati numerici e sperimentali dimostrano una convergenza stabile verso soluzioni fisicamente significative per strutture di coerenza fortemente inomogenee e una robustezza al rumore sperimentale.

Il lavoro suggerisce potenziali estensioni future, come l'incorporazione della propagazione in spazio libero per ricostruire i campi da piani non coniugati di Fourier e l'applicazione del framework alla tomografia quantistica, data l'analogia formale tra la ricostruzione dell'intensità mutua e la ricostruzione della matrice di densità. Tuttavia, il documento mantiene una posizione modesta, concentrandosi sulla capacità dimostrata di ricostruire fase e coerenza in specifici e impegnativi sistemi inomogenei senza estendere eccessivamente le rivendicazioni a campi ottici arbitrari.