Hi Zhanghua, I have the same problem. In my situation, I find the result from meep-0.10.1 is right, while the one given by meep-0.20.3 is wrong, for the result from meep-0.10.1 almost overlaps with MPB's. sincerely, Li The control file and results (the corresponding instantaneous field distribution are attached) are as follows, ***control file*** ---------------------------------------------
(define-param sx 1)
(define-param sy 12)
(define-param sz 3)
(set! default-material (make dielectric (index 3.45)))
(set! geometry-lattice (make lattice (size sx sy sz)))
(define-param r 0.25)
(define thick 0.6)
(define halfthick (/ thick 2))
(define SiO2-height (- (/ sz 2) halfthick))
(define half-SiO2 (/ SiO2-height 2))
(define obj-list1 (make cylinder (center 0.5 0.571577 0) (radius r) (height
thick) (material air)))
(define obj-list2 (make cylinder (center 0.5 -0.571577 0) (radius r) (height
thick) (material air)))
(define obj-list3 (make cylinder (center 0 1.4376 0) (radius r) (height
thick) (material air)))
(define obj-list4 (make cylinder (center 0 -1.4376 0) (radius r) (height
thick) (material air)))
(define obj-list5 (make cylinder (center -0.5 0.571577 0) (radius r) (height
thick) (material air)))
(define obj-list6 (make cylinder (center -0.5 -0.571577 0) (radius r) (height
thick) (material air)))
(set! geometry
(append
(geometric-object-duplicates (vector3 0 (sqrt 3)) 0 3 obj-list1)
(geometric-object-duplicates (vector3 0 (-(sqrt 3))) 0 3 obj-list2)
(geometric-object-duplicates (vector3 0 (sqrt 3)) 0 2 obj-list3)
(geometric-object-duplicates (vector3 0 (-(sqrt 3))) 0 2 obj-list4)
(geometric-object-duplicates (vector3 0 (sqrt 3)) 0 3 obj-list5)
(geometric-object-duplicates (vector3 0 (-(sqrt 3))) 0 3 obj-list6)
(list (make block (material (make dielectric (index 1.45))) (size sx
sy SiO2-height) (center 0 0 (+ halfthick half-SiO2))))
(list (make block (material air) (size sx sy SiO2-height) (center 0 0
(- (- halfthick) half-SiO2))))
(list (make block (material air) (size sx 2 sz) (center 0 6 0)))
(list (make block (material air) (size sx 2 sz) (center 0 -6 0)))
))
(define dpml 0.5)
(set! pml-layers (list (make pml (direction Y) (thickness dpml))
(make pml (direction Z) (thickness dpml))
))
(set-param! resolution 20)
(define fcen 0.245)
(define df 0.2)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ex) (amplitude 1e5) (size 0 0)
(center 0.1 0.2 0.1))
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ex) (amplitude -1e5) (size 0 0)
(center 0.1 -0.2 0.1))
))
(set! symmetries (list (make mirror-sym (direction Y) (phase -1))))
(use-output-directory)
(define-param kx false)
(define k-interp 9)
(if kx
(begin
(set! k-point (vector3 kx 0 0))
(run-sources+ 300
(at-beginning output-epsilon)
(after-sources (harminv Ex (vector3 0.1 0.2 0.1) fcen df)))
(run-until (/ 1 fcen) (at-every (/ 1 fcen 20) output-hfield-z))
)
(run-k-points 300 (interpolate k-interp (list (vector3 0.4) (vector3
0.5)))))
---------------------------------------------------------------
****result from Meep 0.10.1*** (for kx=0.455)
harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error
harminv0:, 0.246201171406986, 9.928119860464e-7, -123991.84078519,
61322.7060466497, -36910.7728140296+48970.084001918i, 1.08583514711351e-7
harminv0:, 0.268271329143659, 8.84533953766205e-4, -151.64558013935,
30.9683519707873, -12.4137304244833+28.3714314185031i, 1.7967056297694e-5
harminv0:, 0.296206106098515, -0.002446469954863, 60.5374501962978,
8972.94922627343, 4310.28322942749+7869.89684172524i, 1.71456122363114e-6
harminv0:, 0.342461339705382, -4.27845597692274e-4, 400.216037692756,
128.310171430319, -12.45306209249-127.704429590357i, 5.04174366078636e-5
---------------------------------------------------------------
****result from Meep 0.20.3*** (for kx=0.455)
harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error
harminv0:, 0.200027330247863, -5.79995633162289e-4, 172.438651957827,
852.986945080003, -820.762756084237+232.239589006533i, 3.99152820955429e-4
harminv0:, 0.244261717858064, 1.86389607229993e-5, -6552.45003968114,
57580.4157695156, -44614.4748421462+36401.2762791333i, 3.69518641021594e-6
harminv0:, 0.288999788766653, 3.91049360446114e-4, -369.518298709093,
1203.42587942411, 291.954794948466+1167.47430163335i, 1.48316971634606e-4
harminv0:, 0.321636684572457, -7.35563680476834e-5, 2186.32793535941,
1295.81755875124, -1211.43911497242-459.954798086586i, 1.45758995846933e-4<<attachment: 0.10.1.png>>
<<attachment: 0.20.3.png>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
