Coexistence of phases in asymmetric nuclear matter under strong magnetic fields
Abstract
The equation of state of nuclear matter is strongly affected by the presence of a magnetic field. Here we study the equilibrium configuration of asymmetric nuclear matter for a wide range of densities, isospin composition, temperatures and magnetic fields. Special attention is paid to the low density and low temperature domain, where a thermodynamical instability exists. Neglecting fluctuations of the Coulomb force, a coexistence of phases is found under such conditions, even for extreme magnetic intensities. We describe the nuclear interaction by using the non–relativistic Skyrme potential model within a Hartree–Fock approach. We found that the coexistence of phases modifies the equilibrium configuration, masking most of the manifestations of the spin polarized matter. However, the compressibility and the magnetic susceptibility show clear signals of this fact. Thermal effects are significative for both quantities, mainly out of the coexistence region.
pacs:
21.65.Cd,26.60.c,97.60.Jd,21.30.xI Introduction
The dense nuclear matter under magnetic fields has been intensively studied (see LAI () and references therein), particularly in relation to astrophysical issues. Investigations of the neutron star structure LATTIMER () as well as the cooling of magnetized stars SHIBANOV (); YAKOVLEV (); BEZCHASTNOV () need the equation of state for magnetized matter as an important input. The presence of very intense magnetic fields in compact stellar objects has been proposed, based on the observational evidence of periodic or irregular radiation from localized sources. According to the energy released and the periodicity of the episodes, these objects have been classified as pulsars, soft gamma ray repeaters and anomalous Xray pulsars. They have been associated with different stages of the evolution of neutron stars. On the star surface the magnetic field could reach values G, as in the case of magnetars and it is expected a growth of several orders of magnitude in its dense interior.
Recent investigations KHARZEEV (); MO (); SKOKOV () have pointed out that matter created in heavy ion collisions could be subject to very strong magnetic fields. As a consequence the particle production would exhibit a distinguishable anisotropy. A preferential emission of charged particles along the direction of the magnetic field is predicted in KHARZEEV () for noncentral heavy ion collisions, due to magnetic intensities MeV. Improved calculations taking care of the mass distribution of the colliding ions MO (), does not modify essentially the magnitude of the produced fields. Furthermore, the numerical simulations performed by SKOKOV () predict larger values MeV.
The effects of magnetic fields on a dense nuclear environment have been described using different models CHAKRABARTY (); BRODERICK (); YUE2 (); SINHA (); RABHI (); RYU (); RYU2 (); ANG (); DONG (); DIENER (); DONG2 (); PGARCIA (); ISAYEV (); UNLP (); UNLP2 (); BORDBAR (); BIGDELI2 (). For instance, covariant field theoretical models have been used to study the role of the magnetic field on hyperonic matter YUE2 (); SINHA (), instabilities at subsaturation densities RABHI (); RYU (), magnetization of stellar matter DONG (), saturation properties of symmetric matter DIENER () and the symmetry energy DONG2 (). Nonrelativistic models have also been used, in the effective interactions of PGARCIA (); ISAYEV (); UNLP () or the microscopic models used in the variational calculations of BORDBAR (); BIGDELI2 (). A comparison of neutron matter results, using different models was presented in UNLP2 ().
It is a well known fact that the nuclear environment experiences thermodynamical instabilities at subnuclear densities and low temperatures. Evidence of this phenomenon can be found in the isospin distillation effect for heavy ion collisions DAS (). These instabilities give rise to a coexistence of phases if the surface tension is low enough. A more complex scenario is obtained when an external magnetic field is added, since there is a competition among two opposite trends. On one hand the magnetic force induces a globally ordered state with aligned spins. On the other hand the nuclear interaction favors the coexistence of phases where two states of different densities and spin polarizations are combined in order to lower the free energy.
In the present work we explore the possibility of a coexistence of phases for nuclear matter under strong magnetic fields, taking as variables the nuclear density, the isospin composition of matter, and the temperature MeV. The possible fields of applications, such as those mentioned before, show a complex scenario where the detailed physical mechanisms are not easy to understand because there is a superposition of effects which can combine to give very different manifestations. Therefore, we aim to present here some of the variables appearing in realistic situations, emphasizing the role of each of these factors, and to understand how they interact in a specific environment. Special attention is paid to the relevant quantities associated with them, as the spin polarization, the isothermal compressibility and the magnetic susceptibility. With this purpose in mind we have analyzed a wide range of isospin composition and we have also reached the extreme value G for the external magnetic field. This selection emphasizes the effects of these variables, which under certain conditions can appear weakened, or hidden by another factors.
The Skyrme model SKYRME (); DOUCHIN (); BENDER1 (); BENDER2 () is appropriate to describe the nuclear interaction under the conditions of interest. This is a nonrelativistic effective model where the in–medium nuclear force is simulated by a density dependent potential. It was successfully used to describe atomic nuclei as well as nuclear matter properties.
This article is organized as follows. We review the Skyrme model for nuclear matter under an external magnetic field in the next section, a brief resume of the Gibbs construction for the coexistence of phases is presented in Section III, the results are shown and discussed in Section IV. A final summary and the main conclusions are given in Section V.
Ii Spin polarized nuclear matter in the Skyrme model
The Skyrme model is an effective formulation of the nuclear interaction which has been employed profusely in the literature BENDER1 (). It consists of a basic Hamiltonian with contact nucleonnucleon potentials including density dependent coefficients,
where represent the Pauli matrices for spin,
is the spin exchange
operator, is the relative momentum
operator and is the total baryonic density.
Note that throughout this article we use units such that ,
.
The
interaction–parameters are fixed to cover a variety of
applications such as exotic nuclei or stellar matter. Using the
HartreeFock approximation, one can find an energy density
functional, which is a convenient way to study thermodynamical
properties of the system.
We are particularly interested in the contributions coming from
terms containing time reversalodd densities and currents, since
they are active when spin states are not symmetrically occupied. A
derivation of these terms can be found in BENDER2 (). We
assume the magnetic field has not dynamics, that is, it behaves as
an external field. There is a direct coupling between nucleons and
the magnetic field, due to their intrinsic magnetic moments. This
implies an additional term to the single
particle spectra of the standard Skyrme model, where is
the Bohr magneton and the Lande factors take account of
the anomalous magnetic moments. They take the values
and for protons and
neutrons, respectively.
Furthermore, the magnetic field induce a quantization of the
energy spectra of charged particles LANDAU (). In the case of
a uniform field, the corresponding Schrödinger equation
exhibits quantized eigenvalues, associated with the motion in
directions orthogonal to the applied field. They are
oscillatorlike levels, depending on a discrete quantum number in
the form , with the cyclotron
frequency of the particle. We can summarize the effects of a
uniform magnetic field over the spectra of nucleons by
(1)  
(2) 
The first two terms in the r.h.s of these equations, are the common results for the Skyrme model, which now have an implicit dependence on the field . The spin index () denotes a spin–up (spin–down) projection, the effective nucleon mass is defined by
(3) 
with the degenerate nucleon mass in vacuum, is
the isospin asymmetry fraction, with standing for
the particle number density of protons and neutrons respectively.
Note that . Since the spin states are not symmetrically
occupied, one can define for each isotopic component the number
density of particles with a given spin polarization .
The spin asymmetry density gives a measure of the spin
polarization , clearly . We have defined for
protons (neutrons).
In Eqs. (1) and (2) we have used the single
particle Skyrme potential energy
(4)  
The expressions for the kinetic energy density will be presented below. The Eqs. (1)–(4) have been written in terms of a set of density dependent coefficients and , which are related to the standard parameters of the Skyrme model by,
We assume baryonic and isospin number conservation, therefore independent chemical potential can be assigned to protons and neutrons. The corresponding distribution functions are the Fermi occupation number for a particle at temperature , with momentum and isospin and spin projections and , respectively.
Now we show explicit expressions for the density of particles with a given spin polarization, the kinetic energy and the isospin asymmetry densities, separately for protons and neutrons
(5)  
(6)  
(7)  
(8)  
(9)  
(10) 
For the proton related quantities, we have taken into account that, assuming along the zaxis, each eigenstate spreads over a bounded region of area in the plane. The component is not bounded and varies continuously. Therefore, the contribution of a charged particle to macroscopic quantities per unit volume has been evaluated by means of the replacement .
The energy density can be split into two terms,
(11) 
one of them depends on . The remaining one is similar to the common contribution of the Skyrme model in a HartreeFock approach, but now it depends implicitly on the magnetic intensity
For completeness we also show the expression for the entropy density in the quasiparticle approach,
The entropy is needed to build up the free energy and the pressure . The magnetization of the system is evaluated in terms of the grand canonical potential according to PATHRIA (),
(14) 
For the system considered, we have . As expected, it can be decomposed into proton and neutron contributions . Finally, the standard relations are used for the isothermal compressibility and magnetic susceptibility ,
Note that in our scheme we were able to develop analytical expressions for the isothermal compressibility and magnetic susceptibility. This gives us some confidence in the evaluation of these magnitudes, as for instance, the susceptibility for low temperatures has fast variations with the density.
In our approach both proton and neutron numbers are conserved separately, therefore the states of polarization of each component are also independent. The global polarization is determined by the condition of minimum free energy . This criterium differs from that in UNLP (), where the Legendre transformed potential was used.
The equilibrium state has a variable spin configuration, depending on and . As will be shown, in the low temperature and low density domain the coexistence of phases imposes a state with a lower degree of polarization than the case which does not consider the phase transition.
Iii Coexistence of phases in nuclear matter
The nuclear interaction gives rise to instabilities in a dense
infinite medium at low temperatures. If the Coulomb interaction is
taken into account and its fluctuations are included, a nucleation
process can be found.
Under the hypothesis assumed in the present work, the system
evolves through a succession of equilibrium states, where it
decomposes spontaneously into two phases of different density and
isospin composition. This phenomenon has been classified as a
non–congruent phase transition SCHRAMM () since there are two
conserved charges, i.e. proton and neutron numbers. Its importance
in the study of the in–medium nuclear interaction has been
emphasized in recent investigations WEISE ().
These two coexisting phases, distinguished in the following by superindices and , have different numbers of proton and neutrons. However they are subject to the thermodynamical equilibrium conditions,
(15)  
(16)  
(17) 
Furthermore, each phase contributes to every intensive additive physical quantity, such as the densities of energy, entropy, etc. Thus, the free energy per unit volume and the density number of nucleons for the whole system can be written as
where the coefficient can be interpreted as the fraction of the partial volume occupied by the state , hence it is bounded by .
In Fig. 1 it is shown how this procedure, commonly known as the Gibbs construction, works for the pressure and the spin asymmetry quotient of each component . The general features of this figure will be discussed in the next section.
Following the standard thermodynamical definitions, the magnetization per unit volume, the magnetic susceptibility and the isothermal compressibility within the coexistence region can be evaluated as
(18)  
(19)  
(20) 
where are the total density of particles and the isospin asymmetry fraction in each phase. The partial contributions to the magnetization , the susceptibility and the compressibility are evaluated in the same way as for a pure single phase.
Iv Results and discussion
For the Skyrme model the SLy4 parametrization is used, for which MeV fm, MeV fm, MeV fm, MeV fm, DOUCHIN (). The saturation density, binding energy, incompressibility and symmetry energy are fm, MeV, MeV and MeV, respectively. Another significative quantity is the inmedium nucleon mass at the saturation density, for which is obtained.
In first place we discuss the effects of the Gibbs construction on the pressure and the spin asymmetry coefficient . These results are shown in Fig. 1, for G and MeV, which is representative for most of the cases studied in this work. The Gibbs construction is shown in dashed lines and replaces, within the coexistence region (CR), the plain results of the model described in Section II. In panel (a) it is shown that the CR includes the instability region where the pressure decreases with density. The Gibbs construction instead, produces a rather linear increasing pressure. The density domain of the CR is reduced and eventually vanishes, by increasing the isospin asymmetry , as well as the temperature (not shown in this figure). For certain values of temperature and magnetic intensity, as for example G, MeV, the CR disappears for asymmetries above a typical value . In these cases and for , a retrograde phase transition takes place. This means that the transition starts and ends at states of similar density and isospin composition, but in between an admixture with states of very different conditions is developed. We illustrate this phenomenon by including afterwards the case G, MeV, .
In the panel (b) of the same figure, we show the spin asymmetry. It can be seen that protons and neutrons are highly polarized at very low densities and they depolarize progressively as the density grows. The coexistence of phases induces an equilibrium state with a significantly reduced degree of polarization, due to the admixture with a higher density and weaker polarization state. For higher magnetic intensities, such as G, the same mechanism causes the frustration of the total neutron magnetization, but it is not able to destroy the magnetic saturation of the proton, as will be discussed subsequently.
The behavior of the pressure as a function of the density for several temperatures and isospin asymmetries is shown in Fig. 2 for G and Fig. 3 for G. The points where a sudden change of slope occurs, correspond to an endpoint of the CR. They do not appear for some particular cases at MeV, because the coexistence does not exist for G and (Fig. 3b), whereas a retrograde transition goes on for , G (Fig. 2b), as explained above. For a given magnetic intensity, thermal effects are more important for lower densities. Furthermore, an increase of the magnetic intensity at constant temperature induces an evident increment of the pressure for , but the opposite trend is observed for . It must be pointed out that the Gibbs construction eliminates all the instabilities for the range of densities and temperatures studied here. In particular there are no regions where the pressure decreases with the density.
The density dependence of the spin asymmetry fraction is shown in Fig. 4 for , several values of , G (Fig. 4a) and G (Fig. 4b). Thermal variations are of no relevance for this quantity. For the lower field intensity, the proton relative polarization is enhanced as decreases, whereas for the neutrons there is only a weak dependence on the isospin composition. For G, the proton component is completely polarized for all the range of and . The effect of the medium polarization is emphasized in this case, as can be seen in the dependence on of the neutron spin asymmetry for . For a fixed total density , the neutron component is completely polarized for and is progressively depolarized as increases. This is a consequence of the dynamical coupling of protons and neutrons (see Eqs. (1)(4)). It must be pointed out that both components, but specially protons in a high sample, tend to recover a high degree of polarization at densities . This feature can be a manifestation of the abnormal spontaneous magnetization described by the Skyrme model at extreme densities.
Results for the free energy per volume as a function of the density, are shown in Fig. 5 ( G) and Fig. 6 ( G). At relatively low densities the kinetic energy and the repulsion between nucleons are small, while the effect of the magnetic field becomes the dominant one. For medium and high densities the repulsive character of the nucleonnucleon interaction and the kinetic energy dominate over the magnetic field and the system increases its energy. This is clearly depicted in Fig. 5, whereas for G in Fig. 6, only the case fits this description. For other values of one should go to higher densities to verify this behavior. From both figures we can see that the addition of protons makes the system to be more bound and so does the increase of the magnetic field. Thermal effects are weak, as opposite contributions tend to cancel each other: the kinetic term increases with , while the entropy contribution does the opposite.
The magnetic susceptibility characterizes the response of the
system to the external field and gives a measure of the energy
required to produce a net spin alignment in the direction of the
field. We have found that this quantity at moderate field
intensities, is sensitive to thermal variations, hence we devote
Figs. 79 to give a more detailed description of the density
dependence of . For all the cases shown, there is a low
density regime where the system has an almost saturated
magnetization (see Fig. 4), therefore the magnetic response is
nearly zero. As it was previously discussed, in most cases the
coexistence of phases frustrates the total magnetization and
consequently enhances the magnetic response within the CR. This
fact can be distinguished by an approximately linear increase,
with a small slope, of the susceptibility as a function of the
total density.
The magnetic susceptibility shows a complex
dependence on the population of the Landau levels. Therefore, in
order to clarify the discussion afterwards, it is worth to make
some general considerations about this relationship. In first
place it must be pointed out that the population of the Landau
levels decreases with both and , but it increases with both
and . However, the distribution function becomes smoother
as the temperature grows, erasing eventually effects due to the
progressive occupation of different Landau levels. In second
place, when very high levels are occupied, further changes in the
quantum number has only imperceptible consequences on the
susceptibility. Finally, the proton component generally shows a
diamagnetic behavior, but it turns to be paramagnetic when the
statistical
occupation is peaked at the lowest Landau level .
Now we focus on the susceptibility for protons in the lower panel
of Fig. 7, where the case of G and MeV, is
illustrated. For increasing densities, the susceptibility changes
from a linear response to an oscillatory behavior. The linear
response turns out from the Gibbs construction, while the
oscillatory behavior is a consequence of the population of the
Landau levels. Note that in the evaluation of the partial
contributions in Eq. (19), the Landau
levels also play a role. However, the narrow range of variation of
the densities and allows the linear behavior for
the susceptibility. For and the system is
diamagnetic, while for it is mainly paramagnetic due to
the fact that only the lowest Landau levels are accessible. In the
upper panel of the same figure, the neutron susceptibility is
depicted. An important change in the scale is observed. The
neutron susceptibility is always paramagnetic and also has a
linear behavior for low densities which changes to oscillatory as
the density grows. Even though neutrons do not have discrete
levels, they are influenced by the dynamical coupling among
protons and neutrons (see Eqs. (1)– (3)).
For the case a sudden decrease at ,
is observed due to a change in the spin configuration of the
system.
A small increase of the temperature, keeping fixed the magnetic intensity, destroys the oscillations at high densities. See Fig. 8 for G, MeV. The magnetic response of protons becomes almost linear for the full range of densities. Two ingredients contribute to these results: the distribution function becomes smoother as the temperature grows and the statistical occupation of the Landau levels increase, being this last effect the dominant one.
For the two figures just described, there is a noticeable difference of magnitude between the susceptibility of protons and neutrons. The proton component is significantly more reactive to changes of the magnetic intensity. The effects of changing the magnetic intensity by keeping fixed the temperature are exposed in Fig. 9, where we present the results for G, MeV. A dramatic change of scale is apparent for both components, but specially for protons. This fact is coherent with the saturation of proton spin alignment (see Fig. 4b), as a consequence the proton component experiences a change of regime from diamagnetic to paramagnetic. In contrast, neutrons which only have a paramagnetic channel, are completely blocked when the saturation of spin alignment is reached. Hence we have , as for . Note that for G, the statistical occupation of the Landau level is strongly dominated by , which explains the paramagnetic behavior of the proton susceptibility.
In Figs. 10 and 11, the isothermal compressibility is presented as
a function of the density for G and G,
respectively. In both figures we show results for , 5 and 10
MeV and isospin asymmetries , 0.4 and 0.8. For all the cases
shown, there is a clear difference between the CR and the domain
of higher densities. In the last case, we found the typical
behavior of an almost incompressible fluid with a slightly
decreasing compressibility. Here the variation of temperature and
isospin composition have weak effects, but they become appreciable
for G. An exceptional behavior is found for
, , MeV, where a local maximum can be seen
around . It seems to be a feature of the model
used and deserves a further investigation.
Within the CR, the compressibility decreases strongly for very low
densities and becomes increasing for moderately higher values. As
a consequence a local maximum of arises at the extreme point
of the CR. This description applies strictly to the lowest
temperature shown here, MeV. As is increased, the effect
becomes less pronounced, but is still perceptible for
G, whereas it is completely missed for G.
As the last matter of analysis we present in Fig. 12 the phase diagram for a fixed value G. Here the closed curves correspond to the isothermal contour of the CR, for both the (Fig. 12a) and (Fig. 12b) planes. The CR for a given temperature is enclosed by the corresponding curve. In the upper panel for each curve there is a maximum value for the pressure, at the left and below that pressure lies the lower density gaseous phase and, at the right the higher density liquid phase. Above that pressure a continuous passage from one phase to another occurs. As the temperature grows the area within the contour is reduced and eventually collapse to the critical point. Note that at very low pressures it is necessary to include states with a small proton excess in order to complete the phase diagram.
V Conclusions
In this work we have addressed the properties of asymmetric nuclear matter under the action of very strong magnetic fields , temperatures below MeV and densities up to twice the saturation density. For the nuclear interaction we have used the non–relativistic Skyrme potential model (SLy4 parametrization) within a Hartree–Fock approximation. We have paid special attention to the low–density and low–temperature domain, where the nuclear interaction gives rise to instabilities, commonly associated with a liquidgas phase transition. If the Coulomb force is not included, the equilibrium state of the system separates spontaneously into coexisting phases. This phase transition was studied in detail in the past SCHRAMM (); WEISE () and it has received renewed attention recently in connection with astrophysical investigations GIULANI () and also in heavy–ion collisions, where the liquid–gas instabilities are related to the formation of fragments occurring in finite nuclear systems DAS ().
Here we introduce an external magnetic field as a new parameter which modifies the properties of the coexisting phases. We propose and analyze the spin polarization, the magnetic susceptibility and the isothermal compressibility as physical quantities that reveal the changes in the global configuration as well as in the internal composition of the equilibrium state. A related investigation, but restricted to neutron matter, has been presented recently in UNLP2 (). There, a comparison of the predictions of very different models was made, which give us some confidence about the common features and warn us about possible singularities of the model chosen.
To obtain the physical equation of state, we have implemented the so called Gibbs construction, which is the appropriate procedure when there are multiple conserved charges.
Referring to the spin polarization, it is higher for very low densities. Neutrons are in general less polarized than protons. At extreme intensities ( G) the effect of the medium polarization is emphasized, due to the dynamical coupling of protons and neutrons (see Eqs. (1)(4)) the relative magnetization of the neutrons is enhanced by increasing the proton fraction. The spin asymmetry depends weakly on the temperature, while it is strongly affected by the magnetic field. The spontaneous separation into independent phases reduces the degree of polarization for both protons and neutrons. This can be understood because for a given total density within the CR, one of such phases has a greater partial density and a lower polarization. This is the reason why neutrons do not reach the magnetic saturation at medium densities under the strongest intensity considered here.
The energy per volume shows that the addition of protons makes the system more bound. And so happens for increasing the magnetic field. At variance, the dependence with temperature is rather small.
The magnetic susceptibility is the most sensitive quantity, as it shows a strong dependence on the isospin asymmetry, the temperature and the magnetic intensity. For protons these dependencies manifest through the population of the Landau levels, which also reflects in the neutron susceptibility due to the dynamical coupling generated by the Skyrme model. The susceptibility in the CR is almost linear.
The isothermal compressibility at very low temperature offers a clear manifestation of the changes in the phase diagram. Out the CR the results corresponds to an almost incompressible fluid, decreasing slowly with density and showing a small dependence on and . Within the CR, the compressibility changes from steeply decreasing at very low densities, to moderately increasing. In such a way a local maximum appears at the end of the CR. This behavior is attenuated by increasing the temperature.
An overview of the general phase diagram shows that the coexistence of phases exists for all the range of magnetic intensities of physical interest. The critical temperature lies above MeV. An increase of the magnetic field has multiple consequences, on one hand it produces an enhancement of the range of pressures involved, while does not modify the density domain. Furthermore the range of isospin asymmetries supported is reduced and for the higher the system evolves through a retrograde phase transition.
It is worthwhile to mention that most of this conclusions become evident because we used an extreme magnetic intensity. For weaker fields, effects such as the temperaturepolarization antagonism or the synergistic combination of neutron excess and spin polarization, are still active but they manifest diffusely.
Acknowledgements
This work was partially supported by the CONICET, Argentina. E. B. acknowledges the support given by the CONICET under contract PIP 0032 and by the Agencia Nacional de Promociones Científicas y Técnicas, Argentina, under contract PICT20102688.
References
 Lai D 2001 Rev. Mod. Phys. 73 629.
 Lattimer J M and Prakash M 2007 Phys. Rep. 442 109.
 Shibanov Y A and Yakovlev D G 1996 Astron. Astrophys. 309 171.
 Yakovlev D G , Kaminker A D, Gnedin O Y and Haensel P 2001 Phys. Rep. 442 1 and references therein.
 Bezchastnov V G and Haensel P 1996 Phys. Rev. D 54 3706; Baiko D A and Yakovlev D G 1999 Astron. Astrophys. 342 192; Chandra D, Goyal A and Goswami K 2002 Phys. Rev. D 65 053003.
 Kharzeev D 2006 Phys. Lett. B 633 260; Kharzeev D E, McLerran L D and Warringa H J 2008 Nucl. Phys. A 803 227.
 Mo Y J, Feng S Q and Shi Y F 2013 Phys. Rev. C 88 024901.
 Skokov V V, Illarionov A Y and Toneev V D 2009 Int. J. Mod. Phys. A 24 5925.
 Chakrabarty S, Bandyopadhyay D and Pal S 1997 Phys. Rev. Lett. 78 2898.
 Broderick A, Prakash M and Lattimer J M 2000 Astroph. J. 537 351.
 Yue P, Yang F and Shen H 2009 Phys. Rev. C 79 025803.
 Sinha M, Mukhopadhyay B and Sedrakian A 2013 Nucl. Phys. A 898 43; Sinha M and Bandyopadhyay D 2009 Phys. Rev. D 79 123001.

Rabhi A, Providência C and Providência J C 2009 Phys. Rev. C 79 015804;
ibid. 2009 80 025806
Rabhi A, Panda P K and Providência C 2011 Phys. Rev. C 84 035803.  Ryu C Y, Cheoun M K, Kajino T, Maruyama T and Mathews G J 2012 Astropart. Phys. 38 25.
 Ryu C Y, Kim K S and Cheoun M K 2010 Phys. Rev. C 82 025804.
 PérezGarcía M A, Providência C and Rabhi A 2011 Phys. Rev. C 84 045803.
 Dong J, Zuo W and Gu J 2013 Phys. Rev. D 87 103010.
 Diener J P W and Scholtz F G 2013 Phys. Rev. C 87 065805.
 Dong J, Lombardo U, Zuo W and Zhang H 2013 Nucl. Phys. A 898 32.
 PerezGarcía M A 2008 Phys. Rev. C 77 065806.
 Isayev A A and Yang J 2009 Phys. Rev. C 80 065801.
 Aguirre R 2011 Phys. Rev. C 83 055804; Aguirre R and Bauer E 2013 Phys. Lett. B 721 136.
 Aguirre R, Bauer E and Vidaña I 2014 Phys. Rev. C 89 035809.
 Bordbar G H and Rezaei Z 2012 Phys. Lett. B 718 1125.
 Bigdeli M 2012 Phys. Rev. C 85 034302.
 Das C B, Das Gupta S, Lynch W G, Mekjian A Z and Tsang M B 2005 Phys. Rep. 406 1.
 Vautherin D and Brink D M 1972 Phys. Rev. C 3 626; Quentin P and Flocard H 1978 Annu. Rev. Nucl. Part. Sci. 28 523.
 Douchin F, Haensel P and Meyer J 2000 Nucl. Phys. A 665 419.
 Bender M and Heenen P H 2003 Rev. Mod. Phys. 75 121.
 Bender M, Dobaczewski J, Engel J and Nazarewicz W 2002 Phys. Rev. C 65 054322.
 Landau L D and Lifschitz E M Quantum mechanics. Nonrelativistic theory, Pergamon Press, Oxford 1991.
 Pathria R K Statistical Mechanics, ButterworthHeinemann, Oxford 1997.
 Hempel M, Dexheimer V, Schramm S and Iosilevskiy I 2013 Phys. Rev. C 88 014906.
 Holt J W and Weise W 2013 Prog. Part. Nucl. Phys. 73 35; Drews M, Hell T, Klein B and Weise W 2013 Phys. Rev. D 88 096011; Wellenhofer C, Holt J W, Kaiser N and Weise W 2014 Phys .Rev. C 89 064009.
 Giulani G, Zheng H and Bonasera A 2014 Prog. Part. Nucl. Phys. 76 116.