reference input parameter library (ripl) for nuclear ...inpp.ohiou.edu/~snp2008/talks/capote.pdf ·...

IAEA Nuclear Data: IAEA Nuclear Data: the the RReference eference IInput nput PParameter arameter LLibrary ibrary (RIPL) for nuclear reaction calculations(RIPL) for nuclear reaction calculations

Roberto Capote, IAEA/NDS, Vienna, AustriaRoberto Capote, IAEA/NDS, Vienna, Austria

SNP2008, July 8-11Athens, OH, USA

Roberto Capote, IAEA Nuclear Data SectionE-mail: [email protected]

M. Avrigeanu Inst. de Fizica si Inginerie Nucleara “Horia Hulubei”, Romania T. Belgya Institute of Isotope and Surface Chemistry, HungaryO. Bersillon Centre d’Etudes Nucleaires de Bruyeres-le-Chatel, FranceR. Capote IAEA Nuclear Data SectionT. Fukahori Nuclear Data Center, JAEA, JapanS. Goriely Institut d’Astrophysique, Université Libre de Bruxelles, BelgiumY. Han China Institute of Atomic Energy, PR China M. Herman National Nuclear Data Center, BNL, USAS. Hilaire DPTA/SPN, CEA/DAM Ile de France, FranceA.V. Ignatyuk IPPE, Obninsk, Russian FederationS. Kailas Bhabha Atomic Research Centre, IndiaA. Koning Nuclear Research and Consultancy Group, The NetherlandsP. Obložinskỷ Brookhaven National Laboratory, USAV. A. Plujko Taras Shevchenko National University, Kiev, Ukraine E. S. Soukhovitskii Joint Institute of Energy and Nuclear Research, Belarus P. Talou Los Alamos National Laboratory, USAP. G. Young Los Alamos National Laboratory, USAG. Zhigang China Institute of Atomic Energy, PR China

RIPL-II and RIPL-III participants

RIPL-III participants



A long time ago before RIPL …

(Recommended/any) inputs for nuclear reaction calculations ?



RIPL BackgroundNuclear reaction theory: sufficiently advanced to meet

most of the requirements for a number of applicationsMajor sources of uncertainty are the input parameters

needed to perform theoretical calculations

Improve the methodology of nuclear data evaluation by increasing predictive power, accuracy and reliability of theoretical calculations by nuclear reaction model codes

RIPL Objective

Improved description of nuclear reactions, easier calculations allowing for a much better understanding



IAEA Nuclear Data Section has addressed these needs through a series of Coordinated Research Projects dedicated to the production of a Reference Input Parameter Library (RIPL)

1994 – 2008The longest running IAEA/NDS project



Reference Input Parameter Library

1994-1997: RIPL-1 starter file (http://www-nds.iaea.org/ripl/ )

Second CRP was initiated on “Nuclear Model Parameter Testing for Nuclear Data Evaluation (Reference Input Parameter Library: Phase II)”, and completed in 2003. Revision, extension and validation of the original RIPL-1 Starter File to produce a consistent RIPL-2 library of recommended input parameters.

Electronic Starter File (known as Reference Input Parameter Library-1) was developed and made available to users throughout the world in 1997 (compilation)

1998-2003: RIPL-2 database (http://www-nds.iaea.org/RIPL-2/)Main goal: Energy applications, E<20MeV



RIPL-3 additional requirementsReactions at high energies for ADS (up to 150

MeV), production of medical radioisotopes (up to 100 MeV) and radiotherapy (up to 250 MeV)

Reactions on nuclei far from stability for ADS and astrophysics

Charged-particle reactions for all non-energy applications

Number of simple routines for the calculation of basic input data from the parameters contained in the library will be provided to reduce a risk of misusing the parameters



Reference Input Parameter LibraryThird (and final) CRP: “Parameters for Calculation of Nuclear Reactions of Relevance to Non-Energy Nuclear Application (Reference Input Parameter Library: Phase III)” started in 2003. The project is close to completion. The update of the RIPL-2 database will be released in September 2008. 2003-2008: RIPL-3 database (http://www-nds.iaea.org/RIPL-3/)



1.- MASSESFRDM file including Audi et al. (2003) massesSkyrme HFB file (HFB-14) including Audi et al. (2003) massesSkyrme HFB spherical density distributions (HFB-14)Gogny HFB spherical density distributions (D1S)

rms(M) = 650-750 keV on 2149 (Z ≥ 8) experimental masses (Audi et al., 2003)

To be compared with- FRDM predictions: rms(M) = 676 keV (2149 Z ≥ 8 nuclei)- Previous HF predictions:

Traditional Skyrme forces: rms(M) >> 2 MeV (120 e-e sph)Ex. Oak Ridge "Mass Table" based on HFB with SLy4

rms(M)=4.7MeV on 570 e-e sph+def nuclei



Comparison with experimental masses(2149 nuclei: Audi, Wapstra & Thibault 2003)

-4

-2

0

2

4

0 20 40 60 80 100 120 140 160

ΔM

[MeV

]

N

M(Exp)–M(HFB-14)

HFB14 model: S. Goriely, M. Samyn, J.M. Pearson, (2007) Phys. Rev. C75, 064312



HFB14 vs. experimental data

2.5

3

3.5

4

4.5

5

5.5

6

6.5

2.5 3 3.5 4 4.5 5 5.5 6 6.5R

exp [fm]

Rth [

fm]

0

0.02

0.04

0.06

0.08

0.1

0 1 2 3 4 5 6 7 8ρ ch

[fm

-3]

r [fm]

32S

1 2 3 4 5 6 7 8r [fm]

208Pb

Charge radii Charge densities



10-2

10-1

100

101

102

0 50 100 150 200 250

Dth/D

exp

BSk9

A50 100 150 200 250

A

BSk13

Impact of the HFB pairing strengthon nuclear level densities at U=Sn

HFB+combinatorial versus experimental s-wave spacings

BSk13=BSk14



0

1

2

3

4

5

6

7

8

0 0.5 1 1.5 2

E coll [M

eV]

β2

240Pu

HFB14: A modified collective correction!! of particular relevance at large deformation --> Fission calculations !!

• a perturbative cranking correction for rotational correlations• a phenomenological correction for “vibrational” correlations

Ecoll = b Erotcrank tanh cβ2( )

Ecoll = b Erotcrank tanh cβ2( )

+ d Erotcrank exp −l(β2 −β2

0)2[ ]

rotational

rotational + vibrational



Fission barriers vs. « experimental » data

52 nuclei with Z ≥ 88

45 nucleirms = 2.0 MeV

rms = 1.2 MeV

NO “VIBRATIONAL” CORRECTION

-4-3-2-101234

B(e

xp)-

B(th

) [M

eV]

-4-3-2-10123

135 140 145 150 155

B(e

xp)-

B(th

) [M

eV]

N

Bin(Exp) – Bin(HFB)

Bout(Exp) – Bout(HFB)



Bin(Exp) – Bin(HFB)52 nuclei with Z ≥ 88

-4-3-2-10123

135 140 145 150 155

B(e

xp)-

B(th

) [M

eV]

N

-4-3-2-101234

B(e

xp)-

B(th

) [M

eV]

Bout(Exp) – Bout(HFB)45 nucleirms = 0.65MeV

rms = 0.67 MeV

WITH “VIBRATIONAL” CORRECTIONFission barriers vs. « experimental » data



2.- DISCRETE LEVELS SCHEME

124328117675Number of newly determined ICC2494523184Number of known ICC

67796058Number of spins chosen from list by spin distributions36783516Number of spins inferred from spin distributions44413581Number of spins inferred from gamma transitions

3250738912Number of unique spins with parenthesis6333968858Number of unique spins

203449200944Number of gammas are processed135406138595Number of levels processed

2113877Number of records read in3020Number of nuclei processed20072005DERIVED FROM ENSDF



3.- RESONANCES

0 50 100 150 200 2500

2

D0(RIPL-2)/D0(Mug2006)s-resonance spacing



Analysis of the resonance parameters for 238UThe set of resonances at the energy region up to 20 keV contains 898 s-wave resonances, 849 p-wave resonances with J=1/2 and 1565 p-wave resonances with J=3/2 [L.Leal et al., Nuclear Data for Science and Technology, Santa Fe, 2004, p.276].

0 50 100 150

10

100

1000

Num

ber o

f res

onan

ces

sqrt(Γn0), eV1/2

U-238s-wave

Nr=897, N0=989+/-2D0=20.2+/-0.2, S0=1.07+/-0.07

0 5 10 15 200

200

400

600

800

1000

Num

ber o

f res

onan

ces

Neutron energy, keV

D0=20.2+/-0.2 eV

0 5 10 15 200

5

10

15

20

Sum

of r

educ

ed w

idth

s, k

eV

Neutron energy, keV

S0=1.07+/-0.07



Average resonance parameters for 238U:D0, eV D1, eV S0, 10-4 S1, 10-4

1965, Gilbert-Cameron 17.7±0.7 -- -- --1979, Rohr et al. 21.5±2.2 -- 1.02 ±0.16 --1984, Mughabghab 20.9±1.1 7.2±0.4 1.2±0.1 1.7±0.31986, Ignatyuk et al. 21.7 ±0.9 7.3±0.5 1.15±0.12 1.7±0.51996, Beijing, RIPL-1 21.0 ±0.05 -- 0.93±0.06 --2002, RIPL-2 (10 keV) 20.8±0.3 7.7±1.0 1.03±0.08 1.6±0.42004, Leal et al., (20 keV) -- -- 1.07±0.07 1.71±0.072006, Mughabghab 20.26±0.72 7.42±0.23 1.29±0.13 2.17±0.192007, present (20 keV) 20.2±0.2 7.59±0.05 1.02±0.02 1.62±0.05



The set of resonances at the energy region up to 6.5 keV contains 203 s-wave resonances and 196 p-wave resonances, which were inserted into the ENDF/B-VII file from the Mughabghab-2006. It is impossible to describe the PT-distributionwith Do=14.9+/-.6 eV /Mug2006/.

0 50 1001

10

100

Num

ber o

f res

onan

ces

sqrt(Γn0), eV1/2

Ag-107s-wave

Nr=203, N0=290+/-5D0=22.4+/-0.3, S0=0.47+/-0.05

0 2 4 60

50

100

150

200

250

Num

ber o

f res

onan

ces

Neutron energy, keV

Ag-107s-wave

D0=22.4+/-0.3 eV0 2 4 60

1

2

3

Sum

of r

educ

ed w

idth

s, k

eV

Neutron energy, keV

Ag-107s-wave

S0=0.47+/-0.05

Analysis of the resonance parameters for 107Ag



D0, eV D1, eV S0, 10-4 S1, 10-4

1965, Gilbert-Cameron 31±6 -- -- --1979, Rohr et al. 24.0±2.8 -- 0.41 ±0.13 --1984, Mughabghab 16±3 -- 0.38±0.07 3.8±0.61986, Ignatyuk et al. 22 ±2 -- 0.42±0.05 3.8±0.51996, Beijing, RIPL-1 22.6 ±0.09 -- 0.54±0.04 --2002, RIPL-2 22.0±0.4 -- 0.40±0.06 3.8±0.82006, Mughabghab 14.9±0.6 8.49±0.25 0.46±0.05 3.76±0.312007, present (B-VII) 22.4±0.5 9.1±0.6 0.47±0.05 2.2±0.3

Average resonance parameters for 107Ag



List of analyzed differences:Nucleus Mug81/84 RIPL-2 Mug06 Present

D0, eV S0 D0, eV S0 D0, eV S0 D0, eV S0

Cm-243 1.1±.2 1.5±.3 .75 ±.15 1.5±.3 1.11±.07 1.20±.22 .69±.06 1.25 ±.15 Am-243m .40±.08 1.3±.2 .40±.08 1.3±.2 .29±.02 1.47±.25 .26±.06 1.3±.1 Pa-232 -- -- .75±.15 .65±.15 .47±.04 1.22±.26 .48±.12 .80±.15Hg-201 -- -- 90±30 1.2±. 5 233±20 .80±.19 70±30 1. 3±.5Hg-198 105±33 -- 105±35 1.3±.5 69±3 -- 100±20 1.3±.3 Dy-156 2.7±.4 -- 4.8±1.6 1.8±.4 2.7±.4 3.3 ±1.3 2.8±.3 1.8±.4 Eu-152 .25±.03 -- .56±.10 1.4±.6 .25±.03 -- .35±.05 1.9±.5Te-130 870±140 .16±.05 1500±500 .2±.1 1130±180* .16±.05 1040±100* --Te-128 260±30 .25±.10 740±150 .2±.1 1510±375 .26±.15 1460±300 .20±.07 I-129 -- -- 30±3 .5±.1 19.0±1.4 .42±.07 13.6±.3 .35±.05 I-127 9.7±.8 .8±.1 15±3 .8±.2 9.7±.8 .6±.1 15.4±.5 .71±.08Ag-107 16±3 .38±.07 22.0±.4 .40±.06 14.9±.6 .46±.05 22.4±.5 .47±.05Pd-110 95±10* .40±.06 150±50 .25±.15 334±43 .47±.17 280±50 .30±.10 Pd-106 67±4* .34±.04 270±90 .6±.3 174±25 .42±.19 144±32 .65±.15



4.- OPTICAL MODEL1. Dispersive CC potentials for nucleon induced reactions

Rigid rotor: Actinides, W-Ta-Hf, Au, Mn, Rh, …Soft rotor: ZrCapote, Soukhovitskii, et al (2005-2008)Kunieda et al (2008)

2. Soft rotor CC OMPs (Soukhovitskii et al, 2004)3. Global dispersive spherical potentials

Neutrons – Morillon & Romain (2005)Protons – Li & Cai (2008)

4. OMPs for complex charged particlesAlphas, Deuteron, Tritons



DISPERSIVE OMPs

Molina, Capote, Quesada and Lozano PRC65(2002) 034616

+ Powerful CC OMP fitting code OPTMANE. Soukhovitskii, S. Chiba, et al



Dispersive Coupled Channels OMP

Dispersive coupled channel analysis of nucleon scattering from 232Th up to 200 MeV

Soukhovitskii, Capote, Quesada and Chiba, PRC72 (2005) 024604




Is a global coupled-channel dispersive optical model potential for actinides feasible?

Capote, Soukhovitskii, Quesada and Chiba, PRC72 (2005) 064210



1 10 100-0.04

-0.02

0.00

0.02

0.04

0.06

Dietrich 2003 Guenther 1982 (shifted by -0.008) Rigid rotor DCC OMP (shifted by +0.006)

[σto

t(W-1

86) -

σto

t(W-1

82)]

/{[σ

tot(W

-186

) + σ

tot(W

-182

)]/2}

ENERGY [MeV]


DCC OMP for tungsten nuclidesCapote, Soukhovitskii, Quesada and Chiba, Varenna 2006; NEMEA-3



A global DCC OMP for actinidesCapote, Soukhovitskii, Quesada, Chiba and Bauge, JNST45 (2008) 333

Complete table in Proceedings of the International Conference on Nuclear Data forScience and Technology, April 22-27, 2007, Nice, France, EDP Sciences, 2008

DCC OMPs for 31 actinides, tungsten and tantalum nuclei derived using approximated Lane consistent formulation with Coulomb corrections



Are DCC OMPs Lane consistent ?

Lane equations

DCC OMPIsospin dependence




Approximate Lane consistency of the dispersive coupled-channels potential for actinidesCapote, Soukhovitskii, Quesada and Chiba, PRC76 (2007) 057602

238U(p,n)

θ(deg.)

0 30 60 90 120 150 180

dσ/d

Ω(b

/sr)

10-5

10-4

10-3

26 MeV, Hansen0+

2+

4+

0++2++4+

Lane, sum of 0++2++4+

232Th(p,n)

θ(deg.)

0 30 60 90 120 150 180

dσ/d

Ω(b

/sr)

10-5

10-4

10-3

26 MeV, Hansen25.8 MeV, Schery0++2++4+

0+

2+

4+

Lane, sum of 0++2++4+



(3387)

(1607)(1928)

(2748)

(3124)

(>3400)

SOFT ROTOR COUPLED SCHEME (7 levels)

--

NUDAT v 2.4

Zr-90 DCCOMP based on soft rotor



DCC OMP - soft rotor couplings

Incident Energy (MeV)

Cro

ss S

ecti

on

(b

arn

s)

1 10 102

10

ENDF/B-VII.0:ZR-90(N,TOT)BROND-2.2:ZR-90(N,TOT)DCCOMP (this work)

1993 Finlay1985 Fedorov1980 Pasechnik1977 De1977 Djumin1975 Guenther1973 Stooksberry1973 Green

n + 90Zr

Soukhovitskii, Capote, Quesada and Chiba, Unpublished



90Zr (Wang, En=24.0 MeV)

θcm(deg.)

0 30 60 90 120 150 180

dσ

/dΩ

(b/s

r)

10-6

10-5

10-4

10-3

10-2

10-1

1000+, Eq=0. MeV

2+, Eq=2.1865 MeV

3-, Eq=2.7479 MeV, data *0.01

Calculated data

90Zr (Dickens, Ep=12.7 MeV)

θcm(deg.)

0 30 60 90 120 150 180

dσ

/dΩ

(b/s

r)

10-7

10-6

10-5

10-4

10-3

10-2

10-1

100

101

102

0+, Eq=0. MeV

02+, Eq=1.75 MeV

2+, Eq=2.18 MeV

3-, Eq=2.74 MeV, data *10-3

Calculated data

n + 90Zr p + 90Zr

DCC OMP - soft rotor couplingsSoukhovitskii, Capote, Quesada and Chiba, Unpublished



5.- LEVEL DENSITIESBased on OBSERVABLES:RIPL-3 discrete levels (2) and Neutron resonances (3)



A NEW global combinatorial NLD formulaS. Hilaire & S.Goriely (2007)

• Particle-hole as well as total parity-, spin- and E-dependent NLD• Deviation from the statistical limit at low energies (discrete counting)

50 100 150 200 25010-2

10-1

100

101

102

A

Dex

p / D

th

http://www-astro.ulb.ac.be/Html/nld_comb_ph.html

292 exp. D0

frms=2.30

s-wavep-wave



-2

-1.5

-1

-0.5

0

0.5

1

1.5

2

50 100 150 200 250

α

A

-4

-2

0

2

4

50 100 150 200 250

δ [M

eV]

A

Renormalization factors to reproduceD0 and cumulative levels

ρrenorm (U) = eα U−δ ρHFB(U −δ)



HFB LD vs OSLO data



Defined local and global systematicsUnpublished (Koning, Hilaire, Goriely)

Impact of LDs on cross section calculations



6.- GAMMA RAY STRENGTH FUNCTIONSLorentzian, EGLO, MLO, SMLO, QRPA-HFB14

(See V. Plujko presentation tomorrow)The E1 gamma-decay strength

function on 144Nd for U=Bn

The E1 photoabsorptioncross section on 144Nd



Comparison of the photoabsorptioncross section calculations on 40Cawith exp.data (V.A. Erokhova et alIzvestiya RAN. Seriya Fiz. 67 (2003) 1479)

Panel (a) shows calculations withGDR parameters from systematics(RIPL); (b) - calculations withGDR parameters obtained fromfitting the exp. data.

HFB-QRPA is microscopicapproach given by S.Goriely et al.

MSA - semi classical methodproposed by V.Abrosimov,O.Davidoskaya.



Comparison of the photoabsorptioncross sections on 208Pb. Panel b shows the low-energy partof the cross sections. Experimentaldata are taken from A. Veyssiere, H.Beil, R. Bergere, P. Carlos, A. Lepretre, Nucl.Phys. A159 (1970) 561 in panel a andfrom V.V. Varlamov, M.E. Stepanov, V.V.Chesnokov, Izvestiya RAN. Seriya Fiz. 67(2003) 656. in panel b. The SLOparameters are taken from theRIPL-2 library.



7.- FISSION

Bin(Exp) – Bin(HFB)52 nuclei with Z ≥ 88

-4-3-2-10123

135 140 145 150 155

B(e

xp)-

B(th

) [M

eV]

N

-4-3-2-101234

B(e

xp)-

B(th

) [M

eV]

Bout(Exp) – Bout(HFB)45 nucleirms = 0.65MeV

rms = 0.67 MeV

HFB14 fission barriers vs. « experimental » data



0

1

2

3

4

5

6

7

8

0 0.5 1 1.5 2 2.5 3

U235

U236

U237

U238

U239

U240

E-E G

S [MeV

]

β2

The U isotopes

Projection of the static path along the quadrupole deformation parameter β2

HFB14

0

1

2

3

4

5

6

7

8

0 0.5 1 1.5 2

Cm242Cm243Cm244Cm245Cm246

Cm247Cm248Cm249Cm250Cm251

E-E G

S [MeV

]

β2

The Cm isotopes



0

2

4

6

8

10

0 0.5 1 1.5 2 2.5

Cm270Cm272Cm274Cm276Cm278Cm280

E-E G

S [MeV

]

β2

280Cm: N=184 shell closureThe Cm isotopes in the very n-rich region: 270 ≤ A ≤ 280



NEW EMPIRE VERSION 3.0



Full damping vs Partial damping.

Incident Energy (MeV)

Cro

ss S

ecti

on

(b

arn

s)

10-1 1

10-6

10-5

10-4

10-3

10-2

10-1

EMPIRE Part.damp.EMPIRE full damp.EMPIRE no class-II

2006 N_tof2001 Shcherbakov1991 Fursov1986 Kanda1986 Goverdovskij1985 Anand1985 Kanda1983 Meadows1982 Behrens1980 Blons1980 D’hondt1978 Nordborg1975 Blons1971 Muir1946 Williams

IMPROVED FISSION MODELLING:BARRIERS + WELLS (absorption)



original HFB

normalized HFB

0.5450.4902.2705.9755.04323492

0.609

0.128

Δgs

-

-

Vc

0.4986.0295.51723692

0.0356.2735.39123592

agsVbVaAZ

0.000.000.5450.4902.2705.9755.44023492

0.609

0.400

Δgs

0.498

-0.080

ags

0.00

0.00

Δsdl

-

-

Vc

0.006.0295.51723692

0.005.8005.54523592

asdlVbVaAZ

RIPL-2

23492

5.6705.00023692

6.0005.25023592

VbVaAZ

235U(n,f)



238U(n,f)RIPL-2

5.5006.30023892

6.1506.40023792

6.0006.45023992

VbVaAZ

original HFB

0.371

0.149

0.253Δgs

3.681

3.808

4.219

Vc

0.5386.5425.99023992

0.4806.4775.92823892

0.2016.4135.55323792

agsVbVaAZ

normalized HFB

0.371

0.149

0.253Δgs

0.538

0.480

0.201

ags

0.00

0.00

0.00Δsdl

3.681

3.745

4.219

Vc

0.006.0686.07423992

0.006.1725.80223892

0.005.9355.92523792

asdlVbVaAZ



238Pu(n,f)RIPL-2

5.2435.96023894

5.0655.54523794

5.3315.96323994

VbVaAZ

original HFB

normalized HFB0.450

0.000

0.000Δgs

0.1825.3315.96323994

0.4005.2435.96023894

0.1345.0655.54523794

agsVbVaAZ

1.000

0.000

0.000Δgs

1.600

1.200

0.134

ags

0.00

0.00

0.00Δsdl

0.005.8406.05023994

0.005.2435.96023894

0.005.0655.34923794

asdlVbVaAZ



CONCLUDING REMARKSOver many years, IAEA staff within the Nuclear Data Section have successfully initiated and overseen the completion of various projects dedicated to satisfying a wide range of user demands for enhancements in the quantification and quality of neutron reaction cross sections.

RIPL-4 hopefully will contain SMMC level densities

The RIPL-3 database represents considerable advancements in the adoption and use of evaluated and highly credible nuclear data both for energy and non-energy applications

One of the most significant database developments have involved important advances in the evolution of a complete and consistent set of input parameters for the calculation of a wide range of nuclear reactions.

Efforts will continue to develop this database further, and monitor all studies that impact and could possibly improve their contents.



RIPL

Experimental data: masses, discrete levels, deformationsModel parameters: OMP, NLD, gamma, fission,etc.

FINAL GOAL: EVALUATION (ENDF-6 formatted file)or NUCLEAR REACTION CALCULATION

Nuclear Data Production

EMPIRE 2.19 (BNL/IAEA) / TALYS (NRG) / GNASH (LANL)



LD: IBM collective states

R. Capote, A. Ventura, F. Cannata , J.M. Quesada, Phys Rev C71, 064320 (2005)“Level densities of transitional Sm nuclei” (Monte Carlo combinatorial + IBM coll)



LD: IBM collective states

R. Capote, A. Ventura, F. Cannata , J.M. Quesada, Phys Rev C71, 064320 (2005)“Level densities of transitional Sm nuclei”

reference input parameter library (ripl) for nuclear ...inpp.ohiou.edu/~snp2008/talks/capote.pdf ·...

Documents