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
