¿Por qué simular?

Las ecuaciones de Maxwell son cuatro líneas. Caben en una camiseta. Describen toda la óptica, toda la radiofrecuencia, toda la luz visible, todo lo que hace un láser, un arcoíris, una fibra óptica o una nanopartícula de oro. El problema no es que no tengamos las ecuaciones — las tenemos desde 1865. El problema es que, salvo en un puñado de casos especiales, no sabemos resolverlas con lápiz y papel.

Cuando hay solución exacta

Hay geometrías con la suficiente simetría para que Maxwell admita soluciones analíticas. Capas planas infinitas: las ondas se reflejan y transmiten, y todo se reduce a multiplicar matrices 2×2 (lo veremos en el artículo 02). Esferas homogéneas: Gustav Mie encontró la solución en 1908, expandiendo los campos en armónicos esféricos. El resultado es una serie infinita que se puede sumar numéricamente con precisión arbitraria.

Mie es la joya de la corona de las soluciones analíticas en óptica. La curva de eficiencia de extinción $Q_{\text{ext}}$ muestra un comportamiento fascinante: para partículas pequeñas ( $x = 2\pi a/\lambda \ll 1$ , el parámetro de tamaño) la dispersión crece como $x^4$ — es la ley de Rayleigh, la razón por la que el cielo es azul. A medida que la partícula crece, aparecen resonancias — los modos propios de la esfera. Y para partículas muy grandes, $Q_{\text{ext}}$ converge a 2, no a 1: la famosa paradoja de extinción, donde la partícula intercepta el doble de luz de lo que su sección geométrica sugiere.

Explora la solución de Mie. Selecciona materiales con distinto índice de refracción y observa cómo cambian las resonancias. Con un metal (índice complejo), parte de la luz se absorbe dentro de la partícula:

Explorar

Material

Sin absorción — toda la luz se dispersa (Q_ext = Q_sca). Las oscilaciones son resonancias de Mie.

Fíjate en lo elaborado del resultado, incluso para la geometría más sencilla posible — una esfera. Las resonancias se estrechan con índices altos (silicio), se amortiguan con absorción (metal), y el límite asintótico $Q_{\text{ext}} \to 2$ aparece siempre. Todo esto sale de una solución exacta. Bonito, ¿verdad?

Ejercicio 1

Selecciona «Vidrio» (n = 1.50) y observa la primera resonancia de Mie. ¿A qué valor de $x$ aparece? Ahora cambia a «Silicio» (n = 3.50). ¿Se mueve la primera resonancia? ¿En qué dirección y por qué?

Solución

Para vidrio, la primera resonancia aparece alrededor de x ≈ 4. Para silicio, baja a x ≈ 1.5. La resonancia ocurre cuando la partícula «cabe» aproximadamente una longitud de onda interna — y como la longitud de onda interna es

\lambda / n

, un índice más alto necesita un tamaño físico menor (menor x) para resonar.

Ejercicio 2

Selecciona «Metal» y observa que $Q_{\text{ext}}$ es mayor que $Q_{\text{sca}}$ . La diferencia es $Q_{\text{abs}}$ . ¿A qué valor de $x$ la absorción domina sobre el scattering? ¿Y cuándo el scattering domina?

Solución

Para partículas metálicas pequeñas (x < 1), la absorción domina:

Q_{\text{abs}} \propto x

mientras que

Q_{\text{sca}} \propto x^4

. Para x grande, el scattering crece más rápido y acaba dominando. El cruce ocurre alrededor de x ≈ 1–2 para este metal (m = 0.2 + 3i).

Ahora la mala noticia.

Cuando no la hay

La solución de Mie funciona porque la esfera tiene simetría esférica perfecta. Los armónicos esféricos forman una base natural, las condiciones de contorno se separan limpiamente, y todo encaja. ¿Pero qué pasa si la partícula no es una esfera? ¿Si es un cubo, una estrella, un dímero de dos esferas, un agregado fractal? ¿Si está cerca de una superficie, dentro de un sustrato, o rodeada de otras partículas?

La simetría se rompe y la solución analítica desaparece. Las ecuaciones de Maxwell siguen siendo válidas — la física no cambia — pero ya no sabemos resolverlas a mano. Necesitamos un ordenador.

Y no es un problema menor. La inmensa mayoría de los problemas reales en fotónica caen en esta categoría:

Cristales fotónicos: periodicidad 2D o 3D, bandas prohibidas, guiado de la luz por defectos en la estructura.
Plasmónica: nanopartículas metálicas con formas arbitrarias, campos enormes en puntas y rendijas, acoplamiento entre partículas.
Fibras y guías de onda: sección transversal no circular, materiales con gradiente de índice, curvatura, empalmes.
Dispositivos integrados: divisores, acopladores, moduladores — geometrías diseñadas para una función específica, sin simetría natural.
Metamateriales: estructuras subwavelength diseñadas para propiedades electromagnéticas imposibles en la naturaleza.

Para todos estos problemas, la estrategia es la misma: discretizar las ecuaciones de Maxwell de alguna forma inteligente y dejar que el ordenador las resuelva. Las formas «inteligentes» de hacerlo constituyen un campo entero de la física computacional — y ese campo es el tema de este módulo.

El paisaje de métodos

No existe un método numérico universal. Cada uno discretiza algo distinto (el volumen, la superficie, las capas), resuelve las ecuaciones de una forma distinta (en el dominio temporal o frecuencial), y tiene sus propias ventajas y limitaciones. Elegir el método correcto para un problema concreto es una de las decisiones más importantes en simulación numérica.

Aquí tienes el mapa de los seis métodos que vamos a estudiar en este módulo. Pulsa en cada uno para ver qué lo hace especial:

Explorar

Pulsa en un método para ver sus detalles

Fíjate en una distinción fundamental: los métodos de volumen (FDTD, FEM, DDA) discretizan todo el espacio — meten puntos o elementos dentro y fuera del objeto. Los métodos de superficie (BEM) solo discretizan el borde — la interfaz entre materiales. Y el Transfer Matrix opera en una dimensión, capa a capa. Cada enfoque tiene su precio y su recompensa.

¿Cómo se elige?

La elección del método depende de tres cosas:

La geometría del problema. ¿Es 1D (capas)? Usa Transfer Matrix. ¿Es una partícula compacta en espacio abierto? BEM o DDA. ¿Es una estructura periódica o una guía de onda? FDTD o FEM.
Lo que quieres calcular. ¿Un espectro de reflexión? Transfer Matrix es instantáneo. ¿La respuesta a un pulso ultracorto? FDTD te da la evolución temporal directamente. ¿Campos cercanos con resolución subnanométrica? BEM maneja las singularidades mejor.
Los recursos disponibles. Algunos métodos escalan brutalmente con el tamaño del problema. FDTD en 3D con buena resolución puede necesitar terabytes de RAM. BEM genera matrices llenas que son caras de invertir. Transfer Matrix corre en milisegundos en un portátil.

No hay un «mejor método». Hay el método adecuado para cada problema. A lo largo de los próximos seis artículos, vamos a construir cada uno desde cero — empezando por el más sencillo (Transfer Matrix) y terminando por el más denso conceptualmente (BEM). Al final del módulo, tendrás la intuición para saber cuál usar y por qué.