250 likes | 425 Views
Consistent nuclear model calculations for deuteron induced reactions on Al and Cu. FENDL-3 2nd RC M NDS/IAEA, Vienna, 23-26 March 2010. Marilena Avrigeanu and Vlad Avrigeanu. http://tandem.nipne.ro/~vavrig/ , http://fp6.cordis.lu/fp6/partners /RCN 49105.
E N D
Consistent nuclear model calculations for deuteron induced reactions on Al and Cu FENDL-3 2nd RC M NDS/IAEA, Vienna, 23-26 March 2010 Marilena Avrigeanu and Vlad Avrigeanu • http://tandem.nipne.ro/~vavrig/, http://fp6.cordis.lu/fp6/partners/RCN 49105 "Horia Hulubei" National Institute for Physics and Nuclear Engineering, P.O. Box MG-6, Bucharest, Romania • Motivation • Elastic Scattering (OMP) • Breakup • Pre-equilibrium • Direct reactions • Evaporation • Activation cross sections calculations 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
IFMIF (International Fusion Material Irradiation Facility) D-Li neutron source for testing fusion reactor candidate materials Motivation: Data Needs • General purpose : Activation/transmutation data library • [U. Fischer: "Nuclear Data Libraries for Advanced Systems: Fusion Devices“, Nov 2007, IAEA, Vienna] • Urgent need for qualified IFMIF reference data library • IEAF/EAF : (n, p, d), to be further developed, improved and validated • Cross section measurements required for neutrons (E >20 MeV): Cr, Co, V, W, Ta, Pb, Bi, Au, Mn deuterons (E < 40 MeV): Cu, Al, Nb protons (E < 12 MeV): Cu, Al, Nb, Ta, W, Au, Pb • Evaluations (n, p, d) to be further developed, improved and validated • Reliable gas production cross-section data (H, He) • Dosimetry data file to be developed for E > 20 MeV (IRDF) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
Nuclear Model Calculations • Pure elastic scattering OP analysis SCAT2 [O. Bersillon] - phenomenological OP + semi-microscopic (DF) OP (local version) • DFOLD [M. Avrigeanu] • - double folding method Optical Model Potential M. Avrigeanu et al., P.R.C62 (2000) 017001; N.P.A693 (2001) 616; E.PJ.A12 (2001) 399 ; N.P.A723 (2003) 104; N.P.A 764 (2006) 246. • Direct reactions -breakup FRESCO-2003 [I.J. Thompson] - elastic transfer among weakly bound systems -direct reactions,(d,p), (d,n), (d,t) • Composite system equilibration • STAPRE-H95 [V. Avrigeanu, M. Avrigeanu] (updated) • - OMP:SCAT2000; preequilibrium:GDH / EXCITON;evaporation:Hauser-Feshbach • TALYS - 1.0 [A. Koning, S. Hilaire, M. Duijvestijn] • - OMP:ECIS’97;preequilibrium:MSD / EXCITON;evaporation: Hauser-Feshbach 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
NO GLOBAL OPTICAL MODEL POTENTIAL (OMP) for d + Nucleus (A<27) • COMPARATIVE ANALYSISofglobal OMPs for d + 27Al, 63,65Cu, 54,56,58Fe, 93Nb: • Lohr-Haeberli (1974): A~40-209, E=8-13 MeV • Perey-Perey (1963,1976): A~40-208, E=12-25 MeV • Daehnick et al. (1980): A~27-238, E=11.8-90 MeV • Bojowald et al. (1988): 27Al, 89Y, 120Sn, and 208Pb at Ed=58.7 and 85 MeV • None of these global OMP describes data at E<15 MeV • Semi-microscopic OMP byusing realistic nucleon-nucleon interaction UDF, & more accurate parameterization for WV, WD, VSO • Phenomenological OMP: parameterization forUR (frozen W&VSO) • Cross-Sections calculations Optical Model Potential 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
Microscopic Real Optical Potential M. Avrigeanu et al., Proc. Int. Conf. on Nuclear Reaction Mechanisms, Varenna, 2006, pp.193 Approximations 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
DF-real &phenomenological imaginary/s.o. potentialsPhenomenological OMP M. Avrigeanu et al., Fusion Eng, Des. 84, 418 (2009) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
OMP Analysis for 63,65,natCu(d,d0) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 63,65Cu:activation cross sectionsfirst step: d-OMP 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d breakup involvement : d+27Al 29Si* elastic 27Al + d 27Al + n+ p inelastic(27Al + n)+p 28Al*+p 28Al + γ+p(γpch. of29Si*) 27Al + n+p(npch. of29Si*) 27Mg + p+p(2pch. of29Si*)24Na + α+p(αpch. of29Si*) inelastic(27Al +p)+n 28Si*+n 28Si + γ+n(γnch. of29Si*) 27Si + n+n(2nch. of29Si*) 27Al + p+n(pnch. of29Si*) 24Mg + α+n(αnch. of29Si*) σBF * σnx /σnTotal σBF * σpx /σpTotal 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d breakup involvement : d+63,65Cu 65,67Zn* elastic 63,65Cu + d 63,65Cu + n+ p inelastic(63,65Cu + n)+p 64,66Cu*+p 64,66Cu + γ+p(γpch. of65,67Zn*) 63,65Cu + n+p(npch. of65,67Zn*) 63,65Ni + p+p(2pch. of65,67Zn*)60,62Co + α+p(αpch. of65,67Zn*) inelastic(63,65Cu +p)+n 64,66Zn*+n 64,66Zn + γ+n(γnch. of65,67Zn*) 63,65Zn + n+n(2nch. of65,67Zn*) 63,65Cu + p+n(pnch. of65,67Zn*) 60,62Ni + α+n(αpch. of65,67Zn*) σBF * σn,x /σnTotal σBF * σp,x /σpTotal 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
1 (a)fpTBU= 0.087 – 0.0066 Z + 0.00163 Z A1/3 + 0.0017 A1/3E – 2E-6 ZE2 • 1 (b)fELBU= 0.031 – 0.0028 Z + 0.00051 Z A1/3 + 0.0005 A1/3E – 1E-6 ZE2 Empirical breakupcomponentsfBUsystematics Q{Al(d,p)} = 5.5 MeV M. Avrigeanu et al., Fusion Eng, Des. 84, 418 (2009) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
Deuteron breakup components 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
BREAKUP contribution to the SECOND particle emission P. Bem,…,M.A.,..Phys. Rev. C 79, 044610 (2009) σBF(p)vsEdσBF(p)vsEdσBF(n)vsEd 27Al(n,p)27MgvsEn27Al(n,α)24NavsEn 27Al(p,n)27SivsEp 27Al(n,p)27MgvsEd 27Al(n,α)24NavsEd27Al(p,n)27SivsEd [via BF] [via BF] [via BF] 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
OMP:deuteron channel: present work neutron / proton channels:Koning-Delaroche (2003) • neutron-proton interaction:<φ(d),φ(p/n)>;V=V0e-(r/rd)2 V0=72.15 MeV; rd=1.484 fm • nucleon bound state:<φ(Target),φ(Residual)>;real-WS V0;r0=1.25; a=0.65 • spectroscopic factors, SIJlsj: experimentalproton/neutronangular distributions ΨJM(ξ,r)~∑lsjAIJlsj[φItarget(ξ) φlsj]JM; |AIJlsj|2 =SIJlsj • 35 levels up to 5.135 MeV for the odd-odd 28Al( 27Al(d,p)28Al) • 24levels up to 11.4 MeV for the even-even 28Si( 27Al(d,n)28Si) • 63 levels up to 3.030 MeV for the odd-odd 64Cu( 63Cu(d,p)64Cu) • 52 levels up to 3.080 MeV for the odd-odd 66Cu( 65Cu(d,p)66Cu) • 119levels up to 6 MeV for the odd-odd 59Co (59Co (d,p) 60Co) STRIPPING calculations details (FRESCO) NO Breakup ? NO Stripping ? ONLY STATISTICAL Hauser-Feshbach Ep (MeV) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
Stripping contribution (FRESCO-code; 35 levels of 28Al) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
Stripping contribution (FRESCO-code; 18 levels of 28Si) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 27Al:activation cross-sections P.Bem,…,M.A.,..Phys. Rev. C 79, 044610 (2009) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 63,65Cu:activation cross-sections M. Avrigeanu, V. Avrigeanu, Eur. Phys. J . Web of Conf. 2, 01004 (2010) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 63,65Cu:STAPRE-HEMPIRETALYS 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 93Nb:Activation cross-sections (preliminary results) 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 93Nb:breakup option inTALYS 1.2 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
d + 93Nb:breakup option inTALYS 1.0 & 1.2 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
CONCLUSIONS • Semi-microscopic OMP analysis • UDF:ρd (charge) & ρAl (charge) & M3Y Paris-NN • WD & VSO phenomenological • improved agreement with data adding the dispersion corrections • Phenomenological OMP analysis for27Al, 63.65,natCu, 93Nb • agreement with all available measured data • good description of (d,d) data vs. TALYS defaultOMPs • improved description of (d,d) data vs. globalOMPs • Completion of d+27Al,63,65Cu activation cross sections • BU, BF - Deuteron break-up mechanism contributions, • DR mechanism considered trough n & p stripping • PEand evaporation mechanisms contributions • Comparison ofSTAPRE-H / TALYS-1.0activation cross sections • Preliminary results for d+93Nb activation cross sections involving TALYS 1.2 • NEED to UPDATE the d-activation libraries 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010
Thank you ! 2nd RCM-FENDL-3, IAEA Vienna, 23-26 March 2010