Dear friend, I think you have a mistake in the supercell-y parameter. Your supercell width is 3 so MPB ensures a Bloch periodicity using a supercell composed by only three holes in y direction. You should probably try with: (define-param supercell-y 10) I'm waiting for hear from you.
2010/12/27 <[email protected]> > Send mpb-discuss mailing list submissions to > [email protected] > > To subscribe or unsubscribe via the World Wide Web, visit > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss > or, via email, send a message with subject or body 'help' to > [email protected] > > You can reach the person managing the list at > [email protected] > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of mpb-discuss digest..." > > > Today's Topics: > > 1. Line-defect waveguide problem (???) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 27 Dec 2010 20:50:38 +0800 > From: ??? <[email protected]> > Subject: [MPB-discuss] Line-defect waveguide problem > To: [email protected] > Message-ID: > <[email protected]> > Content-Type: text/plain; charset="iso-8859-1" > > Dear Prof. Johnson and Meep users, > > Hi all! I am new learner of MPB software and poor in English.So please > forgive me for bad expression. > I am trying to to calculate the waveguide of the square lattice structure > that a line defect of the air holes on the epsilon.I use the > "line-defect.ctl" for reference. like below: > > (define-param supercell-y 3) ; the (odd) number of lateral supercell > periods > (set! geometry-lattice (make lattice > (basis1 1 0) > (basis2 0 1) > (size 1 supercell-y no-size))) > (define-param r 0.2) ; the rod radius in the bulk crystal > (set! geometry (list (make cylinder > (center 0 0 0) (radius r) (height infinity) > (material (make dielectric (epsilon eps)))))) > (set! geometry > (append > (geometric-objects-lattice-duplicates geometry) > (list > (make cylinder (center 0) (radius r) (height infinity) > (material air))))) > > then I got a result like that > o= holes > x=defcet > oooo.... > xxxx.... > oooo.... > oooo.... > xxxx.... > oooo.... > oooo.... > xxxx.... > oooo.... > > Here is the problem,the square lattice structure that I want to calculate > is > ONLY one line-defect in that.That means the line-defect shall not be > periodic.Like below: > > oooo.... > oooo.... > oooo.... > oooo.... > xxxx.... > oooo.... > oooo.... > oooo.... > oooo.... > > Could someone tell me how to write the ".ctl" fle? Really appreciate for > your great help! > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://ab-initio.mit.edu/pipermail/mpb-discuss/attachments/20101227/32467558/attachment-0001.htm > > > > ------------------------------ > > _______________________________________________ > mpb-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss > > End of mpb-discuss Digest, Vol 58, Issue 14 > ******************************************* >
_______________________________________________ mpb-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
