Nanofotónica · Artículo 04

Eigenmodos electrostáticos

La respuesta plasmónica de cualquier nanopartícula se descompone en modos universales que dependen solo de la geometría — no del material ni de la frecuencia.

En M03-01 vimos que una esfera metálica resuena cuando Re(ε)=−2εm\text{Re}(\varepsilon) = -2\varepsilon_m. El «-2» viene de la geometría esférica. Si la partícula fuera un elipsoide, un cubo o una estrella, ese número sería otro. ¿Hay una forma sistemática de encontrar todas las resonancias de cualquier forma? Sí — y la respuesta es un problema de eigenmodos que separa completamente la geometría del material.

Límite electrostático

Cuando la partícula es mucho menor que la longitud de onda (a≪λa \ll \lambda), el campo electromagnético varía lentamente a la escala de la partícula. Podemos ignorar el retardo — la velocidad finita de la luz — y resolver la electrostática: ∇2ϕ=0\nabla^2\phi = 0 dentro y fuera de la partícula, con condiciones de contorno en la superficie.

En este límite, el campo eléctrico total dentro de la partícula se puede escribir como:

E~=E~ext+η M⋅E~\tilde{\mathbf{E}} = \tilde{\mathbf{E}}^{\text{ext}} + \eta\,\mathbf{M}\cdot\tilde{\mathbf{E}}

donde E~ext\tilde{\mathbf{E}}^{\text{ext}} es el campo externo, M\mathbf{M} es un operador integral que depende solo de la geometría de la partícula, y:

η=ε−εmε+εm\eta = \frac{\varepsilon - \varepsilon_m}{\varepsilon + \varepsilon_m}

es el contraste dieléctrico — un número adimensional que codifica toda la información sobre el material y la frecuencia. La separación es exacta: la geometría vive en M\mathbf{M}, el material vive en η\eta.

El operador geométrico M

El operador M\mathbf{M} es un operador integral sobre la superficie de la partícula. Su acción sobre un campo E~\tilde{\mathbf{E}} es calcular las cargas superficiales que ese campo induce, y luego el campo que esas cargas producen. Es la «respuesta electrostática» de la geometría.

Propiedades cruciales:

Problema de eigenvalores

Los modos propios satisfacen:

M⋅en=1ηn en\mathbf{M}\cdot\mathbf{e}_n = \frac{1}{\eta_n}\,\mathbf{e}_n

Cada eigenmodo en\mathbf{e}_n es una distribución de campo (y su correspondiente distribución de carga superficial) que es «natural» para esa geometría. Los eigenvalores ηn\eta_n son números reales entre -1 y +1 que dependen solo de la forma.

La respuesta del campo total se expande en estos modos:

E~=∑ncn1−η/ηn en\tilde{\mathbf{E}} = \sum_n \frac{c_n}{1 - \eta/\eta_n}\,\mathbf{e}_n

donde cn=⟨en∣E~ext⟩c_n = \langle\mathbf{e}_n|\tilde{\mathbf{E}}^{\text{ext}}\rangle es el coeficiente de acoplamiento del campo externo con el modo nn.

Condición de resonancia

La respuesta diverge cuando el denominador se anula:

η(ω)=ηn\eta(\omega) = \eta_n

Es decir, cuando el contraste dieléctrico del material, a la frecuencia ω\omega, coincide con uno de los eigenvalores geométricos. Esto es extraordinario: los modos (las «frecuencias naturales» de vibración) están determinados por la geometría, y el material solo dice a qué frecuencia se excita cada modo.

Cambiar el material (oro por plata, por ejemplo) desplaza las frecuencias de resonancia, pero no cambia los patrones de campo de los modos. Un nanotriángulo de oro y uno de plata tienen los mismos eigenmodos — solo resuenan a frecuencias distintas.

La esfera: eigenmodos analíticos

Para una esfera, los eigenmodos son los armónicos esféricos Ylm(θ,ϕ)Y_l^m(\theta,\phi), y los eigenvalores son:

ηl=12l+1\eta_l = \frac{1}{2l+1}

La condición de resonancia η(ω)=ηl\eta(\omega) = \eta_l da:

ε−εmε+εm=12l+1\frac{\varepsilon - \varepsilon_m}{\varepsilon + \varepsilon_m} = \frac{1}{2l+1}

Para l=1l=1 (dipolar): η1=1/3\eta_1 = 1/3, que da ε=−2εm\varepsilon = -2\varepsilon_m. Ahí está el famoso «-2» de M03-01 — es simplemente el eigenvalor del modo dipolar de la esfera.

Para l=2l=2 (cuadrupolar): η2=1/5\eta_2 = 1/5, que da ε=−1.5εm\varepsilon = -1.5\varepsilon_m.

Para l→∞l \to \infty: ηl→0\eta_l \to 0, que da ε→−εm\varepsilon \to -\varepsilon_m — la frecuencia del plasmón de superficie plana. Los modos de alto orden «ven» una superficie localmente plana.

Eigenmodos de la esfera

Universalidad geométrica

La belleza de este formalismo es su universalidad. Para cualquier forma de partícula:

  1. Diagonaliza M\mathbf{M} → obtienes los eigenmodos y eigenvalores.
  2. Calcula η(ω)\eta(\omega) para tu material → obtienes las frecuencias de resonancia.
  3. Proyecta el campo externo sobre los eigenmodos → obtienes la respuesta completa.

Para formas que no tienen solución analítica (nanotriángulos, nanostars, formas irregulares), M\mathbf{M} se discretiza numéricamente usando BEM. La superficie se divide en NN paneles triangulares, M\mathbf{M} se convierte en una matriz N×NN \times N, y los eigenmodos se obtienen diagonalizando esa matriz.

Resultados típicos para formas comunes:

Reglas de suma para los eigenvalores

Los eigenvalores satisfacen reglas de suma que restringen los valores posibles:

∑nηn=1(traza de M)\sum_n \eta_n = 1 \qquad \text{(traza de M)}

Esto significa que los eigenvalores no son independientes — si un modo tiene un ηn\eta_n grande (resonancia alejada del plasmón de superficie), otros deben compensar. Para la esfera: ∑l(2l+1)12l+1=∑l1\sum_l (2l+1)\frac{1}{2l+1} = \sum_l 1, que diverge — pero cada modo contribuye exactamente 1 a la traza, como debe ser.

Otra regla útil: todos los eigenvalores satisfacen 0<ηn≤10 < \eta_n \leq 1 para partículas convexas. Para partículas con concavidades (dimers, nanoshells), pueden aparecer eigenvalores negativos.

Más allá del electrostático

El formalismo de eigenmodos es exacto en el límite electrostático. Para partículas más grandes (a≳λ/10a \gtrsim \lambda/10), el retardo introduce correcciones:

Pero la imagen conceptual sigue siendo válida: la respuesta plasmónica es una superposición de modos geométricos, cada uno excitado a una frecuencia determinada por el material. El límite electrostático da la intuición correcta; las correcciones retardadas la refinan cuantitativamente.

Ejercicios

Ejercicio 1

Para una esfera de plata en aire (εm=1\varepsilon_m = 1), usa el modelo de Drude (ε(ω)=1−ωp2/ω2\varepsilon(\omega) = 1 - \omega_p^2/\omega^2, con ωp=9.0\omega_p = 9.0 eV) para calcular las frecuencias de resonancia de los modos l=1,2,3l = 1, 2, 3. Verifica que convergen hacia ωp/2\omega_p/\sqrt{2} (plasmón de superficie) cuando l→∞l \to \infty.

Solución
La condición de resonancia es ε=−(l+1)εm/l=−(l+1)/l\varepsilon = -(l+1)\varepsilon_m/l = -(l+1)/l. Con Drude: 1−ωp2/ω2=−(l+1)/l1 - \omega_p^2/\omega^2 = -(l+1)/l, es decir ωl=ωp/1+(l+1)/l=ωp/(2l+1)/l\omega_l = \omega_p/\sqrt{1 + (l+1)/l} = \omega_p/\sqrt{(2l+1)/l}. l=1l=1: ω1=9.0/3=5.20\omega_1 = 9.0/\sqrt{3} = 5.20 eV. l=2l=2: ω2=9.0/5/2=5.69\omega_2 = 9.0/\sqrt{5/2} = 5.69 eV. l=3l=3: ω3=9.0/7/3=5.89\omega_3 = 9.0/\sqrt{7/3} = 5.89 eV. Límite l→∞l \to \infty: ω∞=9.0/2=6.36\omega_\infty = 9.0/\sqrt{2} = 6.36 eV. Los modos convergen al plasmón de superficie desde abajo.
Ejercicio 2

Un nanocubo de oro tiene un eigenmodo de esquina con eigenvalor ηcorner=0.85\eta_{\text{corner}} = 0.85 y un eigenmodo de cara con ηface=0.25\eta_{\text{face}} = 0.25. ¿Cuál de los dos resuena a una frecuencia más baja (más al rojo)? ¿Cuál produce mayor enhancement de campo? Justifica ambas respuestas a partir del valor de η\eta.

Solución
Un η\eta más grande requiere un contraste dieléctrico mayor, es decir ∣ε∣|\varepsilon| más grande, lo que ocurre a frecuencias más bajas (más al rojo) en el modelo de Drude. El modo de esquina (η=0.85\eta = 0.85) resuena a frecuencia más baja que el de cara (η=0.25\eta = 0.25). El modo de esquina también produce mayor enhancement porque su campo se concentra en un volumen más pequeño (punta de la esquina) y resuena donde ∣ε∣|\varepsilon| es grande, dando una respuesta más fuerte. En general, los modos con η\eta cercano a 1 son los más «plasmónicos».
Ejercicio 3

Explica por qué el formalismo de eigenmodos predice que un nanotriángulo de plata y uno de oro tienen exactamente los mismos patrones de campo en los modos, pero a frecuencias distintas. ¿Qué cambia y qué no cambia al sustituir el material?

Solución
El operador M\mathbf{M} depende solo de la geometría. Los eigenmodos en\mathbf{e}_n y los eigenvalores ηn\eta_n son idénticos para un triángulo de oro y uno de plata. Lo que cambia es la función η(ω)\eta(\omega), que depende del material. La condición de resonancia η(ω)=ηn\eta(\omega) = \eta_n se satisface a frecuencias distintas para Au y Ag. Los patrones de campo (dónde se concentra el campo, cuántos nodos tiene) son universales — la geometría los fija. El material solo elige a qué frecuencia cada patrón se enciende.