I'd say your results are pretty close! See: harminv0: *0.26204* -2.85627E-4 *458.71042* 0.14281 0.136049630579314+0.0434083155323208i 1.14604E-8
You could always try different resolutions, or larger cell sizes to minimize the influence of the PML. --Judson On Thu, Sep 6, 2012 at 3:30 AM, Onder Karakilinc <[email protected]>wrote: > Dear Meep users, > > I am trying to calculate resonant frequency and Q factor of structure > mentioned in this paper “Optimal bistable switching in nonlinear photonic > crystals, PHYSICAL REVIEW E 66, 055601, 2002”. > > w0=0.2581 and Q=557 are given, but I have obtained different results as > the following; > > > harminv0: frequency imag. freq. Q |amp| > amplitude error > harminv0: 0.19958 -5.26375E-4 189.58022 0.00137 > -0.00134089084934886+2.9314730010506e-4i 1.43784E-6 > harminv0: 0.23986 -0.00157 76.51416 0.00194 > 3.2124273191714e-4-0.00191486202485204i 2.96654E-6 > harminv0: 0.26168 0.00211 -62.01106 1.27833E-5 > 5.04840831597867e-7-1.27733253595646e-5i 3.40482E-5 > harminv0: 0.26204 -2.85627E-4 458.71042 0.14281 > 0.136049630579314+0.0434083155323208i 1.14604E-8 > harminv0: 0.3192 -3.70434E-4 430.84916 1.54142E-4 > 8.16165519586485e-5+1.30761070218176e-4i 5.32413E-6 > > > What are the key points that I have missed? > > My ctl code for harminv calculation is below > > I'll be very glad of any assitance. Best regards… > > Onder Karakilinc > > > (define-param nbase 1.5) > (define-param epsbase (* nbase nbase)) > > (define-param nrod 3.5) > (define-param epsrod (* nrod nrod)) > > (define-param a 1) > (define-param r 0.25) ; rods radius > > (define-param rsmall (/ r 3)) ; small rods in wg > (define-param rbig (/ (* 5 r) 3)) ; big rods at the center > > (define-param sx 31) > (define-param sy 25) > > (set! geometry-lattice (make lattice (size sx sy no-size))) > > (set! default-material (make dielectric (epsilon epsbase))) > > (set! geometry (list > (make cylinder (center 0 0) (height infinity) (radius r) > (material (make dielectric (epsilon epsrod)))) > )) > > (set! geometry (geometric-objects-lattice-duplicates geometry)) > > (set! geometry (append geometry > (list > (make block (center 0 0) (size sx 1 infinity) > (material (make dielectric (epsilon > epsbase) > ) > ) > ) > ) > ) > ) > > > (set! geometry (append geometry > (geometric-object-duplicates (vector3 1 0) 0 12 > (make cylinder (center 4 0) (radius rsmall) (height infinity) > (material (make dielectric (epsilon epsrod))))))) > > (set! geometry (append geometry > (geometric-object-duplicates (vector3 -1 0) 0 12 > (make cylinder (center -4 0) (radius rsmall) (height infinity) > (material (make dielectric (epsilon epsrod))))))) > > > (set! geometry (append geometry > (list > > (make cylinder (center 0 0 0) (radius rbig) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > (make cylinder (center 1 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > (make cylinder (center 2 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > (make cylinder (center 3 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > (make cylinder (center -1 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > (make cylinder (center -2 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > (make cylinder (center -3 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon epsrod))) > ) > > ))) > > > > > (define-param fcen 0.25) > (define-param df 0.2) > > (set! sources (list > > > (make source > (src (make gaussian-src (frequency fcen) (fwidth > df))) > (component Ez)(center 0.6 0) > ) > > )) > > (set! pml-layers (list (make pml (thickness 2.0)))) > (set! resolution 32) > > > (run-sources+ 300 > (at-beginning output-epsilon) > (after-sources (harminv Ez (vector3 0.6) fcen df))) > > > > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss >
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
