Resonant frequencies in PC slabs!

  Dear Steven and MPB's users, 
    I would like to compute the defect modes (Resonant frequencies) in photonic crystal slab. (According to the document of MPB User tutorial - Finding a point defect state). I do not understand what my .ctl file is wrong! 
   I have a PC slab structure - Triangular lattice (Free standing PC slab with radius of holes: r=0.3a; slab thickness: h =0.4a; dielectric constant: eps =11.56). With this structure, I have calculated the PBG: ~0.28 - 0.36.
   Defect point is formed by removing the center hole, and add another hole in the center of the cell with a dielectric constant of 11.56;
   After making the structure, I have some MPB commands as below:    
   
  (set! resolution 16)
  (set! k-points (list (vector3 (/ -3) (/ 3) 0) ))
  (set! num-bands 50)
  (run-zeven)
  (output-hfield-z 25)
  (get-hfield 25) ; compute the H field for band 25
  (compute-field-energy) ; compute the energy density from H
  (output-efield-z 25)
  (print
   "energy in cylinder: "
   (compute-energy-in-objects (make cylinder (center 0 0 0) (radius r) (height h)
     (material (make dielectric (epsilon eps)))))
      "\n")
  (print "********** 9x9 point defect, targeted solver **********\n")
  (set! num-bands 2) ; only need to compute a single band, now!
  (set! target-freq (/ (+ 0.28 0.36) 2))
  (set! tolerance 1e-8)
  (run-zeven)
  
I really would like to get your suggestions about this trouble!
   
  I thank you so much! and I hope you reply to me soon.
  Sincerely,
  Ngo Quang Minh


 
---------------------------------
Low, Low, Low Rates! Check out Yahoo! Messenger's cheap  PC-to-Phone call rates.

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