Re: [MPB-discuss] fortran won't play nice with C

On Mon, 22 May 2006, Jethro Bell wrote: checking for ctl_get_vector3 in -lctl... no configure: error: Couldn't find the required libctl library. [...] libctl.a does appear in the places where I think it is supposed to (/sw/lib/libctl.a and /sw/lib/libctl.la). You can look in the config.log

Re: [MPB-discuss] Question

On Tue, 23 May 2006, hbchen wrote: For a fixed frequency, I want to picture the corresponding field distribution (Such as E field, H field or energy) in a 2D or 3D photonic crystall. Using find-k function, one can find k at a given frequency. But how to output the corresponding field

Re: [MPB-discuss] what is wrong with my ctl structure??

On Sun, 11 Jun 2006, bldd180 wrote: I am doing some simulations on tetrahedral structure(tetrahedral network of cylindrical bonds connecting all of the nearest neighour points of the diamond lattice).And I wrote the CTL document, but the results I got are totally wrong. we all know this kind

Re: [MPB-discuss] triangular lattice

On Tue, 13 Jun 2006, Kai Lin Woon wrote: I am learning to use the programme. I dont know whether I have input my parameter correctly. Worst still I dont know how to interpret the output data. My structure is basically a triangular lattice air hole array slab. See the data analysis tutorial in

Re: [MPB-discuss] question about output file

in homogenization (effective medium) theory because they give the Wiener bounds [ Aspnes, Am. J. Phys. 50, 704-709 (1982) ]. Cordially, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman

Re: [MPB-discuss] anisotropic

On Sun, 20 Aug 2006, hlren wrote: This is 2D triangular lattice photonic crystal with air holes in GaAs substrate. The backgound dielectric index is 11.4. We calculalte the dispersion relations of line defect waveguides filled with liquid crystal, and the attached file is about the

Re: [MPB-discuss] Re: HDF5/MPB stopped working after installing MEEP

If you make a clean Linux install, I would recommend installing Debian, which has a prepackaged MPB and libctl and everything else you'll need, except for harminv and Meep. (You'll need to use Debian/testing to get the most recent version of libctl.)

Re: [MPB-discuss] Resonant frequencies in PC slabs!

On Thu, 2 Nov 2006, Ngo Quang Minh wrote: I would like to compute the defect modes (Resonant frequencies) in photonic crystal slab. (According to the document of MPB User tutorial - Finding a point defect state). I do not understand what my .ctl file is wrong! You have a fundamental

Re: [MPB-discuss] can MPB output the fields for arbitrary k-point?

) (0,0.153846,0) (0,0.134615,0) (0,0.115385,0) (0,0.0961538,0) (0,0.0769231,0) (0,0.0576923,0) (0,0.0384615,0) (0,0.0192308,0) (0,0,0) Any other suggestion? Thank you Emiliano Steven G. Johnson wrote: On Wed, 8 Nov 2006, Emiliano Di Gennaro wrote: (run-tm (output-at-kpoint (vector3

Re: [MPB-discuss] 3D doubts (partly cleared up)

On Fri, 17 Nov 2006, Vincent Paeder wrote: I've tried the case of a slab of silicon cylinders in n=1.5 medium (the previously attached ctl file), and I also end up with folded bands. So I guess MPB isn't most appropriate for treating slabs, or is there a way to discriminate the bands arising

Re: [MPB-discuss] What can I define the assymmetric photonic crystal slabs?

On Tue, 13 Feb 2007, Nils Jensen wrote: It seems MPB does not support the asymmetric structure, because I check my structure by using this command: h5topng -x (-y) The cross structure that I get, is not correct! MPB absolutely supports asymmetric structures. Perhaps you are using mpbi

Re: [MPB-discuss] MPB vs MPBI discrepancy

On Wed, 7 Mar 2007, MHerrera wrote: When I run the following ctl file using MPBI vs MPB, I get different results between the two, though I believe I should be getting the same output: I can't reproduce your problem. I modified your .ctl file to do only 4 bands, and only at k=0, with a

Re: [MPB-discuss] Point defects in slab structure - resonant frequencies

. For a point defect, there is no translational symmetry so you might as well use a rectangular supercell. Cordially, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb

Re: [MPB-discuss] Question about MPB configure for fftw

On Thu, 17 May 2007, Youri Bae wrote: configure: error: The FFTW libraries could not be found. Make sure FFTW is installed, and that LDFLAGS and CPPFLAGS are set appropriately if necessary. ~/mpb {118}$cd ../fftw-2.1.5 ~/fftw-2.1.5 {119}$pwd /u/home/baeri/fftw-2.1.5 As you can see, I

Re: [MPB-discuss] Error on configure for h5utils.

/index.php/H5utils Cordially, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss

Re: [MPB-discuss] Using MPB in TeraGrid

On Tue, 19 Jun 2007, Lingling Tang wrote: Our group obtained some computing allocation (30,000SUs) from TeraGrid, and we would like to install MPB (and MEEP) on it. The computers have different architecture and operating system, and we are trying to find suitable ones for installing MPB

Re: [MPB-discuss] output-dpwr-in-objects

On Wed, 27 Jun 2007, F.S.F. Brossard wrote: Thanks for that, I confused with (output-dpwr band) which can be called after the run command. You can call output-dpwr-in-objects after the run command too, for that matter, but you have to evaluate the result by passing it a band, e.g.

Re: [MPB-discuss] Asymmetric PC slab, TE/TM modes

to some modes more strongly than others, but computing this coupling coefficient exactly is not possible with MPB (there are some approximations you can use, but an exact calculation requires that you take evanescent modes into account). Regards, Steven G. Johnson

Re: [MPB-discuss] f77 compiler

On Thu, 19 Jul 2007, Dayan Handapangoda wrote: In order to install MPB, i need the f77 compiler to be installed. But couldn't find one. Found a windows version of it, but I need a one to be installed in Fedora 6. Where can I download it? Please help me on this... Install g77 or gfortran

Re: [MPB-discuss] MPB install error

On Sun, 22 Jul 2007, Dayan Handapangoda wrote: When i tried to use ./configure for MPB, the following error message appeared. I've installed g77, do i need to install f77 as well to proceed with MPB installation? I'm quite new to this, someone help me please... checking for g77... g77

Re: [MPB-discuss] Is a double layer(coated) rods 2D analysis possible?

.) Regards, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss

Re: [MPB-discuss] About the lightcone and projected slab mode

) (run-zeven ) (do-stuff (+ k k-inc) k-hi k-inc) ))) (do-stuff 0 0.5 0.1) Best wishes Ryan On Mar 28, 2007 2:20 AM, Steven G. Johnson [EMAIL PROTECTED] wrote: On Tue, 13 Mar 2007, Ryan Hao wrote: In the case of line defect waveguide in photonic crystal slab,I succeed in caculating

Re: [MPB-discuss] About the lightcone and projected slab mode

and project them. This requires some care, however, as described in my paper Linear waveguides in photonic-crystal slabs, in Phys. Rev. B. 62, p. 8212 (2000). It is easier to use the supercell technique above to get an orthogonal unit cell. Regards, Steven G. Johnson

Re: [MPB-discuss] About the lightcone and projected slab mode

On Fri, 16 Nov 2007, Ryan Hao wrote: I have used a non-orthogonal supercell for calvulating triangular lattice. That may be the problem. I chose the supercell like this: o o

Re: [MPB-discuss] negative/complex epsilon

as shown in the output file (petra.out) attached here. Also attached is the comp-eps.ctl file. Complex epsilon matrices are allowed, but they must be hermitian and positive definite -- the epsilon you have supplied violates this condition (it has a negative eigenvalue). Regards, Steven G. Johnson

Re: [MPB-discuss] PERFECT ELECTRIC CONDUCTOR

On Wed, 19 Dec 2007, abiola oladipo wrote: Does anybody know how to define a structure as a perfect electric conductor using the MPB code? MPB does not really support metals; you may need to use Meep (our FDTD code) instead. ___ mpb-discuss mailing

Re: [MPB-discuss] light direction

On Tue, 1 Jan 2008, andika asyuda wrote: My name is Andika Asyuda I am the new user of MPB I want to ask about how to set the direction of incidence of light in MPB There is no incident light in MPB. MPB computes eigenmodes. If you want a transmission spectrum or something like that,

Re: [MPB-discuss] creatin EFS

:_k_as_a_function_of_frequency However, for isofrequency plots I've found it much easier for the most part to just evaluate at a grid of k-points, and then pass it to Matlab to get the contour plot. Regards, Steven G. Johnson ___ mpb-discuss mailing list

[MPB-discuss] our new textbook on photonic crystals available online

-initio.mit.edu/book Hopefully, it will prove useful in your studies, teaching, and research, and in particular will help you to better understand the output (and input) of MPB and Meep. Regards, Steven G. Johnson ___ mpb-discuss mailing list mpb

Re: [MPB-discuss] mpbi: epsilon.h5 output?

On Mar 14, 2008, at 10:04 AM, Mischa Megens wrote: I am having trouble with mpbi (mpb-1.4.2 with inversion symmetry enabled). The layout of its epsilon.h5 output is funny. If I run mpb or mpbi on the example tri-rods.ctl file, it gives different results for epsilon.h5. Mpb gives the

Re: [MPB-discuss] define bock

On Mar 26, 2008, at 10:36 AM, andika asyuda wrote: I try to define triangular lattice structure by MPB, with block as basis but I do not now how to get the usual block. When I define e1,e2,e3 like following (e1 0.5 0.5 0) (e2 0.866 -0.866 0) I get thee rotated block rather than usual

Re: [MPB-discuss] help

On Apr 17, 2008, at 10:36 AM, mahmoud hamoud wrote: Hello dear, i'm trying to install h5utils-1.10.1 on my system gcc 3.4.3 madrakelinux 10.2 and when i tape make i got the following messages: arrayh5.c:382: warning: data definition has no type or storage class arrayh5.c:384: warning: data

Re: [MPB-discuss] output field

define the functions below, then you can just do: (run output-emag) Regards, Steven G. Johnson (define (get-emag which-band) (get-efield which-band) (let ((emag (rscalar-field-make cur-field))) (field-map! emag (lambda (f g) (sqrt (real-part (vector3-cdot f g

Re: [MPB-discuss] output field

On Apr 17, 2008, at 10:37 PM, Mindy Lee wrote: Thanks Dr. Johnson! I tried to run the code but it says ERROR: Unbound variable: get- filename-prefix And I can't find this function output-field-to-file in the manual... Whoops, just replace (get-filename-prefix) with filename-prefix (the

Re: [MPB-discuss] Relationship between primitive cell and reciprocal lattice vectors

in: http://www-math.mit.edu/~stevenj/18.369/spring07/midterm.pdf and the solution in http://www-math.mit.edu/~stevenj/18.369/spring07/midterm-sol.pdf) Your unit cell should be dictated by the periodicity of the structure, not by the propagation direction. Regards, Steven G. Johnson

Re: [MPB-discuss] mode volume

of the maximum in a local variable umax. Then you can just do: (run print-modal-volume) and it will output the modal volume via the above integral for every band and k point. Regards, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab

Re: [MPB-discuss] EM Fields and derivatives

second-order accuracy at best by differentiating planewaves, so you might as well just use a center- difference approximation. Regards, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman

Re: [MPB-discuss] mpb-discuss@ab-initio.mit.edu

) Then just do (run-tm (output-at-band 5 (output-at-kpoint (vector3 ...) output- hfield-z etc))) to output (e.g.) only fields for band 5 at the particular k-point. Regards, Steven G. Johnson ___ mpb-discuss mailing list mpb-discuss@ab

Re: [MPB-discuss] comparison of Guided Mode Expansion, MPB

On May 8, 2008, at 2:34 PM, Chris Kang wrote: My question pertains to the appearance of the bands in the light cone. While the gap seems to remain intact under the light cone, there seems to be some states that MPB is calculating that do not appear on the GME band plot. Is this merely a

Re: [MPB-discuss] need suggestions in creating EFS

In the past I've used the code below to do computations for wavevector diagrams. For example, to get a 2d wavevector diagram for the TM polarization in a square or rectangular lattice, you would do: (wavevector-diagram (kgrid 0 0.5 0 0.5 20 20) TM 0) instead of (run-tm), then grep the

Re: [MPB-discuss] installation with HDF5

On May 15, 2008, at 3:27 PM, Levine, Zachary wrote: I’m also having trouble getting HDF5 recognized.The first problem is that a file called “deflate” isn’t found in “-lz” or zlib.a. But I have a file /usr/local/zlib121-linux2.4/zlib.a which in turn has a file deflate.o in it as

Re: [MPB-discuss] Data analysis question

On May 22, 2008, at 4:30 PM, Mindy Lee wrote: Just one more question. If I use multiple color scheme, how do I print out the colorbar to indicate what color corresponding to what intensity scale? I posted one trick to get a colorbar at:

Re: [MPB-discuss] Data analysis question

On May 22, 2008, at 5:13 PM, Steven G. Johnson wrote: On May 22, 2008, at 4:30 PM, Mindy Lee wrote: Just one more question. If I use multiple color scheme, how do I print out the colorbar to indicate what color corresponding to what intensity scale? I posted one trick to get a colorbar

Re: [MPB-discuss] Implment Do Loop for Objects

On Jul 13, 2008, at 3:47 AM, Ian D. Hosein wrote: I'm trying to simulate the photonic band structure of an FCC lattice with random sphere sizes. I'm working with a supercell of size n (where n is a variable) and then I want to fill the cell with spheres on the lattice sites with

Re: [MPB-discuss] Installation Issue

On Jul 25, 2008, at 12:49 AM, Rakesh Mote wrote: I am facing following problem with MPB installation. I use Cygwin on windows. matrixio.c: In function `write_attr': matrixio.c:91: error: too few arguments to function `H5Dcreate2' matrixio.c:97: error: too few arguments to function

Re: [MPB-discuss] run test error

On Aug 6, 2008, at 3:16 PM, André Gonçalves de Lima wrote: elapsed time for initialization: 1 seconds. mpb: symbol lookup error: /usr/lib/3dnow/atlas/liblapack.so.3gf: undefined symbol: _gfortran_pow_r8_i4 Most likely you have multiple versions of Fortran installed on your computer (e.g.

Re: [MPB-discuss] Defect point inside slab phc

It is rather subtle to properly compute point-defect states for slabs in MPB, because the point-defect states are leaky rather than exponentially confined. This means that the artificial periodicity you impose with the supercell technique is not negligible and has a significant impact on

Re: [MPB-discuss] Slab BandStructure

On Aug 10, 2008, at 9:49 AM, majid sodagar wrote: Thanks for your response. I really have confused about calculating triangular slab bandstructure. I previously asked if the MPB supports PML. and the answer was NO.So how can we simulate slab band structure with it? You don't need PML to

Re: [MPB-discuss] Simulation on a monolayer

A monolayer of particles is mathematically analogous to a photonic- crystal slab. Like a slab, you handle it by using a unit cell in the plane of the particles, and a supercell in the vertical direction, and then looking at the modes only below the light line of the material above/below

Re: [MPB-discuss] extending field solutions with mpb-data

On Sep 14, 2008, at 7:42 PM, Matt Eichenfield wrote: Hi all, It seems like the program mpb-data -m M just tiles field solutions M times, ignorning the plane-wave phase of the solutions as you move from one unit cell to another. It needs to know the k-vector of the solution to extend

Re: [MPB-discuss] PhC from Meep tutorial with MPB

On Sep 15, 2008, at 10:18 AM, Elewout Hallynck wrote: . According to Meep, the first TE bandgap should lie between (approx.) 0.2 and 0.3 . However, when I try to calculate the band structure using MPB, I get the following results: You have to look at the bands lying only underneath the

Re: [MPB-discuss] Propagation constant

On Oct 28, 2008, at 10:19 PM, Yiling Qi wrote: Thanks for your kind response before, I have tried the find-k function, it could work well but still not completely satisfy my concern, since it can give out k(w), however what I would like to get is a relationship between kx/ky and

Re: [MPB-discuss] Propagation constant

On Oct 30, 2008, at 5:35 AM, Yiling Qi wrote: Dear Steven, My simulated structure is a 2D periodic photonic lattice (the contrast of refractive index is very small, the order is only around exp(-4). Then kx and ky are meaningful for the bandstructure, but apparently there will be

Re: [MPB-discuss] about E

On Dec 12, 2008, at 2:17 PM, Chad Husko wrote: 2) Book Look at the (free) online book of Steven and colleagues Ch. 7, where your question is treated in detail: http://ab-initio.mit.edu/book/ Effective areas are in chapter 9 (on fibers). ___

Re: [MPB-discuss] about the source type

Meep questions should go to the meep-discuss mailing list. (However, before posting, look at the archives: the question of the source direction has been discussed multiple times.) Yes, you can rotate and object like a block. In the case of a block, you would do it by rotating the

Re: [MPB-discuss] How to generate the gap map for given index and grid structures?

On Mar 3, 2009, at 7:28 PM, Lingyun Wang wrote: I am asking the same question that someone has asked before (as in this link: http://osdir.com/ml/science.photonic-bands/2004-01/msg2.html) . The link to Steven's example file in his reply post is no longer working. Could someone help me

Re: [MPB-discuss] How can I calculate of group velocity of guided mode?

On Mar 2, 2009, at 9:04 PM, Hassan Aghababaeian wrote: I need to calculate group velocity of guided mode (or Group index) in some geometry like line defects. From froum I understood display-group-velocities command in MPB can be useful. But I couldn't find what is output and how can be

Re: [MPB-discuss] How to set k-points in piperiod pc waveguide?

On Feb 20, 2009, at 4:02 AM, cef...@sina.com wrote: I design a PC waveguide , but the period of waveguide is not equal to the period of crystal .When i calculate the mode of the waveguide , how to set the k-point ? If the periods are commensurate, then the period should be the least

Re: [MPB-discuss] PBG of gyroid structure

On May 4, 2009, at 6:27 PM, Vinodkumar Saranathan wrote: Dear Steven and Mischa, Thanks for your suggestions, I was able to set up the dielectric function. However, I was trying to reproduce the result of Maldovan et al. 2002: http://prola.aps.org/abstract/PRB/v65/i16/e165123 where they

Re: [MPB-discuss] how the bandwidths vary with ky

Your code setting the k points looks fine to me, so it must be something wrong with the rest of your .ctl file. Here is the ctl file I used for that figure (with kx of either 0 or 0.5 for the different curves). (define-param eps 13) (define-param eps2 1) (define-param w 0.2) (set!

Re: [MPB-discuss] error message: Could not find colormap gray

On Jul 12, 2009, at 11:29 AM, andika asyuda wrote: When I used h5topng command, the following message appeared : Could not find colormap gray I saw several discussion about this. They suggested to use h5utils 1.11, but I still cannot find it. I found the newer version instead, version

Re: [MPB-discuss] help!

On Jul 8, 2009, at 10:29 PM, lfwang wrote: Dear Steven and MPB-users: I need to simulate a core-shell fcc structure ( core with dielectric constant epsilon1 + shell with dielectric constant epsilon2) by mpb. But I don't know how to define dielectric constant of spheres(core+shell).

Re: [MPB-discuss] Can't find the libctl header

On Jul 17, 2009, at 4:14 AM, Bin Shao wrote: e.g. what happens if you do gcc -E /usr/include/stdio.h ? It should run and spit out lots of output (the preprocessed stdio.h) not freeze or give an error. It works and spit out lots of output. But I have no idea what I can do next to solve

Re: [MPB-discuss] How to define and calculate the surface mode in 2D-PC of trianular lattice?

One difficulty with your computational cell is that you actually have two surfaces -- one on the left and one on the right. Remember that there are periodic boundary conditions. The solution is to intentionally have two equal surfaces that are mirror symmetric along the y direction, and

Re: [MPB-discuss] define k-path for unconventional waveguide in MPB

On Sep 2, 2009, at 11:44 PM, ryanhao wrote: == but how about the repeated unit takes 2 lattice constant? how to define the K-path? see attached. I have tried to decrease the reciprocal lattice to 1/2 of its previous value, but the

Re: [MPB-discuss] numerical stability of mpb and meep

On Sep 4, 2009, at 1:56 PM, Juan Pablo Vasco wrote: I have a question about mpb... and meep too: ¿how is the behavior of the convergence when i put a very large resolution? It varies between linear and quadratic with resolution, depending on the structure (e.g. sharp corners reduce the

Re: [MPB-discuss] mpb : how to duplicate a geometry object

On Sep 15, 2009, at 2:48 AM, Julien Arlandis wrote: (define air (make dielectric (epsilon 12))) (define whole1 (make cylinder (center -0.2 0 0) (radius 0.06) (height infinity) (material air) )) (define whole2 (make cylinder (center 0.3 0 0) (radius 0.04) (height infinity) (material air)

Re: [MPB-discuss] 4fold symmetry of bands in 6fold symmetry structure...

On Aug 20, 2010, at 8:05 AM, Joachim Kalden wrote: although i defined a hexagonal structure in epsilon, the bands have rather 4fold symmetry (see attachments). I just don't get the clou! (set! geometry-lattice (make lattice (size 4 4 no-size))) Your cell is square, not hexagonal. See the

Re: [MPB-discuss] dispersion diagram slab waveguide; problems with symmetry

On Sep 16, 2010, at 5:44 AM, m...@an3s.net m...@an3s.net wrote: I am trying to calculate the dispersion diagram for a photonic slab waveguide. For this I use a 1-dimensional computational cell in the z-direction. The waveguide itself is in the xy-plane and has a thickness of 1. In the

[MPB-discuss] Meep 1.2.1, MPB 1.5, libctl 3.2.2 are now released

. We have also moved all of the development sources onto Github, and you can use the GIthub issue-tracking system to file bug reports or feature requests: https://github.com/stevengj/mpb https://github.com/stevengj/meep https://github.com/stevengj/libctl Regards, Steven G

Re: [MPB-discuss] MPB usage within new MEEP release

On Apr 9, 2014, at 5:46 AM, Nikola Dordevic n...@zurich.ibm.com wrote: I downloaded fresh MEEP(1.2.1) and MPB(1.5) releases and successfully installed both of then. However, when I tried to use make (eigenmode-source) in MEEP, I got an error that I should configure/build MEEP with MPB.

Re: [MPB-discuss] Seeking help in validating MPB calculation for simple cube lattice

Dear All, I am beginner to Photonic bandgap calculations with MPB. I am trying to validate my MPB program with analytical solution given in Photonic bands: simple-cubic lattice, H. Sami Sözüer and Joseph W. Haus, Vol 10, Issue 2, 1993, but I am unable to replicate the Fig 1 and Fig 2 given

Re: [MPB-discuss] Is it possible to simulate the interface of two lattices? Similar to a 1D lattice with internal "structure"

> On Jul 27, 2017, at 4:40 PM, Chris Flower wrote: > I'm interested in doing bandstructure calculations for a photonic crystal of > two different lattice structures. For instance, imagine (in 2D) a lattice of > rods of radius R1 for all lattice sites with y>0, and radius R2

[MPB-discuss] Meep 1.4, MPB 1.6.1, Libctl 4.0 released

Meep and MPB users, New versions of Meep and MPB are now released. (They require a new version of Libctl, also just released.) The main new Meep feature is an official Python API, thanks to support from Simpetus LLC and funding from the National Science Foundation. We expect that this will

Re: [MPB-discuss] PhC Slab

MPB has periodic boundary conditions, so the modes it calculates above the light line (which are not localized and hence "see" the artificial periodic boundaries in the direction orthogonal to the slab) are essentially artifacts of these boundary conditions. So, when studying waveguides and