The Status of Nuclear Data above 20 MeV Masayoshi SUGIMOTO, Tokio FUKAHORI Japan Atomic Energy Agency IAEA’s Technical Meeting on Nuclear Data Libraries.

The Status of Nuclear Data above 20 MeV

Masayoshi SUGIMOTO, Tokio FUKAHORIJapan Atomic Energy Agency

IAEA’s Technical Meeting on Nuclear Data Librariesfor Advanced Systems: Fusion Devices (NuDL:FD)

31 October – 2 November 2007, Vienna1

1st priority(39)

1H, 12C, 14N, 16O, 27Al, 50,52,53,54Cr, 54,56,57,58Fe, 58,60,61,62,64Ni, 63,65Cu, 180,182,183,184,186W, 196,198,199,200,201,202,204Hg, 204,206,207,208Pb, 209Bi, 235,238U

2nd priority(43)

9Be, 10,11B, 24,25,26Mg, 28,29,30Si, 39,41K, 40,42,43,44,46,48Ca, 46,47,48,49,50Ti, 51V, 55Mn, 59Co, 90,91,92,94,96Zr, 93Nb, 92,94,95,96,97,98,100Mo, 238,239,240,241,242Pu

3rd priority(40)

2H, 6,7Li, 13C, 19F, 23Na, 35,37Cl, 35,38,40Ar, 64,66,67,68,70Zn, 69,71Ga, 70,72,73,74,76Ge, 75As, 89Y, 181Ta, 197Au, 232Th, 233,234,236U, 237Np, 241,242,242m,243Am, 243,244,245,246Cm

4th priority(10)

15N, 18O, 74,76,77,78,80,82Se, 113,115In

Nuclides with red color (66) : Released in March 2004 as JENDL/HE-2004Nuclides with underline are revised for JENDL/HE-2007.

Nuclides with blue color (36): Additionally Release in 2007 as JENDL/HE-2007

JAPAN ATOMIC ENERGY AGENCY

Target Nuclides and Priorites

Neutron & Proton File up to 3 GeV (Total: 132 nuclides)

2

E < 250 MeV

E > 150 MeV

ECISOPTMAN

OM Cal.

JQMDJAM

QMD+SDM

INC+GEM

TOTELA

Systematics

FISCAL

Systematics

GNASH

Pre-equilibriumStatistical Cal.

Isotope Prod. X-sec.Light Particle DDX

Total, ElasticReaction X-sec.

Total, ElasticReaction X-sec.

Transmission C.D.Inela. X-sec.

Isotope Prod. X-sec.Light Particle Spectra

Fraction of MSD

Fission X-sec.

JENDL/HE


Evaluation Methods and Tools

3

No global and wide energy range phenomenological potential parameter set previously

OMP for Cluster Particles

Kalbach’s model made significant underestimation for cluster particle emission spectra.

Cluster Particle Emission

Simple folding potential parameter setsnew

Iwamoto-Harada-Sato Preequilibrium Model (< 100 MeV)

new


Newly Adopted are …

Above are applied for 56Fe, Zr, Nb, W, Pb and Bi.4

Based on methods proposed by Madland

Simple Folding Potential

target

U(E)

Alpha-Particle Potentiale.g. )e.g. )

neutron

proton

=Un(E/4)

Up(E/4)

target

target

neutron × ２

proton × ２

Kunieda-Chiba-Shibata (2007)

Potential for neutron and proton

1 keV-200 MeV, 26 < A < 238


OMP for Cluster Particles

5http://www.jstage.jst.go.jp/article/jnst/44/6/44_838/_article

56Fe 90Zr 120Sn 208Pb

Angle (deg.)

/ r

th.

DeuteronDeuteron Elastic Scattering Elastic ScatteringExp. data PresentBojowald

(`88)

6

58Ni 89Y 208Pb 232Th

Angle (deg.)

/ r

th.

Exp. data PresentBecchetti (`69)

TritonTriton Elastic Scattering Elastic Scattering

7

58Ni 90Zr 120Sn 208Pb

Angle (deg.)

/ r

th.

He-3He-3 Elastic Scattering Elastic ScatteringExp. data Present

8

Total Absorption Cross Section

3He + 118Sn PresentBecchetti

Not Good

Exp.

9

64Zn

92Mo 107Ag

208Pb

Angle (deg.)

/ r

th.

Exp. data PresentAvrigeanu (`94)

AlphaAlpha Elastic Scattering Elastic Scattering

10

Preequilibrium Model for Cluster ParticlesPreequilibrium Model for Cluster Particles

( )( , ) 2 3

(. ,)2 1 ( , , )

( )(

( ),

( ), )

xin

xn l lv mm x

lF

s p h UW

p h E

Emission Rate of Cluster Particles

Formation Factors

Iwamoto-Harada-Sato Model

2p-1h State

Ef

Inverse X-sec（ e.g. Target ）Optical Model

11

Fe(n,xFe(n,x)) Spectrum Spectrum

ENDF/B-VIIJEFF-3.1This work

Sylpen+ (2004)

EEnn= 62.7 MeV = 62.7 MeV

12

Pb(n,xPb(n,x)) Spectrum Spectrum

ENDF/B-VIIJEFF-3.1This work

Blideanu+ (2004)

EEnn= 96 MeV = 96 MeV

13

JENDL-HE2004ENDF/B-VI

JEFF-3.1

Fe(n,xFe(n,x) ) Cross SectionCross Section

Baba+ (1994)

Haight+ (2004)

Incident Neutron Energy (MeV)

(b

arn

)

Present

ENDF/B-VII adopts exp. data. 14

Sample of Evaluated Results

101 102 1030

0.05

0.1

0.15

Incident Energy [MeV]

Cro

ss S

ecti

on [

b]

12C(p,np)11C Hintz, et al. '52 Aamodt, et al. '52 Goebel, et al. '61 Cumming '63 Crandall, et al. '56 Andrews, et al. '68 Measday '66 Parikh '60 Cumming, et al. '58 Horwitz, et al. '60 Kavanagh, et al. '64 Kitwanga, et al. '89 Poskanzer, et al. '64 LA150 JENDL/HE-2004Isotope production Isotope production

cross sections : cross sections : 1212CC

101 102 1030

0.01

0.02

0.03

0.04

0.05


Cro

ss S

ecti

on [

b] exp. JENDL/HE LA150

12C(n,2n)11C

15


101 102 1030

0.01

0.02

0.03

0.04

0.05


Cro

ss S

ecti

on [

b]

12C(p,x)7Be Lange, et al. '95 Bodeman, et al. '93 Lupke, et al. '93 Bodeman, et al. '91 Cumming '63 Willams, et al. '67 Honda, et al. '60 Baker et al. '58 Marquez, et al. '51 Michel, et al. '97 Michel, et al. '95 Bodemann, et al. '93 Schiekel, etal. '96 Fassbender, et al. '98 LA150 JENDL/HE-2004

101 102 1030

0.005

0.01

0.015


Cro

ss S

ecti

on [

b]

27Al(p,x)7Be Grutter '82 Aleksandrov+ '88 Michel+ '86 Dittrich+ '90 Bodemann+ '93 Lupke+ '93 Williams+ '67 Sisterson+ '96 Baker+ '58 Neuzil+ '63 Lafreur+ '66 Taddeucci+ '97 Gruetter '82 Michel+ '97 Michel+ '95 Bodemann+ '93 Schiekel+ '96 Miyano '73 Izosimova+ '78 Fassbender+ '98 Furukawa+ '65 JENDL/HE-2004

Isotope production cross Isotope production cross sections : sections : 1212C, C, 2727AlAl

16


100 101 102 1030

0.1

0.2

0.3

0.4

0.5


Cro

ss S

ecti

on [

b]56Fe(p,n)56Co

Zhuravlev+ '84 Levkovskij '91 Gadioli+ '74 Jenkins+ '70 Remsberg+ '63 Jenkins+ '70 Read '68 Vonach+ '97 Michel+ '83 Barchuk+ '87 Michel+ '86 Aleksandrov+ '90 Michel+ '89 Aleksandrov+ '95 Michel+ '97 Michel+ '95 Schiekel+ '96 Aleksandrov+ '89 Daum '96 LA150 NRG(Petten) JENDL/HE-2004

Isotope production cross sections : Isotope production cross sections : 5656FeFe

17


100 101 102 1030

0.1

0.2

0.3

0.4

0.5

0.6


Cro

ss S

ecti

on [

b]

63Cu(p,n)63Zn Levkovskij '91 Ghoshal '50 Meadows '53 Colle+ '74 Reuland+ '69 LA150 JENDL/HE-2004

100 101 102 1030

0.2

0.4

0.6

0.8

1


Cro

ss S

ecti

on [

b]

65Cu(p,n)65Zn Kopecky '85 Levkovskij '91 Gadioli+ '74 Colle+ '74 Williams+ '67 Greene+ '72 Albouy+ '62 Greenwood+ '84 Pulfer '79 Dmitriev+ '68 LA150 JENDL/HE-2004

Isotope production Isotope production cross sections : cross sections : 63,6563,65CuCu

18

Coupled Channel-Soft Rotator Model

Elastic scattering angular distributions : Elastic scattering angular distributions : 2828SiSiSample of Evaluated Results

19

Elastic scattering angular distributions : FeElastic scattering angular distributions : Fe

100

101

102

103

104

105

106

107

108

0 30 60 90 120 150 180

Mellema et al. (natFe)

Hjort et al. (56Fe)JENDL-HELA150NRG2003

Cro

ss s

ectio

ns [m

b/sr

]

Angle [deg.]

20 MeV x104

22 MeV x 103

24 MeV x 102

26 MeV x101

65 MeV

n + 56Fe

10-2

100

102

104

106

108

1010

0 30 60 90 120 150 180

Comparat et al.Ridley et al.Sakaguchi et al.Brussel et al.Johansson et al.JENDL-HELA150NRG2003

Cro

ss s

ectio

ns [m

b/sr

]

Angle [deg.]

30.3 MeV x104

p + 56Fe

39.8 MeV x103

65 MeV x102

156 MeV x10

182.4 MeV


20

90 MeV90 MeV

120 MeV120 MeV 160 MeV160 MeVDDXs ofDDXs of 90 90Zr(Zr(p,xnp,xn)


21

“Transmission through shields of quasi-monoenergetic neutrons generated by 43- and 68-MeV proton”

Ref.) N. Nakao et al., Nucl. Sci. Eng. 124, 228 (1996). H. Nakashima et al., Nucl. Sci. Eng. 124, 243 (1996).

Neutron spectrum of p+7Li

Experimental arrangement @ TIARA, JAERI

Benchmark Calculation (TIARA Experiment)

22

Shield

IronConcrete

20 cm

20 cm

10- 200 cm

Neutron Beam

by 43 MeV p+7Li

68 MeV p+7Li

Detector

JENDL High Energy File (JENDL/HE)JAPAN ATOMIC ENERGY

AGENCY

Ref.) N. Nakao et al., Nucl. Sci. Eng. 124, 228 (1996). H. Nakashima et al., Nucl. Sci. Eng. 124, 243 (1996).

Integral Experiment at TIARA

23

Ep=68 MeV

t=70 cm

Ep=68 MeV

t=100 cm

Integral Experiment at TIARAJAPAN ATOMIC ENERGY

AGENCY

24

(JENDL/HE-2007)

Benchmark Calculation (TIARA Experiment - Fe)

Ep=68 MeV

t=20 cm

Ep=68 MeV

t=40 cm

25

Benchmark Calculation (TIARA Experiment - Fe)

Ep=43 MeV

t=70 cm

Ep=43 MeV

t=40 cm

26

Benchmark Calculation (TIARA Experiment - Concrete)

Ep=68 MeV

t=100 cm

Ep=68 MeV

t=200 cm

27


Ep=68 MeV

t=25 cm

Ep=68 MeV

t=50 cm

28


Ep=43 MeV

t=25 cm

Ep=43 MeV

t=50 cm

29

Proton Energy: 250 and 350 MeVTarget: Graphite, Al, Fe, PbTarget Thickness ： full stopMeasured Angle: 0 °Flight Distance: 11.4 m and 67.8 m (250MeV proton), 11.4 m and 95 m (350MeV proton)

(NE213 organic liquid scintillator)

RCNP at Osaka Univ.

Integral Experiment at RCNPJAPAN ATOMIC ENERGY

AGENCY

Experimental Set-up

30

Code INC model Evaluated nuclear data Proton and neutron transition energy (MeV)

PHITS Bertini JENDL-HE 150

PHITS Bertini LA150 150

10-6

10-5

10-4

10-3

10-2

0 50 100 150 200 250

Expt.PHITS JENDL-HE + BertiniPHITS LA150 + Bertini

Neu

tron

yie

ld (

n/sr

/MeV

/pro

ton)

Neutron energy (MeV)

250MeV proton on a thick graphite27.5 cm thick ness

10-6

10-5

10-4

10-3

10-2

0 50 100 150 200 250

Expt.PHITS JENDL-HE + BertiniPHITS LA150 + BertiniN

eutr

on y

ield

(n/

sr/M

eV/p

roto

n)


250MeV proton on a thick iron7.5 cm thick ness

JENDL-HE: natC (including 13C)

LA150 : 12C


AGENCY

Case1: JENDL-HE vs. LA150

31

Code INC model Evaluated nuclear data Proton and neutron transition energy (MeV)

PHITS Bertini JENDL-HE 150

PHITS - JENDL-HE -

10-6

10-5

10-4

10-3

10-2

0 50 100 150 200 250

Expt.PHITS JENDL-HE + BertiniPHITS JENDL-HE

Neu

tron

yie

ld (

n/sr

/MeV

/pro

ton)


250MeV proton on a thick graphite27.5 cm thick ness

10-6

10-5

10-4

10-3

10-2

0 50 100 150 200 250

Expt.PHITS JENDL-HE + BertiniPHITS JENDL-HE N

eutr

on y

ield

(n/

sr/M

eV/p

roto

n)


250MeV proton on a thick iron7.5 cm thick ness


AGENCY

Case2: JENDL-HE with Bertini Parameterization vs. no INC

32

30km

～ 12km

Mt. Fuji(3.8km)

Proton (106 GeV) Shower in air

Blue: neutron, proton, etc.

Red: electron, photon, etc.

Green: muon

Proton

Application (Neutron Spectra Due to Cosmic-ray)JAPAN ATOMIC ENERGY

AGENCY

Database for Cosmic-ray Environment

Earth

Atmosphere Model

PHITSPHITS

JENDL/HE

33

Calculated Results

Application (Neutron Spectra Due to Cosmic-ray)JAPAN ATOMIC ENERGY

AGENCY

10–8 10–4 100 10410–4

10–2

100

測定値 (Goldhagen et al.)

シミュレーション (本研究)

高度 11.8 km

rc = 4.3 GV

海面レベルrc = 2.7 GV

中性

子フ

ラッ

クス

(cm

–2 s–

1 leth

arg

y–1 )

中性子エネルギー (MeV)

太陽活動極小期


unregistered

Height 11.8km

See Level

Measured by Nakamura+

Present Simulation

Neutron Energy

Neu

tron

Flu

x

10–8 10–4 100 10410–4

10–2

100

高度 11.28kmrc = 12 GV


測定値 (Nakamura, et al.)

シミュレーション (本研究)

計算式 (本研究)

高度 4.88kmrc = 12 GV

高度 70mrc = 10 GV中

性子

フラ

ック

ス (

cm–

2 s–1 le

tha

rgy–

1 )

中性子エネルギー (MeV)


太陽活動極大期

unregisteredHeight 11.8km

Height 4.88km

See Level

Measured by Nakamura+

Present Simulation

Present Formula

Neutron Energy

Neu

tron

Flu

x

34

Results of some analyses using JENDL/HE were also introduced as those of benchmark calculations. The calculations using JENDL/HE show good performance.

Reviewed were status of method and tools of the evaluation for JENDL High Energy File (JENDL/HE-2007). Especially, improvements of cluster-particle emission spectra and production cross sections.

SummaryJAPAN ATOMIC ENERGY

AGENCY

The first version of JENDL/HE (JENDL/HE-2004 ) has been released in March, 2004 with neutron and proton nuclear data of 66 nuclides up to 3 GeV. Soon, the second version of JENDL/HE (JENDL/HE-2007) will be released with those of totally over 100 nuclides including actinide ones.

35

The Status of Nuclear Data above 20 MeV Masayoshi SUGIMOTO, Tokio FUKAHORI Japan Atomic Energy Agency IAEA’s Technical Meeting on Nuclear Data Libraries.

Documents

proton energy

xa cross sectionbaba

alpha elastic scattering

triton elastic scattering

wide energy

deuteron elastic scatteringexp

64zn92mo107ag208pbangle

58ni90zr120sn208pbangle