Download - Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Transcript
Page 1: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Complete Band Gaps:

You can leave home without them.

Photonic Crystals:Periodic Surprises in Electromagnetism

Steven G. Johnson

MIT

Page 2: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

How else can we confine light?

Page 3: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Total Internal Reflection

ni > no

no

rays at shallow angles > θc

are totally reflected

Snell’s Law:

θi

θo

ni sinθi = no sinθo

sinθc = no / ni

< 1, so θc is real

i.e. TIR can only guidewithin higher indexunlike a band gap

Page 4: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Total Internal Reflection?

ni > no

no

rays at shallow angles > θc

are totally reflected

So, for example,a discontiguous structure can’t possibly guide by TIR…

the rays can’t stay inside!

Page 5: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Total Internal Reflection?

ni > no

no

rays at shallow angles > θc

are totally reflected

So, for example,a discontiguous structure can’t possibly guide by TIR…

or can it?

Page 6: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Total Internal Reflection Redux

ni > no

no

ray-optics picture is invalid on λ scale (neglects coherence, near field…)

Snell’s Law is reallyconservation of k|| and ω:

θi

θo|ki| sinθi = |ko| sinθo

|k| = nω/c

(wavevector) (frequency)

k||

translationalsymmetry

conserved!

Page 7: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Waveguide Dispersion Relationsi.e. projected band diagrams

ni > no

no

k||

ω

light line: ω

= ck / no

light coneprojection of all k⊥ in no

(a.k.a. β)

ω = ck / ni

higher-index corepulls down state

( ∞)

higher-order modesat larger ω, β

weakly guided (field mostly in no)

Page 8: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

J

J

J

J

J

J

J

JJ

J J

0

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

freq

uenc

y (c

/a)

wavenumber k (2π/a)

light cone

Conserved k and ω+ higher index to pull down state

= localized/guided mode.

Strange Total Internal ReflectionIndex Guiding

a

Page 9: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

A Hybrid Photonic Crystal:1d band gap + index guiding

J

J

J

J

J

J

J

JJ

J J

0

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

freq

uenc

y (c

/a)

wavenumber k (2π/a)

light cone

“band gap”

a

range of frequenciesin which there are

no guided modes

slow-light band edge

Page 10: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

A Resonant Cavity

index-confined

photonic band gap

increased rod radiuspulls down “dipole” mode

(non-degenerate)

– +

Page 11: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

J

J

J

J

J

J

J

JJ

J J

0

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

freq

uenc

y (c

/a)

wavenumber k (2π/a)

light cone

“band gap” ω

A Resonant Cavity

index-confined

photonic band gap

– +

k not conservedso coupling to

light cone:

radiationThe trick is to

keep theradiation small…

(more on this later)

Page 12: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Meanwhile, back in reality…

5 µm

[ D. J. Ripin et al., J. Appl. Phys. 87, 1578 (2000) ]

d = 703nmd = 632nmd

Air-bridge Resonator: 1d gap + 2d index guiding

bigger cavity= longer λ

Page 13: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Time for Two Dimensions…

2d is all we really need for many interesting devices…darn z direction!

Page 14: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

How do we make a 2d bandgap?

Most obvioussolution?

make2d patternreally tall

Page 15: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

How do we make a 2d bandgap?

If height is finite,we must couple to

out-of-plane wavevectors…

kz not conserved

Page 16: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

A 2d band diagram in 3d

J

J

J

J

J

J

J

JJJJJ

JJJJ

JJJJJJJ

J

J

J

J

J

J

J

J

JJJJJJJJJJJJJ

JJJJJJJ

JJJJJJ

JJJJJ

JJJJJJJJJJJJJJJJJJJJJJJ

JJJJJ

JJJJJJJ J

JJ

E

E

E

E

E

E

E

E

E

E

EEEEEEEEEEEE

E

E

E

E

E

E

E

E

E

EEEEEEEE

E

E

EEEEEE

E

E

E

EEE

EEEEE

EEEE

EEEE

0

0.1

0.2

0.3

0.4

0.5

0.6

freq

uenc

y (c

/a)

TE bandsTM bands

Square Lattice ofDielectric Rods

(ε = 12, r=0.2a)

Let’s start with the 2dband diagram.

This is what we’d liketo have in 3d, too!

wavevector

Page 17: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

A 2d band diagram in 3d

J

J

J

J

J

J

J

JJJJJ

JJJJ

JJJJJJJ

J

J

J

J

J

J

J

J

JJJJJJJJJJJJJ

JJJJJJJ

JJJJJJ

JJJJJ

JJJJJJJJJJJJJJJJJJJJJJJ

JJJJJ

JJJJJJJ J

JJ

E

E

E

E

E

E

E

E

E

E

EEEEEEEEEEEE

E

E

E

E

E

E

E

E

E

EEEEEEEE

E

E

EEEEEE

E

E

E

EEE

EEEEE

EEEE

EEEE

0

0.1

0.2

0.3

0.4

0.5

0.6

freq

uenc

y (c

/a)

TE bandsTM bands

Square Lattice ofDielectric Rods

(ε = 12, r=0.2a)

Let’s start with the 2dband diagram.

This is what we’d liketo have in 3d, too!

wavevector

3D Structure:

No! When we includeout-of-plane propagation,

we get:

wavevectorfrequency

ω

ω+δω

projected band diagram fills gap!

but this emptyspace looks useful…

Page 18: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Photonic-Crystal Slabs

2d photonic bandgap + vertical index guiding

[ S. G. Johnson and J. D. Joannopoulos, Photonic Crystals: The Road from Theory to Practice ]

Page 19: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Rod-Slab Projected Band Diagram

J

J

J

J

J

J

J

J

JJJJ

JJJJ

JJJJJJJ

J

J

J

J

J

J

J

J

JJ

J

J

J

J

J

J

JJJJJ

JJJJJJJJJ

J

J

J

J

J

J

JJJ

J

J

J

J

J

J

JJJJ

JJJJ

JJJJJJ

J

J

J

J

J

J

J

JJJJ

J

J

JJJJJ

JJJJJJ

JJJJJ

JJ

J

J

J

JJJ

JJ

J

J

J

J

J

JJJJ

JJJJ

JJJJJJ

J

J

J

J

J

JJ

JJJJJJ

JJ

JJJJ

JJJJ

JJJJJJ

J

JJJJJJJJ

JJJJJJJJJJ

JJJJ

JJJ

JJ

JJJJJJJJ

JJJJJJ

JJJJJJ

JJJJJJ

J

J

JJJJJJJJJJJ

JJJJ

JJJJ

JJJJJJJJ

J

JJJJJJJ

JJJJ

JJJJ

JJJJJ

JJJJ J

JJJ

JJJJ

JJ JJJ

JJJJ J

J JJ

E

E

E

E

E

E

E

E

E

EEEEE

EEEE

EEEEEE

E

E

E

E

E

E

E

EE

E

E

E

E

E

E

E

EEEEEEEEE

EEEE

E

E

E

E

E

E

EEE

E

E

E

E

E

E

EEEEE

EEEE

EEEEE

E

E

E

E

E

E

E

EEEE

E

E

E

E

EEEE

EEEEEEEEEEE

E

E

E

E

EEE

EE

E

E

E

E

EEEE

EEEEE

EEEEE

E

E

E

E

E

E

EE

EEEEEEE

EEEEE

EEE

EEE

E

EEE

EEEE

EEEEE

EEEEEEEEE

EEEEEEEEE

E

EEEEEEEE

EEEEEEEE

E

EEEEEEEEE

E

EEEEEEEEEEEEE

EEEEEEE

EEEEE E

EEE EEEEEE

E

EEEE EE

E

E

0

0.1

0.2

0.3

0.4

0.5

0.6

freq

uenc

y (c

/a)

even (TE-like) bandsodd (TM-like) bands

Square Lattice ofDielectric Rods(ε = 12, r=0.2a, h=2a)

light cone

The Light Cone:All possible statespropagating in the air

The Guided Modes:Cannot couple tothe light cone… —> confined to the slab

Thickness is critical.Should be about λ/2 (to have a gap & be single-mode)

Γ X M ΓΓ X

M

Page 20: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Symmetry in a Slab

2d: TM and TE modes

slab: odd (TM-like) and even (TE-like) modes

mirror planez = 0

rE

rE

Like in 2d, there may only be a band gapin one symmetry/polarization

Page 21: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Slab Gaps

J

J

J

J

J

J

J

J

J

JJJJJJJ

JJJJJJ

J

J

J

J

J

J

J

J

J

JJ

J

J

J

J

J

JJJJJJJ

JJJJJJ

JJJ

J

J

J

J

JJJ

J

J

J

J

J

J

J

JJJJJJJ

JJJJJJJ

J

J

J

J

J

J

J

JJJJJJ

JJ

J

JJJJJJJJJ

JJJ

J

JJJJJJJJ

JJJJJJJJJJJJ

JJJJJJJ

J

JJJ

JJJJJ

JJJJJJ

JJ

JJJJJJJ

JJJJJ

J

J

J

JJJJJ

JJJJJJJJJJJJ

JJJ J

J

JJJJJ

JJJJJJJJJ

JJJJJJ J

JJJJ

JJJ J

JJJJJJ

E

E

E

E

E

E

E

E

EEEEE

EEEEE

EEEEE

E

E

E

E

E

E

E

E

EE

E

E

E

E

E

EEEEEEEEEE

EEEEEE

E

E

E

E

E

EEE

E

E

E

E

E

E

EEEEEE

EEEEEEE

EE

E

E

E

E

E

E

EEEEEE

EE

EEEEE

EEEEE

EE

E

E

EEEE

EEEE

EEE

E

EEEE

EEEEE

EEEEEEE

E

EEEE

EEE

EEEEE

EE

EEEEEEE

EEEEE

EE

E

E

EEEEE

EEEEEE

E

EEEEE

EEE E

E

E

EEEEEEEEEEEE

E

EEEEEEEEEE E

EEE

0

0.1

0.2

0.3

0.4

0.5

J

J

J

J

J

J

J

J

JJJJ

JJJJ

JJJJJJJ

J

J

J

J

J

J

J

J

JJ

J

J

J

J

J

J

JJJJJ

JJJJJJJJJ

J

J

J

J

J

J

JJJ

J

J

J

J

J

J

JJJJ

JJJJ

JJJJJJ

J

J

J

J

J

J

J

JJJJ

J

J

JJJJJ

JJJJJJ

JJJJJ

JJ

J

J

J

JJJ

JJ

J

J

J

J

J

JJJJ

JJJJ

JJJJJJ

J

J

J

J

J

JJ

JJJJJJ

JJ

JJJJ

JJJJ

JJJJJJ

J

JJJJJJJJ

JJJJJJJJJJ

JJJJ

JJJ

JJ

JJJJJJJJ

JJJJJJ

JJJJJJ

JJJJJJ

J

J

JJJJJJJJJJJ

JJJJ

JJJJ

JJJJJJJJ

J

JJJJJJJ

JJJJ

JJJJ

JJJJJ

JJJJ J

JJJ

JJJJ

JJ JJJ

JJJJ J

J JJ

E

E

E

E

E

E

E

E

E

EEEEE

EEEE

EEEEEE

E

E

E

E

E

E

E

EE

E

E

E

E

E

E

E

EEEEEEEEE

EEEE

E

E

E

E

E

E

EEE

E

E

E

E

E

E

EEEEE

EEEE

EEEEE

E

E

E

E

E

E

E

EEEE

E

E

E

E

EEEE

EEEEEEEEEEE

E

E

E

E

EEE

EE

E

E

E

E

EEEE

EEEEE

EEEEE

E

E

E

E

E

E

EE

EEEEEEE

EEEEE

EEE

EEE

E

EEE

EEEE

EEEEE

EEEEEEEEE

EEEEEEEEE

E

EEEEEEEE

EEEEEEEE

E

EEEEEEEEE

E

EEEEEEEEEEEEE

EEEEEEE

EEEEE E

EEE EEEEEE

E

EEEE EE

E

E

0

0.1

0.2

0.3

0.4

0.5

0.6

freq

uenc

y (c

/a)

even (TE-like) bandsodd (TM-like) bands

even (TE-like) bandsodd (TM-like) bands

Γ X M Γ Γ M K Γ

Square Lattice ofDielectric Rods(ε = 12, r=0.2a, h=2a)

Triangular Latticeof Air Holes

(ε = 12, r=0.3a, h=0.5a)

light cone light cone

TM-like gap TE-like gap

Page 22: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Substrates, for the Gravity-Impaired

substrate

(rods or holes)

substrate breaks symmetry:some even/odd mixing “kills” gap

BUTwith strong confinement

(high index contrast)

mixing can be weak

superstrate restores symmetry

“extruded” substrate= stronger confinement

(less mixing evenwithout superstrate

Page 23: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Extruded Rod Substrate

S. Assefa, L. A. Kolodziejski

high index

Page 24: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Air-membrane Slabs

[ N. Carlsson et al., Opt. Quantum Elec. 34, 123 (2002) ]

who needs a substrate?

2µm

AlGaAs

Page 25: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Optimal Slab Thickness~ λ/2, but λ/2 in what material?

J

J

J

J

J

JJ

J

J

J

E

E

E

E

E

E

E

E

E

E

E

0

5

10

15

20

25

30

35

0 0.5 1 1.5 2 2.5 3

even (TE-like) gap odd (TM-like) gap

slab thickness (a)

gap

size

(%

)

TM “sees” <ε-1>-1

~ low εTE “sees” <ε>

~ high ε

effective medium theory: effective ε depends on polarization

Page 26: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Photonic-Crystal Building Blocks

point defects(cavities)

line defects(waveguides)

Page 27: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

0.32

0.34

0.36

0.38

0.4

0.42

0.44

0.3 0.35 0.4 0.45 0.5

freq

uenc

y (c

/a)

wavevector k (2π/a)

A Reduced-Index Waveguide

Reduce the radius of a row ofrods to “trap” a waveguide modein the gap.

(r=0.2a)

(r=0.18a)

(r=0.16a)

(r=0.12a)

(r=0.14a)

Still have conservedwavevector—under thelight cone, no radiation

(r=0.10a)

We cannot completelyremove the rods—novertical confinement!

Page 28: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Reduced-Index Waveguide Modes

x →

z →

y →

x →

–1 +1Ez

z →

y →x

→z

y →

x →

–1 +1Hz

z →

y →

Page 29: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Experimental Waveguide & Bend

E

EE

EEEEEEEEEEEEEEEEEEEEEEEEEE

E

EEE

EEEEEEE

EEEEEEEEEEE

EEEEEEEEEEEEEEE

EEEEEEEE

EE

E

EE

EEEEE

EEE

E

EE

EEE

EE

EEEEEEEEEEEEEEEE

EEEEE

EE

E

E

EE

E

EE

EEE

EE

EEEEEEE

EE

EEEE

0

0.2

0.4

0.6

0.8

1

1.2

1200 1300 1400 1500 1600 1700 1800

tran

smis

sion

wavelength (µm)

E experiment

theory

band gap

waveguide mode

[ E. Chow et al., Opt. Lett. 26, 286 (2001) ]

1µm 1µm

GaAs

AlO

SiO2

be

nd

ing

eff

icie

ncy

caution:can easily bemulti-mode

Page 30: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Inevitable Radiation Losseswhenever translational symmetry is broken

e.g. at cavities, waveguide bends, disorder…

k is no longer conserved!

ω(conserved)

coupling to light cone= radiation losses

Page 31: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

All Is Not Lost

Qw

A simple model device (filters, bends, …):

Qr

Q1

Qr1

Qw1= +

Q = lifetime/period = frequency/bandwidth

We want: Qr >> Qw

1 – transmission ~ 2Q / Qr

worst case: high-Q (narrow-band) cavities

Page 32: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Semi-analytical losses

r r t r r r r rE x G x x E x x

defect

( ) ( , ) ( ) ( )= ′ ⋅ ′ ⋅ ′∫ ω ε∆

far-field(radiation) Green’s function

(defect-freesystem)

near-field(cavity mode)

defect

A low-lossstrategy:

Make field insidedefect small

= delocalize mode

Make defect weak= delocalize mode

Page 33: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

(ε = 12)

Monopole Cavity in a Slab

decreasing ε

Lower the ε of a single rod: push upa monopole (singlet) state.

Use small ∆ε: delocalized in-plane, & high-Q (we hope)

Page 34: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

J

J

J

J

J

J

100

1,000

10,000

100,000

1,000,000

0.0001 0.001 0.01 0.1

Qr

∆frequency above band edge (c/a)

Delocalized Monopole Q

ε=6

ε=7

ε=8

ε=9

ε=10

ε=11

mid-gap

Page 35: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Super-defects

Weaker defect with more unit cells.

More delocalizedat the same point in the gap(i.e. at same bulk decay rate)

Page 36: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Super-Defect vs. Single-Defect Q

J

J

J

J

J

J

G

G

G

G

G

G

100

1,000

10,000

100,000

1,000,000

0.0001 0.001 0.01 0.1

Qr

∆frequency above band edge (c/a)

ε=6

ε=7

ε=8

ε=9

ε=10

ε=11

ε=7

ε=8

ε=9

ε=10

ε=11

ε=11.5

mid-gap

Page 37: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Super-Defect vs. Single-Defect Q

J

J

J

J

J

J

G

G

G

G

G

G

100

1,000

10,000

100,000

1,000,000

0.0001 0.001 0.01 0.1

Qr

∆frequency above band edge (c/a)

ε=6

ε=7

ε=8

ε=9

ε=10

ε=11

ε=7

ε=8

ε=9

ε=10

ε=11

ε=11.5

mid-gap

Page 38: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Super-Defect State(cross-section)

∆ε = –3, Qrad = 13,000

(super defect)

still ~localized: In-plane Q|| is > 50,000 for only 4 bulk periods

Ez

Page 39: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

(in hole slabs, too)

Hole Slabε=11.56

period a, radius 0.3athickness 0.5a

Reduce radius of7 holes to 0.2a

Q = 2500near mid-gap (∆freq = 0.03)

Very robust to roughness(note pixellization, a = 10 pixels).

Page 40: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

How do we compute Q?

1excite cavity with dipole source

(broad bandwidth, e.g. Gaussian pulse)

… monitor field at some point

(via 3d FDTD [finite-difference time-domain] simulation)

…extract frequencies, decay rates viasignal processing (FFT is suboptimal)

[ V. A. Mandelshtam, J. Chem. Phys. 107, 6756 (1997) ]

Pro: no a priori knowledge, get all ω’s and Q’s at once

Con: no separate Qw/Qr, Q > 500,000 hard, mixed-up field pattern if multiple resonances

Page 41: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

How do we compute Q?

2excite cavity with

narrow-band dipole source(e.g. temporally broad Gaussian pulse)

(via 3d FDTD [finite-difference time-domain] simulation)

— source is at ω0 resonance,which must already be known (via )1

…measure outgoing power P and energy U

Q = ω0 U / P

Pro: separate Qw/Qr, arbitrary Q, also get field pattern

Con: requires separate run to get ω0, long-time source for closely-spaced resonances

1

Page 42: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Can we increase Qwithout delocalizing?

Page 43: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Semi-analytical losses

r r t r r r r rE x G x x E x x

defect

( ) ( , ) ( ) ( )= ′ ⋅ ′ ⋅ ′∫ ω ε∆

far-field(radiation) Green’s function

(defect-freesystem)

near-field(cavity mode)

defect

Another low-lossstrategy:

exploit cancellationsfrom sign oscillations

Page 44: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Need a morecompact representation

Cannot cancel infinitely many E(x) integrals

Radiation pattern from localized source…

— use multipole expansion & cancel largest moment

Page 45: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Multipole Expansion[ Jackson, Classical Electrodynamics ]

+ + + …

radiated field =

dipole quadrupole hexapole

Each term’s strength = single integral over near field…one term is cancellable by tuning one defect parameter

Page 46: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Multipole Expansion[ Jackson, Classical Electrodynamics ]

+ + + …

radiated field =

dipole quadrupole hexapole

peak Q (cancellation) = transition to higher-order radiation

Page 47: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Multipoles in a 2d example

index-confined

photonic band gap

increased rod radiuspulls down “dipole” mode

(non-degenerate)

– +

as we change the radius, ω sweeps across the gap

Page 48: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

2d multipolecancellation

JJJ

J

J

J

JJ

0

5,000

10,000

15,000

20,000

25,000

30,000

0.25 0.29 0.33 0.37

Q

frequency (c/a)

r = 0.35a

Q =

1,7

73

r = 0.40a

Q =

6,6

24

r = 0.375a

Q =

28,

700

Page 49: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

cancel a dipole by opposite dipoles…

cancellation comes fromopposite-sign fields in adjacent rods

… changing radius changed balance of dipoles

Page 50: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

3d multipole cancellation?

JJ

J

J

J

JJ

J

J

J

0

500

1000

1500

2000

0.319243 0.390536

Q

frequency (c/a) gap topgap bottom

enlarge center & adjacent rods

vary side-rod ε slightlyfor continuous tuning

quadrupole mode

(Ez cross section)

(balance central moment with opposite-sign side rods)

Page 51: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

3d multipole cancellationne

ar f

ield

E zfa

r fie

ld |E

|2

Q = 408 Q = 426Q = 1925

nodal planes(source of high Q)

Page 52: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

An Experimental (Laser) Cavity[ M. Loncar et al., Appl. Phys. Lett. 81, 2680 (2002) ]

elon

gate

row

of

hole

s

Elongation p is a tuning parameter for the cavity…

cavity

…in simulations, Q peaks sharply to ~10000 for p = 0.1a

(likely to be a multipole-cancellation effect)

* actually, there are two cavity modes; p breaks degeneracy

Page 53: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

An Experimental (Laser) Cavity[ M. Loncar et al., Appl. Phys. Lett. 81, 2680 (2002) ]

elon

gate

row

of

hole

s

Elongation p is a tuning parameter for the cavity…

cavity

…in simulations, Q peaks sharply to ~10000 for p = 0.1a

(likely to be a multipole-cancellation effect)

* actually, there are two cavity modes; p breaks degeneracy

Hz (greyscale)

Page 54: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

An Experimental (Laser) Cavity[ M. Loncar et al., Appl. Phys. Lett. 81, 2680 (2002) ]

cavity

quantum-well lasing threshold of 214µW(optically pumped @830nm, 1% duty cycle)

(InGaAsP)

Q ~ 2000 observed from luminescence

Page 55: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

How can we get arbitrary Qwith finite modal volume?

Only one way:

a full 3d band gap

Now there are two ways.

[ M. R. Watts et al., Opt. Lett. 27, 1785 (2002) ]

Page 56: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

The Basic Idea, in 2d

ε1

ε2 < ε1

ε1'

ε2' < ε1'

start with:junction of two waveguides

No radiation at junctionif the modes are perfectly matched

Page 57: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Perfect Mode Matching

ε1

ε2 < ε1

ε1'

ε2' < ε1'

requires:same differential equations and boundary conditions

Match differential equations…ε2 − ε1 = ε2' − ε1' …closely related to separability

[ S. Kawakami, J. Lightwave Tech. 20, 1644 (2002) ]

Page 58: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Perfect Mode Matching

ε1

ε2 < ε1

ε1'

ε2' < ε1'

requires:same differential equations and boundary conditions

Match boundary conditions:field must be TE

(E out of plane, in 2d)

(note switch in TE/TMconvention)

Page 59: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

TE modes in 3d

for

cylindrical waveguides,

“azimuthally polarized”

TE0n modes

Page 60: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

A Perfect Cavity in 3d

R

N layers

N layers

defect

Perfect indexconfinement

(no scattering)

+1d band gap

=3d confinement

(~ VCSEL + perfect lateral confinement)

Page 61: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

A Perfectly Confined Mode

ε1', ε2' = 4, 1

ε1, ε2 = 9, 6λ/2 core

E energy density, vertical slice

Page 62: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Q limited only by finite size

G

G

G

G

G

GG

GG G

J

J

J

J

J

J

J

J

J

J

10

100

1000

10000

100000

1000000

1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6

Q

cladding R / core radius

N = 5

N = 10

~ fixed mode vol.V = (0.4λλλλ)3

Page 63: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Q-tipsThree independent mechanisms for high Q:

Delocalization : trade off modal size for Q

Multipole Cancellation : force higher-order far-field pattern

Qr grows monotonically towards band edge

Qr peaks inside gap

New nodal planes appear in far-field pattern at peak

Mode Matching : allows arbitrary Q, finite V

Requires special symmetry & materials

Page 64: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

Forget these devices…

I just want a mirror.

ok

Page 65: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

ω

k||

modesin crystal

TE

TM

Projected Bands of a 1d Crystal(a.k.a. a Bragg mirror)

incident light

k|| conserved

(normalincidence)

1d b

and

gap

light

line

of a

ir ω =

ck||

Page 66: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

ω

k||

modesin crystal

TE

TM

Omnidirectional Reflection[ J. N. Winn et al, Opt. Lett. 23, 1573 (1998) ]

light

line

of a

ir ω =

ck||

in these ω ranges,there is

no overlapbetween modesof air & crystal

all incident light(any angle, polarization)

is reflectedfrom flat surface

needs: sufficient index contrast & nhi > nlo > 1

Page 67: Complete Band Gaps - Massachusetts Institute of Technologyjdj.mit.edu/photons/tutorial/L4-slabs.pdfrod-slab projected band diagram j j j j j j j j j jjjjjj jjj jjjjj j j j j j j j

[ Y. Fink et al, Science 282, 1679 (1998) ]

Omnidirectional Mirrors in Practice

1

1.2

1.4

1.6

1.8

2

2.2

2.4

2.6

2.8

3

1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2

Inde

x ra

tio,n

2 /n

1

Smaller index, n1

0%

10%

20%

30%

40%

50%

∆λ/λmid

6 9 12 150

5 0

10 0

0

5 0

10 0

0

5 0

10 0

0

5 0

10 0

0

5 0

10 0

normal

450 s

450 p

800 p

800 s

Refle

cta

nc

e (

%)

Wavelength (microns)

Te / polystyrene

Ref

lect

ance

(%

)

contours of omnidirectional gap size