Difracción de Fresnel - Fresnel diffraction

En óptica , la ecuación de difracción de Fresnel para la difracción de campo cercano es una aproximación de la difracción de Kirchhoff-Fresnel que se puede aplicar a la propagación de ondas en el campo cercano . Se utiliza para calcular el patrón de difracción creado por las ondas que pasan a través de una apertura o alrededor de un objeto, cuando se ven desde relativamente cerca del objeto. En contraste, el patrón de difracción en la región de campo lejano viene dado por la ecuación de difracción de Fraunhofer .

El campo cercano se puede especificar mediante el número de Fresnel , F , de la disposición óptica. Cuando se considera que la onda difractada está en el campo cercano. Sin embargo, la validez de la integral de difracción de Fresnel se deduce de las aproximaciones derivadas a continuación. Específicamente, los términos de fase de tercer orden y superiores deben ser insignificantes, una condición que puede escribirse como ${\ Displaystyle F \ gg 1}$

$${\ Displaystyle {\ frac {F \ theta ^ {2}} {4}} \ ll 1,}$$

donde es el ángulo máximo descrito por , a y L el mismo que en la definición del número de Fresnel . ${\ Displaystyle \ theta}$${\ Displaystyle \ theta \ approx a / L}$

Difracción de Fresnel que muestra el punto central de Arago

La difracción de Fresnel múltiple en crestas periódicas poco espaciadas ( espejo estriado ) provoca la reflexión especular ; este efecto se puede utilizar para espejos atómicos .

Tratamientos tempranos de este fenómeno

Algunos de los primeros trabajos sobre lo que se conocería como difracción de Fresnel fueron realizados por Francesco Maria Grimaldi en Italia en el siglo XVII. En su monografía titulada "Luz", Richard C. MacLaurin explica la difracción de Fresnel preguntando qué sucede cuando la luz se propaga y cómo se ve afectado ese proceso cuando una barrera con una rendija o un agujero se interpone en el haz producido por una fuente distante de luz. Utiliza el Principio de Huygens para investigar, en términos clásicos, lo que sucede. El frente de onda que procede de la hendidura y en una pantalla de detección a cierta distancia se aproxima muy de cerca a un frente de onda que se origina en el área del espacio sin tener en cuenta ninguna interacción mínima con el borde físico real.

El resultado es que si el espacio es muy estrecho, solo se pueden producir patrones de difracción con centros brillantes. Si la brecha se ensancha progresivamente, los patrones de difracción con centros oscuros se alternarán con patrones de difracción con centros brillantes. A medida que la brecha se hace más grande, las diferencias entre las bandas oscuras y claras disminuyen hasta que ya no se puede detectar un efecto de difracción.

MacLaurin no menciona la posibilidad de que el centro de la serie de anillos de difracción producidos cuando la luz pasa a través de un pequeño agujero sea negro, pero sí apunta a la situación inversa en la que la sombra producida por un pequeño objeto circular puede, paradójicamente, tener un brillo. centro . (pág.219)

En su Óptica , Francis Weston Sears ofrece una aproximación matemática sugerida por Fresnel que predice las características principales de los patrones de difracción y usa solo matemáticas simples. Al considerar la distancia perpendicular desde el agujero en una pantalla de barrera a una pantalla de detección cercana junto con la longitud de onda de la luz incidente, es posible calcular una serie de regiones llamadas elementos de medio período o zonas de Fresnel . La zona interior es un círculo y cada zona sucesiva será un anillo anular concéntrico. Si el diámetro del orificio circular en la pantalla es suficiente para exponer la primera zona de Fresnel o la central, la amplitud de la luz en el centro de la pantalla de detección será el doble de lo que sería si la pantalla de detección no estuviera obstruida. Si el diámetro del agujero circular en la pantalla es suficiente para exponer dos zonas de Fresnel, entonces la amplitud en el centro es casi cero. Eso significa que un patrón de difracción de Fresnel puede tener un centro oscuro. Estos patrones se pueden ver y medir, y corresponden bien a los valores calculados para ellos.

La integral de difracción de Fresnel

Geometría de difracción, que muestra el plano de apertura (u objeto difractante) y el plano de la imagen, con sistema de coordenadas.

El patrón de difracción de campo eléctrico en un punto (x, y, z) viene dado por:

$${\ Displaystyle E \ left (x, y, z \ right) = {\ frac {1} {i \ lambda}} \ iint _ {- \ infty} ^ {+ \ infty} {E \ left (x ', y ', 0 \ right) {\ frac {e ^ {ikr}} {r}}} {\ frac {z} {r}} dx'dy'}$$

dónde

• ${\ Displaystyle E \ left (x ', y', 0 \ right) \,}$ es el campo eléctrico en la apertura;
• ${\ Displaystyle r = {\ sqrt {\ left (x-x '\ right) ^ {2} + \ left (y-y' \ right) ^ {2} + z ^ {2}}} \,}$;
• ${\ Displaystyle k \,}$es el número de onda ; y${\ Displaystyle 2 \ pi / \ lambda \,}$
• ${\ Displaystyle i \,}$es la unidad imaginaria .

La solución analítica de esta integral es imposible para todas las geometrías de difracción, excepto para las más simples. Por lo tanto, generalmente se calcula numéricamente.

La aproximación de Fresnel

Comparación entre el patrón de difracción obtenido con la ecuación de Rayleigh-Sommerfeld, la aproximación de Fresnel (paraxial) y la aproximación de Fraunhofer (campo lejano).

El principal problema para resolver la integral es la expresión de r . Primero, podemos simplificar el álgebra introduciendo la sustitución:

$${\ Displaystyle \ rho ^ {2} = \ left (x-x '\ right) ^ {2} + \ left (y-y' \ right) ^ {2} \,}$$

Sustituyendo en la expresión de r , encontramos:

$${\ Displaystyle r = {\ sqrt {\ rho ^ {2} + z ^ {2}}} = z {\ sqrt {1 + {\ frac {\ rho ^ {2}} {z ^ {2}}} }}}$$

A continuación, por la expansión binomial,

$${\ Displaystyle {\ sqrt {1 + u}} = (1 + u) ^ {\ frac {1} {2}} = 1 + {\ frac {u} {2}} - {\ frac {u ^ { 2}} {8}} + \ cdots}$$

Podemos expresar como ${\ Displaystyle r}$

$${\ Displaystyle {\ begin {alineado} r & = z {\ sqrt {1 + {\ frac {\ rho ^ {2}} {z ^ {2}}}}} \\ & = z \ left [1+ { \ frac {\ rho ^ {2}} {2z ^ {2}}} - {\ frac {1} {8}} \ left ({\ frac {\ rho ^ {2}} {z ^ {2}} } \ right) ^ {2} + \ cdots \ right] \\ & = z + {\ frac {\ rho ^ {2}} {2z}} - {\ frac {\ rho ^ {4}} {8z ^ { 3}}} + \ cdots \ end {alineado}}}$$

Si consideramos todos los términos de series binomiales, entonces no hay aproximación. Sustituyamos esta expresión en el argumento del exponencial dentro de la integral; la clave de la aproximación de Fresnel es asumir que el tercer término es muy pequeño y puede ignorarse y, en adelante, cualquier orden superior. Para que esto sea posible, debe contribuir a la variación de la exponencial para un término casi nulo. En otras palabras, tiene que ser mucho más pequeño que el período del exponencial complejo; es decir : ${\ Displaystyle 2 \ pi}$

$${\ Displaystyle k {\ frac {\ rho ^ {4}} {8z ^ {3}}} \ ll 2 \ pi}$$

expresando k en términos de la longitud de onda,

$${\ Displaystyle k = {\ frac {2 \ pi} {\ lambda}}}$$

obtenemos la siguiente relación:

$${\ Displaystyle {\ frac {\ rho ^ {4}} {z ^ {3} \ lambda}} \ ll 8}$$

Multiplicando ambos lados por , tenemos ${\ Displaystyle z ^ {3} / \ lambda ^ {3}}$

$${\ Displaystyle {\ frac {\ rho ^ {4}} {\ lambda ^ {4}}} \ ll 8 {\ frac {z ^ {3}} {\ lambda ^ {3}}}}$$

o, sustituyendo la expresión anterior por , ${\ Displaystyle \ rho ^ {2}}$

$${\ Displaystyle {\ frac {1} {\ lambda ^ {4}}} \ left [\ left (x-x '\ right) ^ {2} + \ left (y-y' \ right) ^ {2} \ right] ^ {2} \ ll 8 {\ frac {z ^ {3}} {\ lambda ^ {3}}}}$$

Si esta condición se cumple para todos los valores de x , x ' , y e y' , entonces podemos ignorar el tercer término en la expresión de Taylor. Además, si el tercer término es insignificante, entonces todos los términos de orden superior serán incluso más pequeños, por lo que también podemos ignorarlos.

Para aplicaciones que involucran longitudes de onda ópticas, la longitud de onda λ es típicamente muchos órdenes de magnitud menor que las dimensiones físicas relevantes. En particular:

$${\ Displaystyle \ lambda \ ll z}$$

y

$${\ Displaystyle \ lambda \ ll \ rho}$$

Por lo tanto, en la práctica, la desigualdad requerida siempre será cierta siempre que

$${\ Displaystyle \ rho \ ll z}$$

Entonces podemos aproximar la expresión con solo los dos primeros términos:

$${\ Displaystyle r \ approx z + {\ frac {\ rho ^ {2}} {2z}} = z + {\ frac {\ left (x-x '\ right) ^ {2} + \ left (y-y' \ right) ^ {2}} {2z}}}$$

Esta ecuación, entonces, es la aproximación de Fresnel , y la desigualdad indicada anteriormente es una condición para la validez de la aproximación.

Difracción de Fresnel

La condición de validez es bastante débil y permite que todos los parámetros de longitud tomen valores comparables, siempre que la apertura sea pequeña en comparación con la longitud del camino. Para la r en el denominador, vamos un paso más allá y lo aproximamos solo con el primer término ,. Esto es válido en especial si estamos interesados en el comportamiento del campo sólo en una pequeña zona cerca del origen, donde los valores de x e y son mucho menores que z . En general, la difracción de Fresnel es válida si el número de Fresnel es aproximadamente 1. ${\ Displaystyle r \ approx z}$

Para la difracción de Fresnel, el campo eléctrico en el punto viene dado por: ${\ Displaystyle (x, y, z)}$

$${\ Displaystyle E (x, y, z) = {\ frac {e ^ {ikz}} {i \ lambda z}} \ iint _ {- \ infty} ^ {+ \ infty} E \ left (x ', y ', 0 \ right) e ^ {{\ frac {ik} {2z}} \ left [\ left (x-x' \ right) ^ {2} + \ left (y-y '\ right) ^ { 2} \ right]} dx'dy '}$$

Difracción de Fresnel de apertura circular, trazada con funciones de Lommel

Esta es la integral de difracción de Fresnel; significa que, si la aproximación de Fresnel es válida, el campo de propagación es una onda esférica, que se origina en la apertura y se mueve a lo largo de z . La integral modula la amplitud y fase de la onda esférica. La solución analítica de esta expresión solo es posible en casos excepcionales. Para un caso más simplificado, válido solo para distancias mucho mayores desde la fuente de difracción, consulte Difracción de Fraunhofer . A diferencia de la difracción de Fraunhofer, la difracción de Fresnel tiene en cuenta la curvatura del frente de onda para calcular correctamente la fase relativa de las ondas interferentes.

Formas alternativas

Circunvolución

La integral se puede expresar de otras formas para calcularla usando algunas propiedades matemáticas. Si definimos la siguiente función:

$${\ Displaystyle h (x, y, z) = {\ frac {e ^ {ikz}} {i \ lambda z}} e ^ {i {\ frac {k} {2z}} \ left (x ^ {2 } + y ^ {2} \ right)}}$$

entonces la integral se puede expresar en términos de una convolución :

$${\ Displaystyle E (x, y, z) = E (x, y, 0) * h (x, y, z)}$$

en otras palabras, estamos representando la propagación utilizando un modelo de filtro lineal. Es por eso que podríamos llamar a la función la respuesta al impulso de la propagación en el espacio libre. ${\ Displaystyle h (x, y, z)}$

Transformada de Fourier

Otra forma posible es a través de la transformada de Fourier . Si en la integral expresamos k en términos de longitud de onda:

$${\ Displaystyle k = {\ frac {2 \ pi} {\ lambda}}}$$

y expanda cada componente del desplazamiento transversal:

$${\ Displaystyle {\ begin {alineado} \ left (x-x '\ right) ^ {2} & = x ^ {2} + {x'} ^ {2} -2xx ', \\\ izquierda (y- y '\ right) ^ {2} & = y ^ {2} + {y'} ^ {2} -2yy ', \ end {alineado}}}$$

entonces podemos expresar la integral en términos de la transformada de Fourier bidimensional. Usemos la siguiente definición:

$${\ Displaystyle G (p, q) = {\ mathcal {F}} \ left \ {g (x, y) \ right \} \ equiv \ iint _ {- \ infty} ^ {\ infty} g (x, y) e ^ {- i2 \ pi (px + qy)} \, dx \, dy,}$$

donde p y q son frecuencias espaciales ( números de onda ). La integral de Fresnel se puede expresar como

$${\ displaystyle {\ begin {alineado} E (x, y, z) & = \ left. {\ frac {e ^ {ikz}} {i \ lambda z}} e ^ {i {\ frac {\ pi} {\ lambda z}} \ left (x ^ {2} + y ^ {2} \ right)} {\ mathcal {F}} \ left \ {E (x ', y', 0) e ^ {i { \ frac {\ pi} {\ lambda z}} \ left ({x '} ^ {2} + {y'} ^ {2} \ right)} \ right \} \ right | _ {p = {\ frac {x} {\ lambda z}}, \ q = {\ frac {y} {\ lambda z}}} \\ & = h (x, y) \ cdot G (p, q) {\ big |} _ {p = {\ frac {x} {\ lambda z}}, \ q = {\ frac {y} {\ lambda z}}}, \ end {alineado}}}$$

Es decir, primero se multiplica el campo para ser propagados por una exponencial compleja, calcular su transformada de Fourier bidimensional, en lugar de con y multiplicarlo por otro factor. Esta expresión es mejor que las otras cuando el proceso conduce a una transformada de Fourier conocida, y la conexión con la transformada de Fourier se estrecha en la transformación canónica lineal , que se analiza a continuación. ${\ Displaystyle (p, q)}$${\ Displaystyle \ left ({\ tfrac {x} {\ lambda z}}, {\ tfrac {y} {\ lambda z}} \ right)}$

Transformación canónica lineal

Desde el punto de vista de la transformación canónica lineal , la difracción de Fresnel puede verse como una cizalladura en el dominio del tiempo-frecuencia , correspondiente a cómo la transformada de Fourier es una rotación en el dominio del tiempo-frecuencia.

Ver también

• Difracción de Fraunhofer
• Integral de Fresnel
• Zona de Fresnel
• Número de Fresnel
• Augustin-Jean Fresnel
• Espejo estriado
• Generador de imágenes de Fresnel
• Espiral de Euler

Notas

Referencias

• Goodman, Joseph W. (1996). Introducción a la óptica de Fourier . Nueva York: McGraw-Hill . ISBN 0-07-024254-2.