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.

More Mohammadreza Kazemian's questions See All
Similar questions and discussions