# 18.369 Mathematical Methods in Nanophotonics ```MIT OpenCourseWare
http://ocw.mit.edu
18.369 Mathematical Methods in Nanophotonics
Spring 2008
18.325 Problem Set 3
Problem 1: Perturbation theory
thickness dl and that the ~1 - AE region has thickness a-dl, repeated with
period a, and compute Aw as a function of dl.
(ii) At what dl values does this gap at k =
0 disappear (Aw = O)? What is Aw for
the quarter-wave stack thicknesses?
(a) In class, we derived the 1st-order correction
in the eigenvalue for an ordinary Hermitian
eigenproblem 0 lu) = u lu) for a small per2: Band gaps in MPB
turbation AO. Now, do the same thing for
a generalized eigen~roblemA I u ) = uB Iu)Consider the i d periodic structure consisting of
two alternating layers: ~1 = 12 and ~2 = 1, with
That is, assume we have the solution
thicknesses dl and d2 = a - dl, respectively. To
ace) lU(o)) = , ( o ) ~ ( o ) lU(o))to
help you with this, I've created a sample input
perturbed system and find the firstfile bandgapid.ct1 that is posted on the course
order correction u(') when we change
web page.
both A and B by small amounts AA
and AB. You may assume that u(O) is (a) Using MPB, compute and plot the fractional
TM gap size (of the f i r s t gap, i.e lowest w)
non-degenerate, for simplicity.
vs.
dl for dl ranging from 0 to a. What
(ii) Now, apply this solution to the gendl
gives
the largest gap? Compare t o the
eralized eigenproblem V x V x E =
quarter-wave thicknesses.
\$EE
for a small change A&amp;,and show
that the first-order correction Aw is the (b) Given the optimal parameters above, what
same as the one derived in class using
would be the physical thicknesses in order
the H eigenproblem.
for the mid-gap vacuum wavelength to be
X = 2 ~ c l w= 1.55pm? (This is the wave(b) Recall the problem of the modes in an L x L
length used for most optical telecommunimetal box that we solved in class for the Hz
cations.)
(TE) polarization, and which you solved in
problem set 1 for the Ez (TM) polarization. (c) Plot the Id TM band diagram for this structure, with dl given by the quarter wave
Originally, this box was filled with air ( E =
thickness,
showing the first five gaps. Also
1). Now, suppose that we increase E by some
compute
it
for dl = 0.12345 (which I just
small constant A&amp; in the lower-left \$ x \$
chose
randomly),
and superimpose the two
corner of the box. What is the first-order
plots
(plot
the
quarter-wave
bands as solid
Aw for the first four (lowest w) modes of
lines
and
the
other
bands
as
dashed).
What
the TM polarization? What happens to the
special
features
does
the
quarter-wave
band
degenerate modes?
diagram have?
(c) In class, we calculated the band gap Aw
that appeared in a uniform (Id) material Problem 3: Space group of k
E(X) = ~1 when half of the unit cell a was
changed to ~1 - AE for some small AE, by Consider the structure of problem 2. What are
using perturbation theory. In particular, we the symmetry operations of its space group?
derived the gap between the first two bands How are these reduced (if at all) if we consider
that appears at k = n l a , the edge of the Bk-that is, what is the space group as a function of k? In class, when we derived the band
Brillouin zone.
a be of
gap, we took the modes at k = ~ / to
(i) Now, you should compute the gap that even/odd (cos/sin) form; why was this justified?
Give an example of a Id periodic structure in
appears between band 2 and band 3 at
k = 0, the center of the Brillouin zone, which the modes at 5 = n/a are n o t even or
to first order. However, do it more gen- odd. Is the band diagram of this structure still
erally: assume that the ~1 region has symmetric (i.e. w(- k) = w( k ) ) ?
,
Problem 4: Defect modes in MPB
In MPB, you will create a (TM polarized) defect
mode by increasing the dielectric constant of a
single layer by Ae, pulling a state down into the
gap. The periodic structure will be the same as
the one from problem 2, with the quarter-wave
thickness dl = 1/(1
this, I've created a sample input file defectld.ct1
that is posted on the course web page.
+ m).
(a) When there is no defect (A&amp;) plot out the
band diagram w(k) for the N = 5 supercell,
and show that it corresponds to the band
diagram of problem 2 &quot;folded&quot; as expected.
(b) Create a defect mode (a mode that lies in
the band gap of the periodic structure) by
increasing the e of a single el layer by Ae =
1, and plot the Ez field pattern. Do the
same thing by increasing a single ~2 layer.
Which mode is even/odd around the mirror
plane of the defect? Why?
your conjecture1 with numerical evidence.
(Don't forget that, like in problem set 2, as
the mode becomes less localized you must
increase the computational cell size. Another thing to be careful about is how accurately you know the location of the gap
edge-compute the gap edge using the same
resolution as you are using to compute the
defect, and you may need to increase both
resolutions .)
(f) The mode must decay exponentially far
from the defect (multiplied by an eiEx sign
oscillation and the periodic Bloch envelope,
of course). From the Ez field computed
by MPB, extract this asympotic exponential
decay rate (i.e. K if the field decays e-&quot;&quot;)
and plot this rate as a function of w, for
the first defect mode, as you increase e2 as
above (vary ~2 SO that w goes from the top
of the gap to the bottom).
N
(c) Gradually increase the e of a single ~2 layer,
and plot the defect was a function of AE
as the frequency sweeps across the gap. At
what A&amp; do you get two defect modes in
the gap? Plot the Ez of the second defect
mode. (Be careful to increase the size of the
supercell for modes near the edge of the gap,
which are only weakly localized.)
(d) Via first-order perturbation theory, derive
an exact expression for the rate of change
4 f ~ o)f the defect-mode frequency w of a
single defect layer ~2
AE layer, as AE
is changed. Your expression should be in
terms of the eigenfield E of the localized defect state at the current Ae. Verify your formula numerically by showing, at a couple of
different values of Ae, that it correctly predicts the slope of your w vs. Ae curve above.
(Note that you can export E to Matlab and
compute the necessary integrals there or,
for the Scheme lovers among you, you can
compute the integral directly in MPB using
the compute-energy-in-object s or computefield-integral function .)
+
(e) Is there a minimum Ae (of an e2 layer) t o
create a defect mode, or is a defect mode lo11,
we never merely guess^
calized even for infinitesimal Ae? Support jecture, or perhaps we postulate an ansatz.
we con-
```
{"ft_lang_label":"__label__en","ft_lang_prob":0.9190389,"math_prob":0.9718331,"size":6627,"snap":"2022-40-2023-06","text_gpt3_token_len":1667,"char_repetition_ratio":0.12547185,"word_repetition_ratio":0.012789768,"special_character_ratio":0.2507922,"punctuation_ratio":0.099928625,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.9966573,"pos_list":[0,1,2],"im_url_duplicate_count":[null,1,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2022-09-25T01:44:03Z\",\"WARC-Record-ID\":\"<urn:uuid:e10476a7-78fb-4263-84c9-9ea83dd5c989>\",\"Content-Length\":\"45089\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:60a1c5d1-e1d6-493c-8a95-a0e391751b18>\",\"WARC-Concurrent-To\":\"<urn:uuid:8a5cd319-2dd0-4db0-b938-9d831c311a5a>\",\"WARC-IP-Address\":\"104.21.57.229\",\"WARC-Target-URI\":\"https://studylib.net/doc/13575619/18.369-mathematical-methods-in-nanophotonics\",\"WARC-Payload-Digest\":\"sha1:WQC2SPKH3UCFK2ZAZVDSYRAWWTLGXIUN\",\"WARC-Block-Digest\":\"sha1:KDLYF2OUMKPO2SFIEFTE5ZL4WVNZQMQH\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2022/CC-MAIN-2022-40/CC-MAIN-2022-40_segments_1664030334332.96_warc_CC-MAIN-20220925004536-20220925034536-00762.warc.gz\"}"}