Dynamic Coherent Diffractive Imaging Using Only a Support Constraint in the Complex Plane

Resumen Técnico: Imagen Difractiva Coherente Dinámica Utilizando Solo una Restricción de Soporte en el Plano Complejo

Enunciado del Problema
La caracterización cuantitativa de muestras delgadas, transparentes y variables en el tiempo es un desafío persistente en la metrología. Los métodos interferométricos tradicionales requieren entornos estables y múltiples exposiciones, mientras que la Imagen Difractiva Coherente (CDI) estándar asume típicamente un objeto estacionario. Extender la CDI a muestras dinámicas suele necesitar priores complejos, como estructuras de referencia diseñadas, modelos de movimiento explícitos o restricciones aprendidas basadas en datos. Este artículo aborda el problema de reconstruir un objeto de fase variable en el tiempo a partir de una película de difracción de campo cercano utilizando priores mínimos, evitando específicamente la necesidad de desenrollado de fase explícito incluso cuando la fase acumulada excede ampliamente $2\pi$.

Metodología
Los autores proponen un marco basado en un prior de incremento temporal acotado sobre la dinámica de la muestra, operando bajo la aproximación de película delgada con iluminación espacialmente coherente.

1. Parametrización Multiplicativa Inter-frame:
   En lugar de reconstruir cada cuadro $\psi(x, y, t)$ independientemente, la evolución del campo se parametriza mediante factores de actualización multiplicativos inter-frame, $\Delta\psi(x, y, t)$. El campo en el tiempo $t + \Delta t$ se expresa como $\psi(t + \Delta t) = \psi(t) \cdot \Delta\psi(t)$. Esto transforma la reconstrucción en un problema de factibilidad para la secuencia de factores $\Delta\psi$.

2. Restricción de Sector Circular (Dominio del Objeto):
   El prior central es un límite sobre la tasa de cambio de la longitud de la trayectoria óptica (OPL) de la muestra. Este límite físico se traduce en una restricción sobre el incremento de fase inter-frame, $|\Delta\phi| \leq \alpha/2$. Combinado con una condición de pasividad (la muestra no puede amplificar el campo, por lo que $|\Delta\psi| \leq 1$), los valores admisibles para cada píxel de $\Delta\psi$ se confinan a un sector circular en el plano complejo.

   • Esta restricción generaliza la clásica restricción de no negatividad (un sector degenerado en el eje real).
   • Crucialmente, dado que cada incremento está acotado dentro de $(-\pi, \pi)$, la trayectoria de fase completa se recupera mediante la integración secuencial de estos pequeños incrementos, resolviendo automáticamente los envoltorios de fase sin algoritmos de desenrollado explícitos.

3. Restricción en el Dominio de la Medición:
   La segunda restricción impone la consistencia con la película de intensidad de difracción registrada. Un proyector en el dominio de la medición reemplaza la amplitud calculada del campo propagado con la amplitud medida, manteniendo la fase (proyección de magnitud Gerchberg–Saxton–Fienup).

4. Algoritmo de Reconstrucción:
   El problema se resuelve como un problema de factibilidad utilizando métodos de proyección alternada con reflexiones (Reflejar-Relajar-Reflejar y Reflexiones Alternadas Promediadas Relajadas) para manejar la no convexidad de la restricción de medición.

   • Etapa 1: Reconstruye el campo de sonda estático (condición inicial) a partir del primer cuadro.
   • Etapa 2: Reconstruye los factores inter-frame $\Delta\psi$ en un paso hacia adelante (integrando desde $t=0$) y un paso hacia atrás (integrando desde $t=J$). La estimación final combina los pasos hacia adelante y hacia atrás para equilibrar la acumulación de error direccional.

Contribuciones Clave

• Prior Minimalista: El artículo demuestra que un único prior motivado físicamente (incremento temporal acotado) es suficiente para reconstruir objetos de fase dinámicos sin regularizadores espaciales (por ejemplo, variación total, flujo óptico) o priores aprendidos.
• Desenrollado de Fase Automático: Al recuperar la trayectoria mediante la integración de incrementos acotados, el método maneja inherentemente las fases acumuladas que abarcan muchos ciclos de $2\pi$, eliminando la necesidad de desenrollado de fase explícito.
• Generalizabilidad: El marco se aplica sin modificaciones tanto a dinámicas sintéticas calibradas (SLM) como a procesos fotoquímicos reales y desconocidos.

Resultados Experimentales
El enfoque se validó en tres escenarios distintos:

1. Dinámicas de Reacción-Difusión (SLM): Se reconstruyó un patrón de fase que se expandía espacialmente. Los incrementos inter-frame recuperados y la trayectoria de fase integrada mostraron un buen acuerdo con la verdad fundamental, resolviendo el frente propagante de las especies químicas.
2. Dinámicas de Crecimiento Puro (SLM): Se reconstruyó un soporte espacial fijo con una fase creciendo monótonamente hasta $10\pi$. El método recuperó con éxito los valores de fase y el patrón espacial, confirmando la capacidad de manejar grandes fases acumuladas sin desenrollado. Se observaron artefactos menores en forma de anillos concéntricos debido a la acumulación de error a lo largo del tiempo.
3. Impresión 3D de Fotopolímeros in-situ: El marco se aplicó para monitorear un proceso de fotopolimerización bajo luz azul patroneada. La trayectoria de fase, que representa la evolución espaciotemporal del índice de refracción, se recuperó cuantitativamente. Los resultados rastrearon la expansión de la región polimerizada y el crecimiento de la magnitud de la fase, proporcionando un observable directo para el proceso de impresión.

Significado y Afirmaciones
Los autores afirman que la separación del acoplamiento espacial (manejado por el modelo forward de difracción de Fresnel) y el acoplamiento temporal (manejado por la parametrización inter-frame y la restricción de sector circular) permite una reconstrucción efectiva con muy pocos priores.

El artículo postula que este enfoque ofrece una herramienta práctica para la identificación de sistemas fotoquímicos. La trayectoria de fase recuperada sirve como un observable cuantitativo de la evolución espaciotemporal del índice de refracción, que puede utilizarse para ajustar parámetros cinéticos de ecuaciones de reacción-difusión que gobiernan la polimerización. Los autores señalan que, aunque la acumulación de error en secuencias largas es una limitación, el método proporciona una línea base robusta para la imagen dinámica que no requiere datos de entrenamiento complejos ni estructuras de referencia diseñadas. El trabajo sugiere que las restricciones de soporte acotadas en el plano complejo proporcionan por sí solas una regularización temporal suficiente para la CDI dinámica.