[Meep-discuss] Question about Banddiagram Tutorial
Dear Meep developers and users, I am reading the online tutorial of banddiagram and seeing that the frequency is normailized to a number without unit. Supposing the unit length corresponds to 100nm, what's the actual freq of 0.5? Does it mean f=0.5*c/100 and the angular freq=2*pi*0.5*c/100 with unit of Hz? Thank you Gu ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] how to draw a core-shell cylindrical object?
Dear All, I'd like to draw a core-shell cylindrical object: a cylinder of material A and a circular shell of material B. Could anybody give me some help? Thank you Zongquan Gu ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] A confusion about Gaussian Src
Dear Meepers, I am a littel confused by the Gaussian Src in meep. For example, assuming all the wavelength are normalized to 50nm, I have a Gaussian Src fcen=600nm, ranging from 600nm-700nm, in the code, I write as below, (define-param fcen 50/650=0.077) but now how I could define the fwidth? Thanks Gu___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] Periodicty: Steven said before K=0 at that periodic direction, while tutorial set k=0.4
Dear Steven and Meepers, I red serveral previous mails regarding Periodic Boundary Condition. Some replies from Steven said K should be 0 in that direction, while on website "Band Diagram" tutorial sets K=0.4 to exam the field at K=0.4. Is there any conflict between these setting? Thanks Little Gu___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] how to write k points in non-uniform cell?
Deer Meepers, If there is a 3D cell where X=10, Y=10, Z=100 repeating itself in XY plane, because the cell isn't uniform, how could I write the K points if I want to find the solution from(0,0)>(5,0)>(5,5)>(0,0) at Z=50? Thanks Little Gu ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] How to set Periodic Boundary Condition in XY plane of 2D photonic crystal
Deer Meepers, I'd like to simulate the Ez/Hz distribution ina real 2D photonic crystal, which means it has finite length and a substrate. I set PML in Z direction and want to set Periodic Boundary Condition(PBC) of the cell in XY plane, Here is my question; How to set PBC in XY plane or just leave it as blank like the tutorial on the website? Below is the code: ; all normalized to 50nm (define-param n 4.7) ; index of inner core of cylinder and block substrate (define-param i 4.3) ; index of outer shell of cylinder (define-param w 0.4) ; width of waveguide (define-param r 0.2) ; inner radius of ring (define-param pad 2) ; padding between (cylinder top, block bottom) and edge of PML in Z direction (define-param gap 0.4) ; padding between cylinder and cell edge in XY plane (define-param dpml 1.5) ; thickness of PML (define sxy (* 2 (+ r w gap))) ; cell size in XY plane (set! geometry-lattice (make lattice (size sxy sxy 28))) (set! geometry (list (make cylinder (center 0 0 0.5) (height 20) (radius (+ r w)) (material (make dielectric (index i (make cylinder (center 0 0 0.5) (height 20) (radius r) (material (make dielectric (index n (make block (center 0 0 -10) (size 2 2 2) (material (make dielectric (index n)) (set! pml-layers (list (make pml (direction Z)(thickness dpml (set-param! resolution 20) (define-param fcen 0.0806) ; pulse center frequency (define-param df 0.1) ; pulse width (set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth df))) (component Ez) (center 1 0 0 (run-sources+ 300 (after-sources (harminv Ez (vector3 1 0 0) fcen df))) (run-until 200 (in-volume (volume (center 0 0 0.5) (size 2 0 21)) (at-beginning output-epsilon) (at-end output-efield-z))) Any suggestion is appreciated. Thanks a lot Little Gu ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] Periodic Boundary Condition is default?
Dear Meepers, I red the tutorial "Band diagram" on the website and it specifies PML in Y direction. So the boundary condition in X direction is Periodic Boundary Condtion. I'd like to know if I don't set any boundary condition in this example, the default condition Meep sets in it will be Periodic Boundary Condition in both X and Y direction right? Thanks Gu___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] How to view epsilon from X or Y direction
Greetings, I have one Cylinder with a Block substrate in a 3D cell. (set! geometry-lattice (make lattice (size sxy sxy 14))) (set! geometry (list (make cylinder (center 0 0 0.5) (height 10) (radius r) (material (make dielectric (index n (make block (center 0 0 -5) (size 1 1 1) (material (make dielectric (index n ) ) I am able to view epsilon in different Z plane and it is correct, However, I'd like to view Epsilon from X direction or Y direction (horizental) rather than Z direction(vertical) I write the sentence like this: (in-volume (volume (center 0 0 0) (size sxy sxy 14))(at-beginning output-epsilon) But it doesn't work. Could you help to see how could I do that? Thanks Gu ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] dielectric function of GaAs
Dear all, Is there any reliable resource of introducing the dielectric function of GaAs including dispersion, imaginary part? Any suggestion or recommendation is appreciated. Thanks Little Gu___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
