Re: What can I define the assymmetric photonic crystal slabs? (Correctly!)

matt <[email protected]>
Newsgroups gmane.comp.science.photonic-bands
Message-ID <[email protected]>

It's hard to tell from your email what isn't working.  Also, you didn't 
give your complete source, so it's hard to see if there's some obvious 
mistake (what do you specify for geometry-lattice, for example?)

I was able to get what you described with the following:

(set! num-bands 8)
(define-param eps1 2)
(define-param eps2 3)
(define-param eps0 4)
(define-param r .3)

(set! geometry
       (list
        (make block (material (make dielectric (epsilon eps1))) (center 0 0 -0.2) (size 1 1 0.2))
        (make block (material (make dielectric (epsilon eps0))) (center 0 0 -0.05) (size 1 1 0.1))
        (make block (material (make dielectric (epsilon eps2))) (center 0 0 0.15) (size 1 1 0.3))
        (make cylinder (material air) (center 0 0 0) (radius r) (height 0.6))
        ))
(set! geometry-lattice (make lattice (size 1 1 0.6)))
(set! resolution 64)

(set! k-points (list (vector3 0 0 0)     ; Gamma
                      (vector3 0.5 0 0)   ; X
                      (vector3 0.5 0.5 0) ; M
                      (vector3 0 0 0)))   ; Gamma

(set! k-points (interpolate 1 k-points))
(run-te)
--------------------
% h5topng -x 0 -S 3 epsilon.h5
% display epsilon.png


Kind Regards,
Matt



On Tue, 13 Feb 2007, Nils Jensen wrote:

> Dear Steven and MPB's user,
   I have a problem about MPB would like to discuss with you:

I would like to calculate the band structure of the PC slab: (have a three layer). The first layer is the clad 1 (with refractive index n1, and thickness 0.2a); The second layer is the core (with refractive index n0, and thickness 0.1a); The third layer is the clad 2 (with refractive index n2, and thickness 0.3a). I define my structure as below:   (n1 < n2 < n0)

(set! geometry
(list
    (make block (material (make dielectric (epsilon eps1))) (center 0 0 -0.2) (size 1 1 0.2))
    (make block (material (make dielectric (epsilon eps0))) (center 0 0 -0.05) (size 1 1 0.1))
    (make block (material (make dielectric (epsilon eps2))) (center 0 0 0.15) (size 1 1 0.3))
    (make cylinder (material air) (center 0 0 0) (radius r) (height 0.6)) 
))

It seems MPB does not support the asymmetric structure, because I check my structure by using this command: h5topng -x (-y) ....The cross structure that I get, is not correct!
    Please let me know why it is?

Thanks,
Nils Jensen



____________________________________________________________________________________
No need to miss a message. Get email on-the-go 
with Yahoo! Mail for Mobile. Get started.
http://mobile.yahoo.com/mail
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.