Hello! I am trying to calculate field energy in air before the field passes through a medium and after it passes through. The medium has a refractive index of 2.1609.
I have used two cases: with periodic boundary conditions (using k-point) and without it. Also I have put metal nano cylinder and this leads us to a total of four cases as given below. nc- nano cylinder kp -k-point One - field -energy after the medium two - field - energy before the medium the field energy computations were done at the last time step (after third time step in this case) 1. *no nc and yes kp* : One-3.34630186473115e-5 two-3.5681340908767e-5 2. *yes nc and no kp* : One-6.37494514942021e-8 two-3.80564486417271e-5 3. *yes nc and yes kp* : One-1.75693509082845e-6 two-6.01227187135615e-5 4. *no nc and no kp* : One-3.05763975245008e-5 two-2.52388151348767e-5 Question: How are the results getting affected by k-pojnt condition? Especially if we compare 'yes nc' case then there is a difference of 3 orders of magnitude!! My question is that is my k-point condition correct in the program? ----------------------------------------------------------------------------------------------------------- CTL FILE: (set! geometry-lattice (make lattice (size .5 .09 no-size))) (set! geometry (list (make block (center 0 0) (size .3 .09 0.4) (material(make dielectric (epsilon 2.1609)))))) ;(make cylinder (center 0 0) (radius .035) (height .110) ;(material(make dielectric (epsilon 2.16091) ; (polarizations ;(make polarizability ;(omega 1e-30) (gamma .1561792) ;(delta-epsilon 41.9731539e+60)))))))) (set! k-point (vector3 1 1 0)) (set! sources (list (make source (src (make continuous-src (frequency 2.1777))) (direction X) (component Hz) (center -.2 0) (size 0 .09 0)))) (set! pml-layers (list (make pml (thickness 0.025)(direction X)))) (set! resolution 1000) (run-until 2 (at-beginning output-epsilon)) (run-until 1 (at-every .5 output-efield-x) (at-every .5 output-efield-y)) (print "One" (field-energy-in-box (volume (center 0.15 0) (size 0.01 0.01 0))) " " "two" (field-energy-in-box (volume (center -0.15 0) (size 0.01 0.01 0)))) ---------------------------------------------------------------------------------------------------------------- Thanks! regards -- Ajay Ram S K S M MSc(Physics) Sri Sathya Sai University Prasanthi Nilayam India.
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss