Meep Tutorial/Material dispersion
From AbInitio
(diff) ←Older revision | Current revision | Newer revision→ (diff)
Meep |
Download |
Release notes |
FAQ |
Meep manual |
Introduction |
Installation |
Tutorial |
Reference |
C++ Tutorial |
C++ Reference |
Acknowledgements |
License and Copyright |
In this example, we will perform a simulation with a frequency-dependent dielectric ε(ω), corresponding to material dispersion. (See Dielectric materials in Meep for more information on how material dispersion is supported in Meep.) In particular, we will model a uniform medium of the dispersive material; see also the material-dispersion.ctl
file included with Meep. From the dispersion relation ω(k), we will compute the numerical ε(ω) via the formula:
We will then compare this with the analytical ε(ω) that we specified.
Since this is a uniform medium, our computational cell can actually be of zero size (i.e. one pixel), where we will use Bloch-periodic boundary conditions to specify the wavevector k.
(set! geometry-lattice (make lattice (size no-size no-size no-size)))
We will then fill all space with a dispersive material:
(set! default-material (make dielectric (epsilon 2.25) (polarizations (make polarizability (omega 1.1) (gamma 1e-5) (delta-epsilon 0.5)) (make polarizability (omega 0.5) (gamma 0.1) (delta-epsilon 2e-5)) )))
corresponding to the dielectric function:
To be continued...