Optica: Interferencia de haces multiples

La interferencia de haces multiples se da cuando un número muy grande de ondas mutuamente coherentes se hacen interferir.

El método más común de producir este número de ondas mutuamente coherentes es por división de amplitud. Esta división ocurre por reflexión múltiple entre dos superficies paralelas parcialmente reflectantes.

El primer rayo es parcialmente reflejado y parcialmente transmitido en la primer superficie. La parte transmitida es subsecuentemente reflejada hacia atras y hacia adelante entre las dos superficies.

Comenzemos el análisis matemático calculando la diferencia de camino óptico que se puede obtener de la figura (1) donde:

Archivo:55555.jpg

Figura 1.Diagrama de la diferencia de camino optico entre dos rayos sucesivos

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \sin\theta_{0}=\frac{BA'}{C'B}=\frac{BA'}{CP}=\frac{BA'}{2d\tan\theta}

En donde se puede ver geometricamente que la diferencia entre dos rayos sucesivos transmitidos es $2d\cos \theta$, donde d es la separación entre las dos superficies reflectantes y Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \theta es el angulo entre cualquier rayo interno reflejado y la superficie normal.

$BA'=2d\tan \theta \sin \theta _{0}$

lo cual se pude demostrar que es igual a:

$\Delta \Lambda ={\frac {2d}{\cos \theta }}nt-2d\tan \theta \sin \theta _{0}ni$

factorizando

$={\frac {2d}{\cos \theta }}nt\left[1-\sin \theta \sin \theta _{0}{\frac {ni}{nt}}\right]$

obtenemos la relación de Snell

$={\frac {2d}{\cos \theta }}nt\left[1-\sin ^{2}\theta \right]$

de donde se obtiene:

$\Delta \Lambda =2dnt\cos \theta$

que es la diferencia de fase correspondiente entre dos rayos sucesivos

$\delta ={\frac {2\pi }{\lambda _{0}}}\Delta \Lambda ={\frac {4\pi }{\lambda _{0}}}ntd\cos \theta$

donde Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \lambda es la longitud de onda en el medio y nt es el indice de refracción de el medio entre las superficies reflectantes.

AMPLITUD TOTAL REFLEJADA

---

Ahora definamos las variables.

r es el coeficiente de reflexión de las amplitudes desde el interior del interferometro.

t es el coeficiente de transmisión de las amplitudes para un rayo que viene del exterior.

t' es el coeficiente de transmisión para las amplitudes para un rayo que sale de la cavidad del interferometro.

d es la diferencia de fase entre dos reflexiones consecutivas.

Figura 2. Caminos de los rayos de luz de multiples reflexiones entre dos superficies paralelas.

Tomando la diferencia de fase anterior y sumando las amplitudes de los rayos transmitidos tenemos entonces

$E_{1r}=E_{0}re^{i\omega t}$

$E_{2r}=E_{0}tt'r'e^{i\left(\omega t-\delta \right)}$

$E_{3r}=E_{0}tt'r'^{3}e^{i\left(\omega t-2\delta \right)}$

$E_{Nr}=E_{0}tt'r'^{(2N-3)}e^{i\left(\omega t-\left(N-1\right)\delta \right)}$

donde $E_{0}e^{i\omega t}$ es la onda incidente.

Los términos Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \delta,2\delta,...,(N-1)\delta son las contribuciones a la fase procedentes de una diferencia de longitud de camino óptico entre rayos adyacentes.

$E_{rT}=\sum E_{ir}=E_{0}re^{i\omega t}+\sum _{j=2}^{N}E_{0}tr'^{(2j-3)}t'e^{i\left(\omega t-\left(j-1\right)\delta \right)}$

$=E_{0}e^{i\omega t}\left[r+\left\{\sum _{j=2}^{N}\left(r'^{2}e^{i\delta }\right)^{j-2}\right\}r'tt'e^{-i\delta }\right]$

Si $[r'^{2}e^{i\delta }]<1$, y si el número de términos en la serie se aproxima al infinito, la serie converge. La onda resultante se transforma en:

$E_{r}=E_{0}e^{i\omega t}\left[r+{\frac {r'tt'e^{-i\delta }}{1-r'^{2}e^{-i\delta }}}\right]$

En el caso de absorción cero, cuando no se extrae energía de las ondaas, se puede utilizar las relaciones Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): tt'=1-r^{2}, y r=-r' para volver a escribir la ecuación como

$E_{r}=E_{0}e^{i\omega t}\left[{\frac {r\left(1-e^{-i\delta }\right)}{1-r^{2}e^{-i\delta }}}\right]$

La irradiancia se obtiene tomando el cuadrado complejo de esta amplitud quedando

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): I_{R}=E_{r}E_{r}^{*}=E_{0}^{2}\left[\frac{r\left(1-e^{i\delta}\right)}{1-r^{2}e^{i\delta}}\right]\left[\frac{r\left(1-e^{i\delta}\right)}{1-r^{2}e^{i\delta}}\right]

Que puede transformarse en

$I_{r}=E_{0}^{2}{\frac {2r^{2}\left(1-\cos \theta \right)}{\left(1+r^{4}\right)-2r^{2}\cos \delta }}$

AMPLITUD TOTAL TRANSMITIDA

---

De forma parecida, las amplitudes de las ondas transmitidas proporcionadas por

$E_{1t}=E_{0}tt'e^{i\omega t}$

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): E_{2t}=E_{0}tt'r'^{2}e^{i\left(\omega t-\delta\right)}

$E_{3t}=E_{0}tt'r'^{4}e^{i\left(\omega t-2\delta \right)}$

$E_{Nr}=E_{0}tt'r'^{2\left(N-1\right)}e^{i\left(\omega t-\left(N-1\right)\delta \right)}$

$E_{tT}=E_{0}e^{i\omega t}\left[tt'\sum _{j=1}^{N}r'^{(j-1)}e^{-i\left(j-1\right)\delta }\right]$

Pueden sumarse para dar como resultado

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): E_{t}=E_{0}e^{i\omega t}\frac{tt'}{1-r^{2}e^{-i\delta}}

Multiplicando también por su complejo conjugado

$I_{t}=I_{i}\left(tt'\right)^{2}{\frac {1}{1-r^{2}e^{-i\delta }}}{\frac {1}{1-r^{2}e^{i\delta }}}$

se obtiene la irradiancia del haz transmitido

$I_{t}=I_{i}\left(tt'\right)^{2}{\frac {1}{\left(1+r^{4}\right)-2r^{2}\cos \delta }}$

o bien

$I_{t}=I_{i}{\frac {\left(1-r^{2}\right)^{2}}{\left(1+r^{4}\right)-2r^{2}\cos \delta }}$

Notemos que para dieléctricos

$I_{r}+I_{t}={\frac {\left(1-r^{2}\right)^{2}+2r^{2}\left(1-\cos \delta \right)}{\left(1+r^{4}\right)-2r^{2}\cos \delta }}I_{i}$

expandiendo la ecuación y despejando

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): =\frac{1+r^{4}-2r^{2}-2r^{2}\cos\delta+2r^{2}}{\left(1+r^{4}\right)-2r^{2}\cos\delta}I_{i}=I_{i}

si no se absorbe nada de la energía incidente, la densidad de flujo de la onda incidente será exactamente igual a la suma de la densidad de flujo reflejada por la película más la densidad de flujo total transmitida al salir de la película, es decir

$I_{i}=I_{i}=I_{r}+I_{t}$

si $\delta =2\pi m$ Exixtirá un máximo

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): I_{tmax}=I_{i}\frac{\left(1-r^{2}\right)^{2}}{1+r^{4}-2r^{2}}=I_{i}

$I_{rmin}=0$ (por conservación de energía)

ahora si Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \delta=\left(2m+1\right)\pi se producirá un minimo en la densidad del flujo transmitido

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): I_{tmin}=I_{i}\frac{\left(1-r^{2}\right)^{2}}{(1+r^{2})^{2}}

El máximo correspondiente en la densidad de flujo reflejado es

$I_{rmax}=I_{i}{\frac {4r^{2}}{(1+r^{2})^{2}}}$

Si escribimos

${\frac {I_{t}}{I_{i}}}={\frac {(1-r^{2})^{2}}{1+r^{4}-2r^{2}\left(\cos ^{2}{\frac {\delta }{2}}-\sin ^{2}{\frac {\delta }{2}}\right)}}$

Mediante la identidad trigonometrica $\sin ^{2}\theta +\cos ^{2}\theta =1$

Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \frac{(1-r^{2})^{2}}{1+r^{4}-2r^{2}\left(1-2\sin^{2}\frac{\delta}{2}\right)}

desarrollando el algebra y agrupando terminos

${\frac {\left(1-r^{2}\right)^{2}}{\left(1-r^{2}\right)^{2}+4r^{2}\sin ^{2}{\frac {\delta }{2}}}}$

simplificando se obtiene

${\frac {1}{1+\left[{\frac {2r}{1-r^{2}}}\sin {\frac {\delta }{2}}\right]^{2}}}$

que es llamado coeficiente de fineza

$\digamma =\left({\frac {2r}{1-r^{2}}}\right)^{2}={\frac {4R}{(1-R)^{2}}}$

con lo cual estas ecuaciones se podran escribir como

${\frac {I_{r}}{I_{t}}}={\frac {\digamma \sin ^{2}{\frac {\delta }{2}}}{1+\digamma \sin ^{2}{\frac {\delta }{2}}}}$

y

${\frac {I_{t}}{I_{i}}}={\frac {1}{1+\digamma \sin ^{2}{\frac {\delta }{2}}}}\equiv A\left(\delta \right)$

En donde el término ${\frac {1}{1+\digamma \sin ^{2}{\frac {\delta }{2}}}}\equiv A\left(\delta \right)$ se denomina función de Airy y representa la distribución de la densidad de flujo transmitida

FUNCION DE AIRY

La funcion de Airy representa la distribución de la densidad de flujo transmitida. Cuando Error al representar (SVG (MathML puede ser habilitado mediante un plugin de navegador): respuesta no válida («Math extension cannot connect to Restbase.») del servidor «https://en.wikipedia.org/api/rest_v1/»:): \frac{\delta}{2}=\pi m la función de Airy es igual a la unidad para todos los valores de F y por lo tanto de r. Al aproximarse esta a uno, la densidad de flujo transmitido es muy pequeño, exepto dentro de máximos agudos centrados en los puntos ${\frac {\delta }{2}}=\pi m$. La interferencia de haces multiples conduce a una redistribución de la densidad de energía en comparación con la distribución sinusoidal de dos haces

INTERFEROMETRO DE FABRY-PEROT

---

Este tipo de interferometro forma las franjas con base en multiples reflexiones en dos superficies palnas y paralelas, por lo tanto la interferencia no es solo entre 2 frentes de onda, sino entre un numero muy grande de ellas. A fin de lograr esta reflexiones multiples, las superficies se recubren con una capa reflectora ya sea metalica o dielectrica. Con este tipo de interferencia multiple las franjas ya no tienen perfil senoidal como cuando son solo 2 frentes de onda sino que tienen un perfil muy angosto, aumentando así la presición con que se puede medir su forma y posición. En este interferometro las franjas son de igual inclinación usando una fuente extendida.

las peliculas de metal parcialmente transparentes que se emplean frecuentemente para aumentar la reflectancia (R=r2) absorberan una fracción A de la densidad de flujo, esta fracción se denomina absortancia.

La expresión

$T+R=1$

donde T es la transmitancia, se reescribe como

$T+R+A=1$

Las películas metalicas presentan un desplazamiento adicional de fase ++++, que puede diferir de cero o pi. El desfase entre dos ondas transmitida ssucesivamente es entonces

$\delta ={\frac {4\pi }{\lambda }}n_{t}d\cos \theta +2\delta _{r}$

Para las condiciones que estan siendo consideradas, teta es pequeño y $\pi$ puede considerarse constante.Por lo general, d es tan grande y $\lambda$ tan pequeño que $\pi$ puede omitirse. entonces la ecuación **** puede expresarse como

${\frac {I_{t}}{I_{i}}}={\frac {T^{2}}{1+R^{2}-2R\cos \delta }}$

o de manera equivalente

${\frac {I_{t}}{I_{i}}}={\frac {T^{2}}{(1-R^{2})}}A(\delta )$

Utilizando la ecuación (*) y la definición de la función de Airy, obtenemos

${\frac {I_{t}}{I_{i}}}={\frac {\left(1-R-A\right)^{2}}{\left(1-R\right)^{2}}}A\left(\delta \right)$

simplificando

$=\left(1-{\frac {A}{1-R}}\right)^{2}A\left(\delta \right)$

Dado que la parte absorbida A no es nunca cero, los máximos de la densidad de flujo transmitido siempre serán algo menor que . Por lo tanto, el máximo de transmisión se define como

$I_{i}\left(1-{\frac {A}{1-R}}\right)^{2}\equiv I_{tmax}$

La irradiancia relativa de la distribución de franjas estará determinada por la función de Airy, ya que

${\frac {I_{t}}{i_{tmax}}}=A\left(\theta \right)={\frac {1}{1+\digamma \sin ^{2}{\frac {\delta }{2}}}}$

Los máximos en la transmisión se dan para valores especificos del desfase

$\delta _{max}=2\pi m$

Por tanto, la irradiancia disminuira a la mitad su valor máximo cada vez que 11111111 dado que $\gamma =2{\frac {\delta _{i}}{2}}$

${\frac {1}{2}}={\frac {1}{1+\digamma \sin \left({\frac {{\frac {\delta }{2}}+\delta _{max}}{2}}\right)}}$

de donde se ve que

$\sin \left({\frac {\delta _{\frac {1}{2}}}{2}}+m\pi \right)=\pm sin\left({\frac {\delta _{\frac {1}{2}}}{2}}\right)$

elevando al cuadrado ambos lados de la ecuación

$\sin ^{2}\left({\frac {\delta _{\frac {1}{2}}}{2}}+m\pi \right)=\sin ^{2}\left({\frac {\delta _{\frac {1}{2}}}{2}}\right)$

entonces

$\sin ^{2}\left({\frac {\delta _{\frac {1}{2}}}{2}}\right)={\frac {1}{\digamma }}$

Eliminando el termino cuadrado de la ecuación de lado derecho y despejando obtenemos

$\delta _{\frac {1}{2}}=2\arcsin \left({\sqrt {\frac {1}{\digamma }}}\right)$

Sustituyendo el valor de \gamma

$\gamma =4\arcsin \left({\frac {1}{\sqrt {\digamma }}}\right)$

Dado que F es, por lo general muy grande,

$\arcsin \left({\sqrt {\frac {1}{\digamma }}}\right)\approxeq \left({\sqrt {\frac {1}{\digamma }}}\right)$

entonces se obtiene

${\frac {\delta _{\frac {1}{2}}}{2}}={\frac {1}{\sqrt {\digamma }}}$

y por lo tanto el ancho medio $\gamma =2{\frac {\delta _{i}}{2}}$ es igual a

$\gamma \approxeq \left({\frac {4}{\sqrt {\digamma }}}\right)$

Otro valor interesante es la relacon de la separación de los máximos adyacentes con respecto al ancho medio. Denominada fineza $\Upsilon \equiv {\frac {2\pi }{\gamma }}={\frac {\pi }{2}}{\sqrt {\digamma }}$