Hi, MPB assumes a periodic structure. So if you model a single unit cell and put in a dielectric block in one air hole, it will be repeated periodically. If you want to model a defect in a periodic structure, you need to use a supercell or go for a different simulation tool, which does not require periodicity (e.g. Meep).
Sebastian Dr Sebastian Schulz Postdoctoral Fellow CERC Quantum Nonlinear Optics Department of Physics University of Ottawa www.quantumphotonics.uottawa.ca Tel: 613-562-5800 (7138) > On Oct 25, 2015, at 12:00, mpb-discuss-requ...@ab-initio.mit.edu wrote: > > Send mpb-discuss mailing list submissions to > mpb-discuss@ab-initio.mit.edu > > 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 > mpb-discuss-requ...@ab-initio.mit.edu > > You can reach the person managing the list at > mpb-discuss-ow...@ab-initio.mit.edu > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of mpb-discuss digest..." > > > Today's Topics: > > 1. (no subject) (Ashina Garg) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 25 Oct 2015 01:57:21 +0530 > From: Ashina Garg <ashinagar...@gmail.com> > To: mpb-discuss@ab-initio.mit.edu > Subject: [MPB-discuss] (no subject) > Message-ID: > <CAJ8oNpmFY6D=eU__ri=eBJDdHnN=cT1rMHv7W-etr1=-onw...@mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > I want to create a line defect 1-D PhC waveguide. Below is my ctl file and > the image of epsilon file. I have a default material of 4.41 dielecric > constant and i have inserted blocks of air. To temove a single block, I > have inserted a block of dielectric in the centre of slightly less radius. > I want a single block of dielectric only in the centre which is of same > radius as all other blocks. But whenever the run the code, I'm getting a > dielectric in the centre of each block. Please help me find where am I > going wrong? > > (define-param x-dim 1) > (define-param y-dim 20) > (define-param z-dim 4) > > (set! geometry-lattice (make lattice(size 1 1 no-size))) > (set! default-material (make dielectric (epsilon > 4.41))) > (set! geometry(list(make block (material (make dielectric (epsilon 1))) > (center 0 0 0) (size 0.3 1 infinity)))) > (set! geometry (append geometry > (list (make block (material (make dielectric (epsilon 4.41))) > (center 0 0 0) (size 0.1 1 infinity))))) > > (set! resolution 16) > (set! k-points (list (vector3 0.5 0.5 0))) > (set! num-bands 50) > > > (run) > > Please help > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > <http://ab-initio.mit.edu/pipermail/mpb-discuss/attachments/20151025/970d05fc/attachment-0001.html> > -------------- next part -------------- > A non-text attachment was scrubbed... > Name: epsilon.png > Type: image/png > Size: 1021 bytes > Desc: not available > URL: > <http://ab-initio.mit.edu/pipermail/mpb-discuss/attachments/20151025/970d05fc/attachment-0001.png> > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > mpb-discuss mailing list > mpb-discuss@ab-initio.mit.edu > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss > > ------------------------------ > > End of mpb-discuss Digest, Vol 108, Issue 8 > ******************************************* _______________________________________________ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss