gmane.comp.science.photonic-bands
http://blog.gmane.org/gmane.comp.science.photonic-bands
hourly11901-01-01T00:00+00:00Gmanehttp://gmane.org/img/gmane-25t.png
http://gmane.org
mpb build fails on openSuSE due to ./configure vs.libctl
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2021
<pre>Hi everyone,
I'm reporting a problem (and solution, though) with compiling mpb on
openSuSE 12.3.
Having all required libraries installed, I get the following error on
running ./configure:
checking for libctl dir... /usr/local/share/libctl
checking for gen-ctl-io... gen-ctl-io
checking for ctl_get_vector3 in -lctl... no
And then configure complains about not finding the libctl required.
(Sorry, didn't catch the message exactly, only its meaning above.)
Although mpb finds libctl, as you can see above. Digging into config.log
I found out, that in the check the test program is not linked with
libctl. -lctl is given, but not -L/usr/local/lib. "Patching" ./configure
in line 20233 worked perfectly:
LIBS="-L/usr/local/lib -lctl $LIBS"
I didn't find the place, where LIBS lost -L/usr/local/lib, although it
is given in the environment variable LDFLAGS. Or maybe it is a bug the
openSuSE environment.
Sincerely yours,
Robert Bienert
Universität Kassel
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Robert Bienert2015-02-19T14:51:56Re: MPB basic on non zero kz values
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2020
<pre>Thanks Sir for your reply.
This is solved.
Now, what I'm facing that for a non zero kz value, there will be no TE or
TM polarization. The fields for a particular band at a certain point should
be coupled. is there any way in mpb to map that coupled field? (I tried
that but the results seem to be not reliable)
(eg. if for an inplane condition, if we search for e field distribution for
a band (say band 2) at a point (say K point) of a a hexagonal PC, mpb
returns e field distribution)
Thanks
Subir
Subir Majumder
On Thu, Feb 12, 2015 at 11:30 PM, silviamgr< at >ua.pt <silviamgr< at >ua.pt> wrote:
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Subir Majumder2015-02-13T12:14:21Re: MPB basic on non zero kz values
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2019
<pre>Hello,
if you need to put non-zero values for kz in MPB, for instance in an hexagonal lattice, it can be done like this:
(set! k-points (list (vector3 0 0 kz_mpb) ; gamma
(vector3 0 0.5 kz_mpb) ; m
(vector3 (/ -3) (/ 3) kz_mpb) ; k
(vector3 0 0 kz_mpb))) ; gamma
where, kz_mpb in fact is unitless, it is:
kz_mpb = kz.a/(2pi),
where kz can be measured in [m^-1], and a is the lattice periodicity in [m].
Sílvia M.G. Rodrigues
________________________________________
From: mpb-discuss <mpb-discuss-bounces< at >ab-initio.mit.edu> on behalf of Subir <subir.majumder.23< at >gmail.com>
Sent: Friday, December 26, 2014 10:18 AM
To: mpb-discuss< at >ab-initio.mit.edu
Subject: [MPB-discuss] MPB basic on non zero kz values
Sir, I'm a new user of MPB. I work on photonic crystal fibers. My simple
question is what is the unit of kz values used in the package if need to put
non zero values of kz.
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
</pre>silviamgr< at >ua.pt2015-02-12T18:00:59MPB basic on non zero kz values
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2018
<pre>Sir, I'm a new user of MPB. I work on photonic crystal fibers. My simple
question is what is the unit of kz values used in the package if need to put
non zero values of kz.
</pre>Subir2014-12-26T10:18:36Re: swg
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2017
<pre>Hi,
As far as I know MPB will solve the dispersion curves and the group
velocities when given a propagation direction (k-point) and a structure,
assuming the latter is periodical. So, you'll may end up simulating a
lattice of several segmented waveguides but if they are separated enough
there shouldn't be any problem.
Also the units don't matter that much since you'll have to specify only the
segmentation related to the length of the unit cell (the duty cycle) since
the program uses a normalized unit system. Maybe if you specify more about
what you are trying to simulate I could help more.
Best regards
2014-11-12 14:36 GMT-06:00 Ana Júlia <anajrfoliveira< at >gmail.com>:
Rafael
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Rafael Gonzalez2014-11-12T22:39:02swg
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2016
<pre>hello all
i would like to simulate a 3d SWG segmented waveguide, with a period of
300nm and dutycycle of 50%. How can i (and... is it possible to) obtain the
effective index, the group effective index and dispersion curves?
thanks
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Ana Júlia2014-11-12T20:36:45Supercell si,mulation
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2015
<pre>Hello,
I try to run the 2D hyperuniform structure with rod (raduis 0.2a) in 1600 bands for TE or TM mode. The basis is 22 by 22.
What's the limitation for resolution in the simulation?
Is there any limitation in grid size for the MPB ?
Thank you very much.
?
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Hong Zhu2014-09-21T00:41:18help
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2014
<pre>
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>behrooz rezaei2014-09-16T05:42:54Calculations of dpwr - is it |E|^2 without eps?
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2013
<pre>Hello Prof. and MPB users,
I recently resumed my interest in MPB and began running simple
calculations. I now look at 3D simple cubic grid of dielectric spheres.
I run a ctl file where I calculate the bands, the fields and the power
density.
The title suggests that I might be overlooking something simple, but I will
elaborate on my problem:
;define the cell geometry of spheres in cubic grid
(define-param eps 11.9000000000000000000000000000000000000000000000000)
(define-param r 0.2)
(define-param dielectricon
(make dielectric
(epsilon eps)
)
)
(set! default-material air)
;define the lattice vectors of cubic lattice
(set! geometry-lattice
(make lattice
(basis-size 1 1 1)
(basis1 1 0 0)
(basis2 0 1 0)
(basis3 0 0 1)
)
)
(set! geometry
(list
(make sphere
(center 0 0 0)
(radius r)
(material dielectricon)
)
)
)
;define the k-points
(set! k-points
(list
(cartesian->reciprocal (vector3 0.2 0.3 -0.1))
)
)
;define the simulation settings
(set-param! resolution 32)
(set-param! num-bands 1)
(define-param top-bands 1)
(set-param! tolerance 0.0000000001)
;disable default filename prefix in output files
(set! filename-prefix false)
;special calculation functions to pass to run
(define (output-nonbloch-efield which-band)
(get-efield which-band)
(cvector-field-nonbloch! cur-field)
(output-field-to-file -1 "e.")
)
(define (output-nonbloch-dfield which-band)
(get-dfield which-band)
(cvector-field-nonbloch! cur-field)
(output-field-to-file -1 "d.")
)
(define (output-nonbloch-hfield which-band)
(get-hfield which-band)
(cvector-field-nonbloch! cur-field)
(output-field-to-file -1 "h.")
)
(define (topbands-functions)
(let topbands-loop ((n (- top-bands 1)))
(output-nonbloch-efield (- num-bands n))
(output-nonbloch-dfield (- num-bands n))
(output-nonbloch-hfield (- num-bands n))
(output-tot-pwr (- num-bands n))
(output-hpwr (- num-bands n))
(output-dpwr (- num-bands n))
(if (> n 0) (topbands-loop (- n 1)) 0)
)
)
(run topbands-functions)
I than use matlab to read h5 files of the fields, compute the power density
and write a new h5 file, then I use vis5d to visualize the the power
densities.
U = mpb.epsilon .* reshape(sum(abs(fields.E).^2), size(mpb.epsilon)) ... %
fields.E is 4D matrix: 3xNxNxN, sum() will sum over the 3 components x,y,z
squares Fx^2+Fy^2+Fz^2
+ reshape(sum(abs(fields.H).^2), size(mpb.epsilon)); % epsilon is NxNxN
matrix
filename = [directory '\tot.rpwr.k' sprintf('%02d',k) '.b01.h5']; % or hpwr
or dpwr in order to look at either the elect. or magn. energy
Umpb = hdf5read(filename,'data');
file_new = [filename(1:end-3),'.U',filename(end-2:end)];
copyfile(filename, file_new);
hdf5write(file_new,'/data',U);
What I do not understand is as follows:
1. When I compare in vis5d (1) mpb output of the energy distribution and
(2) my calculated energy distribution, I see that dpwr matches
sum(abs(fields.E).^2) without multiplication in epsilon. It seems that
output-nonbloch-efield gives D instead of E
2. when I check normalization of the fields
using sum(abs(fields.F(:)).^2)*mpb.dV (dV being the spatial resolution of
volume 1/N^3), I get that H is normalized to 1, E is normalized to 0.965
and D is normalized to 1.3333. According to the User Reference, D should
also be normalized to 1, right?
I checked my code again and again, so either there is a bug or I missed
something too simple. Any help will be appreciated.
Thanks in advance for the help.
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Amnon Willinger2014-09-14T20:54:24MPB on OSX 10.9.4
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2012
<pre>Dear Community,
I am interested in installing MPB on my OSX 10.9.4, but currently it is
unsupported. I have, on the other hand, successfully installed meep using
FINK. I used the same procedure in FINK to install MPB but I get the
following error.
Warning! ***HDF5 library version mismatched error***
Thanks,
Ed
*--*
*PhD Candidate*
*Quantum Photonics Lab*
*Research Lab for Electronics, MIT*
*Office: 36-569*
*Scholar Profile
<http://scholar.google.com/citations?user=d5_MmyoAAAAJ&hl=en>*
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Edward H. Chen2014-09-02T11:38:46Is it necessary to fix fields phases when using(output-poynting)?
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2011
<pre>_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Clara Gutierrez Martin2014-08-27T08:54:21MPB replaces y with x
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2010
<pre>Dear all,
I try to use mpb on a server with the following information for mpb:
mpb --version
MIT Photonic Bands 1.4.2, Copyright (C) 1999-2001 by MIT.
Using libctl 3.2.2 and Guile 1.8.8.
The code that I have works well with the mpb that I have on my own ubuntu machine with:
mpb --version
MIT Photonic Bands 1.4.2, Copyright (C) 1999-2001 by MIT.
Using libctl 3.1 and Guile 1.8.8.
When I try to use mpb on the server I see that it replaces x with y for example I define a block with the size of (dx dy dz) but I see in the result that the output is a block with the size of (dy dx dz).
Do you have any idea about this problem?
Thanks in advance,
Ehsan
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Ehsan Saei2014-08-19T13:55:35Group velocity at second and higher bands
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2009
<pre>Dear MPB users and developers.
I tried to check MPB group velocity calculation correctness. At
photonic crystal with negligible contrast group velocity should be
roughly equal light speed at material with changing sings because of
folding. But calculation shows that absolute value of group velocity
depends of band number.
MBP code:
(set! num-bands 7)
(set! geometry-lattice (make lattice (size 1 1 no-size)
(basis1 (/ (sqrt 3) 2) 0.5)
(basis2 (/ (sqrt 3) 2) -0.5)))
(set! default-material air)
(set! geometry (list (make cylinder
(center 0 0 0) (radius 0.2) (height infinity)
(material (make dielectric (epsilon 1.00000001))))))
(set! k-points (list (vector3 0 0 0) ; Gamma
(vector3 (/ -3) (/ 3) 0) ; K
))
(set! k-points (interpolate 64 k-points))
(set! resolution 64)
(run-te display-group-velocities )
; END OF MPB CODE
Example velocity output string:
tevelocity:, 22, #(7.78615212026757e-25 -0.999999999272597 0.0),
#(3.40122036260571e-7 0.805441667563366 0.0), #(-3.40122036718879e-7
0.805441669075089 0.0), #(1.18935285286298e-8 -0.183365906857283 0.0),
#(-1.18935288080422e-8 -0.183365906656533 0.0), #(1.4273767716188e-7
-0.90326523923235 0.0), #(-1.42737677309561e-7 -0.90326524059479 0.0)
First band is ok, but absolute value of second band group velocity
near 0.8 not 1.0.
What is wrong?
Alex Friman.
</pre>Alex2014-08-11T15:37:11Auto Response
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2008
<pre>_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Farshid koohi2014-08-06T16:01:08Is it possible to calculate the band structure of a 3D geometry which is only periodic in 1 or 2 dimensions?
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2007
<pre>Hi,
Is it possible to calculate the band structure of a 3D geometry which is
only periodic in 1 or 2 dimensions?
(i.e. for example a hexagonal lattice of holes in a slab of a given
thickness or a fixed number of layers of a 3D crystal.)
Or is it necessary to use a 3D lattice with a basis3 vector much longer
than the thickness of the structure, i.e. a sort of supercell?
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Mike Taverne2014-08-05T16:12:37Effective index method for solving 1D photoniccrystal slab waveguide
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2006
<pre>Dear Dr. Johnson,
My name is Changyi Li. I am a PhD student at the University of New Mexico.
I would like to calculate the band diagram of a 1D photonic crystal on a
nanobeam (a free-standing rectangular GaN waveguide). The photonic crystal
would be a line of holes.
I have read some papers about the effective index method, which convert the
3D problem to 2D problem. But some paper say that this method is not
accurate enough, at least in some cases. So would you please give me some
suggestion? Should I just solve the 3D problem in mpb or use the effective
index method? And also would you please recommend me some reference if you
have time?
Thank you so much.
Best,
Changyi
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Li .Changyi2014-07-10T16:13:21Elliptical cylinder shape in MPB
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2005
<pre>Hi,
I'm looking to calculate the band structure of an unconventional 3D
photonic crystal which has a simple cubic arrangement. The repeated unit
cell is an open octahedron and I would like to build this structure such
that the edges of the octahedron are made of dielectric elliptical
cylinders. As such, I was wondering whether an elliptical cylinder type
geometric object has been added to the most recent version of MPB. If not,
I would appreciate any advice you can offer with regards to building an
open octahedron unit cell.
Thank you for all your help!
Best,
Viki
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Victoria Chernow2014-07-09T03:09:27Re: simulate hcp inverse opal structure
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2004
<pre>
</pre>pb-student2014-06-18T12:50:01Compillation Error: maxwell_eps.c(307): error: identifier "num_sphere_quad" is undefined
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2003
<pre>_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Lars Liebermeister2014-06-11T14:50:21Re: problem in compiling mpb 1.5
http://permalink.gmane.org/gmane.comp.science.photonic-bands/2000
<pre>Hi
I would recommend to check if the configure script finds the HDF5 library
because when it doesn't the compiling process succeds but with limited
functionality. That might be the problem
Rafael
2014-05-07 0:13 GMT-05:00 Ehsan Saei <e.saei< at >hotmail.com>:
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Rafael Gonzalez2014-05-07T15:37:06problem in compiling mpb 1.5
http://permalink.gmane.org/gmane.comp.science.photonic-bands/1999
<pre>Dear community,
I tried to install mpb 1.5 but when I configure it with "./configure --with-mpi" I get this warning: hdf5.h: accepted by the compiler, rejected by preprocessor
and then I install it with "make install" command. I can find the binary file for mpb-mpi but it doesn't work at all. I tried to test it just --version option but even for that it responds with segmentation fault error
Does this behavior relate to hdf5?
Thanks in advance,
Ehsan
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Ehsan Saei2014-05-07T05:13:44Search EngineSearch the mailing list at Gmanequery
http://search.gmane.org/?group=$group=gmane.comp.science.photonic-bands