1
Metody otrzymywania i właściwości optyczne
materiałów z ujemnym współczynnikiem załamania
Praca dyplomowa magisterska
Milena Dziębaj
Opiekun:
dr hab. inŜ. Włodzimierz Salejda prof. PWr.
Wrocław 2006
Wydział Podstawowych Problemów Techniki
2
Serdecznie dziękuję
Panu profesorowi Włodzimierzowi Salejdzie
za pomoc, cenne uwagi merytoryczne
i wsparcie w trakcie pisania tej pracy
3
Spis treści
Cel pracy ...........................................................................................................................4
Wykaz waŜniejszych oznaczeń i skrótów.........................................................................5
I Wprowadzenie ................................................................................................................6
I.1
Ukryte właściwości równań Maxwella.............................................................. 7
I.2
Ośrodek „dodatni” ........................................................................................... 10
I.3
Ośrodek „ujemny” ........................................................................................... 12
II Wytwarzanie metamateriałów.....................................................................................15
II.1
Pierwsze materiały o ujemnych parametrach .................................................. 15
Powierzchniowy rezonans plazmowy ....................................................... 16
II.2
Tablica długich drutów metalicznych.............................................................. 18
II.3
Rozproszone rezonatory kołowe ..................................................................... 23
II.4
Pierwszy ośrodek o ujemnym współczynniku załamania ............................... 26
II.5
Model linii transmisyjnych.............................................................................. 27
Symulacja rzeczywistego dielektryka ....................................................... 28
Symulacja ujemnego współczynnika załamania ....................................... 34
II.6
Nanostruktury z drutów metalicznych............................................................. 37
II.7
Metamateriały dla zakresu widzialnego .......................................................... 40
III Eksperymenty ............................................................................................................44
III.1
Pierwsze dowody eksperymentalne................................................................. 44
Ujemne załamanie – fikcja czy rzeczywistość? ........................................ 47
III.2
Modulacja transmisji fali EM .......................................................................... 50
III.3
Istnienie fal wstecznych................................................................................... 53
III.4
Ujemne załamanie światła ............................................................................... 54
IV Zastosowania .............................................................................................................58
IV.1
Perfekcyjna soczewka...................................................................................... 58
Idealna soczewka płaska............................................................................ 58
Idealna soczewka sferyczna....................................................................... 61
IV.2
Urządzenia mikrofalowe.................................................................................. 62
IV.3
Najnowsze odkrycia ........................................................................................ 63
V Transmitancja warstwowych układów optycznych ....................................................65
V.1
Ujemne załamanie fali EM .............................................................................. 67
V.2
Polaryzacja typu „s” i „p”................................................................................ 68
V.3
Dyspersja współczynnika załamania ............................................................... 69
V.4
Amplitudowe współczynniki transmisji .......................................................... 69
V Podsumowanie ............................................................................................................74
Dodatek A – Iloczyn wektorowy ....................................................................................76
Dodatek B – Magnetyczny potencjał wektorowy...........................................................76
Dodatek C – Teoria linii transmisyjnych ........................................................................78
Dodatek D – Formalizm macierzowy.............................................................................82
Bibliografia .....................................................................................................................84
4
Cel pracy
Zjawisko ujemnego załamania opisane w 1967 roku przez Viktora Veselago wywołało
oŜywioną dyskusję w świecie naukowym w dziedzinie, której zjawiska i rządzące nimi
prawa wszyscy traktowali juŜ jak dogmat. Opisane przez niego ukryte właściwości równań
Maxwella wymogły ponowną dogłębną analizę bardzo dobrze znanych juŜ obszarów
fizyki.
Celem pracy była analiza syntetyczna osiągnięć naukowych związanych
ze zjawiskiem ujemnego załamania fali elektromagnetycznej oraz z wytwarzaniem
sztucznych materiałów dielektrycznych (metamateriałów) na przestrzeni lat 1967-2006.
Uzasadnione było to niesłabnącym w ostatnich latach zainteresowaniem naukowców tym
zagadnieniem i mnogością publikacji dotyczących tematu. Ze względu na dość liczne
głosy krytyki podwaŜające sam fakt istnienia zjawiska ujemnego załamania, w pracy
przedstawiono wyniki kilku najbardziej przełomowych eksperymentów, potwierdzających
moŜliwość przyjmowania przez współczynnik załamania wartości ujemnych, jak równieŜ
inne ciekawe właściwości metamateriałów.
Omówione zostały dotychczasowe osiągnięcia w dziedzinie zwanej niekiedy „nową
optyką”. Podstawy fizyczne dotyczące zjawiska ujemnego załamania oraz sposób
rozumowania, który doprowadził Veselago do przełomowych wniosków, przedstawione
zostały w rozdziale I. Rozdział II stanowi przegląd technologii wytwarzania
metamateriałów
począwszy
od
pierwszych
prób
uzyskania
takiego
ośrodka,
aŜ do stworzonego w tym miesiącu metamateriału dla zakresu optycznego. Szereg prób
doświadczalnej
weryfikacji
właściwości
projektowanych
ośrodków,
z
których
te najbardziej udane i o kluczowym znaczeniu dla historii zjawiska ujemnego załamania
zebrane zostały w rozdziale III, otworzył drogę do dyskusji na temat potencjalnych
zastosowań sztucznych ośrodków w Ŝyciu mniej lub bardziej codziennym (rozdział IV).
Rozdział V zawiera wyprowadzenie analitycznych wzorów Fresnela na współczynniki
transmitancji i reflektancji dla dwóch warstwowych układów zawierających elementy
ujemne oraz opis zachowania się fali EM na granicy ośrodka dodatniego i ujemnego wraz
ze schematem polaryzacji poszczególnych składowych fali.
5
Wykaz wa
Ŝ
niejszych oznacze
ń
i skrótów
Skróty
LHM – Left-Handed Material; ośrodek „ujemny”;
RHM – Right-Handed Material, ośrodek „dodatni”;
ALMW
–
Array of Long Metallic Wires, tablica długich drutów metalicznych;
SRR – Split-Ring Resonators, rozszczepione rezonatory kołowe;
CSRR – Crossed Split-Ring Resonators, skrzyŜowane rozszczepione rezonatory kołowe;
fala EM – fala elektromagnetyczna;
Oznaczenia
E
r
– wektor natęŜenia pola elektrycznego;
H
r
– wektor natęŜenia pola magnetycznego;
D
r
– wektor indukcji elektrycznej;
B
r
– wektor indukcji magnetycznej;
S
r
– wektor Poyntinga;
k
r
– wektor falowy;
ε
0
– przenikalność elektryczna próŜni;
µ
0
– przenikalność magnetyczna próŜni;
ε
r
– względna przenikalność elektryczna ośrodka;
µ
r
– względna przenikalność magnetyczna ośrodka;
n – współczynnik załamania światła;
ω
– częstość fali elektromagnetycznej;
f – częstotliwość fali elektromagnetycznej;
c – prędkość światła w próŜni;
)
( f
v
– prędkość fazowa fali;
)
( g
v
– prędkość grupowa fali;
d
– grubość warstwy;
r
s
, r
p
– amplitudowe współczynniki odbicia dla polaryzacji s i p;
t
s
, t
p
– amplitudowe współczynniki transmisji dla polaryzacji s i p;
Γ
– macierz charakterystyczna ośrodka;
P – macierz propagacji fali w warstwie dielektrycznej;
D – macierz transmisji fali na granicy ośrodków dielektrycznych;
τ
– ślad macierzy 2 x 2;
σ
– antyślad diagonalny macierzy 2 x 2;
ς
– antysymetryczny antyślad niediagonalny macierzy 2 x 2;
η
– symetryczny antyślad niediagonalny macierzy 2 x 2;
Φ
D
– strumień indukcji elektrycznej;
Φ
B
– strumień indukcji magnetycznej;
j – wektor gęstości prądu.
6
I Wprowadzenie
W ostatnich kilku latach moŜna zauwaŜyć znaczny wzrost zainteresowania nieznanym
dotąd zjawiskiem tak zwanego ujemnego załamania światła. Za inicjatora tej tendencji
powszechnie uwaŜa się rosyjskiego fizyka Victora Veselago, jednakŜe zjawisko to było
przedmiotem zainteresowania naukowców o wiele wcześniej, bo juŜ na początku XX
wieku [1]. Veselago jednak był pierwszym, który swoje przewidywania wyraził otwarcie
dodatkowo popierając je spójnym i pełnym uzasadnieniem oczekiwanych zjawisk.
W 1967 roku
1
Victor Veselago w jednej ze swoich publikacji [1], [2] rozwaŜał
jak zachowywałaby
się
fala
ś
wietlna
padająca
na
wyimaginowany
ośrodek
charakteryzujący się obiema jednocześnie ujemnymi przenikalnościami – elektryczną
i magnetyczną. Rok później praca ta przetłumaczona została na język angielski [4].
Veselago rozwaŜania swe oparł na wnikliwej analizie równań Maxwella, dzięki czemu
odkrył nowe i nieoczekiwane ich właściwości. Konsekwencją tego była teza o istnieniu
nowej grupy materiałów charakteryzujących się nieznanymi dotychczas właściwościami,
które radykalnie zmieniłyby wiele dobrze znanych – jak się wydawało – zjawisk. Przez
wiele lat temat ten nie był poruszany z uwagi na swój jedynie teoretyczny charakter i brak
praktycznych moŜliwości realizacji takiego ośrodka. Jednak od końca XX wieku, kiedy
J.B.Pendry i in. po raz pierwszy opisali obiecujące zastosowania hipotezy Veselago
[5]
−
[7], zjawisko tak zwanego ujemnego załamania nieodmiennie skupia uwagę świata
naukowego.
Niniejszy rozdział zawiera obszerne omówienie teorii wysuniętej przez Victora
Veselago. Podkreślone zostały podstawowe róŜnice między materiałami „dodatnimi”
i „ujemnymi”.
1
W większości publikacji wymieniany jest jednak błędnie rok 1968 jako data pierwszej publikacji Victora
Veselago na ten temat. Obie prace (w języku rosyjskim z roku 1967 jak i w języku angielskim z roku
1968) dostępne są na stronie autora [3].
7
I.1 Ukryte wła
ś
ciwo
ś
ci równa
ń
Maxwella
Podstawowymi równaniami elektrodynamiki są równania zaprezentowane w 1873 r. przez
szkockiego matematyka i fizyka Jamesa Clerka Maxwella [8]
−
[10]. Celem Maxwella było
przedstawienie zjawiska elektromagnetyzmu w jak najprostszy i jednolity sposób.
Równania te – znane dziś jako równania Maxwella – mają następujące postacie
2
Postać
róŜniczkowa
Postać
całkowa
Sens
fizyczny
V
ρ
=
⋅
∇
D
r
Q
dV
d
V
S
=
⋅
=
⋅
∫
∫
ρ
ε
s
E
r
r
0
Prawo Gaussa dla elektryczności –
ź
ródłem pola elektrycznego są ładunki
t
∂
∂
−
=
×
∇
B
E
r
r
dt
d
d
B
L
Φ
−
=
⋅
∫
l
E
r
r
Prawo indukcji Faradaya
– zmienne
w czasie pole magnetyczne wytwarza
wirowe pole elektryczne
0
=
×
∇
B
r
∫
=
⋅
S
d
0
s
B
r
r
Prawo Gaussa dla magnetyzmu
– pole
magnetyczne jest bezźródłowe, linie
pola magnetycznego są zamknięte
t
∂
∂
+
=
×
∇
D
j
H
r
r
r
+
Φ
=
⋅
∫
I
dt
d
d
E
L
0
0
ε
µ
l
B
r
r
Prawo Ampere’a
–Maxwella – zmienne
pole elektryczne oraz przepływający
prąd
wytwarzają
wirowe
pole
magnetyczne
Dodatkowo, równania materiałowe mają postać
E
E
D
r
r
r
r
ε
ε
ε
0
=
=
,
(1.1)
H
H
B
r
r
r
r
µ
µ
µ
0
=
=
.
(1.2)
Na granicy dwóch ośrodków fala elektromagnetyczna musi spełniać następujące warunki
ciągłości
składowych
stycznych
wektorów
natęŜenia
pola
elektrycznego
E
r
i magnetycznego H
r
i składowych normalnych wektorów indukcji elektrycznej
D
r
(1.1)
i magnetycznej
B
r
(1.2)
:
⊥
⊥
⊥
⊥
=
=
=
2
2
2
1
1
1
D
E
E
D
ε
ε
⊥
⊥
=
2
1
B
B
||
2
||
1
E
E
=
||
2
2
||
2
1
||
1
||
1
H
B
B
H
=
=
=
µ
µ
(1.3)
gdzie
⊥
⊥
⊥
⊥
2
1
2
1
,
,
,
D
D
E
E
to składowe wektora natęŜenia i indukcji pola elektrycznego
normalne do kierunku propagacji fali elektromagnetycznej odpowiednio dla ośrodka 1 i 2,
analogicznie
⊥
⊥
2
1
,
B
B
–
składowe
prostopadłe
wektora
indukcji
magnetycznej,
zaś
||
2
||
1
,
E
E
,
||
2
||
1
||
2
||
1
,
,
,
H
H
B
B
– składowe styczne.
2
Spis uŜywanych w pracy oznaczeń znajduje się na stronie 5.
8
Victor Veselago [4] zauwaŜył dwa dodatkowe rozwiązania znanej równości opisującej
związek współczynnika załamania ośrodka i jego przenikalności elektrycznej
i magnetycznej
r
r
n
µ
ε
⋅
=
2
.
(1.4)
Dopuszczając wartości zespolone, uzyskał cztery moŜliwe pierwiastki powyŜszego
równania
r
r
n
µ
ε
⋅
+
=
,
)
(
)
(
r
r
n
µ
ε
−
⋅
−
+
=
r
r
n
µ
ε
⋅
−
=
,
)
(
)
(
r
r
n
µ
ε
−
⋅
−
−
=
(1.5)
Aby sprawdzić, które z powyŜszych moŜliwości są dopuszczalne, Vesalago rozpatrzył
równania Maxwella dla monochromatycznej fali płaskiej o postaci
)
(
0
t
i
e
⋅
−
⋅
⋅
⋅
=
ω
r
k
E
E
r
r
r
r
,
(1.6)
)
(
0
t
i
e
⋅
−
⋅
⋅
⋅
=
ω
r
k
B
B
r
r
r
r
.
(1.7)
Równania Maxwella dla takiej fali prowadzą do równości
H
E
k
r
r
r
⋅
⋅
⋅
=
×
ω
µ
µ
r
0
,
(1.8)
E
H
k
r
r
r
⋅
⋅
⋅
−
=
×
ω
ε
ε
r
0
,
(1.9)
gdzie E
r
– wektor natęŜenia pola elektrycznego, H
r
– wektor natęŜenia pola
magnetycznego, kˆ – wersor na kierunek k
r
, k
r
– wektor falowy
k
k
ˆ
⋅
⋅
=
c
n
ω
r
(1.10)
Uwzględniając wzór
(1.10) otrzymujemy następujące równania:
H
E
k
r
r
⋅
⋅
⋅
=
×
⋅
⋅
ω
µ
µ
ω
r
c
n
0
ˆ
,
(1.11)
E
H
k
r
r
⋅
⋅
⋅
−
=
×
⋅
⋅
ω
ε
ε
ω
r
c
n
0
ˆ
,
(1.12)
z których wynikają następujące cztery moŜliwe rozwiązania:
•
jeŜeli
0
>
r
µ
, to
0
>
n
;
•
jeŜeli
0
<
r
µ
, to
0
<
n
;
9
•
jeŜeli
0
>
r
ε
, to
0
>
n
;
•
jeŜeli
0
<
r
ε
, to
0
<
n
.
Następnie zestawiając wszystkie przypadki otrzymujemy tylko dwie moŜliwości zgodne
z równaniami Maxwella i nie zmieniające ich postaci:
•
gdy
0
>
r
µ
i
0
>
r
ε
, to
0
>
n
(1.13)
•
gdy
0
<
r
µ
i
0
<
r
ε
, to
0
<
n
(1.14)
Oznacza to, iŜ równania Maxwella dopuszczają dwa spośród przytoczonych wcześniej
zespolonych równań na bezwzględny współczynnik załamania ośrodka
r
r
n
µ
ε
⋅
+
=
,
(1.15)
)
(
)
(
r
r
n
µ
ε
−
⋅
−
−
=
.
(1.16)
Pierwsza z powyŜszych zaleŜności (1.15) odpowiada ośrodkowi określanemu przez
Veselago mianem prawoskrętnego, zaś druga (1.16) tak zwanemu ośrodkowi
lewoskrętnemu.
Opisywane tu materiały, ze względu na swoją krótką historię, nie mają jeszcze
jednoznacznie ustalonej nazwy. W głównych pozycjach literaturowych [4] zauwaŜyć
moŜna pewną dowolność w tej kwestii. Ośrodki charakteryzujące ujemne załamanie
nazywane są na przykład materiałami lewoskrętnymi
3
, jednak określenie to zostało
juŜ uŜyte w opisie ośrodków chiralnych [12]. Innym zaproponowanym terminem jest
materiał wsteczny
4
uŜyty przez Lindell’a i in. [13], jednak termin ten z góry narzuca
definicję kierunku wstecznego oraz stwarza problemy przy opisie innych niŜ płaskie czoła
fali. Ziółkowski i Heyman w swojej pracy [14] uŜywają określenia materiał podwójnie
ujemny
5
, które wyraźnie sugeruje jednoczesną ujemność rzeczywistych składowych
przenikalności elektrycznej i magnetycznej, jednak nie podkreśla dostatecznie znaczenia
efektów rozpraszania. Kolejnymi spotykanymi w wielu publikacjach terminami są materiał
o ujemnym współczynniku załamania
6
oraz materiał o ujemnej prędkości fazowej
7
, które
są dość trafnymi nazwami dla tego typu ośrodków. W tej pracy uŜyte zostały dwa ostatnie
określenia, jak równieŜ ośrodek „dodatni” oraz ośrodek „ujemny” intuicyjnie odnoszące
się do dodatniego i ujemnego kąta załamania w omawianych ośrodkach. Ponadto naleŜy
mieć na uwadze, iŜ uŜywane tu określenie metamateriał domyślnie oznacza sztuczny
materiał ujemnie załamujący fale elektromagnetyczne.
3
LHM (left-handed medium)
4
BW (backward medium)
5
DNG (double negative medium)
6
NIM (negative index medium)
7
NPV (negative phase-velocity medium)
10
I.2 O
ś
rodek „dodatni”
Ośrodki, których współczynnik załamania opisać moŜna wzorem (1.15), zwane przez
Veselago ośrodkami prawoskrętnymi, stanowią dobrze poznaną grupę powszechnie
istniejących materiałów. Propagacja fal elektromagnetycznych przez takie ośrodki
jest przedmiotem zainteresowania między innymi optyki czy akustyki i została juŜ
wielokrotnie i wyczerpująco omówiona [8]
−
[10], [15]. Zachowanie fali EM w takim
ośrodku określają prawa odbicia i załamania [10], zaś wektory
E
r
,
H
r
i
k
r
opisujące falę
tworzą prawoskrętną trójkę (Rys. 3).
Szybkość, z jaką fala elektromagnetyczna przemieszcza się w przestrzeni określić
moŜna mierząc jak zmienia się połoŜenie pewnego jej fragmentu, czyli jak szybko
w ośrodku przemieszcza się jej faza [10]. Mówimy wtedy o prędkości fazowej fali
)
(
f
v
i taki opis dobrze sprawdza się w przypadku fali monochromatycznej. JeŜeli mamy
do czynienia z falą modulowaną – jaką jest rzeczywista fala elektromagnetyczna – złoŜoną
z kilku sinusoidalnych fal składowych o róŜnej częstotliwości, to prędkość propagacji
energii moŜe być inna niŜ prędkość fazowa fal składowych. Wtedy mówimy o prędkości
przemieszczania się obwiedni lub paczek falowych, czyli o prędkości grupowej fali
)
(
g
v
(Rys. 1).
Rys. 1 Dwie fale sinusoidalne y
1
i y
2
o zbliŜonych częstotliwościach i długościach fal oraz
obwiednia ich sumy (linia przerywana) rozchodzi się z prędkością grupową [16].
Prędkość fazową moŜna opisać zaleŜnością
π
λ
ω
λ
ω
2
)
(
=
⋅
=
=
f
k
v
f
.
(1.17)
PoniewaŜ dla fal elektromagnetycznych częstość
ω
zaleŜna jest od wektora falowego
(długości fali)
)
(
)
(
k
n
ck
k
=
ω
,
(1.18)
11
więc prędkość fazowa wyraŜa się jako
)
(
)
(
k
n
c
v
f
=
,
(1.19)
gdzie
)
(k
n
jest współczynnikiem załamania dla danej liczby falowej
λ
π
2
=
k
.
Analizując wzór (1.19) moŜna zauwaŜyć, Ŝe gdy n<1 prędkość fazowa moŜe
przekroczyć prędkość światła. Nie oznacza to jednak moŜliwości przekazu informacji
z prędkością większą niŜ prędkość światła
8
. Szybkość przepływu informacji określa
prędkość grupowa, rozumiana jako prędkość przemieszczania się fali modulowanej. Jako
Ŝ
e prędkość grupowa zawiera w sobie informację o tempie transportu energii, więc
z punktu widzenia dynamiki ma ona większe znaczenie niŜ prędkość fazowa.
Prędkość
grupową określa równanie
ω
ω
ω
∂
∂
+
=
∂
∂
=
n
n
c
k
v
g )
(
(1.20)
i jest ona mniejsza od prędkości światła c , czyli teoria względności nie zostaje naruszona.
Prędkość grupowa moŜe osiągać małe wartości
c
v
g
<<
)
(
dla
1
>>
∂
∂
ω
ω
n
. Mówimy
wówczas o świetle spowolnionym, co jest zagadnieniem intensywnie badanym w ostatnich
latach [17].
Poglądowe symulacje dotyczące prędkości grupowej i fazowej zamieszczone są
w Internecie na stronach [18]
−
[20].
8
Fala sinusoidalna ma z góry ustaloną postać na początku i na końcu kanału transmisyjnego, więc nie moŜna
w niej zawrzeć Ŝadnej informacji [16].
12
I.3 O
ś
rodek „ujemny”
Materiał zaproponowany w 1967 r. przez Vesalago [1] był fikcją naukową –
charakteryzował się ujemnym współczynnikiem załamania. Cecha ta powoduje, Ŝe faza
fali przechodzącej przez taki ośrodek zmniejsza się zamiast zwiększać, jak to się dzieje
w ośrodkach prawoskrętnych. Vesalago podkreślał, Ŝe ta podstawowa róŜnica między
ośrodkami prawoskrętnymi i lewoskrętnymi ma decydujący wpływ na niemalŜe wszystkie
zjawiska elektromagnetyczne. Wiele niespotykanych dotąd dla ośrodków prawoskrętnych
efektów badanych jest do dziś [21].
Najłatwiej zauwaŜalnym efektem, wynikającym z ujemnej wartości współczynnika
załamania, jest zmiana kąta ugięcia się fali załamanej na granicy ośrodków o przeciwnych
znakach współczynnika załamania – doznaje ona ugięcia po tej samej stronie normalnej
(Rys. 2). JeŜeli wartości współczynnika załamania obu materiałów są jednakowe, fala
załamana moŜe nie być w ogóle obecna [22].
V
p
(f)
V
p
(g)
V
L
(f)
V
L
(g)
P
L
OP
n
P
> 0
OL
n
L
< 0
Rys. 2 Załamanie fali elektromagnetycznej na granicy ośrodka prawo- i lewoskrętnego.
JeŜeli porównamy przytoczone wcześniej równania Maxwella (1.11) i (1.12)
dla ośrodka prawoskrętnego, któremu odpowiada dodatnia wartość współczynnika
załamania (1.15) oraz dla ośrodka lewoskrętnego o ujemnym współczynniku (1.16),
odkryjemy podstawową róŜnicę między tymi ośrodkami.
13
Analizując równania (1.11) i (1.12) dla pierwszej z dopuszczonych moŜliwości (1.13),
gdzie wszystkie parametry są większe od zera, postać przytoczonych równań nie zmienia
się. Po opuszczeniu wartości skalarnych, moŜna zapisać
H
E
k
r
r
=
×
ˆ
,
(1.21)
E
H
k
r
r
=
×
ˆ
,
(1.22)
co zgodnie z właściwościami iloczynu wektorowego (Dodatek A) oznacza prawoskrętność
trójki wektorów
E
r
,
H
r
i
k
r
jak zostało to pokazane na Rys. 3.
Rys. 3 Wzajemne połoŜenie wektorów natęŜenia pola elektrycznego, magnetycznego,
wektora falowego oraz wektora Poyntinga w ośrodkach „dodatnich”.
JeŜeli natomiast weźmiemy pod uwagę drugą moŜliwość (1.14), gdzie wszystkie
parametry są mniejsze od zera, równania (1.11) i (1.12) przyjmą postać
H
E
k
r
r
⋅
⋅
−
⋅
=
×
⋅
⋅
−
ω
µ
µ
ω
)
(
ˆ
0
r
c
n
,
(1.23)
E
H
k
r
r
⋅
⋅
−
⋅
−
=
×
⋅
⋅
−
ω
ε
ε
ω
)
(
ˆ
)
(
0
r
c
n
.
(1.24)
Analogicznie po opuszczeniu wartości skalarnych zapisać moŜna jako
H
E
k
r
r
=
×
ˆ
,
(1.25)
E
H
k
r
r
−
=
×
ˆ
.
(1.26)
Zgodnie ze wspomnianymi juŜ właściwościami iloczynu wektorowego odpowiada
to lewoskrętnej trójce wektorów
E
r
,
H
r
i
k
r
(Rys. 4).
14
Rys. 4
Wzajemne połoŜenie wektorów natęŜenia pola elektrycznego, magnetycznego,
wektora falowego oraz wektora Poyntinga w ośrodkach
„ujemnych”.
Wektor Poyntinga, definiowany jest za pomocą iloczynu wektorowego
H
E
B
E
S
r
r
r
r
r
×
=
×
=
0
1
µ
,
(1.27)
gdzie E
r
,
B
r
i H
r
są chwilowymi wartościami pola elektromagnetycznego w danym
punkcie. Po raz pierwszy wprowadzony został przez Johna Henry’ego Poyntinga [8]-[10]
i określa kierunek transportu energii przez falę elektromagnetyczną. Jednostką wektora
Poyntinga w układzie SI jest
s
m
J
2
. Za jego pomocą opisać moŜna szybkość przepływu
energii płaskiej fali elektromagnetycznej przez jednostkową powierzchnię. Jak widać
na Rys. 3, w ośrodku dodatnim kierunek i zwrot wektora Poyntinga
S
r
są takie same
jak wektora falowego k
r
, co oznacza Ŝe energia propaguje się zgodnie z kierunkiem
rozchodzenia się fali. W ośrodku ujemnym (
Rys. 4
) kierunki wektorów
S
r
i k
r
są zgodne,
ale ich zwroty przeciwne, co oznacza wsteczną propagację energii fali. Fakt ten stanowi
kolejną istotną róŜnicą między ośrodkiem dodatnim i ujemnym. Fala taka określana jest
mianem fali wstecznej
9
i zjawisko to opisywał juŜ poczynając od roku 1904 H.Lamb [23],
A.Schuster [24], H.C.Pocklington [25], L.I.Mandel’shtam [26], [27] oraz D.V.Sivukhin
[28].
Jak zostało powyŜej zaznaczone, za przenoszenie informacji w fali odpowiada
prędkość grupowa, zatem jej zwrot wskazuje kierunek propagacji energii przy przejściu
fali EM przez granicę ośrodków dodatniego i ujemnego. Energia w takim układzie
propaguje się tak samo jak w przypadku układu zbudowanego tylko z materiałów
dodatnich, zatem zwrot prędkości grupowej przy przejściu przez granicę ośrodków nie
zmienia się. Prędkość grupowa
)
(
1
g
v
i prędkość fazowa
)
(
1
f
v
w ośrodku lewoskrętnym
(izotropowym) są równe co do wartości, lecz antyrównoległe. Przechodząc z ośrodka
dodatniego do ujemnego prędkość fazowa zmienia zwrot na przeciwny (Rys. 2).
9
BW – backward wave
15
II Wytwarzanie metamateriałów
Rozdział stanowi podsumowanie dotychczasowych dokonań naukowych dotyczących
metod wytwarzania metamateriałów o ujemnym współczynniku załamania począwszy
od pierwszych prób budowy kompozytów o moŜliwych do projektowania właściwościach
fizycznych aŜ do najnowszych publikacji związanych z tematem pracy. Omówiono dające
się zauwaŜyć analogie pomiędzy ośrodkami naturalnymi a metamateriałami oraz ewolucję
struktury przestrzennej wytwarzanych kompozytów pozwalającą na rozszerzenie
operacyjnego zakresu częstotliwości aŜ do zakresu optycznego.
W naturalnych materiałach dielektrycznych lokalne (mikroskopowe) oddziaływania
elektromagnetyczne między tworzącymi je atomami i cząstkami powodują efekty
na większą (makroskopową) skalę w postaci parametrów materiałowych znanych jako
przenikalność elektryczna i przenikalność magnetyczna. Aby parametry te miały znaczenie
praktyczne, wykluczone musi być zjawisko dyfrakcji fali EM na strukturze materiału,
co sprowadza się do warunku
λ
>>
d
, gdzie d to wymiary elementarnej komórki
tworzącej sieć krystaliczną danego materiału. Dla naturalnych dielektryków i fali z zakresu
widzialnego warunek ten jest spełniony, bowiem długość fali świetlnej jest o wiele rzędów
większa od rozmiaru atomu i komórki elementarnej. JeŜeli rozwaŜymy jednak przypadek
promieniowania rentgenowskiego, którego
d
≈
λ
10
, obserwować będziemy efekty
dyfrakcyjne. Początkowo wytworzenie jakiegokolwiek sztucznego materiału mogącego
symulować materiał naturalny wydawało się nieosiągalne właśnie ze względu na skalę,
jaką naleŜy osiągnąć, aby wyeliminować dyfrakcję światła. JeŜeli jednak uŜyjemy
promieniowania z zakresu mikrofalowego, których długość fali jest rzędu centymetrów,
stworzenie sieci z komórek o rozmiarach mniejszych od długości takich fal nie jest juŜ
abstrakcją. To właśnie miał na myśli Winston E. Kock wprowadzając po raz pierwszy
w 1949 roku w swojej pracy [29] termin „sztuczny dielektryk”. Opisał on
elektromagnetyczne struktury o moŜliwych praktycznie do wytworzenia rozmiarach,
które naśladowałyby sposób oddziaływania naturalnych kryształów z promieniowaniem
EM. Dwa lata wcześniej w swoich pracach [30], [31] Kock badał wielkogabarytowe
anteny wykorzystujące układ płaskich soczewek metalicznych, początkowo nie zdając
sobie sprawy z analogii zachodzących między jego metalicznymi soczewkami
a naturalnymi dielektrykami. Kiedy jego późniejsze badania sztucznych dielektryków
złoŜonych
z komponentów
o
róŜnych
kształtach
wykazały
liczne
interakcje
z promieniowaniem EM obserwowane w naturalnych kryształach – idea metamateriałów
przestała być jedynie naukową fikcją.
II.1 Pierwsze materiały o ujemnych parametrach
Metamateriał jest to sztucznie wytworzony ośrodek, którego właściwości fizyczne
wynikają nie tylko z rodzaju tworzących go elementów, ale głównie z jego struktury
przestrzennej. Inspiracja metamateriałoznawców i technologów hipotezą Veselago [4]
zaowocowała kilkoma propozycjami takich ośrodków [32]–[37]. Głównym załoŜeniem
było podejście do kaŜdego materiału naturalnie występującego w przyrodzie
jak do kompozytu złoŜonego z atomów i cząstek, których rodzaj i wzajemne ułoŜenie mają
decydujący wpływ na elektromagnetyczne właściwości danego ośrodka i sztuczne
10
Dla kryształu NaCl długość boku komórki elementarnej wynosi 0.562737 nm zaś długość fali
promieniowania X jest rzędu 0.1 nm – dla porównania światło Ŝółte ma długość fali równą 589 nm [16].
16
wytworzenie analogicznego materiału tylko na większą skalę [7]. Struktura taka miałaby
być zbudowana z powtarzających się komórek elementarnych (zwanych takŜe
„fotonicznym
atomem”)
o rozmiarach
znacznie
mniejszych
niŜ długość
fali
elektromagnetycznej, dzięki czemu moŜna by traktować ją jak materiał jednorodny oraz
miała być bardzo lekka, dzięki uŜyciu drobnych i cienkich metalowych elementów.
Jednymi z pierwszych propozycji były ośrodki wytworzone przy uŜyciu tablicy długich
i cienkich drutów metalicznych
11
[5], [6], której
0
<
ε
lub na bazie rozproszonych
rezonatorów kołowych
12
o bardzo wysokiej polaryzowalności magnetycznej [7], których
0
<
µ
(przy częstotliwościach z zakresu mikrofalowego). Zadaniem tak zaprojektowanych
struktur miała być symulacja zachowania plazmy, w której upatrywano klucza
do wytworzenia ujemnych przenikalności elektrycznej i magnetycznej.
Powierzchniowy rezonans plazmowy
Plazma jest czwartym stanem skupienia obok stałego, ciekłego i gazowego, w jakim moŜe
znaleźć się materia [38]. W stan plazmy przechodzi gaz, jeŜeli zostanie podgrzany
do temperatury tak wysokiej, Ŝe elektrony uwalniane są z orbit atomowych, czyli następuje
jego jonizacja. Ten sam efekt uzyskać moŜna oddziałując bardzo duŜym polem
elektrycznym lub zmiennym polem magnetycznym. Obecność w plazmie swobodnych
ładunków elektrycznych jest przyczyną tego, Ŝe plazma jest stanem przewodzącym i silnie
oddziałuje z polem elektromagnetycznym. Cechą charakterystyczną tego stanu jest funkcja
przenikalności elektrycznej, która poniŜej pewnej częstotliwości zwanej częstotliwością
plazmową, przyjmuje wartości ujemne, co skutkuje urojoną wartością stałej propagacji
energii. Ta właściwość plazmy zaowocowała zainteresowaniem naukowców moŜliwością
modelowania ujemnej wartości przenikalności elektrycznej w sztucznych dielektrykach
wykorzystując zjawisko rezonansu plazmowego.
Plazmon jest kwazicząstką powstałą w wyniku kwantyzacji drgań gęstości ładunku
plazmy na powierzchni metalu. Plazmon powierzchniowy jest elektromagnetyczną falą
powierzchniową, o polaryzacji typu p, propagującą się wzdłuŜ powierzchni styku dwóch
materiałów, których stałe dielektryczne mają przeciwne znaki [10]. Związane z nią jest
zanikające wykładniczo pole elektromagnetyczne prostopadłe do kierunku propagacji fali.
Zmieniając kąt padania lub długość fali promieniowania, które są funkcją współczynnika
załamania, moŜliwe jest wzbudzenie powierzchniowego rezonansu plazmowego
13
, czego
rezultatem jest skokowy spadek intensywności promieniowania odbitego. Powierzchniowy
rezonans plazmowy jest zjawiskiem fizycznym mogącym wystąpić, gdy płasko
spolaryzowana fala elektromagnetyczna z zakresu widzialnego lub bliskiego ultrafioletu
pada na powierzchnię metalu przy spełnionych warunkach całkowitego wewnętrznego
odbicia [10]. Wtedy pomimo, iŜ padające promieniowanie jest całkowicie odbijane przez
powierzchnię metalu, pole elektromagnetyczne fotonów powierzchniowych rozciąga się
poza powierzchnię metalu na odległość około ¼ długości fali promieniowania.
Materiały naturalne, na przykład metale, mogą zostać doprowadzone do stanu
plazmy, jednak dla nich częstość plazmowa jest tak wysoka, Ŝe odpowiadające jej
promieniowanie EM charakteryzuje się tak małą długością fali, Ŝe zbudowanie komórki
o rozmiarach jeszcze od niej mniejszych jest właściwie niemoŜliwe. NaleŜało zatem
11
ALMW – Array of Thin Metallic Wires
12
SRR – Split Ring Resonator
13
SPR – Surface Plasmon Resonance
17
znaleźć materiał, dla którego stan plazmy występuje przy niŜszej częstości. Zastosowanie
tutaj znalazły opisane przez Kock’a sztuczne dielektryki [29], będące przedmiotem
zainteresowania takŜe kilku innych prac naukowych [39]
−
[41]. Podsumowania historii
dotyczącej rozwoju sztucznych dielektryków dostarcza praca [42].
Symulacji zachowań plazmowych przy uŜyciu sztucznych dielektryków jako jeden
z pierwszych podjął się w 1962 roku Walter Rotman [43]. Efektem jego pracy był opis
ośrodka dielektrycznego złoŜonego z drutów zorientowanych zgodnie z kierunkiem
wektora falowego przyłoŜonego pola EM oraz ośrodek zbudowany z przewodzących
pasków metalu. Analiza charakterystyk dyspersyjnych potwierdziła, iŜ taki ośrodek jest
dobrą analogią plazmy.
Rys. 5 Schematyczna ilustracja zjawiska powierzchniowego rezonansu plazmowego
(na podstawie: [44], [45]).
Plazmony powierzchniowe traktować moŜna jako powiązane oscylacje przy powierzchni
metalu, których częstotliwość wyznaczona jest przez
0
)
(
)
(
2
1
=
+
s
s
ω
ε
ω
ε
,
(2.1)
gdzie
1
ε
i
2
ε
to funkcje dielektryczne na obu płaszczyznach oddziaływania metalu
z promieniowaniem. JeŜeli pierwszym ośrodkiem jest próŜnia, a drugim metal i jeŜeli
zaniedbamy rozpraszanie, to
2
p
s
ω
ω
=
,
(2.2)
gdzie częstość plazmowa
p
ω
potrzebna do wywołania tego rezonansu zwyczajowo leŜy
w ultrafiolecie i ma postać
e
p
m
ne
0
2
2
ε
ω
=
,
(2.3)
gdzie
e
m
−
efektywna masą elektronu, zaś n to gęstość elektronów w pojedynczym drucie.
18
Częstość plazmowa nie zaleŜy od długości fali padającego promieniowania, więc cechą
charakterystyczną oscylacji plazmowych jest nieskończona prędkość fazowa i zerowa
prędkość grupowa.
Plazmony mają znaczny wpływ na właściwości metalu i jego sposób oddziaływania
z promieniowaniem
elektromagnetycznym,
gdyŜ
przenikalność
dielektryczna
ε
następująco zaleŜy od częstotliwości plazmowej
)
(
1
)
(
2
γ
ω
ω
ω
ω
ε
i
p
+
−
=
,
(2.4)
gdzie
γ
to współczynnik rozpraszania energii plazmonu w układzie
14
. Przenikalność
elektryczna takiej struktury jest silnie ujemna dla częstotliwości mniejszych niŜ plazmowa.
II.2 Tablica długich drutów metalicznych
W 1996 roku J.B.Pendry i in. [5], [6] zaproponowali sposób na przesunięcie częstotliwości
rezonansowej aktywującej powierzchniowy rezonans plazmowy w metalach nawet o sześć
rzędów wielkości (tj. w zakres GHz). Opisywany przez nich materiał składał się z bardzo
cienkich drutów metalicznych o średnicy rzędu 1
µ
m zestawionych w periodyczną sieć
kubiczną o stałej sieci a (Rys. 6). Struktura taka posiada właściwości nie zaobserwowane
dotąd w paśmie GHz. PoniewaŜ druty metaliczne uŜyte do budowy tego materiału
zajmowały znikomą część kaŜdej komórki elementarnej, sieć taka charakteryzowała się
mniejszą koncentracją elektronów, co pozwoliło uzyskać przyrost efektywnej masy
elektronu
e
m .
Rys. 6 Periodyczna struktura złoŜona z cienkich drutów metalicznych, połączonych
na krawędziach i ułoŜonych w kubiczną sieć, (na podstawie: [5])
Równolegle zagadnienie periodycznych struktur sieci metalicznych badała inna
grupa naukowców [46], jednak nie kładli oni szczególnego nacisku na rozmiary
uŜywanych drutów metalicznych, co według Pendry’ego [5], [6] ma kluczowe znaczenie,
bowiem tylko przy takim załoŜeniu, promieniowanie padające na badaną strukturę moŜe
14
Dla prostych metali jest on pomijalnie mały w porównaniu do częstości plazmowej
p
ω
.
19
wnikać w nią wystarczająco głęboko i jednocześnie nie powodować zjawiska wzajemnej
indukcji między poszczególnymi drutami.
Periodyczne struktury elektromagnetyczne budowane na bazie dielektryków były
kilkukrotnie juŜ opisywane [47]–[50]. Metale natomiast, ze względu na obecne w ich
strukturze elektrony walencyjne umoŜliwiające sprawne przekazywanie energii,
charakteryzują się specyficzną odpowiedzią na promieniowanie elektromagnetyczne,
wiąŜącą się z występowaniem na ich powierzchni plazmowego rezonansu gazu
elektronowego. Idealny metal opisać moŜna przy uŜyciu funkcji dielektrycznej
2
2
.
1
ω
ω
ε
p
ideal
−
=
,
(2.5)
i wynikająca z niego idealna zaleŜność dyspersyjna przedstawiona została na Rys. 7.
PowyŜej częstotliwości plazmowej powstają dwa poprzeczne asymptotyczne mody,
dla częstości równej częstości plazmowej występuje jeden mod podłuŜny, zaś poniŜej
częstości plazmowej obecne są tylko zanikające mody związane z urojoną wartością
wektora falowego i promieniowanie wnika w metal bardzo płytko.
w
ek
to
r
fa
lo
w
y
częstość (w jednostkach częstości plazmowej)
Rys. 7 ZaleŜność dyspersyjna dla światła padającego na idealny metal. [6]
Dla metalu rzeczywistego w zaleŜności (2.5) uwzględnić naleŜy uwzględnić tłumienie
wynikające z rezystancji
(
)
γ
ω
ω
ω
ε
i
1
2
.
+
−
=
p
rzecz
.
(2.6)
W większości znanych metali (poza nadprzewodnikami) tłumienie przyjmuje znaczące
wielkości dopiero w podczerwieni.
Celem było wytworzenie kompozytowego materiału, który powieliłby charakterystyczne
dla metali oddziaływanie z falą elektromagnetyczną, ale dla zakresu GHz. Pendry
rozpatrzył propagację promieniowania wzdłuŜ oś OZ i za aktywne uznał tylko druty
do niej równoległe.
20
Rys. 8 Tablica cienkich drutów metalicznych równoległych do osi z i uporządkowanych
w kwadratową sieć w płaszczyźnie x-y (na podstawie [6])
W takim układzie tylko część przestrzeni wypełniona jest metalem, więc efektywna
gęstość elektronów jest mniejsza niŜ w metalu jednorodnym i wynosi
2
2
a
r
n
n
ef
π
=
,
(2.7)
gdzie n – gęstość elektronów w pojedynczym drucie, r – promień drutu, a – stała siatki.
Dominujące znaczenie ma fakt, Ŝe samoindukcja drutów zaleŜy takŜe od ich
wzajemnego ułoŜenia przestrzennego i maleje logarytmicznie wraz z odległością od osi
drutu. NatęŜenie pola magnetycznego wokół kaŜdego drutu w zaleŜności od odległości R
od jego osi opisać moŜna zaleŜnością
R
nve
r
R
I
R
π
π
π
2
2
)
(
2
=
=
H
,
(2.8)
gdzie I to natęŜenie prądu płynącego przez drut, zaś v to średnia prędkość ruchu
elektronów. Wektorowo pole magnetyczne wokół drutu moŜna opisać jako
)
(
)
(
1
0
R
R
A
H
×
∇
=
−
µ
,
(2.9)
gdzie A to magnetyczny potencjał wektorowy (Dodatek B) równy
=
R
a
ve
r
R
A
ln
2
)
(
2
0
π
π
µ
.
(2.10)
a
r
21
Biorąc pod uwagę dodatkowy wkład elektronów w polu magnetycznym do ich momentu
wektorowego, moment na jednostkę długości drutu wynosi
(
)
v
r
m
r
a
v
n
r
e
r
n
r
e
ef
2
2
2
2
0
2
ln
2
)
(
π
π
π
µ
π
=
=
A
,
(2.11)
gdzie
ef
m jest masą efektywną elektronów daną zaleŜnością
=
r
a
n
r
e
m
ef
ln
2
2
2
0
π
π
µ
.
(2.12)
Masa efektywna elektronów w badanej przez Pendry’ego strukturze była równa
p
e
ef
m
m
m
⋅
=
⋅
×
=
×
=
−
83
,
14
10
2,7233
kg
10
4808
,
2
4
26
gdzie
e
m jest masą elektronu a
p
m
masą protonu. Tak duŜa masa efektywna miała wpływ na zmianę wartości częstości
plazmowej
( )
GHz
20
,
8
s
rad
10
15
,
5
ln
2
1
10
2
2
0
0
2
2
=
⋅
×
=
=
=
−
r
a
a
c
m
e
n
ef
ef
p
π
ε
ω
(2.13)
Realizacja geometryczna badanej struktury przedstawiona została na Rys. 9.
Ze względu na łatwość wytworzenia i cenę, model Pendry’ego składał się
ze skrzyŜowanych pod kątem 90
º
arkuszy polistyrenowych, rozsuniętych na odległość
3 mm, wypełnionych cienkimi, powleczonymi złotem, wolframowymi drutami o średnicy
20
µ
m.
Rys. 9 Realizacja geometryczna struktury o ujemnej przenikalności elektrycznej [6].
Na podstawie równania (2.13) wyjaśnić moŜna dlaczego kluczowe znaczenie
w rozumowaniu Pendry’ego odgrywały małe rozmiary poprzeczne drutów. Gdyby druty
nie były cienkie, czynnik
( )
r
a
ln
byłby bliski 1. Wtedy długość fali w próŜni przy
częstotliwości plazmowej wynosiłaby
( )
π
π
π
ω
π
λ
2
ln
2
2
2
2
1
2
2
0
0
0
0
a
r
a
a
c
c
c
p
p
≈
=
=
−
,
22
zaś głębokość penetracji promieniowania w strukturę byłaby równa
( )
π
3
ln
2
r
a
a
. Zatem
gdyby druty tworzące badany materiał nie miały małego promienia, długość fali
p
0
λ
byłaby rzędu stałej siatki i pojawiłyby się efekty dyfrakcyjne, zaś promieniowanie bardzo
płytko wnikałoby do struktury.
Eksperyment
potwierdził
przewidywane
właściwości
trójwymiarowych
struktur
zbudowanych z cienkich drutów metalicznych. Warunkiem koniecznym okazały się
odpowiednie wymiary drutów, które muszą być wystarczająco długie i cienkie. Taka
geometria struktury zaproponowanej przez Pendry’ego [6] pozwoliła na obniŜenie
częstości plazmowej.
23
II.3 Rozproszone rezonatory kołowe
Przedstawiając hipotezę dotyczącą ujemnego załamania fali EM Veselago [2] zdawał sobie
sprawę z tego, Ŝe uzyskanie ujemnej przenikalności magnetycznej będzie stwarzało więcej
trudności niŜ uzyskanie ujemnej przenikalności elektrycznej, czego główną przyczyną jest
brak dowodów na istnienie elementarnej cząstki magnetycznej – monopolu
magnetycznego. Jednak w 1999 roku Pendry i inni [7] opisali sztuczne ośrodki
o niezwykłych właściwościach magnetycznych. Zaprezentowany przez Pendry’ego
posiadał pojemność, która wraz z naturalnie występującą indukcją wynikającą ze struktury
przestrzennej uŜytych pierścieni, dała efekt w postaci odpowiedzi rezonansowej opisanej
przez efektywną przenikalność magnetyczną
(2.14)
gdzie r – promień pierścienia uŜytego w SRR, a – odległość między osiami sąsiadujących
SRR leŜących w jednej płaszczyźnie, l – odległość między płaszczyznami,
ρ
– oporność
metalu, C – pojemność pomiędzy dwoma płaszczyznami metalu.
Z równania (2.14) wynika, Ŝe materiał złoŜony z tablic SRR charakteryzowałby się ujemną
przenikalnością magnetyczną dla częstości bliskich rezonansowym i ograniczony byłby
tylko rezystywnością uŜytego metalu. Częstość rezonansowa wyraŜa się wzorem
3
0
2
0
3
Cr
l
µ
π
ω
=
.
(2.15)
JeŜeli załoŜymy, Ŝe
0
→
ρ
, to ze wzoru (2.14) wynika, Ŝe ujemne wartości
eff
µ
przyjmie
wówczas, gdy drugi czynnik w mianowniku będzie większy od 1, co odpowiada częstości
plazmowej wyraŜonej jako
−
=
2
2
3
0
2
1
3
a
r
Cr
l
mp
π
µ
π
ω
,
(2.16)
gdzie
F
a
r
=
2
2
π
to współczynnik wypełnienia informujący o stopniu wypełnienia komórki
przez SRR.
Umieszczone w powietrzu tablice SRR mają pasmo zabronione w pobliŜu zakresu
częstości pomiędzy
0
ω
a
mp
ω
co sugerować moŜe, Ŝe w tym zakresie przenikalność
magnetyczna przyjmuje wartości ujemne. Jest to zjawisko wąskopasmowe, jednak moŜna
częstość plazmową umieścić w obszarze GHz, powodując tym samym rozszerzenie
wspomnianego wyŜej zakresu.
3
2
0
2
0
2
2
3
2
1
1
Cr
l
r
lj
a
r
eff
ω
µ
π
µ
ω
ρ
π
µ
−
−
−
=
,
,
24
Typowe rozproszone rezonatory kołowe (SRR) formowane są z pary
przewodzących pierścieni nadrukowanych na dielektryku (Rys. 10). Mechanizm
powstawania ujemnej przenikalności magnetycznej jest następujący: zmienne w czasie
pole magnetyczne przyłoŜone wzdłuŜ osi pierścienia indukuje przepływ prądu, który
w zaleŜności od oporności pierścienia wytwarza przeciwne pole magnetyczne
wzmacniające lub przeciwdziałające polu indukującemu [51]. Pierścienie tworzące SRR
posiadają przerwy, dzięki czemu rezonans moŜe zostać osiągnięty przy uŜyciu długości fali
wielokrotnie przekraczających średnicę pierścieni (typowe wymiary pierścieni to około
jedna dziesiąta długości fali padającego promieniowania). Celem takiej geometrii układu
jest generacja moŜliwie największej pojemności magnetycznej w małej przestrzeni
pomiędzy pierścieniami, co pozwala znacząco obniŜyć częstotliwość rezonansu
i skoncentrować pole elektryczne w obszarze SRR [32].
Rys. 10
Pojedynczy rozszczepiony rezonator kołowy zbudowany z dwóch niepełnych
pierścieni. Przerwy w pierścieniach zorientowane są przeciwnie.
Odpowiedź SRR bezpośrednio zaleŜy od oddziałującego na pierścienie pola
magnetycznego [52]. Pole elektryczne takŜe moŜe mieć swój wkład do rezonansu,
ale zaleŜy on od jego orientacji względem przerw w pierścieniach tworzących SRR.
Wynika z tego, iŜ SRR jest w ogólności strukturą anizotropową, a w celu wytworzenia
jednorodnego izotropowego metamateriału, struktury SRR często formowane są
w kubiczne sieci, jak pokazane to zostało na Rys. 11.
Rys. 11
Sześcienna sieć złoŜona z SRR.
25
Olivier J.F. Martin opisał pokrewne struktury SRR o lepszych właściwościach
magnetycznych [52]. Strukturą magnetyczną wykazującą jeszcze większą izotropię
są skrzyŜowane rozszczepione rezonatory kołowe (CSRR)
15
zbudowane z dwóch
prostopadle ustawionych SRR. KaŜdy rozszczepiony rezonator kołowy zbudowany jest
z dwóch aluminiowych pierścieni o średnicy 15 mm i 20 mm, osadzonych na słabo
przewodzącej piance. Zmierzona odpowiedź elektromagnetyczna CSRR okazała się być
idealnie izotropowa w płaszczyźnie poziomej niezaleŜnie od obrotu wokół osi badanej
struktury [52]. Olivier J.F. Martin zasymulował numerycznie trzy róŜne typy geometrii
CSRR dla dwóch skrzyŜowanych elementów SRR Rys. 12 oraz dwa typy geometrii CSRR
dla trzech skrzyŜowanych elementów SRR Rys. 13.
Rys. 12
Trzy typy geometrii podwójnego CSRR: (a) przerwy dwóch skrzyŜowanych SRR
umieszczone są na tym samym biegunie CSRR; (b) na biegunach przeciwnych;
(c) obrócone o 45
o
kaŜda w przeciwne strony
Izotropią w płaszczyźnie XY wykazały się struktury CSRR (a) i (b) przedstawione
na Rys. 12, zaś struktura (c) jest izotropowa tylko dla pewnych kierunków propagacji
i nie kaŜdy kierunek propagacji jest w stanie w ogóle wywołać jej rezonans. Przyczyną są
róŜnice w rozmieszczeniu przerw w pierścieniach tworzących struktury SRR – dla (a) i (b)
ulokowane są one tak, Ŝe przepływ ładunku nie jest zaburzony, w przypadku (c) następuje
zwarcie.
Rys. 13
Dwa typy geometrii potrójnego CRSS: (a) z trzech identycznych elementów SRR
wzajemnie prostopadłych; (b) z trzech elementów SRR o róŜnych rozmiarach.
Potrójne CSRR charakteryzuje idealna izotropia w trzech wymiarach (niezaleŜnie
od kierunku padania fali elektromagnetycznej) dla geometrii przedstawionej na Rys. 13
(b), ale w przypadku (a) kaŜda ze struktur SRR jest zwarta przez pozostałe dwie.
15
CSRR – Crossed Split Ring Resonators
26
II.4 Pierwszy o
ś
rodek o ujemnym współczynniku załamania
Pierwszy
metamateriał
ujemnie
załamujący
fale
EM
zaproponowany
został
przez D.R.Smith’a [32]
−
[35]. Materiał miał postać trójwymiarowej tablicy zbudowanej
z komórek zawierających rozproszone rezonatory kołowe [6] i cienkie druty metaliczne [7]
(Rys. 14). Opisane struktury wykorzystują rozproszone rezonatory kołowe do wytworzenia
ujemnej przenikalności magnetycznej i drutów metalicznych do wytworzenia ujemnej
przenikalności elektrycznej w pewnych – częściowo pokrywających się – pasmach
częstotliwości. W ten sposób otrzymano okno częstotliwości, gdzie przenikalności
magnetyczna
µ
i elektryczna
ε
są jednocześnie ujemne.
UŜyte przez Smitha rozproszone rezonatory miały postać dwóch kwadratowych
miedzianych pierścieni grubości c = 0,25 mm z przerwą równą g = 0,46 mm, odległych
od siebie o d = 0,3 mm. Zewnętrzny wymiar pojedynczego SRR wynosił w = 2,62 mm,
a uŜyte druty metaliczne miały długość 1 cm. Komórka elementarna zbudowana była
ze skrzyŜowanych, osadzonych na podłoŜu krzemowym sześciu rozproszonych
rezonatorów kołowych (SRR) i dwóch cienkich drutów metalicznych (ALMW).
Rys. 14
(a) pojedynczy rezonator kołowy (SRR); (b) komórka elementarna metamateriału
o stałej siatki 0,5 cm [35].
Dla takiego materiału Smith uŜył następujących wzorów opisujących przenikalności
( )
ωγ
ω
ω
ω
ω
µ
j
F
eff
−
−
⋅
−
=
2
0
2
2
0
1
,
(2.17)
( )
2
2
1
ω
ω
ω
ε
ep
eff
−
=
.
(2.18)
W oparciu o powyŜsze równania zostało udowodnione [53], Ŝe taki ośrodek osadzony
w próŜni posiada analogię w postaci modelu opartego na linii transmisyjnej z elementami
indukcyjnymi L, pojemnościowymi C i rezystancyjnymi (Rys. 15).
27
Rys. 15
Model linii transmisyjnych dla ośrodka zbudowanego z komórek elementarnych
zawierających cienkie druty metaliczne i rozszczepione rezonatory kołowe [65].
Odkrycie tego ośrodka było przełomem w historii metamateriałów. Po weryfikacji
eksperymentalnej [54], która potwierdziła oczekiwania (taki materiał posiada ujemny
współczynnik załamania dla promieniowania mikrofalowego z pewnego zakresu
częstotliwości), coraz śmielej zaczęto wierzyć (choć nie wszyscy [55]), iŜ zaskakująca
teoria Veselago znajdzie swe zastosowanie praktyczne.
II.5 Model linii transmisyjnych
Zachowanie się fali EM przy przejściu przez dowolny ośrodek opisują dobrze znane juŜ
równania Maxwella
B
E
r
r
jω
−
=
×
∇
,
(2.19)
D
J
H
r
r
r
ω
j
+
=
×
∇
,
(2.20)
które w jednorodnymm izotropowym ośrodku uzupełnione są przez równania materiałowe,
ś
ciśle zaleŜne od częstości padającego promieniowania
H
B
r
r
)
(
ω
µ
=
,
(2.21)
E
D
r
r
)
(
ω
ε
=
.
(2.22)
Nowe podejście do tych dobrze znanych podwalin współczesnej fizyki
zaprezentowali prawie równocześnie w 1944 roku Gabriel Kron [62], przedstawiając
jednocześnie numeryczna procedurę rozwiązywania równań Maxwella w tej nowej postaci
[63] oraz J.R.Whinnery i S.Ramo [64]. Udowodnili oni, Ŝe dzięki dyskretyzacji
przestrzennej równań Maxwella moŜna zastosować je dla obwodów RLC, gdzie ich pełną
analogią są równania prądowo-napięciowe Kirchhoffa, co było pierwszym krokiem
w kierunku stworzenia modelu naturalnych ośrodków dielektrycznych opartych
na obwodach RLC.
28
Kluczem do wytworzenia sztucznego dielektryka było zbudowanie takiego modelu,
który umoŜliwiałby odnalezienie bezpośrednich analogii do ośrodków spotykanych
w naturze. Podstawowym załoŜeniem, które naleŜy w tym celu poczyniono, było podejście
do kaŜdego materiału – zarówno sztucznego jak i naturalnego – jak do sieci pewnych
podstawowych elementów o bardzo małych wymiarach (w naturalnych materiałach
są to atomy i cząstki). Analogie te powinny dać się takŜe zauwaŜyć w zachowaniu fali
elektromagnetycznej padającej na wytworzony materiał – fale o długości porównywalnej
z wymiarami pojedynczego elementu sieci (atomu, cząstki, komórki elementarnej)
powinny doznawać efektów dyfrakcyjnych takich, jak ma to miejsce w ciałach stałych, zaś
fale o długości odpowiednio większej od wymiarów jednostkowych komórek powinny
załamywać się, a takŜe dawać moŜliwość zdefiniowania odpowiedniego dla tego
przypadku współczynnika załamania fali elektromagnetycznej.
Symulacja rzeczywistego dielektryka
Model linii transmisyjnych (Dodatek C) mogący reprezentować naturalny ośrodek
przy uŜyciu sieci rozproszonych reaktancji moŜe być złoŜony z komórek przedstawionych
na Rys. 16.
{ V
y
, I
x
}
Oś Z
O
ś
Y
Oś
X
{ V
y
, I
z
}
{ V
y
+dV
y
, I
x
+ dI
x
}
{ V
y
+dV
y
, I
z
+ dI
z
}
½ Z
z
½
Z
x
Y
½ Z
z
½
Z
x
Rys. 16
Elementarna komórka modelu linii transmisyjnej dla płaskiego, jednorodnego
ośrodka dielektrycznego.(na podstawie: [65])
Aby podkreślić periodyczność sieci, równania Maxwella (2.19) – (2.22) dla takiego
przypadku rozwiązuje się oddzielnie dla kaŜdej z komórek, dzięki czemu występujące tam
pola elektryczne i magnetyczne moŜna traktować jako statyczne. Zakładamy, Ŝe istnieje
jednostkowa komórka w przestrzeni o wymiarach
∆
x,
∆
y,
∆
z (Rys. 17), której wymiary
są pomijalnie małe w porównaniu do długości fali padającego promieniowania, co pozwala
traktować pole E i H jako statyczne.
29
x
O
ś
X
Rys. 17
Rozkład quasi-statycznych pól w komórce elementarnej (na podstawie: [65])
W przypadku sieci dwuwymiarowej, moŜna przyjąć, Ŝe pole elektromagnetyczne
jest stałe w kierunku y czyli
0
→
∂
∂
y
. W takim przypadku wszelkie zjawiska
elektromagnetyczne zachodzące w komórce opisać moŜna za pomocą kombinacji modów
TE
y
i TM
y
. JeŜeli rozwaŜymy mod TM
y
, dominującymi składowymi pól elektrycznego
i magnetycznego będą E
y
, H
x
i H
z
z
y
H
j
x
E
)
(
ω
ωµ
−
=
∂
∂
,
(2.23)
x
y
H
j
z
E
)
(
ω
ωµ
+
=
∂
∂
,
(2.24)
y
z
x
E
j
x
H
z
H
)
(
ω
ωε
+
=
∂
∂
−
∂
∂
.
(2.25)
Dyskretyzacja przestrzenna równań Maxwella prowadzi do następujących wyraŜeń:
x
H
j
z
x
E
z
x
x
E
z
y
y
∆
−
=
−
∆
+
)
(
)
,
(
)
,
(
0
0
0
0
ω
ωµ
,
(2.26)
z
H
j
z
x
E
z
z
x
E
x
y
y
∆
+
=
−
∆
+
)
(
)
,
(
)
,
(
0
0
0
0
ω
ωµ
,
(2.27)
[
]
[
]
z
x
y
x
E
j
z
z
x
H
z
x
x
H
x
z
x
H
z
z
x
H
y
z
z
x
x
∆
⋅
∆
⋅
⋅
⋅
⋅
+
=
=
∆
⋅
−
∆
+
−
∆
⋅
−
∆
+
)
,
(
)
(
)
,
(
)
,
(
)
,
(
)
,
(
0
0
0
0
0
0
0
0
0
0
ω
ε
ω
(2.28)
RóŜnicę potencjałów w elementarnej komórce definiuje się jako
∫
⋅
−
=
−
'
'
a
a
a
a
dl
V
V
E
r
,
(2.29)
30
zaś natęŜenie prądu
∫
⋅
=
C
dl
I
H
r
,
(2.30)
gdzie a-a’ to dowolna ścieŜka łącząca dolną i górną ścianę sześciennej komórki
elementarnej, zaś C to odpowiednio wybrany zamknięty przekrój przez powierzchnie
dolną i górną. Jako Ŝe pole elektromagnetyczne jest lokalnie niezmienne, otrzymujemy
y
E
V
y
y
∆
=
,
(2.31)
x
H
I
x
z
∆
−
=
,
(2.32)
z
H
I
z
x
∆
=
,
(2.33)
natomiast impedancje i admitancje moŜna zdefiniować jako
z
y
x
j
Z
x
∆
∆
⋅
∆
=
)
(
ω
ωµ
,
(2.34)
x
z
y
j
Z
z
∆
∆
⋅
∆
=
)
(
ω
ωµ
,
(2.35)
y
z
x
j
Y
∆
∆
⋅
∆
=
)
(
ω
ωε
.
(2.36)
W związku z tym poprzednie równania (2.26)–(2.28) przechodzą do postaci
x
x
y
y
I
Z
z
x
V
z
x
x
V
−
=
−
∆
+
)
,
(
)
,
(
0
0
0
0
,
(2.37)
z
z
y
y
I
Z
z
x
V
z
z
x
V
−
=
−
∆
+
)
,
(
)
,
(
0
0
0
0
,
2.38)
[
] [
]
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
0
0
0
0
0
0
0
0
0
0
z
x
YV
z
x
I
z
x
x
I
z
x
I
z
z
x
I
x
z
z
z
z
−
=
−
∆
+
+
−
∆
+
.
(2.39)
Równania (2.37) i (2.38) opisujące róŜnicę potencjałów między odpowiednio przednią
i tylną lub lewą i prawą ścianą sześcianu są wynikiem natęŜenia prądu wywołanego przez
efektywną impedancję Z
x
lub Z
z
. Równanie (2.39) opisujące róŜnicę potencjałów między
górną i dolną ścianą wynika z natęŜenia prądu pochodzącym od admitancji Y.
Wartości impedancji i admitancji dla zadanej częstości promieniowania
0
ω
ω
=
wyraŜają
się zaleŜnościami
z
y
x
j
Z
x
∆
∆
⋅
∆
=
)
(
0
ω
ωµ
,
(2.40)
x
z
y
j
Z
z
∆
∆
⋅
∆
=
)
(
0
ω
ωµ
,
(2.41)
y
z
x
j
Y
∆
∆
⋅
∆
=
)
(
0
ω
ωε
.
(2.42)
PoniewaŜ na obszarze
∆
x
∆
z komórki elementarnej na przestrzeni jej grubości
∆
y
występuje jednorodne pole elektryczne, moŜemy zastosować analogię do kondensatora
31
okładkowego wypełnionego ośrodkiem o
)
(
0
ω
ε
i
)
(
0
ω
µ
, którego pojemność dana jest
zaleŜnością
y
z
x
C
∆
∆
∆
=
)
(
0
ω
ε
.
(2.43)
Obecność lokalnie jednorodnego pola magnetycznego związana jest z przeciwnie
skierowanymi prądami w równoległych płaszczyznach (Rys. 17), których wkład
do przepływu sumuje się w płaszczyźnie
∆
y
∆
z dla prądów płynących wzdłuŜ
∆
x
oraz w płaszczyźnie
∆
x
∆
y dla prądów płynących wzdłuŜ
∆
z. Z tego względu mamy
do czynienia z indukcjami
z
y
x
L
x
∆
∆
⋅
∆
=
)
(
0
ω
µ
,
(dla kierunku x)
(2.44)
x
z
y
L
z
∆
∆
⋅
∆
=
)
(
0
ω
µ
.
(dla kierunku z)
(2.45)
Na uwagę zasługuje fakt, Ŝe rozproszone pojemność i indukcja są silnie
uzaleŜnione od
)
(
0
ω
ε
i
)
(
0
ω
µ
oraz od wymiarów komórki elementarnej. JeŜeli
rozpatrzymy sieć z komórek elementarnych o pomijalnie małych wymiarach
w porównaniu do długości fali promieniowania zawierających rozproszone L’ i C’, moŜe
ona być traktowana jak ośrodek izotropowy. W związku z tym kaŜdy jednorodny
i bezstratny dielektryk moŜe być (przy zadanej częstotliwości promieniowania) traktowany
jako dyskretna sieć z komórek zawierających tylko cewki i kondensatory natomiast
ośrodek rzeczywisty, obarczony pewnymi stratami transmisyjnymi, modeluje się
uwzględniając w obwodzie szeregową rezystancję.
½
L
½
L
O
ś
Y
Oś
X
Rys. 18
Dwuwymiarowy model linii transmisyjnych opisujący ośrodek o parametrach
0
)
(
µ
ω
µ
=
i
r
ε
ε
ω
ε
0
)
(
=
przy uŜyciu rozproszonej szeregowej induktancji
i rozproszonej równoległej pojemności (na podstawie: [65]).
JeŜeli rozwaŜamy idealnie sześcienne komórki elementarne (Rys. 18), czyli jeŜeli
spełnione są warunki
,
,
z
y
x
d
Z
Z
Z
z
x
∆
=
∆
=
∆
=
=
=
(2.46)
32
wspomniane wyŜej impedancje i admitancja będą wyraŜone zaleŜnościami
d
j
Z
)
(
ω
ωµ
=
,
(2.47)
d
j
Y
)
(
ω
ωε
=
,
(2.48)
zaś efektywne stałe materiałowe dla omawianego modelu linii transmisyjnych będą miały
postacie
d
j
Z
ω
ω
ω
µ
)
(
)
(
=
,
(2.49)
d
j
Y
ω
ω
ω
ε
)
(
)
(
=
.
(2.50)
Dla tradycyjnego ośrodka – izotropowego, niemagnetycznego i charakteryzującego się
względną przenikalnością elektryczną
r
ε
otrzymamy
d
j
Z
0
ωµ
=
,
(2.51)
d
j
Y
r
0
ε
ωε
=
,
(2.52)
czyli szeregowa impedancja i równoległa pojemność mają wartości
d
L
0
µ
=
,
[H]
(2.53)
d
C
0
0
ε
ε
=
.
[F]
(2.54)
Rozproszone wartości L’ i C’ przypadające na jednostkę długości są dodatnie, rzeczywiste
i równe
d
L
L
=
=
0
'
µ
,
(2.55)
d
C
C
r
=
=
0
'
ε
ε
.
(2.56)
Z równania falowego dla obwodów
0
2
2
2
2
2
=
+
∂
∂
+
∂
∂
y
y
y
V
z
V
x
V
β
(2.57)
otrzymać moŜna stałą propagacji
β
równą
φ
ω
ε
ε
µ
ω
ω
β
v
C
L
ZY
r
=
=
=
−
±
=
0
0
'
'
.
(2.58)
PowyŜsza relacja dyspersyjna reprezentuje związek stałej propagacji fali z częstością
promieniowania (Rys. 19).
33
Rys. 19
ZaleŜności dyspersyjne dla ośrodka o parametrach
0
)
(
µ
ω
µ
=
i
r
ε
ε
ω
ε
0
)
(
=
(na podstawie: [65])
Wykres
ω
-
β
ilustruje zmiany stałej propagacji wzdłuŜ wyróŜnionej osi
na płaszczyźnie x-z w funkcji częstotliwości. Dostarcza informacji o wartości fazy
i prędkości grupowej w ośrodku. Wykres
β
x
-
β
z
ilustruje zmiany stałej propagacji
w zaleŜności od kierunku propagacji dla zadanej częstotliwości i przedstawia diagram EFS
(equifrequency surface), czyli powierzchnię równej częstotliwości. Dostarcza informacji
o kierunku fazy i prędkości grupowej w ośrodku. Dla kubicznej, elektrycznie obojętnej
komórki elementarnej taka rozproszona sieć modeluje ośrodek izotropowy i diagram EFS
jest kołowy.
Wartość prędkości fazowej w ośrodku odczytać moŜna z wykresu
ω
-
β
jako odcinek
łączący punkty (0,0) i (
β
0
,
ω
0
) natomiast kierunek z wykresu
β
x
-
β
z
jako linię łączącą
punkty (0,0) i (
β
0z
,
β
0x
). Prędkość fazowa definiowana jest jako
β
ω
φ
=
v
.
(2.59)
Analogicznie odczytujemy z wykresów wartość i zwrot prędkości grupowej zdefiniowanej
jako
1
−
∂
∂
=
ω
β
g
v
.
(2.60)
Z analizy wykresu
ω
-
β
wynika, Ŝe stała propagacji typowego ośrodka
prawoskrętnego zamodelowanego przy uŜyciu rozproszonej sieci szeregowych induktancji
i równoległych pojemności jest zaleŜna od częstotliwości padającej fali EM analogicznie
jak ma to miejsce w naturalnych dielektrykach przy niskich częstotliwościach. Na wykres
β
x
-
β
z
widać, Ŝe prędkości fazowa i grupowa (dla ośrodka bezdyspersyjnego) są równe oraz
równoległe i wyraŜają się następująco
g
r
v
C
L
v
=
∂
∂
=
=
=
=
−
1
0
0
1
'
'
1
ω
β
ε
ε
µ
β
ω
φ
.
(2.61)
34
Obie prędkości w tym przypadku są dodatnie, co jest wynikiem wyboru dodatniego
pierwiastka w równaniu (2.61) odpowiadającego jednej z gałęzi wykresu
ω
-
β
,
więc współczynnik załamania, który definiowany jest jako stosunek prędkości światła
w próŜni do prędkości fazowej w ośrodku jest dodatni
r
r
C
L
v
c
n
ε
ε
µ
ε
ε
µ
ε
µ
φ
=
=
=
=
0
0
0
0
0
0
'
'
.
(2.62)
Co więcej – tak jak było to oczekiwane – impedancja falowa ośrodka
η
jest dokładnie
równa uŜytej impedancji w rozproszonym modelu linii transmisyjnych
0
0
0
'
'
Z
C
L
r
=
=
=
ε
ε
µ
η
.
(2.63)
Symulacja ujemnego współczynnika załamania
Dalszym etapem stała się próba modelowania za pomocą komórek elementarnych LC
sztucznych ośrodków dielektrycznych – metamateriałów. Fakt, iŜ przy uŜyciu modelu linii
transmisyjnych o określonej wartości rozproszonych induktancji L’ i pojemności C’ mamy
wpływ na parametry wytwarzanego ośrodka, pozwala uzyskać materiały nie występujące
w naturze. PoniewaŜ L’ i C’ są związane z wartościami przenikalności elektrycznej
i magnetycznej, więc wprowadzenie do komórki elementarnej przedstawionej na Rys. 18
ujemnych wartości L’ i C’ pozwala uzyskać ujemny współczynnik załamania n. Z natury
tych elementów wynika jednak, Ŝe gdy szeregowa impedancja
d
L
j
'
ω
−
i równoległa
admitancja
d
C
j
'
ω
−
mają wartości ujemne, zmienia się ich rola w układzie. W takim
wypadku szeregowa impedancja odgrywa rolę szeregowej admitancji, zaś równoległa
admitancja – równoległej impedancji. Tak zbudowana elementarna komórka
przedstawiona jest na Rys. 20.
O
ś
Y
Oś
X
Rys. 20
Dwuwymiarowy model linii transmisyjnych opisujący ośrodek o jednocześnie
ujemnych, dyspersyjnych parametrach
)
(
ω
µ
µ
−
=
i
)
(
ω
ε
ε
−
=
przy uŜyciu
rozproszonej szeregowej induktancji i rozproszonej równoległej pojemności
(na podstawie: [65]).
35
TakŜe w tym przypadku wymiary komórki elementarnej mają kluczowy wpływ
na zaleŜności dyspersyjne parametrów uzyskanego w ten sposób ośrodka [66], [67].
Korzystając z równań (2.49) i (2.50) efektywną przenikalność elektryczną i magnetyczną
takiego modelu wyznaczyć moŜna jako
( )
d
C
j
Cd
j
2
1
1
ω
ω
ω
ω
µ
−
=
=
,
(2.64)
( )
d
L
j
Ld
j
2
1
1
ω
ω
ω
ω
ε
−
=
=
.
(2.65)
W odróŜnieniu od modelu linii transmisyjnych dla materiałów dodatnich, tutaj parametry
ośrodka są wyraźnie ujemne, a dodatkowo mają charakter dyspersyjny. Uśredniona
po czasie energia elektryczna i magnetyczna zgromadzona w takim materiale jest jednak
dodatnia, więc zasada zachowania energii pozostaje spełniona.
Rozproszone wartości indukcji i pojemności wynoszą dla tego przypadku odpowiednio
d
C
C
⋅
=
'
,
m]
[F
⋅
(2.66)
d
L
L
⋅
=
'
.
m]
[H
⋅
(2.67)
Cechą charakterystyczną stałej propagacji w powyŜszym modelu jest odwrotny związek
z częstością
'
'
1
C
L
ZY
ω
β
−
=
−
−
=
(2.68)
zaś odpowiadające temu wykresy
ω
-
β
oraz
β
x
-
β
z
przedstawione są na Rys. 21.
Rys. 21
ZaleŜności dyspersyjne dla ośrodka o parametrach
)
(
ω
µ
µ
−
=
i
)
(
ω
ε
ε
−
=
(na podstawie: [65])
36
W tym przypadku prędkości fazowa i grupowa są antyrównoległe
g
v
C
L
v
−
=
∂
∂
−
=
−
=
=
−
1
2
'
'
ω
β
ω
β
ω
φ
.
(2.69)
Ze względu na to, Ŝe prędkości grupowe w modelu linii transmisyjnych ośrodka
dodatniego i ujemnego są antyrównoległe, moŜna spodziewać się ujemnego załamania
fali EM przy przejściu przez granicę tych ośrodków.
Współczynnik załamania takiego ośrodka
( ) ( )
0
0
2
0
0
'
'
1
ε
µ
ω
ε
µ
ω
ε
ω
µ
φ
⋅
⋅
−
=
−
=
=
C
L
v
c
n
.
(2.70)
Efektywna impedancja fali
( )
( )
0
'
'
Z
C
L
r
=
=
=
ω
ε
ω
µ
η
.
(2.71)
Opisany tu model został zweryfikowany doświadczalnie poprzez zastosowanie
go do badania zdolności skupiających zbudowanej na jego podstawie płaskiej soczewki
o ujemnym współczynniku załamania [67]. Model znalazł teŜ kilka innych zastosowań
w układach i urządzeniach optycznych [68]
−
[70]. Warto zauwaŜyć, Ŝe takŜe ten model jest
bezstratny, co oznacza, iŜ model linii transmisyjnych pozwala na osiągnięcie ujemnego
współczynnika załamania bez konieczności wprowadzania do układu elementów
rezystywnych.
W modelu ośrodka ALMW/SRR (rozdział II.4) efektywna funkcja przenikalności
przyjmowała ujemne wartości w zakresie częstotliwości między
0
ω
a
mp
ω
, który
w modelu linii transmisyjnych (Rys. 15) odpowiada dokładnie zakresowi częstotliwości,
w którym czynna jest szeregowa gałąź. Analogicznie dla przenikalności elektrycznej
zakres częstotliwości dający ujemną jej wartość odpowiadał zakresowi pracy gałęzi
równoległej w modelu linii transmisyjnych (Rys. 15). W ujęciu linii transmisyjnych
materiał ALMW/SRR wprowadza dodatkowe składowe rezonansowe, co moŜe zostać
wyeliminowane poprzez bezpośrednie zastosowanie indukcji i pojemności. Dzięki temu
moŜliwe jest osiągnięcie bardzo szerokich pasm częstotliwości, w których współczynnik
załamania utrzymuje swoją ujemną wartość.
37
II.6 Nanostruktury z drutów metalicznych
Metamateriał będący kompozytem złoŜonym z par nanodrutów po raz pierwszy opisał
i zasymulował numerycznie Viktor A. Podolskiy [71], [72]. Udowodnił, Ŝe zewnętrzne
pole EM oddziałujące na pojedyncze nanodruty moŜe wzmacniać powstający
na powierzchni materiału plazmon powierzchniowy, w wyniku czego lokalne pola EM
mogą zostać wzrosnąć nawet o trzy rzędy wielkości. Ponadto w kompozycie złoŜonym
z par równoległych nanodrutów występować moŜe rezonans elektryczny i magnetyczny,
co przy
spełnionych
odpowiednich
warunkach,
moŜe
skutkować
ujemnym
współczynnikiem załamania takiego metamateriału.
W odróŜnieniu od metamateriału złoŜonego z cienkich drutów metalicznych
i rozproszonych rezonatorów kołowych, które musiały być projektowane dla konkretnej
częstości padającego promieniowania, zasadniczą zaletą metamateriału z nanodrutów jest
o wiele większa jego uniwersalność. Wzajemne oddziaływania nanodrutów tworzących
kompozyt powoduje powstanie wielu zlokalizowanych plazmonów powierzchniowych
o róŜnych
częstościach
rezonansowych.
Plazmon
powierzchniowy
pochodzący
od pojedynczego nanodrutu jest stosunkowo wąski (50 nm), jeŜeli jednak nanodruty
tworzą kompozyt, zakres częstości rezonansowych jest bardzo duŜy, co umoŜliwia
wykorzystanie takiego materiału w bardzo szerokim przedziale częstości ściśle zaleŜnym
od wymiarów nanodrutów.
Komórka elementarna metamateriału zaproponowanego przez Podolskiy’ego
składała się z pary równoległych nanodrutów metalicznych o długości 2b
1
, średnicy b
2
,
rozsuniętych na odległość d (Rys. 22).
Rys. 22
(a) para równoległych nanodrutów; (b) dwuwymiarowy kompozyt zawierający
pary równoległych nanodrutów (na podstawie [71]).
38
Grupa Podolskiy’ego rozwaŜała przypadek, gdy
1
2
b
d
b
<<
<<
. Promień drutu b
2
musi być duŜo mniejszy od długości fali padającego promieniowania i moŜe być
porównywalny z głębokością penetracji promieniowania w strukturę. Długość drutu 2b
1
moŜe być rzędu długości fali świetlnej. Ze względu na takie wymiary nanodrutów bardzo
trudno było jednoznacznie wyznaczyć rozwiązania równań Maxwella w tym układzie.
Podolskiy zastosował dyskretne przybliŜenie dipolowe [73] zastępując nanodruty duŜą
liczbą zachodzących na siebie, spolaryzowanych sfer ułoŜonych w kubiczną sieć (Rys. 23).
Rys. 23
Model pojedynczego nanodrutu w dyskretnym przybliŜeniu dipolowym.
JeŜeli stała takiej sieci a i promień R pojedynczej sfery są duŜo mniejsze
od długości fali padającego światła, pola EM pochodzące od kaŜdego z dipoli mogą być
traktowane jako kwazistatyczne i są sumą pól zewnętrznych w danym miejscu sieci oraz
pól pochodzących od wszystkich pozostałych dipoli. Podolski udowodnił, Ŝe wyniki
symulacji przy uŜyciu tego przybliŜenia silnie zaleŜą od przyjętego współczynnika
pokrywania się sąsiadujących sfer (a/R).
Padająca fala EM rozchodziła się tak, Ŝe pole elektryczne było równoległe
do osi nanodrutów, podczas gdy pole magnetyczne było do nich prostopadłe. Składowa
magnetyczna padającego promieniowania EM indukuje w drutach przepływ prądu,
będącego sumą prądów płynących w kaŜdym z nanodrutów oraz prądu przesunięcia
pomiędzy drutami. Wartość odpowiedzi rezonansowej silnie zaleŜy od wymiarów drutów
i maleje wraz ze wzrostem promienia drutu lub odległości d między drutami, gdyŜ wtedy
straty stają się znaczące. Efekt jest najmocniejszy dla drutów o promieniach b
2
porównywalnych z głębokością penetracji promieniowania EM (~20nm).
Aby uwzględnić efekt naskórkowy wprowadzono funkcję
( ) ( )
( )
[
]
( )
[
]
∆
+
∆
+
⋅
∆
−
=
∆
i
j
i
j
i
f
1
0
1
1
1
r
r
,
(2.72)
gdzie
j
r
– wektor gęstości prądu,
1
2
2
>>
=
∆
c
b
m
ω
πσ
parametr opisujący stosunek
ś
rednicy nanodrutów do głębokości penetracji promieniowania,
m
σ
przewodność metalu.
39
Przenikalność magnetyczną
µ
kompozytu wyznaczono posługując się momentem
magnetycznym
H
m
dla pojedynczej pary drutów przybliŜonej dwoma równoległymi,
nieskończonymi drutami, dzięki czemu moŜliwe było zastosowanie równania telegrafistów
(Dodatek C). Moment magnetyczny został wyznaczony jako
( )
( )
( )
3
1
1
1
2
2
3
1
tan
2
gb
gb
gb
kd
C
Hb
m
H
−
=
,
(2.73)
gdzie
=
2
2
ln
4
b
d
C
d
ε
pojemność na jednostkę długości,
( )
∆
∆
+
=
2
2
ln
2
b
d
f
i
k
g
d
d
ε
ε
.
Padające pole elektryczne jest równoległe do osi drutów, więc wzbudza w nich
równe prądy, które mogą być traktowane jako niezaleŜne. Całkowity moment dipolowy
dla pary nanodrutów dany jest jako
( )
( )
Ω
+
∆
+
∆
=
cos
1
ln
1
1
3
2
2
1
2
1
2
2
2
1
b
b
b
b
f
E
f
b
b
d
d
m
m
E
ε
ε
ε
,
(2.74)
gdzie
( )
+
+
=
Ω
2
1
1
2
1
2
1
2
1
ln
ln
b
b
kb
i
b
b
k
b
d
d
ε
ε
to bezwymiarowa częstość promieniowania.
Podolskiy wyznaczył efektywną stałą dielektryczną i przenikalność magnetyczną
struktury przedstawionej na Rys. 22 jako
E
d
d
b
b
p
E
2
1
4
1
+
=
ε
,
H
m
d
b
b
p
H
2
1
4
1
+
=
µ
,
(2.75)
gdzie p – gęstość powierzchniowa metalu.
Obszar, w którym taki ośrodek wykazuje ujemną wartość współczynnika załamania
uwarunkowany jest wymiarami nanodrutów, czyli wartościami parametrów b
1
, b
2
oraz d.
Ustalając odpowiednio wartości tych parametrów moŜna przesunąć aktywne pasmo
częstotliwości do zakresu optycznego.
Alternatywnym
sposobem
na
wytworzenie
metamateriału
o
ujemnym
współczynniku załamania jest zastosowanie stosu warstw z nadrukowanymi pojedynczymi
nanodrutami lub uŜycie jednej takiej warstwy umieszczonej nad powierzchnią metalu.
W tym drugim wypadku promieniowanie EM powoduje powstanie obrazu nanodrutów
na powierzchni metalu, co prowadzi do powstania wirtualnej pary równoległych
nanodrutów. Powierzchnię pomiędzy metalem i filmem z nanodrutami wypełnia się
dielektrykiem. Zmieniając grubość i parametry dielektryka moŜna regulować odległość d
pomiędzy nanodrutami.
40
II.7 Metamateriały dla zakresu widzialnego
PoniewaŜ materiały charakteryzujące się ujemnym współczynnikiem załamania nie istnieją
w naturze, w ostatnich latach nieustająco trwają prace badawcze nad wytworzeniem
sztucznych materiałów dających ujemną odpowiedź elektryczną i magnetyczną w pewnym
przedziale częstotliwości. Dwa główne podejścia do tematu realizacji zjawiska ujemnego
załamania obejmują zastosowanie metamateriałów oraz kryształów fotonicznych.
W zakresie promieniowania mikrofalowego ujemny współczynnik załamania został
uzyskany w obu typach struktur, podczas gdy dla zakresu widzialnego ostatnio
zastosowanie znalazły kryształy fotoniczne [74], [75]. W zakresie 10 – 100 THz
zaprezentowane niedawno zostały rezonujące struktury magnetyczne charakteryzujące się
ujemną przenikalnością magnetyczną [76]
−
[78]. JednakŜe, pomijając prace teoretyczne
i symulacje numeryczne [79], ujemne załamanie fal EM dla zakresu widzialnego
potwierdzone doświadczalnie zostało dopiero pod koniec 2005 roku [80].
MoŜliwość zastosowania nanometrowych struktur do uzyskania ujemnego
współczynnika załamania dla zakresu widzialnego opisano teoretycznie [81]
−
[83]
w 2002 roku. Pokazano, Ŝe pary równoległych nanodrutów lub płytek mogłyby
z powodzeniem zastąpić wykorzystywane dotąd rozszczepione rezonatory kołowe (SRR).
Oczekiwano, Ŝe zbudowana z nich struktura będzie charakteryzowała się ujemnym
współczynnikiem załamania nawet, jeŜeli komórka elementarna nie będzie zawierała
dodatkowego metalowego drutu. Analogia pomiędzy typowym SRR i parą nanodrutów
przedstawiona została na Rys. 24 i łatwo ją zrozumieć, jeŜeli potraktujemy SRR
przedstawiony na Rys. 24 (a) jako obwód LC o induktancji L i pojemności przerwy
pierścienia C.
Rys. 24
Analogia pomiędzy rozszczepionym pierścieniem (SRR) a parą nanodrutów:
kolejne fazy przejścia pomiędzy strukturami (a)-(d) [84].
Zmienne pole magnetyczne prostopadłe do powierzchni SRR indukuje w nim
przepływ prądu, który w pobliŜu częstości rezonansowej
LC
LC
1
=
ω
wywołuje
moment magnetyczny prostopadły do płaszczyzny SRR i przeciwstawiający się
zewnętrznemu polu magnetycznemu, co powoduje pojawienie się ujemnej przenikalności
magnetycznej
µ
. JeŜeli przerwa w pierścieniu zostanie powiększona (Rys. 24 (b)),
pojemność C zmniejszy się, co spowoduje wzrost częstości
ω
LC
potrzebnej do wzbudzenia
rezonansu w obwodzie. Gdy drut tworzący pierścień SRR zostanie przerwany równieŜ
z drugiej strony (Rys. 24 (c)), w obwodzie pojawi się druga szeregowa pojemność C,
co dodatkowo zmniejszy wypadkową pojemność obwodu LC. Kolejnym krokiem jest
całkowite otwarcie pierścienia SRR z obu stron (Rys. 24 (d)), prowadzące do otrzymania
pary równoległych drutów. Spadek pojemności C prowadzi do wyŜszej częstości
41
rezonansowej zgodnie ze wzorem
LC
LC
1
=
ω
, co oznacza, Ŝe długość fali
λ
padającego promieniowania EM potrzebnego do wywołania rezonansu w takim obwodzie
jest coraz mniejsza i moŜe osiągnąć nawet wartości z zakresu promieniowania
widzialnego. Zaletą par równoległych nanodrutów w stosunku do wcześniej stosowanych
komponentów są mniejsze wymagania konstrukcyjne pod względem stosunku wymiarów
komórki elementarnej i długości fali
λ
padającego promieniowania. Dla SRR stosunek ten
musiał być co najmniej rzędu
10
=
a
λ
, aby strukturę moŜna było traktować jak
jednorodną, zaś dla par nanodrutów jest on prawie 5 razy mniejszy (
2
≈
a
λ
).
(a)
(b)
(c)
(d)
Rys. 25
(a) Wymiary pary równoległych nanodrutów tworzących komórkę elementarną;
(b) Schemat wzajemnego rozmieszczenia par nanodrutów w próbce; Wykonane
pod mikroskopem elektronowym zdjęcie (c) dwuwymiarowej tablicy i (d) pojedynczej pary
nanodrutów. (na podstawie: [85]).
W roku 2005 grupa badawcza z Purdue University [80] zaprezentowała
metamateriał zbudowany z par równoległych nanodrutów w kształcie trapezu o róŜnicy
szerokości podstaw 20 nm. Opisana przez nich periodyczna dwuwymiarowa tablica
z nadrukowanymi parami złotych drutów była stosunkowo prosta do wytworzenia
w nanoskali i wskazała kierunek dalszych prac nad zaprojektowaniem metamateriału
dla zakresu widzialnego. Nanodruty wytworzone zostały przy uŜyciu litografii
42
elektronowej na szklanym podłoŜu pokrytym cienką warstwą Cr, która zapobiegała
naładowaniu się podłoŜa i została usunięta po zakończeniu procesu. Nanoszenie warstw
powtórzono kilkukrotnie kolejno dla warstw Ti, Au i SiO
2
, aŜ do uzyskania struktury
przedstawionej na Rys. 25.
Dwa równoległe druty tworzyły obwód, który zachowywał się jak linia
transmisyjna (Rys. 26). Prądy przesunięcia na końcach drutów powodowały „zamknięcie”
obwodu. JeŜeli na próbkę padało światło spolaryzowane poprzecznie magnetycznie
względem osi nanodrutów, wzbudzany był jednocześnie rezonans elektryczny
i magnetyczny,
co powodowało
powstanie
ujemnej
odpowiedzi elektrycznej
ε
i magnetycznej
µ
w tak zaprojektowanym metamateriale. Zakres częstotliwości fali
elektromagnetycznej, dla której zjawisko to było moŜliwe uzaleŜniony był od wymiarów
nanodrutów oraz ich wzajemnego rozmieszczenia w tablicy. Współczynnik załamania
w takim materiale powyŜej częstości rezonansowej moŜe przyjąć wartości ujemne [71],
[72]. Rezonans plazmowy moŜe być traktowany jak rezonans w obwodzie LC
o induktancji L, zapewnionej przez metalowe druty, i pojemności C zapewnionej przez
przestrzeń pomiędzy nimi.
Rys. 26
(a) UłoŜenie nanodrutów w metamateriale; (b) Komórka elementarna zawierająca
parę równoległych nanodrutów (na podstawie: [85]).
Doniesienie literaturowe na temat dalszej miniaturyzacji metamateriałów
i przesunięcia operacyjnej częstości w spektrum widzialne pochodzą z sierpnia 2006 roku
[88]. Grupa naukowców G.Dolling, M.Wegener, C.M.Soukoulis oraz S.Linden
wytworzyła metamateriał zbudowany z par równoległych srebrnych nanodrutów
nadrukowanych na szklanym podłoŜu, charakteryzujący się ujemnym współczynnikiem
załamania dla światła o długości fali 780 nm (zob. Rys. 27). Dla zakresu widzialnego
srebro charakteryzuje się znacząco mniejszymi stratami w porównaniu do złota [89],
co pozwoliło po raz pierwszy zademonstrować ujemne załamanie światła dla światła
czerwonego.
43
(c)
Rys. 27
(a) Schemat budowy metamateriału dla zakresu widzialnego; (b) Komórka
elementarna o wymiarach a
x
= a
y
= 300 nm, w
x
= 102 nm, w
y
= 68 nm, t = 40 nm,
s = 17 nm, e
x
= e
y
= e = 8 nm; (c) Zdjęcie pojedynczej pary nanodrutów i
dwuwymiarowej tablicy wykonane pod mikroskopem elektronowym (na
podstawie: [88]).
Stosunek grubości kompozytu 2t+s = 97 nm do szerokości drutu w
y
= 68 nm przekracza 1,
co stwarzało pewne trudności technologiczne. Grubość drutów zorientowanych wzdłuŜ
składowej elektrycznej padającego promieniowania musiała zostać zwiększona
do w
x
= 102 nm. Pozwoliło to na podniesienie efektywnej częstości plazmowej powyŜej
częstości operacyjnej, co było konieczne dla osiągnięcia ujemnej przenikalności
elektrycznej
ε
. Do wytworzenia metamateriału wykorzystano fluorek magnezu MgF
2
,
którego współczynnik załamania wynosi
38
,
1
2
=
MgF
n
oraz szkło o współczynniku
5
,
1
=
szkłz
n
.
44
III Eksperymenty
Naturalną koleją rzeczy była chęć eksperymentalnej weryfikacji zjawiska ujemnego
załamania fali elektromagnetycznej we wzbudzających tyle zainteresowania ośrodkach.
Jako pierwsi dokonali tego R.A.Shelby, D.R.Smith i S.Schultz [54] w 2001 roku. Pomimo
wątpliwości niektórych [55], eksperyment grupy Shelby’ego stał się początkiem nowego
rozdziału w historii metamateriałów. Po okresie teoretycznych dyskusji na temat odkrycia
Veselago [1], metamateriały potwierdziły swą wyjątkowość otwierając tym samym drogę
do licznych nowych zastosowań (rozdział III.1).
Kolejny waŜny eksperyment wykonany został w 2003 roku przez grupę
I.V. Shadrivov’a [90]. Udowodniono, Ŝe wprowadzenie defektów do warstwowej struktury
złoŜonej z naprzemiennie ułoŜonych warstw ośrodka o ujemnym współczynniku załamania
i powietrza pozwala na modyfikowanie zakresu częstości, dla którego obserwuje się
ujemne załamanie fali elektromagnetycznej (rozdział III.2) Istnienie fal wstecznych
w ujemnych metamateriałach zbadał G.V.Eleftheriades [102]. Przeprowadzony przez
niego eksperyment udowodnił, Ŝe model metamateriału wykorzystujący linie transmisyjne
(rozdział II.5) pozwala na łatwe wykorzystanie zjawiska ujemnego załamania do budowy
wydajniejszych urządzeń mikrofalowych (rozdział III.3). Długo oczekiwane ujemne
załamanie dla fal z zakresu widzialnego zrealizowane w praktyce zostało dopiero w 2006
roku [88].
III.1 Pierwsze dowody eksperymentalne
Przeprowadzony w 2001 roku eksperyment [54] dotyczył rozpraszania fali EM
o częstotliwości z zakresu mikrofalowego na próbce metamateriału charakteryzującego się
pasmem częstotliwości, w którym jego efektywny współczynnik załamania n ma wartość
mniejszą od zera. Mierząc kąt odbicia promienia przechodzącego przez pryzmat wykonany
z takiego metamateriału, określano efektywny współczynnik załamania n zgodny
z prawem Snelliusa [10]. Eksperyment potwierdził przewidywania równań Maxwella,
iŜ współczynnik załamania n jest równy ujemnemu pierwiastkowi kwadratowemu
z iloczynu przenikalności
ε
i
µ
dla zakresu częstotliwości, gdzie obie przenikalności
są mniejsze od zera (zob. równanie (1.5)). W doświadczeniu wykorzystano metamateriał
zaproponowany przez D.R.Smitha [32]
−
[35].
Próbka zbudowana była przy wykorzystaniu dwuwymiarowej okresowej tablicy
zawierającej miedziane rozproszone rezonatory kołowe (SRR) i druty metaliczne
(ALMW), wyprodukowanej metodą maski i kwasorytu na podłoŜu krzemowym o grubości
0,25mm. Metamateriał uŜyty w tym doświadczeniu zbudowany był z komórek
elementarnych o wymiarach 5mm, co oznacza, Ŝe dla badanego zakresu częstotliwości
(8-12 GHz) średnia długośc fali promieniowania
λ
=3cm była 6 razy większa od wymiarów
takiej komórki i materiał mógł być traktowany jako jednorodny. Kąt pomiędzy prostymi
normalnymi do powierzchni padania (nr 1 na Rys. 28) i załamania (nr 2 na Rys. 28) był
równy 18,43
o
. Aby określić współczynnik załamania, mierzono odchylenie wiązki
promieniowania mikrofalowego przechodzącej przez pryzmat wykonany z takiego
metamateriału. Próbka i absorbent promieniowania mikrofalowego umieszczone były
pomiędzy dwoma kołowymi płytami aluminiowymi o promieniu 15cm. Na obwodzie
wierzchniej płyty zamontowany był ruchomy detektor promieniowania mikrofalowego
pozwalający na pomiar transmitowanej mocy pod dowolnym kątem. Na ścianę pryzmatu,
45
stanowiąca powierzchnię padania (nr 1), padała wiązka promieniowania mikrofalowego
o poprzecznej polaryzacji magnetycznej. Po przejściu przez pryzmat fala doznawała
ugięcia zgodnego z prawem Snelliusa na płaszczyźnie załamania (nr 2). Aby zmniejszyć
spowodowane dyfrakcją rozmycie kątowe padającego promieniowania, poprowadzone ono
było współosiowym kablem do adaptera falowego, a następnie przechodziło pomiędzy
dwoma aluminiowymi płaszczyznami rozsuniętymi na taką samą odległość jak aluminiowe
koła (1,2 cm), gdzie dodatkowo po bokach ograniczały je rozsunięte na odległość 9,3 cm
płaszczyzny absorbentu promieniowania mikrofalowego. Strzałki na Rys. 28 przedstawiają
bieg fali przy przejściu przez teflonową próbkę kontrolną – kąt
θ
ma tu wartość dodatnią.
Detektor obracany był po obwodzie z krokiem 1,5
o
, a przechodzące widmo mocy mierzone
w funkcji kąta obrotu
θ
detektora od normalnej do powierzchni padania próbki.
Rys. 28
Schemat układu pomiarowego [54].
Pomiary wykonano dla ośmiu róŜnych połoŜeń próbki w układzie. Uśrednione
wyniki przedstawiono na wykresach w funkcji kąta (Rys. 29) (po wcześniejszym
znormalizowaniu mocy transmitowanego przez próbkę promieniowania) oraz w funkcji
częstotliwości (Rys. 30). Maksimum mocy transmitowanej występowało przy bardzo
zbliŜonych kątach dla kaŜdego z połoŜeń próbki.
M
o
c
tr
a
n
sm
it
o
w
a
n
a
Rys. 29
Moc transmitowana jako funkcja kąta [54].
46
Przy częstotliwości 10,5 GHz, mikrofale były uginane pod dodatnimi kątami
w przypadku teflonu i w przeciwną stronę dla próbki z z metamateriału. Maksimum
zarejestrowanej mocy promieniowania dla próbki teflonowej miało miejsce przy kącie
θ
powietrze
= 27º, co odpowiada wartości współczynnika załamania n
teflon
= 1,4
±
0,1 i mieści
się w granicach błędu. Natomiast dla próbki z metamateriału zmierzony kąt wyjściowy
θ
powietrze
=
−
61
o
, przy którym moc osiągała maksimum oznacza, Ŝe współczynnik
załamania tego ośrodka to n
LHM
=
−
(2,7
±
0,1) . Na Rys. 30 przedstawiony został
zmierzony współczynnik w funkcji częstotliwości dla ujemnej próbki (ciągła czarna linia)
i porównany z krzywą teoretyczną (ciągła czerwona linia) oraz wykresem dla teflonu
(ciągła niebieska linia)
LHM
Teflon
LHM (teoria)
3
2
1
0
-1
-2
-3
Częstotliwość [GHz]
8
9
10
11
12
Rys. 30
Współczynnik załamania jako funkcja częstotliwości [54].
Oczekiwano, Ŝe współczynnik załamania będzie szybko dąŜył do wartości mocno
ujemnych przy niskoczęstotliwościowej granicy w badanego obszaru „lewoskrętnego”
(10,2
÷
10,8
GHz),
by
następnie
osiągnąć
wartość
zerową
dla
granicy
wysokoczęstotliwościowej. Wyniki częściowo potwierdziły rozwaŜania teoretyczne,
jednakŜe w pewnych obszarach (przerywana czarna linia) współczynnik załamania n był
albo niemoŜliwy do wykrycia albo zdominowany przez składową urojoną (przerywana
czerwona linia) i nie mógł być wiarygodnie zbadany doświadczalnie. W badanym paśmie
częstotliwości zmierzony współczynnik załamania dla teflonu (linia niebieska na Rys. 30)
jest w przybliŜeniu wartością stałą, zaś współczynnik załamania dla badanego
metamateriału jest ujemny i ma charakter silnie dyspersyjny, co zgadza się z teorią.
Do
wyznaczenia
teoretycznej
zaleŜności
współczynnika
załamania
n
od częstotliwości uŜyto poniŜszych ogólnych równań na zaleŜność przenikalności
ε
i
µ
od częstotliwości
γω
ω
ω
ω
ω
µ
ω
µ
i
mo
mo
mp
+
−
−
−
=
2
2
2
2
0
1
)
(
,
(3.1)
γω
ω
ω
ω
ω
ε
ω
ε
i
eo
eo
ep
+
−
−
−
=
2
2
2
2
0
1
)
(
,
(3.2)
47
przy wykorzystaniu parametrów f
mp
= 10,95 GHz, f
mo
= 10,05 GHz, f
ep
= 12,8 GHz,
f
eo
= 10,3 GHz,
γ
= 10 MHz gdzie
ω
mo
,
ω
mp
– magnetyczna częstość rezonansowa
i plazmowa,
ω
eo
,
ω
ep
– elektryczna częstość rezonansowa i plazmowa,
1
−
=
i
, f=
ω
/ 2
π
.
Istnieją dwa zasadnicze ograniczenia w tym doświadczeniu, które uniemoŜliwiają
zmierzenie efektywnego współczynnika załamania odpowiadającego krawędziom
„lewoskrętnego” pasma częstotliwości. Po pierwsze, kiedy efektywny współczynnik
załamania osiąga zero, długość fali w ośrodku ujemnym staje się bardzo duŜa,
przypuszczalnie większa niŜ wymiary próbki. Nie udało się jednoznacznie określić
współczynnika załamania w zakresie od 10,8GHz do 12GHz, czyli powyŜej
„lewoskrętnego” pasma częstotliwości; uzyskane wyniki odpowiadają raczej urojonej
składowej współczynnika, a nie jego dodatniej wartości. To ograniczenie moŜe zostać
częściowo wyeliminowane poprzez zastosowanie grubszych i szerszych próbek. Po drugie,
poniewaŜ obszar ugięcia jest zawęŜony do około 18,4
o
od normalnej, gdy
3
≥
n
,
promieniowanie ulega raczej całkowitemu wewnętrznemu odbiciu niŜ załamaniu,
co według autorów eksperymentu mogłoby tłumaczyć, dlaczego nie udało się zmierzyć
współczynników załamania o wartościach poniŜej -3 i powyŜej 3.
Ujemne załamanie – fikcja czy rzeczywisto
ść
?
Wyniki eksperymentu [54] zakwestionowane zostały przez N.Garcia i M.Nieto-
Vesperinas’a [55], P.M.Valanju [56] i innych [57], co zaowocowało oŜywioną dyskusją
w środowisku naukowym na temat zjawiska ujemnego załamania i istnienia ośrodków,
w których byłoby ono moŜliwe. PodwaŜona została moŜliwość ujemnego załamania
rzeczywistego promieniowania zawierającego pewne spektrum długości fali [56]. Ujemny
współczynnik załamania występować moŜe jedynie w materiałach dyfrakcyjnych,
co oznacza, Ŝe jest uzaleŜniony od częstości padającego promieniowania, czyli takŜe
od jego długości fali. W związku z tym przy przejściu przez metamateriał wiązka
rzeczywistego promieniowania zostałaby raczej rozszczepiona niŜ ujemnie załamana.
Skrytykowano takŜe sposób przeprowadzenia eksperymentu. Wykorzystany przez
R.A.Shelby’ego, D.R.Smith’a i S.Schultz’a ośrodek [54] charakteryzował się znacznymi
stratami, zaś detektor ustawiony był w niewielkiej odległości od powierzchni pryzmatu,
co było uwaŜane za przyczynę złudzenia, Ŝe fala EM została ugięta pod ujemnym kątem.
Twierdzono [56], Ŝe gdyby detektor znajdował się w większej odległości od próbki, Ŝadne
promieniowanie nie zostałoby zmierzone. Krytyka wyników eksperymentu grupy
Shelby’ego zmobilizowała zwolenników idei ujemnego załamania do szukania
kontrargumentów. Niecały rok później opublikowane zostały uzyskane niezaleŜnie przez
dwie grupy badawcze [58]-[60] wyniki, wyraźnie potwierdzające, przy zastosowaniu
prawa Snelliusa, istnienie ujemnego załamania promieniowania z zakresu mikrofalowego.
A.Houck i in. z MIT
16
[58] wykonali dwa doświadczenia, w których odwzorowali
mikrofalowe pole EM przechodzące przez kilka pryzmatów wykonanych z materiałów
dodatnich i ujemnych oraz płaską płytkę z materiału ujemnego.
Materiałem odniesienia
podobnie jak w [54] był teflon, zaś metamateriał ujemny był kompozytem analogicznym
do opisanego w [35]
17
. Aby ustalić zakres częstotliwości promieniowania, w którym
metamateriał charakteryzuje się ujemnym współczynnikiem załamania, wykonane zostały
16
Massachusetts Institute of Technology
17
Metamateriał uŜywany w doświadczeniach zbudowany był z miedzianych drutów o grubości 50
µ
m,
umieszczonych na podłoŜu o grubości 0,5 mm zestawionych w sieć 6 mm komórek elementarnych.
48
niezaleŜne pomiary właściwości elektromagnetycznych dla struktur zawierających tylko
prostoliniowe przewodniki oraz tylko przerwane pierścienie. Operacyjna częstotliwość
ustalona została na 10,5 GHz. Ponadto udowodniono, Ŝe górna płaszczyzna prowadnicy
falowej nie moŜe dotykać powierzchni próbki, gdyŜ nie obserwowano wtedy zachowania
ujemnego. Przyczyny tego zjawiska upatrywano w zmianę pojemności metalowych
przewodników, gdy odległość płaszczyzny od próbki była mniejsza niŜ 1,5 mm.
W ramach pierwszego eksperymentu A.Houck i in. zaproponowali nowy układ
pomiarowy, wykorzystujący zamkniętą i dokładnie izolowaną prowadnicę falową,
zbudowaną z dwóch równoległych metalowych płaszczyzn (górna płaszczyzna była
dwukrotnie większa od dolnej), w całości pokrytych absorbentem promieniowania
mikrofalowego i mogących przesuwać się względem siebie. UmoŜliwiło to pomiar
rozkładu pola EM w całej przestrzeni wewnątrz prowadnicy. Pomiary wykonane zostały
dla kilku pryzmatów dodatnich (teflon) i ujemnych, o dwóch róŜnych kątach, dla dwóch
połoŜeń górnej płaszczyzny nad próbką. Zmierzona w ten sposób wartość współczynnika
załamania dla teflonu pokrywała się z wartością rzeczywistą
18
, co pozwoliło zweryfikować
poprawność układu pomiarowego. Kąty załamania fali EM obserwowane dla próbki
ujemnej, były wyraźnie umieszczone po przeciwnej stronie normalnej do płaszczyzny
załamania pryzmatu w porównaniu do pomiarów dla próbki teflonowej
19
. Na uwagę
zasługuje fakt, Ŝe w przypadku próbki ujemnej, tłumienie fali miało wartość 25dB
(w porównaniu do 5dB dla pryzmatu z teflonu), co wskazuje na duŜe znaczenie efektów
rozpraszania i absorpcji w metamateriałach.
Drugi eksperyment wykonany przez grupę z MIT dotyczył opisanego wcześniej
teoretycznie przez Pendry’ego [93] i równieŜ skrytykowanego [94] zjawiska super-
soczewkowania. Z wytworzonego metamateriału wykonano płasko-równoległa płytkę
o grubości 6 cm i umieszczono ją w prowadnicy falowej w odległości 2 cm
od mikrofalowej anteny nadawczej. Antena odbiorcza po przeciwnej stronie płytki miała
za zadane rejestrować przechodzące promieniowanie EM. Dla konfiguracji układu z górną
płaszczyzną uniesioną na 2 mm, pochodząca z punktowego źródła promieniowania (anteny
mikrofalowej) fala EM padająca na płytkę została skupiona po jej przeciwnej stronie.
Gdy górna płaszczyzna była obniŜona, skupienie nie występowało, co potwierdziło, Ŝe jest
ono moŜliwe tylko dla płytki z ujemnym współczynnikiem załamania.
PoniewaŜ we wszystkich badanych pryzmatach fala EM propagowała się zgodnie
z prawem załamania Snelliusa, eksperyment ten stanowił wyraźne potwierdzenie istnienia
zjawiska ujemnego załamania. Ponadto wyniki przedstawione dla płasko-równoległej
płytki o ujemnym współczynniku załamania były wstępnym potwierdzeniem słuszności
tezy Veselago [2] i Pendry’ego [93]. Kolejne wnioski z wykonanych przez grupę z MIT
doświadczeń dotyczą wytwarzania metamateriałów. Metamateriał moŜe posiadać
właściwości, których nie mają Ŝadne z tworzących go elementów. Ponadto bardzo waŜna
jest precyzja wykonania modelu. ZauwaŜono [58], Ŝe jeŜeli po cięciu ich laserem płytki
z naniesionymi elementami przewodzącymi nie zostały dokładnie wyczyszczone, obecna
na nich cienka warstwa węgla uniemoŜliwiał obserwację zjawiska ujemnego załamania.
18
Wartość zmierzona to n= 1,52
±
0,07, wartość rzeczywista n
teflonu
= 1,5.
19
Dla pryzmatu o kącie
ϕ
= 18
o
kąt załamania wynosił
θ
= – (6,4
±
2,4º), co odpowiadało wartości
n = – (0,36
±
0,13), zaś dla pryzmatu o kącie
ϕ
= 26º kąt załamania
θ
= – (9
±
2 º) , czyli współczynnik
załamania n = – (0,35
±
0,08).
49
Inne podejście eksperymentalne przedstawione zostało w publikacjach [59], [60].
Podobnie jak w doświadczeniu [58] pryzmaty wykonane były z metamateriałów
zbudowanych z równoległych przewodników i pierścieni, jednak komórka elementarna
zawierała dodatkowy metalowy drut, co pozwoliło na rozszerzenie operacyjnego pasma
częstotliwości. Doświadczenia i symulacje
numeryczne przeprowadzone zostały
dla pryzmatów o kątach 12
°
i 32
°
w zakresie 10–15 GHz, a ich wyniki były zgodne
z teorią. Ponadto uwzględniono rolę strat związanych z rozpraszaniem promieniowania,
które mają bardzo duŜe znaczenie, gdy rozwaŜa się moŜliwości praktycznego zastosowania
metamateriałów. Zidentyfikowano przyczyny strat, zaproponowano sposób wytworzenia
metamateriału,
który
charakteryzuje
się
znacznie
mniejszym
rozpraszaniem
i doświadczalnie sprawdzono jego właściwości [61]. Przyczyną strat była skończona
konduktywność drutów oraz straty pochodzące od dielektrycznego podłoŜa i uŜytego
spoiwa, które odgrywają największą rolę przy duŜej koncentracji pola EM. Symulacje
komputerowe wykazały, Ŝe największa koncentracja pola EM występuje w przerwie
miedzianych pierścieni. Usunięcie dielektryka i spoiwa z tego obszaru spowodowało
poprawę właściwości transmisyjnych metamateriału z 80% do 90% dla próbki o grubości
1cm. Wykazano równieŜ, Ŝe uŜycie metalowych elementów nie grubszych niŜ 3–5
głębokości wnikania promieniowania pozwala zminimalizować straty spowodowane
rozpraszaniem na metalu.
Rys. 31
Komórka elementarna 901 HWD o wymiarach C = 0,025 cm, D = 0,03 cm,
G = 0,046 cm, H = 0,0254 cm, L = 0,33 cm, S = 0,0263 cm, T = 17,0
×
10
-4
cm,
W = 0,025 cm, V = 0,255 cm [61].
Tym samym ujemne załamanie zostało dobitnie potwierdzone praktycznymi
pomiarami. Rozpoczęła się faza nasilonych prac nad miniaturyzacją wytworzonych
struktur i przesunięcia przedziału częstotliwości, w którym współczynnik załamania jest
ujemny do zakresu optycznego.
50
III.2 Modulacja transmisji fali EM
Struktury wielowarstwowe zawierające materiały o ujemnym współczynniku załamania
mogą być traktowane jak sekwencja płasko-równoległych soczewek [90], których
wzajemne efekty pochodzące od płytek materiału ujemnego i przestrzeni materiału
dodatniego (powietrza) znoszą się. Takie struktury periodyczne o uśrednionym
współczynniku załamania równym zero, mają tę właściwość, Ŝe transportowana przez nie
fala jest modyfikowana, co ma wpływ na wzmocnienie lub osłabienie transmisji przez taki
układ. Wprowadzając defekt w postaci warstwy o innej grubości zaburzającej
periodyczność układu udowodniono, Ŝe moŜna zmieniać sposób rozchodzenia się w nim
promieniowania [90]. Metamateriał utworzony był przez trójwymiarową sieć komórek
elementarnych zbudowanych z drutów i rozproszonych rezonatorów kołowych
analogicznych do tych opisanych w pracach [32],[36], [37] (Rys. 32).
Rys. 32
Schemat wielowarstwowej struktury złoŜonej z naprzemiennie ułoŜonych warstw
metamateriału i powietrza o grubości a, zawierający warstwę o grubości b
wprowadzającą defekt do układu; elementarna komórka metamateriału [90].
Główny wkład do efektywnej przenikalności elektrycznej mają druty, zaś do efektywnej
przenikalności magnetycznej rozproszone rezonatory kołowe [5], [92]. Efektywna
przenikalność elektryczna
( )
(
)
ε
γ
ω
ω
ω
ω
ε
i
p
−
−
=
2
1
,
(3.3)
gdzie
⋅
≈
w
p
r
d
d
c
ln
2
π
ω
to efektywna częstość plazmowa,
(3.4)
⋅
⋅
=
w
r
d
S
c
ln
2
2
σ
γ
ε
to tłumienie przenikalności elektrycznej,
(3.5)
51
σ
to przewodność drutu, zaś S to efektywny przekrój drutu określony jako
(
)
<
−
⋅
⋅
⋅
≈
≥
⋅
=
w
w
r
r
S
δ
δ
δ
π
δ
π
dla
r
2
S
r
dla
w
w
2
,
(3.6)
gdzie
πσω
δ
2
c
=
to penetracja promieniowania EM w strukturę materiału.
(3.7)
Aby obliczyć efektywną przenikalność magnetyczną sieci złoŜonej z rozproszonych
rezonatorów kołowych (SRR) ich magnetyzację naleŜy określić jako
r
r
m
I
R
c
n
M
2
2
π
=
,
(3.8)
gdzie
3
3
d
n
m
=
to liczba SRR przypadająca na jedną komórkę elementarną,
(3.9)
r
R
to promień pojedynczego rezonatora kołowego, zaś
r
I
to prąd płynący przez rezonator
kołowy.
Rezonator kołowy moŜna traktować jak efektywny obwód oscylujący o indukcji L
i rezystancji drutu R oraz pojemności C przestrzeni wewnątrz SRR. W takim obwodzie siła
elektromotoryczna (
Π
) pochodzi od zmiennego w czasie pola magnetycznego
rozchodzącej się fali elektromagnetycznej. Przy takich załoŜeniach zmiany prądu
r
I
w pojedynczym SRR opisać moŜna następująco
t
r
I
C
t
r
I
R
t
r
I
L
∂
Π
∂
=
+
∂
∂
+
∂
∂
1
2
2
,
(3.10)
gdzie H’ to lokalna (mikroskopowa) wartość pola magnetycznego, róŜna od wartości
ś
redniej (makroskopowej) pola magnetycznego H, zaś
t
H
c
r
R
∂
∂
⋅
=
Π
'
2
π
.
(3.13)
JeŜeli liczba SRR przypadających na objętość
3
d (czyli pojedynczą komórkę elementarną)
jest wystarczająco duŜa, tablicę SRR traktujemy jako układ dipoli magnetycznych
i zaleŜność między polami magnetycznymi mikroskopowym i makroskopowym opisujemy
jako
M
B
M
H
H
3
8
3
4
'
π
π
−
=
+
=
.
(3.11)
52
W efekcie z równań (3.8)
−
(3.11) otrzymujemy wyraŜenie na efektywną przenikalność
magnetyczną zaleŜną od częstości promieniowania elektromagnetycznego
( )
(
)
µ
ωγ
ω
ω
ω
ω
µ
i
F
F
+
+
+
+
=
3
1
1
2
2
0
3
,
(3.12)
gdzie
2
2
2
=
c
R
L
n
F
r
m
π
π
to współczynnik wypełnienia komórki elementarnej,
(3.13)
LC
1
2
0
=
ω
to częstość rezonansowa obwodu LC,
(3.14)
L
R
=
µ
γ
to tłumienie przenikalności magnetycznej.
(3.15)
Indukcja L, rezystancja R i pojemność C wynoszą odpowiednio
−
=
4
7
8
ln
4
2
r
R
c
R
L
r
r
π
,
(3.16)
r
r
S
R
R
σ
π
2
=
,
(3.17)
g
d
r
C
π
π
4
2
=
,
(3.18)
gdzie r to średnica drutu tworzącego rezonator kołowy,
g
d to szerokość przerwy
w pierścieniu SRR,
r
S
to efektywny przekrój drutu tworzącego SRR, definiowany
analogicznie jak (3.6)
(
)
<
−
≈
≥
=
r
S
r
r
S
r
r
δ
δ
δ
π
δ
π
dla
2r
dla
2
.
(3.19)
Pojemność C jest niezmienna, poniewaŜ
r
d
g
<<
.
Dla materiału o parametrach d = 1 cm,
w
r = 0,05 cm,
r
R
= 0,2 cm, r = 0,02 cm,
g
d = 10
3
cm,
[ ]
s
1
10
2
19
⋅
=
σ
otrzymano zaleŜność rzeczywistych części przenikalności
( )
ω
ε
oraz
( )
ω
µ
od częstotliwości, których wykresy przedstawia Rys. 33. Częstość
rezonansowa w tym wypadku wynosiła 5,82 GHz, zaś pasmo dla którego występują
jednocześnie ujemne przenikalności elektryczna i magnetyczna zawiera się pomiędzy
5,82 GHz a 5,96 GHz.
53
Rys. 33
ZaleŜność rzeczywistej składowej przenikalności elektrycznej i magnetycznej
od częstotliwości [90].
Układ tworzyło siedem warstw metamateriału umieszczonych w powietrzu. Liczba
warstw tworzących układ (Rys. 32) została tak wybrana, aby zminimalizować straty.
Periodyczność struktury zaburzał wprowadzony defekt strukturalny w postaci warstwy
ujemnego materiału o innej grubości. WaŜną cechą badanej struktury periodycznej, mającą
decydujący wpływ na jej właściwości elektromagnetyczne, był jej uśredniony
współczynnik załamania równy zero <n> = 0. Efekt taki uzyskano dzięki zastosowaniu
metamateriału o n =
−
1, a takŜe optymalnemu wyborowi liczby warstw, ich grubości
i szerokości przerw pomiędzy warstwami. Udowodniono, Ŝe przy takiej konfiguracji
w strukturze pojawiają się dodatkowe przerwy wzbronione, które nie występują
w strukturach warstwowych utworzonych wyłącznie z dodatnich materiałów. ZauwaŜono,
Ŝ
e maksimum transmisyjne znajduje się wewnątrz pasma <n> = 0, jeŜeli grubość warstw
jest całkowitą wielokrotnością połowy długości fali padającego promieniowania.
Szerokość pasm transmisyjnych zmniejszała się wraz ze wzrostem liczby lub grubości
warstw. Współczynnik transmisji wewnątrz odkrytych przerw wzbronionych mógł zanikać
nawet dla bardzo cienkich warstw materiału, w czym upatrywano moŜliwości wytworzenia
efektywnych zwierciadeł pracujących w zakresie mikrofalowym.
III.3 Istnienie fal wstecznych
A.Grbic i G.V.Eleftheriades pod koniec roku 2002 zaprezentowali wyniki
przeprowadzonego eksperymentu mającego na celu potwierdzić zjawisko wstecznego
rozchodzenia się fali EM w ośrodkach ujemnych i moŜliwość zastosowania takich
metamateriałów do budowy anten mikrofalowych [102].
Naładowane cząstki poruszające się w ośrodku z prędkością większą niŜ prędkość
fazowa światła emitują spójne promieniowanie zwane promieniowaniem Czerenkowa.
Kąt wypromieniowywanego stoŜkowego frontu falowego zaleŜy od wzajemnego stosunku
prędkości cząstki V i prędkości fazowej fali EM w otaczającym ośrodku (v
(f)
) zgodnie
ze wzorem
V
n
c
V
v
f
0
)
(
/
cos
=
=
θ
54
gdzie n
0
– współczynnik załamania otaczającego ośrodka, c – prędkość światła,
θ
– kąt
pomiędzy wektorem prędkości cząstki a wektorem falowym odbitej fali EM.
Z właściwości funkcji cosinus wynika, Ŝe kąt
θ
w ośrodku ujemnym będzie
rozwarty (
θ
> 90
°
), co oznacza wsteczną propagację frontu falowego emitowanej fali EM.
JeŜeli rozwaŜamy periodyczną strukturę przewodzącą prędkość cząsteczki V zostanie
zastąpiona przez prędkość fazową fali EM w tej strukturze. Wówczas kąt
wypromieniowywania energii wyznaczyć moŜna jako
0
1
1
0
/
/
cos
n
n
n
c
n
c
=
=
θ
gdzie n
1
– współczynnik załamania próbki.
Gdy n
1
<0, emitowane promieniowanie jest skierowane wstecznie do otaczającego strukturę
ośrodka dodatniego, co oznacza istnienie wstecznej fali, której istnienie przewidział
w 1968 roku Veselago [2].
Zaproponowany model ośrodka ujemnego (opisany za pomocą teorii linii
transmisyjnych) został zastosowany w roli anteny odbiorczej i umieszczony został w polu
anteny nadawczej. Obracając próbkę z krokiem 1
°
, zbadano rozkład transmitowanej przez
strukturę mocy w zaleŜności od kąta padania z anteny nadawczej. Uzyskany rozkład
promieniowania elektromagnetycznego wokół próbki wyraźnie potwierdził istnienie
ujemnego promieniowania Czerenkowa. Potwierdzono tym samym kolejną moŜliwość
wykorzystania zjawiska ujemnego załamania fali EM. Obecność fal wstecznych
w ujemnych metamateriałach czyni je dobrymi kandydatami do budowy płaskich,
kompaktowych urządzeń mikrofalowych takich jak mieszacze fal lub soczewki
mikrofalowe. Miniaturyzacja tych urządzeń moŜliwa dzięki zastosowaniu sztucznych
materiałów umoŜliwiłaby ich zastosowanie w komunikacji bezprzewodowej, radarach,
monitoringu i bezprzewodowej transmisji mocy.
III.4 Ujemne załamanie
ś
wiatła
Doświadczalnego potwierdzenia właściwości metamateriałów w zakresie optycznym
dostarczyła grupa badawcza z Purdue University [85] Uzyskane przez nich wyniki
dla metamateriału zbudowanego z par złotych nanodrutów, poparte trójwymiarowymi
symulacjami numerycznymi FDTD
20
. Dzięki uŜyciu róŜnych próbek pokazano, Ŝe rodzaj
podłoŜa i współczynnik wypełnienia metalem maja kluczowy wpływ na właściwości
kompozytu: wartość współczynnika załamania i jego znak. Zaprezentowana w [85] metoda
umoŜliwia bezpośredni pomiar wartości i znaku przesunięcia fazowego dla fali świetlnej
o polaryzacji liniowej „s” i „p” przy przejściu przez badaną próbkę. Czyni to ją bardziej
uniwersalną w porównaniu do metody zaprezentowanej w [86], gdzie mierzono zmiany
fazy bez wykorzystania interferometrów, co wymagało dobrej znajomości próbki przed
pomiarami.
20
FDTD – Finite Difference Time Domain
55
Dla transparentnych materiałów optycznych przesunięcie fazy fali przechodzącej
przez próbkę o grubości d wyraŜa się
λ
π
φ
nd
trans
⋅
=
2
, więc ujemne przesunięcie fazy
0
<
trans
φ
oznacza, Ŝe ośrodek posiada ujemny współczynnik załamania
0
<
n
.
W doświadczeniach interferometrycznych przesunięcie fazy w materiale moŜe być
dokładnie zmierzone dzięki wykorzystaniu kontrolnej warstwy powietrza o tej samej
grubości jako
pow
trans
φ
φ
φ
−
=
, gdzie przesunięcie fazy fali w powietrzu to
λ
π
φ
d
pow
2
=
zatem n jest ujemne w badanym materiale jeŜeli
pow
φ
φ
−
<
. Płytkę zbudowaną z komórek
zawierających metal i dielektryk charakteryzuje absorpcja i odbicia na powierzchniach
granicznych, więc związek przesunięcia fazy
trans
φ
i współczynnika załamania n jest
bardziej skomplikowany. Z tego względu pomiar przesunięcia fazy fali przechodzącej
przez układ został uzupełniony pomiarami amplitudy fali transmitowanej i odbitej.
Zespolony współczynnik załamania warstwy zawierającej pary nanodrutów został
wyznaczony na podstawie równania
(
)
(
)
1
1
1
cos
2
2
−
+
+
+
−
=
p
p
p
n
rt
t
n
t
n
r
nkd
,
gdzie n
p
to współczynnik załamania podłoŜa
21
zaś r,t to amplitudowe współczynniki
reflektanci i transmitancji zmierzone doświadczalnie. W symulacji FDTD do opisania
przenikalności elektrycznej złota wykorzystano model Debye’a oparty na modelu Drudego
(
)
(
)
ω
ω
ε
14
2
2
16
10
0027
,
1
10
3673
,
1
0
,
9
⋅
+
⋅
−
=
i
Au
, gdzie
λ
π
ω
c
2
=
.
Amplitudowe współczynniki transmisji
2
t
T
=
i reflektacji
2
r
R
=
mierzone były
na spektrofotometrze przy uŜyciu światła spolaryzowanego liniowo. Spektrum
transmisyjne zmierzono dla normalnego padania światła, spektrum refleksyjne dla światła
padającego pod niewielkim kątem do normalnej – 8
°
. Ujemny współczynnik załamania
zmierzony został poprzez bezpośrednie pomiary fazy i amplitudy dla promieniowania EM
o częstości z pobliŜa zakresu telekomunikacyjnego (długość fali 1,5
µ
m) przy
zastosowaniu technik interferencyjnych (Rys. 34).
W interferometrze polaryzacyjnym dwa kanały optyczne mają tę samą drogę
geometryczną i róŜnią się polaryzacją światła. RóŜnicę faz pomiędzy ortogonalnie
spolaryzowanymi falami wyznacza się jako
⊥
−
=
∆
φ
φ
φ
||
. Drugi interferometr
22
posiada
dwa kanały optyczne o róŜnych ścieŜkach geometrycznych o tej samej polaryzacji. Wynik
odczytywany jest jako
pow
próbki
φ
φ
δφ
−
=
i opisuje stosunek zmiany fazy wywoływanej
przez próbkę
próbki
φ
w odniesieniu do zmiany fazy wywołanej przez kontrolną warstwę
powietrza
pow
φ
o jednakowej grubości.
21
W doświadczeniu [85] wykorzystano szkło o n
p
= 1,48
22
Określany w literaturze jako „walk-off interferometer” [85]
56
Rys. 34
Schemat pomiaru transmisji i odbicia fali EM przy uŜyciu interferometru
(a) polaryzacyjnego (b) róŜnicowego (na podstawie [85]).
Efekt rozszczepienia („walk-off”) przez kalcyt jest uŜywany do rozdzielenia dwóch
promieni a następnie do ich złączenia ich, Ŝeby wywołać interferencję jak zostało
to pokazana na Rys. 34 (b). Przesunięcia fazy dla obu polaryzacji światła zmierzone
drugim interferometrem
||
δφ
i
⊥
δφ
zostały porównane z otrzymaną z pierwszego typu
interferometru
⊥
−
=
∆
δφ
δφ
φ
||
. Błąd metody dla interferometru polaryzacyjnego określony
został na
±
1,7
°
. ZauwaŜono, Ŝe grubość podłoŜa nie ma wpływu na pomiar róŜnicy faz
⊥
−
=
∆
δφ
δφ
φ
||
, co jest typowe dla interferometrów ze wspólną ścieŜką promienia.
Dla drugiego interferometru zmiana grubości podłoŜa jest źródłem dodatkowego błędu
metody i moŜe zwiększyć go do
±
4
°
.
Metamateriał z równoległych par złotych nanodrutów wykazał ujemny
współczynnik załamania n
≅
−
0,3 dla promieniowania o długości fali 1,3
µ
m
spolaryzowanego poprzecznie magnetycznie (polaryzacja „p”). W tym przypadku światło
wzbudzało jednocześnie rezonans elektryczny i magnetyczny. Dla innych konfiguracji
współczynnik załamania był bliski zeru, ale cały czas dodatni (najniŜszy wynosił n = 0,08
dla 1,1
µ
m dla wielowarstwowej struktury analogicznej do tej opisanej w pracy [87]).
Dla światła o polaryzacji „s” wzbudzony został rezonans elektryczny (ujemna
przenikalność elektryczna) dla promieniowania o długości fali 800 nm, jednak nie udało
się dla tak krótkich fal otrzymać ujemnej przenikalności magnetycznej.
57
Ujemne załamanie światła potwierdzone zostało dopiero niedawno, w sierpniu
2006 roku [88]. Wcześniej dla podobnego zakresu częstotliwości udało się wytworzyć
i zbadać metamateriał o ujemnej przenikalności magnetycznej [89] zbudowany ze słupków
złota o wysokości 85
±
5 nm i średnicy 100 nm, dla których rezonans magnetyczny
wywoływany był przez falę EM o długości 670 nm. W doświadczeniu G.Dolling’a i in.
[88] do wytworzenia badanego metamateriału (zob. rozdział II.7) uŜyto srebrnych
nanodrutów, co pozwoliło na znaczne zmniejszenie strat i przesunięcie uŜytkowego
zakresu częstotliwości w zakres widzialny. W symulacjach FDTD wykorzystano model
Drudego, poniewaŜ dla częstotliwości z zakresu optycznego pozwala on wystarczająco
dokładnie określić właściwości dielektryczne srebra. Aby bez wątpienia określić parametry
metamateriału pomiary interferometryczne róŜnicy faz uzupełniono badaniem czułości
fazy.
Rys. 35 ZaleŜność transmitancji i reflektanci badanego metamateriału uzyskana
(a) doświadczalnie i (b) symulacyjnie [88].
Współczynnik załamania dla rozwaŜanej polaryzacji fali EM o długości 780 nm w badanej
próbce przedstawionej na Rys. 27 wyniósł n =
−
0,6.
58
IV Zastosowania
Ujemne załamanie fali EM, poza swoim czysto teoretycznym charakterem, wzbudza
ogromne zainteresowanie i emocje właściwie od roku 1967. JuŜ Veselago [2] zwrócił
uwagę na pomijane dotąd rozwiązania równań Maxwella i opisał zachowanie się fali
elektromagnetycznej padającej na płytkę o ujemnym współczynniku załamania. Ponad
30 lat później dokładnego i obszernego opisu realizacji takiego materiału dokonała w swej
pracy grupa badawcza z Uniwersytetu Kalifornijskiego w San Diego [32], zaś Pendry
rozwinął koncepcję Veselago opisując sposób działania idealnej soczewki zbudowanej
z płasko-równoległej płytki materiału o ujemnym współczynniku załamania [93]. Artykuł
ten podzielił świat naukowców na dwie grupy: tych, którzy optymistycznie podchodzili
do potencjalnych moŜliwości, jakie daje zastosowanie metamateriałów we współczesnej
optyce oraz sceptyków tej idei [57], [94]
−
[96].
IV.1 Perfekcyjna soczewka
Jednym z najczęściej stosowanych i najdłuŜej znanych elementów optycznych jest
soczewka. MoŜe być uŜywana do skupiania bądź kierowania promieniowania, co pozwala
na wykorzystanie jej w ogromnej ilości układów i dla wielu długości fali – od fal
radiowych po optyczne. JeŜeli uświadomimy sobie jak szeroki jest zakres jej zastosowań,
zrozumiemy dlaczego wokół moŜliwości wykorzystania zjawiska ujemnego załamania fali
elektromagnetycznej w celu wytworzenia soczewki idealnej rozpętała się tak Ŝarliwa
dyskusja pośród naukowców.
Zjawisko załamania stanowi podstawę istnienia soczewek i odwzorowywania
obrazów. KaŜdy materiał o współczynniku załamania innym niŜ współczynnik załamania
otoczenia zmienia kierunek biegu fali elektromagnetycznej przy jej przejściu przez granicę
tych dwóch ośrodków. Kiedy fala pada na powierzchnię styku dwóch materiałów
pod przypadkowym kątem, kierunek rozchodzenia się fali po wejściu do drugiego ośrodka
zmienia się o wartość zaleŜną od współczynników załamania obu ośrodków. Ilościowej
zaleŜności pomiędzy kątem padania
θ
1
i załamania
θ
2
mierzonymi od normalnej
do powierzchni padania, oraz współczynnikami załamania obu materiałów n
1
i n
2
dostarcza
prawo Snelliusa (otrzymane przy Ŝądaniu równości fazy fali padającej i przechodzącej)
2
2
1
1
sin
sin
θ
θ
n
n
=
Idealna soczewka płaska
Veselago przewidział [2], Ŝe materiał o jednocześnie ujemnych względnych
przenikalnościach równych
−
=
−
=
1
1
r
r
µ
ε
,
(3.20)
będzie wykazywał ujemny współczynnik załamania
1
−
=
n
, a materiał taki będzie
zachowywał się jak soczewka skupiająca. W roku 2001 istnienie zjawiska ujemnego
59
załamania fali potwierdzono eksperymentalnie [54] a w kolejnych latach zaprojektowane
zostały materiały, które umoŜliwiają otrzymanie ujemnych przenikalności elektrycznej
[5]
−
[7], [32] i magnetycznej [6], [7], [32], [97], [98], takŜe przy wykorzystaniu struktur
fotonicznych [99]. Zaproponowaną przez Veselago teorię rozwinął Pendry [93], który
wykazał, Ŝe jeŜeli warunek (3.20) zostanie spełniony bardzo dokładnie, rozwaŜana
soczewka będzie skupiała padające promieniowanie idealnie punkt w punkt
23
. Sposób
na realizację takiej soczewki Pendry dostrzegł w metamateriałach wykonanych po raz
pierwszy przez grupę Smith’a [32].
Powszechnie wiadomo, iŜ zdolności skupiające soczewki, a więc rozdzielczość
obrazu, ograniczona jest przez długość promieniowania padającego na układ. śadna
powszechnie stosowana soczewka nie jest w stanie skupić promienia fali EM na obszarze
mniejszym niŜ kwadrat długości fali. WiąŜe się to z faktem, iŜ cała informacja
o przedmiocie przenoszona przez falę EM podzielona jest na dwa obszary – bliski i daleki.
Tradycyjne soczewki są w stanie zrekonstruować tylko dalekie składowe pola EM,
zaś informacja zawarta w bliskim polu w postaci zanikających fal przypowierzchniowych
nie moŜe zostać odtworzona i jest bezpowrotnie tracona, co ma decydujący wpływ
na jakość uzyskanego obrazu.
a)
b)
Rys. 36 (a)
Bieg promieni w płytce z materiału o ujemnym współczynniku załamania;
(b)
Płytka z materiału o współczynniku załamania n=-1 załamuje światło pod ujemnym
kątem względem normalnej. Światło skupiane jest dwukrotnie: wewnątrz płytki i poza nią.
Według Pendry’ego alternatywą dla tradycyjnej soczewki jest płasko–równoległa
płytka z materiału o ujemnym współczynniku załamania (Rys. 36). Światło przechodzące
przez płytkę o grubości d
2
umieszczoną w odległości d
1
od źródła promieniowania
jest skupiane w odległości z = d
2
– d
1
po drugiej stronie płytki. Przenikalność elektryczna
i magnetyczna materiału równe są –1 i współczynnik załamania jest ujemny zgodnie
z równaniem (1.15). Impedancja ośrodka pozostaje w tym przypadku dodatnia
23
W literaturze zjawisko to określane jest mianem „superskupienia” lub „supersoczewkowania”.
60
0
0
εε
µµ
=
Z
,
(4.7)
i jeŜeli taka soczewka umieszczona jest w próŜni, nie występuje odbicie i promieniowanie
jest w całości transmitowane przez układ.
Pendry rozpatrzył punktowe źródło promieniowania elektromagnetycznego
o częstości
ω
, którego pole
)
exp(
)
,
(
)
,
(
t
i
y
ik
x
ik
z
ik
k
k
E
t
r
E
y
x
z
y
x
ω
−
+
+
×
=
∑
,
(4.1)
propaguje się wzdłuŜ osi z soczewki. Zgodnie z równaniami Maxwella mamy
2
2
2
2
y
x
z
k
k
c
k
−
−
+
=
ω
,
2
2
2
2
y
x
k
k
c
+
>
ω
.
(4.2)
Zadaniem soczewki jest korekta fazy kaŜdej ze składowych Fouriera pola
elektromagnetycznego tak, aby w pewnej odległości od soczewki nastąpiło skupienie
i powstał obraz źródła promieniowania.
JeŜeli
2
2
2
2
c
k
k
i
k
y
x
z
ω
−
+
+
=
,
2
2
2
2
y
x
k
k
c
+
<
ω
,
(4.3)
fale przypowierzchniowe zanikają ekspotencjalnie wzdłuŜ osi z i ich amplituda
nie jest odtwarzana. Powstały obraz źródła promieniowania tworzony jest tylko z fal
rozchodzących się wzdłuŜ osi z. PoniewaŜ fale te ograniczone są do
2
2
2
2
c
k
k
y
x
ω
<
+
,
(4.4)
więc największa moŜliwa do uzyskania rozdzielczość to
λ
ω
π
π
=
=
≈
∆
c
k
2
2
max
(4.5)
W przypadku tych fal ich zanikanie ma związek z amplitudą, nie zaś z fazą, zatem aby
temu zapobiec konieczne jest ich wzmocnienie, nie zaś korekta fazy. Pendry udowadnia
[93], Ŝe fale przypowierzchniowe doznają potrzebnego wzmocnienia w procesie transmisji
przez płytkę o ujemnym współczynniku załamania, zatem metamateriały mają zdolność
zapobiegania zanikaniu fal przypowierzchniowych.
61
Idealna soczewka sferyczna
Płytka z materiału o ujemnym współczynniku załamania skupia promieniowanie
elektromagnetyczne, jeŜeli zaś współczynnik ten równy jest –1 zachowuje się jak idealna
soczewka skupiająca produkując idealny obraz odwzorowywanego przedmiotu. Jednak
powstały w ten sposób obraz jest w skali 1:1 i powiększenie go nie jest moŜliwe.
Uzyskanie idealnego i powiększonego obrazu przedmiotu wymaga zastosowania soczewki
skupiającej o zakrzywionej powierzchni, co pociąga za sobą konieczność zmiany definicji
przenikalności elektrycznej
ε
i magnetycznej
µ
materiału tak, aby były one funkcjami
połoŜenia.
W roku 2003 Pendry przedstawił sposób na uzyskanie powiększenia
w perfekcyjnych soczewkach [100]. Aby uzyskać idealną soczewkę powiększającą,
powierzchnia materiału musiała być zagięta, przy jednoczesnym zachowaniu wszystkich
jego właściwości fizycznych. Pendry zaproponował dwie moŜliwości budowy takiej
soczewki (Rys. 37).
Rys. 37
Soczewka o cylindrach współosiowych (a) i stycznych wewnętrznie (b);
2
2
1
2
a
b
r
r
=
(na podstawie [100])
Deformacja kształtu powierzchni pociągała za sobą zmianę charakteru parametrów
takiego materiału. Zastosowana przez Pendry’ego transformacja między tymi geometrią
planarną i cylindryczną szczegółowo opisana została w [100]. Efektem takiej transformacji
była cylindryczna soczewka, która odtwarzała i powiększała bez zniekształceń zawartość
mniejszego cylindra (o promieniu a) na zewnątrz większego cylindra (o promieniu b)
lub wytwarzała pomniejszony idealny obraz obiektu umieszczonego na zwenątrz
większego cylindra wewnątrz cylindra mniejszego. Współczynnik powiększenia soczewki
to
2
2
a
b
. Ponadto jest to soczewka jest krótkoogniskowa i tylko obiekty oddalone od jej
osi o mniej niŜ
a
b
r
2
=
mogą być odwzorowane wewnątrz przestrzeni pomiędzy
cylindrami, zaś obiekty umieszczone w tej przestrzeni jak równieŜ te umieszczone bliŜej
ś
rodka niŜ
b
a
r
2
=
nie wytworzą obrazu na zewnątrz duŜego cylindra. Przenikalności
elektryczna
ε
i magnetyczna
µ
opisanej soczewki muszą być proporcjonalne
do odwrotności kwadratu promienia
2
1
r
[100].
62
IV.2 Urz
ą
dzenia mikrofalowe
Metamateriały o ujemnym współczynniku załamania znalazły zastosowanie takŜe
w budowie róŜnego typu urządzeń mikrofalowych. W tym zakresie jednym z najbardziej
aktywnych metamateriałoznawców jest George V. Eleftheriades. Opisany przez niego
model linii transmisyjnych dla metamateriału o ujemnym współczynniku załamania [53],
[65] okazał się być bardzo uŜyteczny, ze względu na swoją prostotę i łatwość, z jaką
moŜna przy jego uŜyciu budować dwuwymiarowe obwody mikrofalowe.
Kompaktowy, jednowymiarowy modulator fazy zbudowany z naprzemiennie
ułoŜonych warstw materiału ujemnego i dodatniego [101] posiada kilka niewątpliwych
zalet w porównaniu do swojego konwencjonalnego odpowiednika. Zastosowanie warstw
o ujemnym
współczynniku
załamania
pozwoliło
zminiaturyzować
modulator
a jednocześnie uniezaleŜnić jego działanie od wymiarów. Ponadto udowodniono,
Ŝ
e odpowiedni
dobór
wartości
poszczególnych
elementów
pojemnościowych
i indukcyjnych w modelu linii transmisyjnych uŜytego metamateriału pozwala
na osiągnięcie dodatniego, ujemnego lub zerowego przesunięcia fazy fali EM,
przy jednoczesnym zachowaniu niewielkich wymiarów urządzenia. Ze względu na swoje
niewielkie wymiary, planarną geometrię i liniową odpowiedź, modulatory wykorzystujące
ujemne załamanie są łatwe do zintegrowania ich z innymi urządzeniami mikrofalowymi.
MoŜliwość wykorzystania metamateriałów o ujemnym współczynniku załamania
jako anteny mikrofalowej [102], [103] opisana została pod koniec 2002 roku przez
A.Grbic’a i G.V. Eleftheriades’a. Zarówno symulacje jak i eksperyment potwierdziły
dyspersyjny
charakter
odbiorczej
anteny
mikrofalowej
(zob.
rozdział
III.3)
wykorzystującej ujemne materiały: kąt, pod jakim rozchodziła się emitowana fala
wsteczna jest silnie zaleŜny od częstości uŜytego padającego promieniowania
ω
zgodnie
z równaniem
d
LC
B
ω
1
−
≈
, gdzie
d
LC
v
g
2
)
(
ω
≅
.
63
IV.3 Najnowsze odkrycia
Najnowszy sposób wykorzystania metamateriałów przedstawiony został przez Pendry’ego,
Schuringa i Smitha [104] w maju 2006 roku. Zaprezentowali oni sposób na dowolne
kształtowanie linii pola EM przy uŜyciu wydrąŜonej kuli z metamateriału o ujemnym
współczynniku załamania, do której wnętrza promieniowanie EM o określonej częstości
nie wnika (Rys. 38).
Rys. 38
Bieg promieni wewnątrz układu, widok w dwóch (A) i w trzech wymiarach (B)
[104].
Odpowiednio zaprojektowany metamateriał umoŜliwia praktycznie dowolne
zniekształcenie
trajektorii
padającego
promieniowania:
moŜna
spowodować,
by promieniowanie ominęło pewną część przestrzeni, a następnie wróciło na swój
pierwotny tor (Rys. 39).
Rys. 39
Linie pola EM pochodzące od punktowego źródła promieniowania zakrzywiają się
wokół obszaru wewnątrz kuli, a następnie wracają na swój pierwotny tor [104].
Dystorsję pola EM opisano stosując transformacje przestrzenną współrzędnych
do określenia wartości przenikalności elektrycznej
ε
i magnetycznej
µ
przy jednoczesnym
spełnieniu równań Maxwella. Aby skompresować pole, które pierwotnie wypełniało
64
przestrzeń
2
R
r
<
, do przestrzeni
2
1
R
r
R
<
<
pomiędzy współśrodkowymi kulami,
zastosowano współrzędne przestrzenne
=
=
−
+
=
ϕ
ϕ
θ
θ
'
'
2
1
2
1
/
)
(
'
R
R
R
r
R
r
,
(4.6)
które następnie wykorzystano do zdefiniowania przenikalności elektrycznej
ε
’ i magnety-
cznej
µ
’ w tym układzie
'
)
(
)
'
(
'
'
1
2
2
1
2
'
'
r
R
R
R
r
R
r
r
−
−
=
=
µ
ε
,
1
2
2
'
'
'
'
R
R
R
−
=
=
θ
θ
µ
ε
,
1
2
2
'
'
'
'
R
R
R
−
=
=
ϕ
ϕ
µ
ε
.
(4.7)
Dla kuli wewnętrznej
1
R
r
<
wartości
ε
’ i
µ
’ mogą być dowolne i nie mają wpływu
na rozpraszanie
promieniowania
EM,
poza
kulą
zewnętrzną
(
2
R
r
>
)
1
'
'
'
'
'
'
'
'
'
'
'
'
=
=
=
=
=
=
ϕ
ϕ
θ
θ
µ
ε
µ
ε
µ
ε
r
r
, zaś na granicy
2
R
r
=
występuje idealne
dopasowanie opisane warunkami
'
'
'
'
1
'
'
r
ε
ε
ε
ϕ
θ
=
=
oraz
'
'
'
'
1
'
'
r
µ
µ
µ
ϕ
θ
=
=
(4.8)
Dzięki temu osłona wykonana z takiego metamateriału będzie omijana przez
promieniowanie z zakresu częstotliwości, dla którego został on zaprojektowany.
Co więcej, kaŜdy otoczony nią obiekt, będzie niezauwaŜalny z zewnątrz. PowyŜszy opis
wyklucza jednak istnienie pola EM wewnątrz obszaru
1
R
r
<
, co oznacza, Ŝe obserwator
znajdujący się wewnątrz niego obserwator równieŜ nie będzie widział tego, co znajduje się
po drugiej stronie osłony metamateriałowej. Takie osłony mogłyby być uŜywane
do całkowitej ochrony urządzeń wraŜliwych na działanie promieniowania o konkretnej
częstości, a zastosowane dla zakresu widzialnego, stanowiłyby urządzenie kamuflujące.
65
V Transmitancja warstwowych układów optycznych
Metamateriały o ujemnym współczynniku załamania znajdą zapewne zastosowanie
w konstrukcji wielu urządzeń optycznych bazujących na supersieciach. Ze względu na to,
w tym rozdziale prezentujemy wzory Fresnela dla prostych układów wielowarstwowych
zawierających ośrodki ujemne i dodatnie. Supersieć jest nową klasą materiałów
niemoŜliwą do wytworzenia za pomocą powszechnie stosowanych metod wzrostu
kryształów. Składa się z ułoŜonych naprzemiennie co najmniej dwóch rodzajów warstw
materiałów dielektrycznych o grubości kilku atomów. W porównaniu do konwen-
cjonalnych kryształów, supersieci zapewniają większą elastyczność projektowania nowych
materiałów, charakteryzują się lepszą luminescencją kwantową oraz oferują moŜliwość
projektowania struktury pasmowej materiału. Wszystkie te cechy zawdzięczają ściśle
określonemu rozkładowi współczynnika załamania n, przenikalności elektrycznej
ε
i magnetycznej
µ
uzyskiwanym poprzez ustalony porządek ułoŜenia poszczególnych
warstw w strukturze supersieci opisany za pomocą wzorów rekurencyjnych.
Do analizy i obliczeń najprostsze są sieci dwuskładnikowe. Transmitancja
wielowarstwowych układów optycznych omówiona została obszernie w [15], [106].
Tutaj natomiast przedstawiona zostanie supersieć binarna, w której jeden z uŜytych
materiałów jest metamateriałem o ujemnym współczynnikiem załamania n. RozwaŜane są
dwa układy optyczne złoŜone z trzech warstw materiałów przedstawione na Rys. 40 i Rys.
41. Dla uproszczenia zakładamy, Ŝe współczynnik załamania ośrodka wejściowego
i wyjściowego są jednakowe, czyli płytka środkowa otoczona jest dwoma płytkami
wykonanymi z tego samego metamateriału o ujemnym bądź dodatnim współczynniku
załamania. Drugim poczynionym uproszczeniem jest załoŜenie, iŜ uŜywamy
bezdyspersyjnych materiałów dodatnich. Rozpatrzone zostaną dwie konfiguracje.
Pierwsza, przedstawiona jest na (Rys. 40), gdzie pomiędzy półpłaszczyznami ośrodka
ujemnego znajduje się ośrodek dodatni i druga (Rys. 41), gdzie ośrodek dodatni otacza
ośrodek ujemny.
Rys. 40
Ośrodek dodatni otoczony półpłaszczyznami z jednakowego, ujemnego
metamateriału.
66
Rys. 41
Płytka płasko-równoległa wykonana z metamateriału o ujemnym współczynniku
załamania otoczona półpłaszczyznami jednakowego, bezdyspersyjnego ośrodka
dodatniego.
Aby wyeliminować z rozwaŜań przypadki całkowitego wewnętrznego odbicia,
gdy
2
2
π
θ
=
, wybieramy tylko takie kąty padania, które wynikają z zastosowania prawa
załamania
⋅
<
⋅
2
sin
sin
2
1
1
π
θ
n
n
1
sin
2
1
1
<
⋅
n
n
θ
co oznacza, Ŝe kąt padania na pierwszą powierzchnię graniczną musi spełniać warunek
1
2
1
sin
n
n
<
θ
.
PoniewaŜ na wejściu i wyjściu fali elektromagnetycznej uŜyty został ten sam
ośrodek i współczynnik załamania wejściowy i wyjściowy są jednakowe, więc kąt wyjścia
3
θ
fali elektromagnetycznej będzie równy co do wartości kątowi wejścia
1
θ
,
co jest zgodne z prawem Snelliusa. Dla rozpatrywanych układów prawo Snelliusa
moŜna zapisać następująco
Dla układu I:
)
sin(
sin
2
2
1
1
θ
θ
−
⋅
=
⋅
−
n
n
dla pierwszej granicy ośrodków
)
sin(
sin
3
1
2
2
θ
θ
−
⋅
=
⋅
n
n
dla drugiej granicy ośrodków
Dla układu II:
)
sin(
sin
2
1
1
2
θ
θ
−
⋅
−
=
⋅
n
n
dla pierwszej granicy ośrodków
)
sin(
sin
3
2
2
1
θ
θ
−
⋅
=
⋅
−
n
n
dla drugiej granicy ośrodków
Po wykonaniu prostych przekształceń matematycznych wyraźnie widać, iŜ w kaŜdym
z przypadków kąty są sobie równe co do wartości, czyli
3
1
θ
θ
=
.
67
V.1 Ujemne załamanie fali EM
Zachowanie się fali elektromagnetycznej przy przejściu przez granicę ośrodków o róŜnych
co do wartości i przeciwnych co do znaku współczynnikach załamania ilustruje rysunek
Rys. 42
Fala elektromagnetyczna na granicy ośrodków o przeciwnych znakach
współczynnika załamania.
Zakładamy teŜ, Ŝe spełnione są warunki ciągłości linii pola elektromagnetycznego
na granicy dwóch ośrodków (1.3) oraz Ŝe częstość drgań fali elektromagnetycznej
jest stała. Wektory falowe leŜą wtedy w jednej płaszczyźnie zwanej płaszczyzną padania
i mają postać
]
sin
,
0
,
cos
[
]
,
,
0
,
,
[
j
c
j
n
j
c
j
n
j
z
k
j
x
k
j
k
θ
ω
θ
ω
=
=
(5.1)
gdzie: j = 1,2.
68
Wykorzystywana tu jest wielkość zwana liczbą falową i definiowana jako
c
k
ω
=
dla propagacji fali elektromagnetycznej w próŜni. WyraŜa ona całkowitą liczbę długości
fali zawartą w odległości 1 m Zgodnie z prawem odbicia, kąt padania fali EM na granicę
dwóch ośrodków jest równy kątowi odbicia
'
1
1
θ
θ
=
oraz spełnione jest prawo Snelliusa.
V.2 Polaryzacja typu „s” i „p”
W zaleŜności od typu polaryzacji fali elektromagnetycznej rozkład wektorów pól
oraz wektorów falowych fali EM przy przejściu przez granicę ośrodków o przeciwnych
co do znaku współczynnikach załamania przedstawiony jest na Rys. 43 oraz Rys. 44 [15].
Rys. 43
Polaryzacja typu „s” – wektor pola E = [0 , E
y
, 0].
Rys. 44
Polaryzacja typu „p” – wektor pola H = [0 , H
y
, 0].
gdzie:
)
(
1
+
E
to amplituda fali padającej,
)
(
1
−
E
to amplituda fali odbitej,
)
(
2
−
E
to amplituda
fali załamanej,
)
(
2
+
E
to amplituda fali odbitej od drugiej powierzchni granicznej.
69
V.3 Dyspersja współczynnika załamania
Ujemny współczynnik załamania metamateriałów ma charakter silnie dyspersyjny,
co oznacza, Ŝe jego wartość jest bardzo czuła na zmiany częstości promieniowania
elektromagnetycznego padającego na układ. Zjawisko to opisane jest przez następujące
związki dyspersyjne[22]
( )
ε
ε
ε
γ
ω
ω
ω
ω
ω
ε
⋅
⋅
+
−
+
=
i
p
r
2
2
2
1
,
( )
µ
µ
µ
γ
ω
ω
ω
ω
ω
µ
⋅
⋅
+
−
+
=
i
p
r
2
2
2
1
(5.2)
gdzie:
µ
ε
γ
γ
,
– tłumienie dla przenikalności elektrycznej i magnetycznej,
1
−
=
i
.
Dla ułatwienia analizy zakładamy, Ŝe rozpatrujemy tylko częstości bliskie częstości
rezonansowej
ε
ω
i
µ
ω
co oznacza brak tłumienia, czyli
0
=
=
µ
ε
γ
γ
V.4 Amplitudowe współczynniki transmisji
Na podstawie warunków granicznych (1.3) dla fali elektromagnetycznej, stwierdzić
moŜna, iŜ ciągłość składowych stycznych pola elektrycznego i magnetycznego zachowana
jest jedynie w punkcie x = 0. Pozwala to zapisać relacje między natęŜeniami fali padającej
i odbitej oraz padającej i załamanej uŜywając wzorów Fresnela. W ten sposób moŜna
określić amplitudowy współczynnik odbicia r
j,j+1
oraz amplitudowy współczynnik
załamania t
j,j+1
dla obu typów polaryzacji, gdzie indeksy j,j+1 oznaczają numery
sąsiadujących warstw układu optycznego. W rozwaŜanym tutaj przypadku układ składa się
z trzech warstw, co pozwala zapisać wyraŜenia na współczynniki transmisji t i odbicia r
dla obu typów polaryzacji w postaci analitycznej.
JeŜeli rozwaŜamy przejście fali EM z ośrodka ujemnego n
1
<0 do dodatniego n
2
>0
(pierwsza powierzchnia graniczna na Rys. 40 oraz druga powierzchnia graniczna na Rys.
41), amplitudowe współczynnik odbicia r
12
i amplitudowy współczynnik transmisji t
12
dla polaryzacji „s” (poprzeczna polaryzacja elektryczna, Rys. 43) mają postać
Amplitudowy współczynnik odbicia:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
2
2
2
1
2
2
2
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
,
1
1
,
2
2
,
1
1
,
)
(
1
)
(
1
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
µ
n
n
n
n
n
n
c
n
c
n
c
n
c
n
k
k
k
k
E
E
r
r
r
r
r
r
r
x
x
x
x
s
s
+
−
⋅
−
⋅
−
−
−
⋅
−
⋅
−
=
+
−
=
=
+
−
=
+
−
=
=
+
−
(5. 3)
70
Amplitudowy współczynnik transmisji:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
2
2
2
1
1
2
2
2
1
1
1
1
1
1
2
2
2
1
1
1
1
1
1
2
2
,
1
1
,
1
1
,
)
(
1
)
(
2
cos
cos
cos
2
cos
cos
cos
2
cos
cos
cos
2
2
µ
θ
ω
µ
θ
ω
µ
ω
ε
ω
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
n
n
n
n
c
n
c
n
c
n
k
k
k
E
E
t
r
r
r
r
r
r
x
x
x
s
s
+
−
⋅
−
⋅
−
−
⋅
−
⋅
−
⋅
=
+
⋅
=
=
+
⋅
=
+
⋅
=
=
+
−
(5. 4)
Korzystamy w tym miejscu (jak równieŜ w dalszych wyprowadzeniach)
z zaleŜności dyspersyjnych dla przenikalności elektrycznej i magnetycznej określonych
są wzorami (5.2) oraz wzoru (1.16) na ujemny współczynnik załamania.
Transmisję fali elektromagnetycznej spolaryzowanej poprzecznie magnetycznie
(polaryzacja „p”, Rys. 44) opisują natomiast
Amplitudowy współczynnik odbicia:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
2
1
2
2
2
1
2
2
2
1
2
1
2
1
2
1
2
1
2
1
2
1
1
2
2
1
2
2
2
1
2
1
1
2
2
1
2
2
2
1
2
1
,
2
2
1
2
,
2
1
2
1
,
2
2
1
2
,
2
1
1
2
2
1
)
(
1
)
(
1
)
(
1
)
(
1
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
µ
µ
µ
n
n
n
n
n
n
c
n
n
c
n
n
c
n
n
c
n
n
k
n
k
n
k
n
k
n
k
k
E
E
H
H
r
r
r
r
r
r
r
x
x
x
x
p
p
+
−
⋅
−
⋅
−
−
−
⋅
−
⋅
−
=
=
+
−
=
⋅
+
⋅
⋅
−
⋅
=
=
⋅
+
⋅
⋅
−
⋅
=
⋅
⋅
=
=
+
−
+
−
(5.5)
71
Amplitudowy współczynnik transmisji:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
2
1
2
2
2
1
2
1
2
1
2
1
2
1
1
2
1
1
2
2
1
2
2
2
1
2
1
1
2
1
2
1
,
2
2
1
2
,
2
1
2
1
,
2
1
1
2
2
1
)
(
1
)
(
2
)
(
1
)
(
2
cos
cos
cos
2
cos
cos
cos
2
cos
cos
cos
2
2
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
µ
µ
n
n
n
n
c
n
n
c
n
n
c
n
n
n
k
n
k
n
k
n
n
k
k
E
E
H
H
t
r
r
r
r
r
x
x
x
p
p
+
−
⋅
−
⋅
−
⋅
−
⋅
−
⋅
=
=
+
⋅
=
⋅
+
⋅
⋅
⋅
⋅
=
=
⋅
+
⋅
⋅
⋅
⋅
=
⋅
⋅
=
=
+
−
+
−
(5.6)
Dla przejścia fali EM z ośrodka dodatniego n
1
>0 do ujemnego n
2
<0 (druga
powierzchnia graniczna na Rys. 40 oraz pierwsza powierzchnia graniczna na Rys. 41),
poprzez analogię amplitudowe współczynnik odbicia r
12
i amplitudowy współczynnik
transmisji t
12
dla polaryzacji „s” określić moŜna jako
Amplitudowy współczynnik odbicia:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
ω
µ
θ
ω
µ
ω
ε
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
µ
r
r
r
r
r
r
x
x
x
x
s
s
n
n
n
n
n
n
c
n
c
n
c
n
c
n
k
k
k
k
E
E
r
−
⋅
−
⋅
−
+
−
⋅
−
⋅
−
−
=
+
−
=
=
+
−
=
+
−
=
=
+
−
2
1
1
1
2
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
2
1
1
1
2
2
,
1
1
,
2
2
,
1
1
,
)
(
1
)
(
1
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
(5.7)
72
Amplitudowy współczynnik transmisji:
( )
(
)
( )
(
)
( )
(
)
ω
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
r
r
r
x
x
x
s
s
n
n
n
n
n
c
n
c
n
c
n
k
k
k
E
E
t
−
⋅
−
⋅
−
+
⋅
=
+
⋅
=
=
+
⋅
=
+
⋅
=
=
+
−
1
1
1
1
1
1
1
2
2
2
1
1
1
1
1
1
2
2
2
1
1
1
1
1
1
2
2
,
1
1
,
1
1
,
)
(
1
)
(
2
cos
cos
cos
2
cos
cos
cos
2
cos
cos
cos
2
2
(5.8)
Dla fali EM o polaryzacji „p” współczynniki wynoszą
Amplitudowy współczynnik odbicia:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
( )
(
)
ω
µ
θ
ω
µ
ω
ε
µ
θ
ω
µ
θ
ω
µ
ω
ε
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
µ
µ
µ
r
r
r
r
r
r
x
x
x
x
p
p
n
n
n
n
n
n
c
n
n
c
n
n
c
n
n
c
n
n
k
n
k
n
k
n
k
n
k
k
E
E
H
H
r
−
⋅
−
⋅
−
+
−
⋅
−
⋅
−
−
=
=
+
−
=
⋅
+
⋅
⋅
−
⋅
=
=
⋅
+
⋅
⋅
−
⋅
=
⋅
⋅
=
=
+
−
+
−
1
1
2
1
1
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
2
1
1
2
2
1
2
2
2
1
2
1
1
2
2
1
2
2
2
1
2
1
,
2
2
1
2
,
2
1
2
1
,
2
2
1
2
,
2
1
1
2
2
1
)
(
1
)
(
1
)
(
1
)
(
1
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
cos
(5.9)
73
Amplitudowy współczynnik transmisji:
( )
(
)
( )
(
)
( )
(
)
( )
(
)
ω
µ
θ
ω
µ
ω
ε
µ
θ
ω
µ
θ
µ
θ
µ
θ
µ
θ
µ
θ
ω
µ
θ
ω
µ
θ
ω
µ
µ
µ
µ
µ
r
r
r
r
x
x
x
p
p
n
n
n
n
n
c
n
n
c
n
n
c
n
n
n
k
n
k
n
k
n
n
k
k
E
E
H
H
t
−
⋅
−
⋅
−
+
−
⋅
=
=
+
⋅
=
⋅
+
⋅
⋅
⋅
⋅
=
=
⋅
+
⋅
⋅
⋅
⋅
=
⋅
⋅
=
=
+
−
+
−
1
1
2
1
1
1
2
1
2
1
2
1
2
1
1
2
1
1
2
2
1
2
2
2
1
2
1
1
2
1
2
1
,
2
2
1
2
,
2
1
2
1
,
2
1
1
2
2
1
)
(
1
)
(
2
)
(
1
)
(
2
cos
cos
cos
2
cos
cos
cos
2
cos
cos
cos
2
2
(5.10)
Równania (5. 3)–(5.10) moŜna wykorzystać do opisu propagacji fale elektromagnetycznej
w układzie stosując formalizm macierzowy (Dodatek D).
74
V Podsumowanie
Zjawisko ujemnego załamania fali elektromagnetycznej jest zagadnieniem nowym i nie do
końca jeszcze poznanym. Sprzeczne z intuicją właściwości, jakie wykazują ośrodki
charakteryzujące się ujemnym współczynnikiem załamania, wciąŜ wzbudzają wiele
oŜywionych dyskusji w środowisku naukowym.
W ostatnich latach przedstawiono kilka sprzecznych ze sobą opinii [109]. Wydaje
się jednak, Ŝe okres najŜarliwszych dyskusji metamateriały mają juŜ za sobą. Ich
wyjątkowość potwierdzona została kilkoma wiarygodnymi eksperymentami, których
wyniki trudno zakwestionować. Do najbardziej przełomowych naleŜą doświadczenia
wykonane w sierpniu 2006 roku, których wyniki zostaną wkrótce opublikowane [88],
[107]. Po raz pierwszy efekty eksperymentu przy uŜyciu układów zawierających
metamateriały moŜna było obserwować w laboratorium gołym okiem.
Nowy rozdział w historii metamateriałów stanowi moŜliwość ich praktycznego
wykorzystania. Teraz pytanie nie brzmi juŜ „Jak wytworzyć ośrodek o ujemnym
współczynniku załamania?”, ale raczej „Jak go wykorzystać?”. Opisywane w literaturze
zastosowania obejmują zarówno obszary dość dobrze juŜ poznane jak technika
mikrofalowa, poprzez urządzenia wykorzystywane w Ŝyciu codziennym na przykład
w medycynie (MRI
24
– obrazowanie przy uŜyciu magnetycznego rezonansu jądrowego)
aŜ po pomysły z pogranicza science fiction jak powodowanie niewidzialności obiektów
[104] lub ich lewitowania [108]. O ile praktyczne zastosowanie tych ostatnich wciąŜ
jeszcze mieści się poza granicami wyobraźni przeciętnego człowieka, to ulepszenie
metody MRI przyniosłoby bardzo wymierne efekty. Przekroczenie granicy dyfrakcji
(zob. Rozdział IV.1) oznacza moŜliwość obrazowania wnętrza ciała ludzkiego
z dokładnością większą niŜ długość fali uŜywanego promieniowania, co pozwoliłby
na wykrywanie juŜ pojedynczych komórek nowotworowych.
Praca stanowi syntezę dokonań w omawianej dziedzinie. Dość szczegółowo
przedstawiono podstawowe idee leŜące u podstaw koncepcji Veselago (Rozdział I).
Scharakteryzowano stosowaną w literaturze źródłowej terminologię oraz dokonano
klasyfikacji ośrodków na dodatnie i ujemne.
Zaprezentowano zwięzły rys historyczny dotyczący technologii wytwarzania
metamateriałów, a następnie omówiono rozwiązania technologiczne z przełomu XX i XXI
wieku (Rozdział II). DuŜo uwagi poświęcono modelowaniu badanych struktur za pomocą
linii transmisyjnych. Dokładnie omówiono teorię linii transmisyjnych (Dodatek C)
i zasadność stosowania jej w tym przypadku. Praca uwzględnia najnowsze osiągnięcia
technologiczne umoŜliwiające wytwarzanie metamateriałów dla widzialnego zakresu
widma elektromagnetycznego.
Stosunkowo duŜo uwagi poświęcono zebraniu wyników eksperymentów, które
jednoznacznie potwierdziły hipotezę Veselago o istnieniu ośrodków o ujemnym
współczynniku załamania i udowodniły posiadanie przez nie przewidywanych
wyjątkowych właściwości fizycznych (Rozdział III). Zaprezentowano najnowsze
i oczekujące na publikację eksperymenty dla zakresu widzialnego.
24
MRI
−
Magnetic Resonance Imaging
75
Opisany został szereg wybranych zastosowań metamateriałów i układów je
zawierających w dzisiejszej technice (Rozdział IV). DuŜo uwagi poświęcono dokładnemu
wyjaśnieniu teorii Pendry’ego dotyczącej idealnej soczewki płaskiej i sferycznej.
Szczególny nacisk połoŜony został na zastosowania metamateriałów w widzialnym
zakresie widma elektromagnetycznego – zarówno te juŜ wynalezione jak i zupełnie nowe,
wzbudzające wielkie nadzieje naukowców.
Scharakteryzowana została polaryzacja fali elektromagnetycznej na granicy
ośrodków o przeciwnych współczynnikach załamania, co poparte zostało starannymi
rysunkami (Rozdział V). Wyprowadzono wzory Fresnela dla ośrodków dodatnich
i ujemnych, a takŜe omówiono proste wielowarstwowe układy optyczne zawierające
metamateriały. Zdefiniowano pojęcie supersieci optycznych i zwięźle przedstawiono
ich główne zalety. Omówiono takŜe sposób opisu propagacji fali elektromagnetycznej
przez układy warstwowe przy zastosowaniu formalizmu macierzowego (Dodatek D).
Wyprowadzone w rozdziale V wzory Fresnela mogą być przydatne przy numerycznym
modelowaniu propagacji fali elektormagnetycznej przez supersieci optyczne zawierające
ośrodki ujemne.
W pracy zebrano obszerny spis literatury dotyczącej ujemnego załamania fal
elektromagnetycznych. Zawiera on najwaŜniejsze, wydane ostatnio, pozycje ksiąŜkowe,
liczne publikacje naukowe powstałe zarówno w początkowym okresie rozwoju tej
dziedziny jak i oczekujące na wydanie, najnowsze odkrycia oraz adresy najbardziej
wartościowych pod względem merytorycznym stron internetowych dotyczących
metamateriałów. Z tych powodów praca moŜe być wykorzystana jako obszerny materiał
ź
ródłowy do celów naukowych i dydaktycznych.
Według najlepszej wiedzy Autorki, praca stanowi obecnie jedyne tak obszerne
opracowanie w języku polskim zagadnień z zakresu podstaw fizycznych i technologii
otrzymywania metamateriałów wykazujących zjawisko ujemnego załamania fal
elektromagnetycznych.
76
Dodatek A – Iloczyn wektorowy
Iloczynem wektorowym wektorów
]
,
,
[
z
y
x
a
a
a
=
a
r
i
]
,
,
[
z
y
x
b
b
b
=
b
r
jest wektor c
r
określony jako
]
,
,
[
]
,
,
[
x
y
y
x
z
x
x
z
y
z
z
y
z
y
x
b
a
b
a
b
a
b
a
b
a
b
a
c
c
c
−
−
−
=
=
c
r
Wektor wynikowy c
r
ma następujące właściwości:
•
jego wartość jest równa iloczynowi wartości obu wektorów wyjściowych razy sinus
kąta zawartego między nimi
)
,
sin( b
a
b
a
c
r
r
r
r
r
⋅
⋅
=
•
jest on prostopadły do płaszczyzny wyznaczonej przez wektory wyjściowe a
r
i
b
r
•
jego zwrot ustalany jest przy pomocy reguły śruby prawoskrętnej
Wartość iloczynu wektorowego jest równa iloczynowi długości pierwszego wektora przez
długość rzutu drugiego wektora na kierunek prostopadły do pierwszego wektora. Wektor
zerowy otrzymamy, gdy jeden z wektorów wyjściowych jest zerowy lub gdy wyjściowe
wektory są równoległe.
Dodatek B – Magnetyczny potencjał wektorowy
Przypuśćmy, Ŝe mamy wzdłuŜne wzbudzenie plazmowe w strukturze przedstawionej
na Rys. 8. Wektor falowy jest skierowany wzdłuŜ osi z a długość fali padającego
promieniowania jest znacznie większa niŜ stała siatki a.
Wektor indukcji elektrycznej
(
)
[
]
t
kz
i
e
D
ω
−
⋅
=
]
0,
0,
[
0
D
(B.1)
Zgodnie z prawem Ampere’a-Maxwella przepływający prąd oraz zmienne pole
elektryczne wytwarzają wirowe pole magnetyczne
j
D
H
+
∂
∂
=
×
∇
t
(B.2)
JeŜeli zarówno D jak i j mają ten sam rozkład w płaszczyźnie x-y to prawa strona równania
równa jest zero czyli nie ma tam pola magnetycznego. MoŜna się o tym przekonać
obliczając towarzyszący temu gęstość ładunku i gęstość prądu wprowadzone przez
wzdłuŜne pole.
(
)
[
]
t
kz
e
kD
ω
σ
−
⋅
=
⋅
∇
=
l
0
i
D
(B.3)
[
]
(
)
[
]
t
kz
e
D
ω
ω
−
⋅
=
i
0
0,
0,
i
j
(B.4)
JeŜeli wstawimy równania (B.3) i (B.4) do wzoru (B.2), jego prawa strona będzie
równa zero. Tak dzieje się wewnątrz kondensatora rozładowującego się poprzez
77
jednorodny dielektryczny rdzeń – pole magnetyczne nie jest generowane. W naszym
przypadku D i j mają inne rozkłady w płaszczyźnie x-y. Prąd jest ograniczony do bardzo
cienkich metalowych drutów, podczas gdy przy zastosowaniu duŜej długości fali
padającego promieniowania D jest stałe w płaszczyźnie x-y.
Dlatego teŜ w naszym przypadku pole magnetyczne nie jest zerowe. Jego wartość
w pobliŜu drutu obliczyć moŜna stosując następujące przybliŜenie: płaszczyznę x-y
dzielimy na kwadratowe komórki tak, aby w kaŜdej komórce w punkcie przecięcia jej
przekątnych znalazł się cienki metalowy drut. Aby obliczyć pole magnetyczne w otoczeniu
drutu przybliŜamy kaŜdy kwadrat prostopadłu do drutu kołem o takim samym polu
powierzchni, czyli o promieniu równym
π
a
R
k
=
. W kaŜdym punkcie płaszczyzny x-y
badane pole magnetyczne pochodzi tylko od jednego z drutów – tego, który znajduje się
najbliŜej. Pole magnetyczne wewnątrz kaŜdego z kół obliczyć moŜna
R
r
dla
0
R
R
0
dla
2
2
k
k
2
2
>
<
<
⋅
−
=
k
k
R
R
R
j
R
j
H
π
π
(B.5)
zaś odpowiadający temu potencjał wektorowy (skierowany wzdłuŜ osi z)
>
<
<
−
−
=
k
k
2
2
2
0
R
r
dla
0
R
R
0
dla
2
ln
2
)
(
k
k
k
k
R
R
R
R
R
j
R
A
π
µ
(B.6)
Gdzie wykorzystano moŜliwość wyboru stałej integracji tak Ŝe A jest zerem poza
kołem. Powód takiego wyboru był taki, Ŝe wektor A jednego koła nie zachodzi w ten
sposób na inne, więc wykluczona jest induktancja wzajemna pomiędzy drutami,
przynajmniej w tym przybliŜeniu.
MoŜemy więc obliczyć wartość A
+
−
=
−
−
=
2
1
2
ln
2
2
ln
2
)
(
2
2
0
2
2
2
0
a
r
a
r
j
R
R
r
R
r
j
r
A
k
k
k
k
π
π
π
µ
π
µ
(B.7)
Przyjęliśmy
≈
a
r
j
r
A
k
ln
2
)
(
0
π
µ
, poniewaŜ
51
,
8
10
5
10
ln
ln
3
6
−
=
×
=
−
−
a
r
. Resztkowe
pole magnetyczne wynikające z zastosowania przybliŜenia komórek sześciennych kołami,
jest znikome poniewaŜ kwadrat ma poczwórną symetrię, więc korekta pola magnetycznego
byłaby na poziomie R
-4
.
78
Dodatek C – Teoria linii transmisyjnych
W obwodzie elektrycznym długość przewodów łączących ze sobą poszczególne elementy
obwodu nie ma znaczenia tylko w uproszczonych przypadkach, gdy do czynienia mamy
ze stałym napięciem na całej długości przewodu. JeŜeli jednak napięcie w obwodzie
zmienia się z okresem porównywalnym z czasem, jaki potrzebuje sygnał by pokonać
długość przewodu – wpływu jego długości nie da się pominąć i wtedy przewód taki
traktować naleŜy jako linię transmisyjną. Przewód naleŜy traktować jako linię transmisyjną
Mówiąc ogólnie – przewód naleŜy traktować jak linię transmisyjną, jeŜeli obwód zawiera
elementy częstotliwościowe operujące na długości fali porównywalnej z długością drutu.
Linia transmisyjna jest to ośrodek lub struktura słuŜąca do transmisji energii
w róŜnej postaci (np. fali elektromagnetycznej, fali akustycznej, mocy elektrycznej)
z jednego miejsca do drugiego. MoŜe stanowić całość lub część obwodu. Linie
transmisyjne budowane są z drutów, kabli współosiowych, płytek dielektryka, włókien
optycznych lub prowadnic falowych. Opisuje je para równań róŜnicowych zwanych
równaniami telegrafistów opisanych przez Oliviera Heaviside’a – twórcę modelu linii
transmisyjnej. Zasadniczą ideę linii transmisyjnych moŜna ująć w ten sposób,
Ŝ
e przewodnik złoŜony jest z nieskończonej liczby małych segmentów, analogicznych jak
ten na rysunku Rys. 45. Równania telegrafistów moŜna traktować jako uproszczony
przypadek równań Maxwella.
Rys. 45
Schemat elementarnej komórki linii transmisyjnej.
Rozproszona rezystancja R przewodnika reprezentowana jest przez równoległy
rezystor R’, którego wartość wyraŜa się w ohmach na jednostkę długości;
Rozproszona induktancja L (wynikająca z obecności pola magnetycznego wokół
przewodnika z prądem oraz z samoindukcji) reprezentowana jest przez szeregową
cewkę L’ (henry na jednostkę długości);
Pojemność C pomiędzy dwoma przewodnikami reprezentowana jest przez
równoległy kondensator C’ (farad na jednostkę długości);
Konduktancja G dielektryka pomiędzy przewodnikami reprezentowana jest przez
równoległą konduktancję G’ pomiędzy drutem przesyłowym i odbiorczym
(siemens na jednostkę długości).
79
W przypadku, gdy R i G są bardzo małe ich efekt oddziaływania na układ moŜe być
zaniedbany i wtedy linię transmisyjną uznajemy za idealną i bezstratną. W tym przypadku
model opiera się tylko na elementach L i C i otrzymujemy parę równań róŜnicowych –
pierwsze opisujące napięcie V, a drugie prąd I wzdłuŜ linii w zaleŜości od połoŜenia x
i czasu t
( )
( )
t
x
I
t
L
t
x
V
x
,
,
∂
∂
−
=
∂
∂
( )
( )
t
x
V
t
C
t
x
I
x
,
,
∂
∂
−
=
∂
∂
PowyŜsze równania mogą być przekształcone do postaci dwóch równowaŜnych równań
falowych
V
t
LC
V
t
2
2
2
2
1
∂
∂
=
∂
∂
I
t
LC
I
t
2
2
2
2
1
∂
∂
=
∂
∂
W przypadku statycznym równania redukują się do
( )
( )
0
2
2
2
=
⋅
+
∂
∂
x
V
LC
x
x
V
ω
( )
( )
0
2
2
2
=
⋅
+
∂
∂
x
I
LC
x
x
I
ω
gdzie ω jest zadaną częstością
JeŜeli linia jest nieskończenie długa lub gdy jest zakończona swoją impedancją
charakterystyczną, równania te opisują falę elektromagnetyczną poruszającą się
z prędkością
LC
c
1
=
. JeŜeli rozwaŜymy współosiową linię transmisyjną wykonaną
z idealnego przewodnika i próŜni jako dielektryka, wspomniana prędkość jest prędkością
ś
wiatła.
Kiedy R i G nie moŜna zaniedbać równania telegrafistów opisujące pojedynczą
komórkę mają postać
( )
( )
( )
t
x
RI
t
x
I
t
L
t
x
V
x
,
,
,
−
∂
∂
−
=
∂
∂
( )
( )
( )
t
x
GV
t
x
V
t
C
t
x
I
x
,
,
,
−
∂
∂
−
=
∂
∂
RóŜniczkując pierwsze równanie po drodze, a drugie po czasie i dokonując kilku
przekształceń algebraicznych otrzymujemy równania róŜniczkowe z jedną niewiadomą
(
)
GRV
V
t
GL
RC
V
t
LC
V
x
+
∂
∂
+
+
∂
∂
=
∂
∂
2
2
2
2
(
)
GRI
I
t
GL
RC
I
t
LC
I
x
+
∂
∂
+
+
∂
∂
=
∂
∂
2
2
2
2
80
Równania moŜna traktować jak jednorodne równania falowe z tym zastrzeŜeniem,
Ŝ
e tłumienie w nich objawia się obecnością dodatkowych czynników przy V i I.
PowyŜsze
równania
falowe
sugerują
dwa
moŜliwe
rozwiązania
dla przemieszczającej się wzdłuŜ linii transmisyjnej fali
( )
(
)
(
)
kx
t
f
kx
t
f
t
x
V
+
+
−
=
ω
ω
2
1
,
gdzie
v
LC
k
ω
ω
=
=
to liczba falowa (wyraŜona w radianach na metr)
ω
jest częstością kołową (wyraŜoną w radianach na sekundę)
1
f
,
2
f
to dowolna funkcja
LC
v
1
=
prędkość propagacji
Parametr f
1
reprezentuje falę poruszającą się w prawo (w kierunku dodatnich x), zaś f
2
reprezentuje falę poruszającą się w lewo (w kierunku ujemnych x). Napięcie sumaryczne
w dowolnym punkcie linii jest sumą napięć pochodzących od tych dwóch fal.
PoniewaŜ prąd I jest związany z napięciem V poprzez równania telegraficzne, moŜna
zapisać
( )
(
)
(
)
0
2
0
1
,
Z
kx
t
f
Z
kx
t
f
t
x
I
+
−
−
=
ω
ω
gdzie Z
0
jest impedancją charakterystyczną danej linii transmisyjnej, która dla linii
bezstratnej dana jest jako
C
L
Z
=
0
Linia transmisyjna zwarta impedancją charakterystyczną nie będzie posiadała fal stojących
ani odbiciowych, zaś stosunek napięcia do prądu przy zadanej częstotliwości będzie
wartością stałą na całej długości linii. Impedancja charakterystyczna dla liniowego,
jednorodnego, izotropowego ośrodka dielektrycznego dana jest relacją
µ
ε
ε
µ
c
c
Z
=
=
=
1
0
gdzie Z
0
– to impedancja charakterystyczna
ε
– przenikalność elektryczna ośrodka (wyraŜona w faradach na metr)
µ
– przenikalność magnetyczne ośrodka (wyraŜona w henrach na metr)
µε
1
=
c
to prędkość fali w ośrodku
81
JeŜeli rozwaŜanym ośrodkiem jest próŜnia c jest prędkością światła w próŜni rozumianą
jako
0
0
1
ε
µ
=
c
i wtedy impedancja charakterystyczna to
0
0
0
0
ε
µ
µ
=
=
c
Z
, gdzie µ
0
–
stała magnetyczna (przenikalność magnetyczna próŜni), ε
0
– stała elektryczna
(przenikalność elektryczna próŜni).
UŜywając zapisu zgodnego z modelem linii transmisyjnych, ogólne wyraŜenie
na impedancję charakterystyczną ma postać
C
j
G
L
j
R
Z
ω
ω
+
+
=
0
i dla linii bezstratnych
to wspomniane juŜ
C
L
Z
=
0
(bo R i G są pomijalnie małe)
Dla rzeczywistych linii transmisyjnych mamy dwa przypadki:
Dla niskich częstotliwości czyli dla
R
L
<<
ω
oraz
G
C
<<
ω
mamy
G
R
Z
=
0
Dla wysokich częstotliwości czyli dla
R
L
>>
ω
oraz
G
C
>>
ω
mamy
C
L
Z
=
0
Są dwa rodzaje moŜliwych charakterystyk dla linii transmisyjnych. Zazwyczaj G
jest bardzo znikome, więc impedancja charakterystyczna przy niskiej częstotliwości ma
duŜą wartość, zaś dla wysokiej częstotliwości jest niewielka. Punktami przełomowymi dla
zaleŜności impedancja – częstość jest
C
G
=
1
ω
oraz
L
R
=
2
ω
. JeŜeli
C
L
G
R
>>
to oczywistym jest, Ŝe
1
2
ω
ω
>>
. Pomiędzy tymi dwoma częstotliwościami impedancja
charakterystyczna w kablu zmienia się jednostajnie.
82
Dodatek D – Formalizm macierzowy
Propagację fali elektromagnetycznej w wielowarstwowym ośrodku dielektrycznym
o J naprzemiennie ułoŜonych warstwach dwóch typów o róŜnych co do wartości i co do
znaku współczynnikach załamania opisać najłatwiej moŜna posługując się formalizmem
macierzowym [106]
=
Γ
Γ
Γ
Γ
=
Γ
∏
=
+
J
j
j
j
j
in
D
P
D
1
1
,
1
,
22
21
12
11
,
gdzie
1
,
,...,
2
,
1
,
0
+
=
J
J
j
to numery kolejnych ośrodków
Dzięki macierzy charakterystycznej moŜna obliczyć energetyczne współczynniki
transmisji
(transmitancja
ℑ
)
i
odbicia
(reflektancja
)
ℜ
dla
ośrodka
wielowarstwowego.Poszczególne składowe macierzy charakterystycznej to macierz
propagacji w ośrodku j-tym P
j
oraz macierz transmisji D
j,j+1
z ośrodka j do j+1. Definiuje
się je następująco
Macierz propagacji P
j
=
j
j
i
i
j
e
e
P
φ
φ
0
0
,
gdzie
j
x
j
j
j
j
j
j
j
j
k
d
c
n
d
n
d
,
cos
cos
2
⋅
=
=
=
θ
ω
θ
λ
π
φ
,
jest grubością fazową j-tej warstwy układu. Dla warstw dodatnich zakładamy
bezdyspersyjność współczynnika załamania n
j
i jest on z góry ustalony, natomiast
dla warstw ujemnych określamy go w sposób opisany
Macierz transmisji z ośrodka j do j+1
=
+
+
+
+
1
1
1
1
,
1
,
1
,
1
,
j
j
j
j
j
j
j
j
r
r
t
D
.
Elementami macierzy transmisji są amplitudowe współczynniki odbicia i transmisji równe
odpowiednio Transmitancję
ℑ
i reflektancję
ℜ
dla wielowarstwowego ośrodka moŜna
obliczyć ze wzorów
2
11
1
1
3
3
1
cos
cos
Γ
=
ℑ
Γ
θ
θ
n
n
,
2
11
21
Γ
Γ
=
ℜ
Γ
.
83
JeŜeli, tak jak w rozpatrywanym tu przypadku, ośrodek wejściowy i wyjściowy są
identyczne (n
1
= n
3
), układ taki ma unimodularną macierz charakterystyczną
Γ
czyli
1
det
=
Γ
. Wtedy wzór transmitancja wyraŜa się jako
2
11
1
Γ
=
ℑ
Γ
.
Formalizm śladów i antyśladów macierzy charakterystycznej
W takim przypadku transmitancję moŜna wyrazić takŜe uŜywając formalizmu śladów
i antyśladów jako
2
2
4
Γ
Γ
Γ
+
=
ℑ
σ
τ
,
gdzie
22
11
Γ
+
Γ
=
Γ
τ
to ślad macierzy, natomiast
22
11
Γ
−
Γ
=
Γ
σ
jest antyśladem diagonalnym,
który dla niediagonalnych wyrazów macierzy charakterystycznej
Γ
ma postać
12
21
Γ
−
Γ
=
Γ
ς
antysymetryczny antyślad niediagonalny
12
21
Γ
+
Γ
=
Γ
η
symetryczny antyślad niediagonalny
Ś
lady macierzy i symetryczne antyślady niediagonalne mają tylko część
rzeczywistą, zaś antyślady diagonslnr i antysymetryczne antyślady diagonalne tylko część
urojoną.
D
la aperiodycznego układu wielowarstwowego złoŜonego z dwu typów warstw A i B
macierz charakterystyczna
Γ
przyjmuje postać
out
A
A
in
D
Q
D
,
,
=
Γ
,
gdzie Q jest unimodularną macierzą charakterystyczną opisującą propagację fali
elektromagnetycznej w wielowarstwowym ośrodku aperiodycznym umieszczonym między
dwoma jednorodnymi ośrodkami typu A. Macierz Q jest iloczynem macierzy
A
Q
i
B
Q
opisujących propagację fali EM w kaŜdej z warstw
=
=
−
A
A
i
i
A
A
e
e
P
Q
φ
φ
0
0
=
=
−
1
1
1
0
0
1
1
1
BA
BA
BA
i
i
AB
AB
AB
BA
B
AB
B
r
r
t
e
e
r
r
t
D
P
D
Q
A
A
φ
φ
84
Bibliografia
[1]
www.wave-scattering.com/negative.html
;
[2]
V.G.Veselago, Usp. Fiz. Nauk. 92, 517–526 (1967);
[3]
http://zhurnal.ape.relarn.ru/~vgv/
;
[4]
V.G.Veselago, “The electrodynamics of substances with simultaneously negative
values of
ε
and
µ
”, Sov. Phys. Usp. 10 509–514 (1968);
[5]
J.B.Pendry, A.J.Holden, D.J.Robbins, I.Youngs, “Extremely low Frequency
Plasmons in Metallic Mesostructures”, Phys. Rev. Lett. 76 4773–4776 (1996)
[6]
J.B.Pendry, A.J.Holden, D.J.Robbins, W.J.Stewart, “Low frequency plasmons in
thin-wire structures”, J.Phys.: Condens. Matter 10 4785–4809 (1998);
[7]
J.B.Pendry, A.J.Holden, D.J.Robbins, W.J.Stewart, “Magnetism from conductors and
enhanced nonlinear phenomena”, IEEE Trans. Microwave Theory Tech. 47 2075–
2084 (1999);
[8]
D.Halliday, R.Resnick, J.Walker, “Podstawy Fizyki T 2–4”, wydanie pierwsze, PWN
Warszawa (2003);
[9]
R.P.Feynman, R.B.Leighton, M.Sands, “Feynmana wykłady z fizyki T.2, cz.1”,
wydanie piąte PWN Warszawa (2004);
[10]
M.Born, E.Wolf, “Principles of optics”, wydanie siódme rozszerzone, Pergamon
Press, London 1999;
[11]
A.Lakhtakia, M.W. McCall, W.S.Weiglhofer, J.Gerardin, J.Wang, “On mediums
with negative phase velocity: a brief overview”, arXiv:physics/0205027 v1 (2002);
[12]
A.Lakhtakia, wybrane artykuły na temat naturalnej aktywności optycznej (Milestone
Volume 15) SPIE Optical Engineering Press, Bellingham, WA, USA (1990);
[13]
I.V.Lindell, S.A.Tretyakov, K.I.Nikoskinen, S.Ilvonen, „BW media – media with
negative parameters, capable of supporting backward waves”, Microw. Opt.
Technol. Lett. 31 129-133 (2001);
[14]
R.W.Ziolkowski, E.Heyman, “Wave propagation in media having negative
permittivity and permeability”, Phys. Rev. E 64 056625 (2001);
[15]
P.Yeh, “Optical Waves in Layered Media”, Rockwell International Science Center,
Thousand Oaks, California (1988);
[16]
http://home.agh.edu.pl/~kakol/efizyka_pl.htm
;
[17]
A.Figotin, I.Vitebskiy, “Slow light in photonic crystals” arXiv:physics/0504112 v3
(2002);
[18]
http://krypton.mnsu.edu/~7364eb/Math113/groupvelocity.html
;
[19]
http://ocw.mit.edu/index.html
;
[20]
http://www.isvr.soton.ac.uk/SPCG/Tutorial/Tutorial/Tutorial_files/Web-further-
dispersive.htm
;
[21]
J.B.Pendry, D.R.Smith, “Metamorfoza soczewki”, Świat Nauki, nr8 (180), s.46-53,
(2006);
[22]
J.Q.Shen, “Introduction to the theory of left-handed media”, arXiv:cond-
mat/0402213 v1 (2004);
85
[23]
H.Lamb, “On group velocity”, Proc. London Math. Soc. vol.1, 473-479, (1904);
[24]
A.Schuster, “An introduction to the theory of optics”, Edward Arnold, London, 313-
318, (1905);
[25]
H.C.Pocklington, “Growth of a wave-group when the group velocity is negative”,
Nature 71, 607-608, (1905);
[26]
L.I.Mandel’shtam, “Lectures on certain problems in the theory of oscillations”,
(1944) tłumaczenie dostępne na stronie
http://ece-www.colorado.edu/~kuester/
;
[27]
L.I.Mandel’shtam, “Group velocity in a crystal lattice”, (1945) tłumaczenie j.w.;
[28]
D.V.Sivukhin, “The energy of electromagnetic waves in dispersive media”, Opt.
Spektrosk. 3, 308
−
312, (1957);
[29]
W.E.Kock “Metallic delay lenses”, Bell Syst. Tech. J., vol.27, 58
−
82, (1948);
[30]
W.E.Kock, “Radio lenses”, Bell Lab Rec., vol.24, 177
−
216, (1946);
[31]
W.E.Kock, “Metal lens antennas”, Proceedings, IRE and Waves and Electrons,
828
−
836, (1946);
[32]
D.R.Smith, W.J.Padilla, D.C.Vier, S.C.Nemat-Nasser, S.Schultz, “Composite
medium with simultaneously negative permeability and permittivity”, Phys.Rev. Lett.
84
4184
−
4187 (2000);
[33]
D.R.Smith, N.Kroll, “Negative refractive index in left-handed materials”, Phys. Rev.
Lett., vol. 84, nr 14 2933
−
2936, (2000);
[34]
D.R.Smith, D.C.Vier, N.Kroll, S.Schultz, “Direct calculation of permeability and
permittivity for a left-handed metamaterial”, Appl. Phys. Lett., vol.77, nr 14,
2246
−
2248, (2002);
[35]
R.A.Shelby, D.R.Smith, S.C.Nemat-Nasser, S.Schultz, “Microwave transmission
through a two-dimensional, isotropic, left-handed metamaterial”, Appl. Phys. Lett.,
vol.78, nr 4, 489
−
491, (2004);
[36]
M.Bayindir, K.Aydin, E.Ozbay, P.Markoš, C.M.Soukoulis, “Transmission properties
of composite metamaterials in free space”, Appl. Phys. Lett. 81 120–122 (2002);
[37]
K.Li, S.J.McLean, R.B.Greegor, C.G.Parazzoli, M.H.Tanielian, “Free-space focused-
beam characterization of left-handed materials”, Appl. Phys. Lett. 82 2535–2537
(2003);
[38]
http://www.plasmas.org/
;
[39]
S. B. Cohn, “Analysis of the metal-strip delay structure for microwave lenses” J.
Appl. Phys., vol. 20, 257–262, (1949);
[40]
S. B. Cohn, “Experimental verfication of the metal-strip delay-lens theory” J. Appl.
Phys., vol. 24, no. 7,839–841, (1953);
[41]
J. Brown, “Artificial dielectrics” in Progress in dielectrics, vol. 2, 195–225, (1960);
[42]
P.Ikonen, “Artificial dielectrics and magnetics in microwave engineering: A brief
historical
revision”,
http://www.tkk.fi/Yksikot/Sahkomagnetiikka/kurssit/S-
96.4620/reports/artificial_history_pekka.pdf
;
[43]
W.Rotman, “Plasma simulation by artificial dielectrics and parallel-plate media”,
IRE Trans. Antennas Propag., vol.AP-10, nr 10, 82-85, (1962);
[44]
http://www.ifm.liu.se/applphys/sensor/spr.html
;
[45]
http://www.uni-oldenburg.de/biochemie/11906.html
;
86
[46]
D.F.Sievenpiper, M.E.Sickmiller, E.Yablonovitch, “3D wire mesh photonic
crystals”, Phys. Rev. Lett. 76 2480–2483 (1996);
[47]
E.Yablonovitch, T.J.Gmitter, K.M.Leung, “Photonic band structure: The face-
centered-cubic case employing nonspherical atoms” Phys. Rev. Lett. 67 2295–2298
(1991);
[48]
E.Yablonovitch,
“Photonic band-gap crystals” J.Phys.: Condens. Matter 5 2443–
2460 (1993);
[49]
E.Yablonovitch, “Photonic band-gap structures” JOSA B 10 283 (1993);
[50]
J.B.Pendry, “Calculating photonic band structure”, J.Phys.: Condens. Matter 8
1085
−
1108 (1996);
[51]
http://physics.ucsd.edu/~dav/animae.html
;
[52]
http://www.ifh.ee.ethz.ch/~martin/
;
[53]
G.V.Eleftheriades, O.Siddiqui, A.K.Iyer, “Transmission line models for negative
refractive index media and associated implementations without excess resonators”,
IEEE Microwave Wireless Components Lett., vol.13, nr 2, 51
−
53, (2003);
[54]
R.A.Shelby , D.R.Smith , S.Schultz, “Experimental Verification of Negative Index of
Refraction”, Science 292 77–79 (2001);
[55]
N.Garcia, M.Nieto-Vesperinas, „Is there an experimental verification of negative
index of refraction yet?”, Opt. Lett. Vol.27, nr 11, 885
−
887 (2002);
[56]
P.M.Valanju, R.M.Walser, A.P.Valanju, “Wave Refraction in Negative-Index Media:
Always Positive and Very Inhomogeneous”, Phys.Rev.Lett. vol.88, 187401, (2001);
[57]
J.M.Williams, „Some problems with negative refraction”, Phys. Rev. Lett. vol.87,
249703, (2001);
[58]
A.A.Houck, J.B.Brock, I.L.Chuang, “Experimental observations of a left-handed
material that obeys Snell’s law”, Phys.Rev.Lett. vol.90, nr 13, 137401, (2003);
[59]
C.G.Parazzoli, R.B.Greegor, K.Li, B.E.C.Koltenbah, M.Tanielian,“Experimental
verification and simulation of negative index of refraction using Snell's law”, Phys.
Rev. Lett. vol.90, 107401, (2003);
[60]
C.G.Parazzoli, R.B.Greegor, K.Li, B.E.C.Koltenbah, M.Tanielian,“Experimental
determination and numerical simulation of the properties of negative index of
refraction materials”, Opt.Ex. vol.11, nr 7, s.688
−
695, (2003);
[61]
R.B.Greegor, C.G.Parazzoli, K.Li, M.Tanielian, “Origin of dissipative losses in
negative index of refraction materials”, Appl.Phys.Lett. vol.82, nr 14, 2356
−
8,
(2002);
[62]
G. Kron, „Equivalent circuit of the field equations of Maxwell”, Proc. IRE, vol. 32,
nr 5, s. 289
−
299 (1944);
[63]
G.Kron, “Numerical solution of ordinary and partial differential equations by
means of equivalent circuits”, General Electric Company, Schenectady, New York
(1944);
[64]
J.R.Whinnery, S.Ramo, “A new approach to the solution of high-frequency field
problems”, Proc. IRE, vol. 32, nr 5, s.284
−
288 (1944);
[65]
G.V.Eleftheriades,
K.G.Balmain,
“Nagative-refraction
metamaterials
–
Fundamental Principles and Applications”, IEEE Press (2005);
[66]
A.K.Iyer, G.V.Eleftheriades, “Negative refractive index materials supporting 2-D
waves”, IEEE MTT-S Internationam Microwave Symposium Digest, vol. 2, 1067–70
(2002);
87
[67]
G.V.Eleftheriades, A.K.Iyer, P.C.Kremer, “Planar negative refractive index media
using periodically L-C loaded transmission lines”, IEEE Trans. Microwave Theory
Tech., vol.50, no.12, pp.2702-2712, (2002);
[68]
C.Caloz, H.Okabe, H.Iwai, T.Itoh, “Transmission line approach of left-handed
materials”, USNC/URSI National Radio Science Meeting Digest (2002);
[69]
C.Caloz, T.Itoh, “Novel microwave devices and structures based on the
transmission line approach of metamaterials”, IEEE MTT-S Internationam Microwave
Symposium Digest, vol. 1, pp.195 – 198 (2003);
[70]
A.A.Oliner, “A planar negative-refractive-index medium without resonant
elements”, IEEE MTT-S Internationam Microwave Symposium Digest, vol. 1, pp.191
– 194 (2003);
[71]
V.A.Podolskiy, A.K.Sarychev, V.M.Shalaev, “Plasmon modes in metal nanowires
and lefthanded materials”, J. Nonlinear Opt. Phys. Materials 11, 65 (2002);
[72]
V.A.Podolskiy, A.K.Sarychev, V.M.Shalaev, “Plasmon modes and negative
refraction in metal nanowire composites”, Opt.Ex. vol. 11, nr 7, 735, (2003);
[73]
E.M.Purcell, C.R.Pennypacker, “Scattering and absorption of light by nonspherical
dielectric grains”, Astrophys. J. 186, 705, (1973);
[74]
A.Berrier, M.Mulot, M.Swillo, M.Qiu, L.Thylen, A.Talneau, S.Anand, „Negative
refraction at infrared wavelengths in a two-dimensional photonic crystal”, Phys.
Rev. Lett. 93, 073902 (2004);
[75]
E.Schonbrun, M.Tinker, W.Park, J.B.Lee, “Negative Refraction in a Si-Polymer
Photonic Crystal Membrane”, IEEE Phot.Tech.Lett., vol. 17, nr 6, (2005);
[76]
T.J.Yen, W.J.Padilla, N.Fang, D.C.Vier, D.R.Smith, J.B.Pendry, D.N.Basov,
X.Zhang, “Terahertz Magnetic Response from Artificial Materials”, Science 303,
1494-1496, (2004);
[77]
S.Linden, C.Enkrich, M.Wegener, J.Zhou, T.Koschny, C.M.Soukoulis, “Magnetic
Response of Metamaterials at 100 Terahertz”, Science 306, 1351-1353, (2004);
[78]
S.Zhang, W.Fan, B.K.Minhas, A.Frauenglass, K.J.Malloy, S.R.J.Brueck, “Mid-
infrared resonant magnetic nanostructures exhibiting a negative permeability”,
Phys. Rev. Lett. vol.94, 037402 (2005);
[79]
N.-C.Panoiu, R.M.Osgood Jr., “Influence of the dispersive properties of metals on
the transmission characteristics of left-handed materials”, Phys. Rev. E 68, 016611
(2003);
[80]
V.M.Shalaev, W.Cai, U.Chettiar, H.-K.Yuan, A.K.Sarychev, V.P.Drachev,
A.V.Kildishev, „Negative index of refraction In optical metamaterials”, Opt.Lett.
vol.30, nr 24, s.3356-8, (2005);
[81]
A.N.Lagarkov, A. K.Sarychev, “Electromagnetic properties of composites
containing elongated conducting inclusion”, Phys. Rev. B 53, 6318–6336 (1996);
[82]
L.V.Panina, A.N.Grigorenko, D. P. Makhnovskiy, “Optomagnetic composite medium
with conducting nanoelements”, Phys. Rev. B 66, 155411 (2002);
[83]
V.A.Podolskiy, A.K.Sarychev, E.E.Narimanov, V.M.Shalaev, “Resonant light
interactions with plasmonic nanowire systems”, J.Opt.A:Pure Appl.Opt. vol. 7,nr 2,
S32-S37, (2005);
[84]
G.Dolling, C.Enkrich, M.Wegener, J.Zhou, C.M.Soukoulis, S.Linden, “Cut-wire
pairs and plate pairs as magnetic atoms for optical metamaterials”, Opt.Lett., vol.
30
, nr 23, s.3198-3200, (2005);
88
[85]
V.P.Drachev, W.Cai, U.Chettiar, H.-K.Yuan, A.K.Sarychev, A.V.Kildishev,
G.Klimeck, V.M.Shalaev, „Experimental verification of an optical negative-index
material”, Laser Phys.Lett. vol.3, nr 1, s.49-55, (2006);
[86]
S.Zhang,
W.Fan,
N.C.Panoiu,
K.J.Malloy,
R.M.Osgood,
S.R.J.Brueck,
„Demontsration of near-infrared negative-index-materials”, Phys.Rev.Lett., vol. 95,
nr 13, 137404, (2005);
[87]
D.R.Smith, S.Shultz, P.Markos, C.M.Soukoulis, „Determination of negative
permittivity and permeability of metamaterials from reflection and transmission
coefficients”, Phys.Rev.B vol.65, 195104, (2002);
[88]
G.Dolling, M.Wegener, C.M.Soukoulis, S.Linden, “Negative-index metamaterial at
780 nm wavelength”,
http://arXiv:physics/0607135
, (2006);
[89]
A.N.Grigorenko, A.K.Geim, H.F.Gleeson, Y.Zhang, A.A.Firsov, I.Y.Khrushchev,
J.Petrovic, „Nanofabricated media with negative permeability at visible
frequencies”, Nature 438, nr 7066, s.335 (2005);
[90]
I.V.Shadrivov, N.A. Zharova, A.A.Zharov, Y.S.Kivshar, „Defect modes and
transmission properties of left-handed bandgap structures”, Phys. Rev. E vol. 70,
(2004);
[91]
http://www.waves.utoronto.ca/prof/gelefth/publications.html
[92]
A. A. Zharov, I. V. Shadrivov, and Yu. S. Kivshar, “Nonlinear Properties of Left-
Handed Metamaterials”, Phys. Rev. Lett. 91, 037401 (2003)
[93]
J.B.Pendry, “Negative refraction makes a perfect lens”, Phys.Rev.Lett. vol. 85, nr
18, 3966-9 (2000);
[94]
G.W.’t Hooft, “Comment on <Negative refraction makes a perfect lens>”,
Phys.Rev.Lett. vol. 87, nr 24, 249701, (2001);
[95]
N.Garcia, M.Nieto-Vesperinas, „Left-handed materials do not make a perfect lens”,
Phys.Rev.Lett. vol. 88 nr 20, 207403, (2002);
[96]
Long Gen Zheng, Wen Xun Zhang, „Discussion on Negative Refraction and Perfect
Lens”, Progress In Electromagnetics Research Symposium 2005, Hangzhou,
(23-26.08.2005);
[97]
R.F.Broas, D.F.Sieverpiper, E.Yablonovitch, “A high-impedance ground plane
applied to a cell phone handset geometry”, IEEE Trans. Micr. Theory and Tech.
vol.49, 1262-1265, (2001);
[98]
M.C.K.Wiltshire, J.B.Pendry, I.R.Young, D.J.Larkman, D.J.Gilderdale, J.V.Hajnal,
“Microstructured magnetic materials for RF flux guides in magnetic resonance
imaging”, Science 291, 848-851, (2001);
[99]
Chiyan Luo, Steven G. Johnson, J.D.Joannopoulos, J.B.Pendry, “All-angle negative
refraction without negative effective index”, Phys. Rev. Rapid Communications B65,
201104®, (2002);
[100]
J.B.Pendry, “Perfect cylindrical lenses”, Opt.Ex., vol. 11, nr 7, 755, (2003);
[101]
M.A.Antoniades, G.V.Eleftheriades, “Compact linear lead/lag metamaterial phase
shifters for broadband applications”, IEEE Antennas and wireless propagation
letters, vol.2, (2003);
[102]
A.Grbic, G.V.Eleftheriades, “Experimental verification of backward-wave radiation
from a negative refractive index metamaterial”, J.Appl.Phys. vol.92, nr10, 5930-5,
(2002);
89
[103]
G.V.Eleftheriades, “Enabling RF/Microwave devices using negative-refractive index
transmission-line metamaterials”, Radio Science Bulletin nr 312, (2005);
[104]
J.B.Pendry, D.Schuring, D.R.Smith, “Controlling electromagnetic fields”, Science
vol. 312, 1780 (2006);
[105]
R.A.Shelby, D.R.Smith, S.C.Nemat-Nasser, S.Schultz, “Microwave transmission
through a two-dimensional, isotropic, left-handed metamaterial”, Appl. Phys.
Lett. 78, 489–491 (2001);
[106]
A.Klauzer-Kruszyna, „Propagacja światła spolaryzowanego w wybranych
supersieciach aperiodycznych”, praca doktorska przygotowana pod kierunkiem dr
hab. Włodzimierza Salejdy, prof. nadzw. PWr ;
[107]
J.Zhou, T.Koschny, L.Zhang, G.Tuttle, C.M.Soukoulis, „Experimental verification of
negative index of refraction”,
http://arXiv.org/physics/0608301
;
[108]
U.Leonhardt, T.G.Philbin, “Quantum optics of special transformation media”,
http://arXiv.org/quant-ph/0608115
.
[109]
Bliokh K.Yu., Bliokh Yu.P.
What are the left-handed media and what is interesting
about them?
Uspekhi Fizicheskikh Nauk, No 4, 439. (Methodological Notes).