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
