Parameter. Library (RIPL) for nuclear reaction calculations

Transcription

1 IAEA Nuclear Data: the Reference Input Parameter Library (RIPL) for nuclear reaction calculations Roberto Capote, IAEA/NDS, Vienna, Austria

2 RIPL-II and RIPL-III participants M. Avrigeanu Inst. de Fizica si Inginerie Nucleara Horia Hulubei, Romania T. Belgya Institute of Isotope and Surface Chemistry, Hungary O. Bersillon Centre d Etudes Nucleaires de Bruyeres-le-Chatel, France R. Capote IAEA Nuclear Data Section T. Fukahori Nuclear Data Center, JAEA, Japan S. Goriely Institut d Astrophysique, Université Libre de Bruxelles, Belgium Y. Han China Institute of Atomic Energy, PR China M. Herman National Nuclear Data Center, BNL, USA S. Hilaire DPTA/SPN, CEA/DAM Ile de France, France A.V. Ignatyuk IPPE, Obninsk, Russian Federation S. Kailas Bhabha Atomic Research Centre, India A. Koning Nuclear Research and Consultancy Group, The Netherlands P. Obložinskỷ Brookhaven National Laboratory, USA V. A. Plujko Taras Shevchenko National University, Kiev, Ukraine RIPL-III participants E. S. Soukhovitskii Joint Institute of Energy and Nuclear Research, Belarus P. Talou Los Alamos National Laboratory, USA P. G. Young Los Alamos National Laboratory, USA G. Zhigang China Institute of Atomic Energy, PR China 2

3 A long time ago before RIPL (Recommended/any) inputs for nuclear reaction calculations? 3

4 RIPL Background Nuclear reaction theory: sufficiently advanced to meet most of the requirements for a number of applications Major sources of uncertainty are the input parameters needed to perform theoretical calculations RIPL Objective Improve the methodology of nuclear data evaluation by increasing predictive power, accuracy and reliability of theoretical calculations by nuclear reaction model codes Improved description of nuclear reactions, easier calculations allowing for a much better understanding 4

5 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) The longest running IAEA/NDS project 5

6 Reference Input Parameter Library Electronic Starter File (known as Reference Input Parameter Library-1) was developed and made available to users throughout the world in 1997 (compilation) : RIPL-1 starter file ( ) Second CRP was initiated on Nuclear Model Parameter Testing for Nuclear Data Evaluation (Reference Input Parameter Library: Phase II), and completed in Revision, extension and validation of the original RIPL-1 Starter File to produce a consistent RIPL-2 library of recommended input parameters : RIPL-2 database ( Main goal: Energy applications, E<20MeV 6

7 RIPL-3 additional requirements Reactions 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 7

8 Reference Input Parameter Library Third (and final) CRP: Parameters for Calculation of Nuclear Reactions of Relevance to Non-Energy Nuclear Application (Reference Input Parameter Library: Phase III) started in The project is close to completion. The update of the RIPL-2 database will be released in September : RIPL-3 database ( 8

9 1.- MASSES FRDM file including Audi et al. (2003) masses Skyrme HFB file (HFB-14) including Audi et al. (2003) masses Skyrme HFB spherical density distributions (HFB-14) Gogny HFB spherical density distributions (D1S) rms(m) = 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 9

10 Comparison with experimental masses (2149 nuclei: Audi, Wapstra & Thibault 2003) 4 M(Exp) M(HFB-14) ΔM [MeV] N HFB14 model: S. Goriely, M. Samyn, J.M. Pearson, (2007) Phys. Rev. C75,

11 HFB14 vs. experimental data S 208 Pb R th [fm] ρ ch [fm -3 ] R [fm] exp r [fm] r [fm] Charge radii Charge densities 11

12 Impact of the HFB pairing strength 10 2 on nuclear level densities at U=S n BSk9 BSk D th /D exp A BSk13=BSk HFB+combinatorial versus experimental s-wave spacings A 12

13 HFB14: A modified collective correction!! of particular relevance at large deformation --> Fission calculations!! E coll [MeV] a perturbative cranking correction for rotational correlations a phenomenological correction for vibrational correlations 240 Pu rotational ( ) E coll = be crank rot tanh cβ 2 rotational + vibrational E coll = be crank rot tanh( cβ 2 ) [ ] + de crank rot exp l(β 2 β 0 2 ) 2 β 2 13

14 Fission barriers vs. «experimental» data B in (Exp) B in (HFB) 52 nuclei with Z 88 rms = 1.2 MeV B out (Exp) B out (HFB) 45 nuclei rms = 2.0 MeV 14 NO VIBRATIONAL CORRECTION B(exp)-B(th) [MeV] B(exp)-B(th) [MeV] N

15 Fission barriers vs. «experimental» data B in (Exp) B in (HFB) 52 nuclei with Z 88 rms = 0.67 MeV B out (Exp) B out (HFB) 45 nuclei rms = 0.65MeV 15 WITH VIBRATIONAL CORRECTION B(exp)-B(th) [MeV] B(exp)-B(th) [MeV] N

16 2.- DISCRETE LEVELS SCHEME DERIVED FROM ENSDF Number of nuclei processed Number of records read in Number of levels processed Number of gammas are processed Number of unique spins Number of unique spins with parenthesis Number of spins inferred from gamma transitions Number of spins inferred from spin distributions Number of spins chosen from list by spin distributions Number of known ICC Number of newly determined ICC

17 3.- RESONANCES D 0 (RIPL-2)/D 0 (Mug2006) s-resonance spacing

18 Analysis of the resonance parameters for 238 U The 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]. Number of resonances D 0 =20.2+/-0.2 ev Number of resonances Sum of reduced widths, kev N r =897, N 0 =989+/-2 D 0 =20.2+/-0.2, S 0 =1.07+/ S 0 =1.07+/-0.07 sqrt(γ n 0 ), ev 1/2 U-238 s-wave Neutron energy, kev Neutron energy, kev

19 Average resonance parameters for 238 U: D 0, ev D 1, ev S 0, 10-4 S 1, , Gilbert-Cameron 17.7± , Rohr et al. 21.5± ± , Mughabghab 20.9± ± ± ± , Ignatyuk et al ± ± ± ± , Beijing, RIPL ± ± , RIPL-2 (10 kev) 20.8± ± ± ± , Leal et al., (20 kev) ± ± , Mughabghab 20.26± ± ± ± , present (20 kev) 20.2± ± ± ±

20 Analysis of the resonance parameters for 107 Ag 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 It is impossible to describe the PT-distribution with Do=14.9+/-.6 ev /Mug2006/. 250 Number of resonances N r =203, N 0 =290+/-5 D 0 =22.4+/-0.3, S 0 =0.47+/ sqrt(γ n 0 ), ev 1/2 Ag-107 s-wave Number of resonances Ag-107 s-wave D 0 =22.4+/-0.3 ev Sum of reduced widths, kev Ag-107 s-wave S 0 =0.47+/ Neutron energy, kev Neutron energy, kev 20

21 Average resonance parameters for 107 Ag D 0, ev D 1, ev S 0, 10-4 S 1, , Gilbert-Cameron 31± , Rohr et al. 24.0± ± , Mughabghab 16± ± ± , Ignatyuk et al. 22 ± ± ± , Beijing, RIPL ± ± , RIPL ± ± ± , Mughabghab 14.9± ± ± ± , present (B-VII) 22.4± ± ± ±0.3 21

22 List of analyzed differences: Nucleus Mug81/84 RIPL-2 Mug06 Present D 0, ev S 0 D 0, ev S 0 D 0, ev S 0 D 0, ev S 0 Cm ±.2 1.5±.3.75 ± ± ± ±.22.69± ±.15 Am-243m.40± ±.2.40± ±.2.29± ±.25.26± ±.1 Pa ±.15.65±.15.47± ±.26.48±.12.80±.15 Hg ±30 1.2± ±20.80±.19 70± ±.5 Hg ± ±35 1.3±.5 69± ±20 1.3±.3 Dy ± ± ±.4 2.7± ± ±.3 1.8±.4 Eu ± ± ±.6.25± ± ±.5 Te ±140.16± ±500.2± ±180*.16± ±100* -- Te ±30.25± ±150.2± ±375.26± ±300.20±.07 I ±3.5± ±1.4.42± ±.3.35±.05 I ±.8.8±.1 15±3.8±.2 9.7±.8.6± ±.5.71±.08 Ag ±3.38± ±.4.40± ±.6.46± ±.5.47±.05 Pd ±10*.40± ±50.25± ±43.47± ±50.30±.10 Pd ±4*.34± ±90.6±.3 174±25.42± ±32.65±.15 22

23 4.- OPTICAL MODEL 1. Dispersive CC potentials for nucleon induced reactions Rigid rotor: Actinides, W-Ta-Hf, Au, Mn, Rh, Soft rotor: Zr Capote, Soukhovitskii, et al ( ) 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 particles Alphas, Deuteron, Tritons 23

24 DISPERSIVE OMPs Molina, Capote, Quesada and Lozano PRC65(2002) Powerful CC OMP fitting code OPTMAN E. Soukhovitskii, S. Chiba, et al 24

25 Dispersive Coupled Channels OMP Soukhovitskii, Capote, Quesada and Chiba, PRC72 (2005) Dispersive coupled channel analysis of nucleon scattering from 232Th up to 200 MeV 25

26 Dispersive Coupled Channels OMP Capote, Soukhovitskii, Quesada and Chiba, PRC72 (2005) Is a global coupled-channel dispersive optical model potential for actinides feasible? 26

27 Dispersive Coupled Channels OMP Capote, Soukhovitskii, Quesada and Chiba, Varenna 2006; NEMEA-3 DCC OMP for tungsten nuclides 27 [σ tot (W-186) - σ tot (W-182)] / {[σ tot (W-186) + σ tot (W-182)]/2} Dietrich 2003 Guenther 1982 (shifted by ) Rigid rotor DCC OMP (shifted by ) ENERGY [MeV]

28 A global DCC OMP for actinides Capote, Soukhovitskii, Quesada, Chiba and Bauge, JNST45 (2008) 333 Complete table in Proceedings of the International Conference on Nuclear Data for Science 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 28

29 Are DCC OMPs Lane consistent? DCC OMP Isospin dependence Lane equations 29

30 Dispersive Coupled Channels OMP Capote, Soukhovitskii, Quesada and Chiba, PRC76 (2007) Approximate Lane consistency of the dispersive coupled-channels potential for actinides 232 Th(p,n) 238 U(p,n) MeV, Hansen 25.8 MeV, Schery Lane, sum of MeV, Hansen Lane, sum of dσ/dω(b/sr) dσ/dω(b/sr) θ(deg.) θ(deg.) 30

31 Zr-90 DCCOMP based on soft rotor SOFT ROTOR COUPLED SCHEME (7 levels) (3387) -- (1607) (1928) (3124) (2748) (>3400) NUDAT v

32 DCC OMP - soft rotor couplings Soukhovitskii, Capote, Quesada and Chiba, Unpublished Cross Section (barns) 10 ENDF/B-VII.0:ZR-90(N,TOT) BROND-2.2:ZR-90(N,TOT) DCCOMP (this work) 1993 Finlay 1985 Fedorov 1980 Pasechnik 1977 De 1977 Djumin 1975 Guenther 1973 Stooksberry 1973 Green n + 90 Zr Incident Energy (MeV) 32

33 DCC OMP - soft rotor couplings Soukhovitskii, Capote, Quesada and Chiba, Unpublished n + 90 Zr p + 90 Zr 90 Zr (Wang, E n =24.0 MeV) 90 Zr (Dickens, E p =12.7 MeV) dσ/dω(b/sr) ,E q =0. MeV 2 +,E q = MeV ,E q = MeV, data *0.01 Calculated data dσ/dω(b/sr) ,E q =0. MeV ,Eq=1.75 MeV 2 +,E q =2.18 MeV 3 -,E q =2.74MeV,data*10-3 Calculated data θ cm (deg.) θ cm (deg.) 33

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

35 35

36 A NEW global combinatorial NLD formula S. 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) s-wave p-wave D exp / D th exp. D 0 f rms = A 36

37 Renormalization factors to reproduce D 0 and cumulative levels ρ renorm (U) = e α U δ ρ HFB (U δ) α δ [MeV] A A 37

38 HFB LD vs OSLO data 38

39 Impact of LDs on cross section calculations Defined local and global systematics Unpublished (Koning, Hilaire, Goriely) 39

40 6.- GAMMA RAY STRENGTH FUNCTIONS Lorentzian, EGLO, MLO, SMLO, QRPA-HFB14 (See V. Plujko presentation tomorrow) The E1 gamma-decay strength function on 144 Nd for U=B n The E1 photoabsorption cross section on 144 Nd 40

41 Comparison of the photoabsorption cross section calculations on 40 Ca with exp.data (V.A. Erokhova et al Izvestiya RAN. Seriya Fiz. 67 (2003) 1479) Panel (a) shows calculations with GDR parameters from systematics (RIPL); (b) - calculations with GDR parameters obtained from fitting the exp. data. HFB-QRPA is microscopic approach given by S.Goriely et al. MSA - semi classical method proposed by V.Abrosimov, O.Davidoskaya. 41

42 Comparison of the photoabsorption cross sections on 208 Pb. Panel b shows the low-energy part of the cross sections. Experimental data are taken from A. Veyssiere, H. Beil, R. Bergere, P. Carlos, A. Lepretre, Nucl. Phys. A159 (1970) 561 in panel a and from V.V. Varlamov, M.E. Stepanov, V.V. Chesnokov, Izvestiya RAN. Seriya Fiz. 67 (2003) 656. in panel b. The SLO parameters are taken from the RIPL-2 library. 42

43 B in (Exp) B in (HFB) 52 nuclei with Z 88 B out (Exp) B out (HFB) 45 nuclei 43 B(exp)-B(th) [MeV] 7.- FISSION HFB14 fission barriers vs. «experimental» data rms = 0.67 MeV rms = 0.65MeV B(exp)-B(th) [MeV] N

44 HFB14 Projection of the static path along the quadrupole deformation parameter β The U isotopes U235 U236 U237 U238 U239 U The Cm isotopes Cm242 Cm243 Cm244 Cm245 Cm246 Cm247 Cm248 Cm249 Cm250 Cm251 E-E GS [MeV] E-E GS [MeV] β 2 β 2

45 The Cm isotopes in the very n-rich region: 270 A Cm: N=184 shell closure E-E GS [MeV] Cm270 Cm272 Cm274 Cm276 Cm278 Cm β 2 45

46 NEW EMPIRE VERSION

47 IMPROVED FISSION MODELLING: BARRIERS + WELLS (absorption) Full damping vs Partial damping Cross Section (barns) Incident Energy (MeV) EMPIRE Part.damp. EMPIRE full damp. EMPIRE no class-ii 2006 N_tof 2001 Shcherbakov 1991 Fursov 1986 Kanda 1986 Goverdovskij 1985 Anand 1985 Kanda 1983 Meadows 1982 Behrens 1980 Blons 1980 D hondt 1978 Nordborg 1975 Blons 1971 Muir 1946 Williams 47

48 235 U(n,f) RIPL-2 Z A Va Vb original HFB Z A Va Vb Vc Δgs ags normalized HFB Z A Va Vb Vc Δgs ags Δsdl asdl

49 238 U(n,f) RIPL-2 Z A Va Vb original HFB Z A Va Vb Vc Δgs ags normalized HFB Z A Va Vb Vc Δgs ags Δsdl asdl

50 238 Pu(n,f) RIPL-2 Z A Va Vb original HFB Z A Va Vb Δgs ags normalized HFB Z A Va Vb Δgs ags Δsdl asdl

51 CONCLUDING REMARKS Over 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. 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. 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 Efforts will continue to develop this database further, and monitor all studies that impact and could possibly improve their contents. RIPL-4 hopefully will contain SMMC level densities 51

52 Nuclear Data Production Experimental data: masses, discrete levels, deformations Model parameters: OMP, NLD, gamma, fission,etc. RIPL EMPIRE 2.19 (BNL/IAEA) / TALYS (NRG) / GNASH (LANL) FINAL GOAL: EVALUATION (ENDF-6 formatted file) or NUCLEAR REACTION CALCULATION 52

53 LD: IBM collective states R. Capote, A. Ventura, F. Cannata, J.M. Quesada, Phys Rev C71, (2005) Level densities of transitional Sm nuclei (Monte Carlo combinatorial + IBM coll) 53

54 LD: IBM collective states R. Capote, A. Ventura, F. Cannata, J.M. Quesada, Phys Rev C71, (2005) Level densities of transitional Sm nuclei 54