Hi meep users

A warm regard to everybody.
I am trying to study the absorption enhancement of a 2-D thin Si substrate 
using periodic Ag nanoparticles at different periods (0.1- 0.8). Absorption 
peak goes on increases and its width goes on decreases with increasing period. 
Can anybody tell me  where the problem is. I am waiting for any suggestions. My 
ctl is ;


;Characteristic length, a = 1 microm 



(set! geometry-lattice (make lattice (size 0.1 1 no-size)))

(set! geometry (list
       
        (make block (center 0 -0.05) (size infinity 0.1 no-size)
                (material (make dielectric (epsilon 10.1)
        (E-polarizations
        (make polarizability 
        (omega 1.659)(gamma 0.1946)(sigma 0.2445))
        (make polarizability 
        (omega 1.981)(gamma 2.893e-8)(sigma 4.001))))))

        
        (make block (center 0 0.03)(size 0.06 0.06 no-size)
                (material (make dielectric (epsilon 1)
        (E-polarizations

(make polarizability
(omega 1e-20) (gamma 0.038715) (sigma 4.4625e+41))
(make polarizability
(omega 0.65815) (gamma 3.1343) (sigma 7.9247))
(make polarizability
(omega 3.6142) (gamma 0.36456) (sigma 0.50133))
(make polarizability
(omega 6.6017) (gamma 0.052426) (sigma 0.013329))
(make polarizability
(omega 7.3259) (gamma 0.7388) (sigma 0.82655))
(make polarizability
(omega 16.365) (gamma 1.9511) (sigma 1.1133))

))))))


(set! sources (list
               (make source
                 (src (make gaussian-src (frequency 1.21) (fwidth 2.4)))
                 (component Ex)
                 (center 0 0.4) (size 0.1 0))))


(set! k-point (vector3 0 0 0))
(set! ensure-periodicity true)
(set! symmetries (list (make mirror-sym (direction X) (phase -1))))

(set! pml-layers (list (make pml (direction Y) (thickness 0.2))))

(set! resolution 200)

(define planeY1                                  
     (add-flux 1.21 0.88 100
                (make flux-region
                   (center 0 0) (size 0.1 0)(weight -1)(direction Y))
        (make flux-region
                   (center 0 -0.1) (size 0.1 0)(direction Y)))) 

(run-until 200
          (at-beginning output-epsilon)
    (to-appended "ex" (at-every 0.5 output-efield-x))    
    (at-end output-efield-x))

(display-fluxes planeY1)


_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to