equivalent symmetry points, different results
Janos Perczel <[email protected]> Wed, 15 Feb 2017 02:08:12 +0000
| Newsgroups | gmane.comp.science.photonic-bands |
|---|---|
| Message-ID | <76223849952F344E8A155FF5852E880CAB8B68AE@OC11EXPO29.exchange.mit.edu> |
--===============3833275134930253255==
Content-Language: en-GB
Content-Type: multipart/alternative;
boundary="_000_76223849952F344E8A155FF5852E880CAB8B68AEOC11EXPO29excha_"
--_000_76223849952F344E8A155FF5852E880CAB8B68AEOC11EXPO29excha_
Content-Type: text/plain; charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
Dear Professor Steven and MPB users,
Using the code below, I found that two theoretically equivalent X points in=
the BZ of a square lattice give very different results for the output fiel=
ds. I have a lattice of air holes in a PC slab and I am outputting the fiel=
d of the second band at two equivalent X points.
Has anyone else encountered this problem before?
Thank you for considering my question!
Janos
Department of Physics
MIT
(define-param h 0.316) ; the thickness of the slab
(define-param eps 10.56) ; the dielectric constant of the slab
(define-param r 0.4) ; the radius of the holes
(define-param supercell-h 4) ; height of the supercell
(set! geometry-lattice (make lattice (size 1 1 supercell-h)))
(set! geometry
(list (make block (material (make dielectric (epsilon eps)))
(center 0) (size infinity infinity h))
(make cylinder (material air)
(center 0) (radius r) (height supercell-h))))
(set! k-points (list (vector3 0 0.5 0) (vector3 0.5 0 0) ))
(set-param! resolution (vector3 32 32 16))
(set-param! num-bands 2)
(define (output-some-efields b)
(if (=3D b 2)
(output-efield b)))
(begin-time ;
(run-zeven fix-efield-phase output-some-efields))
(display-eigensolver-stats)
--_000_76223849952F344E8A155FF5852E880CAB8B68AEOC11EXPO29excha_
Content-Type: text/html; charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<html dir=3D"ltr">
<head>
<meta http-equiv=3D"Content-Type" content=3D"text/html; charset=3Diso-8859-=
1">
<style type=3D"text/css" id=3D"owaParaStyle"></style>
</head>
<body fpstyle=3D"1" ocsi=3D"0">
<div style=3D"direction: ltr;font-family: Tahoma;color: #000000;font-size: =
10pt;">Dear Professor Steven and MPB users,
<div><br>
</div>
<div>Using the code below, I found that two theoretically equivalent X poin=
ts in the BZ of a square lattice give very different results for the output=
fields. I have a lattice of air holes in a PC slab and I am outputting the=
field of the second band at two
equivalent X points.</div>
<div><br>
</div>
<div>Has anyone else encountered this problem before?</div>
<div><br>
</div>
<div>Thank you for considering my question!</div>
<div><br>
</div>
<div>Janos</div>
<div><br>
</div>
<div>Department of Physics</div>
<div>MIT</div>
<div><br>
</div>
<div><br>
</div>
<div>
<div>(define-param h 0.316) ; the thickness of the slab</div>
<div>(define-param eps 10.56) ; the dielectric constant of the slab</div>
<div>(define-param r 0.4) ; the radius of the holes</div>
<div>(define-param supercell-h 4) ; height of the supercell</div>
<div><br>
</div>
<div>(set! geometry-lattice (make lattice (size 1 1 supercell-h)))</div>
<div><br>
</div>
<div>(set! geometry</div>
<div> (list (make block (material (make dielectric (eps=
ilon eps)))</div>
<div><span class=3D"Apple-tab-span" style=3D"white-space:pre"></span> =
(center 0) (size infinity infinity h))</div>
<div><span class=3D"Apple-tab-span" style=3D"white-space:pre"></span> =
(make cylinder (material air)</div>
<div><span class=3D"Apple-tab-span" style=3D"white-space:pre"></span> =
(center 0) (radius r) (height supercell-h))))</div>
<div><br>
</div>
<div>(set! k-points (list (vector3 0 0.5 0) (vector3 0.5 0 0) ))</div>
<div><br>
</div>
<div>(set-param! resolution (vector3 32 32 16))</div>
<div>(set-param! num-bands 2)</div>
<div><br>
</div>
<div>(define (output-some-efields b)</div>
<div>(if (=3D b 2)</div>
<div>(output-efield b)))</div>
<div><br>
</div>
<div>(begin-time ;</div>
<div>(run-zeven fix-efield-phase output-some-efields))</div>
<div>(display-eigensolver-stats)</div>
</div>
<div> </div>
</div>
</body>
</html>
--_000_76223849952F344E8A155FF5852E880CAB8B68AEOC11EXPO29excha_--
--===============3833275134930253255==
Content-Type: text/plain; charset="utf-8"
MIME-Version: 1.0
Content-Transfer-Encoding: base64
Content-Disposition: inline
X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX18KbXBiLWRpc2N1
c3MgbWFpbGluZyBsaXN0Cm1wYi1kaXNjdXNzQGFiLWluaXRpby5taXQuZWR1Cmh0dHA6Ly9hYi1p
bml0aW8ubWl0LmVkdS9jZ2ktYmluL21haWxtYW4vbGlzdGluZm8vbXBiLWRpc2N1c3M=
--===============3833275134930253255==--