Symmetry Energy of Dilute Warm Nuclear Matter

  • J. B. Natowitz, G. Röpke, S. Typel, D. Blaschke, A. Bonasera, K. Hagel, T. Klähn, S. Kowalski, L. Qin, S. Shlomo, R. Wada, and H. H. Wolter, Cyclotron Institute, Texas A&M University, College Station, Texas 77843-3366, USA, Institut für Physik, Universität Rostock, Universitätsplatz 3, D-18055 Rostock, Germany, Excellence Cluster Universe, Technische Universität München, Boltzmannstraße 2, D-85748 Garching, Germany, GSI Helmholtzzentrum für Schwerionenforschung GmbH, Theorie, Planckstraße 1, D-64291 Darmstadt, Germany, Instytut Fizyki Teoretycznej, Uniwersytet Wrocławski, plac Maksa Borna 9, 50-204 Wrocław, Poland, Bogoliubov Laboratory for Theoretical Physics, JINR Dubna, Joliot-Curie street 6, 141980 Dubna, Russia, Laboratori Nazionali del Sud-INFN, via Santa Sofia 64, 95123 Catania, Italy, Theory Group, Physics Division, Building 203, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, USA, and Fakultät für Physik, Universität München, Am Coulombwall 1, D-85748 Garching, Germany
Phys. Rev. Lett. 104, 202501
Free (a) and internal (b) symmetry energy as a function of the surface velocity. Experimental results are compared with results of theoretical calculations neglecting cluster formation (RMF) and including cluster formation (QS).

The symmetry energy [1] in the nuclear equation of state governs phenomena from the structure of exotic nuclei to astrophysical processes. The structure and the composition of neutron stars depend crucially on the density dependence of the symmetry energy [2]. As a general representation of the symmetry energy coefficient we use the definition

where E(n,δ,T) is the energy per nucleon of nuclear matter with density n, temperature T, and asymmetry δ=(N-Z)/A with Z and N the proton and neutron numbers, and A=N+Z. At low density the symmetry energy changes mainly because additional binding is gained in symmetric matter due to formation of clusters and pasta structures [3].

Our empirical knowledge of the symmetry energy near the saturation density n0 is based primarily on the binding energies of nuclei. The Bethe-Weizsäcker mass formula leads to values of about Esym(n0,0)=2834MeV for the symmetry energy at zero temperature and saturation density n00.16fm-3, if surface asymmetry effects are properly taken into account [4].

In contrast with the value of Esym(n0,0), the variation of the symmetry energy with density and temperature is intensely debated. Many theoretical investigations have been performed to estimate the behavior of the symmetry energy as a function of n and T. A recent review is given by Li et al. [5]; see also [6,7]. Typically, quasiparticle approaches such as the Skyrme Hartree-Fock and relativistic mean-field (RMF) models or Dirac-Brueckner Hartree-Fock (DBHF) calculations are used. In such calculations the symmetry energy tends to zero in the low-density limit for uniform matter. However, in accordance with the mass action law, cluster formation dominates the structure of low-density symmetric matter at low temperatures. Therefore, the symmetry energy in this low-temperature limit has to be equal to the binding energy per nucleon associated with the strong interaction of the most bound nuclear cluster. A single-nucleon quasiparticle approach cannot account for such structures. The correct low-density limit can be recovered only if the formation of clusters is properly taken into account, as has previously been shown in Ref. [8] in the context of a virial expansion valid at very low densities and in Ref. [9].

In this Letter we employ a quantum-statistical (QS) approach which includes cluster correlations in the medium. It interpolates between the exact low-density limit and the very successful RMF approaches near the saturation density. We show that this picture is in agreement with recent experimental findings on Esym at very low densities.

Suitable approaches to account for cluster formation are the nuclear statistical equilibrium (NSE) model [10], cluster-virial expansions [8], and generalized Beth-Uhlenbeck approaches [11]. A thermal Green function approach that allows a generalization of the NSE model by introducing a quasiparticle description also for the bound states was already formulated some decades ago [12], but only recently analyzed with respect to the consequences for nuclear matter [13]. In this QS approach the cluster correlations are described in a generalized Beth-Uhlenbeck expansion. The advantage of this method is that the medium modifications of the clusters at finite density are taken into account. In Ref. [9] the thermodynamic properties of nuclear matter were derived using this approach. The formulation of Ref. [9] is valid in the density and temperature range where the formation of light clusters with A4 dominates and heavier clusters are not yet important. The method requires a sufficiently accurate model for the quasiparticle properties, for which we employ a RMF model with density dependent couplings [14] which gives a good description both of nuclear matter around normal density and of ground state properties of nuclei across the nuclear chart. In order to extend the applicability of this RMF model to very low densities, it has been generalized in Ref. [9] to account also for cluster formation and dissolution.

We note that at very low densities and temperatures below T1MeV new phases may occur. In fact, the formation of a solid phase using Overhauser orbitals including a triple point [15] or Bose-Einstein condensation [16] has been suggested. However, in this Letter we are concerned with experimental data which probe nuclear matter at considerably higher temperatures. The low- T behavior is an interesting issue for future studies.

In the following we focus on finite temperatures and on the subsaturation region n<n0. Experimental information is derived from heavy-ion collisions of charge asymmetric nuclei, where transient states of different density can be reached, depending on the incident energy and the centrality of the collision. In the Fermi energy domain symmetry energy effects have been investigated using judiciously chosen observables [5,17–20].

Recently, the experimental determination of the symmetry energy at very low densities produced in heavy-ion collisions of Zn64 on Mo92 and Au197 at 35 MeV per nucleon has been reported [21]. Results of this study are given in the first four columns of Table I. Note that, as a result of an energy recalibration and reevaluation of the particle yields in different velocity bins, these values are slightly different than those reported in Ref. [21]. The surface velocity vsurf, i.e., the velocity before the final Coulomb acceleration, was used as a measure of the time when the particles leave the source under different conditions of density and temperature. Only values of vsurf<4.5cm/ns are included here, since the system does not reach equilibrium for higher vsurf; see Table I of Ref. [21]. The yields of the light clusters A4 were determined as a function of vsurf. Temperatures were determined with the Albergo method [22] using a H-He thermometer based on the double yield ratio of deuterons, tritons, He3, and He4, and are given in Table I as the average for the two reactions.

Comparisons of the scaled internal symmetry energy Esym(n)/Esym(n0) as a function of the scaled total density n/n0 for different approaches and the experiment. (a) The symmetry energies for the commonly used MDI parametrization of Chen et al. [24] for T=0 and different asy stiffnesses, controlled by the parameter x [dotted, dot-dashed, and dashed (black) lines], for the QS model including light clusters for temperature T=1MeV [solid (green) line], and for the RMF model at T=0 including heavy clusters [long-dashed (orange) line]. (b) The internal scaled symmetry energy in an expanded low-density region. Shown are again the MDI curves and the QS results for T=1, 4, and 8 MeV compared with the experimental data with the NSE entropy (solid circles) and the results of the self-consistent calculation (open circles) from Table .Comparisons of the scaled internal symmetry energy Esym(n)/Esym(n0) as a function of the scaled total density n/n0 for different approaches and the experiment. (a) The symmetry energies for the commonly used MDI parametrization of Chen et al. [24] fo... Show more

The free neutron yield is obtained from the free proton yield and the yield ratio of H3 to He3. To determine the asymmetry parameter of the sources the total proton and neutron yields including those bound in clusters are used. The proton chemical potential is derived from the yield ratio of H3 to He4. The corresponding free proton and free neutron densities are calculated, and the total nucleon density is obtained by accounting also for the bound nucleons according to their respective yields [21]. The total nucleon densities are of the order of 1/100th to 1/20th of saturation density, as seen in Table I.

An isoscaling analysis [23] has been employed (as a function of vsurf) to determine the free symmetry energy Fsym via the expression α=4FsymΔ(Z/A)2/T. Here α is the isoscaling coefficient determined from yield ratios of Z=1 ejectiles of the two reactions and Δ(Z/A)2 is the difference of the squared asymmetries of the sources in the two reactions. With Δ(Z/A)2 and the temperature determined as above, the free symmetry energy is extracted.

Using the free symmetry energy derived in this way from the measured yields, the internal symmetry energy can be calculated if the symmetry entropy is known. The values of the symmetry entropy SsymNSE for given parameters of temperature and density within the NSE model are shown in Table I, column 5. They are calculated with the equivalent expression of Eq. (1) as the difference between the entropies of pure proton or neutron and symmetric nuclear matter. In contrast with the mixing entropy that leads to a larger entropy for uncorrelated symmetric matter in comparison with pure neutron matter, the formation of correlations, in particular, clusters, will reduce the entropy in symmetric matter; see also Fig. 9 of Ref. [9]. For parameter values for which the yields of free nucleons in symmetric matter are small, the symmetry entropy may become positive, as seen in Table I for low temperatures. The fraction of nucleons bound in clusters can decrease, e.g., due to increasing temperature or the dissolution of bound states at high densities due to the Pauli blocking. Then, the symmetric matter recovers its larger entropy so that the symmetry entropy becomes negative, as seen in Table I also in the QS and self-consistent (SC, see below) calculations.

The results obtained in this way for the internal symmetry energy Esym=Fsym+TSsymNSE are shown in Table I, column  6. We note that in Ref. [21] the symmetry entropy was estimated using results of the virial expansion of Ref. [8] leading to different internal symmetry energies. However, this approximation is unreliable at the densities considered here.

In Table I, we also give results of the QS model [9] for the free and internal symmetry energies (columns 7 and 8) at given T and n. In Fig. 1(a) the experimentally obtained free symmetry energy is compared to the results of the RMF calculation without clusters and the QS model with clusters [9]. There are large discrepancies between the measured values and the results of calculations in the mean-field approximation when cluster formation is neglected. On the other hand, the QS model results correspond nicely to the experimental data. In Fig. 1(b) we compare the internal symmetry energy derived from the experimental data with the RMF and QS results. Again, it is clearly seen that the quasiparticle mean-field approach (RMF without clusters) disagrees strongly with the experimentally deduced symmetry energy while the QS approach gives a rather good agreement with the experimental data.

In Fig. 2 we present results for different approaches of extracting the internal symmetry energy and compare with the experimental values. In Fig. 2(a) we show theoretical results for T at or close to zero. A widely used momentum-dependent parametrization of the symmetry energy (MDI) at temperature T=0MeV was given in Refs. [5,24] and is shown for different assumed values of the stiffness parameter x. For these parametrizations the symmetry energy vanishes in the low-density limit. We compare this to the QS result at T=1MeV (at lower T crystallization or Bose condensation may occur as discussed above). In this approach the symmetry energy is finite at low density. The T=1MeV curve will also approach zero at extremely low densities of the order of 10-5fm-3 because the temperature is finite. The RMF, T=0 curve is discussed below. Also note that the underlying RMF model for the quasiparticle description with n0=0.149fm-3, Esym(n0)=32.73MeV gives a reasonable behavior at high density similar to the MDI, x=0 parametrization. We thus see that our approach successfully interpolates between the clustering phenomena at low density and a realistic description around normal density.

In Fig. 2(b) we compare the theoretical results to the experimental ones, full (red) circles (Table I, column 6), in an expanded low-density region. In addition to the MDI parametrization, we show the QS results [9] for T=1, 4, and 8 MeV, which are in the range of the temperatures in the experiment. The QS results including cluster formation agree well with the experimental data points, as seen in detail in Fig. 1. We conclude that medium-dependent cluster formation has to be considered in theoretical models to obtain the low-density dependence of the symmetry energy that is observed in experiments.

The temperatures and densities of columns 2 and 3 in Table I will be modified if medium effects on the light clusters are taken into account [25]. We have carried out a self-consistent determination of the temperatures TSC and densities nSC taking into account the medium-dependent quasiparticle energies as specified in Ref. [13] (columns 9 and 10 of Table I). Compared to the Albergo method results [21], the temperatures TSC are about 10% lower. Significantly higher values are obtained for the inferred densities nSC which are more sensitive to the inclusion of medium effects. We have also calculated the free and internal symmetry energies corresponding to these self-consistent values of TSC and nSC according to Ref. [9] (columns 11 and 12 of Table I). These results are also shown in Fig. 2(b) as open (purple) circles. The resultant internal symmetry energies are 15%–20% higher than the QS model values for T and n given in columns 2 and 3 in Table I.

We have restricted our present work to that region of the phase diagram where heavier clusters with A>4 are not relevant. The generalization of the given approach to account for clusters of arbitrary size would lead to an improvement in the low-density, low-temperature region when nuclear statistical equilibrium is assumed. Alternatively, one can introduce the formation of heavier nuclei in the presence of a nucleon and cluster gas, cf. Refs. [26,27].

The simplest approach to model the formation of heavy clusters is to perform inhomogeneous mean-field calculations in the Thomas-Fermi approximation assuming spherical Wigner-Seitz cells. In Fig. 2(a) preliminary results for the zero-temperature symmetry energy of such a calculation is shown by the long-dashed line using the same RMF parametrization as for the QS approach introduced above; for details see Ref. [28]. The symmetry energy in this model approaches a finite value at zero density in contrast with the behavior of the MDI parametrizations and conventional single-nucleon quasiparticle descriptions.

In conclusion, we have shown that a quantum-statistical model of nuclear matter, which includes the formation of clusters at densities below nuclear saturation, describes quite well the low-density symmetry energy which was extracted from the analysis of heavy-ion collisions. Within such a theoretical approach the composition and the thermodynamic quantities of nuclear matter can be modeled in a large region of densities, temperatures, and asymmetries that are required, e.g., in supernova simulations.


  1. We use the historical term “symmetry energy” although the community should invoke renaming it “asymmetry energy” as it characterizes the dependence of the nuclear binding energy on the asymmetry δ=(NZ)/A.

  2. J. M. Lattimer and M. Prakash, Phys. Rep. 442, 109 (2007).
  3. G. Watanabe et al., Phys. Rev. Lett. 103, 121101 (2009).
  4. P. Danielewicz, Nucl. Phys. A727, 233 (2003).
  5. B. A. Li et al., Phys. Rep. 464, 113 (2008).
  6. C. Fuchs and H. H. Wolter, Eur. Phys. J. A 30, 5 (2006).
  7. T. Klähn et al., Phys. Rev. C 74, 035802 (2006).
  8. C. J. Horowitz and A. Schwenk, Nucl. Phys. A776, 55 (2006).
  9. S. Typel, G. Röpke, T. Klähn, D. Blaschke, and H. H. Wolter, Phys. Rev. C 81, 015803 (2010).
  10. J. P. Bondorf et al., Phys. Rep. 257, 133 (1995).
  11. M. Schmidt, G. Röpke, and H. Schulz, Ann. Phys. (N.Y.) 202, 57 (1990).
  12. G. Röpke et al., Nucl. Phys. A379, 536 (1982); ; A424, 594 (1984).
  13. G. Röpke, Phys. Rev. C 79, 014002 (2009).
  14. S. Typel, Phys. Rev. C 71, 064301 (2005).
  15. H. R. Jaqaman, Phys. Rev. C 38, 1418 (1988).
  16. H. Takemoto et al., Phys. Rev. C 69, 035802 (2004).
  17. V. Baran et al., Phys. Rep. 410, 335 (2005).
  18. M. B. Tsang et al., Phys. Rev. Lett. 102, 122701 (2009).
  19. S. Wuenschel et al., Phys. Rev. C 79, 061602(R) (2009).
  20. C. Sfienti et al., Phys. Rev. Lett. 102, 152701 (2009).
  21. S. Kowalski et al., Phys. Rev. C 75, 014601 (2007).
  22. S. Albergo et al., Nuovo Cimento Soc. Ital. Fis. 89A, 1 (1985).
  23. M. B. Tsang et al., Phys. Rev. Lett. 86, 5023 (2001).
  24. L. W. Chen et al., Phys. Rev. Lett. 94, 032701 (2005); ; Phys. Rev. C 76, 054316 (2007).
  25. S. Shlomo et al., Phys. Rev. C 79, 034604 (2009).
  26. J. M. Lattimer et al., Nucl. Phys. A535, 331 (1991).
  27. H. Shen et al., Nucl. Phys. A637, 435 (1998).
  28. S. Typel et al. (to be published).

About the Authors

Image of J. B. Natowitz
Image of G. Röpke
Image of S. Typel
Image of D. Blaschke
Image of A. Bonasera
Image of K. Hagel
Image of T. Klähn
Image of S. Kowalski
Image of L. Qin
Image of S. Shlomo
Image of R. Wada
Image of H. H. Wolter

Related Articles

Focus: Better Microwaves from a Metamaterial

Focus: Better Microwaves from a Metamaterial

A system that produces radiation by sending electrons through a metamaterial waveguide could potentially generate high-power microwaves with high efficiency. Read More »

Arts &amp; Culture: Science and Filmmaking in an Equal Relationship

Arts & Culture: Science and Filmmaking in an Equal Relationship

Imagine Science—a yearly film festival that takes place in New York, Paris, and Abu Dhabi—is forging new ways to bring together scientists and filmmakers. Read More »

Synopsis: Explaining Aftershock Clustering

Synopsis: Explaining Aftershock Clustering

A study of bursting phenomena like earthquakes suggests that events appear to cluster in time because of the way that small events like aftershocks are identified. Read More »

More Articles