need suggestions in creating EFS

----- Forwarded Message ----
From: andika asyuda <[email protected]>
To: [email protected]
Sent: Friday, March 28, 2008 6:57:33 AM
Subject: need suggestions in creating EFS



----- Forwarded Message ----
From: andika asyuda <[email protected]>
To: [email protected]
Sent: Thursday, March 27, 2008 11:41:26 PM
Subject: need suggestions in creating EFS



----- Forwarded Message ----
From: andika asyuda <[email protected]>
To: [email protected]
Sent: Thursday, March 27, 2008 10:46:42 PM
Subject: need suggestions in creating EFS

Hi every MPb user.

My name is Andika. I get an assignment to construct EFS in photonic crystal case. Usually, I use autocad for manual plotting for this task, I found my previous method is time consuming. . But now, I need to plot EFS by matlab, or something else. I use following ctl file

(set! num-bands 4)

(set! geometry-lattice (make lattice (size 1 1 no-size)
                         (basis1 (/ (sqrt 3) 2) 0.5)
                         (basis2 (/ (sqrt 3) 2) -0.5)))

(set! geometry (list (make block
                       (center 0 0 0) (size 0.7 0.7 infinity) (e1 0.5 0.5 0) (e2 0.866 -0.866 0)
                       (material (make dielectric (epsilon 12.96))))))

(set! k-points (list 
    (vector3 0 0 0)
    (vector3 -0.2887 0.2887 0)
    (vector3 -0.236209 0.314945 0)
        (vector3 0 0 0)
    (vector3 -0.183718 0.341191 0)
        (vector3 0 0 0)
    (vector3 -0.131227 0.367436 0)
        (vector3 0 0 0)
    (vector3 -0.0787364 0.393682 0)
        (vector3 0 0 0)
    (vector3 -0.0262455 0.419927 0)
         (vector3 0 0 0)
    (vector3 0.0262455 0.446173 0)
         (vector3 0 0 0)
    (vector3 0.0787364 0.472418 0)
         (vector3 0 0 0)
    (vector3 0.131227 0.498664 0)
         (vector3 0 0 0)
    (vector3 0.183718 0.524909 0)
         (vector3 0 0 0)
    (vector3 0.236209 0.551155 0)
         (vector3 0 0 0)
    (vector3 0.2887 0.5774 0)
         (vector3 0 0 0)
                     
))
            
            
(set! k-points (interpolate 100 k-points))

(set! resolution 32)

(run-tm (output-at-kpoint (vector3 (/ 3) (/ 3) 0)
              fix-efield-phase output-efield-z))


then, I use contour command in matlab. But, I do not get closed loop curve.  . Please give me some suggestions or guidance in crrating EFS. It wll be a great help if I get a detail explanation. I will use this EFS to do wave vector formalisn analysis. I found my previous method is time consuming. in ploting and converting k-point from reciprocal lattice to x-y basis



      Looking for last minute shopping deals?  
Find them fast with Yahoo! Search.




      Be a better friend, newshound, and know-it-all with Yahoo! Mobile.  Try it now.




      Never miss a thing.   Make Yahoo your homepage.





      ____________________________________________________________________________________
Looking for last minute shopping deals?  
Find them fast with Yahoo! Search.  http://tools.search.yahoo.com/newsearch/category.php?category=shopping

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