hi. I'm trying to simulate a pretty simple geometry including nanoparticles by mcnp. my whole cod is as below:1
project nano
C Cell cards
1 2 -1 -7 +8 +9 -10 -11 +12 fill=1
2 1 -0.00125 -1 +2 -3 +4 -5 +6 #1
3 2 -1 14 -15 16 -17 fill=2 lat=1 u=1
4 3 -19.3 -13 u=2
5 2 -1 13 u=2
6 0 1:-2:3:-4:+5:-6
C Surface cards
1 px 5
2 px -5
3 py 5
4 py -5
5 pz 5
6 pz -5
7 px 0.5
8 px -0.5
9 py 1.5
10 py 2.5
11 pz 0.5
12 pz -0.5
13 sy 2 7E-6
14 pz -7E-6
15 pz 7E-6
16 px -7E-6
17 px 7E-6
C DATA cards
mode e
IMP:e 1 1 1 1 1 0
SDEF pos=0.1 0.1 0.1 ERG=d1 PAR=3 CELL=2
SI1 L 0.36 0.71 0.81
SP1 0.3 0.49 0.2
M1 6012 -0.000150
7014 -0.784431
8016 -0.210748
18000 -0.004671
M2 1001 2
8016 1
M3 79197 1
*f8:e 1
NPS 1000000
C END of data cards
it is about a tumor ( cell 1) that is filled by gold nano-particles. nano-particles have a sphere shape with a radius equal to 70 nanometers. in order to fill the tumor by nano-particles I created a cube with each side equal to 150 nanometers ( cell 3 ) and I filled it by nano-spheres ( cell 4 & 5). And at the end cell, 1 ( the tumor ) is a lattice of cell 3 ( which includes nano-spheres ).
the problem is when I run the code, dose tally at cell 1 is always zero! and "track entering" and "population" in every cell except 2 and 5 are zero too! In another word, no particle enters cells 1 & 3 and 4! specifically, cell 1 is more important because I want to calculate dose in tumor. I can't see where is my mistake, so please let me know if you noticed anything.