tag:blogger.com,1999:blog-852400474287196648.comments2015-02-27T07:58:09.364-08:00Majid's ResearchMajid al-Dosarihttps://plus.google.com/111269600501787951096noreply@blogger.comBlogger21125tag:blogger.com,1999:blog-852400474287196648.post-4043886135301678112015-02-27T07:58:09.355-08:002015-02-27T07:58:09.355-08:00I too had the same feeling after using MEEP-Python...I too had the same feeling after using MEEP-PythonVikraman Knoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-43261879207571059652014-05-03T13:44:42.621-07:002014-05-03T13:44:42.621-07:00Sorry John but my priorities have shifted signific...Sorry John but my priorities have shifted significantly away from this work. Did you go through the math and code of this page?Majid al-Dosarihttp://about.me/majidaldosarinoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-84912923554893363812014-04-29T06:26:31.870-07:002014-04-29T06:26:31.870-07:00Hi Majid,
I have a question about how to define a...Hi Majid,<br /><br />I have a question about how to define a material by pure Drude model with Lorentzian oscillating term. Let say a material epsilon(omega)= eps_inf + sigam/(omega*(omega+i*gamma)) with eps_inf=12,sigma=0.3(2*pi*c/a)²,gamma=1.6(2*pi*c/a). I am looking forward to your reply! Thanks in advance for your kindly help!James Greennoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-22767108657448451302013-12-06T05:04:56.796-08:002013-12-06T05:04:56.796-08:00Hi Majid! Thank you so much for this information.Hi Majid! Thank you so much for this information.James Greennoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-7915958180299773502013-12-05T20:19:24.147-08:002013-12-05T20:19:24.147-08:00Sorry I'm not involved in this work any more. ...Sorry I'm not involved in this work any more. But I think found the code that made this to give you start.<br /><br /><br />----<br /><br />(set! geometry-lattice (make lattice (size 10 5 no-size)))<br /><br />(set! pml-layers (list (make pml (thickness 1))))<br /><br />(set-param! resolution 30)<br /><br />(define myAg (make dielectric (epsilon 1)<br /><br />(E-polarizations<br /><br />(make polarizability<br /><br />(omega 1e-20) (gamma 0.0038715) (sigma 4.4625e+39))<br /><br />; (make noisy-lorentzian-susceptibility<br /><br />; (noise-amp 10) (omega 1e-20) (gamma 0.0038715) (sigma 4.4625e+39))<br /><br />; (make polarizability<br /><br />; (omega 0.065815) (gamma 0.31343) (sigma 7.9247))<br /><br />; (make polarizability<br /><br />; (omega 0.36142) (gamma 0.036456) (sigma 0.50133))<br /><br />; (make polarizability<br /><br />; (omega 0.66017) (gamma 0.0052426) (sigma 0.013329))<br /><br />; (make polarizability<br /><br />; (omega 0.73259) (gamma 0.07388) (sigma 0.82655))<br /><br />; (make polarizability<br /><br />; (omega 1.6365) (gamma 0.19511) (sigma 1.1133))<br /><br />)))<br /><br />(set! geometry (list<br /><br /> (make block (center -2.5 0.5 0.0) (size 5 infinity infinity)<br /><br /> (material myAg ))))<br /><br /> <br /><br />(set! sources (list<br /><br /> ;(make source<br /><br /> ; (src (make gaussian-src (frequency 2.5) (fwidth 5)))<br /><br /> ; (component Ex)<br /><br /> ; (center 0 0 0))<br /><br /> (make source<br /><br /> (src (make gaussian-src (frequency 2.5) (fwidth 5)))<br /><br /> (component Ex)<br /><br /> (center -3 0 0))<br /><br /> ; (make source<br /><br /> ; (src (make continuous-src (frequency 2.5) ))<br /><br /> ; (component Ex)<br /><br /> ; (center -3 0 0))<br /><br /> )<br /><br /> )<br /><br /> <br /><br />;why don't i get an output when i put this ?? <br /><br />;(This must be done after specifying the<br /><br />; geometry, sources, resolution, etcetera,<br /><br />; because all of the field parameters are initialized when flux planes are created.) <br /><br />(define trans ; transmitted flux <br /><br /> (add-flux 5.0 5.0 40<br /><br /> (make flux-region (center 4 0.0 0.0) (size 0 5 0))<br /><br /> ))<br /><br /> <br /><br />(run-until 1000 (at-every 0.6 (output-png Ex "-Zc bluered")))<br /><br />; (run-until 10<br /><br /> ; (at-beginning output-epsilon)<br /><br /> ; (to-appended "ex" (at-every 0.6 output-efield-z)))<br /><br />; (run-sources+ <br /><br /> ; (stop-when-fields-decayed 50 Ex (vector3 4 0 0.0) 1e-3)<br /><br /> ; (output-png Ex "-Zc bluered")<br /><br /> ; )<br /><br /> <br /><br />(display-fluxes trans)Majid al-Dosarinoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-54662448564789278452013-10-01T11:22:52.259-07:002013-10-01T11:22:52.259-07:00I could post my code if you insist but I'm not...I could post my code if you insist but I'm not sure it would be that useful to you (my input data was not arranged as nicely with a timestamp!) Your main difficulty is in arranging the data, and for that the solution is in pandas timeseries indexing http://pandas.pydata.org/pandas-docs/stable/timeseries.html<br />In my case, once I had my hourly data, I was able to iterate by week.<br /><br /><br />As for the plot, I've used matplotlib to plot a million points without much trouble.Majid al-Dosarinoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-85439825744840740652013-10-01T03:26:32.289-07:002013-10-01T03:26:32.289-07:00Dear Majid
This sort of plot looks great! Recently...Dear Majid<br />This sort of plot looks great! Recently I began working with Python (used Matlab before) and tried to create a plot similar to this one. But I ran into problems (array size) when I tried to group my data in days and hours of the day with the lambda x-function.<br />So my question is how you managed to arrange your data, especially which type of data did you use? I used a dateframe that I fetched from a MySQL-Server. Could you give me some advise or even share your code as a iPython Notebook?<br />Thanks in advance<br />JDBJabbaDabbaDeeBeenoreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-71958820980525072342013-09-23T11:03:34.450-07:002013-09-23T11:03:34.450-07:00Sorry I'm not involved with this work anymore ...Sorry I'm not involved with this work anymore so I can't help you. However, if you leave your post here someone might help you. Thanks for your interest. Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-8332163406611119582013-09-23T05:02:20.874-07:002013-09-23T05:02:20.874-07:00Hi Majid,
I am trying to find the dispersion plot ...Hi Majid,<br />I am trying to find the dispersion plot for a air silver air interface.(IMI structure). The thickness of the metal film is 50 nm.The code is given below but i do not get expected dispersion curve(for even and odd modes).kindly tell me if my code is incorrect in logics<br /><br />;DRUDE MODEL<br /><br />(define myAg(make dielectric (epsilon 1)<br /> (polarizations (make polarizability (omega 1e-20) (gamma 0.0038715) (sigma 4.4625e+39) <br /> )<br /> )<br />))<br /><br />(define-param eps2 1) ; dielectric constant of waveguide2<br />(define-param eps3 1) ; dielectric constant of waveguide3<br />(define-param w1 .5) ; width of waveguide1<br />(define-param w2 1) ; width of waveguide2<br />(define-param w3 2) ; width of waveguide3<br /><br /><br />(define-param sy 12) ; size of cell in y direction (perpendicular to wvg.)<br />(define-param dpml 1) ; PML thickness (y direction only!)<br /><br />(set! geometry-lattice <br /> (make lattice (size 1 sy no-size))<br />)<br /><br />(set! geometry<br /> (list (make block (center 0 0) (size infinity w1 infinity) (material myAg))<br /> (make block (center 0 -0.75) (size infinity w2 infinity) (material (make dielectric (epsilon eps2))))<br /> (make block (center 0 1.25) (size infinity w3 infinity) (material (make dielectric (epsilon eps3))))<br /> <br /> )<br />)<br /><br />(set-param! resolution 20)<br /><br />(set! pml-layers (list (make pml (direction Y) (thickness dpml))))<br /><br />(define-param fcen 0.5) ; pulse center frequency <br />(define-param df 1.5) ; pulse freq. width: large df = short impulse<br /><br />(set! sources (list<br /> (make source<br /> (src (make gaussian-src (frequency fcen) (fwidth df)))<br /> (component Ez) (center 0 -0.9))))<br /><br />;(set! symmetries (list (make mirror-sym (direction Y) (phase -1))))<br /><br /><br />(define-param k-points <br /> (list <br /> (vector3 0.0 0 0) <br /> (vector3 1 0 0))) <br /><br />; define a series of k-points<br />(define-param k-interp 19) <br />(set! k-points (interpolate k-interp k-points))<br /><br />; calculates frequencies for these k-points<br />(run-k-points 300 k-points)<br /> <br />Sir<br />kindly reply as i am stuck and i will be looking forward to your reply<br /> <br /><br /><br />comet88http://www.blogger.com/profile/12060556969297970319noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-51235571121737635982012-12-18T13:35:12.494-08:002012-12-18T13:35:12.494-08:00Your concern should go to a Python MEEP expert.Your concern should go to a Python MEEP expert.Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-36571240100919853692012-12-14T03:10:01.961-08:002012-12-14T03:10:01.961-08:00i am getting segmentation error while running any ...i am getting segmentation error while running any of the example files of python meep in ubuntu 12.10 , the latest release of ubuntu. i have amd turion x2 processor , 4 gb ram. i have tried increasing the memory limits to unlimited by using the "ulimited" command . but still getting the same error. please help. My email id is ravishaw_kol@yahoo.co.in I am in dire need of python meep for my project , please help.Ravi Shawhttp://www.blogger.com/profile/07322838771453002422noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-18357468334803213742012-12-06T17:41:40.301-08:002012-12-06T17:41:40.301-08:00I'm now thinking I should have used a fast (wh...I'm now thinking I should have used a fast (which would be non-Python) regexp processor. I used Python functions that process text supposedly in C. Could an expert chime in on the fastest way?Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-37150609638618754072012-10-24T14:51:21.891-07:002012-10-24T14:51:21.891-07:00Check this post http://juluribk.com/2011/04/27/pla...Check this post http://juluribk.com/2011/04/27/plasmonic-materials-in-meep/. Part 1 only deals with conversions.Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-34621198406194123982012-10-13T01:00:14.780-07:002012-10-13T01:00:14.780-07:00Hi, I like your blog. Thanks for the explanations....Hi, I like your blog. Thanks for the explanations. one question, how do you simulate the permittivity of SiC in meep? I know one can simulate dispersive materials in meep. But just want to know the parameters.kevvu kekahttp://www.blogger.com/profile/10215841596934332795noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-8842254518102104322012-09-01T17:48:46.138-07:002012-09-01T17:48:46.138-07:00Thanks and great! I hope you continue to develop i...Thanks and great! I hope you continue to develop it.Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-47095566155272413382012-08-24T03:01:18.682-07:002012-08-24T03:01:18.682-07:00You might be interested to try my Python bindings ...You might be interested to try my Python bindings to Meep. They are object-oriented, but much less awkward to use than Python-MEEP. I reached much the same conclusion as you, that it was more trouble to use Python-MEEP than it was worth.<br /><br />My bindings are available on GitHub: https://github.com/ptomato/Pythonic-Meep<br /><br />Right now it is more of a proof of concept because only a few functions are implemented, but I have written very thorough documentation so far ;-)Philip Chimentohttp://www.blogger.com/profile/16545910914624009114noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-9173014639372160882012-06-29T00:14:13.291-07:002012-06-29T00:14:13.291-07:00i haven't checked yet but perhaps it only appe...i haven't checked yet but perhaps it only appears there are extra modes because of numerical artifacts. increasing the resolution might help. i remember it was much worse on a lower resolution.Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-22716289072608626002012-06-13T12:00:46.313-07:002012-06-13T12:00:46.313-07:00Thanks. it's unintended that i have a gap b/w ...Thanks. it's unintended that i have a gap b/w the silver and pml. i intended /that/ gap to be zero with no overlap. i thought i checked that by substituting (epsilon somenumber) and outputting the epsilon as a function of distance. i'll check again.Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-50697410087541564162012-06-13T11:46:50.204-07:002012-06-13T11:46:50.204-07:00Majid,
Excellent example! One question: it seems ...Majid,<br /><br />Excellent example! One question: it seems that the PML layer is not overlapping with the silver layer. In fact, there is an air layer between the PML and silver, causing unwanted modes to exist. Is there a reason why PML is not overlapping with Ag?Xufeng Wanghttp://www.blogger.com/profile/00622577336074268776noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-53983315593585419862012-04-30T08:31:23.686-07:002012-04-30T08:31:23.686-07:00that's what i used.that's what i used.Majid al-Dosarihttp://www.blogger.com/profile/16192275403547721607noreply@blogger.comtag:blogger.com,1999:blog-852400474287196648.post-65575314132808065692012-04-30T02:50:34.599-07:002012-04-30T02:50:34.599-07:00Hi, dear
i used your code, but it seems a lot of f...Hi, dear<br />i used your code, but it seems a lot of frequeicies is missing by harminv, i could not get the graph correctly, can you send me a full version of your codesnowlandhttp://www.blogger.com/profile/05439434645891374075noreply@blogger.com