Recent Activities Related to Nuclear
Data Evaluation and Calculations inChina
Sun
Sun Weili
WeiliInstitute
Institute ofof Applied
Applied Physics
Physics andand Computational
Computational
Mathematics
Mathematics Beijing
Beijing 100088,100088, China
China

China Nucl. Data Committee

China Nucl. Data Center
Evaluation
BenchmarkMeasurement

NetworkSichuan U.
Lan Zhou U.
Tsing Hua U.
Nan Kai U.
IAPCM
……

Chinese Evaluated Nuclear Data Library(CENDL)
No.of Nuclides of CENDL
250
20015010050020036
CENDL-1

5567

CENDL-2

CENDL-2.1

CENDL-3.0

1995--2000

CONTENT OF CENDL- 3.0
140115120CENDL-3.0
CENDL-2.1
CENDL-2100
80
60
40
2033
22 22 21191010 1012818
5890
Light

StructuralFission Products& Medium

Heavy

Actinides

Nucl.

Content
1,2,3H,3,4He,6,7Li,9Be,10,11B,12C,14N,16O,19F,23Na,natMg,27Light Elements

Al, natSi, 31P,natS,natCl,natK,natCa
natTi,natV,50,52,53,54,natCr,55Mn,54,56,57,58,natFe,59Co,58,60,61,62,6

Structural
Materials

4,natNi,63,65,natCu,natZn
69,71,natGa,83,84,85,86Kr,85,87,natRb,88,89,90Sr,89,91Y,90,91,92,93,94,95,96,natZr,93,95Nb,95,97,98,100,natMo,99Tc,99,100,101,102,103,104,105

Ru,103Rh,105,108Pd,107,109,natAg,113,natCd,115,natIn,natSn,121,12
Fission Products & 3,natSb,130Te,127I,124,129,131,132,134,135,136Xe,133,134,135,137Cs,130,1
Medium Elements 32,134,135,136,137,138,natBa,139La,140,141,142,144Ce,141Pr,142,143,144,1
45,146,147,148,150,natNd,147,148,149Pm,144,147,148,149,150,151,152,154,nat
Sm,151,153,154,155,natEu,152,154,155,156,157,158,160,natGd,164Dy
Heavy Elements
Actinides

natLu,natHf,181Ta,natW,197Au,natHg,natTi,204,206,207,207,natPb,233,234,235,236,238,239U, 237Np,238,239,240,241,242Pu, 241,242Am,
249Bk,249Cf

Table 1. The Nuclides of CENDL-3.1 (>300 )1. Covariance study
Code:
EXPOV:SPC :Covariance Matrix for exp. dataSpline fitting
58,60,61,62,64,nat
58,60,61,62,64,natNi, 63,65,nat
63,65,natCu,

RAC :

27
27Al

R matrix theory for light nuclide
6，76，7Li , 10，11
10，11B, 16
16O

2. Data file for ADS
Code: MEND
- Neutron:

20 -- 250MeV

50,52,53,54Cr, 54,56,57,58Fe, 90,91,92,94,96Zr,
180,182,183,184,186W, 204,206,207,208Pb,238U.

- Proton: thresholds – 250MeV
54,56,57,58Fe, 180,182,183,184,186W,204,206,207,208Pb, 209Bi, 238U.

3. Calculation
for Light Nuclei Reaction(< 20 MeV)

multi-particle
recoil effect
3-body Break-up

n+A(A+1)*

R1*

R2*

Methods
 Unified Hauser-Feshbach and Exciton model

Zhang J. S et al， NSE. 114 (1993) 55

2- body process;3-body break-up process;multi-particle emission; recoil effect;
kinematics & dynamics;

 Model theory for light nucleiZhang J. S. et al， NSE. 133 (1999) 218

Data file (CENDL-3.1)
 6 Li, 7 Li, 9Be, 10B, 11B,

12C, 16 O, 19 F

( 1 p-shell)
Particularly, File 6 (DDX for all emittedparticles ) is established, where the energybalance is strictly considered.
Code: LUNF

n + 6 Li

Neutron DDX

n + 16 O

Neutron DDX

n + 19 F

Neutron DDX

Energy Spectra of (n, 2n) reaction

ENDF/B6.8

11Spectra of B(n, np)

reaction

Energy Spectra of 11B(n, nα) reaction

Neutron spectrum for the FNS C12CNeutron spectrum for the FNS LiO2
7Li

5He emission
 Generally,n, p, d, t, 3 He, 4He emission

 5He emission,

possibility ?

Experimentally,
Turk et al. Nucl. Phys.A431( 1984) 31

Theoretically,Zhang J.S. Sci. in China Ser.
G. Phy. 47 (2004) 137

Threshold3He vs 5He

95He ⎯⎯→ n + α∗

核素9Be

3He

5He

24.04

3.74

10B

17.34

5.89

11B

25.27

10.43

12C

21.1

8.95

14N

18.62

13.41

16O

15.64
17.08

8.56
5.17

19F

Open Channels (<20 MeV)
¾ n + 10B⎧0⎪⎪1⎪⎪2⎪3⎪⎪4⎪⎪5⎪
10n + B = ⎨6⎪⎪7⎪8⎪⎪9⎪⎪10⎪⎪11⎪12⎩

Eth = 0.0000 Q = 11.453
n'+10BEth = 7.9077 Q = −7.1840
p +10Be
Eth = 0.0000 Q = 0.225
α + 7LiEth = 0.0000 Q = 2.790
5He+ 6Li
Eth = 5.8944 Q = -5.3550
d + 9Be
Eth = 4.8014 Q = −4.362
t + 8Be
Eth = 0.0000 Q = 0.230
2n + 9B
Eth = 9.2858 Q = −8.436(np + pn )+ 9Be Eth = 7.2494 Q = −6.586(nα + αn )+ 6Li Eth = 4.9082 Q = −4.459(nd + dn )+8Be Eth = 6.6341 Q = −6.027( pα + αp )+ 6He Eth = 7.9077 Q = −7.184
2α + tEth = 0.0000 Q = 0.323
γ +11B

With 5He emission

No 5Heemission

Nuclear model calculations(1) Study of Nucleus–nucleus microscopic opticalpotential and applications(2) Study of Folding model and application

Optical potential without free parameter
Application to unstable nuclei

(1) Study of Nucleus–nucleus microscopic
optical potential and applications
Ma Yin-Qun and Ma Zhong-Yu
China Institute of Atomic Energy

Outline
1.The isospin dependent MOP of the nucleonnucleus scatterings have been obtained in theframework of the DBHF approach, where thelocal density approximation is used.
ieff

VM i 1i2i
2= U + Us +[U s − (U 0 − Vc ) ]E2Ei0

2. Extending this MOP to the case for a composite
projectile based on folding modelv
VFM ( R) =1r i r r
ρ (ki )Veff (ki )dki
3 ∑ ∫ i(2π ) i = p ,nIn the framework of momentum represensitivety In order to compare with experimental data, twonormalization factor are introduced. Then, theoptical potential is written as:rrrr
U opt ( R) = N R ReVFM ( R) + N I ImVFM ( R) + UC ( R)

6He+ 12C,

elastic scattering6He+ Be
9Ye Y L,PRC71 014604

6He+ Bi209

6 Li

12 C, 28 Si, 40 Ca, 58 Ni, 90 Zr, 208 Pb

30 – 318 MeVneutron-rich nuclei
11 Be
12 C,
11 Li
12 C, 28 Si,
17 F
12 C, 14 Ni,4 He

12 C,16O,542.3 MeV
319, 637 MeV
170 MeV40 Ca, 58 Ni, 90 Zr, 208 Pb

17F+ C
1217F + Ni
14

α +12C

α + 16 O

α + Zr90(2) Study of Folding Modeland Application
Nucleon-Nucleus interaction

V ( E , r ) = Vd ( E , r ) + Vex ( E , r )Direct term

Exchange term

Khoa: Phys. ReV. C 63, 034007 (2001)

Isospin dependent interaction
CDM3Y{+ ∫ dr′{⎡⎣ρ ( r, r′) + ρ ( r, r′) ⎤⎦ v ( ρ, E, s)± ⎡⎣ρ ( r, r′) − ρ ( r, r′) ⎤⎦ v ( ρ, E, s)} exp( ik ( E, r ) s)}

V (E, r) = ∫ dr′ ⎡⎣ρP ( r′) + ρn ( r′) ⎤⎦ v00d ( ρ, E, s) ± ⎡⎣ρP ( r′) − ρn ( r′) ⎤⎦ v01d ( ρ, E, s)
PP

ex00nnex
01

One-bodyDensity matrix

One-body density Matrix
Campi and Bouyssy: “Local approximation”

ρ ( R , R + s)

ρ (R + s 2) ˆj1 ⎡⎣ k F ( R + s 2 ) s ⎤⎦

ˆj1 ( x ) = 3 ( sin x − x cos x ) x3

Calculation of one – body desity matrix
based on shell modelA

ρ ( r , r ′ ) = ∑ ∑ ϕ i* ( r , σ ) ϕ i ( r ′, σ )i =1σWoods-Saxon potential

Comparison of folding potential bydifferent one-body density matrix
03253MeV p+ S

V/MeV-10-20PRESENT

LA-30-40(V exact-V)/V exact*100-50
50

12

34

5678

0

10%

-5-10-150

1234R/fm

5678

Optical potential
Normalization factor

U ( E, r ) = Nr[V ( E, r )] + iW ( E, r ) +Vso ( E, r )( l ⋅σ ) +VC
Folding model

Phenomenal
Koning-Delarochepotential

Folding vs JLM
40Ca(p, p)

Folding vs KD

Folding vs JLM

Folding vs KD

58 Ni (n, n)

90Zr (p, p)

7.76 -- 61.4 MeV

Folding vs JLM

Thank you !!!

Structural materials (33)Comparing to CENDL-2, the maindevelopment of CENDL-3.0 for thestructured materials is that in addition to
the data of natural elements, the data of
their isotopes were also included, and madethe data consistent between the naturalelements and their isotopes.

90ZrThe evaluated total cross section of n + 90

Fission products (115)most fission products nuclidesFrom MF-1 to MF-5 were provided in
CENDL-3.0 for, MF-1 to MF-6 are variable
for others.

101 fission products nuclides were sent
to join in the international comparison ofFP and coordinated by WPEC Subgroup21, and 38 of them were selected as thedata file of release 7 of ENDF/B-VI.

140Ce(n,2n) reaction cross section of n + 140

Sn (n, γ) capture cross section

Actinides (18)15 actinides were evaluated or reevaluated.
MF-1~6, 12~15 were included forimportant actinides (i.e. U, Pu isotopes),and MF-1~5 for others.
The results of the benchmark testing
havebeen
consideredduringtheevaluation process.

240Pu (n, tot)

Fig.8
Fig.8

240Pu (n, f)

240Pu (n, inl) (MT=51)

1.3 The benchmark testing of CENDL-3.0
In order to test the reliability of the data from CENDL-3.0,
benchmark testing for some important nuclides of CENDL-3.0
has been done, the calculations and analyses of benchmarks
were done with Monte Carlo code MCNP and transformation
codes. The data processing was carried out by using theinternationally used code system NJOY97. The results werecompared with ENDF, JENDL and JEFF validated CENDL-3.0
to identify the source of the discrepancies with theexperimental results.
Table 1 and 2 show the results of calculation with different
233
233U and 238
238U as the examples.
libraries Keffofeff
Some
Some ofof results
results based
based onon benchmark
benchmark FNS
FNS andand OKTAVIANOKTAVIAN forfor
light
light andand medium-heavy
medium-heavy nuclides
nuclides ofof CENDL-3.1
CENDL-3.1 areare presented
presented
here
here (Fig.11-13).
(Fig.11-13).

Keff for 233U
Assembly

Experiment

CENDL-3.0

ENDF/B-VI

JENDL-3.2

233U Jezebel

1.000(～0.001)

0.998220.99235

1.01460

233U-F-002a

1.000(～0.001)

0.99518

0.99544

1.00799

233U-F-002b

1.000(～0.001)

0.995830.999431.00959

233U-F-003a

1.000(～0.001)

0.99474

0.99755

1.01168

233U-F-003b

1.000(～0.001)

0.995230.99989

1.00929

233U-F-004a

1.000(～0.001)

1.00715

1.00494

1.01747

233U-F-004b

1.000(～0.001)

1.00808

1.00940

1.01721233U-F005a

1.000(～0.001)

0.99951

0.994091.00852

233U-F-005b

1.000(～0.001)

1.00204

0.99803

1.00763

Flattop-23

1.000(～0.001)

0.992181.00146

1.01140

Averaging

1.000(～0.001)

0.998820.99926

1.01154

Keff for 238U
Benchmarks of Thermal Reactor
Assembly

CENDL-3.0

JEF-2

JENDL-3.2

ENDF/B-6

TRX-1

0.99750.99520.99010.9908TRX-2

0.99980.9972

0.99200.9924

BAPL-UO2-1

1.00101.00200.99750.9952

BAPL-UO2-2

1.00031.0014

0.99720.9951

BAPL-UO2-3

1.00071.0014

0.9976

0.9960Conclusion
1. CENDL-3.0 was finished in 2000, and it has been improved
since then. CENDL-3.1 is be preparing for its release.
2. Nuclides and data files of CENDL-3.1 were increased andextended compared with CENDL-2.1.
3. All evaluations performed based on the new experimental
data and model calculations carefully. They were improved
compared with CENDL-2.1.
4. The data of the most important nuclides were validated bythe benchmark testing. Some of them are better than otherevaluated libraries according to the results of thebenchmark testing .

2.4 Nuclear structure and decay data
CNDC has taken permanent responsibility for evaluating andupdating NSDD for A=51, and 195-198 mass chain. The data
have been revised using available experimental decay andreaction data for mass chain A=197, and are being updated formass chain A=196. Updated evaluation of A=197 has been sentto NNDC, USA
USA andand was
was publishedpublished in NDS in 2004. The
evaluations of mass chain A=52-56 were being updated at Jilin
233U were being evaluated on theUniversity. The decay data of 233
basis of the new measured data.

2.3 The systematic study of fission yield data
Based on the mass distribution data up to 200 MeV measured
by Zoller, the systematic on dependence of chain yield on
incident neutron energy for each mass number A was studied.
And also the systematics of mass distribution on mass A andincident neutron energy was investigated by using 5 (or 3)
Gaussian model. The calculated results could reproduce theexperimental data used well. The investigation also shows thatthe correlation between the parameters of the systematic andthe yields calculated with the systematics is quite complicatedand, in general, is quite strong.

Reaction CrossSection

CDM3Y nucleon-nucleon Interactionvd(ex )
00(01)v( ρ , E, s ) = F ( ρ ) g ( E ) v

d(ex )
00(01)

3d ( ex )0 0 ( 0 1)(s ) = ∑ Y
υ =1d ( ex )0 0 ( 0 1)

ex p ( − Rυ s )Rυ s(s)Isospin dependent interaction{+ ∫ dr′{⎡⎣ρ ( r, r′) + ρ ( r, r′) ⎤⎦ v ( ρ, E, s)± ⎡⎣ρ ( r, r′) − ρ ( r, r′) ⎤⎦ v ( ρ, E, s)} exp( ik ( E, r ) s)}

dd′′′′′V (E, r) = ∫ dr ⎡⎣ρP ( r ) + ρn ( r ) ⎤⎦ v00 ( ρ, E, s) ± ⎡⎣ρP ( r ) − ρn ( r ) ⎤⎦ v01 ( ρ, E, s)
PPex
00nnex01

Code:
LUNF :

Light elements

SUNF :

Structural Materials, fission products
Medium elements
Heavy elements & actinides

FUNF :

APMN/APOM94:
DWUCK + ECIS:INPUT:

spherical optimal potential
inelastic / deformed nuclei
RIPL + adjustNeutron spectrum for the OKTAVIAN Ti

TiTi

α+208Pb

Plan
2006 --- 2010
- More than 100 nuclides will be evaluated.
- Revision and validation.
- covariance data.

CDM3Y nucleon-nucleon Interactionv

d(ex )
00(01)v( ρ , E, s ) = F ( ρ ) g ( E ) v

d(ex )
00(01)

3d ( ex )0 0 ( 0 1)(s ) = ∑ Y
υ =1d ( ex )0 0 ( 0 1)

ex p ( − Rυ s )Rυ s(s)Comparison of elastic angular distribution
41032

p+ S, E =53 MeVlab3(dσ/dΩ)/(mb/sr)

10Khan et al.
PresentFold-LA

2101

100

100

102030Θ40o/( )

c.m.

5060

II．Progresses on nuclear data evaluations
2.1 Nuclear reaction model study
The model was improved and completed for 1p shell light
nuclides, which contains the dynamics and kinematics ofnuclear reactions.
A method to set up file-6 of light nuclei for evaluated neutron
data in ENDF/B-6 format below 20 MeV has been established andthe energy balance was strictly considered. This method has12C.
been used in the calculation of n +12
The possibility of 55He emission has been investigated in the lightnuclear reactions, and the formulation for calculating of the 55Heemission was developed. Also the double differential crosssections of 55He emission, as well as the spectra of neutron andalpha particle from the breakup of 55He, were set up and used in
the series codes of LUNF.

(2) Study of CDCC approach andapplication
An Haixia and Cai Chonghai
NanKai University

outline A new code is developed based on CDCC

approach and P3C5 algorithm.
 Contribution of break-up process is considered.

 Applied to elastic scattering by Deutron incident

on target ranging from 12 C to 208 Pb

Deutron elastic
angular distribution
Compared toSpherical Optical
Model calc. withglobal potential

Deutron elastic
angular distribution

Application to elastic scattering
for stable and unstable nuclei
1.

6He incident on 12C, 9Be, 208Pb, 209Bi

22.5 –229.8 MeV