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>wrote: > Send meep-discuss mailing list submissions to > 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 > > You can reach the person managing the list at > 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> > To: "meep-discuss@ab-initio.mit.edu" <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> > 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 > 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 http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
