Primary Isotope Yields and Characteristic Properties of the Fragmenting Source in Heavyion Reactions near the Fermi Energies
Abstract
For central collisions of Ca Ca at 35 MeV/nucleon, the density and temperature of a fragmenting source have been evaluated in a selfconsistent manner using the ratio of the symmetry energy coefficient relative to the temperature, , extracted from the yields of primary isotopes produced in antisymmetrized molecular dynamics (AMD) simulations. The values are extracted from all isotope yields using an improved method based on the Modified Fisher Model (MFM). The values of obtained, using different interactions with different density dependencies of the symmetry energy term, are correlated with the values of the symmetry energies at the density of fragment formation. Using this correlation, the fragment formation density is found to be . Using the input symmetry energy value for each interaction temperature values are extracted as a function of isotope mass . The extracted temperature values are compared with those evaluated from the fluctuation thermometer with a radial flow correction.
pacs:
25.70.PqGBK
I. Introduction
A wealth of nuclear phenomena, including dynamical processes in nuclear reactions, nuclear structure and nuclear astrophysics are intimately affected by the density dependence of the symmetry energy term in the nuclear equation of state (EOS) Baoan2008 (). Investigations of the density dependence of the symmetry energy have been conducted using many observables such as isotopic yield ratios Tsang2001 (), isospin diffusion Tsang2004 (), neutronproton emission ratios Famiano2006 (), giant monopole resonances Li2007 (), pygmy dipole resonances Klimkiewicz2007 (), giant dipole resonances Trippa2008 (), collective flow zak2012 () and isoscaling Xu2000 (); Tsang2001_1 (); Huang2010 (). The ratio of the symmetry energy coefficient relative to the temperature, , extracted from isotope yield ratios Huang2010_1 (); Huang2010_2 (), has been studied within the frame work of the Modified Fisher Model (MFM) Fisher1967 (); Minich1982 (); Hirsch1984 (); Bonasera2008 ().
While temperature is one of the key variables in characterizing nuclear reactions, it is very difficult to determine the temperature of hot nuclear matter in a dynamical process. Several nuclear thermometers have been proposed. These include the slope of energy spectra Westfall1982 (); Jacak1983 (), momentum fluctuations Wuenschel2010 (), double isotope yield ratios Albergo1985 () and excited state distributions Morrissey1984 () among others. However they may not be generally applicable in all circumstances and even for a given system the extracted temperature values from these thermometers may be quite different from each other Guo2013 (). In our recent work, the isotopic yield ratio method was applied to extract values from the experimentally reconstructed primary fragment yields. These ratios were compared to those calculated from AMD primary generated fragment yields obtained using Gogny interactions with different density dependencies of the symmetry energy Lin2014 (). In the analysis, we found that the extracted values change according to the interactions used. From the dependence on interaction, the density, symmetry energy and temperature at the time of fragment formation were determined in a selfconsistent manner.
In order to perform a more detailed investigation of the relationship between the experimentally extracted observables and input parameters such as the density dependent symmetry energy term in the model, we apply a similar procedure to AMD primary events, using AMD as an event generator. There are three major reasons to use AMD as the event generator for this work. One is its capability to reproduce the experimental isotope yields. AMD results, such as multiplicity, angular distribution and energy spectra, have often been compared with those from the experimental data for intermediate energy heavy ion collisions and reproduce them reasonably well Ono96 (); Ono99 (); Ono02 (); Ono04 (); Wada98 (); Wada00 (); Wada04 (); Hudan06 (). In one of our recent works in Ref.Lin2014 (), the yields of the experimentally reconstructed primary hot isotopes are well reproduced by those of the AMD simulations. Second is to eliminate the secondary cooling effect. As shown in Ref.Huang2010_1 (); Lin2014 (), the sequential decay of the primary hot isotopes significantly alters the yield distribution and distorts the information inherent in the primary hot fragment yields. Third is to simplify the initial conditions, using zero impact parameter to eliminate effects of transverse flow and neck emission among others. The AMD events are generated for central collisions ( fm) of at 35 MeV/nucleon, using interactions having different density dependencies of the symmetry energy term, i.e., the standard Gogny interaction which has an asymptotic soft symmetry energy (g0), the Gogny interaction with an asymptotic stiff symmetry energy (g0AS) and the Gogny interaction with an asymptotic superstiff symmetry energy (g0ASS) Ono99 (); Ono2003 (). The extracted temperature values are compared to those obtained using a fluctuation thermometer.
II. Improved MFM Model and Extraction of
In Ref. Huang2010_1 () only yields of isotopes with and 3 were used to extract the ratio of the symmetry energy coefficient relative to the temperature, . In this article, an improved method is used to extract the values. In the improved method, all available isotope yields are employed. The improvement is made possible by taking into account the mass dependent temperature in the free energy in an iterative manner.
In the framework of MFM, the yield of an isotope with mass and ( neutrons and protons) produced in a multifragmentation reaction, can be given as Minich1982 (); Hirsch1984 (); Bonasera2000 (); Bonasera2008 (); Huang2010 (); Huang2010_1 ()
(1) 
Using the generalized WeiszckerBethe semiclassical mass formula Weizsacker1935 (); Bethe1936 (), can be approximated as
(2) 
In Eq.(1), and originate from the increases of the entropy and the mixing entropy at the time of the fragment formation, respectively. () is the neutron (proton) chemical potential. is the critical exponent. In this work, the value of is adopted from the previous studies Bonasera2008 (). In general coefficients, , , , and the chemical potentials are temperature and density dependent. In this formulation a constant volume process at an equilibrium, is assumed in the free energy, and therefore the term ”symmetry energy” is used throughout this work along Ref. Marini2012 (). If one assumes a constant pressure at the equilibrium process Sobotka2011 (), the term ”symmetry enthalpy” should be used. Experimentally, whether the equilibrium process takes place at constant pressure or volume can not be determined, and thus we use ”symmetry energy” throughout the paper, keeping in mind the ambiguity Marini2012 ().
In the previous analyses Huang2010_1 (); Huang2010_2 (); Chen2010 (); Lin2014 (); Lin2014_2 (), the temperature in Eq.(1) was assumed to be identical to the temperature of the fragmenting source and treated as a constant for all isotopes. However as seen below, this temperature turns out to be (fragment) mass dependent. This mass dependence on the temperature was not recognized in the previous analyses, just because the mass dependence was masked by the larger error bars. In this improved method, the error bars become small and the mass dependence becomes evident. In order to take into account the mass dependence of the temperature in Eq.(1), the temperature is replaced by an apparent temperature . We attribute this mass dependence to the system size effect as discussed in Sec. IV.1. In the improved MFM formulation, therefore, this system size effect is empirically realized by reducing the apparent temperature as increases as . is the temperature of the fragmenting source and is a constant quantifying the mass dependence.
In order to study the density, temperature and symmetry energy in the fragmenting source, the improved MFM of Eq.(1) is utilized to extract the value from the available isotope yields. Since the value in Eqs.(1) and (2) depends on 5 parameters, , , , and (defined by ), the optimization process of these parameters is divided into the following three steps to minimize the ambiguity for each parameter. For a given value,

Optimize and values from mirror isobars and fix these parameter values.

Optimize , and values from isotopes.

Using extracted parameters in step (1) and step (2), values are extracted from all available isotopes. Comparing the extracted values with three different interactions, the density of the fragmenting source is extracted. Using this density, the value of the symmetry energy coefficient, , for each interaction is determined. The temperature is then calculated as the ratio of to .
It is expectable that if the value is properly selected which means the mass dependence is well considered, a constant is obtained. Since the value is small as seen below, we perform the optimization of the parameter in an iterative manner in the following analysis, that is, in the first round is set in and calculate the temperature as a function of , using steps (1)(3). From this plot a new value, , is extracted from the slope. In the second round, is used for the steps (1)(3) and a new value, , is extracted. If is 0 within a given error range, the iteration stops and the value is fixed as the mass dependent parameter of the apparent temperature and value is determined. Otherwise the iteration continues.
Details of steps (1)(3) are first described below for a given value. In the step (1), following Ref. Huang2010_1 (), the isotope yield ratio between isobars with and , , is utilized, which is
(3)  
where is the yield of isotopes with and , and . When the above equation is applied for a pair of mirror nuclei of odd mass isotopes with and , the symmetry energy term, pairing term and mixing entropy terms drop out and the following equation is obtained.
(4) 
For all available mirror isobars, and are optimized in Eq.(4). The values and the fit result for is shown in Fig.1.
In the step (2) we apply Eq.(1) to the isotopes with with the extracted and values in the step (1). For isotopes, the ratio of the free energy relative to the temperature can be calculated from Eq.(1) and Eq.(2) without the symmetry energy term as
(5) 
where . The value of on the right of the second equation can be calculated from the simulated or experimental values when the value is fixed. Non zero values show the deviation of the mass distribution of isotopes from the power law distribution determined by the critical exponent Bonasera2008 (). When other values are used, the parameter values change accordingly. In order to eliminate the constant , all isotope yields are normalized to the yield of Bonasera2008 (); Huang2010 (); Huang2010_1 (). For the first round with , the renormalized values of from the AMD events with the g0 interaction are plotted as a function of the isotope mass using solid points in Fig. 2(a). The values of , and are used as free parameters to fit the given values, employing Eq.(5). A typical search result is shown by open circles in Fig. 2(a) for the case of the g0 interaction. Similar quality results are obtained for the events generated using the g0AS and g0ASS interactions. One should note that the value of make a small contribution and the contribution is evident as a staggering in the vs plot. Therefore the essential free parameters in this step are and . The extracted parameter values are summarized in Table 1 for the first round and the final round .
k=0.0  

g0  1.77  2.74  
g0AS  1.76  2.66  
g0ASS  1.77  2.72  
k=0.007  
g0  1.55  2.29  
g0AS  1.53  2.20  
g0ASS  1.55  2.29 
In step (3) Eq.(1) is applied to yields of all isotopes with and . From Eq.(1), and values can be related to the modified free energy, as
(6) 
where is the free energy relative to the temperature, , subtracted by the calculated contributions of the volume, surface, Coulomb and pairing terms, using the parameters in Table 1. Resultant values are shown by symbols in Fig. 2(b). They exhibit quadratic relationships with minimum values close to zero. The minimum values are at or near isotopes and therefore reflect approximately the difference between the data and fits point in Fig. 2(a). In this step, the and the values are optimized. Since the values are extracted from the step(1), the optimization is made around the values in the fifth column of Table 1 in a small margin. The values are extracted from the quadratic curvature of the isotope distribution for each given and plotted in Fig. 2(c) separately for the g0, g0AS and g0ASS interactions. As one can see for the first round with , the extracted values increase as increases in all cases, and they more or less parallel each other.
III. Selfconsistent Determination of Density and Temperature
In order to determine the density and temperature at the time of the fragment formation, the parallel behavior of the observed values in Fig. 2(c) is utilized. As suggested in Ref. Ono2003 (), the observed differences are attributed to the difference of the symmetry energy at the density at the time of the fragment formation. The ratios between g0/g0AS and g0/g0ASS of for the first round are shown in Fig. 3(a). The ratios show flat distributions as a function of for both cases. The extracted average ratio values are shown by lines in the figure and the values are given in the first column of Table 2. In Fig. 3(b) the symmetry energy coefficient is plotted as a function of the density for the three interactions used in the calculations and in Fig. 3(c) their ratios, and , are plotted. Using the ratio values determined from Fig. 3(a) and the density dependence of the values in Fig. 3(c), the implied densities of the fragmenting sources are indicated by the shaded vertical areas shown in Fig. 3(c). The extracted density values for each case are given in the second column of Table 2. Assuming that the nucleon density should be same for the three different interactions used, the nucleon density of the fragmenting source is determined from the overlap of the extracted values. This assumption is reasonable for the central collisions because the nucleon density is mainly determined by the stiffness of the EOS and not by the density dependence of the symmetry energy term. From the overlapped density area in Figs. 3(c), is extracted as the density at the time of the fragment formation. Using this density value, the corresponding symmetry energy values at that density are extracted for the three different interactions from Fig. 3(b). They are given in the third column of Table 2.
k  int  (MeV)  

k=0.0  g0  26.60.3  
g0/g0AS  1.140.02  0.710.05  
g0AS  23.71.3  
g0/g0ASS  1.410.03  0.660.02  
g0ASS  18.70.7  
k=0.007  g0  26.80.2  
g0/g0AS  1.140.02  0.720.05  
g0AS  24.01.3  
g0/g0ASS  1.400.03  0.660.02  
g0ASS  18.70.7 
Once the symmetry energy value is determined for a given interaction, the temperature, , can be calculated as / (. The extracted values are shown as a function of by open symbols for the first round in Fig. 4. The larger errors of , comparing to those in Fig. 2(c), originate from the errors of and extracted for each interaction which are shown in the third column of Table 2 and Fig. 2(c), respectively. The temperature values extracted from the three different interactions agree with each other very well and show a monotonic decrease as increases from MeV at to MeV at . From this slope the extracted temperature as function of , , is determined for the first round, assuming .
The iteration are repeated four times in this work. The same plots as Fig.2, but with the value for the fourth (final) round, , are shown in Fig.5 and the extracted parameters are also given in Table 1. A very similar quality of results to those of the first round with are obtained, even though the optimized parameter values are quite different between those of the first round and of the fourth round . The extracted values parallel each other and show a rather flat distribution as a function of for up to 15 in Fig. 5(c). As seen in Fig. 4, in which the extracted values are shown by closed symbols as a function of for the fourth round with , the extracted values are consistent with 5.5 MeV within the error bars. Since values show a flat distribution as a function of , the iteration is stopped at this round and MeV is taken as the temperature of the emitting source.
The extracted parameter values of , density and symmetry energy values for the fourth round are very similar to those of the first round as shown in Table 2. The values and errors of these parameters are essentially determined by the ratios and their errors of the values between different interactions as discussed in Fig.3(a). These ratio values are stable between the first and fourth rounds, even though the optimized parameters values in Table 1 are quite different between these rounds. The extracted value of the mass dependent factor has some errors, but the above fact ensures that the extracted density, temperature and symmetry values and their errors in Table 2 are rather stable independent of the choice of the values within its error bar, when the parameters in Table 1 are optimized for the given value.
The decreasing trend of the temperature as increases is often observed in heavy ion collisions and normally attributed to variations in impact parameter Trautmann07 (). The heavier fragments tend to be produced in more peripheral collisions and therefore show lower temperature. However in this study all events analyzed are generated in the same class of events, central collisions with fm. Therefore we attribute the decreasing trend may originate to the different fragment formation processes rather than to the centrality of the events.
IV. Discussion
IV.1 Temperature and Fluctuation thermometer
The temperature at the time of the fragment formation has been studied, using different thermometers Guo2013 () as mentioned earlier. Here we will investigate the applicability of the fluctuation thermometer using different isotopes, which gives the temperature values as a function of fragment mass A.
For a given variable , the fluctuation is given as . In our present work, focusing on momentum fluctuations, is used, following Ref. Wuenschel2010 (). Different relationships have been proposed to derive the temperature from observed momentum fluctuations. In Refs. Zheng2011 (); Zheng2012 (), Zheng et al. reported that there may be significant effects on the fluctuation temperature resulting from the quantum nature of the particles examined. When the temperature is determined from the fermions, such as n, p, t etc., the momentum distribution is expected to follow the FermiDirac distribution Zheng2011 (). On the other hand, when bosons, such as d or particles, are used, the BoseEinstein distribution is expected to be appropriate Zheng2012 (). In this paper the fluctuation temperature is calculated from IMF momentum distributions including those in the excited states from the primary fragments of AMD simulations without the afterburner. Even for a given IMF isotope, therefore, fermions and bosons are mixed and we assume the classical momentum distribution for these IMFs in the following analysis.
Under the assumption of the MaxwellBoltzmann distribution, the momentum distribution is,
(7) 
and the fluctuation of , , can be given as
(8) 
From above one can show that for a given particle with mass , and therefore the fluctuation temperature is obtained. The fluctuation temperature values calculated from the AMD events with the g0 interaction are represented by open circles in Fig. 6(b). Very similar results of the fluctuation temperature values are obtained for the other two interactions. As seen in the figure, over the entire range of isotope mass, the fluctuation temperature values are quite different from the values in the previous section. shows a broad peak MeV around , and then a monotonic decrease as increases. For comparison the fluctuation temperature values from the AMD events without the Coulomb interaction are also calculated and shown by open squares. They show a similar trend to those with the Coulomb interaction included, but about 1.0 to 1.5 MeV lower temperature in the entire mass range.
In order to understand the behavior of the extracted fluctuation temperature, simple MonteCarlo calculations are performed. In this model, a percolation model is utilized to generate fragments from 80 lattice points () Staufer1979 (). No charge and internal excitation energy are considered in the percolation model. The motion of particles was assigned as the sum of the thermal motion and radial expansion. One can use a more sophisticated method to generate the fragments, such as those used in statistical multifragmentation model (SMM) Bondorf85 () or microcanonical Metropolitan Monte Carlo model (MMMC) Gross90 (), but here we just use a simpler model available for our purpose, because we are only concerned with the distribution of thermal and collective motions assigned to these fragments. In the central AMD collisions ( fm), collective contributions to the particle energy occur only in the radial direction. Both radial flow and Coulomb repulsion contribute to the radial expansion motion. The Coulomb acceleration is treated as a part of the radial flow and the radial velocity is defined as , where is the flow parameter and is the coordinate vector. Fragments generated by the percolation model are distributed uniformly inside a sphere of radius under the condition of , where is the radius of the fragment with mass A. They are calculated as and , where is adopted from the above analysis. This condition allows some overlap of two spheres but for the purpose of the restriction of the distribution, this condition is good enough to see the global effect of the fluctuation temperature. Under a thermal equilibrium condition, the thermal velocity , where , is given by a MaxwellBoltzmann velocity distribution as
(9) 
where is the input parameter in the model. The total velocity is the sum of these two velocity vectors, . For a given parameter set of and values, which are determined below, the temperature is evaluated using the fluctuation thermometer and the extracted temperature values are shown in Fig. 6(a) for the thermal motion only, (solid curve), and for the sum of the thermal motion and radial motion, (dashed curve). As expected, the temperature values for the thermal motion show a constant temperature. The temperature with the radial flow does not increase linearly as increases because of the condition on the radial distance of the fragments . In order to apply this model to the reaction events, momentum conservation, for all fragments in one event), has to be considered. In the actual application, among more than a hundred million events used in Fig. 6(a), only those with MeV/c are selected as an approximation of . When the momentum conservation of the system is required, the momenta of the larger fragments tend to be limited in smaller side relative to those of used in Fig. 6(a). The results are shown by solid and dashed lines in Fig. 6(b). The parameters and are adjusted to reproduce the fluctuation temperature, extracted from the AMD events in Fig. 6(b). As seen in the figure, and show additional significant decreasing trends for heavier fragments, compared to those shown in Fig. 4(a). Using the parameter sets of MeV and c/fm ( c/fm without Coulomb), the values from both AMD calculations with and without Coulomb are reproduced reasonably well by . At the same time, , the radial flow corrected temperature, reproduces well the temperature values of from the selfconsistent symmetry energy analysis of the AMD events.
IV.2 Radial flow
From the above MonteCarlo simulation, the average radial flow energy, excluding the Coulomb contribution, is evaluated as MeV/nucleon from Fig. 6(b). The average Coulomb repulsive energy is evaluated as MeV/nucleon. Since the radial flow energy is not modified by the secondary decay cooling process, the extracted radial flow energy in the above analysis can be compared to the existing experimental data as shown in Fig. 7. Most of the existing data are distributed along the line and, as seen in the figure, the extracted radial flow energy is consistent with that trend.
IV.3 Volume and surface contribution of Symmetry energy
The increase of seen in Fig.2(c) can be attributed to the surface contribution of the symmetry energy in finite nuclei, assuming a constant temperature. In the WeiszckerBethe semiclassical mass formula used in Eq.(2), the coefficient, , is independent of the nuclear size, which originates from the volume nature of the symmetry energy. However advanced mass formulas have introduced size dependence on reflecting the surface effect Mayer1966 (); Danielewicz2003 (). Following the formulation of Danielewicz et al. Danielewicz2003 (); Danielewicz2009 (), the symmetry coefficient can be expressed as where and the suffixes V,S denote volume and surface contributions, respectively. One should note that the value is independent of the temperature of the emitting source. has been obtained from ground state nuclei Jiang2012 (). When we apply this formula to fit the results in Fig. 2(c) for the g0 interaction, one can get and . The fit result is shown in Fig.8. If we assume and use the value MeV from Table 2, we get MeV, which is much lower than that extracted above from the fluctuation thermometer. The extracted temperature independent value, is also significantly smaller than values extracted from the ground state nuclei. Since the mass dependence of the temperature extracted above is consistent to the fluctuation temperature after the radial flow correction, we conclude that the mass dependence of the values in Fig. 2(c) originates from the mass dependent apparent temperature and that the surface contribution of the symmetry energy is small, if any. This conclusion is also consistent with our recent results of Zn + Sn at 40 MeV/nucleon, a slightly larger system, when the improved MFM method is applied to the reconstructed isotope yields Liu2014 (). In the analysis a slightly smaller mass dependence of values is observed. If the mass dependence originates from the surface contribution, the mass dependence should remain same. On the other hand this smaller mass dependence is consistent with the MonteCarlo analysis in the previous section for the system size of .
V Conclusions and Summary
The following conclusions are drawn from the above analysis from central collision events of AMD for Ca + Ca at 35 MeV/nucleon:

The density of the fragmenting source is evaluated as , using the extracted values, based on the improved MFM model and the comparisons between the results from different density dependent symmetry energy terms.

The apparent temperature values show a monotonic decrease as the fragment decreases from MeV at to MeV at . This gives the mass dependent apparent temperature .

The fluctuation thermometer, uncorrected for flow effects, gives larger temperature values and a different shape compared to the apparent temperatures extracted from above. They show a structure, i.e., a broad peak at and a monotonic decrease for larger . Those from the AMD events without the Coulomb interaction show a similar trend but exhibit about MeV lower values.

After the correction for the radial flow requiring momentum conservation in a MonteCarlo model, the behavior of the fluctuation thermometer values are well understood as the combination of the thermal and radial collective motions. The bell shape structure of the temperature values is well reproduced by the model. The mass dependence of the apparent temperature observed from the selfconsistent symmetry energy analysis originates from the system size effect through momentum conservation. This indicates that the fluctuation thermometer can be a reasonable probe to measure the temperature of the fragmenting source in a multifragmentation process, if one can properly correct the contribution of the collective motion. The evaluated expansion energy from the fluctuation thermometer is consistent with existing experimental values.

Surface contribution of the symmetry energy to the values is examined. We conclude that the surface contribution in is small, if any.
Summarizing, an improved method is proposed for the extraction of the ratio of the symmetry energy coefficient relative to the temperature, , taking into account the mass dependence of the apparent temperature, based on the MFM model. This method is applied for the central collisions of the AMD events generated for at 35 MeV/nucleon. Gogny interactions, g0, g0AS and g0ASS, with three different density dependencies of the symmetry energy are employed. As a function of IMF charge , the ratios of the extracted values from the different interactions are essentially constant and reflect the differences of the symmetry energy at the density at the time of the fragment formation. Using this correlation, is evaluated as the density at the time of fragmenting source and the symmetry energy value at that density are extracted for each interaction. The temperature values are then determined. The extracted temperatures show a monotonic decrease as the fragment increases, changing from 5 MeV to 3 MeV when increases from 4 to 18. The extracted temperature values are compared to those of the fluctuation thermometer, assuming . The fluctuation thermometer shows significantly higher temperature values and a different trend as a function of fragment mass. After the correction for the radial flow the trends and the differences are well understood as resulting from the combination of the thermal and radial collective motions under momentum conservation. The mass dependence of the apparent temperature observed from the improved method originates from the system size effect under this momentum conservation. This indicates that the fluctuation thermometer can provide a reasonable probe for the temperature of the multifragmentation source if the contribution of collective flows is properly corrected.
Acknowledgments
X.Liu and W. Lin thank Prof. J. B. Natowitz for his hospitality during the stay in the Cyclotron Institute, TAMU. This work is supported by the National Natural Science Foundation of China (Grants No. 11075189) and 100 Persons Project (0910020BR0, Y010110BR0) and ADS project 302 (Y103010ADS) of the Chinese Academy of Sciences, the U.S. Department of Energy under Grant No. DEFG0393ER40773 and the Robert A. Welch Foundation under Grant A0330.
References
 (1) BaoAn. Li et al., Phys. Rep. 464, 113 (2008).
 (2) M. B. Tsang et al., Phys. Rev. Lett. 86, 5023 (2001).
 (3) M. B. Tsang et al., Phys. Rev. Lett. 92, 062701 (2004).
 (4) M. A. Famiano et al., Phys. Rev. Lett. 97, 052701 (2006).
 (5) T. Li et al., Phys. Rev. Lett. 99, 162503 (2007).
 (6) A. Klimkiewicz et al., Phys. Rev. C76, 051603 (2007).
 (7) L. Trippa et al., Phys. Rev. C77, 061304 (2008).
 (8) Z. Kohley et al., Phys. Rev. C85, 064605 (2012).
 (9) H. S. Xu et al., Phys. Rev. Lett. 85, 716 (2000).
 (10) M. B. Tsang et al., Phys. Rev. C64, 054615 (2001).
 (11) M. Huang et al., Nucl. Phys. A847, 233 (2010).
 (12) M. Huang et al., Phys. Rev. C81, 044620 (2010).
 (13) M. Huang et al., Phys. Rev. C82, 054602 (2010).
 (14) M. E. Fisher, Rep. Prog. Phys. 30, 615 (1967).
 (15) R. W. Minich et al., Phys. Lett. B118, 458 (1982).
 (16) A. S. Hirsch et al., Nucl. Phys. A418, 267c (1984).
 (17) A. Bonasera et al., Phys. Rev. Lett. 101, 122702 (2008).
 (18) G. D. Westfall et al., Phys. Lett. B116 118, (1982).
 (19) B. V. Jacak et al.. Phys. Rev. Lett. 51, 1846 (1983).
 (20) S. Wuenschel et al., Nucl Phys A843, 1, (2010).
 (21) S. Albergo et al. Nuovo Cimento A89, 1, (1985).
 (22) D. J. Morrissey et al., Phys. Lett. B148, 423 (1984).
 (23) C. Guo and F. Zhang, Nuclear Science and Techniques 24, 050513 (2013).
 (24) W. Lin et al., Pjys. Rev. C89, 021601(R) (2014).
 (25) A. Ono and H. Horiuchi, Phys. Rev. C53, 2958 (1996).
 (26) A. Ono, Phys. Rev. C59, 853 (1999).
 (27) A. Ono, S. Hudan, A. Chbihi, J. D. Frankland, Phys. Rev. C66, 014603 (2002).
 (28) A. Ono and H. Horiuchi, Prog. Part. Nucl. Phys. 53, 501 (2004).
 (29) R. Wada et al., Phys. Lett. B422, 6, (1998).
 (30) R. Wada et al., Phys. Rev. C62, 034601 (2000).
 (31) R. Wada et al., Phys. Rev. C69, 044610 (2004).
 (32) S. Hudan, R. T. de Souza, A. Ono, Phys. Rev. C73, 054602 (2006).
 (33) A. Ono, P. Danielewicz et al., Phys. Rev. C68, 051601(R) (2003). For g0AS, x=1/2 and for g0ASS x=2 are used in Eq.(2) in Ref. Ono2003 ().
 (34) A. Bonasera et al., Riv. Nuovo Cimento 23, 1 (2000).
 (35) C. F. von Weizscker, Z. Phys. 96, 431 (1935).
 (36) H. A. Bethe, Rev. mod. Phys. 8, 82 (1936).
 (37) P. Marini et al., Phys. Rev. C85, 034617 (2012).
 (38) L. G. Sobotka, Phys. Rev. C 84, 017601 (2011).
 (39) Z. Chen et al., Phys. Rev. C81, 064613 (2010).
 (40) W. Lin et al., to be published.
 (41) W. Trautmann et al., Phys. Rev. C76, 064606 (2007).
 (42) H. Zheng and A. Bonasera, Phys. Lett. B696, 178, (2011).
 (43) H. Zheng, G. Giuliani and A. Bonasera, Nucl. Phys. A892, 43, (2012).
 (44) D. Staufer, Phys. Rep. 54, 1 (1979), p=0.28 is used. The other choice of p deos not change the results.
 (45) J. P. Bondorf, R. Donangelo, I. N. Mishustin, Nucl. Phys. A443, 321 (1985).
 (46) D. H. E. Gross, Rep. Prog. Phys. 53, 605 (1990).
 (47) B. Borderie and M. F. Rivet, Prog. Part. Nucl. Phys. 61, 551 (2008) and see references therein.
 (48) W. D. Myers and W. J. Swiatecki, Nucl. Phys. A81, 1 (1966).
 (49) P. Danielewicz, Nucl. Phys. A727, 233 (2003).
 (50) P. Danielewicz and Jenny Lee, Nucl. Phys. A818, 26 (2009).
 (51) H. Jiang, G. J. Fu, Y. M. Zhao and A. Arima, Phys. Rev. C85, 024301 (2012) and references therein.
 (52) X. Liu et al., to be published.