Montecarlo simulation of the decay of warm superdeformed nuclei

Similar documents
Warm superdeformed nuclei:

Study of shell structure and order-to-chaos transition in warm rotating nuclei with the radioactive beams of SPES

Rotational motion in thermally excited nuclei. S. Leoni and A. Bracco

4. Rotational motion in thermally excited nuclei *

GDR FEEDING OF THE HIGHLY-DEFORMED BAND IN 42 Ca

Damping mechanisms and order-to-chaos transition in the warm rotating 163 Er nucleus

The rotational γ -continuum in the mass region A 110

Structure of hot nuclear states

Nuclear Structure at Extreme conditions through γ spectroscopy measurements

COLLECTIVITY IN LIGHT NUCLEI AND THE GDR

Isospin Character of Low-Lying Pygmy Dipole States via Inelastic Scattering of 17 O ions

p(p,e + ν)d reaction determines the lifetime of the sun

Recent Results from Ridge Analysis- Evidence for Hyperdeformation?

Investigation of Pygmy Dipole Resonance in neutron rich exotic nuclei

Evidence for the wobbling mode in nuclei

Order-to-Chaos Transition in Rotational Nuclei

Silvia M. Lenzi University of Padova and INFN. Silvia Lenzi, 10th Int. Spring Seminar on Nuclear Physics, Vietri sul Mare, May 21-25, 2010

Heavy-ion sub-barrier fusion reactions: a sensitive tool to probe nuclear structure

Evolution Of Shell Structure, Shapes & Collective Modes. Dario Vretenar

Collective Excitations in Exotic Nuclei

High-spin studies and nuclear structure in three semi-magic regions of the nuclide chart High-seniority states in Sn isotopes

Probing neutron-rich isotopes around doubly closed-shell 132 Sn and doubly mid-shell 170 Dy by combined β-γ and isomer spectroscopy.

Coexisting normal and triaxial superdeformed structures in 165 Lu

Gamma-ray spectroscopy I

Structure of halo nuclei and transfer reactions

Coherent and incoherent nuclear pion photoproduction

New Trends in the Nuclear Shell Structure O. Sorlin GANIL Caen

Giant Resonances Wavelets, Scales and Level Densities

One- and two-phonon wobbling excitations in triaxial 165 Lu

Evidence for the inuence of reaction dynamics on the population of compound nuclei

Fusion of light halo nuclei

Submitted to Acta Physica Polonica, Proceedings of the XXXI Zakopane School of Physics

Monte Carlo Simulation for Statistical Decay of Compound Nucleus

Photopion photoproduction and neutron radii

The many facets of breakup reactions with exotic beams

Beyond mean-field study on collective vibrations and beta-decay

The two-step (and multiple-step) γ cascade method as a tool for studying γ-ray strength functions. Milan Krtička

Angular-Momentum Projected Potential Energy Surfaces Based on a Combined Method. Jianzhong Gu. (China Institute of Atomic Energy, Beijing, China)

Fission fragment mass distributions via prompt γ -ray spectroscopy

Nuclear Isomerism. Phil Walker. University of Surrey. on the occasion of the 70th birthday of Geirr Sletten

Fitting Function for Experimental Energy Ordered Spectra in Nuclear Continuum Studies

Sunday Monday Thursday. Friday

Chapter 6. Summary and Conclusions

arxiv:nucl-th/ v1 28 Sep 1995

K-isomer decay rates. Phil Walker. - introduction - exploring the different degrees of freedom - exotic nuclei

Shape Coexistence and Band Termination in Doubly Magic Nucleus 40 Ca

Physics of neutron-rich nuclei

Microscopic description of fission in the neutron-deficient Pb region

Radiative Capture Reaction

Low-spin structure of 210 Bi

Shape coexistence in light Krypton isotopes

Cranked-RPA description of wobbling motion in triaxially deformed nuclei

Microscopic Fusion Dynamics Based on TDHF

SUB-BARRIER FUSION of HEAVY IONS

Dipole Response of Exotic Nuclei and Symmetry Energy Experiments at the LAND R 3 B Setup

The Nuclear Many-Body Problem

arxiv:nucl-ex/ v1 9 Apr 2001

Probing shell evolution with large scale shell model calculations

Giant Dipole Resonance - New Experimental Perspectives

Heavy-ion fusion reactions for superheavy elements Kouichi Hagino

Subbarrier fusion of carbon isotopes ~ from resonance structure to fusion oscillations ~

Spectroscopic Quadrupole Moment in 96,98 Sr : Shape coexistence at N=60. E.Clément-GANIL IS451 Collaboration

arxiv:nucl-ex/ v1 27 May 2002

Medium polarization effects and pairing interaction in finite nuclei

Structure of neutron-rich Mg isotopes explored by beta-decay of spin-polarized Na isotopes

Shell model description of dipole strength at low energy

Nuclear structure aspects of Schiff Moments. N.Auerbach Tel Aviv University and MSU

SUPPLEMENTARY INFORMATION

Spin Cut-off Parameter of Nuclear Level Density and Effective Moment of Inertia

Oblate nuclear shapes and shape coexistence in neutron-deficient rare earth isotopes

Study of multi-nucleon transfer reactions with light nuclei

DSAM lifetime measurements for the chiral bands in 194 Tl

How to do C.C. calculations if there is only limited experimental information on intrinsic degrees of freedom?

Nonstatistical fluctuations for deep inelastic processes

Application of quantum number projection method to tetrahedral shape and high-spin states in nuclei

Microscopic description of shape coexistence and shape transition

Linking nuclear reactions and nuclear structure to on the way to the drip lines

Particle excitations and rotational modes in nuclei with A 70-90

First results with PARIS array

Proximity Decay and Tidal Effects

A Comparison between Channel Selections in Heavy Ion Reactions

Testing the validity of the Spin-orbit force Nuclear forces at the drip-line O. Sorlin (GANIL, Caen, France) PART 1:

Stability of heavy elements against alpha and cluster radioactivity

Symmetries and collective Nuclear excitations PRESENT AND FUTURE EXOTICS IN NUCLEAR PHYSICS In honor of Geirr Sletten at his 70 th birthday

(Multi-)nucleon transfer in the reactions 16 O, 3 32 S Pb

Observation and Quadrupole Moment Measurement of the First Superdeformed Band in the A~60 Mass Region

Compound and heavy-ion reactions

New theoretical insights on the physics of compound nuclei from laser-nucleus reactions

Pygmy dipole resonances in stable and unstable nuclei

Studying the nuclear pairing force through. Zack Elledge and Dr. Gregory Christian

Collective excitations of Λ hypernuclei

Mapping Fission in Terra Incognita in the neutron-deficient lead region

Radioactivity at the limits of nuclear existence

Conversion Electron Spectroscopy in the Second Minimum of Actinides

Report on the benchmarking of the event generator for fusion-evaporation reactions

Mapping Low-Energy Fission with RIBs (in the lead region)

Nature of low-energy dipole states in exotic nuclei

Spectroscopy of fission fragments using prompt-delayed coincidence technique

First results from the AGATA Demonstrator. Francesco Recchia Università di Padova

Statistical-Model and Direct-Semidirect-Model Calculations of Neutron Radiative Capture Process

Transcription:

Montecarlo simulation of the decay of warm superdeformed nuclei E. Vigezzi INFN Milano Understanding the dynamics in the SD well: probing γ strength functions, energy barriers, level densities, residual interactions S. Leoni et al., PRL 11, 14252 (28); PRC; to be published

Two High Statistics EUROBALL-IV Experiments 151 Tb and 196 Pb HECTOR BGO INNER BALL 21 detectors 8 Large BaF 2 for GDR analysis 3 Si + 17 Er 196 Pb + 4n Thin target, E beam = 148 MeV, L max 53 Statistics: 8 days 2 1 8 <F Ge > = 3 <F > = 1 GDR in A = 19 (D. Montanari et al., in preparation) 27 Al + 13 Te 151 Tb + 6n Thin target, E beam = 155 MeV, L max 8 Statistics: 17 days 9 1 9 <F Ge > = 5 <F > = 22 discrete spectroscopy in SD well (J. Robin et al., PRC77, 1438(28)) search for linking transitions SD ND (J. Robin et al., PRC78, 34319(29))

Quasi-Continuum Spectroscopy: Warm Rotation RIDGES Discrete Unresolved Bands E* 168 Yb 168 Yb I-2 I I+2 T= T I-2 I I+2 E2 bump I+1 spectra I I+8 I+6 I+4 I+2 I

Compound Nucleus Chaos Warm Rotation strongly interacting bands I I-2 I+2 SD T= Regular bands Order L m ND ND SD tunneling SD to ND minima

151 Tb Analysis SD gated spectrum E2 Bump SD gated 1 st Ridge Rot-Plane 2 nd Ridge Total 1 st Ridge counts counts [x1 3 ] counts [x1 3 ] counts [x1 3 ] 3 25 2 15 1 5 75 125 175 15 1 5 2 15 1 E E = 1 kev 25 E = 2 kev 7 5 3 E = 1 kev -1 1 E 1 -E 2 [kev] Intensity [%] Intensity [%] Intensity [%] Intensity [%] Intensity 1 3 E2 SD gated 1 1 3 SD gated Ridge 1 1 3 2. Ridge 1 1 3 1 8 independent observables Yrast I yrast = 2 % 1. Ridge 3 4 5 6 7 N path N path N path Number of Discrete Bands Counts Fluctuations 75 5 25 1 75 5 25 1 75 5 25 N (2) path 14 ± 2 12 ± 5 27 ± 8 Neve P 2 1 1 SD gated Ridge 2. Ridge 1. Ridge 1 observed 2 3 4 5 6 7

196 Pb Analysis 8 independent observables SD gated spectrum E2 Bump SD gated 1 st Ridge Rot-Plane 2 nd Ridge Total 1 st Ridge counts [x1 3 ] counts [x1 3 ] counts [x1 3 ] Counts [x1 3 ] 8 6 4 2 8 6 4 6 5 4 15 1 5 5 75 1 125 E E = 8 kev E =16 kev E = 8 kev -8 8 E 1 -E 2 [kev] Intensity [%] Intensity [%] Intensity [%] Intensity [%] 1 3 1 1 3 1 1 3 1 1 3 1 Intensity E2 SD gated E2 + M1 SD gated Ridge 2. Ridge 1. Ridge Yrast I yrast = 1.3 % 3 4 5 N path N path N path 2 75 Number Discrete Bands Counts Fluctuations N (2) path Neve P 2 1 1 4 SD gated Ridge 5 25 75 5 25 2. Ridge 1. Ridge 3 4 5 8 ± 3 12 ± 4 3 observed 2 ± 9

Intensity [%] Intensity [%] Intensity [%] Intensity [%] Intensity 1 3 E2 SD gated 1 1 3 SD gated Ridge 1 1 3 2. Ridge 1 1 3 1 Yrast I yrast = 2 % 3 4 5 6 7 1. Ridge 151Tb 4% of the yrast intensity: The rest must come from damped bands 2-3 times higher than the yrast intensity: bands decay to ND well

at T Vnn = 34 MeV fm 3 Vnp= 5 MeV fm 3 H = H def - J x + V res Configuration mixing

E2 fragmentation rotational damping Width [kev] 35 3 25 2 15 1 5 rot I-2 I rot n branch 1 1 1 151 Tb n branch = 2 n branch = 2 onset of damping x 1/1 I = 5 + I = 4 + I-2 I 3 4 5 6 7 1 2 3 4 U [MeV] K. Yoshida, M. Matsuo, NPA612(1997)126

E2 ND E1 E2 SD E1 extrapolated SD, B(E2), barriers 2 E* [MeV] ND tunneling SD rot microscopic levels, B(E2), barriers 15 rot I-2 I I+2 SD 1 ND -decay in SD well: microscopic + extrapolated 3 4 5 6

Microscopic calculations for decay-out at T K. Yoshida, M. Matsuo and Y. Shimizu, NPA696, 85(21) ND SD M t 1 2 DSD 1 2 1 e action integral S ( E ) ds 2M ( V ( q( s)) E ) path along tunneling path 2S driven by pairing S S 2 15 1 5 12 1 151 Tb (,+) U =. MeV U =. MeV.4 196 Pb.2 (,+) 8.4.6 6.8 1. 4 1.2 1.4 2 1.6 1.8 3 4 5 6.8 1.6 2.4 2.8 3.6 4.4 5.2 extension of model used for SD yrast 152 Dy

Ingredients of the Monte Carlo simulation: Microscopic levels in SD well and actions: calculated up to 5 MeV above yrast (151 Tb) and 3 MeV (196Pb); extrapolated above. Inertial mass multiplied by.7 (151Tb) and 2.5 (196Pb) E2 transitions: calculated in the SD well; schematic parametrization in the ND well E1 transitions: GDR tail, with a scaling factor to reproduce yrast intensity Entry distribution: calculated Level density in the ND well: P. Demetriou and S. Goriely, NPA 695 (21) 95

151 Tb: simulated -spectra Counts [arb. un.] counts [x1 3 ] 6 4 2 7 5 3 Total matrix E = 1 kev E = 1 kev SIMUL DATA -1 1 E 1 -E 2 [kev] Counts [arb. un.] counts [x1 3 ] 15 1 5 15 1 5 SD gated matrix E = 1 kev E = 1 kev -1 1 E 1 -E 2 [kev] counts counts 15 1 5 375 125 175 25 2 15 1 5 E2 SD gated 75 125 175 E SIMUL DATA

Exp. versus Theory Analysis of Intensities 151 Tb 196 Pb exclusive quantities SD gated decay-flow SD gated γ-spectrum E2 Bump SD gated γ-γ 1st Ridge inclusive quantities average decay-flow Rot-Plane 2nd Ridge Total γ-γ 1st Ridge

SD gated quantities are sensitive to E1/E2 balance Enhanced dipole transitions have been observed in A= 15 and A = 19 196 Pb, 19 Hg 194 Hg D. Rosbach et al. PLB513(21)9 T. Lauritsen et al., PRL89(22)28251 B(E1) ~ 1-4 -1-3 W.u. 1 to 1 times stronger than GDR tail Theory: J. Kvasil et al., PRC75 (27)3436 T. Nakatsukasa et al., PLB343(1995)19

Results are sensitive to the absolute value and to the spin-dependence of the calculated actions

Why to consider fluctuations? T. Dossing et al, Phys. Rep. 268 (1996) 1 Because they give a measure of the number of paths populating a certain region of the γ-γ spectrum. N path = 1/ Σ W i 2 Npath along the ridge gives a measure of the number of discrete bands as a function of angular momentum. But: each band is weighted by the square of its intensity. When all paths have the same probability, W i = 1/Ntrue and Npath= Ntrue

n mix 14 12 1 8 6 151 Tb I = 5 + -48 + -46 + 4 2 14 15 16 17 18 19 E [MeV] 194 Hg 194Hg: Bands survive in the mixing region. They are based on mixed states but are not damped (ergodic bands). Can they be populated and affect the fluctuation spectra?

Exp. versus Theory I Cranking Onset of Damping n branch < 2 N path 75 5 25 Number of Discrete Bands 151 Tb 196 Pb SD gated Ridge SD gated Ridge 75 5 25 I-2 + decay-out N path 1 75 5 25 2. Ridge 2. Ridge 125 1 75 5 25 + Population (N path = 1/w i2 ) N path 1 75 5 25 1. Ridge 1. Ridge 125 1 75 5 25 2 3 4 5 6 7 3 4 5 Strong suppression due to uneven population of the bands

194Hg : 1 Regular bands of chaotic nature A. Lopez-Martens et al., PRL1(28)1251

Conclusions Warm Rotation in Superdeformed nuclei is a test bench for cranked shell model at finite temperature & tunneling through potential barrier Experimental analysis 151 Tb & 196 Pb: Intensities and Fluctuations of Quasi-continuum spectra provide strong constraints on theoretical models Possibility of testing decay-out models over a wide spin range Data interpretation: Microscopic Monte Carlo simulation, reducing free parameters to the minimum Good overall agreement with data Indication of nuclear structure effects in the E1 strength at low excitation energy (~ 1 MeV) Collaboration

Participants to the Experiments Milano University & INFN: S. Leoni, A. Bracco, G. Benzoni, N. Blasi, S. Brambilla, F. Camera, F. Crespi, A. DeConto, P.Mason, D. Montanari, B. Million, M. Pignanelli, O. Wieland IRES, Strasbourg: G.Duchene, J.Robin, D. Curien, Th.Bysrki, F.A.Beck et al., Krakow, Poland: A.Maj, M. Kmiecik, P.Bednarczyk, W. Meczynski, J. Styczen, et al. NBI, Copenhagen: B. Herskind, G. Hagemann, G. Sletten et al. Oliver Lodge Laboratory, University of Liverpool: P.J.Twin KTH, Stockholm: A.Odahara, K.Lagergren + EUROBALL collaborations Theory: M.Matsuo (Niigata University), Y.R.Shimizu (Kyushu University), E.V.

Outlook: detailed analysis of other cases 1. The SD NUCLEUS 152 Dy So far investigated in details via a parameter dependent model T.Lautitsen et al., PRC75(27)6439 2. The PECULIAR case of 194 Hg 1 Regular bands of chaotic nature A. Lopez-Martens, PRL1(28)1251 chaotic intrinsic-motion regular rotational decay chaotic rotational decay Can theory reproduce this? ergodic bands rotational damping

Yb (Z=7) I=4h,U=2MeV highly aligned orbits E2 strength 16 176 Yb rot 98 168 Yb I no highly aligned orbits I-2 I-2 neutron rot for U 2 MeV N P

151 Tb: simulated discrete SD bands population 1 4 = + YRAST 1 3 even spins 1 4 = - 1 3 even spins Intensity 1 Intensity 1 1-1 1-1 1-2 25 3 35 4 45 5 55 6 65 7 1-2 25 3 35 4 45 5 55 6 65 7 75 1 4 = + 1 4 = - 1 3 odd spins 1 3 odd spins Intensity 1 Intensity 1 1-1 1-1 1-2 25 3 35 4 45 5 55 6 65 7 1-2 25 3 35 4 45 5 55 6 65 7 75

Discrete Spectroscopy Info s 151 Tb 196 Pb Intensity (%) 14 12 1 8 6 4 2 % 8,x1 6 6,x1 6 4,x1 6 2,x1 6, Yrast 8 1 12 14 16.8-1.6 MeV Intensity (%) 14 12 1 8 6 4 1.3 % 6 4 2 E Yrast 1 2 3 4 5 6 7 8.2-.8 MeV 2 2 2 4 6 8 2 14 Number of SD band 2 4 6 8 2 14 Number of SD band 151 Tb tentative assignment SD SD 196 Pb firm assignment 3698 462 ND ND J. Robin et al., submitted to PRC A.N. Wilson et al., Phys. Rev. Lett. 95 (25)

Compound Nucleus Chaos SD Warm Rotation strongly interacting bands I I-2 I+2 T= Regular bands Order Superdeformation at Finite Temperature Challenging topic experiment: focus on ~1% -decay theory: cranking at T coupled to tunneling to ND well L m ND ND SD tunneling SD to ND minima present status experiment: partial infos on 143 Eu, 152 Dy, 194 Hg theory: schematic OR many parameters

Aim of this work: Step forward in Experiment & Theory EXPERIMENT: Study of warm rotation in 151 Tb and 196 Pb Evaluation of several independent experimental oservables stringent test of -decay flow THEORY: Development of new Monte Carlo model Based on microscopic calculations Towards a parameter free model

196 Pb microscopic calculations Cranked shell model at T E2 fragmentation U (MeV) U (Mev) 2. 1.5 1..5. 2. 1.5 1..5. H = H def J x + V res = - yrast (,+) 3 4 5 6 red even green odd = + n branch Width [kev] 1 1 1 196 Pb n branch = 2 n branch = 2 rot I-2 onset of damping I = 3 +..5 1. 1.5 2. 2.5 U [MeV] 25 2 15 1 5 I x 1/1 rot I = 2 + 3 4 5 K. Yoshida, M. Matsuo, NPA612(1997)126 I-2 I

The Entry Distribution: is NOT a parameter, it can be measured or calculated n d/di [arb. un.] 1 27 Al + 14 Te @ E lab =152 MeV 8 6 4 2 CN fusion (Grazing code) 151 Tb (Cascade) n 2 4 6 8 1 Entry 151 Tb E * [MeV] 4 3 2 1 ND SD E2 E1 E2 E1 tunneling SD yrast counts [arb.un.] 6 4 2 <I> = 58 ± 8 <U> = 8.3 ± 2.5 MeV 2 I=3 I=4 15 1 5 ND yrast 5 5 2 U [MeV] 5 5 2 U [MeV] 3 4 5 6 7 8 A. Winther, Nucl. Phys. A594, 23(1995) F. Pulhofer, Nucl. Phys. A28, 267 (1977)

Exp. versus Theory Analysis of Intensities 151 Tb 196 Pb exclusive quantities SD gated decay-flow inclusive quantities average decay-flow Intensity [%] Intensity [%] Intensity [%] Intensity [%] 1 3 E2 SD gated 1 1 3 SD gated Ridge 1 1 3 2. Ridge 1 1 3 1 1. Ridge Yrast I yrast = 2 % 3 4 5 6 7 Intensity [%] Intensity [%] Intensity [%] Intensity [%] 1 3 1 1 3 1 1 3 1 1 3 1 E2+M1 E2 SD gated SD gated Ridge Yrast I yrast = 1.3% 2. 2Ridge 1. Ridge 3 4 5 SD gated spectrum E2 Bump SD gated 1 st Ridge Rot-Plane 2 nd Ridge Total 1 st Ridge

Exp. versus Theory Analysis of Intensities 151 Tb 196 Pb exclusive quantities SD gated decay-flow inclusive quantities average decay-flow Intensity [%] Intensity [%] Intensity [%] Intensity [%] 1 3 2. Ridge 1 1 3 1 1 3 E2 SD gated 1 1 3 SD gated Ridge 1 1. Ridge Yrast I yrast = 2 % 3 4 5 6 7 Intensity [%] Intensity [%] Intensity [%] Intensity [%] 1 3 1 1 3 1 3 1 1 3 1 E2+M1 Yrast I yrast = 1.3% E2 SD gated SD gated Ridge 2. 2Rid Ridge 1. Ridge 3 4 5 SD gated spectrum E2 Bump SD gated 1 st Ridge Rot-Plane 2 nd Ridge Total 1 st Ridge

SD gated quantities are sensitive to E1/E2 balance Enhanced octupole vibrations have been observed in A= 15 and A = 19 196 Pb, 19 Hg 194 Hg T(E1) [s -1 ] 4x9 3x9 2x9 1x9 SD 5 5 2 25 3 E [MeV] 7 6 enhanced GDR T. Lauritsen et al., PRL89(22)28251 D. Rosbach et al. PLB513(21)9 B(E1) ~ 1-4 -1-3 W.u. 1 to 1 times stronger Theory: Nilson + QRPA J. Kvasil et al., PRC75 (27)3436 T. Nakatsukasa et al., PLB343(1995)19 T(E1) [s -1 ] 5 4 3 2 Standard GDR 151 Tb.5 1. 1.5 2. 2.5 E [MeV]

The E1 decay strength: the tail of the GDR f GDR = Sum of 3 Lorentzian ND SD T 3 ( E1) H1 n K E1 B( E1) E E2 E1 E2 E1 T(E1) [s -1 T(E1) [s -1 ] 4x9 19 3x9 19 2x9 19 1x9 19 tunneling ND SD Hindrance factor ~ 1-2 (tuned to reproduce the intensity of the yrast band) K.E.G. Lobner, Phys. Lett. 26B, 369(1968) G. Leander, PRC38, 728(1988) 55 5 2 25 3 E [MeV]

n mix 14 12 1 8 6 151 Tb I = 5 + -48 + -46 + 4 2 14 15 16 17 18 19 E [MeV] 4 151 Tb 196 Pb 194 Hg 194 Hg % 2 <2 2-4 n mix 4-6 >6

SD gated quantities are sensitive to E1/E2 balance Enhanced octupole vibrations have been observed in A= 15 and A = 19 196 Pb, 19 Hg 194 Hg T(E1) [s -1 ] 4x9 3x9 2x9 1x9 SD 5 5 2 25 3 E [MeV] 7 6 enhanced GDR T. Lauritsen et al., PRL89(22)28251 D. Rosbach et al. PLB513(21)9 B(E1) ~ 1-4 -1-3 W.u. 1 to 1 times stronger Theory: Nilson + QRPA J. Kvasil et al., PRC75 (27)3436 T. Nakatsukasa et al., PLB343(1995)19 T(E1) [s -1 ] 5 4 3 2 Standard GDR 151 Tb.5 1. 1.5 2. 2.5 E [MeV]