Respuesta EM · Artículo 06

Teorema óptico

La extinción — cuánta energía una partícula saca del haz incidente — es exactamente proporcional a la amplitud de dispersión en la dirección forward. Es el resultado más profundo de la teoría de dispersión.

Imagina una onda plana que llega a una nanopartícula. Parte de la energía se absorbe (se convierte en calor). Parte se dispersa (se redirige en todas direcciones). La suma de ambas es la extinción: σext=σabs+σscat\sigma_{\text{ext}} = \sigma_{\text{abs}} + \sigma_{\text{scat}}. Medir la extinción es sencillo — mides cuánta luz falta en la dirección forward. Pero calcularla requiere sumar la absorción (que depende del calentamiento interno) y la dispersión (que requiere integrar sobre todas las direcciones). El teorema óptico dice que hay un atajo: basta con conocer la amplitud de dispersión en una sola dirección — la dirección forward.

Partícula pequeña: σ_ext desde Im(α)

Empecemos por el caso más simple. Una partícula mucho menor que λ\lambda responde como un dipolo con polarizabilidad α(ω)\alpha(\omega). El momento dipolar inducido es p=ε0αE0\mathbf{p} = \varepsilon_0\alpha\mathbf{E}_0.

La potencia extraída del haz incidente (extinción) es el trabajo del campo incidente sobre el dipolo:

Pext=ω2Im ⁣[pE0]=ωε02E02Im[α]P_{\text{ext}} = \frac{\omega}{2}\,\text{Im}\!\left[\mathbf{p}^*\cdot\mathbf{E}_0\right] = \frac{\omega\varepsilon_0}{2}|E_0|^2\,\text{Im}[\alpha]

Dividiendo por la intensidad incidente I=ε0cE02/2I = \varepsilon_0 c|E_0|^2/2:

σext=PextI=ωcIm[α(ω)]=kIm[α]\sigma_{\text{ext}} = \frac{P_{\text{ext}}}{I} = \frac{\omega}{c}\,\text{Im}\left[\alpha(\omega)\right] = k\,\text{Im}[\alpha]

O equivalentemente:

σext=4πωcIm ⁣[α(ω)]\sigma_{\text{ext}} = \frac{4\pi\omega}{c}\,\text{Im}\!\left[\alpha(\omega)\right]

(El factor 4π4\pi depende del sistema de unidades; aquí usamos Gaussianas para la polarizabilidad). La extinción depende solo de Im[α]\text{Im}[\alpha] — la parte de la respuesta que está en fase con la velocidad del campo, no con el campo mismo.

El teorema óptico general

Para un dispersor de tamaño arbitrario, la onda dispersada en el campo lejano se escribe:

Escatreikrrf(r^)\mathbf{E}_{\text{scat}} \xrightarrow{r\to\infty} \frac{e^{ikr}}{r}\,\mathbf{f}(\hat{r})

donde f(r^)\mathbf{f}(\hat{r}) es la amplitud de dispersión en la dirección r^\hat{r}. El teorema óptico establece:

σext=4πkIm ⁣[e^f(k^i)]\sigma_{\text{ext}} = \frac{4\pi}{k}\,\text{Im}\!\left[\hat{e}^*\cdot\mathbf{f}(\hat{k}_i)\right]

donde k^i\hat{k}_i es la dirección de la onda incidente (forward) y e^\hat{e} es su polarización. La extinción total — absorción más dispersión a todos los ángulos — está codificada en la amplitud de dispersión en una sola dirección.

Derivación por conservación de energía

El campo total es E=Ei+Escat\mathbf{E} = \mathbf{E}_i + \mathbf{E}_{\text{scat}}, donde Ei=E0e^eikr\mathbf{E}_i = E_0\hat{e}\,e^{i\mathbf{k}\cdot\mathbf{r}} es la onda incidente. El flujo de energía a través de una esfera lejana de radio RR es:

SdA=Si+Sscat+SintdA\oint \langle\mathbf{S}\rangle\cdot d\mathbf{A} = \oint \langle\mathbf{S}_i + \mathbf{S}_{\text{scat}} + \mathbf{S}_{\text{int}}\rangle\cdot d\mathbf{A}

El término Si\mathbf{S}_i (Poynting de la onda incidente) se anula al integrar sobre la esfera cerrada (tanta energía entra como sale). El término Sscat\mathbf{S}_{\text{scat}} da PscatP_{\text{scat}}. El término de interferencia Sint\mathbf{S}_{\text{int}} da Pext-P_{\text{ext}}.

La evaluación del término de interferencia requiere calcular:

SintRe ⁣[Ei×Hscat+Escat×Hi]\mathbf{S}_{\text{int}} \propto \text{Re}\!\left[\mathbf{E}_i \times \mathbf{H}_{\text{scat}}^* + \mathbf{E}_{\text{scat}} \times \mathbf{H}_i^*\right]

Al integrar sobre la esfera lejana, las fases ei(krkr)e^{i(\mathbf{k}\cdot\mathbf{r} - kr)} oscilan rápidamente y se cancelan en todas las direcciones excepto la forward (r^=k^i\hat{r} = \hat{k}_i), donde la fase es estacionaria. La evaluación por fase estacionaria da exactamente el teorema óptico.

¿Por qué la dirección forward?

El teorema óptico no es evidente — ¿por qué la extinción total depende solo de la dispersión forward? La razón es la interferencia. Lo que mide un detector en la dirección forward es la suma del campo incidente y el campo dispersado:

Eforward=Ei+Escat(k^i)\mathbf{E}_{\text{forward}} = \mathbf{E}_i + \mathbf{E}_{\text{scat}}(\hat{k}_i)

La intensidad detectada es Ei+Escat2|\mathbf{E}_i + \mathbf{E}_{\text{scat}}|^2. La reducción respecto a Ei2|\mathbf{E}_i|^2 — la «sombra» — viene del término de interferencia, que es proporcional a Re[EiEscat]\text{Re}[\mathbf{E}_i^*\cdot\mathbf{E}_{\text{scat}}]. Pero la extinción (potencia total extraída) es proporcional a Im[f]\text{Im}[\mathbf{f}], no a Re.

¿Cómo se reconcilia? La onda dispersada forward tiene un desfase de π/2\pi/2 respecto a la incidente (viene de la condición de radiación eikr/re^{ikr}/r). Así, Re[EiEscat]Im[f]\text{Re}[E_i^* \cdot E_{\text{scat}}] \propto \text{Im}[\mathbf{f}]. La sombra está hecha de Im(f).

Descomposición: absorción y dispersión

La extinción se descompone en dos canales:

σext=σabs+σscat\sigma_{\text{ext}} = \sigma_{\text{abs}} + \sigma_{\text{scat}}

Para una partícula dipolar:

El albedo a=σscat/σexta = \sigma_{\text{scat}}/\sigma_{\text{ext}} mide cuánto de la extinción es dispersión. Para partículas muy pequeñas (ka1ka \ll 1), σscat(ka)4\sigma_{\text{scat}} \propto (ka)^4 mientras σabs(ka)\sigma_{\text{abs}} \propto (ka), así que a0a \approx 0 — la absorción domina. Para partículas grandes, la dispersión domina.

Límite de la sección eficaz: 3λ²/(2π)

¿Cuál es la máxima extinción que un dispersor puede tener? Para un dispersor dipolar, la respuesta tiene una forma sorprendentemente simple. En resonancia, la polarizabilidad satisface una cota impuesta por la conservación de energía (unitariedad):

Im[α]k36πα2\text{Im}[\alpha] \geq \frac{k^3}{6\pi}|\alpha|^2

La igualdad se da cuando no hay absorción (dispersión pura). El máximo de σext\sigma_{\text{ext}} para un dipolo se alcanza cuando esta cota se satura:

σextmaˊx=3λ22π\sigma_{\text{ext}}^{\text{máx}} = \frac{3\lambda^2}{2\pi}

Este resultado es notable:

Explorar
Amortiguamiento γ (eV) 0.15
Secciones eficaces: extinción = dispersión + absorción
Sin pérdidas (γ → 0): σ_ext = σ_sca → 3λ²/(2π). Con pérdidas: la absorción aparece y la dispersión cae.

Consecuencias para el diseño de absorbedores

El límite 3λ2/(2π)3\lambda^2/(2\pi) tiene implicaciones directas para el diseño de absorbedores perfectos:

  1. No puedes hacer un absorbedor más pequeño que λ\sim\lambda: la sección eficaz máxima es λ2\sim\lambda^2, independientemente de cuánto esfuerzo pongas en la nanoestructura. La difracción pone un piso.
  2. Para absorción pura (a = 0): la máxima absorción dipolar es σabsmaˊx=3λ2/(8π)\sigma_{\text{abs}}^{\text{máx}} = 3\lambda^2/(8\pi), que se alcanza cuando σscat=σabs\sigma_{\text{scat}} = \sigma_{\text{abs}} (condición de matched scattering). Un absorbedor perfecto no puede eliminar la dispersión — necesita dispersar exactamente lo mismo que absorbe.
  3. Conexión con láseres (Módulo 05): la sección eficaz de absorción de un átomo de dos niveles en resonancia es exactamente 3λ2/(2π)3\lambda^2/(2\pi) — el límite dipolar. No es coincidencia: un átomo en resonancia es el dispersor dipolar más eficiente posible.

Conexión con Mie

El espectro de extinción de Mie (Módulo 03, Art. 01) se puede escribir como:

σext=2πk2l=1(2l+1)Re(al+bl)\sigma_{\text{ext}} = \frac{2\pi}{k^2}\sum_{l=1}^{\infty}(2l+1)\,\text{Re}(a_l + b_l)

donde ala_l y blb_l son los coeficientes de Mie. El máximo de cada término es (2l+1)λ2/(2π)(2l+1)\lambda^2/(2\pi) — cuando Re(al)=1\text{Re}(a_l) = 1 e Im(al)=0\text{Im}(a_l) = 0 (resonancia unitaria). Esto es exactamente el teorema óptico término a término.

En la práctica, para una nanopartícula de Au de 20 nm, solo el término dipolar l=1l = 1 importa (ka0.25ka \approx 0.25). La eficiencia de extinción en resonancia (~520 nm) es Qext5Q_{\text{ext}} \approx 5 — la partícula extingue un área 5 veces mayor que su sección geométrica. Pero está lejos del límite dipolar (Qextmaˊx102Q_{\text{ext}}^{\text{máx}} \sim 10^2) porque las pérdidas del oro (Im(ε) alto) reducen la amplitud de dispersión.

Cota unitaria para la polarizabilidad

La conservación de energía exige σabs=σextσscat0\sigma_{\text{abs}} = \sigma_{\text{ext}} - \sigma_{\text{scat}} \geq 0. Sustituyendo las expresiones dipolares:

kIm[α]k46πα20k\,\text{Im}[\alpha] - \frac{k^4}{6\pi}|\alpha|^2 \geq 0
Im[α]k36πα2\text{Im}[\alpha] \geq \frac{k^3}{6\pi}|\alpha|^2

Esta es la cota unitaria. Para saturarla, escribimos α=6πi/(k3(S1))\alpha = 6\pi i/(k^3(S - 1)) donde S1|S| \leq 1 es el elemento de la matriz S. La extinción es máxima cuando S=1S = -1:

σextmaˊx=kIm ⁣[6πik3(2)]=3πk2=3λ22π\sigma_{\text{ext}}^{\text{máx}} = k\,\text{Im}\!\left[\frac{6\pi i}{k^3(-2)}\right] = \frac{3\pi}{k^2} = \frac{3\lambda^2}{2\pi}

Ejercicios

Ejercicio 1

Una nanopartícula de oro de 10 nm de radio tiene una polarizabilidad en resonancia (~520 nm) de aproximadamente α(3+12i)×104\alpha \approx (3 + 12i) \times 10^4 nm³. Calcula σext\sigma_{\text{ext}}, σscat\sigma_{\text{scat}} y σabs\sigma_{\text{abs}}. ¿Cuál es el albedo? ¿Domina la absorción o la dispersión?

Solución
k=2π/520=0.0121k = 2\pi/520 = 0.0121 nm⁻¹. σext=kIm[α]=0.0121×1.2×105=1450\sigma_{\text{ext}} = k\,\text{Im}[\alpha] = 0.0121 \times 1.2\times10^5 = 1\,450 nm². α2=(32+122)×108=1.53×1010|\alpha|^2 = (3^2 + 12^2)\times10^8 = 1.53\times10^{10} nm⁶. σscat=k4α2/(6π)=(0.0121)4×1.53×1010/(6π)17\sigma_{\text{scat}} = k^4|\alpha|^2/(6\pi) = (0.0121)^4 \times 1.53\times10^{10}/(6\pi) \approx 17 nm². σabs=145017=1433\sigma_{\text{abs}} = 1\,450 - 17 = 1\,433 nm². Albedo a=17/14500.012a = 17/1\,450 \approx 0.012. La absorción domina completamente (~99%). Esto es típico de nanopartículas metálicas pequeñas: la dispersión escala como (ka)4(ka)^4 y es negligible cuando ka1ka \ll 1.
Ejercicio 2

El límite de la sección eficaz dipolar es 3λ2/(2π)3\lambda^2/(2\pi). Para λ=1550\lambda = 1550 nm (telecomunicaciones), ¿cuánto vale en μm²? Si una antena óptica de InP tiene un área física de 1 μm², ¿puede capturar más luz que su área geométrica?

Solución
σextmaˊx=3(1.55)2/(2π)=1.15\sigma_{\text{ext}}^{\text{máx}} = 3(1.55)^2/(2\pi) = 1.15 μm². Sí, una antena óptica dipolar resonante a 1550 nm puede capturar ~1.15 μm² de flujo — más que su área física de 1 μm². La antena «concentra» el campo de un área mayor que ella misma. Esto es posible porque el campo cercano distorsiona las líneas de flujo de energía, canalizándolas hacia la partícula. Es la base de las antenas ópticas para fotodetección.
Ejercicio 3

Un absorbedor perfecto dipolar tiene σabs=σscat\sigma_{\text{abs}} = \sigma_{\text{scat}} (condición de matched scattering). Demuestra que en este caso σabs=3λ2/(8π)\sigma_{\text{abs}} = 3\lambda^2/(8\pi). ¿Es mejor o peor que una partícula que solo absorbe (σscat=0\sigma_{\text{scat}} = 0)?

Solución
Si σabs=σscat\sigma_{\text{abs}} = \sigma_{\text{scat}}, entonces σext=2σabs\sigma_{\text{ext}} = 2\sigma_{\text{abs}}. Además, σext3λ2/(2π)\sigma_{\text{ext}} \leq 3\lambda^2/(2\pi), con el máximo cuando la cota unitaria se satura con S=0|S| = 0 (que da exactamente σabs=σscat\sigma_{\text{abs}} = \sigma_{\text{scat}}): σabs=3λ2/(8π)\sigma_{\text{abs}} = 3\lambda^2/(8\pi). Una partícula que «solo absorbe» (sin dispersión) tendría σabs=σext3λ2/(2π)\sigma_{\text{abs}} = \sigma_{\text{ext}} \leq 3\lambda^2/(2\pi), pero eso viola la cota unitaria salvo que α=0|\alpha| = 0. En realidad, maximizar la absorción requiere σscat=σabs\sigma_{\text{scat}} = \sigma_{\text{abs}} — dispersar tanto como absorbes. Es contraintuitivo pero es consecuencia directa de la unitariedad.