Progress of implementation of the GDH model in TALYS

KIT – University of the State of Baden-Württemberg and

Large-scale Research Center of the Helmholtz Association www.kit.edu

A.Yu.Konobeyev, U.Fischer, P.E.Pereslavtsev

Institute for Neutron Physics and Reactor Technology, KIT

A.Koning

Nuclear Data Section, International Atomic Energy Agency

M.Blann

75-5865 Walua Rd. D423, Kailua-Kona, Hi. 967409, US

Progress of implementation of the GDH model

in TALYS

KIT – University of the State of Baden-Württemberg and

Large-scale Research Center of the Helmholtz Association www.kit.edu

Objective

Discussion of implementation of the GDH model in

the TALYS-1.74 code and results of test calculations

JEFF MeetingNovember 30 – December 4, 2015

3

The hybrid and GDH model

J.J.Griffin, 1966 : exciton model

G.D.Harp, J.M.Miller, B.J.Berne, 1968 : master-equation approach

M.Blann, 1971 : “hybrid” model

M.Blann, 1972 : “geometry dependent hybrid” model

M.Blann, H.K.Vonach, 1983 :modification and basic improvements

J.Bisplinghoff, 1986 : critical analysis. Configuration mixing

1994 : cluster emission in GDH

M.Blann, 1996 : hybrid Monte Carlo simulation (HMS)

2005 : improvement of cluster emission calculation in GDH

C.A.Soares Pompeia, B.V.Carlson, 2006 : new analysis

2010-2014 : further improvements

JEFF MeetingNovember 30 – December 4, 2015

4

Emission and transition rate

The GDH model

P U,C K O D

t t t t

d d d d

d d d d

0

e2 x

l n x nel 0 n nx x x

d ( p 1,h,U)( 2l 1)T X gD

d ( p,h,E )

inve x x x x xx 2 3

x

( 2S 1) ( ),

g

x 0 x l

V ( )

.

Cluster emission

0 10 20 30 40 50

100

101

Deuteron energy (MeV)

Bertrand (73)

Ep=61.5 MeV

D

K-O

F(1,1)

F(2,0)

EQ

TOTAL

54Fe(p,d)x

d

/d d

(m

b/M

eV

)

JEFF MeetingNovember 30 – December 4, 2015

5

GDH + TALYS

TALYS-1.0, 2009

TALYS-1.2, 2011

TALYS-1.74

Exciton and hybrid model

different, sometimes close results

Uncertainty of calculated cross-sections, particle distributions etc.

TALYS-1.0(p,x) reactions, targets Z=12-83,

Ep up to 150 MeV, N = 16,045

Factors default models

GDH

20.3 20.9

1.24 1.27

2.33 1.98

1.33 1.33

5.02

N

1iexp

i

calc

i

exp

i

N

1H

N

1iexp

i

calc

iCE

N

1R

N

1icalc

i

exp

iEC

N

1R

2/11

N

1i

2

expi

expi

N

1i

2

expi

expi

calci

expi

/

)lg()lg(

10S

JEFF MeetingNovember 30 – December 4, 2015

6

GDH in TALYS-1.74

GDH: mpreeqmode=5

different GDH options: subroutine gdhinput

added subprograms: 55

modified ALICE/ASH subprograms: 30

new subroutines: 10

modified TALYS subroutines: 6

Memory load: +0.15 %

JEFF MeetingNovember 30 – December 4, 2015

7

Nucleon energy distributions

JEFF MeetingNovember 30 – December 4, 2015

8

JEFF MeetingNovember 30 – December 4, 2015

9

JEFF MeetingNovember 30 – December 4, 2015

10

JEFF MeetingNovember 30 – December 4, 2015

11

JEFF MeetingNovember 30 – December 4, 2015

12

JEFF MeetingNovember 30 – December 4, 2015

13

JEFF MeetingNovember 30 – December 4, 2015

14

Deuteron energy distributions

JEFF MeetingNovember 30 – December 4, 2015

15

JEFF MeetingNovember 30 – December 4, 2015

16

JEFF MeetingNovember 30 – December 4, 2015

17

JEFF MeetingNovember 30 – December 4, 2015

18

JEFF MeetingNovember 30 – December 4, 2015

19

a-particle energy distributions

JEFF MeetingNovember 30 – December 4, 2015

20

JEFF MeetingNovember 30 – December 4, 2015

21

JEFF MeetingNovember 30 – December 4, 2015

22

JEFF MeetingNovember 30 – December 4, 2015

23

JEFF MeetingNovember 30 – December 4, 2015

24

Conclusion

The GDH model was implemented in TALYS-1.74

First test calculations were performed