Practical Vectorial Mode Solver for Dielectric Waveguides Based on Finite Differences

Resumen Técnico: Solutor Vectorial Práctico para Guías de Onda Dieléctricas Basado en Diferencias Finitas

Autor: Ergun Simsek
Fecha: 15 de abril de 2025

1. Planteamiento del Problema

El análisis modal de guías de onda dieléctricas es fundamental en la fotónica moderna y las telecomunicaciones. Aunque existen métodos numéricos establecidos (como el Método de Elementos Finitos - FEM, y el Método de Diferencias Finitas - FDM), el tratamiento de la ecuación de onda vectorial completa presenta desafíos significativos:

• Soluciones espurias: La aparición de modos no físicos en las simulaciones.
• Ineficiencia computacional: Dificultades para resolver sistemas complejos con alta precisión.
• Problemas en interfaces: La presencia de discontinuidades en la permitividad ($\varepsilon$) en las interfaces de materiales genera errores numéricos, especialmente al calcular la componente longitudinal del campo eléctrico ($E_z$) mediante derivadas espaciales de $\varepsilon$.
• Enfoque tradicional: Muchos métodos convencionales se basan en la formulación del campo magnético o en simplificaciones que no garantizan un cumplimiento estricto de las condiciones de frontera en todas las componentes del campo eléctrico.

2. Metodología

El autor propone un solutor numérico basado en diferencias finitas diseñado específicamente para la formulación del campo eléctrico en guías de onda dieléctricas lineales ópticas y no magnéticas.

• Formulación Vectorial Completa:

  • Se parte de la ecuación de onda vectorial para el campo eléctrico $\mathbf{E}$.
  • A diferencia de enfoques que calculan $E_z$ post-procesando $E_x$ e $E_y$, este método incorpora las tres componentes del campo eléctrico ($E_x, E_y, E_z$) en una única formulación autoconsistente.
  • Esto evita la necesidad de derivar $E_z$ directamente de las derivadas de la permitividad en las interfaces, reduciendo así los artefactos numéricos.

• Discretización y Ecuaciones:

  • Se asume una variación de campo armónico en el tiempo $e^{j(\omega t - \beta z)}$, reduciendo el problema 3D a 2D.
  • Las ecuaciones acopladas se transforman en un problema de autovalores generalizado.
  • Para manejar la dependencia cuadrática de la constante de propagación $\beta$, se introduce una variable auxiliar $F = \beta E$, convirtiendo el sistema no lineal en un sistema lineal aumentado:
    $$\begin{pmatrix} M & L \\ 0 & I \end{pmatrix} \begin{pmatrix} E \\ F \end{pmatrix} = \beta \begin{pmatrix} 0 & I \\ I & 0 \end{pmatrix} \begin{pmatrix} E \\ F \end{pmatrix}$$
  • Se utiliza un esquema de diferencias finitas de cuarto orden (estencil central) para aproximar las derivadas espaciales, mejorando la precisión respecto a esquemas de primer orden.
  • La permitividad relativa ($\varepsilon_r$) se evalúa en el centro de las celdas de la malla, mientras que los campos se evalúan en los bordes.

• Implementación:

  • El problema de autovalores se resuelve utilizando la función eigs de MATLAB, optimizada para matrices dispersas.
  • Se filtran los modos físicamente significativos ($n_{eff} < n_{ring}$) y se normalizan los campos para una densidad de potencia de 1 W/m².

3. Contribuciones Clave

• Formulación Autoconsistente de Campo Eléctrico: La inclusión simultánea de las tres componentes de $\mathbf{E}$ en la matriz del sistema asegura un cumplimiento robusto de las condiciones de frontera en las discontinuidades de permitividad, mitigando errores comunes en esquemas simplificados.
• Equilibrio entre Simplicidad y Precisión: El método no busca la precisión de máquina (como lo hacen los solutores de alto orden complejos), sino ofrecer una herramienta práctica, interpretable y eficiente para análisis exploratorio y su integración con enfoques basados en datos.
• Reducción de Artefactos Numéricos: Al evitar el cálculo directo de $E_z$ a partir de derivadas de $\varepsilon$ en las interfaces, se minimizan las inexactitudes discretas típicas en los bordes de los materiales.

4. Resultados

La precisión del solutor se validó comparando sus resultados con un solutor comercial de referencia (Tidy3D) en dos escenarios distintos a una longitud de onda de 1550 nm:

• Caso 1: Guía de onda rectangular de $Si_3N_4$ (1.6 µm x 0.7 µm) en $SiO_2$.

  • Se compararon los primeros 4 modos resonantes.
  • La diferencia en los índices de refracción efectivos ($n_{eff}$) fue menor al 0.23% para todos los modos.
  • Los perfiles de campo eléctrico y magnético mostraron una distribución de intensidad concentrada en el núcleo, coincidiendo con la física esperada.

• Caso 2: Guía de onda trapezoidal sobre sustrato con película delgada (simulando niobato de litio).

  • Geometría compleja con ancho variable (1.2 µm a 1.0 µm).
  • La diferencia en $n_{eff}$ fue menor al 0.88% para los primeros 4 modos.
  • El solutor capturó correctamente la simetría del modo fundamental y las variaciones de campo en las interfaces.

En ambos casos, el tiempo de cálculo fue de aproximadamente 2-3 segundos para obtener los primeros 10 modos, demostrando alta eficiencia computacional.

5. Significado e Impacto

Este trabajo presenta una herramienta computacional valiosa para el diseño y análisis de guías de onda dieléctricas. Su principal valor radica en:

• Accesibilidad: Ofrece una implementación simplificada basada en diferencias finitas que es más fácil de entender e integrar que los métodos de elementos finitos complejos o los solutores de alto orden.
• Robustez: Proporciona soluciones modales fiables sin requerir la complejidad algorítmica extrema de los métodos de precisión de máquina, siendo suficiente para la mayoría de las aplicaciones de ingeniería fotónica.
• Versatilidad: Su estructura flexible permite futuras extensiones, como la incorporación de técnicas de promediado de permitividad o tratamientos de frontera más avanzados si se requiere mayor precisión en casos específicos.

En conclusión, el solutor propuesto llena un nicho importante entre los métodos aproximados rápidos y los solutores vectoriales completos de alta complejidad, ofreciendo un equilibrio óptimo para el análisis modal en estructuras dieléctricas modernas.