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

Technische Zusammenfassung: Phasen- und Kohärenzrekonstruktion aus Nahfeld- und Fernfeldintensitäten

Problemstellung
Die Quantifizierung der Kohärenzeigenschaften komplexer optischer Felder ist entscheidend für Anwendungen, die von hochleistungsfähigen Laserarrays bis hin zu kohärenten Quantensystemen reichen. Während iterative Phasenrekonstruktionsalgorithmen, wie das Gerchberg–Saxton (GS)-Framework, der Standard für die Rekonstruktion komplexer Felder aus Intensitätsmessungen sind, setzen sie typischerweise voll kohärente Felder voraus. Viele wichtige physikalische Systeme, darunter große Laserarrays, Röntgen-/Elektronenstrahlen und strukturierte Quantenzustände, weisen jedoch nur eine partielle räumliche Kohärenz auf. In diesen Fällen kann das Feld nicht durch eine einfache komplexe Amplitude beschrieben werden, sondern erfordert eine vierdimensionale Mutual-Intensitätsfunktion, $J(r, r') = \langle E^*(r)E(r') \rangle$, welche die erstklassige räumliche Kohärenz repräsentiert. Bestehende Rekonstruktionsmethoden stützen sich oft auf restriktive Annahmen (z. B. homogene Kohärenz, bekannte Modenbasen oder niedrigdimensionale Dichtematrizen) oder scheitern an der Rekonstruktion der vollständigen Mutual-Intensität aus Intensitätsmessungen in nur zwei Ebenen (Nahfeld und Fernfeld).

Methodik
Die Autoren schlagen ein neues Paradigma zur Rekonstruktion der vollständigen Mutual-Intensität partiell kohärenter Felder aus Intensitätsmessungen in zwei Fourier-konjugierten Ebenen vor. Der Ansatz passt den GS-Algorithmus an, um die Mutual-Intensität iterativ zwischen der Nahfeld- und der Fernfeldebene zu propagieren. Es werden zwei komplementäre Implementierungen entwickelt:

1. Tensor-GS-Algorithmus: Diese Methode rekonstruiert die Mutual-Intensität direkt durch Anwendung alternierender Projektionen auf den vierdimensionalen Mutual-Intensitäts-Tensor. Der Prozess umfasst:

   • Die Propagation der Mutual-Intensität zwischen den Ebenen mittels vierdimensionaler Fourier-Transformationen.
   • Das Erzwingen von Zielintensitäts-Constraints durch Projektion der Diagonalelemente der Mutual-Intensität auf die gemessenen Nahfeld- und Fernfeldintensitäten.
   • Die Aktualisierung der Off-Diagonal-Elemente, um die zugrunde liegende Kohärenzstruktur zu bewahren und gleichzeitig die Hermitesche Symmetrie sowie die physikalische Konsistenz (Sicherstellung von $|\tilde{g}_1(k, k')| \leq 1$) zu erzwingen.
   • Iteration bis zur Konvergenz.

2. Monte-Carlo (MC) GS-Variante: Um die Rechenkosten im Zusammenhang mit vierdimensionalen Fourier-Transformationen zu reduzieren, stellt diese Methode das partiell kohärente Feld als ein Ensemble von $P$ unabhängigen, zufälligen komplexen Realisierungen dar.

   • Die Felder werden zwischen den Ebenen unter Verwendung von Standard-zweidimensionalen Fourier-Transformationen propagiert.
   • Die Feldamplituden werden in jedem Schritt mit einer gemeinsamen Funktion multipliziert, um sicherzustellen, dass die Ensemble-gemittelte Intensität der Zielintensität entspricht.
   • Die Mutual-Intensität wird durch Mittelung über das Ensemble der $P$ Realisierungen geschätzt.
   • Die Genauigkeit skaliert als $1/\sqrt{P}$, was eine kontrollierte Approximation ermöglicht. Die Komplexität skaliert als $PM \log(M)$ (wobei $M$ die Anzahl der räumlichen Moden ist), im Vergleich zu $M^2 \log(M)$ für die Tensor-GS-Methode.

Wesentliche Beiträge

• Allgemeines Rekonstruktions-Framework: Die Arbeit präsentiert eine allgemeine Methode zur Rekonstruktion der vollständigen räumlichen Mutual-Intensität aus Intensitätsmessungen in nur zwei Ebenen, ohne die Notwendigkeit einer vordefinierten Modenbasis, einer niedrigdimensionalen Parametrisierung oder mehrerer Detektionsebenen.
• Duale algorithmische Ansätze: Die Einführung sowohl eines hochpräzisen Tensor-GS-Verfahrens als auch einer recheneffizienten MC-GS-Variante bietet Flexibilität für unterschiedliche Systemgrößen und Genauigkeitsanforderungen.
• Umgang mit Inhomogenität: Die Methoden sind speziell für räumlich inhomogene partiell kohärente Felder konzipiert, ein Regime, in dem traditionelle Inverse-Fourier-Transformations-Ansätze (IFT) versagen.

Ergebnisse
Die Autoren validierten beide Methoden durch numerische Simulationen und experimentelle Demonstrationen:

• Numerische Simulationen:
  • Lineare Arrays: Simulationen an linearen Arrays von bis zu $N=600$ Gaußschen Strahlen mit einer vorgegebenen Gaußschen Abkling-Kohärenz zeigten, dass IFT-Methoden die Kohärenz zwar für kurze Kohärenzlängen ($\xi \ll Nd$) korrekt rekonstruieren, sowohl die Tensor-GS- als auch die MC-GS-Methoden jedoch eine Genauigkeit beibehalten, selbst wenn die Kohärenzlänge die Array-Größe überschreitet ($\xi > Nd$), wo die Mutual-Intensität stark inhomogen wird.
  • Ring-Arrays: Simulationen an Ring-Arrays ($N=50$) zeigten, dass die IFT-Methode aufgrund des Zusammenbruchs der Annahme der Translationsinvarianz qualitativ versagt. Im Gegensatz dazu rekonstruierten beide vorgeschlagenen Methoden die korrekte Struktur der Mutual-Intensität präzise, wobei die Tensor-GS-Methode eine hochtreue Rekonstruktion lieferte, während die MC-GS-Methode eine recheneffiziente Approximation bot.
• Experimentelle Demonstration:
  • Die Tensor-GS-Methode wurde auf ein dreieckiges Array aus 130 gekoppelten Lasern mit kontrollierter Kopplung und gestaffelter Chiralität angewendet.
  • Das System erzeugte Zustände mit nicht-monotonem Zerfall der Mutual-Kohärenz und nicht-trivialen Phasen (gestaffelte Anti-Vortex- und Gleichgewichtssuperposition von Vortex/Anti-Vortex-Zuständen).
  • Die Rekonstruktion zeigte eine exzellente Übereinstimmung mit der Theorie und reproduzierte die erwarteten Phasenwindungen (z. B. $2\pi/3$) sowie den Gaußschen Zerfall der Kohärenzamplitude.
  • Der quadratische Mittelwert (RMS) des Phasenfehlers wurde mit $2\pi/250$ für den gestaffelten Anti-Vortex-Zustand und $2\pi/65$ für den Gleichgewichtszustand gemessen.

Bedeutung und Ansprüche
Die Arbeit behauptet, dass diese Ergebnisse darauf hindeuten, dass der Ansatz praktikabel für die Charakterisierung großer partiell kohärenter optischer Systeme ist. Die Methoden ergänzen bestehende tomographische Ansätze, indem sie die Kenntnis der Kohärenzstruktur oder mehrerer Detektionsebenen überflüssig machen. Die Autoren merken an, dass die theoretische Eindeutigkeit der Rekonstruktion der Mutual-Intensität aus zwei Fourier-konjugierten Ebenen eine offene Frage bleibt, ihre numerischen und experimentellen Ergebnisse jedoch eine stabile Konvergenz hin zu physikalisch sinnvollen Lösungen für stark inhomogene Kohärenzstrukturen und eine Robustheit gegenüber experimentellem Rauschen demonstrieren.

Die Arbeit deutet auf potenzielle zukünftige Erweiterungen hin, wie etwa die Einbeziehung der Freiraumpropagation zur Rekonstruktion von Feldern aus nicht-Fourier-konjugierten Ebenen sowie die Anwendung des Frameworks auf die Quantentomographie, angesichts der formalen Analogie zwischen der Rekonstruktion der Mutual-Intensität und der Rekonstruktion der Dichtematrix. Die Autoren wahren jedoch eine bescheidene Haltung und konzentrieren sich auf die demonstrierte Fähigkeit, Phase und Kohärenz in spezifischen, anspruchsvollen inhomogenen Systemen ohne Überdehnung der Ansprüche auf beliebige optische Felder zu rekonstruieren.