On ring and bubble formations in heavy ion collisions

On ring and bubble formations in heavy ion collisions 1 2 3 corresp first-author Corresponding author. E-mail address: Konstantin.Cherevko@gmail.com Corresponding author. E-mail address: Konstantin.Cherevko@gmail.com Konstantin.Cherevko@gmail.com 1 2 3 3 1 2 4 corresp Corresponding author. E-mail address: fszhang@bnu.edu.cn Corresponding author. E-mail address: fszhang@bnu.edu.cn fszhang@bnu.edu.cn College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China Beijing Radiation Center, Beijing 100875, China Beijing Radiation Center, Beijing 100875, China Physics Department, Taras Shevchenko National University of Kyiv, 4 Glushkova av., Kyiv 03022, Ukraine Physics Department, Taras Shevchenko National University of Kyiv, 4 Glushkova av., Kyiv 03022, Ukraine Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator of Lanzhou, Lanzhou 730000, China Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator of Lanzhou, Lanzhou 730000, China The work is devoted to the implementation of the hydrodynamic laws to the head-on heavy ion collisions within the energy range 50–100 MeV/ A . The hydrodynamic mechanisms of the bubble and ring structures formation are investigated. It is shown that there is a possible hydrodynamic explanation of the different structures being formed in the case of soft ( K =200 MeV) and stiff ( K =400 MeV) equations of state. Within the suggested approach the final geometry of the system is defined in the initial stage of the collision and is very dependent on the sound velocity in the nuclear matter. The obtained results are in a good correspondence with the Boltzmann-like transport theory calculations and the experimental data for the selected energy range. Nuclear incompressibility Head-on collisions Hydrodynamics 1 Introduction 1 The knowledge of the nuclear equation of state (EOS) is one of the fundamental goals in nuclear physics [ 1 1 , 2 2 , 3 3 , 4 4 ] which has not yet been achieved [ 5 5 , 6 6 ] . The possibility of extracting information on the EOS at high baryon density, is restricted to two fields of research that are observations of astrophysical compact objects and studies of hot nuclear systems, created in high energy proton induced reactions or in heavy-ion collisions (HIC). Probably this explains the intensive studies of heavy ion collisions during the last several decades. The very exciting parts of them are the studies devoted to the nuclear multifragmentation (MF) when a number of intermediate mass fragments 3 < Z < 20 (IMF) are being formed from the compressed and excited system [ 7 7 , 8 8 , 9 9 ] . IMFs may provide a unique probe to study the reaction mechanism and hot nuclear matter properties [ 10 10 ] . In spite of such an interest to the topic, there are still many questions to be answered. The serious challenge in the studies is the inclusiveness of the heavy ion collisions experiments [ 11 11 ] . The nowadays detectors are very sophisticated and allow the high detection rate but still some data are missing [ 12 12 ] . At this step a powerful instrument of computer simulations is of high importance. There is a number of transport theories used for simulating collisions [ 13 13 , 14 14 , 15 15 ] and statistical decay models [ 16 16 ] that provide interesting microscopic results and give predictions on possible MF mechanisms. Within the Boltzmann-Uehling-Uhlenbeck (BUU) model it was shown that at the energy range 60–75 MeV/ A "bubble" or "doughnut" structures can be formed depending on the stiffness of the nuclear EOS [ 17 17 , 18 18 ] . Later the formation of "doughnut" structures was confirmed experimentally [ 19 19 ] . Unfortunately up to now there is no model describing in details the underlying physical mechanisms of the exotic shapes formation in the head-on HIC. There are different possible ways to look for the explanation. One is to continue with the computer simulations. In that case one should bear in mind that we are still far from having models that are formally well founded, practically applicable, and sufficiently realistic to be quantitatively useful [ 20 20 ] . Another way is to consider the qualitative Van der Waals and nucleon-nucleon interactions similarity [ 21 21 ] and to look for the explanation in thermodynamics [ 22 22 , 23 23 ] or hydrodynamics well developed for ordinary liquids. In that case a number of possible physical mechanisms being responsible for the phenomena can be suggested. Among them are the mechanism that includes the shockwave formation and interference and the one that studies the thermodynamic properties of the system (temperature, surface tension) together with the influence of the angular momentum on the system behavior [ 24 24 ] . The other possibilities are the possible capillary effects with "capons" [ 25 25 ] being formed on the nuclei surface and the spinodal decomposition of the system. Therefore, some deep analysis of the phenomena is needed to choose the appropriate decay channel. The first step toward a better understanding of the process is the analysis of the qualitative picture that comes from the BUU calculations and the available experiments [ 17 17 , 18 18 , 19 19 ] . First to mention is the fact that the expanding velocity of the outer surface of the system is much smaller than the expanding velocity of the inner surface that suggests the importance of the surface effects. Second is the system behavior during the final stage. BUU calculations predict the simultaneous breakup with few fragments of similar masses and low kinetic energies. This suggests that the dynamic process might be responsible for the system evolution but at the breakup stage Coulomb force seems to be of great importance. Third is the angular distribution of fragments. In the case of soft EOS it is almost isotropic when for the stiff EOS fragments are in the plane perpendicular to the beam direction. From all the above the hydrodynamic approach seems to be quite attractive as it allows explanation of the observed behavior of the system. In this work, we report an attempt to develop the hydrodynamic approach with the shockwave mechanism [ 17 17 , 26 26 ] that can explain the observed phenomena and reveal its physical nature. 2 Model description 1 When studying the system from the hydrodynamic point of view the question arises whether such a macroscopic theory works for the nuclear systems. Originally the idea to apply such a macroscopic theory to nuclei was introduced in 1970s–1980s by Siemens [ 27 27 ] and Stocker [ 28 28 ] . The domain of their studies was in a higher energy range than we are interested in, nevertheless, the existing confirmation for the applicability of hydrodynamics for nuclear systems suggests that it should be possible to use it for the energy range in focus [ 29 29 ] . The symmetry of the system allows the simplification of the model. Namely from the collision of two nuclei it is possible to change to the collision of a spherical nucleus of radius R with a rigid wall that moves toward it with the velocity υ 0 . The difference is in the viscous effects, therefore, the slip boundary condition on the wall surface should be added. In our studies we start from the acoustic approximation developed in the theory of liquid droplet collisions [ 30 30 , 31 31 ] . Within this model according to Huygens principle the expanding nucleus edges emit wavelet that propagate with the sound velocity in all the directions and are the only source of the surface distortion. In that case the compression stage when the nucleus remains spherical lasts until the contact point velocity <math id="M1"><mrow><msub><mi>υ</mi><mi>c</mi></msub><mo>=</mo><msub><mi>υ</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>R</mi><mo>−</mo><msub><mi>υ</mi><mn>0</mn></msub><mi>t</mi><mo stretchy="false">)</mo><mo>/</mo><msqrt><mrow><mn>2</mn><mi>R</mi><msub><mi>υ</mi><mn>0</mn></msub><mi>t</mi><mo>−</mo><msubsup><mi>υ</mi><mn>0</mn><mn>2</mn></msubsup><msup><mi>t</mi><mn>2</mn></msup></mrow></msqrt></mrow></math> [ 32 32 ] (can be found from the geometry of the system) is higher than the sound velocity υ s inside the nucleus ( Fig.1 Fig.1 ). That first compression stage is characterized by the uniform pressure along z -axis and the pressure near the contact edge being higher than the pressure in the central part [ 31 31 ] . After that lateral jetting occurs and the pressure gradually decreases while the form of the nuclei is changing ( Fig.1 Fig.1 (a)). In all the calculations a simple approximate relation between the sound velocity υ s and the particle velocity υ p is used: <math display="block" id="M2"><mrow><msub><mi>υ</mi><mi>s</mi></msub><mo>=</mo><msub><mi>υ</mi><mrow><mi>s</mi><mn>0</mn></mrow></msub><mo>+</mo><mi>k</mi><msub><mi>υ</mi><mi>p</mi></msub><mo>,</mo></mrow></math> where k is a coefficient that is used as an adjusting parameter in the model and υ s 0 is the ambient sound velocity. The compression part of the nuclear matter EOS [ 33 33 ] <math display="block" id="M3"><mtable columnalign="left"><mtr><mtd><mi>U</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>A</mi><mrow><mo>(</mo><mrow><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mi>B</mi><msup><mrow><mo>(</mo><mrow><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mi>γ</mi></msup></mtd></mtr><mtr><mtd><mi>P</mi><mo>=</mo><msup><mi>ρ</mi><mn>2</mn></msup><mfrac><mrow><mo>∂</mo><msub><mi>E</mi><mi>B</mi></msub></mrow><mrow><mo>∂</mo><mi>ρ</mi></mrow></mfrac><mo>,</mo><mi/><mi/><mi/><msub><mi>E</mi><mi>B</mi></msub><mo>=</mo><mfrac><mn>3</mn><mn>5</mn></mfrac><msub><mi>E</mi><mi>F</mi></msub><mo>+</mo><mfrac><mn>1</mn><mi>ρ</mi></mfrac><mstyle displaystyle="true"><mrow><mo>∫</mo><mrow><mi>U</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow></mrow></mstyle><mi>d</mi><mi>ρ</mi></mtd></mtr></mtable></math> together with the continuity equations and momentum and energy conservation <math display="block" id="M4"><mtable columnalign="left"><mtr><mtd><mfrac><mrow><mo>∂</mo><mi>ρ</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mi>d</mi><mi>i</mi><mi>v</mi><mo stretchy="false">(</mo><mi>ρ</mi><mover accent="true"><mi>υ</mi><mo stretchy="true">→</mo></mover><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi>ρ</mi><msub><mi>υ</mi><mi>i</mi></msub><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi>ρ</mi><msub><mi>υ</mi><mi>i</mi></msub><msub><mi>υ</mi><mi>k</mi></msub><mo>+</mo><mi>p</mi><msub><mi>δ</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><msub><mi>x</mi><mi>k</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><mo>∂</mo><mo stretchy="false">(</mo><mi>ρ</mi><mi>E</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mi>d</mi><mi>i</mi><mi>v</mi><mo stretchy="false">(</mo><mi>υ</mi><mo stretchy="false">(</mo><mi>E</mi><mo>+</mo><mi>p</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo><mo>=</mo><mn>0</mn></mtd></mtr></mtable></math> are used for phenomena analysis. Parameters A , B and γ are chosen to give the different values of incompressibility [ 34 34 , 35 35 ] . U ( ρ ), E B , ρ , P and υ are the Skyrme type mean field potential, binding energy per nucleon, density, pressure and particle velocity, respectively. In addition simple geometrical considerations are used to find the relations between the normal and tangential components of sound and particle velocities in the disturbed and undisturbed parts of the system [ 32 32 ] . The introduced acoustic model is applicable during the first stage of the collision until there is no lateral jetting [ 31 31 , 32 32 ] . This fact allows calculating the jetting time from the anomalous behavior of the contact edge velocity ( Fig.2 Fig.2 ). One may see that <math id="M5"><mrow><msub><mi>U</mi><mi>c</mi></msub></mrow></math> decreases from infinity to some final value and after that starts increasing. Such a picture seems to be unphysical and the minimum should indicate the boundaries of the model applicability and therefore it corresponds to the jetting time. When jetting occurs the shockwave starts to travel up along the free surface of the nucleus. This allows using the model of a jet formation between the rigid wall and the shockwave due to the pressure difference and the existing particles flow inside the nucleus. In this case assuming that the matter in the jet is at normal density ρ 0 it is possible to start with the hydrodynamic equations for the jet parameters [ 36 36 ] that take into account work against the surface forces L with the surface tension coefficient σ : <math display="block" id="M6"><mrow><mi>d</mi><msub><mi>L</mi><mi>s</mi></msub><mo>=</mo><mn>8</mn><mi>π</mi><mi>σ</mi><mi>R</mi><mi>d</mi><mi>R</mi><mo>,</mo></mrow></math> <math display="block" id="M7"><mrow><mi>L</mi><mo>=</mo><mfrac><mi>π</mi><mn>2</mn></mfrac><msubsup><mi>R</mi><mrow><mtext>jet</mtext></mrow><mn>2</mn></msubsup><msub><mi>ρ</mi><mrow><mtext>bound</mtext></mrow></msub><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mi>t</mi></munderover><mrow><msubsup><mi>υ</mi><mrow><mtext>jet</mtext></mrow><mn>3</mn></msubsup><mi>d</mi><mi>t</mi></mrow></mrow></mstyle><mo>.</mo></mrow></math> From Eqs.(1)–(5) it is possible to calculate the main characteristics of the system geometry depending on the density at contact point ρ bound : <math display="block" id="M8"><mrow><msub><mi>R</mi><mrow><mtext>jet</mtext></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mrow><mn>32</mn></mrow></mfrac><mfrac><mrow><msubsup><mi>ρ</mi><mrow><mtext>bound</mtext></mrow><mn>2</mn></msubsup><msubsup><mi>υ</mi><mrow><mtext>jet</mtext></mrow><mn>4</mn></msubsup><msup><mi>t</mi><mn>2</mn></msup></mrow><mrow><msub><mi>ρ</mi><mn>0</mn></msub><mi>σ</mi><mn>,</mn></mrow></mfrac></mrow></math> where R jet and υ jet are the size of the jet (Figs.3 and 4) and the particle velocity in the origin of the jet. To check the hydrodynamic approach we have chosen 93 Nb+ 93 Nb system with the impact velocity V =0.34 c that is the same as in the existing BUU calculations [ 17 17 ] . The parameter k =0.77 is adjusted from the comparison of jetting time in nuclear fluid dynamics calculations for 80 Kr+ 80 Kr at 400 MeV/ A [ 26 26 ] and in our model. 3 Results and discussion 1 All the calculations are provided for the two types of the EOS. The results are presented in Tables 1 and 2. 1 1 1 1 This work 1 1 BUU [ 17 17 ] 1 1 Jetting time, τ jet , fm/ c 1 1 3.5 1 1 1 1 Shape at τ jet , Δ D/D 1 1 0.2 1 1   1 1 Pressure release time, fm/ c 1 1 63 1 1 40~80 1 1 Maximum density, ρ / ρ 0 1 1 ~1.76 1 1 ~1.5 1 1 1 1 Soft EOS 1 1 Stiff EOS 1 1 Shape change, Δ D/D 1 1 0.5 1 1 0.4 1 1 Geometry, ( D +2 L )/ D 1 1 1.1 1 1 1.7 1 1 Time from collision, fm/ c 1 1 24 1 1 21 First thing to mention is the good correspondence of the pressure release time (time when the expansion shockwaves reach the center of the nuclei and therefore the density of the system drops down) obtained in our calculations with that in the BUU ( Table 1 Table 1 ). Our calculations have shown that from the hydrodynamic point of view cases of "soft" and "stiff" EOS differ not only quantitatively but different mechanisms are involved in the system evolution. The geometry of the system for the two cases at a moment of time when the shockwaves reach the top of the nuclei and its evolution are shown in Figs.3 and 4. From the comparison one may see difference in the qualitative pictures for different EOS. Soft EOS ( Fig.3 Fig.3 ): Such geometry allows the shockwaves produced on impact to be reflected from the free boundary and focused on the symmetry axis of the system. In that case rarefaction with cavity formation along the axis occurs ( Fig.3 Fig.3 (b)). Next steps in such a scenario are the cavity collapse and bubble entrapment inside the system. This mechanism shows the possibility of the bubble formation in the collision. Stiff EOS ( Fig.4 Fig.4 ): In that case system looks not like a spherical nuclei with small side parts (as it is for soft EOS) but rather like an expanding "pancake" . In our opinion this picture has some similarity with the problem of expanding liquid sheet [ 37 37 ] . Therefore we suggest that in this case the rarefaction along the symmetry axis does not play that important role but one rather observes the "pancake" becoming thinner and the thicker rim being formed due to the surface and viscous forces. Both mechanisms give low values of the expanding velocities at the final state corresponding to the low kinetic energies of the fragments obtained in the BUU calculations. We would also like to mention the possibility of the other geometry such as uniform system for different values of the impact energy or intermediate values of the incompressibility coefficient. Its occurrence as well defining the energy range where the exotic topologies can be observed require some further studies of the possible system splashing or recoil. Intuitively, the difference between the observed qualitative pictures has the origin in the initial stage when the compression pressure and temperature are different for different incompressibility coefficients. At this stage the only governing parameter is the sound velocity in the system. Viscous and surface forces come to play later on when the geometry has been defined already and govern the fragment formation processes. Therefore, we suggest that the exotic structures formation in the head-on heavy ion collisions can be explained from the hydrodynamic point of view. In order to explain the breakdown of such structures and to predict the fragments mass distribution it is necessary to consider the possible Rayleigh-Plateau instabilities, capillary waves and the influence of the long range Coulomb forces at the final stage of the system evolution. 4 Conclusion 1 Suggested approach allows for a simple physical picture of the exotic structure formation in the head-on heavy ion collisions. The qualitative picture obtained within our model is in a good correspondence with the one observed within the BUU calculations and in the experiment. The straightforward link between the EOS and the exotic shapes of different topology is explained from the hydrodynamic point of view. Combination of the introduced hydrodynamic approach together with the microscopic calculations can reveal the physical nature of the multifragmentation phenomena and give the possibility to extract the data on the EOS from the head-on heavy ion collisions. References Bertch G F and Das Gupta S . Phys Rep , 1988 , 160 : 189 - 233 . Phys Rep Abe Y , Ayik S , Reinhard P G , et al . Phys Rep , 1996 , 275 : 49 - 196 . Phys Rep Baran V , Colonna M , Greco V , et al . Phys Rep , 2005 , 410 : 335 - 466 . Phys Rep Li B A , Chen L W , Ko C M . Phys Rep , 2008 , 464 : 113 - 281 . Phys Rep Baldo M and Maieron C . J Phys G: Nucl Part Phys , 2007 , 34 : R243 - R283 . J Phys G: Nucl Part Phys Xie W J , Su J , Zhu L , et al . Phys Lett B , 2013 , 718 : 1510 - 1514 . Phys Lett B Kwiatkowski K , Fiedman W A , Woo L W , et al . Phys Rev C , 1994 , 49 : 1516 - 1524 . Phys Rev C Pan J and Gupta S . Phys Rev C , 1995 , 51 : 1384 - 1392 . Phys Rev C Zhang F S , Chen L W , Jin G M , et al . Phys Rev C , 1999 , 60 : 064604 . Phys Rev C Rodrigues M R D , Wada R , Hagel K , et al . J Phys: Conf Ser , 2011 , 312 : 082009 . J Phys: Conf Ser Hűfner J . Phys Rep , 1985 , 125 : 129 - 185 . Phys Rep De Souza R , Le Neindre N , Pagano A , et al . Eur Phys J , 2006 , 30 : 275 - 291 . Eur Phys J Fuchs C and Wolter H H . Eur Phys J A , 2006 , 30 : 5 - 21 . Eur Phys J A Zhang F S and Suraud E . Phys Lett B , 1993 , 319 : 35 - 40 . Phys Lett B Zhang F S and Suraud E . Phys Rev C , 1995 , 51 : 3201 - 3210 . Phys Rev C Tsang M B , Bougault R , Charity R , et al . Eur Phys J A , 2006 , 30 : 129 - 139 . Eur Phys J A Bauer W , Bertsch G F , Schulz H . Phys Rev Lett , 1992 , 69 : 1888 - 1891 . Phys Rev Lett Xu H M , Gagliardi C A , Tribble R E , et al . Phys Rev C , 1993 , 48 : 933 - 936 . Phys Rev C Stone N T B , Bjarki O , Gualtieri E E , et al . Phys Rev Lett , 1997 , 78 : 2084 - 2087 . Phys Rev Lett Ono A and Randrup J . Eur Phys J A , 2006 , 30 : 109 - 120 . Eur Phys J A Shlomo S , Kolomietz V M . Rep Prog Phys , 2005 , 68 : 1 - 76 . Rep Prog Phys Cherevko K V , Bulavin L A , Sysoev V M . Phys Rev C , 2011 , 84 : 044603 . Phys Rev C Su J and Zhang F S . Phys Rev C , 2013 , 87 : 017602 . Phys Rev C Wong C Y . Phys Rev Lett , 1984 , 52 : 1393 - 1396 . Phys Rev Lett Khodel V A . Sov J Nucl Phys , 1974 , 19 : 404 . Sov J Nucl Phys Stöcker H and Greiner W . Phys Rep , 1986 , 137 : 277 - 392 . Phys Rep Siemens P J and Rasmussen J O . Phys Rev Lett , 1979 , 42 : 880 - 883 . Phys Rev Lett Stocker H , Cusson R Y , Maruhn J A , et al . Z Phys A , 1980 , 294 : 125 - 135 . Z Phys A Baumgardt H G , Schott J U , Sakamoto Y , et al . Z Phys A , 1975 , 273 : 359 - 371 . Z Phys A Heynmann F J . J Appl Phys , 1969 , 40 : 5113 - 5122 . J Appl Phys Lesser M B . Proc R Soc Lond A , 1981 , 377 : 289 - 308 . Proc R Soc Lond A Haller K K , Poulikakos D , Venticos Y , et al . J Fluid Mech , 2003 , 490 : 1 - 14 . J Fluid Mech Aichelin J . Phys Rep , 1991 , 202 : 233 - 360 . Phys Rep Zhang F S . Z Phys A , 1996 , 356 : 163 - 170 . Z Phys A Zhang F S and Ge L X . Nuclear multifragmentation . Beijing : Science Press , 1998 . Nuclear multifragmentation Kutateladze S S and Styrikovich M A . Hydrodynamics of Gas-Liquid Systems , 1976 , Moscow : Energy Publ . Hydrodynamics of Gas-Liquid Systems Taylor G . Proc R Soc Lond A , 1959 , 253 : 296 - 312 . Proc R Soc Lond A 10.13538/j.1001-8042/nst.2013.05.015 Supported by National Natural Science Foundation of China (NSFC) projects (Nos.11025524 and 11161130520) and National Basic Research Program of China (2010CB832903) and Ministry of Education, Science Youth and Sports of Ukraine project (No. M/175-2012) . 2013-06-27 2013-10-01 2021-04-11 2013 Nuclear Science and Techniques 5 24