Re: [Meep-discuss] error when defining geometry
On Nov 2, 2010, at 8:42 AM, Silviu Popescu wrote: Hello, I have an error in my code but I can't figure it out. (set! geometry (append (list (make block (center 0 0 (- -1.85 grosimeRetea)) (size infinity infinity 3.4) (material (make dielectric (epsilon eps_SiO2))) ) ) (geometric-object-duplicates (vector3 0.5 0 0) 0 (/ Sx 0.5) (make block (center (0 0 (/ grosimeRetea -2) ) ) (size umplere infinity grosimeRetea) (material (make dielectric (epsilon eps_Si))) ) ) ) ) :16:80: In expression (0 0 (/ grosimeRetea -2)): :16:80: Wrong type to apply: 0 ABORT: (misc-error) Can someone give me a tip ? The expression (0 0 foo) in Scheme means to call the function "0" with arguments 0 and foo. Since 0 is not a function, this gives an error. Fix this by getting rid of your extra parens: (center 0 0 ...) Moral: parentheses in Scheme are important! --SGJ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
Re: [Meep-discuss] strange behaviour when setting geometry-lattice parameters
On Nov 2, 2010, at 11:43 AM, Silviu Popescu wrote: Dear meep users, I don't understand the following behaviour. When I define geometry lattice like this: (define-param Sx 8) (define-param Sz 8 ) (set! geometry-lattice (make lattice (size Sx no-size Sz))) (set-param! resolution 10) \ This is because Sx and Sz are predefined constants in Meep referring to components of the Poynting vector (which happen to have values 100 and 104). The define-param command does nothing if the variable is already defined (in order to handle the case where the variable was defined on the command line), and so it uses the predefined values in this case. --SGJ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] strange behaviour when setting geometry-lattice parameters
Dear meep users, I don't understand the following behaviour. When I define geometry lattice like this: (define-param Sx 8) (define-param Sz 8 ) (set! geometry-lattice (make lattice (size Sx no-size Sz))) (set-param! resolution 10) At the output I have: Initializing structure... Working in 3D dimensions. Computational cell is 100 x 0 x 104 with resolution 10 block, center = (0,0,-2.15) . And for (set! geometry-lattice (make lattice (size 8 no-size 8))) (set-param! resolution 10) . At the output I have: Initializing structure... Working in 3D dimensions. Computational cell is 8 x 0 x 8 with resolution 10 block, center = (0,0,-2.15) . Can someone explain this behaviour ? I use: Meep 1.1.1, Copyright (C) 2005-2009 Massachusetts Insitute of Technology. Using libctl 3.0.3 and Guile 1.8.7. Best Regards, Silviu P. ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] gaussian source output propagation
Hi All, Thanks for MEEP, I have been finding it very useful, and am currently in the process of porting it over to our cluster. Looking forward to it! However, this time I believe I am using MEEP for a possibly inappropriate purpose; I was hoping I could describe the situation and find out whether I can believe the results. I am attempting to simulate the light propagation of a VCSEL, out of the aperture, through a small air gap, substrate, cladding material, and into a waveguide core. Sparing you the details (unless they are needed), I have modeled this as a 2D structure that would look exactly as a cross-section of this stack up would be drawn. The VCSEL is modeled as a 1-D gaussian source, with the length equal to the aperture diameter. All of these materials are transparent at the wavelength of interest, I am using only a real epsilon for each material. There are several hard index discontinuities. I am getting essentially no propagation through the stack. Is the source a reasonable way to represent the VCSEL? I know the divergence angle of the VCSEL, but don't know how to put collimation into the source. Is this whole situation something that I cannot expect MEEP to help me with? Thanks, Chris ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] Verification of material (Au) dispersion
Dear MEEP users, I get a Drude model for Au (the model itself has been proved to be accurate with applicable wavelength range: 500~1000nm), and want to verify it through simulaltion as the MEEP Tutorial of Material dispersion. I made a slight change to the given tutorial .ctl file. However, the result achieved is not what I have expected. The .ctl command and partial simulation results are as follows. Any statement about my problem will be highly appreciated. With best regards. *** The ctl command: ; Define lengthscale a = 200nm ; Wavelength range : 500~1000nm ; Corresponding frequency range in MEEP : 0.2 ~ 0.4 (set! geometry-lattice (make lattice (size no-size no-size no-size))) (set-param! resolution 20) (set! default-material (make dielectric (epsilon 5.9673) (E-polarizations (make polarizability (omega 0.67e-10) (gamma 0.0106133) (sigma 4.4673e20)) (make polarizability (omega 0.43338) (gamma 0.0699067) (sigma 1.09)) ) ) ) (define-param fcen 0.3) (define-param df 0.2) (set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth df))) (component Ez) (center 0 0 0 (define-param kmin 0.2) (define-param kmax 0.4); I am not sure if the values for kmin and kmax are right. However, even if I change these values, the results and the problem are the same. (define-param k-interp 99) (define kpts (interpolate k-interp (list (vector3 kmin) (vector3 kmax (define all-freqs (run-k-points 200 kpts)) ; a list of lists of frequencies (map (lambda (kx fs) (map (lambda (f) (print "eps:, " (real-part f) ", " (imag-part f) ", " (sqr (/ kx f)) "\n")) fs ) ) (map vector3-x kpts) all-freqs ) * Partial simulation results: ... harminv0:, frequency, imag. freq., Q, |amp|, amplitude, error run 0 finished at t = 250.0 (1 timesteps) freqs:, 1, 0.2, 0, 0 freqs-im:, 1, 0.2, 0, 0 harminv1:, frequency, imag. freq., Q, |amp|, amplitude, error run 1 finished at t = 250.0 (1 timesteps) freqs:, 2, 0.202, 0, 0 freqs-im:, 2, 0.202, 0, 0 ... ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] error when defining geometry
Hello, I have an error in my code but I can't figure it out. (set! geometry (append (list (make block (center 0 0 (- -1.85 grosimeRetea)) (size infinity infinity 3.4) (material (make dielectric (epsilon eps_SiO2))) ) ) (geometric-object-duplicates (vector3 0.5 0 0) 0 (/ Sx 0.5) (make block (center (0 0 (/ grosimeRetea -2) ) ) (size umplere infinity grosimeRetea) (material (make dielectric (epsilon eps_Si))) ) ) ) ) :16:80: In expression (0 0 (/ grosimeRetea -2)): :16:80: Wrong type to apply: 0 ABORT: (misc-error) Can someone give me a tip ? Thank you, Silviu Popescu___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] About cylinder of geometric-object.
To come straight to the point,About cylinder of geometric-object. The result becomes a jet-black screen if it analyzes it by wanting to use a cylinder thing, and using cylinder. Are not there any solutions? My best regards. SHO ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] About cylinder of geometric-object
Delivery to the following recipient failed permanently ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] About cylinder of geometric-object
To come straight to the point,About cylinder of geometric-object. The result becomes a jet-black screen if it analyzes it by wanting to use a cylinder thing, and using cylinder. Are not there any solutions? My best regards. SHO ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] About cylinder of geometric-object.
Hello. To come straight to the point,About cylinder of geometric-object. The result becomes a jet-black screen if it analyzes it by wanting to use a cylinder thing, and using cylinder. Are not there any solutions? My best regards. SHO ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
[Meep-discuss] ,About cylinder of geometric-object.
Hello. To come straight to the point,About cylinder of geometric-object. The result becomes a jet-black screen if it analyzes it by wanting to use a cylinder thing, and using cylinder. Are not there any solutions? My best regards. SHO ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
Re: [Meep-discuss] Question about (run-k-point s...) and nonlinear material.
Thank you, Prof. Steven. I've got some ideas about this problem: I see a small down movement of the branch of the dispersion diagram in the Kerr-nonlinear periodic structure. At this moment with (run-k-points...) method I just see this movement when the ampliude of the signal is not big enough. And when I started to simulate the structure without periodical boundary condition I started to observe gap solitons, as you write in your letter. Now I tried to understand where is the scope of applicability of my simulation. Of course it is very important question as well. Thank you, Alexandr. ___ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss