NZ equilibration in two and three bodies dynamically deformed nuclear systems ( 7 Zn+ 7 Zn @ 35 MeV/nucleon) Dr. 1 1Cyclotron Institute, Texas A&M University, TAMU DNP 218, Kona 1
Systems studied: 7 Zn+ 7 Zn at 35 MeV per nucleon using NIMROD DNP 218, Kona!2
n p DNP 218, Kona!3
Velocity gradient and surface tension amplifies instabilities DNP 218, Kona!4
DNP 218, Kona!5
TLF* E HF LF 3F α - Simultaneous or double rupture? - 3F separated first? - Fragments aligned? String of pearls? - Fragments aligned with the Vcm vector? DNP 218, Kona!6
Previously on NZ equilibration Time-dependence of NZ equilibration examining composition of HF-LF emitted from PLF* vs the breakup alignment angle (serves as clock for equilibration) PRL 118, 6251 (216) PRC 95, 44 (217) Variation of composition vs alignment angle shows exponential behavior for LF & HF, suggesting 1st-order kinetics, for all systems studied Small systematic effects in composition for reactions of relatively npoor projectile with nrich target No significant differences in rate cts between systems of different initial composition DNP 218, Kona = (N-Z)/A Yield and measured composition, used to extract an estimate for purely dynamical component.12 7 Zn +.12 7.1.8.8.6.6.4.4.2.2.1 2 4 6 8 1 12 Zn + Zn ZH=12, ZL =7 14 16.12 18 6 8 1 12 14 16 18 6 8 1 12 14 16 18 Ni.8.6.4.4.2.2 2 4.1.6 2 Ni Zn + HF LF.8 Ni +.1.12 Zn 4 6 8 1 12 14 16 18 2 4 α (degrees)!7
Symmetrize Dalitz plot.7.6.5 Preliminary HF LF rd 3 F ZH ZL Z3 ZH ZL Yield.4.3 ZL Z3.2 ZH Z3.1 ZH 5 1 15 2 25 Atomic Number ZL ZH ZT ZL ZT Z3 ZT Z3 Preliminary 1 Simmetrized Dalitz-Z.8.6.4 5 4.2 3 2Z s are approximately equal 1Z is different (and bigger).2.4.6.8 2 1 DNP 218, Kona 1 1.8.6.4.2.2.4.6.8 1!8
Velocity distributions of HF, LF and 3F Dashed lines = beam velocity (i.e..27c) and mid-velocity (i.e..13c). Preliminary The three of them peaked above midvelocity (indicating that HF, LF and 3F originate from the PLF*). LF (3F) produced at velocities > than mid-velocity and < than the HF (LF). HF produced closer to beam velocity. Normalized Yield.6.5.4.3 Z H =12, Z =7, Z =3 L 3 HF LF rd 3 F.2 Hierarchy in the velocity distributions (correlated to the charge sorting): HF forward from LF, and both are forward respect to 3F..1.5.1.15.2.25.3.35.4.45.5 Velocity (c) DNP 218, Kona!9
Rupture scenarios: Who is separating first? Example of the rupture scenario where Z3 separates first First rupture Second rupture Vcm Z3 snap Vrel α ZL ZH snap ZH ZL ZH ZL Vrel Vcm β DNP 218, Kona!1
.1.8 Preliminary Rupture scenarios: Angular distributions Rupture I ZHZL-Z3 ZH-ZL Normalized Yield.6.4 α.2 β Normalized Yield Normalized Yield Normalized Yield 2 4 6 8 1 12 14 16 18 angle (degrees).1 Rupture II ZLZ3-ZH.8 ZL-Z3.6.4.2 2 4 6 8 1 12 14 16 18 angle (degrees).1.8.6.4.2 Rupture III ZHZ3-ZL ZH-Z3 Statistical contribution (color hatches) ZH forward from ZL, and both forward from Z3 Peak in the angular yield when Z3 and ZL separate first When ZH separates first, the dynamical contribution seems to be smaller that when Z3 or ZL separate first DNP 218, Kona 2 4 6 8 1 12 14 16 18 angle (degrees) angle (degrees)!11
Rupture scenarios: correlations between the angles of both breakings Another way to show the angular correlation More dyn-dyn and dyn-stad combinations ZHZL 18 16 14 12 1 8 6 4 7 6 5 4 3 2 ZLZ3 18 16 14 12 1 8 6 4 7 6 5 4 3 2 More stad-dyn combinations 2 1 2 1 2 4 6 8 1 12 14 16 18 (ZHZL_Z3) Preliminary 2 4 6 8 1 12 14 16 18 (ZLZ3_ZH) More dyn-dyn and stad-dyn combinations ZHZ3 18 16 14 12 1 8 6 4 8 7 6 5 4 3 2 18 16 14 12 1 8 6 4 dyn stat dyn dyn stat stat dyn stat 2 1 2 2 4 6 8 1 12 14 16 18 (ZHZ3_ZL) 2 4 6 8 1 12 14 16 18 Looks like HF is not separating first DNP 218, Kona!12
Angles between the relatives velocities & Angle between the projections of HF, LF and 3F velocities, in a plane perpendicular to Vcm Cartoon of a plane perpendicular to Vcm and a plane formed by the three fragments Plane perpendicular to Vcm asymmetric angular projection plane perpendicular to Vcm ZH ZH Z3 + symmetric angles on the relatives velocity plane of the 3 fragments Z3 Vcm? Vcm ZL ZL DNP 218, Kona Plane of the relatives velocities!13
Angles between the projections of HF, LF and 3F velocities, in a plane perpendicular to Vcm The majority of events populate combinations of (18 o, o,18 o ) 1.8 Preliminary Dalitz-VplanePerpToVcm LF3F 6 Plane perpendicular to Vcm.6.4 5.2.2.4 4 3 2 Vcm Vcm.6.8 HFLF HF3F 1 1 1.8.6.4.2.2.4.6.8 1 Two fragments are closer to each other, and apart from the other one. DNP 218, Kona!14
Angles between the relatives velocities Plane formed by the 3 fragments Plane of the relatives velocities Preliminary LF3F-HF3F VrelHF3F VrelHFLF Vcm VrelLF3F HFLF-HF3F HFLF-LF3F It looks like red, blue and green angles are similar (slightly following r > b > g) DNP 218, Kona!15
the show down 1.8 Preliminary Dalitz-VplanePerpToVcm LF3F 6 Preliminary LF3F-HF3F.6.4 5.2 4.2 3.4 2.6.8 HFLF HF3F 1 HFLF-HF3F HFLF-LF3F 1 1.8.6.4.2.2.4.6.8 1 Two fragments seemingly closer to each other and apart from the third one The three angles are approximately similar Plane perpendicular to Vcm Combining both:. Relative velocities plane is oblique to Vcm. The three fragments don t seem to be aligned Vcm Vcm DNP 218, Kona Plane of the relatives velocities!16
Angle between Vrels and Vcm: oblique angle Preliminary hangle_oblique_vrels_vcm Counts 18 16 14 12 1 8 6 The angular momentum is controlling the rupture, because it is occurring in preferential planes. Plane perpendicular to Vcm 4 2 ZH ZH Z3 2 4 6 8 1 12 14 16 18 Angle(degrees) Z3 Vcm The rotation is perpendicular to the direction of the movement. The system breaks due to rotation and in a double rupture, where the plane of the Vrels is the plane perpendicular to the angular momentum, (very oblique to Vcm) Vcm ZL Plane of the relatives velocities ZL DNP 218, Kona Plane of the relatives velocities!17
Comparison to Boltzmann-Langevin One Body (BLOB): Something else we want to explore BLOB calculation: Violent nuclear reaction dynamics, neck fragmentation Dynamical transport model Unifies in a common approach the description of fluctuations in nuclear matter, and a predictive description of the disintegration of nuclei into nuclear fragments Characterize the very fast early stages of the collision process which are out of equilibrium 7Zn+ 7 Zn 35AMeV (b=5.7) producing 45 Ca, 4 Cl, 6 Li,d (at 4fm/c) https://www.youtube.com/watch?v=l5ul-iux1ik Paolo Napolitani Sizable neck fragments work in progress DNP 218, Kona!18
Summary Studied time dependence of the NZ equilibration, examining composition of fragments from the PLF* vs angle (serves as clock for the equilibration). Composition follows exponential behavior suggesting first-order kinetics. We extracted a purely dynamical composition component, found target effects and no significant differences in rate constants among HF-LF or different systems. Hierarchy in the velocity distributions (strongly correlated to the charge sorting) Sizes: one big, two smaller and similar. The relative velocities plane is oblique to Vcm and the three fragments don t seem to be aligned. The physics is dominated by the angular momentum, so the system breaks due to rotation and in a double rupture scenario. Comparing with dynamical transport models (BLOB), with sizable neck fragments DNP 218, Kona!19
Thank you! SJY group and collaborators: A. Rodriguez Manso, A.B. McIntosh, J. Gauthier, K. Hagel, L. Heilborn, A. Jedele, Z.W. Kohley, Y.W. Lui, L.W. May, A. Zarrella, R. Wada, A. Wakhle and S.J. Yennello Department of Energy DE-FG2-93ER4773 Welch Foundation A-1266 DNP 218, Kona!2
Backup Slides DNP 218, Kona!21
Neutron Ion Multidetector for Reaction Oriented Dynamics (NIMROD) Multidetector array for reactions between massive target and projectiles: Total of 228 detector modules arranged in 14 annular rings (2-3 detectors/module) Projectile energies go from 2MeV to 4GeV 4π coverage (3.6 to 167. ) = nearly complete geometrical coverage excellent isotopic ID DNP 218, Kona Detector Module ΔE E!22
Angular (α) distributions: 7Zn.9.8 Lifetime of PLF* correlated with rotation angle: + 7Zn ZH=14, ZL=5 Dynamical Statistical.7.6.5.4.3.2.1.9 Excess yield largest and most strongly aligned for most asymmetric splits..8 4 6 8 1 12 14 16 18.7 Normalized Yield Angular distribution peaked for most aligned configuration. Decreases in yield with decreasing alignment. 2 ZH=14, ZL=7 Less aligned decays represent longer decay times..6.5.4.3.2.1.9.8 2 4 6 8 1 12 14 16 18 ZH=12, ZL=7.7 YT = Yd + Ys.6.5.4 Dynamical yield dominates at small angles and decreases as α increases..3.2.1 DNP 218, Kona 2 4 6 8 1 12 14 16 18 α (degrees)!23
Composition vs decay alignment: Overall composition Dynamical composition.14.12 Z H =14, Z =5 L.14.12 Z H =14, Z =5 L ΔObservable = fsδs + fdδd.1.1.8.8 As angle of rotation increases, Δ HF (Δ LF ) start off very n-poor (n-rich), then evolve towards each other. Exponential fit: a + b e -c(α) a = equilibrium value c = rate constant for equilibration First-order kinetics. (N-Z)/A =.6.4.2.14.12.1.8.6.4 Z H =14, Z =7 L HF LF = (N-Z)/A dyn.6.4.2.14.12.1.8.6.4 Z H =14, Z =7 L HF LF Dynamical composition generally follows the same trend as the overall composition (more extreme)..2.14.12.1 Z H =12, Z =7 L 7Zn + 7 Zn.2 2 4 6 8 1 12 14 16 18.14.12.1 Z H =12, Z =7 L 7Zn + 7 Zn Rate of change of the composition unaffected by correction.8.6.8.6.4.4 DNP 218, Kona.2 2 4 6 8 1 12 14 16 18 α (degrees).2 2 4 6 8 1 12 14 16 18 α (degrees)!24
Composition vs decay alignment: for different systems Δ/α ( 7 Zn, Ni) correlations essentially the same..12.1 7 Zn + 7 Zn Ni + Ni Δ/α ( Zn) correlation shifted to lower values (i.e. lower equilibrium composition). Same rate constants and change from initial to final value. Comparing the symmetric and the asymmetric systems: (N-Z)/A =.8.6.4.2.12.1.8.6 Zn + Z H =12, Z HF LF Zn =7 L Zn + Ni Δ/α correlations for the symmetric system are systematically lower target effects..4.2 2 4 6 8 1 12 14 16 18 2 4 6 8 1 12 14 16 18 α (degrees) DNP 218, Kona!25
Equilibrium composition <Δ> = a + b e -c(α) ZL, ZH values fairly clustered equilibrium value for LF (HF) depends on ZL (ZH) but not on ZH (ZL). Comparing 7 Zn ( ) and Ni ( ) systems (similar n-rich system composition): ~ same equilibrium composition. Zn ( ): consistently less n-rich equilibrium compositions Equilibrium Composition.18.16.14.12.1.8.6.4.2 nd 2 Heaviest Frag.(LF) 7 Zn + Zn Zn + Zn Ni + Ni Zn + Ni 7 Heaviest Frag.(HF) 7 7 Zn + Zn Zn + Zn Ni + Ni Zn + Ni The asymmetric Zn + Ni system ( ): slightly + n-rich daughters than the n-poor syst. slightly - n-rich daughters than the n-rich syst..2 2 4 6 8 1 12 14 16 18 Atomic Number target effect DNP 218, Kona!26
Rate constant for equilibration: <Δ> = a + b e -c(α) Exponential slope Average rate constant zs -1 : LF 4 ± 1 HF 4 ± 2 Relevant parameter to calculate the equilibration times. Describes how fast the equilibration occurs within the PLF*. ) -1 Rate Constant (zs 18 16 14 12 1 8 6 4 2 2 18 16 14 12 1 8 6 4 2 7 Zn + HF LF Zn + 7 Zn Zn 18 16 14 12 1 8 6 4 2 Ni + Zn + 2 2 2 4 6 8 1 12 14 16 1 2 4 6 8 1 12 14 16 18 Ni Ni Atomic Number Agreement of rates force driving the equilibration is independent of the size of both partners only depends on the difference in asymmetry DNP 218, Kona!27
Angular (α, φin and Θout) distributions: why α? Consistent with a significant amount of dynamical decay. Statistical decay α and φin: Symmetric distribution about 9 ~ sinusoidal for low spin Squashed out symmetrically toward edges as spin increases. Distributions strongly peaked toward and asymmetric large yield of non-statistical decay. The distributions: φin (Θout near 7 or 11 ) looks like the total φin φin (Θout near 3 or 15 ) looks much flatter. dynamical yield is preferentially closer to the plane. Normalized Yield Normalized Yield.9 Z =7 H =12, Z L.8.7.6.5.4.3.2.1.9.8.7.6.5.4.3.2.1 5 1 15 α (degrees) 5 1 15 θ out (degrees) Normalized Yield (degrees) θ out.9.8.7.6.5.4.3.2.1 15 1 5 5 (degrees) 1 15 ϕ in ϕ in 5 1 15 (degrees) 45 4 35 3 25 2 15 1 5 φin angle should describe the deformation alignment in a similar manner as it does the α angle. DNP 218, Kona!28
Composition vs decay alignment, for α, φin : Rotation of the PLF* as it decays into HF and LF is predominantly around an axis perpendicular to the RP. (N-Z)/A.1.9.8.7 Z H =12, Z =7 L HF (vs α) HF (vs ϕ LF (vs α) LF (vs ϕ ) in ) in =.6 Composition as a function of these two angles expected to be similar..5.4 2 4 6 8 1 12 14 16 18 angle (degrees) Rotation axis can be canted from perpendicular, we expect that α would provide a truer measure of the time. DNP 218, Kona!29
Secondary decay (SD): using GEMINI++ open markers: <Δ> of the final state fragments Varying the initial excitation energy: In both cases ( ), exponential dependence maintained. Shift to lower composition and muting of the amplitude (stronger for higher excitation energy).14.12.1.8 GEMINI++ Decay to Z=7 Primary E* = 1AMeV Secondary E* = 1AMeV Primary E* = 2AMeV Secondary E* = 2AMeV Varying the starting <Δ>: System with initially larger asymmetry shifted down strongly by SD. (System farther from valley of stability feels a stronger force driving it back toward the valley) After SD, the more n-rich system remains more n-rich. (N-Z)/A =.6.14 =.1.12.1 =.1 =.8 =.8.8 The trend is not destroyed or created, and the characteristic rate of the exponential is retained..6 2 4 6 8 1 12 14 16 18 α (degrees) DNP 218, Kona!3
Rate of equilibration: Average rate constant degree -1 : Symmetric systems: LF.3 ±.1 HF.2 ±.1 Asymmetric system: LF.2 ±.1 HF.1 ±.2 Relevant parameter to calculate the equilibration times. Describes how fast the equilibration occurs within the PLF*. ) -1 Rate Constant (degrees.16.14.12.1.8.6.4.2.2.16.14.12.1.8.6.4.2.2 7 Zn + HF LF Zn + 7 Zn Zn.16.14.12.1.8.6.4.2 Ni + Zn +.2 2 4 6 8 1 12 14 16 1 2 4 6 8 1 12 14 16 18 Ni Ni Atomic Number Agreement of rates force driving the equilibration is independent of the size of both partners only depends on the difference in asymmetry DNP 218, Kona!31
Time-Scale t = α/ω ω(angular frequency), J (angular momentum) where ω = Jħ/I eff Fit The J is determined using the width of the out-of-plane α particle distribution GEMINI++ simulations: reproducing this width can be done with spin from 1ħ (E*/ A=.8MeV) to 5ħ (E*/A=1.2MeV). We can take J=22ħ with a factor of 2.2 uncertainty. Moment of inertia is calculated using a 2 touching spheres model: Ieff = m ZH r CM,ZH 2 + ⅖m ZH r ZH 2 + m ZL r CM,ZL 2 + ⅖m ZL r ZL 2 Using a complete rotational period: t = (1-4) zs DNP 218, Kona!32