Solo el borde importa: BEM

El BEM parte de una observación profunda: para medios homogéneos por trozos, las ecuaciones de Maxwell se pueden reformular como ecuaciones integrales sobre las superficies que separan los materiales. El volumen desaparece. Solo queda el borde.

La idea

En FEM y FDTD discretizas todo el volumen — cada punto del espacio donde el campo puede existir. En BEM, usas el hecho de que dentro de cada región homogénea, el campo satisface la ecuación de Helmholtz. Y la solución de la ecuación de Helmholtz en una región homogénea se puede escribir completamente en términos de los campos en su frontera, usando el teorema integral de Green:

donde $G(\mathbf{r}, \mathbf{r}')$ es la función de Green de espacio libre — la respuesta del medio a una fuente puntual. En 3D:

Si discretizas la superficie en $N$ elementos y aplicas las condiciones de contorno (continuidad de los campos tangenciales), obtienes un sistema lineal. La clave: las incógnitas viven solo en la superficie, no en el volumen.

Volumen vs. superficie

Aquí está la diferencia visual entre discretizar el volumen (FEM/DDA) y discretizar solo el borde (BEM). Ajusta la resolución y observa cómo escala el número de incógnitas:

En 2D, el volumen tiene $O(N^2)$ incógnitas y el borde $O(N)$. En 3D, la diferencia es aún mayor: $O(N^3)$ vs $O(N^2)$. Para una esfera de 100 nm resuelta con celdas de 2 nm, un método volumétrico necesita ~65.000 incógnitas; BEM necesita ~8.000.

Acelerando BEM: el método de multipolos rápido

La matriz densa de BEM es su talón de Aquiles. Para $N = 10^4$ elementos de superficie, la matriz tiene $10^8$ entradas. Almacenarla y resolver el sistema con un método directo es prohibitivo para problemas grandes.

La solución es el FMM (Fast Multipole Method): las interacciones entre elementos lejanos se agrupan en expansiones multipolares. Solo las interacciones cercanas se calculan exactamente. El resultado: la multiplicación matriz-vector pasa de $O(N^2)$ a $O(N \log N)$ o incluso $O(N)$, sin almacenar la matriz completa. Con FMM + solucionador iterativo (GMRES), BEM escala a millones de elementos.

Ventajas y limitaciones

• Reducción dimensional: un problema 3D se reduce a un problema en superficies 2D. Menos incógnitas, menos memoria.
• Dominios abiertos gratis: la función de Green incorpora la condición de radiación al infinito automáticamente. No hay PML, no hay truncamiento artificial. La luz se va al infinito de verdad.
• Campos cercanos precisos: BEM calcula los campos directamente a partir de las corrientes superficiales. No hay interpolación de una rejilla — ideal para hotspots plasmónicos donde el campo varía en fracciones de nanómetro.
• Matrices densas: cada elemento de superficie interactúa con todos los demás (vía la función de Green). La matriz es $N \times N$ y llena. Sin FMM, almacenarla cuesta $O(N^2)$ y resolver el sistema cuesta $O(N^3)$ con métodos directos.
• Integrales singulares: cuando $\mathbf{r} \to \mathbf{r}'$, la función de Green diverge. Estas integrales singulares requieren técnicas especiales (integración analítica, regularización, sustracción de singularidades). Es la parte más difícil de implementar BEM correctamente.
• Solo medios homogéneos por trozos: cada región debe tener $\varepsilon$ constante. Gradientes de índice o materiales inhomogéneos no encajan en BEM puro. Para esos casos, se combina BEM con FEM (métodos híbridos FEM-BEM).

¿Dónde brilla BEM?

BEM es el método natural para problemas de scattering en espacio abierto con interfaces bien definidas:

• Plasmónica: nanopartículas metálicas con campos enormes en puntas y rendijas. BEM da precisión subnanométrica en los campos cercanos.
• SERS: Surface-Enhanced Raman Scattering necesita calcular el campo local con precisión extrema — el enhancement va como $|E|^4$.
• Antenas ópticas: dímeros, bow-ties, nanoantenas de forma compleja. BEM captura el acoplamiento entre estructuras cercanas con precisión.
• Secciones eficaces: BEM calcula naturalmente $Q_{\text{ext}}$, $Q_{\text{sca}}$ y $Q_{\text{abs}}$ a partir del teorema óptico, sin necesidad de integrar campos lejanos en una esfera.

BEM en la práctica

El software de referencia en nanofotónica es MNPBEM (Metal NanoParticle BEM), un toolbox de MATLAB/Octave desarrollado por Hohenester y Trügler (Universidad de Graz). Usa la formulación PMCHWT, incluye materiales dispersivos (Drude, datos tabulados de Johnson & Christy), y calcula campos cercanos, secciones eficaces, EELS y espectros SERS. Otros: SCUFF-EM (open-source, C++, FMM incluido, problemas grandes) y Bempp (Python, genérico).

Un ejemplo concreto: el campo cercano de un dímero de nanopartículas de oro con gap de 1 nm — el hotspot del que hablamos en el módulo de plasmónica. MNPBEM con ~2.000 elementos de superficie por esfera (4.000 total) resuelve el problema en minutos y da el enhancement con precisión de ~1%. Con FDTD, necesitarías una rejilla de 0.25 nm de resolución y un dominio enorme — horas de cálculo para peor precisión.

Cerrando el módulo

Con BEM cerramos el recorrido por los métodos numéricos. Cada uno tiene su nicho:

• TMM: capas planas — instantáneo.
• Rayleigh/Mie: partículas pequeñas o esféricas — analítico.
• FDTD: respuesta temporal, geometría arbitraria — intuitivo.
• DDA: partículas no esféricas en espacio abierto — flexible.
• FEM: geometrías complejas en dominio cerrado — adaptativo.
• BEM: interfaces nítidas en espacio abierto — preciso.

El arte no está en conocer un solo método, sino en saber cuál usar para cada problema. Y la mejor forma de desarrollar esa intuición es resolver problemas — muchos problemas — con herramientas diferentes y comparar.

Ejercicios