Optical Fourier Architecture for Universal Nonlinear Functions

Resumen Técnico: Arquitectura de Fourier Óptica para Funciones No Lineales Universales

Planteamiento del Problema
El artículo aborda un cuello de botella fundamental en la computación óptica: la dificultad de integrar transformaciones no lineales en procesos dinámicos estrictamente lineales sin depender de conversiones entre los domios óptico y electrónico, las cuales consumen muchos recursos y generan latencia. Mientras que los procesadores ópticos coherentes sobresalen en operaciones de matriz-vector lineales con un paralelismo masivo y latencias de sub-nanosegundo, carecen de un método determinista y universal para sintetizar funciones no lineales arbitrarias. Los enfoques previos para compilar circuitos ópticos parametrizados para la síntesis de funciones han dependido de la optimización numérica no lineal o búsquedas heurísticas, que sufren de mínimos locales, cuellos de botella de convergencia y una falta de garantía de existencia para las configuraciones de circuito requeridas.

Metodología
Los autores proponen una arquitectura algebraica exacta que evalúa una serie de Fourier arbitraria finita utilizando un circuito óptico lineal de dos modos ($2 \times 2$). La innovación central radica en tratar la variable de entrada no como un estado a ser codificado, sino como un parámetro codificado directamente en la fase de desplazadores de fase de un solo modo sintonizables.

La metodología procede a través de los siguientes pasos:

1. Representación de Fourier: Una función objetivo $f(x)$ definida en un intervalo compacto se aproxima mediante una serie de Fourier truncada $\hat{f}(x)$. Esta se mapea a un polinomio de Laurent $\hat{f}(z)$ en el círculo unitario complejo, donde $z = e^{i\pi x/\tau}$.
2. Restricción Unitaria: Dado que las redes ópticas lineales son estrictamente unitarias, el polinomio objetivo $\hat{f}(z)$ no puede realizarse solo; debe ir acompañado de un polinomio auxiliar $g(z)$ tal que la intensidad total $|\hat{f}(z)|^2 + |g(z)|^2$ permanezca constante para todo $z$ en el círculo unitario.
3. Factorización Espectral: Los autores demuestran la existencia de tal polinomio complementario $g(z)$ utilizando el teorema de Fejér-Riesz. Los coeficientes de $g(z)$ se determinan de manera determinista mediante la factorización espectral, implementada específicamente a través de la descomposición de Cholesky de una matriz de Toeplitz. Esto evita por completo la optimización iterativa.
4. Síntesis de Circuito Recursiva: Una vez establecidos los vectores de coeficientes para $\hat{f}(z)$ y $g(z)$, los autores utilizan un algoritmo de reducción algebraica (análogo al algoritmo de Schur clásico). Este algoritmo identifica recursivamente vectores de extremo ortogonales para determinar una secuencia de unitarias $SU(2)$ estáticas (realizadas como interferómetros Mach-Zehnder) y desplazadores de fase. Este proceso descompone sistemáticamente la operación objetivo en una cascada de $2N$ unitarias estáticas y desplazadores de fase para una expansión de $N$ armónicos.

Contribuciones Clave

• Prueba de Existencia: El artículo proporciona una prueba rigurosa de que existe un circuito óptico exacto y sin pérdidas para cada serie de Fourier finita, eliminando la incertidumbre inherente a los métodos de búsqueda heurística.
• Compilación Determinista: Los autores derivan un método analítico para construir el circuito. La síntesis no es iterativa, requiriendo solo una única pasada de factorización espectral y reducción recursiva, con una complejidad computacional que escala como $O(N)$.
• No Linealidad Pasiva: La arquitectura demuestra que se puede lograr una no linealidad efectiva en un sistema óptico estrictamente lineal codificando el argumento de la función en los desplazadores de fase, integrando así la computación no lineal dentro del proceso lineal que actúa sobre los modos ópticos espaciales.
• Escalabilidad: El sistema óptico resultante exhibe una profundidad de $O(N)$ para una expansión de $N$ armónicos, permitiendo que las evaluaciones de funciones ocurran dentro del tiempo de vuelo óptico pasivo.

Resultados
Los autores validaron sus afirmaciones mediante simulaciones numéricas:

• Funciones Suaves: La arquitectura sintetizó con éxito la función $\text{sinc}(x)$ suave con un error acotado utilizando una profundidad de circuito modesta. Se demostró que el error surge únicamente de la truncación de la serie de Fourier, no de la compilación del circuito.
• Funciones Discontinuas: El método fue probado con una función $\text{sinc}(x)$ escalonada por tramos con discontinuidades de salto discretas. A pesar de los desafíos planteados por las oscilaciones de Gibbs y la necesidad de altos cortes armónicos, el algoritmo se mantuvo numéricamente estable.
• Eficiencia: Para series de Fourier y polinomios auxiliares que constan de más de mil términos, el programa completo (incluyendo configuración, factorización, síntesis y graficación) se completó en menos de cinco segundos en una CPU de una computadora portátil estándar.

Significancia y Reclamaciones
El artículo afirma proporcionar una "base universal y determinista" para la computación óptica no lineal de una sola variable en plataformas fotónicas integradas. Al tender un puente entre la teoría de aproximación de funciones universales y el hardware fotónico, este trabajo ofrece una vía para operaciones no lineales programables sin la necesidad de materiales ópticos no lineales desafiantes, fuentes de un solo fotón o no linealidades inducidas por medición.

Los autores enfatizan que, si bien el algoritmo se inspira en el Procesamiento de Señales Cuánticas (QSP) y en los marcos de trabajo de re-carga de datos (data re-uploading), la computación opera enteramente sobre amplitudes ópticas coherentes clásicas. Esto hace que el enfoque sea accesible para la fotónica integrada de corto plazo. El alcance del trabajo es modesto, señalando explícamente que las estrategias algebraicas actuales están restringidas a polinomios univariados. La extensión a funciones multivariables o la concatenación de circuitos ópticos sin conversión electro-óptica sigue siendo un problema abierto.