Hi Filip,
I can't wait to read your post. please give us the link when you finish it.
Best regard,
P. Tapsanit
Department of Physics,
Tohoku unibersity
________________________________
From: "[email protected]"
<[email protected]>
To: [email protected]
Sent: Thursday, May 16, 2013 1:00 AM
Subject: meep-discuss Digest, Vol 87, Issue 10
Send meep-discuss mailing list submissions to
[email protected]
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
[email protected]
You can reach the person managing the list at
[email protected]
When replying, please edit your Subject line so it is more specific
than "Re: Contents of meep-discuss digest..."
Today's Topics:
1. Calculating complex epsilon in MEEP 1.2 (Kahli Burke)
2. Re: Calculating complex epsilon in MEEP 1.2 (Filip Dominec)
----------------------------------------------------------------------
Message: 1
Date: Tue, 14 May 2013 23:48:08 -0700
From: Kahli Burke <[email protected]>
To: "[email protected]" <[email protected]>
Subject: [Meep-discuss] Calculating complex epsilon in MEEP 1.2
Message-ID: <[email protected]>
Content-Type: text/plain; charset=us-ascii
Hi,
I was looking at some metal material properties created by Aaron Webster
(http://falsecolour.com/aw/meep_metals/meep-metals.pdf), he referenced a way to
have MEEP show the complex epsilon for the material using a 1x1x1 pixel grid
and the 'meep-fields-analytic-chi1' function. This would let you specify a
frequency for the field and very quickly print out the complex dielectric at
any frequency.
This function doesn't exist anymore, it was in 1.1.1 but not in 1.2. I'm
curious if there is a similar low level function that can do something similar,
I'd like to check that these material properties are still lining up with other
analytic models like the Drude and Debye models.
It looks like the polarization code was rewritten for 1.2, perhaps to support
more complex materials like anisotropic dielectrics. I looked through the SWIG
wrapper code and there are some functions that reference 'chi1' (e.g.
'meep-material-function-eff-chi1inv-row') but I'm not sure how to use them or
if they would be of any use for my goal. Any insight would be appreciated.
Thanks,
Kahli Burke
University of Oregon
------------------------------
Message: 2
Date: Wed, 15 May 2013 10:33:00 +0200
From: Filip Dominec <[email protected]>
To: Kahli Burke <[email protected]>
Cc: "[email protected]" <[email protected]>
Subject: Re: [Meep-discuss] Calculating complex epsilon in MEEP 1.2
Message-ID:
<calfwsr5yrth+u2by38nbnfpjnywnbav+0oxqjkarxhxa8vt...@mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1
Hi, Kahli,
few months ago, I needed to run some time-domain and frequency-domain
simulations on the same structure, which contained dispersive
materials. As MEEP does not accept dispersive model in
frequency-domain computations yet, I had to compute the value of
complex epsilon for each frequency manually.
The formula turned out to be the simplest one, known from textbooks,
something like:
eps_complex(omega) = eps_infinity + sigma/(omega_t**2- i*gamma*omega -
omega) + ...
... where omega_t, gamma, and sigma are the three parameters that are
passed to MEEP to define the resonance. You may add as many lorentzian
terms as needed.
I verified that the points (of s-parameters calculation) obtained from
the frequency-domain simulation match the time-domain curve very well.
I am going to put the whole procedure on the web soon. The conclusion
is that you may get along without the meep-fields-analytic-chi1
function and use the Lorentz formula to get correct results.
Filip
2013/5/15, Kahli Burke <[email protected]>:
> Hi,
>
> I was looking at some metal material properties created by Aaron Webster
> (http://falsecolour.com/aw/meep_metals/meep-metals.pdf), he referenced a way
> to have MEEP show the complex epsilon for the material using a 1x1x1 pixel
> grid and the 'meep-fields-analytic-chi1' function. This would let you
> specify a frequency for the field and very quickly print out the complex
> dielectric at any frequency.
>
> This function doesn't exist anymore, it was in 1.1.1 but not in 1.2. I'm
> curious if there is a similar low level function that can do something
> similar, I'd like to check that these material properties are still lining
> up with other analytic models like the Drude and Debye models.
>
> It looks like the polarization code was rewritten for 1.2, perhaps to
> support more complex materials like anisotropic dielectrics. I looked
> through the SWIG wrapper code and there are some functions that reference
> 'chi1' (e.g. 'meep-material-function-eff-chi1inv-row') but I'm not sure how
> to use them or if they would be of any use for my goal. Any insight would be
> appreciated.
>
> Thanks,
>
> Kahli Burke
> University of Oregon
> _______________________________________________
> meep-discuss mailing list
> [email protected]
> http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
>
------------------------------
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
End of meep-discuss Digest, Vol 87, Issue 10
********************************************
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
