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

_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss

Reply via email to