Full-Wave Harmonic Balance Framework for Dispersive Time-Varying Photonic Structures

Resumen Técnico: Marco de Equilibrio Armónico de Onda Completa para Estructuras Fotónicas Dispersivas y Variantes en el Tiempo

Planteamiento del Problema
Las estructuras fotónicas que varían en el tiempo ofrecen capacidades únicas, tales como la conversión de frecuencia, la no reciprocidad y el control dinámico del frente de onda, al romper la simetría de traslación temporal. Sin embargo, modelar con precisión estos sistemas sigue siendo un desafío significativo cuando los materiales involucrados son fuertemente dispersivos (por ejemplo, óxidos conductores transparentes, grafeno, materiales con épsilon cercano a cero). En tales plataformas, la modulación temporal no puede tratarse como una simple perturbación instantánea de un parámetro independiente de la frecuencia. En su lugar, la causalidad dicta que la respuesta del material depende de la historia de la excitación, lo que significa que cada armónico de Floquet generado experimenta una respuesta material distinta, permaneciendo al mismo tiempo acoplado a todos los demás armónicos a través de la modulación. Los métodos computacionales existentes enfrentan compromisos: los métodos en el dominio del tiempo (como FDTD) a menudo requieren simulaciones largas para alcanzar el estado estacionario y tienen dificultades con las actualizaciones causales y dispersivas, mientras que los métodos analíticos especializados o de ondas acopladas a menudo carecen de la flexibilidad geométrica para manejar estructuras 3D complejas, medios inhomogéneos o espectros de modulación arbitrarios.

Metodología
Los autores proponen un marco de equilibrio armónico de onda completa que resuelve directamente la respuesta electromagnética en estado estacionario en el dominio de la frecuencia, evitando el avance temporal. El núcleo de la metodología consiste en reformular el problema variante en el tiempo como un conjunto acoplado de subproblemas en el dominio de la frecuencia dentro de un espacio de Floquet extendido.

Los elementos clave de la formulación incluyen:

• Representación de Fuentes Secundarias: Los componentes variantes en el tiempo no se representan como condiciones de contorno dependientes del tiempo, sino como fuentes secundarias inducidas. Se utilizan densidades de polarización volumétrica ($P$) para medios volúricos modulados, y densidades de corriente superficial ($J_s$) para láminas conductoras moduladas.
• Expansión de Floquet: Todos los campos y fuentes inducidas se expanden en una serie de Floquet truncada ($F(r, t) = \sum F_n(r)e^{j\omega_n t}$), donde $\omega_n = \omega_0 + n\Omega$.
• Formulación de Sistema Acoplado: El problema se resuelve como un sistema lineal global $(L - SC)X = Y_{ext}$.
  • $L$ representa el operador de Maxwell en el dominio de la frecuencia de bloque para el fondo estático, evaluado en la frecuencia física de cada armónico ($\omega_n$) para considerar la dispersión del material.
  • $S$ representa los operadores de inyección de fuentes.
  • $C$ representa los operadores de generación de fuentes derivados de la expansión de Fourier temporal de los parámetros del material. Estos operadores introducen un acoplamiento fuera de la diagonal en el espacio armónico, representando el mecanismo de conversión de frecuencia.
• Implementación: El marco se implementa utilizando el Método de Elementos Finitos (FEM). El vector de incógnitas estándar de FEM se eleva a un vector de dominio armónico, permitiendo que cada grado de libertad espacial porte múltiples coeficientes armónicos. Esto preserva la flexibilidad geométrica del modelado de onda completa mientras captura el acoplamiento inter-armónico.

Contribuciones Clave

1. Marco Unificado: El método proporciona una estructura algebraica única para tratar tanto medios volúricos modulados como láminas conductoras moduladas, soportando espectros de modulación temporal arbitrarios (monocromáticos, multitono o no sinusoidales).
2. Manejo de la Dispersión: A diferencia de los métodos que asumen una dispersión débil o utilizan aproximaciones de envolvente, este marco evalúa la dispersión del material y los operadores de radiación en la frecuencia específica de cada armónico generado.
3. Observables Resueltos por Armónicos: La solución proporciona directamente amplitudes, fases y perfiles de campo en estado estacionario para todos los armónicos de Floquet retenidos, permitiendo el cálculo de la dispersión, absorción, flujo de potencia y distribuciones de campo específicos de cada armónico sin necesidad de post-procesamiento de datos en el dominio del tiempo.

Resultos
El marco fue validado y demostrado a través de dos ejemplos representativos:

1. Cilindro de Grafeno Enrollado con Bombeo Paramétrico:

   • Configuración: Un cilindro dieléctrico envuelto en una lámina de grafeno con un nivel de Fermi modulado temporalmente ($\Omega = 2\omega_0$), operando en un régimen de amplificación paramétrica degenerada.
   • Validación: Los resultados numéricos para la dispersión y absorción resueltas por armónicos mostraron una excelente concordancia con una solución analítica de dispersión de Floquet derivada para la geometría cilíndrica.
   • Hallazgos: La simulación capturó con éxito la redistribución de energía desde la portadora hacia el armónico de conversión descendente ($n=-1$), confirmando la capacidad del marco para modelar la fuerte dispersión y la modulación rápida en sistemas de plasmones superficiales resonantes.

2. Metasuperficie de Espacio-Tiempo Óptico (basada en ITO):

   • Configuración: Una metasuperficie reflectante que consiste en una rejilla de $TiO_2$ binaria sobre una capa de ITO (óxido de indio y estaño) con épsilon cercano a cero, respaldada por oro. La frecuencia de plasma del ITO se moduló para impulsar cambios ultrafast en la permitividad.
   • Aplicación: El marco se utilizó para diseñar y optimizar la metasuperficie para la deflexión de haz por conversión de frecuencia. El objetivo era convertir una onda incidente en bandas laterales y redirigirlas hacia un orden de difracción espacial específico ($m=1$).
   • Hallazgos: El diseño optimizado redirigió con éxito las bandas laterales generadas ($\omega_{\pm 1}$) hacia el primer orden de difracción espacial. Los resultados demostraron la interacción entre la dispersión ENZ, el acoplamiento Drude dependiente del tiempo y la difracción espacial, validando la utilidad del método para el diseño de dispositivos fotónicos complejos y realistas.

Significancia
El artículo afirma que este marco establece una plataforma general para modelar y diseñar estructuras fotónicas activas y dispersivas. Al combinar la transparencia física de los enfoques de equilibrio armónico con la flexibilidad geométrica del FEM de onda completa, el método supera las limitaciones de los resolvedores existentes para sistemas variantes en el tiempo. Permite a los investigadores cuantificar directamente cómo la modulación redistribuye la energía electromagnética entre canales desplazados en frecuencia, proporcionando observables esenciales (potencia, absorción y difracción resueltas por armónicos) para el diseño de sistemas fotónicos de espacio-tiempo de próxima generación basados en grafeno, óxidos conductores transparentes y metasuperficies resonantes. Este trabajo posiciona al marco no solo como una herramienta de análisis, sino también como una herramienta para el diseño inverso de dispositivos que requieren conversión de frecuencia, amplificación paramétrica y control dinámico del frente de onda.