A systematic TDDFT data for nuclear fission analysis 岩田順敬 (Yoritaka Iwata) 東工大 科学技術創成研究院 (Tokyo Tech) with T. Nishikawa (Nuclear Engineering, Ltd) K. Kean (Tokyo Tech. and JAEA) S. Chiba (Tokyo Tech.) Introduction on nuclear chart Fission mechanism is studied for this region Nucleus in our scope Thermal Neutron Fission CS We have less information about the fissility BE/A Halflife 236 U Motivation ◇ Nuclear fission is one of theoreticallytheoretically difficult processes in sub-atomic physics. DifficultDifficult toto treattreat microscopicallymicroscopically HereHere wewe taketake lessless phenomenologicalphenomenological (ab-initio(ab-initio like)like) treatmenttreatment suchsuch asas TDDFTTDDFT (time-dependent(time-dependent densitydensity functionalfunctional theory).theory). ◇◇ InIn termsterms ofof fission,fission, twotwo topicstopics areare focusedfocused onon …… __ many-bodymany-body forceforce __ frictionfriction andand dissipationdissipation Fission dynamics mapped on Eq. of Motion Γ(t) :random force macroscopic Potential V Friction γ … by CHF … by TDDFT microscopic TDDFT calc. R Contents ◇ CHF for FissionFission barrierbarrier heightheight SLy4SLy4 byby Kean-Iwata-ChibaKean-Iwata-Chiba …… effectiveeffective many-bodymany-body forceforce effecteffect ◇ TDDFT + Langevin model for Friction by Nishikawa-Iwata-Chiba …… effectiveeffective dissipationdissipation includedincluded inin TDDFTTDDFT SV-basSV-bas ◇ Result Theoretical Nuclear data … provided by theoretical systematics ◇ Summary Density functional theory WHY microscopic ? ➔ Nucleon degree of freedom （ Fermi statistics ） ➔ Nucleon-nucleon interaction （ incl. spin-orbit force ） ➔Mean-field type description ➔ Coulomb force nucleus ➔ Less phenomenological Less free parameter (no free parameter for dynamics) ーーーーーー ➔TDDFT and CHF use the same functional form proton neutron Fission dynamics to refine fissility ... mapped on Eq. of Motion Γ(t) :random force macroscopic Potential V Friction γ … by CHF … by TDDFT FissilityFissility microscopic TDDFT calc. R ~ 50 ◇ CHF for FissionFission barrierbarrier heightheight FissionFission barrierbarrier heightheight byby constraintconstraint Hartree-FockHartree-Fock withwith pairingpairing (density(density dependentdependent pairing)pairing) 15 236UU ] ] V V 10 e e M M α = 1/6 [ [ E E 5 SLy4SLy4 ββ=Q=Q20 0 0 0.2 0.4 0.6 0.8 1 1.2 SinceSince thethe barrierbarrier heightheight isis examinedexamined byby β,β, onlyonly symmetricsymmetric fissionfission cancan bebe identifiedidentified ([Thr.]([Thr.] barrierbarrier heightheight becomesbecomes ratherrather higher).higher). Fission barrier height We need a reliable fission barrier inner barrier (outer barrier) Fission barrier height is 4 to 7 MeV (exp.), while they are calculated from 7 to 14 MeV (SLy4). … almost twice larger The improved interaction is necessary ! _ why are most of the Skyrme interactions so large ? _ is there any way to re-fit the Skyrme interaction ? in the corresponding part of nuclear chart Medium effect force effect Hamiltonian density of density functional theory Analogous to Alpha term: density dependent term on the fission barrier height ... Perturbation analysis for alpha-effect 15 236UU 10 α = 1/6 5 alpha value +10 % SLy4SLy4 ββ=Q=Q20 0 20 0 0.2 0.4 0.6 0.8 1 1.2 [Perturbation analysis] -10 % We examined 10 % different cases starting from Original alpha = 1/6 = 0.1666... Sensitivity of fission barrier height on alpha We have confirmed Only 10 % change of a … can contribute to … ~40% difference of fission barrier height This fact is confirmed in U-isotopes 231 to 239 Pu-isotopes 237 to 245 using SkM*, SLy4, SV-bas interactions. alpha value -10 % Alpha-perturbation theory This fact is our motivation to change alpha slightly to have a better description of fission barrier height What to change [Strategy] +10 % Also change this one to keep the total contribution as a whole term Validity of the proposed method Binding energy [MeV] Original calc (Proposed calc) Original Proposed Radius [fm] RIPL3 r_p, r_n, r_tot (r_p, r_n, r_tot) Only small deference in these values (quality is kept) Correlation between experimental and theoretical results for SLy4 and SkM* interactions ● original ● proposed Theoretical database I using substantially refined Skyrme forces …m-SLy4 m-SkM* m-SV-bas in terms of fission ◇ TDDFT + Langevin model Friction along the collision axis Process of extracting friction effect Model time evolution Duration time ~ 1.5 [10-21 s] 296116 = 296Lv symmetric fission by TDDFT R Y.Iwata and S. Heinz, CERN Proceedings 002 (2012) Averaged friction coefficient Microscopic inputs R, dR/dt : extracted ◇ TDDFT + Langevin model It should play a essential role in fission Friction along the collision axis Validity for TDDFT + 4D Langevin method 236UU This is good enough (more than expected), now we engage to calculate the other cases Linearly scaled energy-dependence ← Root-mean squared line This is not only for the U236, but for the other isotopes. Indeed ... [systematics] Linearly scaled energy-dependence N=150 120 Uranium isotopes 100 80 Friction coeff. g 60 40 N=146 120 20 100 0 Energy 1 2 3 4 5 6 807 8 9 10 11 60 40 N=142 20 120 0 100 1 2 3 4 5 6 7 8 9 10 11 80 60 40 20 0 1 2 3 4 5 6 7 8 9 10 11 [systematics] Linearly scaled energy-dependence N=152 140 Plutonium isotopes 120 100 80 Friction coeff. g 60 40 N=148 20 120 0 100 1 2 3 4 5 6 7 8 9 10 11 Energy 80 60 40 N=144 20 140 0 120 1 2 3 4 5 6 7 8 9 10 11 100 80 60 40 20 0 1 2 3 4 5 6 7 8 9 10 11 [systematics] Linearly scaled energy-dependence N=154 120 100 Curium isotopes 80 Friction coeff. g 60 N=150 40 120 20 100 0 2 3 4 5 6 807 8 9 10 11 Energy 60 40 20 N=146 0 120 2 3 4 5 6 7 8 9 10 11 100 80 60 40 20 0 2 3 4 5 6 7 8 9 10 11 [Systematics] Linearly scaled energy-dependence Modeled by ... γ = a (E/A) + b b 115 120 125 130 135 140 145 -4.5 -5 U Pu [As a physical law] [As a physical law] Cm Cf This would be the first time -5.5 Fm to show the energy a dependence of friction -6 using microscopic calculations -6.5 Statistical law found by Nishikawa-san -7 Neutron number Energy dependence is rather moderate compared to mass dependence in which Z dependence is larger than N dependence Friction along the collision axis 236U Similarity can be found in BE How does the friction behave ? Thermal Neutron Fission CS Similarity can be found in BE/A BE/A Halflife 236 U Theoretical Database II (SV-bas and more) SkM* , SLy4 calcs will follow before long. SLy4 0.58 100 98 0.54 96 A 94 / γ 0.5 0 92 0.46 1.48 1.52 1.56 1.6 1.64 N/Z The model eq. for exploring Z=120 Stochastic 1D eq. to refine the fissility Γ(t) :random force + Coulomb pot Potential V Friction γ … by CHF … by TDDFT microscopic TDDFT calc. R Summary SkM*SkM* ◇ FissionFission barrier/Potentialbarrier/Potential SLy4SLy4 (mainly(mainly byby K.K. KeanKean sansan)) … Alpha-perturbation theory SV-basSV-bas ◇ Friction (mainly by Nishikawa san) … Friction formula (depending on mass and energy) These are substituted into both _ 1D Langevin type equation (for refining fissility) _ Usual 3D and 4D Langevin equation .