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

Resumo Técnico: Recuperação de Fase e Coerência a partir de Intensidades de Campo Próximo e Campo Distante

Definição do Problema
A quantificação das propriedades de coerência de campos ópticos complexos é crítica para aplicações que variam desde arranjos de lasers de alta potência até sistemas quânticos coerentes. Embora algoritmos de recuperação de fase iterativos, como o arcabouço de Gerchberg–Saxton (GS), sejam o padrão para reconstruir campos complexos a partir de medições de intensidade, eles tipicamente assumem campos totalmente coerentes. Muitos sistemas físicos importantes, incluindo grandes arranjos de lasers, feixes de raios X/elétrons e estados quânticos estruturados, exibem apenas coerência espacial parcial. Nesses casos, o campo não pode ser descrito por uma simples amplitude complexa, mas requer uma função de intensidade mútua quadridimensional, $J(r, r') = \langle E^*(r)E(r') \rangle$, que representa a coerência espacial de primeira ordem. Métodos de reconstrução existentes frequentemente dependem de suposições restritivas (ex: coerência homogênea, bases modais conhecidas ou matrizes de densidade de baixo posto), ou falham em recuperar a intensidade mútua completa a partir de medições de intensidade em apenas dois planos (campo próximo e campo distante).

Metodologia
Os autores propõem um novo paradigma para reconstruir a intensidade mútua completa de campos parcialmente coerentes a partir de medições de intensidade em dois planos Fourier-conjugados. A abordagem adapta o algoritmo GS para propagar e projetar iterativamente a intensidade mútua entre os domínios de campo próximo e campo distante. Duas implementações complementares foram desenvolvidas:

1. Algoritmo Tensor GS: Este método reconstrói diretamente a intensidade mútua aplicando projeções alternadas ao tensor de intensidade mútua quadridimensional. O processo envolve:

   • Propagar a intensidade mútua entre os planos usando transformadas de Fourier quadridimensionais.
   • Impor restrições de intensidade alvo projetando os elementos diagonis da intensidade mútua sobre as intensidades medidas no campo próximo e no campo distante.
   • Atualizar os elementos fora da diagonal para preservar a estrutura de coerência subjacente, garantindo a simetria hermitiana e a consistência física (garantindo que $|\tilde{g}_1(k, k')| \leq 1$).
   • Iterar até a convergência.

2. Variante MC GS de Monte Carlo: Para reduzir o custo computacional associado às transformadas de Fourier quadridimensionais, este método representa o campo parcialmente coerente como um conjunto de $P$ realizações complexas aleatórias independentes.

   • Os campos são propagados entre os planos usando transformadas de Fourier bidimensionais padrão.
   • As amplitudes de campo são multiplicadas por uma função comum em cada etapa para garantir que a intensidade média do conjunto corresponda à intensidade alvo.
   • A intensidade mútua é estimada por meio da média do conjunto sobre as $P$ realizações.
   • A precisão escala como $1/\sqrt{P}$, permitindo uma aproximação controlada. A complexidade computacional escala como $PM \log(M)$ (onde $M$ é o número de modos espaciais), comparada a $M^2 \log(M)$ para o método Tensor GS.

Principais Contribuições

• Arcabouço de Reconstrução Geral: O artigo apresenta um método geral para reconstruir a intensidade mútua espacial completa a partir de medições de intensidade em apenas dois planos, sem exigir bases modais predefinidas, parametrização de baixo posto ou múltiplos planos de detecção.
• Abordagens Algorítmicas Duais: A introdução de ambos o método de alta precisão Tensor GS e a variante computacionalmente eficiente MC GS proporciona flexibilidade para diferentes tamanhos de sistema e requisitos de precisão.
• Tratamento de Inhomogeneidade: Os métodos são especificamente projetados para lidar com campos parcialmente coerentes espacialmente inhomogêneos, um regime onde as abordagens tradicionais de transformada inversa de Fourier (IFT) falham.

Resultados
Os autores validaram ambos os métodos através de simulações numéricas e demonstração experimental:

• Simulações Numéricas:
  • Arranjos Lineares: Simulações em arranjos lineares de até $N=600$ feixes Gaussianos com decaimento Gaussiano prescrito de coerência mostraram que, enquanto os métodos IFT recuperam a coerência com precisão apenas para comprimentos de coerência curtos ($\xi \ll Nd$), os métodos Tensor GS e MC GS mantêm a precisão mesmo quando o comprimento de coerência excede o tamanho do arranjo ($\xi > Nd$), onde a intensidade mútua torna-se fortemente inhomogênea.
  • Arranjos em Anel: Simulações em arranjos em anel ($N=50$) demonstraram que o método IFT falha qualitativamente devido à quebra das suposições de invariância de translação. Em contraste, ambos os métodos propostos recuperaram com precisão a estrutura correta da intensidade mútua, com o Tensor GS fornecendo uma reconstrução de alta fidelidade e o MC GS oferecendo uma aproximação computacionalmente eficiente.
• Demonstração Experimental:
  • O método Tensor GS foi aplicado a um arranjo triangular de 130 lasers acoplados com acoplamento controlado e quiralidade alternada.
  • O sistema gerou estados com decaimento de coerência mútua não monotônico e fases não triviais (anti-vórtex alternado e superposição igual de estados de vórtice/anti-vórtex).
  • A reconstrução alcançou excelente concordância com a teoria, reproduzindo as rotações de fase esperadas (ex: $2\pi/3$) e o decaimento Gaussiano da magnitude de coerência.
  • O erro de fase RMS (raiz-média quadrática) foi medido em $2\pi/250$ para o estado anti-vórtex alternado e $2\pi/65$ para o estado de superposição igual.

Significância e Alegações
O artigo afirma que estes resultados indicam que a abordagem é prática para caracterizar sistemas ópticos parcialmente coerentes de grande escala. Os métodos complementam as abordagens tomográficas existentes ao remover a necessidade de conhecimento prévio da estrutura de coerência ou de múltiplos planos de detecção. Os autores observam que, embora a unicidade teórica da reconstrução da intensidade mútua a partir de dois planos Fourier-conjugados permaneça uma questão aberta, seus resultados numéricos e experimentais demonstram convergência estável para soluções fisicamente significativas para estruturas de coerência fortemente inhomogêneas e robustez ao ruído experimental.

O trabalho sugere potenciais extensões futuras, como a incorporação de propagação em espaço livre para reconstruir campos a partir de planos não Fourier-conjugados e a aplicação do arcabouço à tomografia quântica, dada a analogia formal entre a reconstrução da intensidade mútua e a reconstrução da matriz de densidade. No entanto, o artigo mantém uma postura modesta, focando na capacidade demonstrada de reconstruir fase e coerência em sistemas inhomogêneos específicos e desafiadores, sem estender excessivamente as alegações para campos ópticos arbitrários.