Benchmark on neutron and charged particle production with MCNPX

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

1

projectile target

energy (MeV)

Model combination isa cem ind x x x x x x x x

p p p p

U U U U

585 660 800 1200

ina x x x x

isd

Berd

x0 x0 x0

X X X

p p p p n

W W W W W

800 1200 1800 2500 34  173 (10)

x x x x x

x

x x

x0

X X

p p p p p p n

Ta Ta Ta Ta Ta Ta Ta

585 800 1000 1200 1800 2500 34  173 (10)

x x

x x

x x

x0

X

x x x

x x x

X X X

p

Si

1200

x

x

X

p p p p p p n

Pb Pb Pb Pb Pb Pb Pb

585 600 800 1200 1800 2500 34  173 (10)

x

x

x

x x x x

x

x x

p p p p p p

Fe Fe Fe Fe Fe Fe

585 800 1200 1600 2500 2600

x x x x x x

x x x

p p p p

C C C C

585 600 660 800

x x x x

x

p p

Be Be

660 800

x

⇨Benchmark on particle production with MCNPX (neutron differential cross section, neutron multiplicity, light charged particle production etc...), for several targets of importance for radioprotection task.

● ● ●

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

x0

X X

x x x x x x

x0 x0 x0 x0 x0

X X X X X

x x

x x x x

x0 x0

X

x

x

x0

X

Eurisol October 28, 2005

●

●

MCNPX 2.5f Thin target No transport INC models : ● INCL4 ● CEM2K (CEM) ● BERTINI ● ISABEL Fission-Evaporation models : ● ABLA ● Dresner (RAL,ORNL) ● CEM2K (GEM)

2

Neutron double differential cross section ⇨neutron double differential cross section computed for different targets.

INCL4-ABLA — simulation ● data Data from : W.B. Amian et al., Nucl. Sci. Eng., vol.112, 1992, p.78

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

3

Neutron double differential cross section 5

5

θ = 30o

2

2

0.5

0.5

0.2

0.2

5

θ = 120

o

2

5

⇨Ratio between simulation and data for differential neutron cross section, plotted for several INC/Fission-Evaporation model combinations. θ = 150o

2

0.5

0.5

0.2

0.2

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

θ = 60o

— — — --—

Eurisol October 28, 2005

INCL-ABLA CEM ISABEL-ABLA ISABEL-RAL BERTINI-RAL

4

Neutron multiplicity p(1200 MeV) + U

✔ The simulated neutron multiplicity distribution (dashed line) has been folded with BNB detector efficiency (full line) to be compared with the data taken by NESSI collaboration (black points).

Data from : NIMB 170 (2000) 299 or Julich annual report 2001 (Herbach)

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

5

Light charged particle production p(x MeV) + Fe

⇨Light charged particle production cross section as a function of incident proton energy.

Data from : Final Report HINDAS p.89-90

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

6

Computing time Reaction Number of Events INCL4 ISABEL ISABEL ISABEL CEM BERTINI

ABLA RAL ABLA ORNL GEM RAL

Reaction Number of Events INCL4 ISABEL ISABEL ISABEL CEM BERTINI

p(800 MeV) + U p(800 MeV) + Pb p(800 MeV) + W 1000000 1000000 2000000 769 min 140 min 149 min 171 min 37 min 36 min

409 min 124 min 142 min

945 min 406 min 414 min

27 min 26 min

64 min 45 min

p(800 MeV) + Fe p(800 MeV) + C 2000000 1500000

ABLA RAL ABLA ORNL GEM RAL

p(800 MeV) + Be 1660000

103 min 177 min 181 min

12 min 123 min 124 min

10 min 173 min 166 min

42 min 24 min

33 min 7 min

39 min 7 min

INTEL Pentium IV 2.6 GHz

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

7

Conclusion Neutron DDXS ABLA

CEM

GEM

ISABEL

ABLA

ISABEL

RAL

BERTINI

RAL

    

    

LCP

    

Computer Time

Pr el im in ar y

INCL4

Neutron Mult.

    

⇨All models are giving reasonable results for neutron double differential cross section and neutron multiplicity, however CEM2k used with default parameters not reproduce so well the data ⇨INCL4 is much slower than other models for heavy nuclei, ISABEL requires a much greater execution time for light nuclei ⇨Perspectives :

● ● ●

To continue benchmark, especially with light charged particles Benchmark with thick targets Simulations with real geometry

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

8

Neutron DDXS INCL-ABLA

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

9

Neutron DDXS CEM

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

10

Neutron DDXS ISABEL-ABLA

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

11

Neutron DDXS ISABEL-RAL

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

12

Neutron DDXS BERTINI-RAL

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

13

Neutron DDXS θ = 30°

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

14

Neutron DDXS θ = 60°

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

15

Neutron DDXS θ = 120°

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

16

Neutron DDXS θ = 150°

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

17

Neutron multiplicity

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

18

d2σ/dEdΩ (mb/MeV/sr)

Light charged particle production p(600 MeV) + C INCL-ABLA α

45° 75° 105°

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

19

Light charged particle production p(600 MeV) + C

RAPP Benjamin CEA/DSM/DAPNIA/SPhN

Eurisol October 28, 2005

20