Dear all i am tring to simulate Ag nano particles. this requires a complex eplsion value as e1+ ie2 where e1 and e2 are positive. i am putting a single nano particle in the center of the region of consideration. thecthl file is attached. after some time steps the solution just blows out and results in NaN and as a result in a black png image.the png image of this stucture shows faults along the boundary in the image as white dots. i am not able to grasp the problem. why the sub pixel averaging not working.is this the problem with ctl coding or a bug in the negative epsilon of the meep. can Dconductivity be used or i should go for polarizability.ctl program:;resolution 1nm;grid size 120 by 120;epsilon = 10.84 + 0.762i;wavelength lambda = 500nm;plane polarized monochromatic light polarised along y ( Ey ) and traveling;along xdirection(set! geometrylattice (make lattice (size 120 120 nosize)))(set! geometry (list(make cylinder (center 0 0) (radius 30) (height 110)(material (make dielectric (epsilon 10.84)(Dconductivity (/(* 2 pi 0.002 .7620) 10.84)))))))(set! sources (list (make source (src (make continuoussrc (frequency .002))) (component Ey) (direction X) (center 54 0) (amplitude 1) (size 0 100 0))))(set! Courant .5)(set! pmllayers (list (make pml (thickness 5.0))))(set! resolution 1)(rununtil 1000 (atbeginning outputepsilon) (atevery 50 outputefieldy) (atevery 50 outputefieldx))
<<attachment: newtrial-eps-000000.00.png>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

