arxiv: v2 [nucl-ex] 25 Apr 2018

Transcription

1 Evidence for high excitation energy resonances in the 7 alpha disassembly of 28 Si arxiv:8.766v2 [nucl-ex] 25 Apr 28 X. G. Cao,,2 E. J. Kim, 2, K. Schmidt, 4,2 K. Hagel, 2 M. Barbui, 2 J. Gauthier, 2 S. Wuenschel, 2 G. Giuliani, 2,5 M. R. D. Rodriguez, 2,6 S. Kowalski, 4 H. Zheng, 2,5 M. Huang, 2,7 A. Bonasera, 2,5 R. Wada, 2 G. Q. Zhang,,2 C. Y. Wong, 8 A. Staszczak, 9 Z. X. Ren, Y. K. Wang, S. Q. Zhang, J. Meng,, and J. B. Natowitz 2 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 28, China 2 Cyclotron Institute, Texas A&M University, College Station, Texas 7784 Division of Science Education, Chonbuk National University, 567 Baekje-daero Deokjin-gu,Jeonju 54896, Korea 4 Institute of Physics, University of Silesia, 4-7 Katowice, Poland. 5 Laboratori Nazionali del Sud, INFN, via Santa Sofia, 62, 952 Catania, Italy 6 Instituto de Física, Universidade de São Paulo, Caixa Postal 668, CEP , São Paulo, SP, Brazil 7 College of Physics and Electronics Information, Inner Mongolia University for Nationalities, Tongliao, 28, China 8 Physics Division, Oak Ridge National Laboratory, Oak Ridge, USA 9 Institute of Physics, Maria Curie-Sk lodowska University, Lublin, Poland State Key Laboratory of Nuclear Physics and Technology, School of Physics, Peking University, Beijing 87, China Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto , Japan (Dated: April 26, 28) The excitation function for the 7 alpha de-excitation of 28 Si nuclei excited to high excitation energies in the collisions of 5 MeV/nucleon 28 Si with 2 C reveals resonance structures that may indicate the population of toroidal high-spin isomers such as those predicted by a number of recent theoretical calculations. This interpretation is supported by extended theoretical analyses. PACS numbers: 25.7.Pq Keywords: Intermediate energy heavy ion reactions, Clusters I. INTRODUCTION Light nuclei in their ground states in the valley of stability usually have spherical or near-spherical geometries []. With increasing excitation energy and/or angular momentum, these nuclei may exhibit more exotic shapes [2 5]. Theoretical investigations of the possible existence of extremely exotic nuclear shapes have a long history. Wheeler suggested that, under certain conditions, nuclei could assume toroidal shapes []. Pursuing this suggestion, Wong explored possible toroidal and bubble nuclei and predicted exited toroidal states in the mass region of 4 < A < 7 and A < 25 [6 8]. Large shell effects in light nuclei and large Coulomb energies in heavy nuclei favor exotic toroidal configurations. A recent search for heavy toroidal systems indicated that the probability of planar fragmentation configurations in the experimental data was much greater than predicted by quantum molecular dynamics calculations [9]. Extending such studies Wong and collaborator predicted that, toroidal configurations were also possible for nuclei with a sufficiently high angular momenta [ 2]. They defined the region of mass and angular momentum in which such configurations might be realized. More recent theoretical studies have used microscopic techniques to address this question of light toroidal nuclei. In particular, Zhang et al. [], Ichikawa et al. [4, 5] and Staszczak and Wong [6 8], using varied approaches, have predicted the existence of toroidal isomers in light nuclei. We can understand the origin of possible light-a toroidal isomers in the following simple way. In a nucleus with a toroidal shape, there are toroidal magic numbers 2(2m + ), with integer m, arising from large energy gaps between single-particle levels in the light mass region (Fig. of [6]). Extra stability [9] associated with toroidal magic numbers leads to an excited local energy minimum that is stable against the expansion and contraction of the toroidal major radius (see Fig. 2 of [6], Fig. 8 of [7], and Fig. of []). Such an excited state residing in a local energy minimum under a toroidal shape constraint will be called a diabatic state, and the corresponding constrained calculation a diabatic calculation []. This is in contrast to an adiabatic state of the lowest energy minimum in an adiabatic calculation without a shape constraint [4 8]. Relative to a diabatic toroidal local energy minimum core, Bohr- Mottelson spin-aligning particle-hole excitations [2] can be constructed to yield a toroidal nucleus with a spin, I=I z, by promoting nucleons with angular momentum aligned opposite to a chosen symmetry axis to populate orbitals with angular momentum aligned along the symmetry axis [4 8]. The spinning toroidal nucleus possesses an effective rotational energy that tends to expand the toroid, whereas the energy associated with the nuclear bulk properties tends to contract the toroid.

2 2 TABLE I: Calculated parameters for the predicted toroidal isomer in 28 Si. Left to right: Spin I=I z, quadrupole moment Q 2, the cranking rotation frequency hω, excitation energy E x, toroidal major radius R, minor radius d, aspect ratio R/d, and maximum density ρ max I Q 2 hω E x R d R/d ρ max [ h] [b] [MeV] [MeV] [fm] [fm] [fm ] 28 Si The balance between the two energiesgives rise to a local toroidal energy minimum [8]. For small values of I, the toroidal minimum occurs as an excited (diabatic) state above the sphere-like ground states. As I increases, the crossing of the toroidal and sphere-like energy surfaces takes place, and the toroidal high-spin energy minimum switches to become the lowest energy (adiabatic) state. Adiabatic self-consistent calculations have located an extensive region of toroidal high-spin isomer (THSI) states that are stable against the expansion and contraction of the toroids [4 8]. Motivated by the predictions of toroidal isomeric states in a number of light nuclei, we have undertaken searches for evidence of their existence. For 28 Si, the nucleus investigated in the present work, the Staszczak and Wong calculations indicate that toroidal shapes with I= become possible at excitation energies greater than 85 MeV. The existence of a stabilized state with angular momentum of 44 h and excitation energy of 4.8 MeV [6] was predicted. See Table. In the calculation the toroidal shape of this short-lived isomer 28 Si is characterized by a radius of 4. fm and a cross sectional radius of.45 fm for the cylindrical ring containing the nucleons. Thus the aspect ratio is The predicted toroidal states, although expected at very high energies, are analogous to yrast traps already observed in more conventionally shaped nuclei [22]. Should such very highly excited stabilized toroidal states of light nuclei exist, their lifetimes should be short. They may de-excite or undergo shape relaxation rather quickly. In either case the most-likely de-excitation modes are particle or cluster emission and fragmentation [2]. It is well documented that macroscopic toroids fragment as a result of the development of Plateau-Rayleigh instabilities [24 28]. In the basic Rayleigh description [25] the dominance of a single mode of symmetric fragmentation leads to disassembly into equal size pieces, the number of which is of the order of the aspect ratio of the toroid [24 28]. This correspondence is most accurate for large aspect ratios [27]. Modern numerical simulations taking into account the viscosity of the fluid and the surrounding medium indicate that more complicated symmetric breakups involving different size fragments are possible [28]. As has been already discussed in the literature, nuclear tori might also manifest Plateau-Rayleigh instabilities [6, 7, ]. The aspect ratio predicted by Staszczak and Wong for the 4.8 MeV, 44 h state in 28 Si suggests that, while other fragmentations are possible [27], the dominant instability of that toroid would lead to a break-up into fragments. However the actual dominant mode will be affected by the temperature dependent viscosity of the disassembling nucleus [29]. In the nuclear case the discreteness of the nucleons, the existence of Coulomb forces, shell effects and variations in the fragment binding energies may also modify the fragmentation pattern of the torus. While strongly reduced Coulomb energies might be expected to provide a signature for the disassembly of predicted heavy or super-heavy nuclear toroids in their meta-stable ground-states [8, 9], this is probably not so for the predicted high spin light tori. These light tori are predicted to have very high excitation energies. Thus, in their disassembly, a lowering of the Coulomb repulsion between fragments is probably not observable as the release of large deformation and rotational energies will normally lead to large kinetic energies for the observed decay fragments. Finally, given the high excitation energies involved, it is likely that the initial fragments will often be excited and undergo subsequent de-excitations, smearing the signature of a Plateau-Rayleigh instability. All of these considerations suggest that judicious choices of reaction mechanism, exit channels and observables will be necessary to probe the possible existence of these very exotic and very interesting nuclei. Two notable features of the Staszczak-Wong and Ichikawa calculations are that. the cross-sectional radii of the cylindrical rings containing the nucleons are.5 fm, essentially equal to the alpha particle radius. Indeed, in their search for isomeric states in 4 Ca, Ichikawa et al. used a ring of alpha particles as an initial configuration for their cranked Hartree Fock calculations [4]. 2. the toroidal rings corresponding to the predicted isomers have matter densities approximately 2/ of ρ where ρ is the normal central density of such nuclei [5 7]. Since the general importance of α-like correlations in the structure and properties of light nuclear systems at normal and reduced densities is now well established [ 5], the features noted here suggest that the disassembly into alpha particles or alpha-conjugate nuclei might be favored for light nuclear toroids. In 986 Wilkinson specifically suggested that spinning rings of alpha particles might be stabilized by circumnavigating neutrons. He discussed their stabilities toward electromagnetic or fragmenting de-excitations and the possibility of producing them in heavy ion collisions [6]. The idea of stabilization by surrounding neutrons is also a feature of the extended Ikeda diagram systematics [7].

3 We have searched for toroidal isomers in the disassembly of 28 Si produced in near-fermi-energy collisions of 28 Si with 2 C. While some nucleon transfer and nucleon-nucleon collisions leading to early (preequilibrium)particle emission occur [8 44], such collisions, many of which lead to essentially binary exit channels, are capable of producing projectile like nuclei in the A=28 mass region with high excitation energy and high angular momentum. Calculations[4 44]indicate that an angular momentum in the range of 4 h and an excitation energy as high as 7 MeV can be reached. These calculations should be considered as only indicative of possible angular momentum range as they do not have the ingredients to explore detailed quantum structure at such high excitation and angular momentum. When the excitation energy E and the angular momentum I of an emerging 28 Si corresponds to that of a toroidal high-spin isomer, the collective cranking motion and the rearrangement of the single-particle motion of the nucleons may eventually lead to the toroidal high-spin isomer. Here we report results for an investigation of 28 Si, focusing on the 7α decay channels of excited projectile-like fragments produced in the reaction 5 MeV/nucleon 28 Si + 2 C. In this reaction the energy available in the center of mass is 294 MeV, the maximum angular momentum, L max, is 94 h (a reaction cross section of 247 mb), L crit for fusion is 26 h and the rotating liquid drop limiting angular momentum is 4 h [45]. These parameters indicate that the bulk of the reaction cross section will lead, not to fusion, but to initially binary configurations of excited projectile-like and target-like nuclei. This is consistent with experimental results reported for similar collisions in this energy region [8, 9]. EXPERIMENTAL DETAILS The experiment was performed at Texas A&M University Cyclotron Institute. A 5 MeV/nucleon 28 Si beam produced by the K5 superconducting cyclotron impinged on a 2 C target. The reaction products were measured using a 4π array, NIMROD-ISiS (Neutron Ion Multidetector for Reaction Oriented Dynamics with the Indiana Silicon Sphere)[46, 47] which consisted of 4 concentric rings covering from.6 to 67 in the laboratory frame [46]. In the forward rings with θ lab 45, two special modules were set having two Si detectors (5 and 5 µm) in front of a CsI(Tl) detector ( cm), referred to as super-telescopes. The other modules (called telescopes) in the forward and backward rings had one Si detector (one of 5, or 5 µm) followed by a CsI(Tl) detector. The pulse shape discrimination method was employed to identify the light charged particles with Z in the CsI(Tl) detectors. Intermediate mass fragments (IMFs), were identified with the telescopes and super-telescopes using the E E method. In the forward rings an isotopic resolution up to Z = 2 and an elementalidentificationup toz = 2wereachieved. In the backward rings only Z = 2 particles were identified, because of the detector energy thresholds. In addition, the Neutron Ball surrounding the NIMROD-ISiS charged particle array provided information on average neutron multiplicities for different selected event classes. Further details on the detection system, energy calibrations, and neutron ball efficiency can be found in [8, 47, 48]. It is important to note that, for symmetric collisions in this energy range, the increasing thresholds with increasing laboratory angle lead to a condition in which the efficiencies strongly favor detection of projectile like fragments from mid- peripheral events. Modeling of these collisions using an Antisymmetrized Molecular Dynamics (AMD) code [4] and applying the experimental filter demonstrates that this is primarily an effect of energy thresholds. Nucleons are emitted during the initial phase of the collision. They are often modeled as emission from a virtual mid-rapidity source having a velocity close to that of the nucleon-nucleon collision frame [4]. While these nucleons, ejected prior to equilibration of the remaining system, are not from the de-excitation of the primary exit channel products, some of them appear in the projectile velocity region. They should not be included in the calorimetric determination of the thermalized excitation energy. ANALYSIS Forthe 28 Si+ 2 Creaction,atotalof7.5millionevents were recorded and a significant proportion of events had significant alpha-like (AL) mass emission(i.e. alpha particles or alpha-conjugatenuclei). About.9 5 events had AL=28. Of these 82 detected events had 7 alpha particles. The AMD calculations indicate that for detected A=28 channels 9% of the nucleons originate from the projectile. Some mixing is consistent with damped collisions. A very careful inspection of experimental invariant velocity plots for each reaction channel and each emitted species confirmed that, projectile like sources strongly were strongly dominant in the detected events, except for Z= particles. These Z= particles showed some contribution from both target-like and mid-rapidity sources, the latter characteristic of pre-equilibrium emission. The observed source velocities for the selected A=28 events decrease slowly with exit channel complexity but are always above the center of mass velocity. Complete fusion is a rare process and target like fragments are generally not detected as indicated above. Searching for 8 alpha events we found that their yield was a factor of less than for 7 alpha events. This indicates small possibility of contamination from such events. All of these features

4 4 were taken into account in the following analysis. To characterize the source excitation energies involved we have used calorimetric techniques to determine the excitation energy, E x, of the primary projectile-like fragments. E x, is normally defined as the sum, of the kinetic energies of ejected particles and fragments in the frame of the total projectile-like nucleus, (determined by reconstruction of the mass and velocity of the primary excited nucleus from its de-excitation products) minus the reaction Q-value. See equation (). E x = Σ Mcp i= K cp(i)+m n < K n > Q () where M cp is the multiplicity of charged particles, K cp is the kinetic energy of a charged particle in the source frame, M n is the neutron multiplicity and < K n > is the average neutron kinetic energy in the source frame. The events initially selected typically had a few Z= particles (and neutrons) and, in very rare cases, a heavier fragment, associated with them. In our event selection, which focusses on the alpha-conjugate exit channels, we have allowed Z= particles and neutrons in the events but we have determined the excitation energies excluding the energies of Z= particles and neutrons as they are primarily pre-equilibrium particles, representing energy dissipation but not energy deposition into the projectile-like fragment [4]. Invariant velocity plots for the alpha particles indicated that the alpha particles resulting from mid-rapidity pre-equilibrium emission were negligible and that the small number detected from the target like source (given the thresholds and geometry of the NIMROD detector) could be effectively removed by rejection of alpha particles with energies greater than 4 MeV in the projectile-like fragment source frame. In Figure we present the excitation functions derived in this manner for the alpha-conjugate exit channels. For comparison, results from calculations using the phenomenological event generator HIPSE [44] and the AMD transport model with a Gemini afterburner are also presented [4, 42]. While the agreement between the data, the AMD and the HIPSE results is generally good, we note that the 7-alpha exit channel results differ considerably. The experimental 7-alpha results appear to have structure at the higher excitation energies. In the following we focus on this exit channel. Two hybrid codes, AMD-GEMINI [4, 42, 4] and HIPSE-GEMINI [42 44] were also used to calculate the 7-alpha excitation energy spectrum. The modeling of light particle emission using the code GEMINI was thoroughly examined in reference [4]. In that work, formulations of the barriers and transmission coefficients, the level density and the yrast line were carefully explored. We have used the default parameter prescriptions derived from this work. Results of both calculations, filtered through the NIMROD acceptance and normalized to the data, are also presented in the figure. It is not expected that high energy resonances could appear in these models and, indeed, both are structure-less. The AMD results are somewhat broader in energy and shifted to slightly higher energy than the HIPSE results. As will be seen in this figure, the 7 alpha distribution spans the energy region in which toroidal configurations are predicted and the 4.8 MeV stabilized state is predicted to exist. After reachinga maximumat MeV the excitation energy distribution shows some structure at 26 and 8 MeV. The granularity and angular resolution of NIMROD-ISiS are not ideally suited to searches suchasthis, asthe transformationtothe sourceframerelies upon the angle of detection. Through simulations we have determined that the observed experimental width in excitation energy of an initially sharp state at 4 MeV will have a standard deviation, σ, of 4 MeV resulting from the angular uncertainty. Taking this into account, the broad structures apparent in the excitation energy spectrum are consistent with much narrower resonances in the excitation energy distribution. We have checked this by adding a 7 alphas delta function with Ex=4 MeV to the uncorrelated 7 alpha events and observing that the resultant filtered spectrum is consistent with our observed spectrum. Clearly an experiment with much better angular resolution, allowing better resolution for the excitation energy spectrum, would be very desirable. In figure 2, on the left hand side, we compare the experimental spectrum for the 7 alpha events to an uncorrelated 7 alpha spectrum. The spectrum of uncorrelated events was constructed by randomly selecting 7 alpha particles from 7 different events. The random selection was done many times to assure that the statistical fluctuations for this spectrum would be much lower than those of the correlated event spectrum. This uncorrelated spectrum is taken to represent the 7 alpha phase space in this excitation energy region. Subtracting the normalized uncorrelated spectrum from the data results in the difference spectrum depicted in the right of Figure 2. Some excess is seen at 4 MeV. The peaks at 26 and 8 MeV are quite prominent and there is a tailing toward higher energies. The under-shoot at MeV and below may suggest that the normalization of the uncorrelated spectrum is too conservative. If that spectrum is lowered the peaks in the difference spectrum will be even more prominent. For comparison we have also employed the AMD calculated spectrum as a background spectrum to be subtracted from the experimental spectrum. For this purpose the AMD-GEMINI spectrum was shifted to agree with the experimental uncorrelated spectrum at the lower edge of the experimental spectrum. The difference spectrum obtained by subtracting the normalized AMD- GEMINI spectrum from the experimental data, is very similar to, but not exactly the same as that obtained when subtracting the experimental uncorrelated spectrum from the experimental data. The structure in the experimental 7α spectrum ap-

5 Mg+α Ne+2α O +α C +4α α C+α O + C Data AMD+Gemini HIPSE+Gemini FIG. : Excitation functions for the alpha conjugate exit channels in the de-excitation of 28 Si. The shapes of the experimental data are compared with results of both AMD and HIPSE calculations. pears to reflect mechanisms not encompassed in the dynamic reaction models or normal statistical decay treatments. Relative to the uncorrelated background derived from the experiment the statistical significance of the difference peak at 4 MeV is 5.σ, at 26 MeV is 8.σ and at 8 MeV is 7.2σ [49]. We take these to be the minimum values for the statistical significance since the construction of the uncorrelated spectrum includes contributions from the peak region and the number of uncorrelated alphas may, therefore, be overestimated in that region. Evaluation of the statistical significance of the observed peaks is sensitive to the background assumed. If we base the test for the statistical significance on the use of the AMD GEMINI result we find a statistical significance of 4, 26 and 8 MeV peaks to be 4.2σ, 6.σ, and 6.6σ, respectively. CROSS SECTION AND ANGULAR MOMENTUM To determine the cross section for the 7 alpha channel we have assumed that the total events detected with our minimum bias trigger ( particle or fragment detected), corrected for detector efficiency, represent a total reaction cross section of 247 mb [45]. The overall detector efficiency was determined from the ratio of numbers of AMD-GEMINI generated events before and after the detector filter. The specific efficiency for the 7 alpha 5 45 Data Mixed 4 event data (a) 5 AMD+Gemini simulation Data - (mixed event data) Data - (AMD+GEMINI simulation) (b) FIG. 2: Excitation energy distribution leading to observed 7-alpha events. Left panel-the data are represented by a solid black line. An uncorrelated spectrum derived from event mixing is represented by a long dashed line. The filtered result from an AMD-GEMINI calculation is indicated by the dashed and dotted line (see text). The last two are normalized to the experimental spectrum at the lower edge of the spectrum. On the right the differences between the experimental spectrum and the others are presented. Relative to the uncorrelated background derived from the experiment the statistical significance of the difference peak at 4 MeV is 5.σ, at 26 MeV is 8.σ and at 8 MeV is 7.2σ. See text channel was also determined in a similar fashion. Because double hit corrections can be large, these results are very sensitive to the number of, and excitation energies of intermediate 8 Be nuclei produced. Thus there may be a significant systematic uncertainty in the resultant cross section. Based on the estimated uncertainties

6 6 we estimate the cross section for the 7 alpha channel to be.9 mb with a systematic uncertainity of 25%. We estimate the cross sections for the 26 and 8 MeV peaks, respectively, to be 5µb and 28µb with similar uncertainities. If the state observed at 8 MeV corresponds to the predicted 4.8 MeV toroidal state the angular momentum would be 44 h. Our results give no direct information on the angular momentum. As already noted, the AMD and HIPSE calculations, employing semi classical techniques, indicate that angular momenta in the range of 4 h are reached but they do not have the ingredients necessary to explore detailed structure at such high excitation energy and angular momentum. MOMENTUM SPACE SHAPE ANALYSES We have utilized a shape analysis technique [5, 5] to diagnose the momentum space source shape for the 7α events. This type of analysis is a popular method to study emission patterns of sources, dynamical aspects of multi-fragmentation and collective flows of particles in heavy ion collisions at intermediate and relativistic energies. Although we are not able to observe directly the geometric shape of these sources, the momentum space correlations among the detected fragments provide hints of shapes and information on the disassembly dynamics. Should nuclei with exotic shapes undergo simultaneous fragmentation into equal sized cold fragments the observed momentum space distributions of products would be directly related to the initial geometric configuration of the de-exciting system. In contrast, if the primary fragments are excited or the emission is sequential, deexcitation of the initially produced fragments could significantly modify the initial momentum space distribution. For the light alpha-conjugate nuclei explored in this study, alpha emission from excited heavier primary fragments would produce large perturbations of that distribution. The analysis employs a tensor constructed from the product momenta, written as: T ij = Σ N ν= pν i pν j where N is the total nucleon number, p ν i is the momentum component of ν th nucleon in the center-of-mass and i refers to the Cartesian coordinate. The tensor can be diagonalized to reduce the event shape to an ellipsoid. The eigenvalues of the tensor: λ, λ 2 and λ, normalized by: λ +λ 2 +λ = and orderedaccording to: λ λ 2 λ, can quantitatively give shape information of the events. The sphericity is defined as: S = 2 ( λ ), and coplanarity is defined as: C = 2 (λ 2 λ ). In the sphericitycoplanarity plane, the ideal rod, disk and sphere events are exactly located at the three vertices of the triangle: (,) (/4, /4) and (,), respectively. Two fragment events will always appear at, while three fragments events will appear on the rod-disk axis. To appear as a disk at the apex of this triangle requires a simultaneous and symmetric fragmentation into or more equally sized pieces. The results of the momentum shape analysis are shown in Fig. (a). While there are some events located in the disk region, the bulk of the events are elsewhere. Also depicted in Figure are sphericity-coplanarity plots of the AMD-GEMINI results. Figure (b) portrays the filtered final results. Figure (c) presents the plot for the fm/c freezeout momentum distribution predicted by AMD. Clearly this freeze-out distribution is much more rod to disk like than that observed after de-excitation and filtering. The model results suggest that the observed final distribution normally results from processes in which an initial breakup into larger excited fragments is followed by alpha particle de-excitation of those fragments. COMPARISON WITH TOROIDAL SHELL MODEL The 7α de-excitation of 28 Si projectile-like fragments encompasses the excitation energy range for which toroidal configurations are predicted to exist. The resonance structure indicates the existence of stabilized configurations which may correspond to toroidal high-spin states such as those predicted in several theoretical calculations [ 8]. If the observed state at 8 MeV corresponds to the 4.8 MeV state predicted to exist by Staszczak and Wong [6], the angular momentum would be 44 h. Resonances at such high excitation energies are rather unusual,andweknowofnomodelexceptthetoroidalisomer model that predicts 28 Si resonances at these high excitation energies. We have explored whether the experimental data may be described in terms of the toroidal shell model in which nucleons move in the toroidal potential V(ρ,z) = 2 mω2 (ρ R) mω2 z 2 [6, 7, 6]. Upon neglecting the small spin-orbit interaction, the toroidal single-particle energy ǫ(nλω) for the nλω state, in the I= toroidal core with a major radius R and R d, is approximately ǫ(nλω) hω (n+)+ h2 Λ 2 2mR2, (2) where n=n z +n ρ is the harmonic oscillator quantum number, ω the oscillator frequency, m the nucleon mass, Λ the orbital angular momentum about the symmetry z- axis, Ω z =Λ+s z, and s z the intrinsic nucleon spin. Relative to the toroidal I= core at E occupying the lowest toroidal single-particle states, the spinaligning Bohr-Mottelson particle-hole excitations leading to toroidal high-spin isomers can be constructed by following the crossings of Routhian energy levels as a function of increasing cranking frequency hω (Fig. (b) of

7 7.5.4 Data (a).5.4 AMD+Gemini (b) (c) Coplanarity..2. Coplanarity..2. Coplanarity Sphericity Sphericity Sphericity FIG. : Shape analysis of alpha-conjugate exit channels in the de-excitation of AL = 28, Z = 4 projectile-like nuclei produced in the 5 MeV/nucleon 28 Si + 2 C reaction. Left-Experimental data. Middle-Filtered results from AMD-GEMINI calculation. Right- Results for AMD primary fragments at fm/c. [6]). The contribution to I z and (E I E ) from a particle-hole excitation can be easily obtained from the changes in particle-hole state quantum numbers (nλω) and ǫ(nλω). For 28 Si, for example, the p-h excitation contributes I z =8 h, and (E I E )=7 h 2 /2mR 2 bypromoting a nucleon from ( 7/2) to 4(9/2). Similar contributions can be obtained for the 2p-2h excitation by additional promotion from ( 5/2) to 4(7/2)), p-h from 2( 5/2) to 5(/2), and 4p-4h from 2( /2) to 5(9/2). From such calculations, we obtain the spin I=I z and the relative energy (E I E ) (in units of ( h 2 /2mR 2 )), for various 28 Si toroidal high-spin isomer states as the signature for toroidal 28 Si in Table II. TABLE II: Toroidal high-spin isomers (THSI) of 28 Si in the toroidal shell model. The spin-aligning (n particle)-(n hole) excitations, for neutrons (ν) and protons (π), relative to a toroidal core with I= and energy E, lead to the THSI state of spin I=I z, and excitation energy E I. Configurations I (E I E ) E I in h 2 /2mR 2 (MeV) (p-h) ν(p-h) π 9.82 (p-h) ν(p-h) π (p-h) ν(2p-2h) π+(2p-2h) ν(p-h) π (2p-2h) ν(2p-2h) π (2p-2h) ν(p-h) π+(p-h) ν(2p-2h) π (p-h) ν(p-h) π (p-h) ν(4p-4h) π+(4p-4h) ν(p-h) π (4p-4h) ν(4p-4h) π We can identify the I=44 h, 6 h and 28 h toroidal highspin isomers in the toroidal shell model as the resonances E A =8.7 MeV, E B =25.4 MeV, and E C =2.7 MeV, respectively, with a high degree of confidence based on the following grounds: (i) the theoretically predicted I=44 h state at E 44 = 4 MeV in Ref. [6] is close to the observedexperimentalresonanceenergye A =8.7MeV, (ii) the theoretical energy spacing between the I=44 h an I=6 h THSI states and between the I=6 h and the I=28 h THSI states are equal, and likewise the experimental spacing between resonances E A and E B and betweenresonancese B ande C arealsoapproximateequal, and (iii) above the I=44 h THSI state, there exist, theoretically, additional I=5 h and 56 h THSI states at higher energies, and there remains significant experimental cross section above the resonance E A =8.7 MeV. Expressing E I as E I = E +[(E I E )/( h 2 /2mR 2 )] ( h 2 /2mR 2 ), () we find by using the quantity [[(E I E )/( h 2 /2mR 2 )] in TableII thate =9.82MeVandR=5.56fm, withwhich we can determine E I of all other THSI states presented in Table II. Table II reveals that up to 4p-4h excitations, the toroidal shell model predicts THSI yrast states built on the toroidal I= core. They occur within the range of excitation energies of the present measurement. Note that some of the states are doubly degenerate. Based upon the above THSI energy levels E I as primary ingredients, we have constructed a simple phenomenological semi-empirical formula to estimate the relative THSI production cross sections by including the width parameters σ I. We start by assuming that the 7α cross section obtained after subtracting the uncorrelated cross section from the correlated cross section in Fig. 2(b) arises dominantly from toroidal configurations. We assume further that the reaction cross section above the observed threshold is proportional to the distribution of deposited angular momentum I which is governed by the impact parameter and the collision dynamics. For toroidal THSI production, the energy and angular momentum must match that of a THSI state. Hence the sum over I with (di)= is carried out for THSI states. Neglecting other unknown factors, we developed the following semi-empirical cross section formula for the 7α

8 8 channel from toroidal configurations σ toroidal (E x,7α) = A I=I toroid g I I +exp{(i I max )/a} exp{ (E x E I ) 2 /2σI 2 }, (4) 2πσI where gi is the state degeneracy factor and I max and the diffusion parameter a are introduced phenomenologically to describe initial-state dynamical and/or finalstate structural limitations. The most important primary ingredients in the above formula are the THSI spin I, energies E I, and degeneraciesg I. Upon using the toroidalshell model oftableii to fix these primary ingredients, we extract the secondary quantities of the widths and other parameters. We find that the gross features of the excitation function can be well described by the semi-empirical formula (4) with extracted widths and fitting parameters as shown in Fig 4. The extracted widths for the sharp I=44 h, 6 h and 28 h resonances are small, of the order of the experimental bin size. The extracted widths for the I=5 h and 56 h states are large, which may indicate that the particle-hole excitations for these highest lying THSI states may involve promoting particles to populate states with large intrinsic widths close to the particle drip line. However, low statistical uncertainities and uncertainities in background in that region may contribute to the apparent broadening. The present analysis indicates that the a spinning toroidal 28 Si expresses itself as a set of THSI states with a unique signature listed in Column of Table II. In the excitation energy spectrum we find different facets of the toroidal signature,i.e, the presence of sharp resonances at appropriate energies, the spacing between some of the resonances, the apparent presence of THSI states at energies higher than the predicted I=44 h state. The approximate matching of the experimental excitation energy spectrum with the semi-empirical cross section formula containing all THSI states up to 4p-4h excitations as exhibited in Fig. (4), provides strong evidence for the possible production of toroidal high-spin isomers in the present experiment. ADDITIONAL THEORETICAL SUPPORT FROM RELATIVISTIC MEAN-FIELD CDFT THEORY Our experimental results and the phenomenological analysis suggest the existence of more than one stabilized state. Previous non-relativistic adiabatic calculations without a toroidal shape constraint give only the I=44 h state [6], and cannot be used to calculate diabatic THSI states such as those with I<44 h. The phenomenological description of the toroidal shell model predicts a total of THIS states built on the toroidal I= core. These toroidal isomers states can be searched for counts 8 Data - (mixed event data) Fit A= 7.,±.4 Imax= 57.9,± 2.4 a I = 6.,±.6 σ 4,6 = 6.7,± 2.52 σ 28 = 2.5,±. σ 6 =.64,±.6 σ 44 =.4,±.74 σ 5 =.2,±.7 σ 56 = 5.25,± FIG. 4: The experimental (correlated data) - (mixed event data) (in solid points) are compared with the results of semiempirical formula (4) containing the signature of the toroidal high-spin isomers as primary ingredients (in solid curve) with the list of extracted widths and other parameters. and examined by the covariant density functional theory (CDFT) [52], which exploits basic properties of QCD at low energies, in particular, symmetries and the separation of scales. It is worth pointing out that the CDFT theory has provided an excellent description of ground states and excited (diabatic) states for nuclei all over the periodic table with a high predictive power [5 55]. Using a universal density functional and without assuming the existence of clusters a priori, CDFT provides a high degree of confidence in the investigation of nuclear toroidal structures. With the most successful density functionals PC- PK [2] and DD-ME2 [56], the newly-developed cranking CDFT in D lattice space [57, 58] has been applied to investigate the toroidal states in 28 Si. In these calculations, the z axis is chosen as the symmetry axis. Grid points are respectively taken for x, y and z with a step size.8 fm. The self-consistency of calculations is achieved with an accuracy of 4 MeV for the single-particle levels. The pairing correlations are neglected. By choosing a trial initial wave function with ring-like configuration of seven alpha-particles on the plane with z =, a toroidal state with I = 44 h is obtained at rotational frequency hω = 2.5 MeV, corresponding to p-h configurations for both neutrons and protons. The excitation energy of this toroidal state is 47.9 MeV for PC-PK and 45.8 MeV for DD-ME2. Local energy minima are found for the toroidal state with I=28 h, corresponding to 2p-2h excitations for both neutrons and protons, and the toroidal state with I=6 h,

9 9 Density [fm - ] (a) (b) I =28, PC-PK Total neutron density Total neutron density Neutron single-particle density r [fm] FIG. 5: Neutron radial density distributions (z direction is integrated) of the occupied single-particle levels (blue and thin lines) as well as the total density distribution (black and thick lines) in toroidal state with I = 28 h. corresponding to the 2p-2h excitation for neutrons and p-h for protons. The toroidal state with I = 44 h can be easily obtained by the adiabatic CDFT calculations. For the toroidal states with I = 28 h and 6 h, diabatic or configuration-fixed CDFT calculations are necessary. Additional toroidal states with I= h, 4 h, 6 h, 5 h and 56 h have also been located. They are shown in Table III and IV, respectively. In the column of hω, the mid-points of a range of the cranking rotational frequencies are given for toroidal diabatic (excited) state configurations arising from particle-hole excitations without specifying hω. In order to investigate the stability of these highspin torus isomers against particle emission, we have examined the radial density distributions of the occupied single-particle levels as well as the total density distributions for these toroidal states. As an example, neutron densities in toroidal state with I = 28 h are shown in figure 5. It can be clearly seen that all radial density distributions are well localized and the stability of the toroidal isomer against particle emission is hereby demonstrated. Similar conclusions can be drawn for other toroidal states as well. Results from the CDFT theory confirm the previous theoretical result predicting a toroidal high-spin isomer state with I=44 h at an excitation energy of E x MeV, and support the identification of the resonance observed experimentally at E x =8 MeV as this possible toroidal state. In addition all THSI states obtained in the toroidal shell model have also been located in the relativistic mean-field CDFT theory, supporting the use TABLE III: Properties of toroidal states in 28 Si obtained by the covariant functional PC-PK [2]. In the table, hω is rotational frequency, β 2 is quadrupole deformation, related to the quadrupole moment Q 2 of [6] by β 2= Q 2 /[(/ 5π)A 5/ r] 2 with r =.2 fm, E is excitation energy, and R, d, and ρ max are determined by fitting the calculated density distributions to the Gaussian function ρ(x,y,z) = ρ maxe [( x 2 +y 2 R) 2 +z 2 ]/(d 2 /ln2). I Configuration hω β 2 E R d ρ max [ h] [MeV] [MeV] [fm] [fm] [fm ] (ph) ν(ph) π (ph) ν(2p2h) π (2p2h) ν(ph) π (ph) ν(ph) π (2p2h) ν(2p2h) π (2p2h) ν(ph) π (ph) ν(2p2h) π (ph) ν(ph) π (ph) ν(4p4h) π (4p4h) ν(ph) π (4p4h) ν(4p4h) π TABLE IV: Same as Table. III but for the covariant functional DD-ME2 [56]. I/ h Configuration hω β 2 E R d ρ max [MeV] [MeV] [fm] [fm] [fm ] (ph) ν(ph) π (ph) ν(2p2h) π (2p2h) ν(ph) π (ph) ν(ph) π (2p2h) ν(2p2h) π (2p2h) ν(ph) π (ph) ν(2p2h) π (ph) ν(ph) π (ph) ν(4p4h) π (4p4h) ν(ph) π (4p4h) ν(4p4h) π of the THSI states in the toroidal shell model as the signature for toroidal high-spin isomers. There is however one notable difference that may reveal new physics associated with the toroidal shape. The spacings between energy levels in the relativistic CDFT theory appear to be significantly greater than their corresponding spacings in the toroidal shell model, or in the experimental data. One interesting possibility may be that the toroidal THSI nuclei under consideration have such a distorted shape and low densities that they may probe the nuclear energy density functional in a new regime for which the extrapolation from the normal nuclear matter in the CDFT theory may not be adequate. The possibility of using toroidal nuclei to probe the nu-

10 clear density functional at lower density may add an interesting dimension to the study of toroidal high-spin isomers. SUMMARY AND CONCLUSIONS In conclusion, the excitation function of the 7α channel in deep-inelastic collisions of 28 Si on 2 C reveals resonance structures at high excitation energies. The features of these structures appear to coincide with those predicted by a number of theoretical calculations in which the toroidal shell effects stabilize the nucleus in a diabatic I= state against major radius variations while spinaligning particle-hole expectations lead to many highspin toroidal isomers. From the theoretical perspective, if the present results are confirmed by further studies, a very large number of diabatic and adiabatic toroidal high-spin isomers in a very large light-mass region may be opened up for future investigations. Finally, we note that recent experimental and theoretical works provide indications that clustering effects are important in the collisions of alpha conjugate nuclei [59, 6]. The latter work, published as this paper was being prepared for submission provides theoretical evidence that toroidal alpha substructures can be quite commonly produced in such collisions. We strongly encourage further experimental work on collisions of light alpha-conjugate systems, both for the production of exotic clustered states and for the investigation of the dynamical evolution during such collisions. A higher granularity detector system and addition of gamma ray detection could offer significant improvements for such studies. ACKNOWLEDGEMENTS This work was supported by the United States Department of Energy under Grant # DE-FG-9ER477 and under grant # DE-AC5-OR22725 with UT- Battelle, LLC (Oak Ridge National Laboratory) and by The Robert A. Welch Foundation under Grant # A. Partial support by the National Natural Science Foundation of China under Contracts No. 4255, No. 52, No. 62, and No. 529 and the Youth Innovation Promotion Association CAS (No. 279) are acknowledged. Travel support for C. Y. Wong under the CUSTIPEN Program is thankfully acknowledged. We appreciate useful communications from A. Ono, J. Maruhn, T. Ichikawa and S. Umar. We also greatly appreciate the efforts of the staff of the TAMU Cyclotron Institute. [] N. J. Stone, Atom. Dat. Nucl. Dat. Tab. 9, 75 (25). [2] S. Cohen, F. Plasil and W.J. Swiatecki, Ann.Phys.(NY) 82, 557 (974). [] J. A. Wheeler, Nucleonics Notebook, 95(unpublished), see also p. 297 in G. Gamow, Biography of Physics, Harper & Brothers Publishers, N.Y. 96; Princeton University Graduate Course Physics 576 Take-Home Examination Problem 2, May 22, 96 (unpublished). [4] R. A. Brown and L.E. Scriven, Proc. R. Soc. Lond. A 7, (98). [5] S. Åberg, H.Flocard and W. Nazerewicz, Ann.Rev.Nuc. Part. Sci, 4, 49 (99). [6] C. Y. Wong, Phys. Lett. B4, 446 (972). [7] C. Y. Wong, Ann. Phys. (NY) 77, 279 (97). [8] C. Y. Wong, in SuperheavyElements, ed. M. A. K. Lodhi (Pergamon Press,New York, 978) p [9] R. Najman et al., Phys. Rev. C 92, 6464 (25). [] C. Y. Wong, Phys. Rev. C 7, (978). [] C. Y. Wong, Phys. Rev. Lett. 55, 97 (985). [2] J. Y. Zhang and C. Y. Wong, Phys. Rev. C 4, 94 (986). [] W. Zhang, H. Z. Liang, S. Q. Zhang, J. Meng, Chin. Phys. Lett. 27, 2 (2). [4] T. Ichikawa et al.,phys. Rev. Lett. 9, 225 (22). [5] T. Ichikawa et al., Phys. Rev. C 89, 5(R) (24); ibid., 9, 44 (24). [6] A. Staszczak and C. Y. Wong, Phys. Lett. B78, 4 (24). [7] A. Staszczak and C. Y. Wong, Acta Phys. Pol. B 46, 675 (25). [8] A. Staszczak and C. Y. Wong, Phys. Scr. 9, 46 (25). [9] M. Brack, J. Damgaard, A. S. Jensen, H.C. Pauli, V. M. Strutinsky, C. Y. Wong, Rev. Mod. Phys (972). [2] P. W. Zhao, Z. P. Li, J. M. Yao and J. Meng, Phys. Rev. C 82, 549 (2). [2] A. Bohr and B. Mottelson, Nucl. Phys. A45, c (98). [22] M. Ploszajczak et al., Nucl. Phys. A, 477(978). [2] M. B. Tsang et al., Eur. Phys. J. A, 29 (26). [24] J. Plateau, Annual Report of the Board of Regents of the Smithsonian Institution, 27, (86); ibid., 286 (865); ibid., 255 (866). [25] Lord Rayleigh, Phil. Mag. 28, 6 (94). [26] E. Pairam and A. Fernández-Nieves, Phys. Rev. Lett. 2, 245 (29). [27] A. K. Nurse, S. R. Coriell and G. B. McFadden, J. Res. NIST 2, 74 (25). [28] M. A. Fontelos, V. J. Garcia-Garrido, and U. Kindelán, SIAM J. Appl. Math. 7 (6), 94 (2). [29] H. Hofmann, The Physics of Warm Nuclei: With Analogies to Mesoscopic Systems, Oxford University Press, USA (28). [] Clusters in Nuclei, ed. C. Beck, Springer International Publishing, (2); ibid., 2 (22); ibid., (24). [] K. Ikeda, N. Tagikawa and H. Horiuchi, Prog. Theor. Phys. Suppl. E 68, 464 (968). [2] W. von Oertzen, M. Freer and Y. Kanada-En yo, Phys. Rep. 42, 4 (26). [] M. Freer, Rep. on Prog. in Phys (27). [4] Y. Kanada-En yo, H. Horiuchi, Prog. Theor. Phys. Suppl.

11 42, 25 (2). [5] Y.Funaki, T. Yamada, H.Horiuchi, G. Röpke, P. Schuck, and A. Tohsaki, Phys. Rev. Lett (28). [6] D. H. Wilkinson, Nucl. Phys. A452, 296 (986). [7] W. von Oertzen and M. Milin, Covalent Binding on the Femtometer Scale: Nuclear Molecules, Clusters in Nuclei - Vol. - Lecture Notes in Physics Springer Heidelberg, 875, 47 (24). [8] K. Hagel et al., Phys. Rev. C 5, 27 (994). [9] Y. Larochelle et al., Phys. Rev. C 5, 82 (996). [4] A. Ono, Phys. Rev. C 59, 85 (999). [4] W. Cassing, Z. Phys. A27, 87 (987). [42] R. J. Charity, Phys. Rev. C 82, 46 (2). [4] D. Mancusi, R. J. Charity and J. Cugnon, Phys. Rev. C 82, 446 (2). [44] D. Lacroix, A. Van Lauwe, and D. Durand, Phys. Rev. C 69, 5464 (24). [45] W. W. Wilcke et al., Atom. Dat. Nucl. Dat. Tab. 25, 89 (98). [46] S. Wuenschel et al., Nucl. Instr. Meth. Phys. Res. A64, 578 (29). [47] R. Wada et al., Phys. Rev. C 69, 446(24). [48] W. Lin et al., Phys Rev C 9, 446 (24). [49] S. I. Bityukov and N. V. Krasnikov, arxiv:physics 9897v2 (998). [5] G. Fai and J. Randrup, Nucl. Phys. A44, 55 (98). [5] J. P. Bondorf, C. H. Dasso, R. Donangelo and G. Pollarolo, Phys. Lett. B24, 28 (99). [52] J. Meng, ed., Relativistic Density Functional for Nuclear Structure, International Review of Nuclear Physics, Vol. (World Scientfic, Singapore, 26). [5] J. Meng, J. Peng, S. Q. Zhang, P. W. Zhao, Front. Phys. 8, 55 (2). [54] J. Meng, H. Toki, S.-G. Zhou, S. Zhang, W. Long, and L. Geng, Progress in Particle and Nuclear Physics 57, 47 (26). [55] H. Z. Liang, J. Meng, and S.-G. Zhou, Physics Reports 57, (25). [56] G. A. Lalazissis, T. Ni sć and D. Vretenar, P. Ring, Phys. Rev. C 7, 242 (25). [57] Z. X. Ren, S. Q. Zhang, J. Meng, Phys. Rev. C 95, 24 (27). [58] Z. X. Ren et al., to be published. [59] K. Schmidt et al., Phys. Rev. C 95, 5468 (27). [6] B. Schuetrumpf and W. Nazerewicz, Phys. Rev. C 96, 6468 (27).