2. Propagación a través de una interfase dieléctrica y película delgada.

En este capítulo vamos a analizar la evolución de una onda electromagnética al

propagarse a través de un medio dieléctrico. Se analiza la propagación del campo
primeramente sin considerar pérdidas y posteriormente se incluyen pérdidas debido a la
componente imaginaria de la función dieléctrica. Hacemos el ejercicio de considerar
funciones dieléctricas realistas tomadas de la literatura.8
Se espera observar que los campos reflejados y transmitidos tiendan a la solución
estacionaria, la cual también es calculada por medio de una aplicación de condiciones de
frontera por medio del Método de Matriz de Transferencia.9
Finalmente, reformulamos el MDFDT para adaptarlo al caso cuando el medio
material incluya pérdidas.

2.1 Reflexión y transmisión por una interfase; caso estacionario

Consideramos el caso en el cual una onda plana, de frecuencia ω y vector de onda
ki se aproxima desde el medio 1 de permitividad ε1 y permeabilidad μ1 e incide a la
interfase con un medio 2 (ε2,μ2). Esta condición se ilustra en la Fig. 2.1 y se expresa
mediante las funciones discontinuas en z =z0 como

⎧ ε1 z < zo
⎪
ε ( z) = ⎨
⎪ε (2.1-1)
⎩ 2 z ≥ z0

⎧ μ1 z < zo
⎪
μ ( z) = ⎨
⎪μ (2.1-2)
⎩ 2 z ≥ z0

Cuando una onda incide a una interfase, una fracción del campo incidente es
reflejada y una fracción es transmitida en el medio, tal y como lo ilustran la función

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

22
a la simulación del Campo Electromagnético.
 ⎧ei ( k1 z −ωo t ) + rei ( − k1 z −ωo t ) z < zo
⎪⎪
E x ( z, t ) = ⎨
⎪ (2.2)
⎪⎩ tei ( k 2 z −ωo t ) z ≥ zo

(ε1,μ1) ( ε 2 ,μ 2 )

z0 z
Figura. 2.1 Interfase entre dos medios

Por medio de la aplicación de las condiciones de frontera, la continuidad del

campo eléctrico paralelo (Ex) y campo magnético paralelo (Hy) es posible encontrar los
coeficientes de Fresnel8
k − k2z
r = 1z (2.3-1)
k1z + k 2 z

2k1z
t= (2.3-2)
k1z + k 2 z

En el caso de incidencia normal, estas expresiones se reducen a

ε 1 μ1 − ε 2 μ 2
r= (2.4-1)
ε 1 μ1 + ε 2 μ 2

2 ε 1 μ1
t= (2.4-2)
ε 1 μ1 + ε 2 μ 2

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

23
a la simulación del Campo Electromagnético.
2.2 Reflexión y transmisión por una interfase: aplicación del MDFDT

Para simular la propagación de una onda a través de una interfase entre dos
medios, es necesario describir la interfase del medio material de una forma discreta. El
punto zo es relacionado con el punto discreto k_2 por medio de la relación
z o = k _ 2 * Δz . La discontinuidad de la función dieléctrica se describe por medio de las
relaciones:

⎧ ε1 k <k_2
⎪
ε (k ) = ⎨
⎪ε (2.5-1)
⎩ 2 k ≥k_2

⎧ μ1 k <k_2
⎪
μ (k ) = ⎨
⎪μ (2.5-2)
⎩ 2 k ≥k_2

Para introducir estas relaciones dentro del algoritmo del MDFDT retomamos las
ecuaciones (1.5-1) y (1.5-2)

E xn +1 / 2 (k ) = E xn −1 / 2 (k ) −
0.5
ε
[H n
y (k + 1 / 2) − H yn (k − 1 / 2) ] (2.6-1)

H yn +1 (k + 1 / 2) = H yn (k + 1 / 2) −
0.5
μ
[E n +1 / 2
x (k + 1) − E xn +1 / 2 (k ) ] (2.6-2)

Estas ecuaciones las podemos escribir en código computacional mediante las

relaciones
ex(k) = ex(k) + cb(k)*( hy(k-1) – hy(k) ) (2.7-1)
hy(k) = hy(k) + ca(k)*( ex(k) – ex(k+1) ) (2.7-2)
en las cuales la cantidad que afecta a los paréntesis es
cb(k) = 0.5/ ε
ca(k) = 0.5/ μ
y donde los valores k especifican la coordenada del material.
En estas ecuaciones es posible introducir diferentes valores de permeabilidad
magnética μ , tal y como es el caso de los materiales izquierdos10. Sin embargo, en este
trabajo nos restringiremos al caso en que los medios tienen μ = 1 .

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

24
a la simulación del Campo Electromagnético.
 En la Fig. 2.2 mostramos nuestra simulación de la incidencia de una onda
sinusoidal de un medio de ε 1 = 1 a un medio de función dieléctrica ε 2 = 4 , la longitud de
onda λ1 es del orden de 632.8 x10-9. Se muestran también los coeficientes de Fresnel para
la reflexión y la trasmisión, calculados por las relaciones (2.4).

1981 pasos temporales

2,0

ε1=1 ε2=4
1,5
r=0.336
1,0

0,5
Ex(z,t)

t=0.664

0,0

-0,5

λ2
-1,0

-1,5
λ1

-2,0
0 1 2 3 4 5 6 7 8
-6
z(metros) x10

Figura 2.2.- Propagación en un medio dieléctrico.

La longitud de onda en el medio esta definida por el vector de onda

2π 2π
k 2 z = n2 =
λ λ2

Donde λ y λ2 son la longitud de onda en el vacío y la longitud de onda en el

medio 2. La relación entre ellas es
λ
λ2 =
n2

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

25
a la simulación del Campo Electromagnético.
2.3 Propagación en un medio con pérdidas

Existe una gran variedad de materiales cuya función dieléctrica (ε) no es

puramente real, sino que es un número complejo (ε=εr + iεi), es decir, contiene una parte
real (εr) y una parte imaginaria (εi). Esta parte imaginaria es responsable de las pérdidas
por absorción en el material. Para este tipo de medios, no es posible aplicar directamente
las ecuaciones (2.6) ya que la implementación numérica de estas ecuaciones con
constante dieléctrica compleja da como resultado campos eléctricos no realistas. Es por
ello que para incorporar la parte imaginaria de la función dieléctrica es necesario hacer
una reformulación del algoritmo.
Una de las características mas importantes del MDFDT es la posibilidad de
simular las propiedades ópticas de medios dependientes de las frecuencias,11 los cuales
por lo general incluyen pérdidas.
Usualmente la parte imaginaria de la función dieléctrica es introducida en el
MDFDT por medio de una conductividad, utilizando el termino J de la ecuación de
Ampere Maxwell (ec. 1.1-1), esta es una estrategia probada1,2, sin embargo en el lenguaje
de propiedades ópticas se utiliza la parte imaginaria de la función dieléctrica. Por lo
anterior en esta sección proponemos una forma alternativa de presentar la formulación
del algoritmo para medios con parte imaginaria de la función dieléctrica, o expresado en
otra forma, del índice de refracción que es la forma en que comúnmente es reportado en
libros de parámetros materiales para óptica.8
Nuestra derivación comienza con el argumento de que el termino J de la ecuación
de Ampere-Maxwell (1.1-1) corresponde a una corriente externa y en nuestro caso no
existe, por lo tanto la ecuación de Maxwell-Ampere queda
1 ∂D(x, t )
∇ × H ( x, t ) = (2.8)
c ∂t

El término del lado derecho lo podemos escribir de la forma:

∞
1 ∂D(x, t ) 1 ∂
∫ ε (ω )E(ω)e
1 − iωt
= dω
c ∂t c ∂t 2π (2.9)
−∞

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

26
a la simulación del Campo Electromagnético.
 Se elige una onda que oscila en el tiempo de la forma

E x (x, t ) = Eo e −iω o t (2.10)

Para introducir la ecuación (2.10) en la ec. (2.9) aplicamos la transformada

∞

∫ E(t)e ω dt
1 i t
E (ω ) =
2π
−∞

∫
1
E (ω ) = Eoe −iωot eiωt dt
2π
−∞

∫
Eo
E (ω ) = e( −iωo +ω )t dt
2π
−∞
Recordando la definición de la función δ
∞

∫
1
δ ( x − a) = e α ( x − a ) dα
2π
−∞
Así

E (ω ) = 2π E oδ (−ω o + ω ) (2.11)

Sustituyendo (2.11) en (2.9) tenemos

∞
1 ∂D(x, t ) 1 ∂
∫ ε (ω ){ 2π E δ (−ω
1 − i ωt
= o o + ω ))e dω
c ∂t c ∂t 2π
−∞

1 ∂D( x, t ) Eo ∂
= ε (ωo )e −iωot dω
c ∂t c ∂t

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

27
a la simulación del Campo Electromagnético.
La forma compleja de la función dieléctrica:
ε (ωo ) = ε r (ωo ) + iε i (ωo )

Considerando lo anterior tenemos

1 ∂D( x, t ) ε r (ωo ) ∂ ⎡ ω ε (ω )
= Eoe −iωot ⎤ + o i o ⎡ Eoe −iωot ⎤
c ∂t c ∂t ⎢⎣ ⎥⎦ c ⎢⎣ ⎥⎦

Al Segundo término del lado derecho le ha sido aplicada la derivada temporal.

Retomando la ec. (2.10), la ec. (2.8) se puede escribir
ε (ω ) ∂ ω ε (ω )
∇ × H ( x, t ) = r o E (t ) + o i o E (t )
c ∂t c
Si retomamos la polarización unidimensional donde H = ĵH y y E = iˆE x

tenemos la ecuación escalar

∂ ε ∂ ω ε (ω )
− Hy = r Ex + o i o Ex (2.12)
∂z c ∂t c
Ahora escribimos esta ecuación en forma de diferencias finitas y obtenemos:

H yn (k + 1 / 2) − H yn (k −1 / 2) ε E n +1 / 2 (k ) − Exn −1 / 2 (k ) ωoε i (ωo ) Exn +1 / 2 (k ) + Exn −1 / 2 (k )

− = r x +
Δz c Δt c 2
(2.13)
Nótese que el último término lo hemos aproximado como el promedio entre dos
pasos temporales. Ahora despejamos el término E xn +1 / 2 (k )

⎛ ω ε (ω )Δt ⎞
E xn+1/ 2 (k )⎜⎜1 + o i o
⎝ 2ε r ⎟
⎠
⎛ ω ε (ω )Δt ⎞ 1 Δt
⎟ = E xn−1/ 2 (k )⎜1 − o i o
⎜
⎝ 2ε r
⎟−
⎟
⎠ ε r Δz
(
c H yn (k + 1 / 2) − H yn (k − 1 / 2) )
Retomando el paso cΔt / Δz = 1 / 2 llegamos a la ecuación
⎛ ω ε (ω )Δt ⎞
⎜1 − o i o ⎟
( )
⎜ 2ε r ⎟ 0.5
n +1 / 2
Ex (k ) = ⎝ ⎠ E n −1 / 2 (k ) − H ny (k + 1 / 2) − H yn (k − 1 / 2)
x
⎛ ω ε (ω )Δt ⎞ ⎛ ω ε (ω )Δt ⎞
⎜1 + o i o ⎟ ε r ⎜⎜1 + o i o ⎟
⎜ 2ε r ⎟ 2ε r ⎟
⎝ ⎠ ⎝ ⎠
(2.14-1)

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

28
a la simulación del Campo Electromagnético.
que junto con la condición de Faraday

H yn +1 (k + 1 / 2) = H yn (k + 1 / 2) −
0.5
μ
[E n +1 / 2
x (k + 1) − E xn +1 / 2 (k ) ] (2.14-2)

resuelven completamente el problema de un medio con parte imaginaria de la función

dieléctrica. Estas ecuaciones se implementan mediante las líneas de programas
ex(k) = ca(k) * ex(k) + cb(k) * ( hy(k-1) – hy(k) ) (2.15-1)
hy(k) =hy(k) + 0.5 * ( ex(k) –ex(k-1) ) (2.15-2)
donde
eaf = dt * epsilon_i*ω0 / (2 * epsilon)
ca(k) = (1 – eaf) / (1 + eaf)
cb(k) = 0.5 / (epsilon * (1 + eaf))

En la Fig. 2.3 se ilustra el resultado de la simulación para un material dieléctrico

con pérdidas, para un valor de la función dieléctrica ε=4+i

2381 pasos temporales

2,0

1,5 ε1=1 ε2=4+i

r=0.341
1,0

0,5
Ex(z,t)

t=0.658

0,0

-0,5

-1,0 λ2

-1,5 λ1

-2,0
0 1 2 3 4 5 6 7 8
-6
z(metros) x10

Figura 2.3.- Propagación en un Material dieléctrico con pérdidas

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

29
a la simulación del Campo Electromagnético.
 Los coeficientes de Fresnel están dados por r=3.41+.0542i y t=6.58+0.0542i,
donde ahora estas expresiones son determinadas por números complejos, ya que el índice
de refracción del medio 2 es complejo. En la figura ilustramos la parte real de estos
coeficientes de Fresnel.
El decaimiento de la onda electromagnética puede determinarse analíticamente
por medio de la forma del campo. Según la ec. 2.2 tenemos
Ex ( z , t ) = tei ( k 2 z ( z − z o ) −ωt ) z > zo

Donde k 2 z esta dado por

k2z =
ω
c
(n2 r + in2 i ) (2.16)

la onda tiene las siguientes componentes

⎛ 2π ⎞ ⎛ 2π ⎞
E x ( z , t ) = t exp⎜⎜ i n r ( z − z o )⎟⎟ exp⎜⎜ − ni (z − z o )⎟⎟ exp(− iωt ) (2.17)
⎝ λo ⎠ ⎝ λo ⎠
La primera componente corresponde a una onda oscilatoria y la segunda a un
decaimiento exponencial. Este decaimiento puede ser caracterizado por una distancia de
penetración δ que esta dada por
⎛ 2πni δ ⎞
E x (δ ) = exp(−1) = exp⎜⎜ − ⎟⎟
⎝ λo ⎠
de lo cual obtenemos que δ esta dada por
λo
δ= (2.18)
2πni
Esta relación indica que la distancia de decaimiento es inversamente proporcional
a la parte imaginaria del índice de refracción.

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

30
a la simulación del Campo Electromagnético.
 (ε1,μ1) (ε2,μ2) ( ε 3 ,μ 3 )

0 d z

Figura 2.5 Ilustración de la reflexión en una película delgada.

2.4 Reflexión y transmisión a través de una película delgada

Como un ejemplo más realista o mas complicado, ahora consideramos la reflexión

sobre una película delgada. Esto es, una placa de material limitada en ambos lados por
medios semi-infinitos que es ilustrada en la Fig. 2.5, que pueden tener diferentes
propiedades uno de otro. A la izquierda el de la línea z = 0 consideramos el medio 2 y a
la derecha de z = d consideramos el medio 3. Entre ellos se encuentra el medio 2. La
aplicación de las condiciones de frontera para los campos estacionarios da lugar a la
expresión para los coeficientes de reflexión y transmisión para la película delgada
r12 + r23 e ik 2 z d
r=
1 + r12 r23 e ik 2 z d
Donde r12 y r23 son los coeficientes de Fresnel entre el medio 1 y 2 y 2 y 3
respectivamente. En la Fig. 2.6 mostramos la variación de la reflexión ( R = rr * ) en
función del espesor de la distancia óptica, que es definida como la distancia por el índice
de reflexión de la película delgada. Los medios semi-infinitos al lado izquierdo y derecho
son ε1 = 1.0 y ε3 = 2.25. Consideramos tres valores diferentes para la película delgada ε2
= 1.44, ε2 = 4.0 y ε2 = 9.0.

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

31
a la simulación del Campo Electromagnético.
 ε2 = 9.0
0.5

0.4
Reflexión

0.3

ε2 = 4.0
0.2

0.1

ε2 = 1.44

0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0
Distancia óptica (n2d/λo)

Figura 2.6 Variación de la reflexión en función de la distancia óptica de una película delgada.

ε1=1.0
0.5
ε2=9
ε3=2.25
0.4
Reflexión

0.3

0.2

0.1

0.0
0.0 0.5 1.0 1.5 2.0 2.5 3.0
Frecuencia (ω/ωo)

Figura 2.7 Variación de la reflexión de una película delgada en función de la frecuencia.

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

32
a la simulación del Campo Electromagnético.
La reflexión es periódica según lo indica el factor

ω 2π
β = k 2 z d = n2 d = n2 d
c λo

Esta relación indica que la reflexión será periódica ante cambios tanto en la
variación del espesor de la película (Fig. 2.6), como en la variación de la frecuencia (Fig.
2.7).
En la Fig. 2.8 mostramos un ejemplo de la simulación de la transmisión de una onda
sinusoidal a través de una película delgada de espesor λo/8. La parte (a) ilustra la
variación espacial del índice de refracción. La parte (b) ilustra la propagación de la onda
con el MDFDT. Finalmente la parte (c) ilustra la distribución del campo que es predicha
por la solución estacionaria que es obtenida por medio de aplicación de las condiciones
de frontera.

2,5
Indice de refracción

λo/8
(a)
2,0

1,5

1,0

0,5

1,50
(b)
0,75
Ex(z,t)

0,00

-0,75

-1,50

1,50
(c)
Ex(z,t)

0,75

0,00

-0,75

-1,50

0 1 2 3 4 5 6 7 8
-6
z (metros) x10

Figura 2.8 Propagación del campo electromagnético a través de una película delgada. (a)
variación del índice de refracción. (b) distribución obtenida a partir del MDFDT. (c)
Distribución a partir de la solución del campo estacionario.

Aplicación del Método de Diferencias Finitas en el Dominio del Tiempo

33
a la simulación del Campo Electromagnético.