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
Re: 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:44Outputting Fields at One k-point for One Band
http://permalink.gmane.org/gmane.comp.science.photonic-bands/1990
<pre>Hi everyone,
I'm sorry if this question is basic, but after searching for a while I really couldn't figure out an answer.
My question is this. Given a dielectric structure for which we find the band structure over a given path over the first Brillouin Zone (e.g. Gamma-X-M-Gamma), how do we get the fields at just one k-point and just one band? For example, if I want the fields from the 33rd band at the X-point.
I have tried:
(run (output-at-k-point (vector3 0.5 0 0) (output-efield-x 33) (output-efield-y 33) (output-efield-z 33)))
but this doesn't work. I saw a post from Steven Johnson at https://www.mail-archive.com/mpb-discuss< at >ab-initio.mit.edu/msg00062.html explaining that we need to give the run function a band function, so I tried making a function as follows:
(define (MyBandFunction)
(begin
(output-at-k-point (vector3 0.5 0 0)
(output-efield-x 33)
(output-efield-y 33
(output-efield-z 33)
)
)
)
(run MyBandFunction)
but still no dice. By my understanding of how MPB works is that it takes your dielectric structure, takes a band and outputs the electric field for all the k-points. So, I figured I could use the "current-k" global variable. Something like
(define (MyBandFunction2)
(if (eq? current-k (vector3 0.5 0 0))
(
(begin
(output-efield-x 33)
(output-efield-y 33
(output-efield-z 33)
)
)
)
)
(run MyBandFunction2)
but it still doesn't work. Would someone please be able help me figure out how to get the output for one particular k-point at one band?
Thank-you very much!
Pranai
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>Pranai Vasudev2014-03-14T02:57:25Re: Photonic crystal slab waveguide
http://permalink.gmane.org/gmane.comp.science.photonic-bands/1989
<pre>Hi,
I am trying to find the line defect mode for square lattice of size 7X7. I
am not understanding how to take the k-points
This is what I took:
(set! k-points (list (vector3 0 0 0) ; Gamma
(vector3 0.5 0 0))) ; X
(set! k-points (interpolate 16 k-points))
Is it correct?
</pre>Akash2014-02-14T06:19:40Re: mpb ./configure fails with Fedora 20 x86_64
http://permalink.gmane.org/gmane.comp.science.photonic-bands/1988
<pre>
Hi
I managed to make mpb 1.4.2 on Fedora 20:
- install compat-guile18 and compat-libf2c-34
- replace guile-config with guile1.8-config (configure must find script for
1.8 first)
- make your on libctl:
./configure --prefix=/opt/mpb-1.4.2/libctl LDFLAGS="-lm"
- config mpb like this
export gen_ctl_io_ok=/opt/mpb-1.4.2/libctl/bin/gen-ctl-io
./configure --with-libctl=/opt/mpb-1.4.2/libctl/share/libctl
--prefix=/opt/mpb-1.4.2 CPPFLAGS="-I/usr/include/ctl
-I/opt/mpb-1.4.2/libctl/include -DH5_USE_16_API=1"
LDFLAGS="-L/opt/mpb-1.4.2/libctl/lib"
I found in
http://thread.gmane.org/gmane.comp.science.photonic-bands/1902/focus=1904
a hint about a prerelease version 1.5. Which I could build with fewer
tweaks:
./configure --prefix=/opt/mpb-1.5 CPPFLAGS="-I/usr/include/ctl"
Im not the end user, my customers are currently testing the prerelease, so
I cannot talk about the output of the prerelease version.
Regards
Bengt
_______________________________________________
mpb-discuss mailing list
mpb-discuss< at >ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss</pre>tim.5.bigi_gigi< at >xoxy.net2014-02-05T11:54:14Re: mpb ./configure fails with Fedora 20 x86_64
http://permalink.gmane.org/gmane.comp.science.photonic-bands/1987
<pre>Hi,
guile as well as quile-devel are installed by the script.
</pre>René Genz2014-02-02T17:46:46Search EngineSearch the mailing list at Gmanequery
http://search.gmane.org/?group=$group=gmane.comp.science.photonic-bands