Problem in 1-D waveguide
Ashina Garg <[email protected]> Tue, 20 Oct 2015 13:09:12 +0530
| Newsgroups | gmane.comp.science.photonic-bands |
|---|---|
| Message-ID | <CAJ8oNpniVCm3czR4Vg8CoBWJ5ypbzv9+AF53KBKGz2Naqf1K3Q@mail.gmail.com> |
--001a114e5252acc4dd0522845bb8
Content-Type: multipart/alternative; boundary=001a114e5252acc4d20522845bb6
--001a114e5252acc4d20522845bb6
Content-Type: text/plain; charset=UTF-8
I want to create a line defect 1-D PhC waveguide. Below is my ctl file and
the image of epsilon file. I have a default material of 4.41 dielecric
constant and i have inserted blocks of air. To temove a single block, I
have inserted a block of dielectric in the centre of slightly less radius.
I want a single block of dielectric only in the centre which is of same
radius as all other blocks. But whenever the run the code, I'm getting a
dielectric in the centre of each block. Please help me find where am I
going wrong?
(define-param x-dim 1)
(define-param y-dim 20)
(define-param z-dim 4)
(set! geometry-lattice (make lattice(size 1 1 no-size)))
(set! default-material (make dielectric (epsilon
4.41)))
(set! geometry(list(make block (material (make dielectric (epsilon 1)))
(center 0 0 0) (size 0.3 1 infinity))))
(set! geometry (append geometry
(list (make block (material (make dielectric (epsilon 4.41)))
(center 0 0 0) (size 0.1 1 infinity)))))
(set! resolution 16)
(set! k-points (list (vector3 0.5 0.5 0)))
(set! num-bands 50)
(run)
--001a114e5252acc4d20522845bb6
Content-Type: text/html; charset=UTF-8
Content-Transfer-Encoding: quoted-printable
<div dir=3D"ltr">I want to create a line defect 1-D PhC waveguide. Below is=
my ctl file and the image of epsilon file. I have a default material of 4.=
41 dielecric constant and i have inserted blocks of air. To temove a single=
block, I have inserted a block of dielectric in the centre of slightly les=
s radius. I want a single block of dielectric only in the centre which is o=
f same radius as all other blocks. But whenever the run the code, I'm g=
etting a dielectric in the centre of each block. Please help me find where =
am I going wrong?<br><br>(define-param x-dim 1)<br>(define-param y-dim 20)<=
br>(define-param z-dim 4)<br><br>(set! geometry-lattice (make lattice(size =
1 1 no-size)))<br>=C2=A0(set! default-material (make dielectric (epsilon 4.=
41)))=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=
=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0 <br>(set! g=
eometry(list(make block (material (make dielectric (epsilon 1)))<br>=C2=A0=
=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0 (c=
enter 0 0 0) (size 0.3 1 infinity))))<br>(set! geometry (append geometry<br=
>(list (make block (material (make dielectric (epsilon 4.41)))<br>=C2=A0=C2=
=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0=C2=A0 (cent=
er 0 0 0) (size 0.1 1 infinity)))))<br>=C2=A0<br>(set! resolution 16)<br>(s=
et! k-points (list (vector3 0.5 0.5 0)))<br>(set! num-bands 50)<br><br><br>=
(run)<br><br></div>
--001a114e5252acc4d20522845bb6--
--001a114e5252acc4dd0522845bb8
Content-Type: image/png; name="epsilon.png"
Content-Disposition: attachment; filename="epsilon.png"
Content-Transfer-Encoding: base64
X-Attachment-Id: f_ifz27u6i0
iVBORw0KGgoAAAANSUhEUgAAAUAAAAFACAIAAABC8jL9AAADxElEQVR4nO3TsQmAABDAwMfJf3Ot
rCxFJHA3QLrMzOzuedvdeedZ09fX/6h/vMwBPzIwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjC
DAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEh
zMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkY
wgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iB
IczAEGZgCDMwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZ
GMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OY
gSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCE
GRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMAND
mIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMw
hBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjAD
Q5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgz
MIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYw
A0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEhzMAQZmAI
MzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZgCDMwhBkYwgwMYQaG
MANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEGhjADQ5iBIczAEGZg
CDMwhBkYwgwMYQaGMANDmIEhzMAQZmAIMzCEGRjCDAxhBoYwA0OYgSHMwBBmYAgzMIQZGMIMDGEG
hjADQ5iBIczAEGZgCDMwhBkYwgwMYQaGMANDmIEh7AJubPHL8uJP9QAAAABJRU5ErkJggg==
--001a114e5252acc4dd0522845bb8
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline
X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX18KbXBiLWRpc2N1
c3MgbWFpbGluZyBsaXN0Cm1wYi1kaXNjdXNzQGFiLWluaXRpby5taXQuZWR1Cmh0dHA6Ly9hYi1p
bml0aW8ubWl0LmVkdS9jZ2ktYmluL21haWxtYW4vbGlzdGluZm8vbXBiLWRpc2N1c3M=
--001a114e5252acc4dd0522845bb8--