Óptica de Fourier · Artículo 02

La difracción es una transformada de Fourier

Cuando la luz pasa por una rendija, el patrón que aparece en la pared es la transformada de Fourier de la forma de esa rendija. La naturaleza calcula — y nosotros aprendemos a leer el resultado.

En el artículo anterior descompusimos imágenes en sumas de rejillas sinusoidales y vimos que la transformada de Fourier es el mapa que nos dice qué frecuencias hay y con qué peso. Ahora viene algo extraordinario: la naturaleza hace exactamente ese cálculo cada vez que la luz pasa por una abertura. El patrón de difracción que aparece en una pantalla lejana es la transformada de Fourier de la forma de esa abertura. Ni más, ni menos.

La luz no viaja en línea recta (del todo)

Parece que la luz viaja en línea recta — y en la mayoría de las situaciones cotidianas funciona perfectamente suponer eso. Pero si miras con cuidado, cuando la luz pasa por una rendija muy estrecha, no se limita a proyectar una franja del mismo ancho al otro lado. Se esparce. El borde de la sombra no es nítido: aparecen franjas claras y oscuras que no deberían estar ahí si la luz fuera sólo rayos rectos.

Christiaan Huygens propuso en 1678 una idea elegante: cada punto de un frente de onda se comporta como una fuente puntual que emite su propia ondita esférica. El frente de onda un instante después es la envolvente de todas esas onditas. Si el frente de onda pasa por una abertura, sólo los puntos dentro de ella contribuyen — y esas onditas interfieren entre sí al propagarse. Esa interferencia es lo que produce el patrón de difracción.

Augustin Fresnel formalizó esto un siglo más tarde con la integral de difracción de Huygens-Fresnel:

U(P)=1iλΣU(Q)eikrrcosθ  dSU(P) = \frac{1}{i\lambda} \iint_{\Sigma} U(Q)\, \frac{e^{ikr}}{r}\, \cos\theta \; dS

Cada punto QQ de la abertura Σ\Sigma emite una ondita eikr/re^{ikr}/r que se propaga hasta el punto PP de la pantalla. El campo total es la suma (integral) de todas esas contribuciones. Parece complicada, pero el punto clave es que esta integral, en el caso más relevante para la óptica, se simplifica a algo que ya conocemos.

La rendija que calcula

Supongamos que la pantalla está lejos de la abertura — lo que los ópticos llaman campo lejano o aproximación de Fraunhofer. En ese límite, la integral de Huygens-Fresnel se reduce, sorprendentemente, a una transformada de Fourier.

De Huygens-Fresnel a la transformada de Fourier

Partimos de la integral de difracción. La distancia desde un punto Q=(x,y)Q = (x', y') de la abertura hasta un punto P=(xP,yP)P = (x_P, y_P) de la pantalla a distancia zz es:

r=z2+(xPx)2+(yPy)2r = \sqrt{z^2 + (x_P - x')^2 + (y_P - y')^2}

Para zz grande, expandimos en serie:

rz+(xPx)2+(yPy)22zr \approx z + \frac{(x_P - x')^2 + (y_P - y')^2}{2z}

Expandiendo el cuadrado: (xPx)2=xP22xPx+x2(x_P - x')^2 = x_P^2 - 2x_P x' + x'^2. El término xP2/2zx_P^2/2z es una fase global (no depende de la posición en la abertura). El término x2/2zx'^2/2z es la fase cuadrática de Fresnel. En la aproximación de Fraunhofer (za2/λz \gg a^2/\lambda), este término es despreciable. Solo queda el término cruzado:

rz+xP22zxPxzr \approx z + \frac{x_P^2}{2z} - \frac{x_P x'}{z}

Sustituyendo en la integral y usando 1/r1/z1/r \approx 1/z:

U(xP)t(x)e2πixPλzxdx=F{t} ⁣(xPλz)U(x_P) \propto \int t(x') \, e^{-2\pi i \frac{x_P}{\lambda z} x'} dx' = \mathcal{F}\{t\}\!\left(\frac{x_P}{\lambda z}\right)

La frecuencia espacial es u=xP/(λz)u = x_P / (\lambda z): posiciones lejanas del centro corresponden a frecuencias altas.

Si llamamos t(x)t(x) a la función de transmisión de la rendija (1 donde la luz pasa, 0 donde está bloqueada), entonces la intensidad en la pantalla lejana es:

I(u)t(x)e2πiuxdx2=F{t}(u)2I(u) \propto \left| \int t(x)\, e^{-2\pi i u x}\, dx \right|^2 = \left|\mathcal{F}\{t\}(u)\right|^2

Para una rendija de ancho aa, esa transformada nos da exactamente la función sinc al cuadrado:

I(u)sinc2 ⁣(πauλf)I(u) \propto \operatorname{sinc}^2\!\left(\frac{\pi\, a\, u}{\lambda f}\right)

Un máximo central brillante, lóbulos secundarios que decaen rápidamente, y ceros donde las contribuciones de las distintas partes de la rendija se cancelan exactamente. Pero no confíes sólo en la fórmula: compruébalo. Mueve el deslizador y observa qué le pasa al patrón de difracción cuando cambias el ancho de la rendija:

Explorar · Rendija simple
Apertura
Intensidad |FT|² Patrón de difracción
Ancho de rendija 40 px

Fíjate en algo contraintuitivo: cuanto más estrecha es la rendija, más ancho es el patrón de difracción. Y al revés: una rendija amplia produce un pico central estrecho. Esto es exactamente la propiedad de reciprocidad de la transformada de Fourier — una función estrecha en el espacio tiene un espectro ancho en frecuencia, y viceversa. La naturaleza nos lo está mostrando con luz.

Dos rendijas: el experimento de Young

Thomas Young, en 1801, hizo pasar luz por dos rendijas paralelas y observó algo que los defensores de la teoría corpuscular no podían explicar: franjas de interferencia. Zonas claras y oscuras alternadas, como si las ondas que salen de cada rendija se reforzaran en unos puntos y se cancelaran en otros.

Desde nuestra perspectiva de Fourier, el experimento de Young no es misterioso — es exactamente lo que esperarías. La función de transmisión ahora son dos rectángulos separados. Su transformada de Fourier es el producto de un sinc (por el ancho de cada rendija) y un coseno (por la separación entre ellas). El coseno produce las franjas de interferencia; el sinc modula su envolvente.

Juega con los dos parámetros: el ancho de cada rendija controla la envolvente (lo rápido que decae), y la separación controla la frecuencia de las franjas de interferencia:

Explorar · Doble rendija
Apertura t(x)
Intensidad |FT|²
Ancho de rendija 12 px
Separación 60 px

Si separas mucho las rendijas, las franjas se aprietan. Si las acercas, se ensanchan. Es la misma reciprocidad de siempre: distancia grande en el espacio produce frecuencia alta en el patrón de difracción. La transformada de Fourier conecta los dos mundos.

De dos rendijas a muchas: la rejilla de difracción

Si dos rendijas producen franjas, ¿qué pasa con 5? ¿Con 20? ¿Con 50? La respuesta es una de las herramientas más útiles de la óptica: la rejilla de difracción.

A medida que aumentas el número de rendijas, los máximos de interferencia se hacen más estrechos y más brillantes. Con pocas rendijas, los picos son anchos y vagos. Con muchas, se convierten en líneas afiladas. La posición de esos máximos la da la ecuación de la rejilla:

dsinθ=mλd\,\sin\theta = m\lambda

donde dd es el periodo de la rejilla, θ\theta el ángulo, mm el orden de difracción, y λ\lambda la longitud de onda. Cada longitud de onda se desvía un ángulo diferente — por eso una rejilla de difracción separa colores, igual que un prisma pero con más control.

Explora cómo cambia el patrón al variar el número de rendijas, su ancho individual y el periodo entre ellas:

Explorar · Rejilla de difracción
Apertura t(x)
Intensidad |FT|²
Rendijas N 8
Ancho de rendija 6 px
Periodo 40 px

Observa con atención: entre los máximos principales aparecen N2N - 2 mínimos secundarios. Con 2 rendijas no los ves; con 50, la estructura fina es evidente. La rejilla de difracción es la base de la espectroscopía — la herramienta que nos dice de qué están hechas las estrellas, qué contaminantes hay en el aire, y qué energías emiten los átomos.

Cualquier apertura

Todo lo anterior es 1D — rendijas infinitamente largas en una dirección. Pero la misma idea se extiende a dos dimensiones sin ningún cambio conceptual. Si la abertura tiene una función de transmisión t(x,y)t(x,y), el patrón de difracción en campo lejano es:

I(u,v)t(x,y)e2πi(ux+vy)dxdy2=F{t}2I(u,v) \propto \left| \iint t(x,y)\, e^{-2\pi i(ux+vy)}\, dx\, dy \right|^2 = \left|\mathcal{F}\{t\}\right|^2

Un cuadrado produce un patrón sinc en dos direcciones (una cruz). Un círculo produce los anillos de Airy — la función de Bessel de primer orden, que es el equivalente circular del sinc. Un triángulo produce un patrón hexagonal. Un anillo produce anillos concéntricos más marcados. Selecciona diferentes aperturas y observa la relación entre la forma y su patrón de difracción:

Explorar · Aperturas 2D
Apertura
Apertura t(x,y)
Patrón |FT|²

Cada patrón te está contando algo sobre la geometría de la abertura. Un patrón alargado en una dirección significa que la abertura es estrecha en esa dirección. Simetrías en la abertura se reflejan en simetrías del patrón. La transformada de Fourier traduce entre forma y frecuencias — y la difracción es esa traducción hecha visible.

La matemática: la integral de difracción

Recapitulemos lo que hemos construido desde la intuición. La integral de difracción de Huygens-Fresnel nos dice el campo en un punto PP de la pantalla:

U(P)=1iλΣU(Q)eikrrcosθ  dSU(P) = \frac{1}{i\lambda} \iint_{\Sigma} U(Q)\, \frac{e^{ikr}}{r}\, \cos\theta \; dS

En la aproximación de Fraunhofer (pantalla lejana), hacemos dos simplificaciones: 1/r1/r es prácticamente constante sobre la abertura, y la fase krkr se puede aproximar de forma lineal en las coordenadas de la abertura. Con esas dos simplificaciones, la integral se reduce a:

U(u,v)t(x,y)e2πi(ux+vy)dxdy=F{t}(u,v)U(u,v) \propto \iint t(x,y)\, e^{-2\pi i(ux + vy)}\, dx\, dy = \mathcal{F}\{t\}(u,v)

donde u=xP/λzu = x_P / \lambda z y v=yP/λzv = y_P / \lambda z son las coordenadas angulares escaladas. La intensidad medida en la pantalla es el módulo al cuadrado del campo:

I(u,v)=U(u,v)2=F{t}(u,v)2I(u,v) = \left|U(u,v)\right|^2 = \left|\mathcal{F}\{t\}(u,v)\right|^2

La condición de Fraunhofer se cumple cuando la distancia a la pantalla zz es mucho mayor que a2/λa^2/\lambda, donde aa es el tamaño de la abertura. En la práctica, esto se puede lograr también con una lente — lo veremos en el artículo siguiente.

¿Y esto para qué sirve?

La conexión entre difracción y transformada de Fourier no es una curiosidad elegante. Es la base de tecnologías que usamos todos los días y de descubrimientos que cambiaron la ciencia:

En el próximo artículo veremos cómo una simple lente convergente calcula la transformada de Fourier de la luz que la atraviesa — sin electricidad, sin silicio, sin algoritmos. La óptica de Fourier deja de ser matemáticas sobre papel y se convierte en algo que puedes construir sobre un banco óptico.

Ejercicios

Ejercicio 1

Usa el simulador de doble rendija de arriba. Fija el ancho de cada rendija y duplica la separación entre ellas (dd). ¿Qué le pasa a la separación entre las franjas de interferencia? Demuestra que la separación angular entre máximos consecutivos es Δθλ/d\Delta\theta \approx \lambda / d, y por tanto se reduce a la mitad cuando dd se duplica.

Solución

Los máximos de interferencia de la doble rendija ocurren cuando dsinθ=mλd\,\sin\theta = m\lambda. La separación angular entre máximos consecutivos (mm y m+1m+1) es:

Δθλd\Delta\theta \approx \frac{\lambda}{d}

Si duplicas dd, Δθ\Delta\theta se reduce a la mitad: las franjas se aprietan. Puedes verificarlo en el simulador: al aumentar la separación, las franjas se hacen más estrechas y más juntas. Es la reciprocidad de la transformada de Fourier: una distancia grande en el espacio produce una frecuencia alta (franjas apretadas) en el patrón de difracción.

Ejercicio 2

Usa el simulador de rendija única. Empieza con un ancho aa grande y observa el patrón de difracción. Ahora reduce aa a la mitad. ¿Qué le pasa al ancho del máximo central? Sabiendo que los primeros ceros del patrón sinc están en sinθ=±λ/a\sin\theta = \pm\lambda/a, calcula el ancho angular del máximo central y confirma tu observación.

Solución

Los primeros ceros del sinc ocurren en sinθ=±λ/a\sin\theta = \pm\lambda/a, así que el ancho angular del lóbulo central es:

Δθcentral=2λa\Delta\theta_{\text{central}} = \frac{2\lambda}{a}

Si reduces aa a la mitad, el ancho del lóbulo central se duplica. Cuanto más estrecha la rendija, más ancho el patrón de difracción. En el simulador se ve claramente: al reducir el ancho de la rendija, la luz se esparce más. Es la propiedad fundamental de reciprocidad de la transformada de Fourier: función estrecha en el espacio = espectro ancho en frecuencias.