Benemerita Universidad Autonoma de Puebla

Facultad de Ciencias Fsico-Matematicas

Luz atrapada en una estructura construida con
metamateriales
Tesis presentada al
Colegio de Fsica
como requisito parcial para la obtencion del grado de
Licenciado en Fsica
por
Esperanza Juarez Perez

asesorada por
Dra. Martha Palomino Ovando

Puebla , Pue. ,
Abril de 2014.

Benemerita Universidad Autonoma de Puebla

Facultad de Ciencias Fsico-Matematicas
Luz atrapada en una estructura construida con
metamateriales
Tesis presentada al
Colegio de Fsica
como requisito parcial para la obtencion del grado de
Licenciado en Fsica
por
Esperanza Juarez Perez

asesorada por
Dra. Martha Palomino Ovando

Puebla, Pue.,
Abril de 2014.

Ttulo: Luz atrapada en una estructura construida con

metamateriales
rez Pe
rez
Estudiante:Esperanza Jua

COMITE

Dr. Rosendo L. Lozada Morales

Presidente

Dra. Rosario Pastrana Sanchez

Secretario

M. C. Pedro Tolentino Eslava

Vocal

Dr. Martn R. Palomino Merino

Suplente

Dra. Martha Palomino Ovando

Asesor

Indice general
Agradecimientos

VII

Resumen

IX

Introducci
on

XI

Objetivos

XV

1. METAMATERIALES
1.1. Definici
on y caracteristicas de un metamaterial
1.1.1. Definici
on e historia . . . . . . . . . . . . . . .
1.1.2. Propiedades . . . . . . . . . . . . . . . . . . . .
1.2. Aplicaciones . . . . . . . . . . . . . . . . . . . . . . .
1.3. Tipos de metamateriales . . . . . . . . . . . . . . .
1.3.1. Dielectricos artificiales . . . . . . . . . . . . . .
1.3.2. Materiales magneticos artificiales . . . . . . . .
1.3.3. Materiales quirales . . . . . . . . . . . . . . . .
1.3.4. Materiales con recubrimiento . . . . . . . . . .
1.4. Modelo de Gupta . . . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

2. DESARROLLO TEORICO
2.1. Ecuaciones de Maxwell . . . . . . . . . . . . . . . . . . .
2.1.1. Fundamentos de las ecuaciones de Maxwell . . . . .
2.1.2. Obtenci
on de las ecuaciones de Maxwell . . . . . . .
2.2. M
etodo de la matriz de transferencia . . . . . . . . .
2.2.1. Ecuaci
on de onda . . . . . . . . . . . . . . . . . . . .
2.2.2. Polarizaci
on P (Polarizacion Transversal Magnetica)
2.2.3. Polarizacion S (Polarizacion Transversal Electrica) .
2.3. C
alculo de tiempo de fase . . . . . . . . . . . . . . . . .

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.
.
.

1
1
1
2
4
5
5
6
6
7
7

.
.
.
.
.
.
.
.

9
9
10
13
15
15
16
20
23

3. RESULTADOS
27
3.1. Espectros de transmisi
on . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.2. Tiempos de fase . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
Conclusiones

35

iii

Indice de figuras
1.1. El agua con ndice de refraccion negativo si se introduce un popote se vera doblar en
direcci
on contraria a la norma . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.2. Efecto Doppler en (a) un material convencional y (b) un material izquierdo . . . . . . .
1.3. El telefono movil de la figura (a) (modelo BL40 de la marca LG) fue el primer dispositivo

3
4

de este tipo en integrar una antena basada en metamateriales (modelo MTM-C de la

empresa norteamericana Rayspan) (b). Esta tecnologa permite reducir a una quinta
parte el tama
no de antena que sera necesario en caso de utilizar tecnologa convencional

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
1.4. Rayos de luz atravesando objetos invisibles. . . . . . . . . . . . . . . . . . . . . . .
1.5. Grafica de () y () modelo de Gupta . . . . . . . . . . . . . . . . . . . . . . .
de antenas impresas.

2.1.
2.2.
2.3.
2.4.
2.5.
2.6.
2.7.

. .
.
Ley de Ampere . . . . .
~ . .
Lneas de fuerza de B
Polarizaci
on P . . . . . .

.
.
.
.
.
Onda incidiendo en la red .
Polarizaci
on S . . . . . . .

Dispositivo de Faraday

Fen
omeno de inducci
on

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

5
7
8
10
11
12
13
17
18
21

3.1. Forma que van a tener las multicapas a) caso periodico y b) caso no periodico . . . . . . 27
3.2. Comparacion de espectros de transmision entre estructuras periodicas y aquellas donde
la capa de metamaterial se incrementa, en la figura para ambos casos n=6, db=2 cm, a)
caso peri
odico da=2.25 cm; b) caso no peri
odico di=1 cm y df=4 cm . . . . . . . . . . 28
3.3. Comparacion entre los espectros de transmision entre estructuras multicapas donde en un
caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene constante
para esta figura n=8, db=2 cm; a) caso peri
odico da= 2.31 cm b) caso no peri
odico di =
1 cm , df = 4 cm

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29

3.4. Espectros de transmision para estructuras multicapas de metamateriales y dielectricos,

en un caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene
constante para esta figura n=8, db=2 cm; a) caso peri
odico da= 2.46 cm b) caso no
peri
odico di = 0.5 cm , df = 5 cm

. . . . . . . . . . . . . . . . . . . . . . . . . . . 29

3.5. Espectro de transmision para estructuras multicapas de metamateriales y dielectricos,

en un caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene
constante para esta figura n=10, db=2 cm; a) caso peri
odico da= 2.35 cm b) caso no
peri
odico di = 1 cm , df = 4 cm

. . . . . . . . . . . . . . . . . . . . . . . . . . . 30
v

INDICE DE FIGURAS

vi

3.6. Espectro de transmision para estructuras multicapas de metamateriales y dielectricos,

en un caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene
constante para eta figura se tiene n=8, db=3cm; a) caso peri
odico da=2.31cm b) caso no
peri
odico di=1cm, df=4cm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.7. Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=6, db=2
cm, a) caso peri
odico da=2.25 cm; b) caso no peri
odico di=1 cm y df=4 cm. Para esta
estructura se va a obtener un tiempo en el vacio de 0,849ns

. . . . . . . . . . . . . . 31

3.8. Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=8, db=2
cm; a) caso peri
odico da= 2.31 cm b) caso no peri
odico di = 1 cm , df = 4 cm. Para esta
estructura se va a obtener un tiempo en el vacio de 1,14ns

. . . . . . . . . . . . . . . 32

3.9. Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=8, db=2
cm; a) caso peri
odico da= 2.46 cm b) caso no peri
odico di = 0.5 cm , df = 5 cm. Para
esta estructura se va a obtener un tiempo en el vacio de 1,18ns . . . . . . . . . . . . . 32
3.10. Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=10, db=2
cm; a) caso peri
odico da= 2.35 cm b) caso no peri
odico di = 1 cm , df = 4 cm. Para esta
estructura se va a obtener un tiempo en el vacio de 1,44ns

. . . . . . . . . . . . . . . 33

Agradecimientos
A mis padres Cris
oforo y Mary Cruz por apoyarme en toda la carrera tanto moral
como econ
omicamente, a ustedes padres les dedico esta tesis que con tanto esfuerzo realice
gracias a su apoyo y
animo.
En especial a Ud. Dra. Martha Palomino quien tuvo la paciencia para orientar esta
tesis ya que todos sus comentarios fueron muy constructivos, gracias por ense
narme tanto.
Mis hermanas Jessenia y Cecilia que tuvieron que soportar todos mis estados de
animo a
un m
as cuando no sabia como realizar alg
un calculo.
A ti Ruben que est
as conmigo sin condiciones, en las buenas y en las malas, en
especial en las desveladas apoy
andome y comprendiendo mi trabajo.
A mis grandes amigos que encontre en el transcurso de la carrera y que mencionarlos
en una lista sera difcil ya que a todos los aprecio y los quiero por igual.

A todos ustedes gracias

vii

Resumen
Se ha mostrado recientemente en la literatura que la luz puede ser retardada en una
estructura de capas alternadas de dielectricos y metales, donde el ancho de la pelcula
del metal crece regularmente, y se cumple la condicion a da + b db = 0, donde da y db
representan los anchos de las pelculas del metal y del dielectrico respectivamente y a y
b las permitividades dielectricas correspondientes, esta ecuacion puede cumplirse ya que
para ciertos rangos de frecuencia en los metales la permitividad puede ser negativa, en
este trabajo se analizar
an estructuras para na da + nb db = 0, donde na y nb son los ndices
de refracci
on de los materiales, para cumplir con esta ecuacion es necesario que uno de
los ndices de refracci
on sea negativo, esto solo es posible cuando se tienen estructuras
donde uno de los componentes es un metamaterial, se estudiara la propagacion de luz en
una estructura donde se alternan capas de dielectricos y metamateriales las pelculas de
metamateriales crecen de una capa a otra, se calculara el tiempo de transito con el fin de
determinar si la luz sufre un retardo en estas estructuras y se haran las comparaciones
con el caso peri
odico.

ix

Introducci
on
Desde el principio de los tiempos la humanidad no hubiera subsistido, de no haber
tenido la habilidad de fabricar herramientas. As la civilizacion se ha desarrollado a
medida que se han dominado las propiedades de los materiales.
Al principio la luz solo era algo que diferenciaba el da de la noche, pero ahora
despues de muchos milenios ya se sabe mucho acerca de este tema. En la actualidad se
ha llegado a saber que con ayuda de la luz y una lente se pueda ver incluso los atomos.
Para poder trabajar con la luz se tuvo que dejar de lado a los materiales creados por
la naturaleza, y se tuvieron que inventar nuevos materiales: los metamateriales, estos
son materiales creados de forma artificial, en laboratorios, a partir de materiales naturales.
Cuando la luz llega a un objeto, una parte se refleja. Pero otra parte puede penetrar en
el y, al hacerlo, choca con los atomos que lo componen y se desva. Existe una magnitud
fsica llamada ndice de refracci
on, que indica cuanto se desva la luz al atravesar un material. Este ndice sale de dividir la velocidad de la luz en el vaco por la velocidad de la luz
dentro del material. En el vaco, la luz no choca con nada por lo que alcanza su velocidad
maxima, mientras que, dentro de cualquier material, la luz sera frenada por los millones
de atomos que lo componen y su velocidad sera menor que en el vaco, para la luz visible,
el ndice de refracci
on de cualquier material es siempre positivo y mayor que 1. As que
para poder controlar la luz se tendra que manipular el ndice de refracion de los materiales.
Hace aproximadamente unos 50 a
nos, el fsico ruso Vctor Veselago [1] especulo con
la posibilidad de crear materiales con ndice de refraccion negativo, lo cual en un inicio
solo result
o ser una curiosidad cientfica, pero fue hasta inicios de la decada pasada
cuando el ingles John Pendry, se le ocurrio una manera de lograrlo [2]. Los materiales de
ndice negativo (NIM por sus siglas en ingles), tambien se denominan como materiales
izquierdos (LHM por sus siglas en ingles), pueden tener simultaneamente para un mismo
rango de frecuencias la permitividad electrica y la permeabilidad magnetica negativas por
lo que se les asocia un ndice de refraccion tambien negativo. Los NIM no existen en la
naturaleza, por lo que deben ser fabricados artificialmente, y hasta ahora ha sido posible
crearlos con estas caractersticas en rangos de frecuencia de las microondas e infrarrojo [3].
Estos nuevos materiales han atrado recientemente mucha atencion debido a sus
aplicaciones sin precedentes, por ejemplo en la actualidad, las Fuerzas Armadas de
distintos pases dedican importantes esfuerzos de investigacion y recursos economicos
para conseguir reducir la detectabilidad de sus plataformas frente a los sensores utilizados
xi

xii

Introduccion

por el oponente, con el objetivo de aumentar sus probabilidades de supervivencia en

los escenarios de operaciones. Esta reduccion de detectabilidad se plasma en distintas
estrategias de reducci
on en todas sus variantes: radar, infrarroja, visible, ac
ustica, etc.
Entre las aplicaciones m
as revolucionarias y futuristas podran se
nalarse la proteccion
de las ciudades contra los efectos de los terremotos, utilizando estructuras basadas en
el concepto de metamaterial para desviar las ondas mecanicas ssmicas causadas por el
movimiento de tierras. Utilizando metamateriales se observa un efecto interesante y es
conocido como la inversion de efecto doppler tal que la frecuencia de las ondas disminuye
en lugar de aumentar al acercarse la fuente. Uno de los primeros ambitos donde se
empezo a aplicar con exito las tecnologas de metamateriales fue el de los dispositivos y
componentes que forman parte de los sistemas de RF: antenas, filtros, amplificadores, etc.
La introducci
on de las tecnicas de metamateriales en el dise
no de estos componentes ha
permitido obtener importantes mejoras en dichos dispositivos, as como desarrollar nuevos
tipos de componentes. Los sistemas de radiofrecuencia constituyen, una parte esencial de
los equipos y plataformas militares, proporcionando capacidades militares consideradas
crticas: comunicaciones t
acticas inal
ambricas, deteccion a larga distancia de objetivos y
vigilancia 24h mediante radar, guiado de proyectiles, obtencion de inteligencia, etc. En
los proximos a
nos se espera que los sistemas de RF sigan desempe
nando un papel de
primer orden en los distintos escenarios operativos, mejorando las capacidades actuales y
proporcionando nuevas capacidades (como, por ejemplo, las armas de energa dirigida de
radiofrecuencia) [4].
La diferencia principal entre los materiales y metamateriales es que en los primeros
sus propiedades dependen de su composicion, es decir, que dependen de los atomos y
moleculas que los forman. En cambio las propiedades de los metamateriales dependen
mas de su estructura que de su composicion, as que dependiendo de la distribucion de los
matateriales que conforman a un metamaterial, sera el comportamiento de la luz cuando
los atravieza.
Los metamateriales tienen propiedades opticas facinantes, ya que pueden retrasar la
luz hasta pararla. Se ha logrado por primera vez atrapar la luz de un arcoris en el rango
de frecuencias visible, de 457 a 633 nm [5]. Este nuevo dispositivo fue fabricado por los
fsicos de las universidades de Towson y Purdue, EEUU, que consta de una lente convexa
de cristal recubierta de ambos lados con una pelcula de oro. La luz es almacenada en la
capa de aire que queda en medio. En estos lugares se almacena la luz (momentaneamente),
la posicion donde se almacena la luz depende de su color. De esta forma se lograra parar
la luz de un arcoris y atraparlo en el dispositivo.
Otros dispositivos que han dado pauta a la posibilidad de detener la luz son cristales
fotonicos unidimensionales contruidos con capas alternadas de metales y dielectricos en
donde la luz puede ser atrapada gracias al cambio de signo de la constante dielectrica en
la interface metaldielectrica, donde los espesores de las losas metalicas aumentan linealmente y cumplen con la condici
on a da + b db = 0 donde da y db representan los anchos
de las pelculas del metal y del dielectrico respectivamente y a , b las permitividades
dielectricas correspondientes. La diferencia de grosores de los periodos del metal ofrecen
diferentes longitudes de onda de corte. Cuando se construye una estructura de este tipo

xiii
donde los grosores del metal aumentan linealmente, la luz de diferentes longitudes de onda
se localiza en diversas posiciones en el espacio. Esto es conocido como el arcoris atrapado
de almacenamiento de la luz [5]. Como una extension a este fenomeno resulta un objeto
interesante de estudio la construccion de un cristal fotonico unidimensional compuestos
de capas alternadas de dielectricos y metamateriales, dado que los metamateriales tienen
la caracterstica de que para ciertas regiones de frecuencia el ndice de refraccion es
negativo es posible que se cumpla con la condicion na da + nb db = 0 siendo na y nb los
ndices de refracci
on de los dos medios.
En esta tesis se pretende determinar el retraso que tiene la luz en una estructura de
metamaterial-dielectrico, esto es, se analizara bajo que circunstancias es posible atrapar la
luz. Para ello se calcular
a mediante el metodo de la fase estacionaria el tiempo de transito
de la luz, en la estructura periodica de capas alternadas de metamaterial y dielectrico
y se har
an las comparaciones con estructuras dielectrico-metalicas, por otra parte se
analizar
an los casos periodicos, con quellas estructuras donde la capa de metamaterial
crece uniformemente.
Para lograr estos objetivos la tesis se desarrolla de la siguiente manera: en el primer
captulo se describen los metamateriales, su estructura, composicion, propiedades y aplicaciones. En el segundo captulo se presenta el desarrollo teorico que lleva a establecer
las ecuaciones que determinar
an los tiempos de transito, esto es, partiendo de las ecuaciones de Maxwell se construye la matriz de transferencia que nos permite calcular los
coeficientes de transmisi
on y reflexion, necesarios para los calculos del tiempo de fase. En
el tercer captulo se presentan los resultados as como su discusion, finalmente se dan como
captulos aparte las conclusiones y la bibliografa.

Objetivos
Objetivo principal
Calcular e interpretar fsicamente el tiempo de transito de una onda electromagnetica
que atraviesa una estructura metamaterial-dielectrica, donde las capas de metamateriales
crecen uniformemente.

Objetivos especficos

1. Calcular el tiempo de tr
ansito de una onda electromanetica en una estructura de
capas alternadas de metamaterial-dielectrico.
2. Comparar los tiempos de transito de una onda electromagnetica en una estructura
de capas alternadas de metamaterial-dielectrico donde las capas de metamaterial
crecen uniformemente, con los casos donde estas estructuras son peeriodicas.

xv

Captulo 1

METAMATERIALES
1.1.
1.1.1.

Definici
on y caracteristicas de un metamaterial
Definici
on e historia

Los metamateriales son tambien conocidos como: materiales izquierdos o materiales

de ndicie de refracci
on negativo. El termino de metamaterial fue creado por Rodger M.
Walson de la Universidad de Texas en Austin en 1999[6]. El prefijo meta ( en griego)
significa m
as all
a, as que en este sentido el nombre Metamaterialessignifica los
sistemas que est
an m
as all
a de los materiales convencionales. La palabra Metamaterial
apareci
o por primera vez en la literatura en 2000 cuando Smith publico su artculo
semanal sobre un material estructurado con permeabilidad y permitividad negativa en el
rango de frecuencias de las microondas [8].
J. B. Pendry fue el primero en teorizar una forma practica de hacer un metamaterial
izquierdo (LHM). Izquierdo en este contexto significa un material en el que la regla
de la mano derecha no es obedecida, lo que permite que una onda electromagnetica
transmita energa (con una velocidad de grupo) en la direccion opuesta a su velocidad
de fase. La idea inicial de J. B. Pendry, era que una distribucion de cables metalicos
alineados a lo largo de la direccion de propagacion de la onda dan lugar a una permitividad efectiva negativa ( < 0). Sin embargo, existen materiales naturales (como
Ferroelectricos) con permitividad negativa: el reto era construir un material que tuviera
para el mismo rango de frecuencias una permeabilidad negativa ( < 0). En 1999, Pendry
demostr
o que un anillo (en C) con el eje a lo largo de la direccion de propagacion podra
proporcionar esa permeabilidad negativa. De esa manera, una distribucion periodica
de esos cables y anillos poda dar lugar a un ndice de refraccion efectivo negativo [10] [11].
Los metamateriales se definen en sentido amplio como estructuras artificiales, homogeneas desde el punto de vista electromagnetico, que exhiben propiedades fsicas que
no est
an presentes en la naturaleza. Seg
un el Instituto Europeo Virtual para Materiales
Electromagneticos Artificiales y Metamateriales son una disposici
on artificial de elementos estructurales, dise
nada para conseguir propiedades electromagneticas ventajosas e
inusuales. S
olo se pueden obtener artificialmente en laboratorios por medio de una serie
de metodos de fabricaci
on y se caracterizan por tener unas propiedades macroscopicas
1

CAP
ITULO 1. METAMATERIALES

1.1. DEFINICION Y CARACTERISTICAS DE UN METAMATERIAL

distintas a las de sus materiales constituyentes basicos, las cuales dependen de su
estructura, normalmente ordenada, y poco de su composicion [4].
Los metamateriales poseen capacidades extraordinarias para desviar las ondas
electromagneticas. Mientras que todos los materiales encontrados en la naturaleza tienen
un ndice de refracci
on positivo (el ndice de refraccion es una medida de cuanto se
desvan las ondas electromagneticas al pasar de un medio a otro), los metamateriales
presentan un ndice de refracci
on negativo, lo que implica que las ondas incidentes se
desven hacia el otro lado de la perpendicular a la superficie normal que separa los
medios. Por este motivo, tambien se les denomina como materiales izquierdos. Esto se
debe a que los metamateriales tienen simultaneamente una permitividad electrica () y
una permeabilidad magnetica () negativas, por lo que el ndice de refraccion, n, toma
un valor negativo [4].

1.1.2.

Propiedades

Un metamaterial cumple con varias propiedades electromagneticas inusuales, estas

propiedades provienen del dise
no y de la forma en la que se distribuyen y estructuran los
componentes con los que se realiza el metamaterial. Las propiedades de los metamateriales
son principalmente debido a la arquitectura celular. Se pueden dise
nar materiales con
caractersticas no observadas en los materiales que existen en la naturaleza, esto a partir
de que las celdas unitarias que definan su estructura tengan propiedades especficas de tal
manera que den una respuesta efectiva. Este es la mayor ventaja de los metamateriales[6].
Refracci
on negativa: Ley de Snell inversa. La respuesta electromagnetica de un
material se caracteriza por los par
ametros macroscopicos que provienen de una respuesta
efectiva dada por la permitividad electrica y la permeabilidad magnetica. Valores
negativos en ambos par
ametros para un determinado rango de frecuencias producen un
ndice de refracci
on negativo.
El cambio de signo en el ndice de refraccion afecta de manera directa a la ley de Snell.
Este fenomeno, supone la intervenci
on de dos medios distintos. El angulo de refraccion
(ley de Snell de refracci
on) se invierte con respecto a la situacion convencional cuando
una onda electromagnetica cruza la interfase entre dos medios que poseen signos opuestos
en sus ndices de refracci
on. El
angulo de reflexion, por otra parte, permanece inalterado,
ya que solo afecta la propagaci
on dentro de un medio.
Una expresi
on general de la ley de Snell teniendo en cuenta la existencia de medios
con ndices negativos se puede escribir como:
s1 |n1 | cos 1 = s2 |n2 | cos 2

(1.1)

en la que los subndices numericos hacen referencia a cada uno de los distintos medios, s
representa el signo del tipo de medio y los angulos de incidencia y refraccion (los ndices
numericos se asignan dependiendo de la relacion del medio con el sentido de propagacion
de la onda). Para dar una idea m
as gr
afica del comportamiento de un material izquierdo,
2

CAP
ITULO 1. METAMATERIALES

1.1. DEFINICION Y CARACTERISTICAS DE UN METAMATERIAL

se podra decir que si el agua tuviese un ndice de refraccion negativo, un pez que estuviese
nadando, parecera estar haciendolo en el aire por encima de la superficie del agua. En
la figura 1.1 vemos otro ejemplo de como se vera un objeto al introducirlo en agua que
tuviera un ndice de refracci
on negativo [4] [3].

Figura 1.1: El agua con ndice de refraccion negativo si se introduce un popote se vera doblar en direccion
contraria a la norma

Inversi
on del efecto Doppler: El efecto Doppler es el cambio aparente en la
frecuencia de una onda que percibe un observador que se encuentra en movimiento
relativo con respecto a la fuente de dicha onda. El observador distingue los cambios
de frecuencia fijando una referencia de fase definida. La frecuencia percibida por el
observador se puede expresar como w
~ = w~0 + ~k~v , siendo w0 la frecuencia real, ~k el vector
de onda y ~v la velocidad de la fuente con respecto al detector. El signo de la velocidad es
positivo cuando la fuente se aproxima al observador. En este caso, dado que el vector de
onda posee la misma direcci
on y sentido, la frecuencia percibida es mas alta. Lo contrario
sucede cuando la fuente se aleja del observador.

Sin embargo, cuando las ondas se propagan a traves de un medio que posee un ndice
de refracci
on negativo, el vector de onda cambia su signo con respecto a la situaci
on
considerada normal, alterando el resultado percibido por parte del observador. Es decir, los
vectores de onda y velocidad son antiparalelos cuando la fuente se aproxima al observador.
La figura 1.2 ilustra la diferencia existente entre ambos casos. Los metamateriales exhiben
esta propiedad de invertir el efecto Doppler como consecuencia del cambio de signo relativo
que se produce entre el vector de propagacion y el vector de onda [4].
3

CAP
ITULO 1. METAMATERIALES
1.2. APLICACIONES

Figura 1.2: Efecto Doppler en (a) un material convencional y (b) un material izquierdo
Inversi
on del efecto Cherenkov: La radiacion Cherenkov es la radiacion proveniente de la energa electromagnetica emitida por una partcula cargada que cruza un
aislante a una velocidad m
as grande que la correspondiente a la velocidad de la luz
en el medio. Esta radiaci
on electromagnetica se emite dentro de un cono, con angulos
de salida originados en la partcula cargada. El angulo del cono corresponde a la
condicion de fotones emitidos que interfieren constructivamente, y puede determinarse
mediante la relaci
on cos = 1 , donde es la velocidad de la partcula normalizada
por la velocidad de la luz en el vaco y n es el ndice de refraccion. Dado que el
angulo de salida depende del signo del ndice de refraccion, el uso de metamateriales
puede invertir el resultado. Concretamente, la radiacion Cherenkov viaja hacia atras (en
sentido contrario a la direcci
on de propagacion) en un medio o material de ndice negativo.
El fenomeno de radiaci
on Cherenkov es similar al fenomeno de fugas de radiacion
presente en algunos tipos de antenas [4].
Una caracterstica adicional importante en metamateriales es su homogeneidad a gran
escala. Este problema est
a conectado con el hecho de que las propiedades electromagneticas de los medios de comunicaci
on en general se asignan con referencia a una determinada
frecuencia de la radiaci
on electromagnetica [9].

1.2.

Aplicaciones

A tan solo 12 a
nos de haberse creado los primeros metamateriales, su campo de
investigacion y aplicaci
on ha crecido r
apidamente, con la participacion de investigadores
en diversas disciplinas como la fsica b
asica, la optica, ciencia de materiales, la mecanica y
la ingeniera electrica. No olvidando que estos surgen del campo del electromagnetismo [8].
Una de las aplicaciones m
as conocida es el uso de los metamateriales como medios
para crear la invisibilidad, este efecto se logra al refractar la luz en otra direccion lo que
tiene como consecuencia el no ser visto.

CAP
ITULO 1. METAMATERIALES
1.3. TIPOS DE METAMATERIALES
En la
optica una aplicaci
on atractiva es crear una superlente, que puede ser ampliamente utilizada en el tratamiento de imagenes medicas de super-resolucion, formaci
on
de im
agenes
opticas. La primer superlente fue creada en el 2004 para el regimen de las
microondas, y mostr
o tener una resolucion tres veces mejor que el lmite de difracci
on.
Mas tarde, Fang propuso las primeras superlentes opticas usando una pelcula delgada
de plata lo que rompe el lmite de difraccion y de esta manera es posible producir una
super-resoluci
on de im
agenes [6] [7].
Los metamateriales prometen no solo mejorar las prestaciones de un gran n
umero de
dispositivos y sistemas en los campos de la electronica y de la optoelectronica, sino introducir cambios verdaderamente revolucionarios en un gran n
umero de sectores relacionados
con estas tecnologas. De esta manera los metamateriales se dise
nan teoricamente para
controlar el camino de propagacion de todo tipo de ondas (principalmente electromagneticas, pero tambien ac
usticas, ssmicas, etc.). Los metamateriales ya estan siendo utilizados
en telefonos m
oviles como el representado en la figura 1.3. Otras aplicaciones como la invisibilidad se encuentran todava lejos del alcance de la tecnologa actual, siendo necesarios
grandes avances en el
area de la nanotecnologa para conseguir fabricar metamateriales
con las caractersticas requeridas. En la actualidad, la investigacion sobre los efectos de la
invisibilidad u ocultamiento de objetos con metamateriales es extraordinariamente intensa
debido a las posibles aplicaciones de estos [4].

Figura 1.3: El telefono movil de la figura (a) (modelo BL40 de la marca LG) fue el primer dispositivo de
este tipo en integrar una antena basada en metamateriales (modelo MTM-C de la empresa norteamericana
Rayspan) (b). Esta tecnologa permite reducir a una quinta parte el tama
no de antena que sera necesario
en caso de utilizar tecnologa convencional de antenas impresas.

1.3.
1.3.1.

Tipos de metamateriales
Diel
ectricos artificiales

Los dielectricos artificiales son estructuras heterogeneas construidas a una escala de

una determinada longitud de onda de tal manera que se pueda controlar la permitividad
efectiva de la macroestructura. Winston E. Kock demostro experimentalmente que un
arreglo de placas met
alicas planas, con periodicidad espacial p 2 , podra formar el
patr
on del haz de una antena en el mismo sentido que una lente dielectrica homogenea
5

CAP
ITULO 1. METAMATERIALES
1.3. TIPOS DE METAMATERIALES
cuyo ndice de refracci
on n sea menor que la unidad: 0 n 1
s
 2

1
2p

(1.2)

Este resultado se obtiene teniendo en cuenta que cuando un campo electrico E es

paralelo a las placas, el modo fundamental de la gua de ondas paralela a las placas tiene
una velocidad de fase mayor que la velocidad de la luz en el vaco. Desde que se trabaja con
dispositivos en el rango de longitudes de onda milimetricas y centimetricas, la investigacion
de dielectricos artificiales en las microondas y frecuencias de radio ha aumentado [4].

1.3.2.

Materiales magn
eticos artificiales

La permitividad electrica de un dispositivo se puede controlar mediante la eleccion

del material en funci
on de la frecuencia de excitacion . Por ejemplo, los cristales de
silicio exhiben permitividades positivas a frecuencias opticas y se utilizan en una amplia
variedad de dispositivos, como microscopios opticos, telescopios o simplemente lentes de
gafas comerciales. Los metales nobles, que se comportan como gases de electrones libres,
tienen una permitividad negativa a frecuencias de excitacion debajo de su frecuencia
intrnseca electrica de plasma pe . Para obtener permitividad negativa sin que se tengan
perdidas de transmisi
on de potencia elevadas a traves de metales opacos, varios grupos
de investigaci
on utilizan arreglos de hlos finos.
Mientras que los compuestos electricamente sensibles son relativamente faciles de encontrar, los materiales magneticos con permeabilidad magnetica que se encuentran en
la naturaleza son mucho m
as raros y consisten esencialmente en ferritas (materiales ferromagneticos) excitados a una frecuencia de plasma magnetica pm . Ademas, no se conocen
materiales con propiedades magneticas en frecuencias de terahertz o mayores. Por este
motivo, obtener materiales con la permeabilidad magnetica a altas frecuencias es muy
complicado y requiere estructuras dise
nadas especficamente para generar una respuesta
magnetica. En el libro Magnetism from conductors and enhanced nonlinear phenomena
de John Pendry sugiri
o el uso de resonadores de anillos abiertos (Split Rings Resonators SRRs), b
asicamente de tama
no de longitudes de onda electomagneticas del orden
de milmetros, para generar una permeabilidad magnetica efectiva. En estas estructuras,
cuando el campo electromagnetico incidente tiene una polarizacion perpendicular al plano
del arreglo de anillos, se inducen corrientes resonantes en lazo que emulan a un dipolo
magnetico [2] [4].

1.3.3.

Materiales quirales

En los materiales quirales la geometra de la estructura interna produce una respuesta

macroscopica y se observa la rotaci
on de la polarizacion del vector de propagacion. Los
materiales con una importante capacidad de rotacion o dicrosmo circular se utilizan para
la fabricacion de dispositivos optoelectronicos, como polarizadores circulares en aplicaciones como la fotografa. Las estructuras giroscopicas en los materiales quirales se estan
investigando como otra alternativa para obtener materiales NRI (negative refractive index). La quiralidad se obtiene mediante estructuras helicoidales y geometras retorcidas.
6

CAP
ITULO 1. METAMATERIALES
1.4. MODELO DE GUPTA
Se han demostrado peque
nos valores NRI en las microondas mediante materiales quirales
basados en estructuras helicoidales. Tambien se ha realizado el mismo experimento con
estructuras m
as peque
nas en el visible y el infrarrojo, y en estas longitudes de onda se
ha demostrado una significativa polarizacion circular pero no se ha obtenido refracci
on
negativa [4].

1.3.4.

Materiales con recubrimiento

Los materiales con recubrimiento se investigaron en las decadas de los 80 y de los 90 y

se concluy
o te
oricamente que los objetos invisibles no pueden existir. No se consideraron
otras opciones hasta que aparecieron dos artculos en la revista Science en 2006, lo que
implic
o que se reabrieran las investigaciones. Ambos artculos explicaban de una manera
similar c
omo hacer un objeto invisible a la luz incidente: el objeto debe tener un ndice de
refracci
on con un gradiente de distribucion de tal manera que todos los rayos incidentes
sean guiados alrededor de su superficie y recuperen su trazado original en el otro extremo.
La figura 1.4 muestra un objeto al que se le ha puesto un recubrimiento de metamaterial
que produce una desviaci
on de la luz te tal suerte que el objeto en el centro no se le puede
ver [4].

Figura 1.4: Rayos de luz atravesando objetos invisibles.

1.4.

Modelo de Gupta

Las propiedades inusuales de las ondas electromagneticas en un medio capaz de poseer

simult
aneamente constante dielectrica negativa  y la permeabilidad magnetica negativa
fueron descubiertos por Veselago a finales de los sesenta. Estos medios son referidos
como materiales de mano izquierda LH.
A continuaci
on se muestra el modelo de Gupta que esta descrito por las siguientes
ecuaciones, este modelo se deriva de los modelos propuestos por Pendry .
( 2 p2 )
() =
,
2

() =

( 2 b2 )
( 2 02 )

(1.3)

donde p = 212GHz, b = 26GHz, 0 = 24GHz. El ndice de refraccion sera positivo

para > p y negativo para 0 < < b . En otros rangos de , el ndice de refracci
on
7

CAP
ITULO 1. METAMATERIALES
1.4. MODELO DE GUPTA
sera puramente imaginario y las ondas correspondientes representan modos evanescentes.
La figura 1.5 muestra la permeabilidad magnetica y la permitividad electrica de acuerdo
al modelo usado por Gupta, se observa el peque
no rango de frecuencias en donde ambas
funciones son negativas [12] [13] [14].

Figura 1.5: Grafica de () y () modelo de Gupta .

Captulo 2

DESARROLLO TEORICO
2.1.

Ecuaciones de Maxwell

La Teora Electromagnetica del fsico escoces James Clerk Maxwell (1831-1879) es

una de las obras intelectuales mas importantes en la historia de las ciencias. Su aparici
on
se inicia en 1861 (On Physical Lines of Force) y se completa en un tercer trabajo en
1865 (A Dynamical Theory of the Electromagnetic Field ). Para esta epoca ya se conocan
varias leyes individuales sobre el comportamiento de la electricidad y el magnetismo, pero
no se tena una teora formal que explicara los fenomenos de naturaleza electromagnetica.
Maxwell supo seleccionar cuatro fenomenos basicos fundamentales como principios,
con los cuales arm
o un modelo fsico matematico capaz de explicar la totalidad de las
leyes en esa disciplina y predecir fenomenos desconocidos. Esta teora es considerada
el nacimiento de la Fsica Moderna y las consecuencias incidieron en todas las ramas
de la fsica, adem
as de conformar un modelo completo para los fenomenos clasicos del
electromagnetismo, explic
o la optica ondulatoria y, en parte, la naturaleza de la luz.
Predijo la existencia de ondas electromagneticas y demostro que el campo es un ente
fsico real e independiente de la materia.
Cabe se
nalar que la Teora de Relatividad Especial esta implcita en las ecuaciones
de Maxwell pues ellas se cumplen con rigor en todos los sistemas inerciales, lo que
permite deducir naturalmente las Transformaciones de Lorentz como relaciones u
nicas de
transformaci
on de coordenadas entre sistemas inerciales.
La formulaci
on moderna del electromagnetismo fue elaborada en 1884 por el cientfico
Olivier Heaviside (1850-1925), para lo cual estructuro el analisis vectorial y replante
o la
formulaci
on de Maxwell, llev
andola a la forma que trata la bibliografa actual mediante
ecuaciones diferenciales en derivadas parciales.

CAP
ITULO 2. DESARROLLO TEORICO
2.1. ECUACIONES DE MAXWELL

2.1.1.

Fundamentos de las ecuaciones de Maxwell

Los cuatro fen

omenos b
asicos tomados como postulados del electromagnetismo son:
1. Ley de Faraday sobre la fuerza electromotriz inducida
Esta ley fue descubierta por Michael Faraday en 1831, quien se desempe
naba como
encargado del pa
nol del laboratorio (ordenanza) de la Royal Institution de Inglaterra,
usando un dise
no muy simple, como muestra la figura:

Figura 2.1: Dispositivo de Faraday

Al mover el im
an dentro del carton, que tena enrollado un alambre de cobre, las
laminas met
alicas del electroscopio se abran, indicando la acumulacion de cargas
electricas en ambas hojuelas como consecuencia de una corriente electrica por el
alambre de cobre, simult
anea con el movimiento. Esto indica que en el conductor de
cobre existe un campo electrico, condicion que solo se cumple cuando hay movimiento
relativo entre el im
an y el conductor.
De esta manera Faraday descubrio que la electricidad y el magnetismo se relacionaban funcionalmente si los campos eran variables en el tiempo. La forma matematica
de la ley de Faraday es:
I

~ d~s = d
E
dt
C

Z Z

~ d~a
B

(2.1)

El primer miembro (circulaci

on del campo electrico) es la definicion de la denominada fuerza electromotriz inducida en el conductor, siendo C la curva definida por
el alambre de cobre. El segundo miembro es la variacion temporal (debida al movimiento del im
an) del flujo magnetico a traves de la superficie que tiene por borde a
la curva C.
2. Ley de Gauss-Faraday sobre inducci
on el
ectrica
Los experimentos de inducci
on electrica realizados por Faraday (antes del a
no 1831)
mostraron que si una carga Q es encerrada por un recipiente conductor inicialmente
neutro, pero sin establecer contacto directo con el cuerpo cargado (ver figura 2.2), el
recipiente conductor reordena sus cargas (fenomeno de induccion) de tal manera que
las superficies interior y exterior del recipiente quedan cargadas con signo opuesto.
10

CAP
ITULO 2. DESARROLLO TEORICO
2.1. ECUACIONES DE MAXWELL
La carga total inducida en cada superficie resulta de magnitud exactamente igual a
la de la carga encerrada.

Figura 2.2: Fenomeno de induccion

Los variados experimentos de Faraday sobre induccion permitieron comprender que

los medios conductores poseen una cantidad inmensa de cargas libres en su interior
que pueden reordenarse, y mostraron que la carga neta de un conductor permanece constante ante fen
omenos inductivos, confirmando la conservacion de la carga.
As como tambien se verifico que para cuerpos en reposo el interior de los conductores
es neutro, sin campo electrico, a
un en presencia de cuerpos externos cargados.
La expresi
on matem
atica de esta ley fue dada por Gauss y reformulada por Heaviside
~ que fuera definido
con la actual forma vectorial, utilizando el campo de inducci
on D,
y medido por Faraday, cuyo modulo en un punto cualquiera del espacio representa la
densidad de carga inducida maxima que podra obtenerse si se pusiera una plaquita
met
alica (transversal al campo).

~ d~a =
D

Z Z Z
dV

(2.2)

~ a traves de cualquier superficie que

El primer miembro es el flujo del campo D
encierre la carga Q, mientras que el segundo miembro representa la carga total
encerrada.
3. La ley de Ampere
Hasta el a
no 1820 se pensaba que la electricidad y el magnetismo eran fenomenos no
relacionados. En una conferencia que daba el dinamarques Oersted (para conseguir
fondos para sus proyectos), justamente mientras intentaba mostrar dicha independencia, pos
o una br
ujula sobre un conductor con corriente provocando que la aguja
se orientara de manera transversal al conductor. As, de casualidad, descubrio que
una corriente electrica esta rodeada por un campo magnetico (ver figura 2.3).
11

CAP
ITULO 2. DESARROLLO TEORICO
2.1. ECUACIONES DE MAXWELL

Figura 2.3: Ley de Ampere

Luego, Oersted repiti
o el experimento ante sus alumnos y, aunque no logro dar una
explicaci
on satisfactoria, lo publico. Fue el fsico matematico frances A. Ampere
(1775-1836) quien interpret
o y dio la expresion matematica del fenomeno (que lleva
su nombre), adem
as de proponer a las corrientes como u
nica causa del magnetismo,
propuesta conocida como la Hip
otesis de Ampere.
Hoy se sabe que las corrientes electricas y el campo magnetico asociado no son
causa y efecto ya que ambos, corriente y campo, aparecen simultaneamente con el
movimiento de cargas. Matem
aticamente la ley de Ampere se expresa:
I
~ d~s = I
H
(2.3)
C

~ siendo C cualquier curva cerrada que

El primer miembro es la circulacion de H,
rodee a la corriente I. Esta ley es valida solo para corrientes constantes. La ley de
Ampere puede ser expresada usando el vector densidad de corriente, cuya relacion
con la corriente est
a dada por:
Z Z
I=

J~ d~a

(2.4)

Siendo S la secci
on del conductor donde circula la corriente. Dado que el contorno C
de la ley de Ampere encierra la corriente y que fuera del conductor el vector J~ es nulo
(notar que J~ no depende del tiempo), podemos extender el recinto de integracion
hasta el borde C, quedando:
I
Z Z
~
H d~s =
J~ d~a
(2.5)
C

12

CAP
ITULO 2. DESARROLLO TEORICO
2.1. ECUACIONES DE MAXWELL
4. No existencia de monopolos magn
eticos
La experiencia mostr
o que no existen polos magneticos aislados. Si un iman se parte
al medio se obtienen dos imanes de menor intensidad. Esto muestra una particu~ cuyas lneas de fuerza son necesariamente
lar propiedad del campo magnetico (B),
cerradas pues no tienen ni fuentes ni sumideros.

~
Figura 2.4: Lneas de fuerza de B

2.1.2.

Obtenci
on de las ecuaciones de Maxwell

Dado que las expresiones anteriores estan dadas por integrales y no son validas en un
punto ya que describen un fen
omeno extenso, por lo cual no siempre es posible encontrar
una relaci
on funcional v
alida punto a punto entre las magnitudes que intervienen en una
ecuaci
on integral.
As que para pasar estas ecuaciones a su forma diferencial se hace uso del procedimiento
que realizo Heaviside utilizando dos teoremas del analisis vectorial:
Teorema de Gauss
I

~ d~a =
A

Z Z Z

~ dV
A

(2.6)

~ d~a
A

(2.7)

Teorema de Stokes
I

~ d~s =
A

Z Z

As que partiendo de la Ley de Faraday sobre la fuerza electromotriz inducida, usando

el Teorema de Stokes y realizando los calculos apropiados se obtiene la primera ecuaci
on
de Maxwell.
~
~ = B
E
t

(2.8)

~ y B,
~ con sus
Esta ecuaci
on muestra que en un punto cualquiera pueden coexistir E
formas funcionales relacionadas por la ecuacion dada.

13

CAP
ITULO 2. DESARROLLO TEORICO
2.1. ECUACIONES DE MAXWELL
Ahora partiendo de la ley de Gauss-Faraday sobre induccion electrica y usando el
teorema de Gauss se obtine la segunda ecuacion de Maxwell.
~ =
D

(2.9)

~ son las cargas positivas y los

Esta ley escalar nos indica que las fuentes del campo D
sumideros las cargas negativas. El campo electrico asociado a una carga nace en ella (si
es positiva) o muere en ella (si es negativa).
Para obtener la tercera ecuaci
on de Maxwell se parte de la ley de Ampere y se usa el
teorema de Stokes y la segunda ecuaci
on de Maxwell obteniendo:
~
~ = J~ + D
H
(2.10)
t
Finalmente para obtener la cuarta ecuacion de Maxwell se acepta que las lneas de
fuerza del campo magnetico son cerradas, hecho verificado experimentalmente, la expresion
~ no tiene fuentes ni sumideros. En
matematica es inmediata pues el campo magnetico B
consecuencia, su divergencia es nula.
~ =0
B

(2.11)

En resumen las cuatro ecuaciones de Maxwell, descritas por Heaviside, son consideradas
los Principios de la Teora Electromagnetica, que corresponden a cuatro fenomenos basicos
que no tienen demostraci
on te
orica. Es importante recalcar que de estas ecuaciones se
deducen todas las leyes conocidas del electromagnetismo, conformando una teora clasica
completa [15][16].
Ellas son:
~
~ = B
E
t
~
D =

(2.12)

~
~ = J~ + D
H
t
~
B =0
donde:
E: Campo electrico existente en el espacio, debido a las cargas
B: Campo de inducci
on magnetica, debido a las corrientes
H: Campo magnetico que resume los efectos magneticos de la materia
D: Campo de desplazamiento que resume los efectos electricos de la materia
J: Densidad de corriente electrica
: Densidad de carga electrica
con
H=

B
,

D = E
14

(2.13)

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

y
: permeabilidad magnetica
: permitividad electrica

2.2.

M
etodo de la matriz de transferencia

La forma en que una onda se propaga al atravesar una interfase obedece a la ley de
refracci
on, y el c
alculo de la transmision y la reflexion de dicha onda se resuelve aplicando
las Ecuaciones de Maxwell y las condiciones de frontera en forma adecuada, dado que
se va a estudiar una superred compuesta de dos tipos de materiales, para calcular los
coeficientes de transmisi
on
o reflexion, se va hacer uso de el metodo de la matriz de
transferencia (MT Matriz de Transferencia), se hara el desarrollo para dos casos, cuando la
luz incide sobre las interfaces con polarizacion Transversal Magnetica (TM o polarizaci
on
P) y cuando la radiaci
on incide con polarizacion transversal electrica (TE o polarizacion S).
El Metodo de MT permite relacionar el campo electrico (magnetico) a lo largo de toda
la estructura, donde esta se construye repitiendo una celda unitaria n veces, el Metodo de
MT se realizar
a para dos casos especiales de polarizacion conocidas como Polarizacion P
y Polarizaci
on S para hallar los coeficientes de transmision y reflexion.

2.2.1.

Ecuaci
on de onda

Partiendo de la Ecuaciones de Maxwell para un medio isotropico D=E, con permeabilidad =1 medio no magnetico y suponiendo que no esisten cargas o correintes libres
externas en la regi
on de interes es decir que se toma =0 y J=0. As que en este caso las
Ecucaciones de Maxwell (Ecuaciones 2.12) nos quedan como sigue:
~ =0
E

(2.14)

~ =0
B

(2.15)

~
~ + 1 B = 0
E
c t

(2.16)

~
~ =  E
(2.17)
B
c t
Combinando las cuatro ecuaciones mencionadas anteriormente se va a obtener la ecuancion de onda. Ahora aplicando el rotacional a la ecuacion 2.16 se obtiene que.
~
~ + 1 B = 0
( E)
c
t
~ (B
~ C)
~ = B(
~ A
~ C)
~ C(
~ A
~ B)
~ entonces se tiene que:
recordando que A
~ = ( E)
~ ( )E
~ = 2 E
~
( E)
15

(2.18)

(2.19)

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

reemplazando la ultima parte de la igualdad anterior en la ecuacion (2.18) y sustituyendo
(2.17) se obtiene:

~
2 E

se define v 2 =

~
 2 E
=0
2
2
c t

(2.20)

c2
para obtener finalmente la ecuacion de onda

~
2 E

~
1 2E
=0
v 2 t2

(2.21)

Las soluciones a esta ecuaci

on son de la forma

~ = E~0 e(~k~rt)
E
donde

k2 =

(2.22)

2

c2

Esta es la ecuaci
on de onda para el campo electrico en un medio lineal, homogeneo e
isotropico. An
alogamente se obtiene la ecuacion de onda para el campo magnetico
aplicando el rotacional a la ecuaci
on (2.17)

~
2 B

~
1 2B
=0
2
2
v t

(2.23)

con soluciones de la forma

~ = B~0 e(~k~rt)
B

2.2.2.

(2.24)

Polarizaci
on P (Polarizaci
on Transversal Magn
etica)

Se considera la propagaci
on de ondas en una superred creciendo en en direccion z y las
capas son paralelas al plano xy, con periodicidad en una dimension formada por un medio
dielectrico homogeneo alternandose con metamaterial que va creciendo exponencialmente
a una unidad, libre de cargas y corriente donde no hay carga ni corrientes sobre la superficie
( = 0 y J = 0), con incidencia normal, como se muestra acontinuacion.[17][18]
16

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

Figura 2.5: Polarizacion P

de la figura es facil ver que:
~k = (kx , 0, kz )
~ = (Ex , 0, Ez )
E

(2.25)

~ = (0, By , 0)
B

(2.27)

(2.26)

Las ecuaciones del campo eletrico y del campo magnetico quedan de la siguiente forma:
~ = Ex ei(~k~rt) + Ez ei(~k~rt) k
E

(2.28)

~ = Bei(~k~rt)
B
para relacionar las componentes del campo electrico E se tiene que
~ = ikx Ex ei(kx x+kz zt) + ikz Ez ei(kx x+kz zt) = Ex eiA + Ez eiA = 0
E

(2.29)

Ahora se calcula la amplitud del campo electrico E

~ = ikx Ex ei(kx x+kz zt) + ikz Ez ei(kx x+kz zt) = 0
E

(2.30)

iei(kx +kz zt) [kx Ex + kz Ez ] = 0

(2.31)

entonces
as que se obtine
kx
Ex
(2.32)
kz
Se calcula la amplitud del campo magnetico B usando. Aplicando el rotacional se tiene





~ =
 =
[ikx Ez eiA ikz Ex eiA ]
E
(2.33)
y
z
 x

 E eiA 0 E eiA 
x
z
Ez =

ahora realizando

~
1 B
i
~
= By ei(k~rt)
c t
c
17

(2.34)

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

igualando las ecuaciones (2.33) y (2.34)
kz Ex kx Ez =

By
c

(2.35)

sustituyendo la ecuaci
on (2.32) se obtiene que

Ex (
Recordando que

kz2 + kx2

) = By
kz
c

k 2 = kx2 + kz2 =

 2
c2

By =

(2.36)

donde = 1, entonces

Ex
ckz

(2.37)


ckz

(2.38)

Se define la impedancia como

Y =
entonces

By = Y E x

(2.39)

las condiciones de frontera que satisfacen los campos en la interfase son:

Ex+ = Ex

Bx+ = Bx

(2.40)

Figura 2.6: Onda incidiendo en la red

Las componentes de las ondas incidentes son :
Exj = A1j eikzj zj A2j eikzj zj
Byj = Yj Exj = Yj (A1j eikzj zj + A2j eikzj zj )
Se puede escribir de forma matricial
18

(2.41)

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

Exj
Byj


=
zj

1
Yj

1
Yj



Exj
Y Exj

1
Yj

=
zj



1
Yj

A1j eikzj zj
A2j eikzj zj

= Bj

Llamamos Bj a la matriz

eikzj (zj zj )

1
Y

1
Y

eikzj (zj zj )


(2.42)
zj
0

A1j eikzj zj

!
0

A2j eikzj zj

zj

y, a d = zj zj .
zj

se define

T (d) =

eikzj d
0
ik
0
e zj d


(2.43)

entonces
0

Exj
Y Exj


= Bj T (d)
zj

A1j eikzj zj

!
(2.44)

A2j eikzj zj

zj

Por otro lado para la frontera en zj que es la j-esima capa se tiene que:
0

Exj
Y Exj

A1j eikzj zj


0

= Bj

!
(2.45)

A2j eikzj zj

zj

zj

despejando de la ecuaci
on (2.45) y sustituyendo en (2.44) se obtiene que.




Exj
Exj
1
= Bj T (d)Bj
Y Exj z 0
Y Exj z
j

(2.46)

Ahora se define la matriz mj = Bj T (d)Bj1 donde mj es llamada la matriz de

transferencia en cada capa. Despues de realizar las operaciones correspondientes
calculando Bj1 se encuentra que:

mj =

cos(kzj d)
iYj sin kzj d

i
Yj

sin kzj d
cos kzj d

!
(2.47)

Esta matriz relaciona los campos electrico y magnetico de la onda incidente entre capas,
finalmente la matriz que relaciona la amplitud de la onda incidente con las ondas reflejada
y transmitida M resulta del producto de todas las matrices mj de cada una de las capas
por lo que se tiene M = mN mN 1 ...mj ...m1 ; esto es:


ExN
Y ExN


=M

zn

19

Ex0
Y Ex0


(2.48)
z=0

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

Para obtener los coeficientes de reflexi
on y transmision se considera que la onda incidente
esta normalizada, es decir que A1 = 1 y A2 = r para z=0, de las ecuaciones 2.41 obtenemos:

Ex0 (0) = 1 r

(2.49)

By0 (0) = Y0 Ex = Y0 (1 + r)

(2.50)

En z = L, donde L es la longitud total de la estructura A1 = t y A2 = 0

ExN +1 (L) = teikzN +1 L

(2.51)

ByN +1 (L) = YN +1 teikzN +1 L

(2.52)

Sustituyento en la ecuaci
on 2.48 se tiene la relacion matricial :

teikN +1 L
tYN +1 eikN +1 L


=M

1r
Y0 (1 + r)


(2.53)

De donde se obtiene el siguiente sistema de ecuaciones:

teikN +1 L = m11 (1 r) + m12 Y0 (1 + r)

tYN +1 e

ikN +1 L

= m21 (1 r) + m22 Y0 (1 + r)

(2.54)
(2.55)

Cuyas soluciones son:

t=

(m11 m12 Y0 )(m21 + m22 Y0 ) (m11 + m12 Y0 )(m21 m22 Y0 )

eikN +1 L [(m11 m12 Y0 )YN +1 (m21 m22 Y0 )]

(2.56)

(m11 + m12 Y0 )YN +1 (m21 + m22 Y0 )

(m11 m12 Y0 )YN +1 (m21 m22 Y0 )

(2.57)

r=

Donde t y r son los coeficientes de reflexion y transmision respectivamente.

2.2.3.

Polarizacion S (Polarizaci
on Transversal El
ectrica)

Se analiza ahora la propagaci

on de la luz sobre el mismo medio con polarizacion transversal electrica.[17][18]
20

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

Figura 2.7: Polarizacion S

Ahora se considera la polarizacion transversal electrica, donde los vectores estan dados
por:
~ = (0, Ey , 0)
E
~ = (Bx , 0, Bz )
B

(2.59)

~k = (kx , 0, kz )

(2.60)

k = kx2 + kz2

(2.61)

(2.58)

se observa que

Las ecuaciones del campo electrico y magnetico estan dadas por:

~ = Ey ei(~k~rt)
E

(2.62)

~ = Bx ei(~k~rt) + Bz ei(~k~rt) k
B

(2.63)

Analogamente realizando las operaciones que se mostraron en la seccion 2.2 se encuentra

que:

Bz =

Bx kx
kz

(2.64)

~ y realizando operaciones similares y mediante el uso de la

Aplicando el rotacional a B
ecuaci
on (2.17) se llega a :

Bx = Y E y
zc
Donde Y = k
es la impedancia.
Una onda incidente se tienen las siguientes componentes tangenciales:

21

(2.65)

CAP
ITULO 2. DESARROLLO TEORICO

2.2. METODO DE LA MATRIZ DE TRANSFERENCIA

Ey = A1 eikz z A2 eikz z
ikz z

Bx = Y (A1 e

ikz z

A2 e

(2.66)

(2.67)

Expresado en forma matricial :



Ey
Bx


=

zj

A1 eikz z
A2 eikz z



1 1
Yj Yj


= Bj
zj

A1 eikz z
A2 eikz z


(2.68)
zj

(2.69)
con

Bj =

1 1
Yj Yj

que puede escribirse como :

0

Eyj
Bxj


= Bj

eikz (zj zj )

eikz (zj zj )

A1 eikz zj

zj

!
(2.70)

A2 eikz zj

z0 j

sea
0

d = zj z j y

T (d) =

eikz d
0
0
eikz d

Para la frontera en la j-esima capa se tiene:



Eyj
Bxj


= Bj
z0 j

A1 eikzj zj
0
A2 eikzj zj

(2.71)
z0 j

De donde obtenemos:


Eyj
Bxj

Bj T (d)Bj1

zj

Eyj
Bxj


(2.72)
z0 j

se encuentra que la matriz de transferencia es :

mj =

cos(kz d)
Yij sin kz d
iYj sin kz d
cos kz d

!
(2.73)

Para el arreglo completo de N capas se tiene la siguiente relacion:



Eyn
Bxn


=M

zN

22

Ey1
Bx1


(2.74)
z0

CAP
ITULO 2. DESARROLLO TEORICO

2.3. CALCULO
DE TIEMPO DE FASE
Siguiendo la hip
otesis de que la onda incidente esta normalizada para z = 0 tenemos que
A1 = 1 y A2 = r, entonces para la onda incidente y reflejada :
Ey (0) = 1 r

(2.75)

Bx (0) = Y0 (1 r)

(2.76)

Para la onda trasmitida z = n + 1, y A1 = t y A2 = 0

Ey (L) = teikN +1 L

(2.77)

ikN +1 L

(2.78)

Bx (L) = YN +1 te
Sustituyendo en la ecuaci
on (2.74) se obtiene:


teikN +1 L
Yn+1 teikN +1 L


=m

1 r
Y0 (1 r)


(2.79)

Se obtiene el siguiente sistema de ecuaciones

teikN +1 L = m11 (1 + r) + m12 Y0 (1 + r)

(2.80)

YN +1 teikN +1 L = m21 (1 + r) + m22 Y0 (1 + r)

resolviendo el sistema de ecuaciones se llega a las soluciones t y r coeficientes de transmisi
on
y reflexi
on respectivamente:

t=

(m11 + m12 Y0 )(m22 Y0 m21 ) (m12 Y0 m11 )(m21 Y0 + m22 )

(m11 + m12 Y0 )YN +1 eikN +1 L + (m21 + m22 )eikN +1 L

(2.81)

(m12 Y0 m11 )YN +1 + (m22 Y0 m21 )

(m12 Y0 + m11 )YN +1 (m22 Y0 + m21 )

(2.82)

r=

2.3.

C
alculo de tiempo de fase

El tiempo fase es el tiempo que tarda el paquete de onda en atravesar una region del
espacio. Este se obtiene a partir de un analisis de la fase de la funcion que describe un
paquete de ondas incidente.
En varios casos se representa a una partcula a travez de un paquete de onda, y el
tiempo de fase de los paquetes incidente, reflejado y transmitido se les asocia como los
tiempos de incidencia, reflexi
on y transmision.
Los paquetes de onda se obtienen dependiendo de la energa E que tenga la particula,
parte del paquete se refleja y otra parte se penetra en la region y emerge en el otro lado
como otro paquete de onda. Como condicion inicial se pone una partcula incidente a la
23

CAP
ITULO 2. DESARROLLO TEORICO

2.3. CALCULO
DE TIEMPO DE FASE
izquierda sobre el eje z de la estrucura z = l y este se propaga hacia la derecha con
velocidad de grupo (
v ). Para representar a esta partcula, desde el punto de vista de la
mecanica cu
antica, se considera un paquete de onda con distribucion de probabilidad de
impulso |(k)|2 , as como una distribucion de probabilidad espacial |(z)|2 , Alrededor de
m
v
l
la media de valores k =
. El paquete de onda llegara a z = 0 en t0 = .
}
v
Una vez que el paquete de onda inicial (z) = (z, 0) es definido, se puede obtener la
distribucion de impulso mediante una transformada de fourier.
k0

1
(z) = (z, 0) =
2

dk(k)eikz

(2.83)

de modo que
1
(k) =
2

eikz (z)dz

(2.84)

Para fines pr
acticos se considera inicialmente un paquete de onda gaussiano centrado en
z = l

(z) =

2
z02

1
4

eik(z+l) e

(x+l)2
(z0 )2

(2.85)

que corresponde a la distribuci

on de impulso centrado en k
1

(k) = g k k eikl =

ikl

2k2
con

1 e e

(kk)2
4 2
k

(2.86)

k2 = hk 2 i hki2 y hki = k

La ecuacion (2.85) representa un paquete de onda gaussiano centrado en z = l a t = 0

y la ecuacion (2.86) representa la funcion de distribucion del paquete de onda incidente
en el espacio de los vectores de onda, mientras que para los paquetes de onda reflejado y
transmitido son las siguientes :

r (k) = |r(k)|

(2.87)

t (k) = |t(k)|

(2.88)

dada la funci
on
1

g k k =

1 e
2

2k

(kk)2
4 2
k

(2.89)

La funcion de distribuci
on de paquete de onda incidente, reflejado y transmitido espacial
son:
24

CAP
ITULO 2. DESARROLLO TEORICO

2.3. CALCULO
DE TIEMPO DE FASE

g k k ei[k(z+l)wt] dk

(2.90)

rg k k ei[k(z+l)+wt] dk

(2.91)

tg k k ei[k(z+l)wt] dk

(2.92)

(z, t) =

Z
r (z, t) =

t (z, t) =

donde r y t son los coeficientes de reflexion y transmision.

Los paquetes de onda incidente, reflejado y transmitido son.
2
1
(z, t) = a |(z, t)|
2


2
1
r (z, t) = a |r (z, t)|
2
2

1
t (z, t) = a |t (z, t)|
2

(2.93)

(2.94)

(2.95)

donde a es el periodo de la estructura.

El tiempo de fase se calcula con el metodo de fase estacionaria propuesto por Zohta
[19] [20] . El tiempo de fase se obtiene aplicando la ecuacion de Zohta,
=

Im(T )
+
tan1
v
Re(T )

(2.96)

donde T es la amplitud de transmision asociada a la onda , d es el ancho de la estructura y

v es la velocidad de grupo asociada al paquete la cual aproximamos con la relacion v = /k

Al termino
td =

Im(T )
tan1

Re(T )

Zohta lo define como tiempo de retardo.[13] [21] [22] [23] [24]

25

(2.97)

Captulo 3

RESULTADOS
Los resultados que se van a mostrar a continuacion seran obtenidos utilizando un
sistema de multicapas donde la celda unitaria esta formada por una capa de metamaterial
que est
a descrita con el modelo de Gupta y una capa de aire como dielectrico, al
construir la superred el ancho del dielectrico se mantiene constante, en tanto que el
ancho de la capa del metamaterial se va incrementando en una unidad de tal manera
que la n-esima capa tiene un ancho de espesor b + n (ver figura 3.1), el medio de
incidencia y transmisi
on es el aire. Un sistema como este hace que se retenga la luz
para algunas frecuencias, para ello se calculara la transmision de la luz a traves de la
estructura mediante el metodo de la matriz de transferencia descrito en el captulo 2, y
se calcular
a el tiempo de tr
ansito mediante el metodo de fase estacionaria descrito por
el modelo de Zhota. Estos resultados se compararan con una estructura periodica donde se alternan metamaterial y dielectrico pero el metamaterial se va a mantener constante.

Figura 3.1: Forma que van a tener las multicapas a) caso periodico y b) caso no periodico

27

CAP
ITULO 3. RESULTADOS

3.1. ESPECTROS DE TRANSMISION

3.1.

Espectros de transmisi
on

En esta secci
on se muestran espectros de transmision para estructuras descritas de
acuerdo a la figura 3.1, los anchos de los materiales se eligen de tal manera que se cumpla
la condicion d/2c 1 donde debe estar en el rango donde existe la respuesta metamaterial. A continuaci
on se muestra una serie de graficas donde se comparan los espectros
de transmisi
on para estructuras peri
odicas donde el ancho del dielectrico y del metamaterial es constante con aquellas estructuras donde la capa de metamaterial va creciendo en
una unidad, las estructuras est
an hechas de tal forma que se cumple en ambos casos la
condicion de que la cantidad de dielectrico y metamaterial son iguales para los dos tipos
de estrutura. En la descripci
on de cada grafica el n
umero de capas se denota por la letra
(n), el espesor de las capas de metamaterial para el caso periodico se denotan con (da), y
en el caso en que var
a el ancho del metamaterial, el espesor inicial se denota por (di), el
espesor final del metamaterial (df) y el espesor de la capa de dielectrico por (db), el cual
es el mismo en toda la estructura.

Figura 3.2: Comparacion de espectros de transmision entre estructuras periodicas y aquellas donde la
capa de metamaterial se incrementa, en la figura para ambos casos n=6, db=2 cm, a) caso peri
odico
da=2.25 cm; b) caso no peri
odico di=1 cm y df=4 cm

28

CAP
ITULO 3. RESULTADOS

3.1. ESPECTROS DE TRANSMISION

Figura 3.3: Comparacion entre los espectros de transmision entre estructuras multicapas donde en un
caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene constante para esta figura
n=8, db=2 cm; a) caso peri
odico da= 2.31 cm b) caso no peri
odico di = 1 cm , df = 4 cm

Figura 3.4: Espectros de transmision para estructuras multicapas de metamateriales y dielectricos, en

un caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene constante para esta
figura n=8, db=2 cm; a) caso peri
odico da= 2.46 cm b) caso no peri
odico di = 0.5 cm , df = 5 cm

29

CAP
ITULO 3. RESULTADOS

3.1. ESPECTROS DE TRANSMISION

Figura 3.5: Espectro de transmision para estructuras multicapas de metamateriales y dielectricos, en un

caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene constante para esta figura
n=10, db=2 cm; a) caso peri
odico da= 2.35 cm b) caso no peri
odico di = 1 cm , df = 4 cm

Figura 3.6: Espectro de transmision para estructuras multicapas de metamateriales y dielectricos, en un

caso el ancho de los metamateriales va cambiando y en el otro caso se mantiene constante para eta figura
se tiene n=8, db=3cm; a) caso peri
odico da=2.31cm b) caso no peri
odico di=1cm, df=4cm

Al analizar las gr
aficas de las figura 3.2 a la 3.5 se observa que las resonancias de
Fabry-Perot se incrementan con el aumento de n
umero de capas en ambos casos tanto
el periodico como el no peri
odico, al incrementar el n
umero de capas los espectros de
tranmision coinciden a bajas frecuencias, sin embargo al rededor de las frecuencias donde
se cumple la condici
on n1 dm + n2 dd = 0 donde dm es el ancho de la cantidad total de
metamaterial y dd es la cantidad total de dielectrico, los espectros de transmision difieren
uno del otro considerablemente, as mismo los gaps no quedan bien definidos en los casos
no periodicos.
30

CAP
ITULO 3. RESULTADOS
3.2. TIEMPOS DE FASE

3.2.

Tiempos de fase

En las siguientes gr
aficas se muestran los tiempos de fase calculados mediante la
ecuaci
on (2.96) donde nuevamente se comparan las curvas obtenidas para los casos
periodico y no peri
odico, cabe mencionar que en las graficas mostradas en cada figura
la cantidad de metamaterial y dielectrico que atraviesan las ondas electromagneticas
son iguales, esto es, si la estructura no repercutiese, el tiempo de transito de las ondas
electromagneticas para cada caso deberan ser los mismos . Al igual que en la secci
on
anterior en la descripci
on de cada grafica se da el n
umero de capas (n), el espesor de las
capas de metamaterial para el caso periodico se denotan con (da), el espesor inicial del
metamaterial (di), espesor final del metamaterial (df) y el espesor de la segunda capa
que es el dielectrico (db), tambien se adiciona a cada cada grafica el tiempo que tardara
la luz en atravesar si estuviera en el vacio para una estructura de dielectrico.

Figura 3.7: Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=6, db=2 cm, a)
caso peri
odico da=2.25 cm; b) caso no peri
odico di=1 cm y df=4 cm. Para esta estructura se va a obtener
un tiempo en el vacio de 0,849ns

31

CAP
ITULO 3. RESULTADOS
3.2. TIEMPOS DE FASE

Figura 3.8: Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=8, db=2 cm;
a) caso peri
odico da= 2.31 cm b) caso no peri
odico di = 1 cm , df = 4 cm. Para esta estructura se va a
obtener un tiempo en el vacio de 1,14ns

Figura 3.9: Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=8, db=2 cm;
a) caso peri
odico da= 2.46 cm b) caso no peri
odico di = 0.5 cm , df = 5 cm. Para esta estructura se va a
obtener un tiempo en el vacio de 1,18ns

32

CAP
ITULO 3. RESULTADOS
3.2. TIEMPOS DE FASE

Figura 3.10: Tiempo de fase para las estructuras periodicas y no periodicas, para ambos n=10, db=2
cm; a) caso peri
odico da= 2.35 cm b) caso no peri
odico di = 1 cm , df = 4 cm. Para esta estructura se va
a obtener un tiempo en el vacio de 1,44ns

En correspondencia con las curvas de transmision las graficas del tiempo de fase, muestran similares comportamientos, los tiempos de retardo en las regiones de baja transmisi
on
tienden a cero, en tanto que las resonancias de Fabry-Perot implican tiempos grandes de
transmisi
on, el resultado interesante que se observa en todas las curvas es que para las
estructuras no peri
odicas los tiempos de retardo para algunas frecuencias son muy grandes
comparadas con el caso peri
odico, por ejemplo e la figura 3.6 para una frecuencia de 5.48
GHz el tiempo de tr
ansito de la estructura periodica es 10 veces mayor que el tiempo
de tr
ansito de la estructura periodica, por citar otro ejemplo en la figura 3.8 para una
frecuencia de 5.55 GHz, el tiempo para el caso periodico es de 7.8 ns en tanto que en el
caso no peri
odico es de 140 ns, 20 veces mayor el tiempo de transito. Por otra parte en las
regiones cercanas al gap cero, esto es, donde se cumple la condicion n1 dm + n2 dd = 0 los
tiempos de tranmisi
on son reducidos incluso negativos, como ocurre en los gaps de Bragg.

33

Conclusiones
En esta tesis se ha estudiado el tiempo que le toma a la luz atravesar una estructura
construida con capas alternadas de dielectricos y metamateriales, esto es, se trata de
cristales fot
onicos unidimesionales, para ello se revisaron algunos principos de la teora
electromagnetica, partiendo de las ecuaciones de Maxwell y mediante el uso de la
matriz de transferencia es posible calcular los coeficientes de tranmision y reflexion de
la luz, posteriormente usando el metodo de Zhota para el calculo del tiempo de fase, se
encontr
o el tiempo que le toma a la luz atravesar la estructura. Los cristales fotonicos
discutidos en este trabajo fueron hechos a base de alternar dielectricos y metamteriales
en donde el ancho de las capas del metamaterial va creciendo a lo largo de la estructura,
la motivaci
on surgi
o a raz de resultados publicados donde se encuentra que en cristales
fotonicos de capas alternadas de dielectricos y metales, cuando el ancho de metal crece
en la estructura, se ha encontrado que la luz queda atrapada momentaneamente en el
cristal fot
onico y esto ocurre cerca de la region de frecuencia donde se cumple la condici
on
a da + b db = 0. En el caso de los metamateriales se estudio este comportamiento donde se
cumple la condici
on na da + nb db = 0 efecto solo observable en metamateriales porque s
olo
estos pueden tener un ndice de refraccion negativo. Con el fin de comparar el efecto que
produce el crecimiento suscesivo de las capas de metamaterial, se hicieron comparaciones
con cristales fot
onicos donde la estructura es periodica, esto es, el ancho de las capas
de metamaterial y de dielectrico se mantienen constantes. Las comparaciones fueron
hechas entre estructuras peri
odicas y estructuras donde el ancho del metamaterial se va
incrementando, con la condici
on de que la cantidad total de metamaterial y dielectrico sea
la misma para los dos casos. Los resutlados observados fueron los siguientes siguientes:

Los especros de transmisi

on muestran una regularidad para el caso periodico aun en
la regi
on de frecuencias donde se cumple la condicion na da + nb db = 0 a diferencia
de los espectros de transmision del caso no periodico que reducen la transmisi
on
cerca de las frecuencias del gap cero, como puede verse en las figuras 3.2 a 3.5.
Al parecer el rompimiento de la periodicidad, genera estados localizados, en estas
frecuencias la transmisi
on es reducida.

Las resonancias de Fabry Perot son similares para ambos tipos de estructuras peri
odicas o no peri
odicas, incluso en la medida que el n
umero de capas se incrementa,
estas tienden a coincidir en regiones fuera del gap cero esto puede verse en la figura
3.5. En tanto que para frecuencias cercanas al gap cero la diferencia de los espectros
35

de transmisi
on se hace m
as evidente.
Si se parte del hecho que las estructuras mostradas en las figuras 3.6 a 3.9 contienen
la misma cantidad de dielectrico y de metamaterial, esperaramos resultados
similares al medir el tiempo de tr
ansito, sin embargo se puede observar que mientras
en las regiones de baja frecuencia los tiempos de fase coinciden, en la medida
que se acercan a frecuencia correspondientes al gap cero los tiempos de fase de
la estructura no peri
odica se reducen, llegando a valores incluso negativos, este
resultado puede ser atribuido a una reduccion de la amplitud de tranmison, tal como ocurre en los gaps de Bragg, y que ha sido ampliamente discutido en la literatura.
En el rango de frecuencias cercanas al gap cero, existen tiempos de fase muy grandes comparados con el caso periodico, incluso son mayores hasta en dos ordenes
de magnitud, para estos valores se puede decir que la luz se detiene momentaneamente, generando un efecto similar al reportado en las redes construidas a base de
metales y dielectricos donde se habla del arcoris atrapado. Se puede decir que los
metamateriales tambien son factibles de atrapar la luz.

Bibliografa
[1] V. Veselago, Electrodynamics of substances with simultaneously negative values of
and , Sov. Phys. Usp., Vol.10, 509, 1968.
[2] J. Pendry; A. Holden; D. Robbins; W. Stewart, Magnetism from Conductors and
Enhanced Nonlinear Phenomena, IEEE Trans. Microwave Theory Tech., Vol. 47,
2075, 1999.
[3] G. Eleftheriades; K. Balmain, Negative-Refraction Metamateriales: Fundamental
Principle and Applications, John Wiley and Sons, Inc., Hoboken , New Jersey, 2005
[4] 09 Monografas del SOPT: Los metamateriales y sus aplicaciones en defensa,
075-11-024-3, marzo 2011
[5] K. Tsakmakidis; A. Broardman; O. Hess, Trapped rainbow storage of light in
metamateriales, Nature, Vol 459 Noviembre 2005
[6] T. Jun Cui; D. Smith; R. Liu, Metamaterials: Theory, Design and Applications,
Springer New York Dordrecht Heidelberg London, 2009
[7] F. Perez; I. Solis; M. Palomino, Elastic superlattices with simultaneously negative
effective mass density and shear modulus, Journal of Applied Physics, ISSN:
0021-8979, 113, 093508, (2013)
[8] W. Cai; V. Sholaev, 0ptical Metamaterials: Fundamentals and Applications, Springer New York Dordrecht Heidelberg London, 2010
[9] A. Sihvola, Metamaterials; A Personal View, Radioengineering Vol 18 No. 2 June
2009
[10] G. Stekolschik, Los metamateriales: Luz obediente, Tecnologa Exacta mente, A
no
18 No. 48, septiembre 2011

37

[11] J. Pendry; A. Halden; D. Robbins; W. Stewart, Magnetism from Conductors and

Enhanced Nonlinear Phenomena, 0018-9480/99, IEEE Vol. 47 No. 11, noviembre
1999
v
[12] S. Dutta Gupta, R. Arum; G. Agarwal, Subluminal to superluminal propagation in
a left-handed medium, 10.1103/PhysRevB.69.113104, marzo 2004
[13] H. K. Ramrez Estudio de efectos de absorci
on en la propagaci
on de ondas electromagneticas en medios estratificados que contienen metamateriales, DOCTORADO
EN FISICA APLICADA (FCFM-BUAP), 9 de diciembre de 2011.
[14] H. Kinto; M. Palomino; F. Ramos, Photonic modes in dispersive and lossy superlattices containing negative-index materials, Progress In Electromagnetics ResearchPIER B, Vol. 35, 133-149, (2011) ISSN, 1559-8985 EMW Publishing, direccion
electronica de la revista: http://www.jpier.org/PIERB/pier.php?paper=11062911
[15] J. David Jackson ,Classical Edectrodynamics, 3rd edition, John Wiley and Sons,
Hoboken , New Jersey, 1999
[16] P, Gomez; E. Gonz
ales, Las ecuaciones de Maxwell, Creative commons, 2012
[17] O. J. Bressan, 0ndas, The General Secretarial of the organization of America States
Washintong D.C., 1975
[18] R. K. Wangsness, Campos Electromagneticos, LIMUSAS.A. de C.V. Grupo Noriega
Editores Balderas 95, Mexico D.F., 1998
[19] Y. Zohta, Solid State Commun, Vol 73, 845, 1990
[20] Y. Zohta, Phys, Rev. B41, 7879, 1990
[21] G. O. Mendoza, Tiempo de retardo en estructuras peri
odicas construidas con
metamateriales, LICENCIATURA EN FISICA (FCFM-BUAP), 29 de agosto de
2007
[22] C. Aguilar; M. Palomino, Transmisi
on de paquetes de ondas electromagneticas a
traves de estructuras peri
odicas con respuesta
optica no lineal, Superficies y Vaco
17(2), 9-14, junio 2004

[23] V. de Aquino; V. Aguilera; M. Goto; H. Iwamoto, Tunneling time through a

rectangular barrier, Physical Review Vol 58 Num. 6, diciembre 1998
[24] C. Aguilar; M. Palomino; F. RAmos, Transit time of electromagnetic Gaussian packets traversing non linear, periodic layered systems, Optical Material, 29, 1667-1671.
(2007). ISSN: 0925-3467