Nuclear symmetry energy in a modified quark meson coupling model
We study nuclear symmetry energy and the thermodynamic instabilities of asymmetric nuclear matter in a self-consistent manner by using a modified quark-meson coupling model where the confining interaction for quarks inside a nucleon is represented by a phenomenologically averaged potential in an equally mixed scalar-vector harmonic form. The nucleon-nucleon interaction in nuclear matter is then realized by introducing additional quark couplings to , , and mesons through mean-field approximations. We find an analytic expression for the symmetry energy as a function of its slope . Our result establishes a linear correlation between and . We also analyze the constraint on neutron star radii in matter with equilibrium.
pacs:26.60.+c, 21.30.-x, 21.65.Qr, 95.30.Tg
One of the major focuses in the study of nuclear matter has recently been to understand the equation of state (EOS) of asymmetric nuclear matter and the density dependence of the nuclear symmetry energy. The nuclear symmetry energy is a fundamental quantity which determines several important properties of very small entities such as the atomic nuclei as well as very large objects such as neutron stars prakash (). In fact, the behavior of nuclear symmetry energy is most uncertain among all properties of dense nuclear matter. Furthermore, the symmetry energy is important for modeling nuclear matter by probing the isospin part of nuclear interactions. Recent tili () experimental studies of isospin-sensitive observables in intermediate-energy nuclear reactions involving radioactive beams have been quite useful in providing some constraints on the density dependence of nuclear symmetry energy at subsaturation densities. The effects of symmetry energy and its slope on neutron star properties is an important area of study. Another area of relevance in the study of asymmetric nuclear matter is the instabilities associated with possible liquid-gas phase transitions at subsaturation densities. Such liquid-gas phase transition plays an important role in the description of the crust of compact star matter at densities between fm and saturation density ( fm). Here, we would like to address these two relevant aspects in the study of asymmetric nuclear matter in a phenomenological model that we have used in our earlier work for symmetric nuclear matter.
There has been a proliferation of phenomenological models to describe infinite nuclear matter and also properties of finite nuclei. These are in fact essential steps in the development of this area of study for which realistic first-principles theoretical descriptions as well as adequate experimental or observational data are not available. Variations in different phenomenological approaches stretching from the nonrelativistic to the relativistic are tried incorporating some further aspects of theoretical requirements in the model. All these models are usually set in terms of parameters that are fit to reproduce the properties of either finite nuclei nl3 () or bulk nuclear matter. As a result, most of the models behave more or less similarly as far as the equation of state is concerned around the saturation density and at zero temperature. However, when these models are used to describe nuclear matter at subsaturation densities to explain the liquid-gas phase transition or at high densities to explain neutron star matter, they yield very different results. Therefore, it has been seen as essential to incorporate some constraints related to symmetry energy and its derivatives by using up-to-date theoretical and experimental information. Most of the relativistic-mean-field (RMF) models are attempts in these directions. However, in these models nucleons are treated as structureless point objects. Therefore, as a next step in the requirement of incorporating the quark structure of the nucleon with meson couplings at the basic level; quark-meson-coupling (QMC) models have been proposed guichon () and properties of nuclear matter have been studied in great detail in a series of works ST (); recent (); temp (); phase (). In these models nucleons are described as a system of nonoverlapping MIT bags which interact through effective scalar and vector meson exchanges at the quark level. However, it has been argued that the hadronic structure described by the MIT-bag model suffers from some theoretical inadequacy due to the sharp bag boundary in breaking chiral symmetry, which is a good symmetry of strong interactions within the partially conserved axial current (PCAC) limit. Therefore, more sophisticated versions such as the Cloudy Bag Model (CBM) have been proposed for the study of hadronic structure. So, to further include this aspect of the physics requirement, it would be more appropriate to develop a quark-meson coupling model where nucleon structure is described by models like the CBM instead of by MIT bags. As an alternative approach barik (); frederico89 (); batista () to the CBM, the relativistic independent quark model with a phenomenologically averaged confining potential in equally mixed scalar-vector harmonic form in the Dirac frame work has been used extensively with remarkable consistency in the baryonic as well as the mesonic sector bdd (). This model has provided a very suitable alternative to the otherwise successful cloudy bag model in describing hadronic structure with its static properties and various decay properties.
We therefore proposed in our earlier work rnm () a modified quark-meson coupling model (MQMC), which is based on a suitable confining relativistic independent quark potential rather than a bag to address the nucleon structure in vacuum as an alternative approach to QMC for the study of symmetric nuclear matter. This attempt was not so much as to plead superiority of MQMC over QMC at this level. It has only incorporated an essential aspect of the physics requirement missing in MIT-bag model to stand as an alternative to the more appropriate CBM. Further investigations are necessary to check its consistency and its predictability for any new physical features. In the MQMC model, we have studied the bulk nuclear properties such as the compressibility, the structure of EOS, and also discussed some implications of chiral symmetry in nuclear matter along with the nucleon and nuclear term and the sensitivity of nuclear matter binding energy with variations in the light quark mass. The results obtained in such a picture for symmetric nuclear matter were quite encouraging. In the present attempt, we study the bulk properties of asymmetric nuclear matter and also the low-density instabilities of the system in such a model. To treat the asymmetric nuclear matter, we incorporate in our model the contribution of the isovector vector meson in addition to those of the isoscalar scalar meson and isoscalar vector meson considered earlier for symmetric nuclear matter rnm (). Such studies are also useful to discuss the systems such as neutron stars with .
A correlation between the symmetry energy and its slope has been verified recently by Ducoin et al. ducoin () for a set of effective relativistic and nonrelativistic nuclear models. Such a study was based on numerical results for and obtained from different parametrizations. Theoretically and are constrained stone (); chen (). In a recent paper, Santos et al. correl () have established an analytic relationship between these quantities. In this context QMC-based models have not been studied. We have made an attempt to set up a relationship between these two quantities analytically.
The paper is organized as follows: In Sec. II, a brief outline of the model describing the nucleon structure in vacuum is discussed. The nucleon mass is then realized by appropriately taking into account the center-of-mass correction, pionic correction, and gluonic correction. The EOS is then developed. In Sec. III, we discuss the nuclear symmetry energy, its slope and incompressibility, and observe its density dependence. The thermodynamic instabilities of the system are analyzed in Sec. IV. We establish the analytic relationship between and and discuss the results in Sec. V.
Ii Modified quark meson coupling model
Recently, the modified quark-meson coupling model was adopted for symmetric nuclear matter where the interaction was realized in a mean-field approach through the exchange of effective mesonic fields coupling to the quarks inside the nucleon rnm (). We now extend this model to asymmetric nuclear matter and include the contribution of the isovector vector meson, , in addition to and mesons. In view of this, we briefly present the outlines of our approach rnm () in the present context.
We first consider nucleons as a composite of constituent quarks confined in a phenomenological flavor-independent confining potential, in an equally mixed scalar and vector harmonic form inside the nucleon rnm (), where
Here are the potential parameters. The confining interaction here provides the zeroth-order quark dynamics of the hadron. In the medium, the quark field satisfies the Dirac equation
where , and ; where , , and are the classical meson fields, and , , and are the quark couplings to the , , and mesons, respectively. is the quark mass and is the third component of the Pauli matrices. In the present paper, we consider nonstrange and quarks only. We can now define
where the effective quark energy, and effective quark mass, . We now introduce and as
The ground-state quark energy can be obtained from the eigenvalue condition
The solution of equation (5) for the quark energy immediately leads to the mass of the nucleon in the medium in zeroth order as
We next consider the spurious center-of-mass correction , the pionic correction for restoration of chiral symmetry, and the short-distance one-gluon exchange contribution to the zeroth-order nucleon mass in the medium. The center-of-mass correction and the pionic corrections in the present model are found, respectively, as rnm ()
with the axial vector nucleon form factor given as
The pseudovector nucleon pion coupling constant can be obtained from the familiar Goldberg Triemann relation by using the axial-vector coupling-constant value in the model, as discussed in Ref. rnm ().
The color-electric and color-magnetic contributions to the gluonic correction which arises due to one-gluon exchange at short distances are given as:
and due to color-magnetic contributions, as
where and are the numerical coefficients depending on each baryon. The color-electric contributions to the correction of baryon masses due to one gluon exchange are calculated in a field-theoretic manner rnm (). It can be found that the numerical coefficient for color-electric contributions such as , and comes out to be zero. From calculations we have and for the nucleons. The quantities and are given in the following equation
In the calculation we have taken as the strong-coupling constant in QCD at the nucleon scale barik (). The color-electric contribution is zero here, and the gluonic corrections to the mass of the nucleon are due to color-magnetic contributions only.
Finally, treating all these corrections independently, the mass of the nucleon in the medium becomes
The total energy density and pressure at a particular baryon density for the nuclear matter becomes
where is the spin degeneracy factor for nuclear matter. The nucleon density becomes
Therefore, the total baryon density becomes and the (third component of) isospin density . The proton fraction, is defined as
where and are the proton and neutron densities.
The vector mean-fields and are determined through
where and . Finally, the scalar mean-field is fixed by
The iso-scalar scalar and iso-scalar vector couplings and are fit to the saturation density and binding energy for nuclear matter. The isovector vector coupling is set by fixing the symmetry energy. For a given baryon density, , , and are calculated from Eqs. (20) and (21), respectively.
Iii The symmetry energy
We may define the neutron-excess parameter so that the nuclear symmetry energy can be obtained as the difference between the total energy per baryon of pure neutron matter and that of isopspin-symmetric matter at baryon density . Here we consider the nuclear matter consisting of protons and neutrons only with as the proton fraction. An expansion of the total energy per baryon, , with respect to the neutron-excess parameter, becomes steiner ()
where , etc. are the first-, second- and third- order derivatives with respect to in a Taylor’s expansion. However, neglecting Coulomb contributions near the isospin symmetry of QCD, demands the total energy of pure neutron matter to be same as that of pure proton matter, for which the odd powers in are to be forbidden in the above expansion. Again for densities near or below the saturation density (), truncation of this expansion to quadratic terms in is considered to be a good approximation. In view of that, the coefficient of the quadratic term in can be identified as the symmetry energy
This may be considered to be a good approximation even for small proton fraction , which can be valid for finite nuclei. But for nuclear matter at densities in excess of the saturation density , effects of higher order than quadratic in the expansion may be important. Therefore, in order to study the density dependence of , one may expand this as a function of around saturation density in terms of a parameter to obtain
so as to consider the symmetry-energy parameters as follows:
Apart from the quantities in Eq. (LABEL:cesym), the following quantities calculated from pressure and energy density for the consideration of constraints and correlations studies are
and the volume part of the iso-spin incompressibility
We have assumed since the volume term is dominant dutraetal (). These parameters characterize the density dependence of nuclear symmetry energy around normal nuclear matter density and thus provide important information on the behavior of nuclear symmetry energy at both high and low densities. Also, the curvature parameter distinguishes the different parametrizations. A more significant measurement would be the evaluation of the shift of the incompressibility with asymmetry, which is given by
because this value can be correlated to experimental observations of the giant monopole resonance (GMR) of neutron-rich nuclei. Recent observations of the GMR tili () on even- Sn isotopes give a quite stringent value of MeV. In the present model we determine this value for three quark masses of , , and MeV and observe that they are consistent with the GMR measurements.
The compressibility at saturation density can be determined analytically, from Eq. (LABEL:ceng):
The study of the correlation between symmetry energy and its slope can be performed analytically. For this purpose we use the Eqs. (16) and (23) to find . In this model, we get the closed-form expression
The correlation function exhibits the dependence on the different bulk parameters , and .
Iv Stability Conditions
Nuclear forces have an attractive long-range part and a repulsive hard core similar to a Van der Waals fluid. It is expected to present a liquid and a gas phase characterized by the respective densities. Nucleons can be either protons or neutrons. Such a two-component system undergoes liquid-gas phase transition. The asymmetric nuclear matter (ANM) shows two types of instabilities chomaz (): a mechanical instability conserving the proton concentration and a chemical instability occurring at constant density.
We consider asymmetric nuclear matter characterized by proton and neutron densities and transform these into a set of two mutually commuting charges unique (). In infinite matter the extensivity of free-energy implies that it can be reduced to a free energy density: which at reduces to energy density only. Since, we deal with a two-component nuclear medium, spinodal instabilities are intimately related to phase equilibria and phase transitions. Although it consists of unstable states, the spinodal region of the phase diagram can be addressed by standard thermodynamics.
The condition for stability implies that the free energy density is a convex function of the densities . A local necessary condition is the positivity of the curvature matrix:
Here we used , where the effective chemical potentials in the present context are given by
Since we consider a two-fluid system, is a symmetric matrix with two real eigenvalues ms (). The two eigenvalues are given by,
and the eigenvectors by
The largest eigenvalue is always positive whereas the other can take on negative value. We are interested in the latter, because it defines the spinodal surface, which is determined by the values of , and . The smallest eigenvalue of becomes negative. The associated eigenvector defines the instability direction of the system, in isospin space.
V Results and Discussion
We set the model parameters by fitting the nucleon mass MeV and charge radius of the proton fm in free space. Taking standard values for the meson masses; namely, MeV, MeV and MeV and fitting the quark-meson coupling constants self-consistently, we obtain the correct saturation properties of nuclear matter binding energy, MeV, pressure, , and symmetry energy MeV at fm. The values of , , and obtained this way and the values of the model parameters at quark masses , , and MeV are given in Table 1.
|MQMC (5 MeV)||-15.7||0.151||0.93||159||32.0||84.7||-27.7||-535.9||103.2||-590.8|
|MQMC (40 MeV)||-15.7||0.151||0.91||208||32.0||84.9||-28.4||-537.6||94.2||-575.9|
|MQMC (300 MeV)||-15.7||0.151||0.76||349||32.0||89.1||-14.5||-549.0||-15.6||-545.1|
|DD typel ()||-16.0||0.149||0.56||239||31.6||55.9||-95.3||-431.1||576.8||-462.57|
|QMC alex (); alextemp ()||-15.7||0.150||0.77||291||33.7||93.5||-10.0||-570.8||29.4||-580.24|
|FSUGold todd ()||-16.2||0.148||0.61||229||32.6||60.4||-51.4||-414.0||425.7||-276.07|
|BKA24 bka24 ()||-15.9||0.147||0.60||227||34.2||84.8||-14.9||-523.7||112.4||-421.55|
|BSR12 bsr12 ()||-16.1||0.147||0.61||232||34.0||77.9||-44.2||-511.6||324.2||-414.30|
In Figs. 1(a) and 1(b), we plot the binding energy per nucleon for nuclear matter as a function of density corresponding to MeV and MeV, respectively, for different values. In Fig. 2, we compare the variation of the binding energy per nucleon for quark mass and MeV with that of QMC and observe that, for MeV, MQMC compares well with that of QMC. It is observed from Figs. 1 and 2 that, at low quark mass, the equation of state is softer. In Table 2, we compare the nuclear matter properties at saturation for quark masses , , and MeV, respectively, in the present model to QMC alex (); alextemp (), and some of the approved models as suggested in Ref. dutraetal ().
The value of the compressibility is determined to be , , and MeV respectively, for quark masses , , and MeV. A recent calculation stonemos () has predicted to be in the range from the experimental GMR energies in even-even Sn and Cd. Furthermore, the value of the effective mass calculated in the present model at quark mass MeV is which compares well with the empirical value of the effective mass, which is mahaux ().
We observe that the symmetry energy shows an extremely linear behavior. This is further justified from the plot for the slope parameter . This is based on equation (23) for . However, if we consider terms higher than the quadratic one in defining the relation in Eq. (22), it would be more appropriate to use the expression as in Eq. (24) to show the density dependence of for higher nuclear densities . This has been shown in Fig. 4 in comparison with several other models, as noted there, including QMC.
v.1 Correlation between the symmetry energy and its slope
We study the correlation function with the variation of quark masses in Fig. 5. We observe that the function increases with quark masses.
The established value of binding energy and the saturation density in the nuclear mean-field models cannot be applied to incompressibility and effective mass since the latter are found as output in this model where the coupling constants are fixed in a self consistent manner by taking into consideration the binding energy and saturation density. Therefore, we have taken the variation of with different quark masses at the same and . We observe that, because there is only one isovector parameter in the expression for and , the variation is linear. Such linearity in the behavior was also observed in nonrelativistic models correl (). It indicates one of the limitations of the model parameters. We expect a nonlinear behavior between and for the models with more than one isovector parameter. In the Fig. 6, we have shown the slope of symmetry energy, at various quark masses. It is interesting to note that there is a linear relationship of the slope of the symmetry energy with quark mass. This is a direct consequence of the dependence of the symmetry energy on .
We next study the mechanical instability and its dependence on the isospin asymmetry of the system by plotting the pressure as a function of density and the asymmetry parameter . In Fig. 7, we show that the mechanical instability occurs in the region where the slope of the pressure with respect to density is negative. We observe that the mechanical-instability region shrinks when the isospin-asymmetry increases.
The system is stable under separation of two phases if the free energy of a single phase is lower than the free energy in all two-phase configurations.
In the spinodal area we can get the signature of the mechanical instability by finding the velocity of sound determined in the model as
where , is the velocity of sound and is the speed of light.
In Fig. 8 we show the sound velocity as a function of density by changing the asymmetry parameter. There is a reduction in the instability when we move away from the symmetric nuclear matter. Moreover, the velocity becomes imaginary when we enter into the spinodal area unique ().
The positivity of the local curvature matrix is equivalent to the condition that both the trace and the determinant are positive. In the present model the above condition is violated and the system is in the unstable region of a phase transition. Further it is to be pointed out that for a two component, n-p thermodynamical system, the stability parameter is given by the condition;
but in charge symmetric matter the isoscalar (total density) and isovector (concentration) oscillations are not coupled and there are two separate conditions for instability asym (). These conditions are for mechanical instability
and for chemical instability
In the ANM the isoscalar and isovector modes are coupled and the two separate inequalities do not select the nature of instability. Moreover, we observe in Fig. 9, a large difference in the behavior of the stability parameter in Eq. (37) inside the instability region. For higher asymmetry, the range of the stability parameter is smaller than at lower asymmetries. To understand this effect we follow the Landau-dispersion-relation approach for small-amplitude oscillations in Fermi liquids landau (). For a two component matter, the interaction is characterized by the Landau parameters which is defined by the relationship
where and represents the single-particle level density at the Fermi energy. At zero temperature it has the simple form
The dispersion relations give the properties of density (isoscalar) modes and gives the concentration (isovector) modes.
In the unstable region of dilute asymmetric nuclear matter we have isoscalar-like unstable modes, hence , while the combination . In the Fig. 10 we plot the values obtained from the calculation of these two quantities in the unstable region at zero temperature for . An important observation we make from the comparison of Figs. 9 and 10 is the shift in the maximum instability density region. In Fig. 9 the largest instability (the most negative value) is at . However, the most negative Pomeranchuk condition , which corresponds to the fastest unstable mode, is present in more dilute matter at .
In the following we study the direction of instability of the system. In Fig. 11, we show the ratio of the proton versus neutron density fluctuations corresponding to the unstable mode which defines the direction of instability of the system. We plot the results for different proton fractions and observe that the instabilities tend to restore the isospin symmetry for the dense (liquid) phase leading to the fractionation of the ANM.
Figure 12 shows the proton-neutron density fluctuation ratio as a function of the isospin asymmetry for a fixed nuclear density, and compares it to the QMC and Brueckner-Hartree-Fock (BHF) calculations. The relativistic models give larger fluctuation ratios than the corresponding value of . We also observe that the fluctuation ratio in the present model is larger compared to the nonrelativistic BHF model. A pure mechanical disturbance would occur chomaz () if the instability preserves the ratio between protons and neutrons, i.e., . Conversely if then we should observe pure chemical disturbance. In the present case we observe that the disturbance along the unstable eigen direction conserves neither nor but has mixed character with both chemical and mechanical contents.
v.3 Constrain on neutron star radii
The symmetry energy plays an important role in describing the mass-radius relationship in neutron stars. Neutron stars are compact objects maintained by the equilibrium of gravity and the degenerecy pressure of the fermions together with a strong nuclear repulsion force due to the high density reached in their interior. The slope of the symmetry energy, , constrains the neutron star radii. It is confirmed that the radii for the neutron stars with canonical mass are not affected by the symmetry energy at saturation density menezes (). However, in some cases the radii increase with , while in others, there is a decrease. In fact the radii are correlated with a variation of the slope . The radii increase up to a maximum value, then drop again. This behavior can be associated with a maximum theoretical value of , and provide a possible constraint to nuclear matter. In the present model the value of comes out to be MeV which is very close to the experimental observation tili (). The most direct connection between the astrophysical observations and the nuclear symmetry energy concerns neutron star radii () which are highly correlated with neutron star pressures near . It is to be noted that Lattimer and Prakash constrain () found the radii of neutron stars for masses near the canonical mass , obey a power-law relation:
where is the radius of a star, is the pressure of neutron star matter at density , and is a constant for a given density and mass. Considering the structure of a neutron star as pure neutron matter, the value of this constant at quark mass MeV in our model comes out to be
which is very near to km predicted by Lattimer et al. constrain (). It is to be noted that the mass of the neutron star for matter with equilibrium comes out to be 2.7 with quark mass MeV and 1.64 with quark mass MeV. The details of such calculations incorporating the hyperons in the composition of neutron stars is in progress.
In the present work we have studied the EOS for asymmetric nuclear matter by using a modified quark-meson coupling model (MQMC). Self-consistent calculations were made by using a relativistic quark model with chiral symmetry along with the spurious center-of-mass correction, pionic correction for restoration of chiral symmetry, and short-distance correction for one-gluon exchange to realize different bulk nuclear properties. The instability in the two-component nuclear system is then analyzed. In asymmetric matter the isoscalar and isovector modes are coupled and the two separate inequalities for density oscillations and concentration oscillations no longer maintain a physical meaning for the selection of the nature of the instabilities.
The symmetry energy, its slope , and curvature parameter were found in reasonable agreement with experimental values. Without considering self interactions in the scalar field, we found an analytic expression for the symmetry energy as a function of its slope . Our result establishes a linear correlation between and . We also study the variation of correlation function with the variation of quark masses. The symmetry energy is correlated with neutron star radii. In this model we observe that, at twice the saturation density , the constant is found 5.90 km in the canonical-mass region of .
The authors would like to acknowledge the financial assistance from BRNS, India for the Project No. 2013/37P/66/BRNS.
- (1) M. Prakash, I. Bombaci, M. Prakash, P. J. Ellis and J. M. Lattimer, Phys. Rep. 280, 1 (1997).
- (2) T. Li, et al., Phys. Rev. Lett. 99 162503 (2007); Bao-An Li, Lie-Wen Chen, Che Ming Ko, Phys. Reps. 464, 113 (2008).
- (3) G.A. Lalazissis, J. König, P. Ring, Phys. Rev. C 55, 540 (1997).
- (4) P. A. M. Guichon, Phys. Lett. B 200, 235 (1988); P.A.M. Guichon, K. Saito, E. Rodionov, A.T. Thomas, Nucl. Phys. A 601 (1996) 349.
- (5) K. Saito and A.W. Thomas, Phys. Lett. B 327, 9 (1994); 335, 17 (1994); 363, 157 (1995); Phys. Rev. C 52, 2789 (1995); K. Saito, K. Tsushima, and A.W. Thomas, Nucl. Phys. A 609, 339 (1996); Phys. Rev. C 55, 2637 (1997); Phys. Lett. B 406, 287 (1997).
- (6) P. G. Blunden and G.A. Miller, Phys. Rev. C 54, 359 (1996); N. Barnea and T.S. Walhout, Nucl. Phys. A 677, 367 (2000); H. Shen and H. Toki, Phys. Rev. C 61, 045205 (2000); P.K. Panda, R. Sahu, C. Das, Phys. Rev. C 60, 38801 (1999); P.K. Panda, M.E. Bracco, M. Chiapparini, E. Conte, and G. Krein, Phys. Rev. C 65, 065206 (2002); P.K. Panda, and F.L. Braghin, Phys. Rev. C 66, 055207 (2002).
- (7) P.K. Panda, A. Mishra, J.M. Eisenberg, and W. Greiner, Phys. Rev. C 56, 3134 (1997); I. Zakout, and H.R. Jaqaman, Phys. Rev. C 59, 962 (1999).
- (8) G. Krein, D.P. Menezes, M. Nielsen, and C. Providencia, Nucl. Phys. A 674, 125 (2000), P.K. Panda, G. Krein, D.P. Menezes and C. Providência, Phys. Rev. C 68, 015201 (2003); P.K. Panda, D.P. Menezes, C. Providência. Phys. Rev. C 69, 025207 (2004).
- (9) T. Frederico, B. V. Carlson, R. A. Rego, M. S. Hussein. J. Phys. G 15, 297 (1989).
- (10) N. Barik and B.K. Dash, Phys. Rev. D 33, 1925 (1986), ibid, Phys. Rev. D 34, 2092 (1986); N. Barik and R.N. Mishra, Phys. Rev. D 61, 014002 (2000),
- (11) E.F. Batista, B.V. Carlson and T. Frederico, Nucl. Phys. A 697, 469 (2002).
- (12) N.Barik, B.K.Dash and M.Das, Phys. Rev. D 31, 1652 (1985) ibid 32, 1725 (1985);34, 2803 (1986). N.Barik, P.C.Dash, A.R.Panda, Phys. Rev. D 46, 3856 (1992); N.barik, P.C. Dash, Phys. Rev. D 49, 299 (1994).
- (13) N. Barik, R. N. Mishra, D. K. Mohanty, P. K. Panda and T. Frederico, Phys. Rev. C 88 015206 (2013).
- (14) C. Ducoin, J. Margueron, C. Providenĉia and I. Vidãna, Phys. Rev. C 83, 045810 (2011).
- (15) J. R. Stone and P.G. Reinhard, Prog. Part. Nucl. Phys. 58, 587 (2007).
- (16) L. W. Chen, C.M. Ko, Bao-An Li, and J. Xu, Phys. Rev. C 82, 024321 (2010).
- (17) B. M. Santos and M. Dutra,O. Lourenco and A. Delfino Phys. Rev. C 90 035203 (2014).
- (18) A. W. Steiner, Phys. Rev. C 74 045808 (2006)
- (19) M. Dutra et al, Phys. Rev. C 90, 055203 (2014)
- (20) P. Chomaz, M. Colonna and J. Randrup, Phys. Rep. 389 263 (2004).
- (21) J. Margueron and P.Chomaz, Phys. Rev. C 67 041602 (2003).
- (22) H. Müller and B. D. Serot, Phys. Rev. C 52 2072 (1995).
- (23) S.Typel, Phys. Rev. C 71,064301(2005)
- (24) A.M. Santos C. Providência, P.K. Panda, Phys. Rev. C 79, 045805 (2009).
- (25) B.G.Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett. 95, 122501 (2005)
- (26) B.K.Agrawal, Phys. Rev. C 81, 034323(2010)
- (27) S. K. Dhiman, R. Kumar and B. K. Agrawal, Phys. Rev. C 76, 045801 (2007)
- (28) P.K. Panda, A.M. Santos, D.P. Menezes, C. Providência, Phys. Rev. C 85, 055802 (2012).
- (29) J.R.Stone, N. J. Stone and S. A. Moszkowski, Phys. Rev C 89, 044316 (2014)
- (30) M. Jaminon and C. Mahaux, Phys. Rev. C 40, 354 (1989)
- (31) B. Liu, V. Greco, V. Baran, M. Colonna, and M. Di Toro, Phys. Rev. C 65, 045201 (2002).
- (32) L.D. Landau and E.M. Lifshitz, Statistical Physics (Pergamon Press Ltd., Oxford, England, 1981) page 208; G. Baym and C. Pethick, Landau Fermi-liquid theory (John Wiley and Sons, New York, 1991); L.Scalone, M.Colonna and M.Di Toro, Phys. Lett. B 461, 9 (1999).
- (33) L.L.Lopes and D.P.Menezes, Braz. J. Phys. 44, 774 (2014).
- (34) J. M. Lattimer and Y. Lim, Astrophys. J. 771 51 (2013).