Dear Mischa, thanks for your response. I wouldn't be updating the structure every time step, more like every optical (half-)cycle. There is a little problem with the routine though. In the MEEP source code it says: // FIXME: how to handle changes in susceptibilities? which is exactly what I need to do. So I will look into the source code a bit further, but I'm not optimistic.
best, Dries ________________________________ From: Mischa Megens [mischa.meg...@gmx.com] Sent: Wednesday, July 17, 2013 11:05 PM To: Oosten, D. van (Dries) Subject: Re: [Meep-discuss] meep-discuss Digest, Vol 89, Issue 7 Hi Dries, There is fields::phase_in_material that is supposed to help you do what you want, but I have a suspicion that meep is not going to be very fast if you try changing structures every time step. But maybe taking a look at the routine points you in the right direction. Best, Mischa On 7/17/2013 1:41 PM, D.vanOosten wrote: Dear Mehdi, thanks for your response. I know the paper and the authors very well. However, they use a very different approach. They determine the field distribution using MPB and then in a postprocessing step, determine how each k vector is influence by their small perturbation. In my case, I need to calculate the effects dynamically and the effects cannot be treated perturbatively, as the material actually becomes metallic during the pulse, due to the high carrier density. Best, Dries -------- Original message -------- From: Mehdi Kamali <m.kamali4...@gmail.com><mailto:m.kamali4...@gmail.com> Date: 17/07/2013 22:21 (GMT+01:00) To: meep-discuss@ab-initio.mit.edu<mailto:meep-discuss@ab-initio.mit.edu> Subject: Re: [Meep-discuss] meep-discuss Digest, Vol 89, Issue 7 Dear Dries Look at the following article. It may be useful; PRL 108, 033902 (2012) On Wed, Jul 17, 2013 at 7:30 PM, <meep-discuss-requ...@ab-initio.mit.edu<mailto:meep-discuss-requ...@ab-initio.mit.edu>> wrote: Send meep-discuss mailing list submissions to meep-discuss@ab-initio.mit.edu<mailto:meep-discuss@ab-initio.mit.edu> To subscribe or unsubscribe via the World Wide Web, visit http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss or, via email, send a message with subject or body 'help' to meep-discuss-requ...@ab-initio.mit.edu<mailto:meep-discuss-requ...@ab-initio.mit.edu> You can reach the person managing the list at meep-discuss-ow...@ab-initio.mit.edu<mailto:meep-discuss-ow...@ab-initio.mit.edu> When replying, please edit your Subject line so it is more specific than "Re: Contents of meep-discuss digest..." Today's Topics: 1. dynamically updating the susceptibility during a MEEP simulation (Oosten, D. van (Dries)) ---------------------------------------------------------------------- Message: 1 Date: Wed, 17 Jul 2013 07:26:35 +0000 From: "Oosten, D. van (Dries)" <d.vanoos...@uu.nl<mailto:d.vanoos...@uu.nl>> To: "meep-discuss@ab-initio.mit.edu<mailto:meep-discuss@ab-initio.mit.edu>" <meep-discuss@ab-initio.mit.edu<mailto:meep-discuss@ab-initio.mit.edu>> Subject: [Meep-discuss] dynamically updating the susceptibility during a MEEP simulation Message-ID: <fa281bb491411842a76ebd8ca731d3a4230ac...@ictsc-w-s202.soliscom.uu.nl<mailto:fa281bb491411842a76ebd8ca731d3a4230ac...@ictsc-w-s202.soliscom.uu.nl>> Content-Type: text/plain; charset="us-ascii" Dear fellow MEEP users, I have been struggling for quite some time with this problem. I run a MEEP simulation in which a femtosecond laser pulse impinges on a structure containing a semiconductor material (say silicon). In our experiment, we use such high light intensities that during the pulse, an extremely high carrier density is generated in the semiconductor. A this happens during the pulse, we have to take the changing carrier density, and therefore the changing plasma frequency, into account. We can do this in our own FDTD code relatively easily, but we would like to make use of MEEPs readily available MPI support. So the question is, is there a way to change the susceptiblity of the structure during a simulation. What we have in mind is the following. We will step() MEEP forward for half an optical cycle, calculate the intensity distribution inside the semiconductor and use it to update the carrier density (and perhaps carrier temperature). From this we will calculate new parameters for the Drude susceptiblity of the semiconductor. We would then need a way to update these. Simply writing them back into the array using C does not work. Upon inspection of the source code of MEEP I understand that this is because the array you specify is actually only used as a input to generate the chunks. So I would need a way to trigger the distribution of the new parameters to the chunks. My understanding of the code is however by far not good enough to know how to achieve that. Is there a way of doing this without changing the source code of MEEP? Some undocumented function that when used wisely can achieve this? I would be very thankful for any advice! kind regards, Dries ------------------------------ _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu<mailto:meep-discuss@ab-initio.mit.edu> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss End of meep-discuss Digest, Vol 89, Issue 7 ******************************************* _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu<mailto:meep-discuss@ab-initio.mit.edu> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
