Good day everyone, I am simulating scatering properties of some particle. For a size smaller than 35nm calculations run very well and I obtain good results. However when I increase the size of Nanoparticle to more than 40nm, at the end of calculations I obtain '+nan of -nan' as result (see this result below). I really do not know what can be the matter. Has anyone experienced this kind of problem? if it is the case, How can I handle the problem?
King Regards, Maximin. ###############results############################### Computational cell is 6.2 x 6.2 x 6.2 with resolution 20 block, center = (0,0,0) size (3.2,3.2,3.2) axes (1,0,0), (0,1,0), (0,0,1) dielectric constant epsilon diagonal = (1,1,1) sphere, center = (0,0,0) radius 0.4 dielectric constant epsilon diagonal = (1.00001,1.00001,1.00001) time for set_epsilon = 0.00224357 s lorentzian susceptibility: frequency=5.7743, gamma=168.685 lorentzian susceptibility: frequency=0.27441, gamma=0.0145899 lorentzian susceptibility: frequency=0.244131, gamma=0.0761473 lorentzian susceptibility: frequency=0.292098, gamma=0.0257822 lorentzian susceptibility: frequency=0.338166, gamma=1e-20 lorentzian susceptibility: frequency=1e-20, gamma=0.555016 ----------- on time step 1 (time=0.025), 108.367 s/step field decay(t = 2.025): 2.54885897962345e-50 / 2.54885897962345e-50 = 1.0 field decay(t = 3.05): 7.6748148994115e-15 / 7.6748148994115e-15 = 1.0 field decay(t = 4.05): 1.31630221258758e-11 / 1.31630221258758e-11 = 1.0 on time step 162 (time=4.05), 0.0248906 s/step field decay(t = 5.075): 9.44556429017964e-11 / 9.44556429017964e-11 = 1.0 field decay(t = 6.1): 1.62836044305629e-9 / 1.62836044305629e-9 = 1.0 . . . field decay(t = 28.625): 0.138646609982224 / 0.275371699173149 = 0.50348895837348 field decay(t = 29.65): 0.213473641248043 / 0.275371699173149 = 0.775219973181829 field decay(t = 30.675): 0.131922108486818 / 0.275371699173149 = 0.479069232179402 field decay(t = 31.7): 0.0954612933808837 / 0.275371699173149 = 0.346663414096374 on time step 1307 (time=32.675), 0.0243314 s/step field decay(t = 32.725): 0.0940538122071552 / 0.275371699173149 = 0.341552209212377 field decay(t = 33.75): 0.0196737053708999 / 0.275371699173149 = 0.0714441804657978 field decay(t = 34.775): 0.025208020379212 / 0.275371699173149 = 0.0915417977043517 field decay(t = 35.8): 0.0114590958545844 / 0.275371699173149 = 0.0416131936905366 on time step 1472 (time=36.8), 0.0243567 s/step field decay(t = 36.825): 0.00415199935839779 / 0.275371699173149 = 0.0150777998278868 field decay(t = 37.85): 0.00349274026407682 / 0.275371699173149 = 0.0126837299350818 field decay(t = 38.875): 5.05190094437579e-4 / 0.275371699173149 = 0.00183457521580649 field decay(t = 39.9): 4.95306600649438e-4 / 0.275371699173149 = 0.00179868375049681 on time step 1636 (time=40.9), 0.0244389 s/step field decay(t = 40.925): 1.02314595843005e-4 / 0.275371699173149 = 3.7155087523599e-4 field decay(t = 41.95): 9.65202405693598e-5 / 0.275371699173149 = 3.50508933413195e-4 field decay(t = 42.975): 8.15482036517457e-5 / 0.275371699173149 = 2.96138651490361e-4 field decay(t = 44.0): 6.70990812661126e-5 / 0.275371699173149 = 2.43667310284931e-4 on time step 1797 (time=44.925), 0.024944 s/step field decay(t = 45.025): 8.52356861711356e-5 / 0.275371699173149 = 3.09529579209012e-4 field decay(t = 46.05): 7.0618746192027e-5 / 0.275371699173149 = 2.56448815924338e-4 field decay(t = 47.075): 7.911590079658e-5 / 0.275371699173149 = 2.87305852540908e-4 field decay(t = 48.1): 7.03635288690586e-5 / 0.275371699173149 = 2.55522005639422e-4 . . . on time step 7472 (time=186.8), 0.0249913 s/step field decay(t = 187.5): 2.10763513530665e-7 / 0.275371699173149 = 7.65378265680601e-7 field decay(t = 188.525): 2.60825953182775e-7 / 0.275371699173149 = 9.47177774498794e-7 field decay(t = 189.55): 1.69770375757494e-7 / 0.275371699173149 = 6.16513520696785e-7 field decay(t = 190.575): 2.84247479109823e-7 / 0.275371699173149 = 1.03223199756302e-6 on time step 7636 (time=190.9), 0.0245153 s/step field decay(t = 191.6): 2.8709873678423e-7 / 0.275371699173149 = 1.04258621218627e-6 field decay(t = 192.625): 2.81378832234225e-7 / 0.275371699173149 = 1.02181463483398e-6 field decay(t = 193.65): 2.63798543898019e-7 / 0.275371699173149 = 9.57972604629015e-7 field decay(t = 194.675): 3.30769015732812e-7 / 0.275371699173149 = 1.20117287552063e-6 on time step 7796 (time=194.9), 0.0251093 s/step field decay(t = 195.7): 2.57771844758281e-7 / 0.275371699173149 = 9.3608691645614e-7 field decay(t = 196.725): 2.95754091439449e-7 / 0.275371699173149 = 1.07401774520585e-6 field decay(t = 197.75): 3.52576435602714e-7 / 0.275371699173149 = 1.28036554468518e-6 field decay(t = 198.775): 4.33952266121196e-7 / 0.275371699173149 = 1.57587823085746e-6 on time step 7957 (time=198.925), 0.0248622 s/step field decay(t = 199.8): 2.56442035606597e-7 / 0.275371699173149 = 9.31257774043621e-7 . . . field decay(t = 219.275): 0.0723841190338135 / 0.275371699173149 = 0.262859688381773 on time step 8777 (time=219.425), 0.0244014 s/step field decay(t = 220.3): 0.0600824356079102 / 0.275371699173149 = 0.218186675639937 field decay(t = 221.325): 0.306595802307129 / 0.306595802307129 = 1.0 field decay(t = 222.35): 1.60675430297852 / 1.60675430297852 = 1.0 field decay(t = 223.375): 1.18775939941406 / 1.60675430297852 = 0.73922901417613 on time step 8941 (time=223.525), 0.0245355 s/step field decay(t = 224.4): 1.40089416503906 / 1.60675430297852 = 0.871878271893941 field decay(t = 225.425): 4.59901428222656 / 4.59901428222656 = 1.0 field decay(t = 226.45): 6.646728515625 / 6.646728515625 = 1.0 field decay(t = 227.475): 22.711181640625 / 22.711181640625 = 1.0 on time step 9104 (time=227.6), 0.0245667 s/step field decay(t = 228.5): 51.885009765625 / 51.885009765625 = 1.0 field decay(t = 229.525): 69.09765625 / 69.09765625 = 1.0 field decay(t = 230.55): 210.25 / 210.25 = 1.0 field decay(t = 231.575): 1489.4775390625 / 1489.4775390625 = 1.0 on time step 9265 (time=231.625), 0.0248858 s/step field decay(t = 232.6): 1168.78515625 / 1489.4775390625 = 0.784694717172876 field decay(t = 233.625): 2268.140625 / 2268.140625 = 1.0 field decay(t = 234.65): 3555.140625 / 3555.140625 = 1.0 field decay(t = 235.675): 13514.0625 / 13514.0625 = 1.0 on time step 9428 (time=235.7), 0.0245732 s/step field decay(t = 236.7): 43472.25 / 43472.25 = 1.0 . . . field decay(t = 1156.1): 1.447139343293e304 / 1.62976459153031e304 = 0.887943786982249 field decay(t = 1157.125): 5.00945341949369e304 / 5.00945341949369e304 = 1.0 field decay(t = 1158.15): 9.41755611784805e304 / 9.41755611784805e304 = 1.0 on time step 46365 (time=1159.12), 0.0256248 s/step field decay(t = 1159.175): 8.11024864779285e304 / 9.41755611784805e304 = 0.861184 field decay(t = 1160.2): 3.99126241090812e305 / 3.99126241090812e305 = 1.0 field decay(t = 1161.225): 1.15101521553237e306 / 1.15101521553237e306 = 1.0 field decay(t = 1162.25): 2.96460715401868e306 / 2.96460715401868e306 = 1.0 on time step 46526 (time=1163.15), 0.024989 s/step field decay(t = 1163.275): 4.52712386536198e306 / 4.52712386536198e306 = 1.0 field decay(t = 1164.3): 9.3874440472146e306 / 9.3874440472146e306 = 1.0 field decay(t = 1165.325): 1.16787580528217e307 / 1.16787580528217e307 = 1.0 field decay(t = 1166.35): 6.71837611119561e307 / 6.71837611119561e307 = 1.0 on time step 46688 (time=1167.2), 0.0248351 s/step field decay(t = 1167.375): 1.42595336796086e308 / 1.42595336796086e308 = 1.0 field decay(t = 1168.4): +inf.0 / +inf.0 = +nan.0 run 0 finished at t = 1168.4 (46736 timesteps) flux1:, 0.08032, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.0834185507246377, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.0865171014492754, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.089615652173913, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.0927142028985507, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.0958127536231884, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.0989113043478261, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0 flux1:, 0.102009855072464, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0, +nan.0
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
