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

Resumen Técnico: Recuperación de Fase y Coherencia a partir de Intensidades en Campo Cercano y Lejano

Planteamiento del Problema
La cuantificación de las propiedades de coherencia de campos ópticos complejos es crítica para aplicaciones que van desde matrices de láseres de alta potencia hasta sistemas cuánticos coherentes. Si bien los algoritmos de recuperación de fase iterativos, como el marco de Gerchberg–Saxton (GS), son estándares para reconstruir campos complejos a partir de mediciones de intensidad, estos típicamente asumen campos totalmente coherentes. Muchos sistemas físicos importantes, incluyendo grandes matrices de láseres, haces de rayos X/electrones y estados cuánticos estructurados, exhiben solo una coherencia espacial parcial. En estos casos, el campo no puede describirse mediante una simple amplitud compleja, sino que requiere una función de intensidad mutua tetradimensional, $J(r, r') = \langle E^*(r)E(r') \rangle$, que representa la coherencia espacial de primer orden. Los métodos de reconstrucción existentes suelen depender de suposiciones restrictivas (por ejemplo, coherencia homogénea, bases modales conocidas o matrices de densidad de bajo rango) o fallan al intentar recuperar la intensidad mutua completa a partir de mediciones de intensidad en solo dos planos (campo cercano y campo lejano).

Metodología
Los autores proponen un nuevo paradigma para la reconstrucción de la intensidad mutua completa de campos parcialmente coherentes a partir de mediciones de intensidad en dos planos de Fourier conjugados. El enfoque adapta el algoritmo GS para propagar e iterar la proyección de la intensidad mutua entre los dominios de campo cercano y campo lejano. Se desarrollan dos implementaciones complementarias:

1. Algoritmo Tensor GS: Este método reconstruye directamente la intensidad mutua aplicando proyecciones alternas al tensor de intensidad mutua tetradimensional. El proceso implica:

   • Propagar la intensidad mutua entre los planos utilizando transformadas de Fourier tetradimensionales.
   • Imponer restricciones de intensidad objetivo mediante la proyección de los elementos diagonales de la intensidad mutua sobre las intensidades medidas en el campo cercano y el campo lejano.
   • Actualizar los elementos fuera de la diagonal para preservar la estructura de coherencia subyacente mientras se asegura la simetría hermítica y la consistencia física (asegurando que $|\tilde{g}_1(k, k')| \leq 1$).
   • Iterar hasta la convergencia.

2. Variante de Monte Carlo (MC) GS: Para reducir el costo computacional asociado con las transformadas de Fourier tetradimensionales, este método representa el campo parcialmente coherente como un conjunto de $P$ realizaciones complejas aleatorias independientes.

   • Los campos se propagan entre los planos utilizando transformadas de Fourier bidimensionales estándar.
   • Las amplitudes de campo se multiplican por una función común en cada paso para asegurar que la intensidad promedio del conjunto coincida con la intensidad objetivo.
   • La intensidad mutua se estima mediante el promedio del conjunto sobre las $P$ realizaciones.
   • La precisión escala como $1/\sqrt{P}$, lo que permite una aproximación controlada. La complejidad computacional escala como $PM \log(M)$ (donde $M$ es el número de modos espaciales), en comparación con $M^2 \log(M)$ para el método Tensor GS.

Contribuciones Clave

• Marco de Reconstrucción General: El artículo presenta un método general para reconstruir la intensidad mutua espacial completa a partir de mediciones de intensidad en solo dos planos, sin requerir bases modales predefinidas, parametrización de bajo rango o múltiples planos de detección.
• Enfoques Algorítmicos Duales: La introducción de tanto un método de alta precisión (Tensor GS) como un método computacionalmente eficiente (MC GS) proporciona flexibilidad para diferentes tamaños de sistema y requisitos de precisión.
• Manejo de la Inhomogeneidad: Los métodos están diseñados específicamente para manejar campos parcialmente coherentes espacialmente inhomogéneos, un régimen donde los enfoques tradicionales de transformada inversa de Fourier (IFT) fallan.

Resultados
Los autores validaron ambos métodos mediante simulaciones numéricas y demostración experimental:

• Simulaciones Numéricas:
  • Matrices Lineales: Las simulaciones en matrices lineales de hasta $N=600$ haces gaussianos con una coherencia prescrita de decaimiento gaussiano mostraron que, mientras los métodos de IFT recuperan la coherencia con precisión solo para longitudes de coherencia cortas ($\xi \ll Nd$), los métodos Tensor GS y MC GS mantienen la precisión incluso cuando la longitud de coherencia excede el tamaño de la matriz ($\xi > Nd$), donde la intensidad mutua se vuelve fuertemente inhomogénea.
  • Matrices de Anillo: Las simulaciones en matrices de anillo ($N=50$) demostraron que el método IFT falla cualitativamente debido a la ruptura de las suposiciones de invarianza de traslación. En contraste, ambos métodos propuestos recuperaron con precisión la estructura correcta de la intensidad mutua, siendo el Tensor GS el que proporcionó una reconstrucción de alta fidelidad y el MC GS una aproximación computacionalmente eficiente.
• Demostración Experimental:
  • El método Tensor GS se aplicó a una matriz triangular de 130 láseres acoplados con acoplamiento controlado y quiralidad escalonada.
  • El sistema generó estados con decaimiento de coherencia mutua no monotónico y fases no triviales (anti-vórtice escalonado y superposición igual de estados de vórtice/anti-vórtice).
  • La reconstrucción logró una excelente concordancia con la teoría, reproduciendo los giros de fase esperados (por ejemplo, $2\pi/3$) y el decaimiento gaussiano de la magnitud de la coherencia.
  • El error de fase raíz cuadrática media (RMS) se midió en $2\pi/250$ para el estado de anti-vórtice escalonado y $2\pi/65$ para el estado de superposición igual.

Significancia y Reivindicaciones
El artículo sostiene que estos resultados indican que el enfoque es práctico para caracterizar sistemas ópticos parcialmente coherentes de gran escala. Los métodos complementan los enfoques tomográficos existentes al eliminar la necesidad de conocimiento previo de la estructura de coherencia o de múltiples planos de detección. Los autores señalan que, si bien la unicidad teórica de la reconstrucción de la intensidad mutua a partir de dos planos de Fourier conjugados sigue siendo una cuestión abierta, sus resultados numéricos y experimentales demuestran una convergencia estable hacia soluciones físicamente significativas para estructuras de coherencia fuertemente inhomogéneas y robustez ante el ruido experimental.

El trabajo sugiere posibles extensiones futuras, como la incorporación de la propagación en espacio libre para reconstruir campos a partir de planos no conjugados de Fourier y la aplicación del marco a la tomografía cuántica, dada la analogía formal entre la reconstrucción de la intensidad mutua y la reconstrucción de la matriz de densidad. Sin embargo, el artículo mantiene una postura modesta, centrándose en la capacidad demostrada para reconstruir la fase y la coherencia en sistemas inhomogéneos específicos y desafiantes, sin extender excesivamente sus afirmaciones a campos ópticos arbitrarios.