HiI have three questions. I get the result in dielectric structure. But I do not get the same result in metal-dielectric structure compare to the text book. My questions are
First, how to get the 2D metal-dielectric w-k diagram.
Second, 3D dielectric-metal-dielectric slab structure w-k diagram.
Third, U-shape Transmission spectrum. (negative index material)
I think that First case is very simple problem but very important. My result (w-k diagram) is a little same compared to the text book, but do not exactly match.My idea is
First 2D metal-dielectric w-k diagram,
x-direction-size=1*a(lattice constant) , boundary condtion :bloch
y-direction-size=10*a(lattice constant), boundary conditoin : PML
upper half side(y>0) region: metal, wp= 1.0e16, gamma=0.0 therefore delta-epsilon=1.0e40
lower half side(y<0) region: dielectric epsilon: ed=1.0(air)
Source: Gaussian pulse, component: Hz.
In this case, wsp must content to be this value wsp=wp/(1+ed)^0.5=1/(1+1)^0.5= 0.707. But at w-k diagram wsp is about 0.9. When I simulate in bulk metal, wp is exactly 1. I think that meep do not supported in metal-dielectric case. Or I do not know how to make the ctl file.
Second, 3D dielectric-metal-dielectric slab structure w-k diagram.
x-direction-size=1, boundary condition: bloch
y-direction-size=1, boundary condition: bloch
z-direction-size=5, boundary condition: PML
Thinkness 0.3 metal is implemented at the center. background epsilon=1.0 (air)
I think that if First problem is solved, Secon problem is sloved, too. But if First problem is not solved, Second problem is not solved, either.
Third, U-shape Transmission spectrum. (negative index material, metamaterials)
SCIENCE vol 36 19 November 2004 Magnetic response of metamaterials at 100 Terahrtz. Figure 3. C, D, G, H. When I simulate the structure, do not appear the resonance peak.
lattice constant a=100nm
x-direction-size: 7*a=700nm, boundary condition: PML, size: 0.5
y-direction-size: 7*a=700nm, boundary condtion: PML, size: 0.5
z-direction-size: 4*a=400nm, boundary condition: PML, size: 0.5
wp=2pi2.175*1.65e15, gamma=2pi6.5*1.65e12
delta-epsilon=(wp/wn/omega)^2=(0.725 /1e-20)^2=5.25625e39
gamma=0.002167.
Thanks.First.(set! geometry-lattice (make lattice (size 1 10 no-size)))
(set! eps-averaging? false)
(set! geometry
(list
(make block (center 0 2.5) (size 1 5)
(material (make dielectric (epsilon 1.0))))
(make block (center 0 - 2.5) (size 1 5)
(material (make dielectric (epsilon 1.0)
(polarizations
(make polarizability (omega 1.0e-20) (gamma 0.0) (delta-epsilon 1.0e+40))))))))(set-param! resolution 50)
(set! pml-layers (list (make pml (direction Y) (thickness 0.5))))
(define-param k-interp 15)
(define-param fcen 0.5)
(define-param df 1.0)
(set! sources
(list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df) ))
(component Hz) (center 0 0))))(run-k-points 300 (interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0)) ))
Second
(set! geometry-lattice (make lattice (size 1 1 5)))
(set! eps-averaging? false)
(set! geometry
(list
(make block (center 0 0 0) (size 1 1 5)
(material (make dielectric (epsilon 1.0))))
(make block (center 0 0 0) (size 1 1 0.3)
(material (make dielectric (epsilon 1.0)
(polarizations (make polarizability (omega 1.0e-20) (gamma 0.0) (delta-epsilon 1.0e+40 ))))))))(set! pml-layers (list (make pml (direction Z) (thickness 0.5))))
(set-param! resolution 50)
(define-param k-interp 15)
(define-param fcen 0.5)
(define-param df 1.0)(set! sources
(list (make source (src (make gaussian-src (frequency fcen) (fwidth df)))
(component Hx) (center 0 0 0) )))(run-k-points 300
(interpolate k-interp (list (vector3 0 0 0) (vector3 0.5 0 0))))Third
(set! geometry-lattice (make lattice (size 7 7 4)))
(set! eps-averaging? false)
(define-param a_value 100e-9)(define wn (/ (* 2 pi 2.99792458e8) a_value))
(define-param e0 1)
(define-param g1 (* 2 pi 6.5 1.65e12))
(define-param wp (* 2 pi 2.175 1.65e15))(define material1
(make dielectric (epsilon e0)
(polarizations
(make polarizability
(omega 1e-20) (gamma (/ g1 wn)) (delta-epsilon (/ (* wp wp) (* wn wn 1e-20 1e-20)))
))))(set! geometry
(list
(make block (center -0.975 -1.15 0) (size 1.25 0.9 0.2)
(material material1))
(make block (center 0.975 -1.15 0) (size 1.25 0.9 0.2)
(material material1))
(make block (center - 1.15 0 0) (size 0.9 1.4 0.2)
(material material1))
(make block (center 1.15 0 0) (size 0.9 1.4 0.2)
(material material1))
(make block (center 0 1.15 0) (size 3.2 0.9 0.2)
(material material1))
(make block (center 0 0 0.125) (size 7 7 0.05)
(material (make dielectric (epsilon 3.8))))))
(set! pml-layers (list (make pml (thickness 0.5))))
(set-param! resolution 20)
(define-param fcen 0.5)
(define-param df 0.9)
(define-param nfreq 500) ; number of frequencies at which to compute flux(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey)
(center 0 0 -0.7)
(size 6 6 0)
)))(define trans ; transmitted flux
(add-flux fcen df nfreq
(make flux-region
(center 0 0 0.7) (size 6 6 0) (direction Z)
)))
(run-sources+ 500
(at-beginning output-epsilon))
(display-fluxes trans)
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
