Hi,all.
When I used MEEP to find modes, following the tutorial, the output was a
little different from the tutorial. In the tutorial, the 6th column was
error and both real and imaginary were given. But in the outputs I got , as
below, the 5th column looks more like error. Did I misunderstood?
harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error
harminv0:, 0.200068338343141, -6.44246544433111e-7, 155273.117156714,
0.181773665809152, -0.00285709430770058+0.181751210707974i,
1.38429090645002e-7
harminv0:, 0.396517414512765, 3.12197776468352e-8, -6350420.22076926,
6.01863473151113, 5.94258839734845+0.953733280876562i, 2.07446847341345e-9
harminv0:, 0.451456811625371, -4.53066538043226e-8, 4982235.20959188,
6.21148573078014, -2.60806998247918-5.63742192408696i, 3.6554550330908e-9
harminv0:, 0.541228570195831, -1.24345057501898e-9, 217631718.167637,
2.72191028281296, 1.98859335997502+1.85857258032777i, 5.63431336862739e-10
harminv0:, 0.647353541562013, -1.5294418123807e-8, 21163065.3850883,
3.24407420463741, -3.14294479500851+0.803688659073113i, 1.09432171289535e-9
harminv0:, 0.701592821418967, 7.18355839346233e-9, -48833237.1640132,
1.07549227402891, 0.769211650564473+0.751662868666374i, 7.19555520626556e-10
my ctl file:
(define-param a 1)
(define-param eps 13.022)
(define-param thick 0.6)
(define-param r (* 0.25 a))
;The cell dimensions
(define-param sz 5);
(define-param dpml 1);
(set! geometry-lattice (make lattice (size (/ a 1) (* a (sqrt 3)) no-size)))
;(set! geometry-lattice (make lattice (size a a no-size)))
(set! geometry
(list ;(make block (center 0 0) (size infinity infinity infinity)
(material (make dielectric (epsilon eps))))
; (make cylinder (center 0 0) (radius r) (height infinity) (material (make
dielectric (epsilon eps)) ))
; (make cylinder (center 0 0) (radius r) (height infinity) (material (make
dielectric (epsilon eps)) ))
(make cylinder (center (/ a 4) (* (/ (sqrt 3) 4) a ))
(radius r) (height infinity) (material (make dielectric (epsilon
eps)) ))
(make cylinder (center (/ a -4) (* (/ (sqrt 3) -4) a ))
(radius r) (height infinity) (material (make dielectric (epsilon
eps)) ))
;(make cylinder (center 0 0) (radius r) (height infinity) (material (make
dielectric (epsilon eps)) ))
; (make cylinder (center (/ a 4) (* (/ (sqrt 3) 4) a)) (radius r) (height
infinity) (material air))
; (make cylinder (center (/ a -4) (* (/ (sqrt 3) -4) a)) (radius r) (height
infinity) (material air))
; (make cylinder (center (/ a 2) (* (/ (sqrt 3) 2) a)) (radius r) (height
infinity) (material (make dielectric (epsilon eps))))
; (make cylinder (center (/ a 2) (* (/ (sqrt 3) -2) a)) (radius r) (height
infinity) (material (make dielectric (epsilon eps))))
; (make cylinder (center (/ a -2) (* (/ (sqrt 3) 2) a)) (radius r) (height
infinity) (material (make dielectric (epsilon eps))))
; (make cylinder (center (/ a -2) (* (/ (sqrt 3) -2) a)) (radius r) (height
infinity) (material (make dielectric (epsilon eps))))
))
;(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(set-param! resolution 30)
(define-param fcen 0.4); pulse center frequency
(define-param df 0.8); pulse freq. width:large df=shot impulse
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hz) (center 0.12 0 0))))
;(set! symmetries (list (make mirror-sym (direction Z) (phase -1))))
(define-param kx false); if true, du run at specified kx and get fields
(define-param k-interp 9)
(if kx
(begin
(set! k-point (vector3 0.25))
(run-sources+ 300 (at-beginning output-epsilon)
(after-sources (harminv Hz (vector3 0.1234 0) fcen df)))
(run-until (/ 1 fcen) (at-every (/ 1 fcen 20) output-hfield-z)))
(run-k-points 1200 (interpolate k-interp (list (vector3 0) (vector3 0.5)
(vector3 0.5 (/ (sqrt 3) 6) 0) (vector3 0)))))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
