Since i did not subscribe this mail-list before and my email address was not in the list, the server rejected this e-mail. But Steven maybe saw the mail from the server log, and replied it. I am so grateful to him, and i think the original mail may be interesting to someone.
Zhichao
---------- Forwarded message ----------
Date: 2006-11-3 下午4:55
Subject: calculation with symmetry
Hi, all
I calculate a two-dimensional structure which is periodic along the
y-direction and mirror-symmetric about z=0 plane. The following ctl file is
used to obtain the resonate frequency of the Bloch mode with ky= 0.5*2pi/a.
If I did not add the symmetry, the results are correct. It is very
surprising that Meep does not give any resonate frequency when I add the
symmetry restriction.
Something wrong in my ctl file? Or the translational symmetry and the
mirror symmetry can not be applied simultaneously in Meep?
Best regards
Zhichao Ruan
(define-param a 1)
(define-param wy (/ 1 16))
(define-param wx 1.2)
(define-param thick (/ 25 4))
(define-param zhigh 15)
; The cell dimensions
(define-param dpml 1) ; PML thickness (z direction only!)
(set-param! resolution 32)
(set! geometry-lattice (make lattice (size no-size a zhigh)))
(set! geometry
(list (make block (center 0 0 0) (size infinity 1.2 thick)
(material (make perfect-metal)))
(make block (center 0 0 0) (size infinity wy (+ thick 0.1))
(material air))))
(set! pml-layers (list (make pml (direction Z) (thickness dpml))))
(define-param fcen 0.25) ; pulse center frequency
(define-param df 0.49) ; pulse freq. width: large df = short impulse
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hx) (center 0 (/ 1 16) 0))
))
;;;;;;;;;;
(set! symmetries (list (make mirror-sym (direction Z))))
;;;;;;;;;;
(run-k-points 300 (list (vector3 0 0.5 0)))
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
