Warm dense matter and cooling of supernovae remnants

We study the thermal effects on the nuclear matter (NM) properties such as binding energy, incompressibility, free symmetry energy and its coefficients using NL3, G3 and IU-FSU parameter sets of relativistic mean-field models. These models being consistent with the properties of cold NM, have also been used to study the effect of temperature by incorporating the Fermi function. The critical temperature for the liquid-gas phase transition in the symmetric NM is found to be 14.60, 15.37 and 14.50 MeV for NL3, G3 and IU-FSU parameter sets respectively, which is in excellent agreement with previous theoretical and experimental studies. We inspect that the properties related to second differential coefficient of the binding energy and free symmetry energy at saturation density ( i.e. K0(n,T)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_{0}(n,T)$$\end{document} and Qsym,0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Q_{sym,0}$$\end{document} ) exhibit the contrary effects for NL3 and G3 parameters as the temperature increases. We find that the prediction of saturated curvature parameter ( Ksym,0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_{sym,0}$$\end{document} ) for G3 equation of state at finite temperature favour the combined analysis of Ksym,0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_{sym,0}$$\end{document} for the existence of massive pulsars, gravitational waves from GW170817 and NICER observations of PSR J0030+0451. Further, we investigate the cooling mechanism of newly born stars through neutrino emissivity controlled by direct Urca process and instate some interesting remarks about neutrino emissivity. We also deliberate the effect of temperature on the M-R profile of Proto-Neutron star.


Introduction
One of the most prominent energetic events of the universe is manifested by the core-collapse supernovae (CCSN) explosion of the giant stars having mass in the range of 8-40 times a e-mail: ankit.k@iopb.res.in b e-mail: harish.d@iopb.res.in (corresponding author) c e-mail: patra@iopb.res.in that of the mass of sun [1]. The concept that the gravitational instability during the evolution of a massive star results in a rapid compression and then thermonuclear explosion, was first suggested by Burbidge, Burbidge, Fowler, and Hoyle (designated as "B 2 F H") [2]. The energy emerge during this explosion is carried off by the photons and the neutrinos, which are billion trillion trillion in numbers and hauled the most of the energy released. The observation of Einstein Observatory (HEAO-2) first supported the fact that high neutrino emissivity is mainly responsible for the rapid cooling of newly born dense star [3]. Studies suggest that the neutrinos also play a significant role in the dynamics of supernovae explosion and control many important aspects of the collapsed core, generally known as newly born neutron star [4,5]. In this exploring era of every tiny and massive object through modern science, neutron stars still hold the mystery in itself and are poorly known objects. These are the ideal play field in the observable universe to inspect the theories of dense matter physics and unfold new opportunities. A large part of our perception about the dynamics of supernovae explosion and the composition of neutron stars is based on the equation of state (EoS), which makes the EoS a key ingredient of our study. Since it is believed that the matter evolved during the core-collapse phenomenon is unable to attain β−equilibrium condition for few seconds [6], due to the promptness of the CCSN explosion, so, usually the EoS for a dense matter with either same amount of protons and neutrons or a proton fraction of ∼ 0.2 is taken into account to cipher the dynamics of collapse event [7]. Many other important aspects of core-collapse events like how many protons are converted into neutrons, quenching rate, the mass-radius profile of newly born neutron star and its composition are determined by the EoS [8].
Another important dimension which makes the EoS of nuclear matter (NM) at finite temperature more interesting is the study of the dynamics of heavy ion collision reaction and structure of exotic nuclei. Under extreme conditions of density, isospin and temperature the EoS can be explored through many experimental facilities such as GANIL -SPIRAL2 [9] facility in France, CSR in China [10] and the FRIB technique in the United States [11]. Nuclear collisions can produce hot and dense hadronic NM in terrestrial laboratories momentarily, which can procure a lot of information about the thermodynamical properties of NM like incompressibility, symmetry energy and its derivatives [12]. In heavy-ion collision experiments, the symmetry energy plays an important role in understanding the isoscaling behaviour of multifragmentation phenomena. The observations of the experiments at Texas A&M University and National Superconducting Cyclotron Laboratory suggest a connection between the nuclear symmetry energy and primary fragment yield distribution of statistical multi-fragmentation model [13][14][15]. On the other hand, symmetry energy at finite temperature also plays an important role in the cooling mechanism of newly born hot astrophysical objects [16][17][18] and has effective impact on the cooling rate through direct Urca and modified Urca processes [19,20].
In the recent years, the density dependence of the symmetry energy and its derivatives (L sym , K sym , Q sym ) have been used to constrain the EoS near saturation density, which make them more vital to decipher properly. It is well known that the slope parameter (L sym ) and the size of the neutron skin in super-heavy nuclei are connected by a strong linear correlation. The slope (L sym ) and curvature parameter (K sym ) also control the location of the neutron drip line, core-crust transition density and gravitational binding energy of neutron star [21]. The isovector skewness parameter (Q sym ) is the most ambiguous quantity, due to the large fluctuations in its value obtained from different models. Recent studies show that Q sym is related to the incompressibility [22] of the system and also suggest an important role of Q sym in the cooling of newly born proto-neutron star [8]. We study in detail the effects of temperature on all these nuclear parameters and the correlations among them using the most familiar NL3, IU-FSU and recently developed G3 parameter sets. With thermal effects being the main focus of our study, we also explore the dependence of thermal index (Γ th ) on the density for all the parameter sets. In the extant work, we also investigate the dependence of cooling mechanism of a hot dense matter on the EoS and the variation in the mass-radius profile of a proto-neutron star with temperature. The temperature in the interior of a newly born dense star just after the supernova explosion can vary from 10 to 100 MeV [23]. The analyzation of the x-rays emitted during the stellar evolution of the young star ensure the fact that it losses most of its energy through extremely rapid neutrino emission enhanced by the direct Urca process (named after a casino in Rio) [24][25][26]. However, the direct Urca process dominates only during the initial stage of the cooling and can no longer operate once the proton fraction inside the core reaches a threshold value. The proton fraction in a newly born proto-neutron star depends on how the nuclear symmetry energy scales with density-which is an important aspect of this work [24]. In the mean-field approximation, the symmetry energy and the magnitude of the neutrino emissivity (Q) seem to be controlled by the EoS of different parameter sets.
The paper is organized as: In Sect. 2, we explain the formalism for temperature dependent quantum hadron-dynamics model. Section 3 is divided into two parts, where 3.1 is devoted to the theoretical structure of the relativistic meanfield formalism at finite temperature and acquaintance with the NM parameters. In Sect. 3.2 we presented the results for various NM parameters of symmetric NM using NL3, G3 and IU-FSU parameter sets. The detailed framework and the results of neutrino emissivity through direct Urca process are extended in Sect. 4. The mass-radius profile of protoneutron star is discussed in Sect. 5. Finally, the discussion and the concluding remarks are outlined in the Sect. 6.

Temperature dependent QHD model
Temperature dependent QHD model is based on relativistic covariant field theory, proposed to obtain the expediential in-medium nuclear properties around the saturation density. The nucleons (neutrons and protons) are represented by the Dirac spinor ψ in the Lagrangian density. We take in account the σ , ω, ρ and δ mesons in our Lagrangian, represented by the σ , ω μ , ρ μ and δ fields respectively. The mesons act as the mediators and represent the effective nuclear interaction of the nucleons through the meson fields σ , ω μ , ρ μ and δ. The basic Lagrangian density used in the present work, which includes the interaction of nucleon fields ψ with a scalar field (σ ), a vector field (ω) and isovector fields (ρ μ and δ) and the cross-coupled interactions of these meson fields up to fourth order [27], is given by where M is the nucleon mass; m σ , m ω , m ρ and m δ are the masses of mesons and g σ , g ω , g ρ and g δ are the coupling con-stants for the σ , ω, ρ and δ mesons respectively; κ 3 (or κ 4 ) and ζ 0 are the self-interacting coupling constants of the σ and ω mesons respectively; η 1 , η 2 , η ρ and Λ ω are the coupling constants of non-linear cross-coupled terms. The quantities F μν and R μν being the field strength tensors for the ω and ρ mesons respectively, defined as F μν = ∂ μ ω ν − ∂ ν ω μ and R μν = ∂ μ ρ ν −∂ ν ρ μ . The τ are the Pauli matrices and behave as the isospin operator, which carry the isospin component of the nucleons. In the further calculations of the present work, we will consider relativistic mean-field approximation and isotropy in isospin space. Within the assumed approximation, the meson fields are replaced by their expectation values and for rotationally invariant systems only the time component (μ = 0) of the isovector field survives [28,29]. The third component of the isospin operator (τ 3 ) when operates on neutron and proton gives, τ 3 | p >= (+1)| p > and τ 3 |n >= (−1)|n >. Applying the relativistic mean-field approximation, the Lagrangian density yields the following Dirac equation for the nucleon field Assuming NM as an uniform system, where the fields of σ , ω and ρ mesons are independent of position and taking the thermodynamic argument of free gas into consideration, the baryon, scalar and isovector densities for the appropriate Fermi momentum at finite temperature can be evaluated as [30][31][32] n = α= p,n ψ α ψ α = n p + n n = α= p,n where n p and n n are the proton and neutron densities respectively, k α is the nucleon Fermi momentum, T is the temperature, f α (μ * α , T ) andf α (μ * α , T ) are the thermal Fermi distribution function of the nucleon and the anti-nucleon, and μ * α is the effective chemical potential of the nucleon. The familiar Fermi distribution function for the particle and the anti-particle can be written as [31] where k B is the Boltzmann constant and E * α is the effective energy of the nucleon which can be written as The effective chemical potential of the proton and the neutron can be derived as [31] where μ p and μ n are the usual the chemical potential of proton and neutron relative to free nucleon mass. Another important thing we would like to mention for the formalism is the contribution of loop correction for the renormalized theory. Recently, Prakash et al. [33] reported that the twoloop approximation improves the energy density functional by adding density-dependent contributions to the Hartree terms of mean field theory from the exchange of isoscalar, isovector and pseudoscalar mesons. It is found that the thermal properties of proto-neutron stars differ significantly if we include the two-loop correction in the mean-field theory . However, it has also been reported in the references ( [34][35][36]) that the loop contribution in the mean field theory can be mocked up by refitting the parameters using experimentally available data. We calculate thermal properties with the mean-field approach without loop correction because in our work the coupling constants of the Lagrangian are obtained by fitting the experimental and observational data.

Theoretical formalism
In the following section, the formalism used for the calculations of symmetric and asymmetric NM properties has been discussed. To describe the asymmetric NM, we introduce the asymmetric parameter t, which is defined as, t = n n − n p n n + n p .
We can assign the desired asymmetry in the NM by varying the value of t. For symmetric NM (SNM), t = 0 and t = 1 for pure neutron matter. In the present work, we extend our calculations for SNM (t = 0). The energy density and the pressure for the NM can be calculated from the Lagrangian using the expression for energy-momentum tensor [37], which is where φ j includes all the fields present in the Lagrangian. Using this expression of the energy density and pressure for a warm nuclear system can be naively derived as [27] Since we deal with the temperature dependent NM, we also did some study regarding phase coexistence, which unfolds the liquid-gas phase transition in thermal NM. The NM system can remain only in one phase (gaseous) once it reaches the critical point [31]. For symmetric NM, the inflection point of the pressure curve with respect to the total nucleon density determines the critical point [32], that is T C being the critical temperature. One of the basic and fundamental quantity of the NM is incompressibility, also known as the isoscalar incompressibility or compression modulus (K ) [8]. K directly influences the curvature of the equation of state and gives adequate information about the nature of the equation of state. Higher the value of K, more stiff the equation of state will be [38]. The incompressibility is related to the equation of state through [37] K (n, Another important quantity which controls the equation of state of NM is symmetry energy. Symmetry energy has a significant contribution to the pressure of the astrophysical objects which is responsible for gravitational attraction [39]. In case of the hot NM, we can define two forms of symmetry energy, one we call as the nuclear symmetry energy (E sym ) and the other as free symmetry energy (F sym ). In order to calculate the free symmetry energy of the hot NM, we need entropy density S and free energy density F of that system. So the free energy density is given by and the entropy density for a hot NM is calculated by [37] Various theoretical studies have shown that we can calculate the temperature and density dependence of nuclear symmetry energy in hot symmetric NM using empirical parabolic approximation, which can be estimated as the difference of the energy per nucleon of pure neutron matter and symmetric NM [38]. In the similar way, we can find the free symmetric energy for a hot NM, which is of utmost importance for astrophysical phenomena. The free symmetric energy of a symmetric NM can be defined as [18,38] which can be interpreted as the remnant of free energy per nucleon for pure neutron matter and symmetric NM. The free symmetry energy can be expanded in a Taylor series expansion around χ . It is a dimensionless variable which describes the density deviation from saturation density and for mathematical convention and simplification of the expressions of higher order derivatives defined as, χ = (n − n 0 )/3n 0 , n 0 being the saturation density. So the expansion can be written as, where L sym , K sym and Q sym are the slope parameter, curvature parameter and skewness parameter respectively.
Recently it is convinced that these parameters are of utmost importance in nuclear and astrophysics by showing the correlation of these parameters with different nuclear and astrophysical properties [8]. The expressions for these parameters can be extracted as [40] L sym (n, Although we extract the results from this expansion at high densities, but we should not exclude the fact that at very high nuclear densities the error bar in higher order coefficients extracted through series expansion is quite high [8]. Thermal index commit an important aspect on core-collapse supernovae simulations [41]. The thermal index can be calculated with the help of simple formula discussed in [42] Γ th = 1 + E th P th (31) where E th and P th are the thermal energy density and pressure respectively. E th and P th are given by and P(T ) are the energy density and pressure at temperature T.

Results
Throughout our calculations in this paper, we use NL3, G3 and IU-FSU parameter sets, where NL3 endue us the stiff EoS and the other two (G3 and IU-FSU) facilitate us to examine the softer region. The obtained results for the variation in the properties of the NM with temperature and density are discussed in the present section. The coupling constants and the empirical values of nuclear properties at saturation of cold EoS for the assumed parameter sets (NL3, G3 and IU-FSU) are given in Table 1 [43,44]. We start our discussion with the energy and pressure of hot NM and procurement of critical temperature for liquid-gas phase transition. In Fig. 1, we depicted the calculated results of binding energy per nucleon and pressure density from T = 0 to 20 MeV as a function of density. We observed an increase in the binding energy and saturation density of the SNM with temperature for all the three supposed parameter sets. The marked point on each curve of the upper panel of the Fig.1 represents its minima. It has been already reported that NM saturates around 0.148 fm −3 with binding energy per nucleon around −16 MeV for both the parameter sets (NL3 and G3) at zero temperature [30] while the saturation density for IU-FSU parameter set is 0.154 fm −3 . Here we found that as we increase the temperature from 0 to 20 MeV, n 0 increases linearly from 0.148 to 0.196 fm −3 and similarly the binding energy also goes on increasing which is clear from the Fig.1. The increase in the binding energy per nucleon with temperature indicates that the system becomes more loosely bound at higher temperature. Lower panel of Fig.1 shows the variation of pressure with temperature as a function of density. This adaptation of pressure is utterly important in determining the critical parameters of liquid-gas phase transition specially the critical temperature, T C . Various theoretical and experimental studies predict the value of T C for SNM in the range of 10 − 20 MeV [32,[50][51][52][53]. We also found the value of T C using Eq. 22 as 14.60, 15.37 and 14.50 MeV for NL3, G3 and IU-FSU parameter sets respectively. Other important characteristics quantities of the liquid-gas phase transition like P C , ρ C can also be determined using T C with the help of correlations derived in reference [32].
The variation of the incompressibility modulus has also been depicted in Figs. 2 and 3. As the value of n increases, the chart of K exhibits an anomalous behaviour. We observed that the magnitude of K shows a maxima around 4 − 5 times of the saturation density for all the three parameter sets. However, the magnitude of K is quite high for NL3 parameter set in comparison to G3 and IU-FSU, which indicates a strong dependence of K on the nature of the EoS. We also observed an unfamiliar proneness in the K values at saturation density (denoted as K 0 ) for different temperature which are shown in Fig. 3.For NL3 parameter set K 0 shows an increment with increase in temperature while for G3 and IU-FSU parameter sets the case is reversed. This contemplation of K 0 confirms the sensitiveness of incompressibility on the choice of parametrization and the cross-coupling of field variables responsible for the softening of EoS. This fact has also been supported by the study in reference [54].
We also analysed that the G3 and IUFSU parameter sets indulge all the constraints of the elliptic flow heavy-ion collision experiment [55] on the value of K 0 , which is K 0 = 220±40 MeV. NL3, on the other hand, being the stiff parameter set depicts the higher value of K 0 and does not support Table 1 The coupling constants and the NM properties at saturation for cold EoS of NL3 [43], G3 [44] and IU-FSU [42] parameter sets. The nucleon mass (M) is 939.0 MeV. All of the coupling parameters are dimensionless and the NM parameters are in MeV, except k 3 and n 0 which are in fm −1 and fm −3 respectively. The NM parameters are given at saturation point and T = 0 K for NL3, G3 and IU-FSU parameter sets in the lower panel.  metry energy at saturation density (F sym,0 ) increases as we increase the temperature. Another important aspect that we observed here is the difference in the magnitude of the F sym for NL3 and G3 parameter sets. The magnitude range of F sym for NL3 is quite high in comparison to G3 and IU-FSU, which indicates that the stiff EoS allocate higher value of symmetry energy. Again the determination of higher value of F sym,0 along with K 0 by NL3 forces suggested that we need a better parameterized EoS, so that we can match the prediction of nuclear parameters made by various theoretical and experimental studies. On the other hand, the spectrum of F sym at finite temperature determined by our recently developed G3 parameter and the popular IU-FSU, fall in the expected range which is verified by other studies also [42,[56][57][58][59]. L sym changes slightly as we change the temperature for NL3 parameter set, while for G3 and IU-FSU it remains unaffected. The change in slope parameter at saturation density (L sym,0 ) with tem- perature is shown in Fig. 5. The effect of temperature on the second and third derivatives of F sym , which are known as curvature and isovector skewness parameter respectively (denoted by K sym and Q sym ) is shown in Figs. 6 and 7 for both the parameter sets. The variation of K sym and Q sym with density somehow reflects the sight of the sinusoidal wave, with extremum of the curve near twice the saturation density as can be seen in Fig. 6. K sym , being the higher order parameter, can be constrained with the help of neutron star observations [49]. Recently, a study done by Josef Zimmerman et al. merged the data reported by two most vital experiments (PSR J0030+0451 by NICER Collaboration [60] and GW170817 by LIGO/Virgo [61,62]) and derive a joint 1-σ constraint on curvature parameter at saturation density (K sym,0 ). It is believed to be the most reliable bound on K sym,0 till date and reported within the 90 % confident bounds as K sym,0 = −102 +71 −72 MeV [49]. In Fig. 7 we can see that the value of K sym,0 estimated by G3 parameter set for entire temperature spectrum falls into the range of 1-σ constraint discussed above. NL3 and IU-FSU parameter sets, on the other hand, not only inconsistent with the estimated constraint but surprisingly falls entirely in the opposite magnitude range. This observation clearly reject the NL3 and IU-FSU phenomenological model in determining the NM parameters for astrophysical observations and hints for G3 EoS, best suited for the study of astrophysical phenomena. Studies show that there is a strong correlation between K sym,0 , tidal deformability and R 1.4 (radius of NS with mass 1.4 M ) [63][64][65]. By obligating K sym,0 in the befitting range we can measure tidal deformability theoretically, which is quite a complicated quantity to measure independently, precisely upto a certain level of accuracy. Q sym being the most ambiguous quantity also shows entirely different behaviour for the NL3 and G3 parameter set. For G3 it falls in the positive magnitude range while for NL3 and IU-FSU it beholds the negative magnitude as the temperature increases. Although there is no reported constraint on Q sym,0 till now, but we believe that since G3 parameter set satisfy all the desired results for other NM parameters, so it is obvious to assume it of the right kind. However, some predictions made on the basis of skyrme interactions supports the negative magnitude of skewness parameter [66], so it is difficult to state anything about Q sym with certainty. Figure 8 enounce the variation of thermal index for the avowed parameter sets. The nature of the EoS of the astrophysical phenomenon is also dictated in some manner by the density dependence of thermal index. As we can see in Fig.  8 that the thermal index in the low density region approaches the value of 5/3 for all temperatures, which indicates that the system described by our model for all the parameter sets almost behaves like a non-relativistic ideal gas [42]. We also observe that as the temperature goes on increasing the max- imum value of thermal index decrease which is consistent with the fact that at high temperature, the system approaches relativistic behaviour. The value of Γ th around 2 has been appraised as an important aspect in the neutron-star merger simulations, which has also been reflected in our calculations [67]. Thermal index curve for NL3 is comparatively more steep than G3 and IU-FSU parameter sets and approaches towards the relativistic gas behaviour more rapidly. The variations in the thermal index with density for different EoS solidify the frequency of gravitational wave oscillations and also decides the delay time of black hole collapse [42,67].

Cooling through neutrino emission
In this section, we provide an approach to study the thermal evolution of the newly born neutron star through EOS. It is widely believed that the dense neutron star formed just after the supernova explosion contains tantamount nucleons i.e. almost equal number of protons and neutrons [68]. After that, the process of cooling and neutronization [68,69] happens slowly and finally it achieves the thermal stabilization and the beta equilibrium, which we call cold neutron star. The direct Urca processes n −→ p + e − +ν e and p + e − −→ n + ν e , governs the cooling of newly born neutron stars [3] through neutrino emission and requires a super threshold proton fraction inside the super-dense star core to operate. To realize the effect of neutrino emissivity on cooling mechanism and EOS, we consider a degenerate dense matter containing nucleons (neutrons and protons) and the electrons at finite temperature, for which the Lagrangian is redefined as where L is the Lagrangian defined in Eq. (1), m e being the the mass of electron and the last term represents the electrons contribution in the matter, which are considered to be non-interacting particles. Following the same procedure as defined in Sect. 2, we will get the total energy density and pressure for this system (33) (34) where E and P is the energy density and pressure defined in Eq. 20, k e is the Fermi momentum of electron, f e is the Fermi distribution of electron for finite temperature and E e is the energy of electron given by We maintained the same number density for proton and electron i.e. n p = n e , to achieve charge neutrality in the described system. The expression for the neutrino emissivity (Q) in high dense (neutron star) system was first estimated by Lattimer et al. in non-relativistic manner [70]. Since the movement of nucleons in neutron star cores is relativistic, so the relativistic expression for neutrino emissivity was later calculated by L. B. Leinson and A. Perez. The non-relativistic emissivity is quite small than what is predicted by the relativistic formalism. The detailed explanations and the derivation for the formula of neutrino emissivity (Q) in relativistic framework and mean-field approximation, which is used here can be found in the reference [71,72]. The formula for neutrino emissivity is given by [72] where G F = 1.166 × 10 −11 MeV −2 is the Fermi weak coupling constant, C = 0.973 is the Cabibbo factor, C V = 1 and C A = 1.26 are the vector and axial-vector constants respectively and the constant C M = 3.7 represents the weak magnetism effects. The condition inevitable for Urca processes to go is represented by Θ(k e + k p − k n ) = 1, if k e + k p − k n ≥ 0 and zero otherwise, where k e , k p and k n are the Fermi momenta of electron, proton and neutron respectively. Since, our defined system of nucleons and leptons (n-p-e) clearly satisfy the above necessary condition required for the initialisation of direct Urca process, so, we calculated the neutrino emissivity for both NL3, G3 and IU-FSU parameter sets which is displayed in Fig. 9. The dashed lines in the Fig. 9 represent the results for G3, the solid lines stand for NL3 and the dotted lines for IU-FSU parameter set. Some very interesting remarks can be concluded based on these curves of Q depicted at different temperatures. The most vital observation of this work is that the neutrino emissivity is responsible for cooling of newly born neutron star only in its initial stage i.e. when the temperature is quite high and the neutrons are enough thermally excited to trigger the direct Urca process. As the star cools down, the magnitude of Q decreases substantially and then cooling mainly takes place through photon emission. This behaviour can be seen clearly in Fig. 9 for all the assumed parameter sets. The magnitude of Q for 40 MeV (represented by black line) temperature is negligible in comparison to 80 MeV curve (blue line). Also, initially the diffusion rate of neutrinos is so high that the matter cools down within a fraction of seconds and it is believed that it is so effective that it lowers the temperature to about 1 billion Kelvin during this momentary span [24,73]. Another important dimension that we observe in this curve is the difference in the magnitude of the Q values for NL3, G3 and IU-FSU parameter sets. This tells us about the dependence of cooling property of newly born NS on the EoS. As we know that NL3 provides the stiffest EoS, which means the high mass of the proto-neutron star and G3 being the softer EoS in comparison to NL3 predicts lower mass. Keeping this fact in mind we concluded that the cooling of proto-neutron star through direct Urca process is slow in heavier star and fast in the lighter one.

Proto-neutron star
To understand the complete penchant of temperature on a dense matter system, we extend our calculations to obtain the mass-radius (M-R) profile of proto-neutron star (PNS) in this section. The newly born proto-neutron star is hot, leptonrich and has a core in evolutionary stage which ultimately grow either in a cold neutron star or crash into a black hole. As we discussed in the previous section that a PNS losses a significant amount of energy through neutrino emission and neutrino emissivity plays an important role in the evolution of the newly born star. However, some of the neutrinos with small mean free path unable to escape the core of the neutron star and got trapped inside the star. These trapped neutrinos affects the early evolutionary illation and properties of the PNS [74][75][76]. We modify our Lagrangian mentioned in Eq. (32) and added the terms responsible for neutrino effects to derive the EoS for the proto-neutron star system envisaged in the present section. The modified Lagrangian is where φ ν e is the wave-function of electron neutrino. The relevant energy density and pressure for the above mentioned Lagrangian can be derived as [75] where the summation is over the total number of trapped neutrinos; E total and P total are the energy density and pressure defined in eq. 33 and μ ν e is the chemical potential of electron neutrino. We also maintained the necessary β− equilibrium and charge neutrality conditions i.e. μ n = μ p + (μ e − μ ν e ), where μ n , μ p , μ e and μ ν e are the chemical potentials of the neutron, proton, electron and neutrino respectively; n p and n e are the number densities of proton and electron. We can also define the lepton fraction as Y L = n e + n ν e n (41) where n e , n ν e and n are the number density of electrons, neutrinos and baryons respectively. The effect of trapped neutrinos on the early evolutionary stages of PNSs can be represented by fixing Y L in a certain defined range i.e. Y L ∼ 0.1−0.4. We fix Y L = 0.4 in our calculations, which is the best precedence to explore the properties of newly born PNS [75][76][77]. Also, the entropy per baryon (S) of the star matter can be derived as [77] S = E P N S + P P N S − α= p,n n α μ α nT (42) For a fixed entropy, the effect of the neutrino trapping is to keep the electrons concentration high so that matter is more proton-rich in comparison to the case in which the neutrinos are not trapped [16]. We study the M-R profile for the PNS using different kind of EoS, first by placing the temperature constant throughout the star (fixed temperature EoS) [78][79][80] and another by fixing the certain defined values of entropy per baryon (S = 1&2) for the star matter (fixed entropy EoS) [81][82][83][84]. Although we have used the constant temperature EoS to describe the mass-radius profile of proto-neutron star, but this assumption is valid only for low temperature spectrum. If the matter's temperature is significantly larger than the corresponding critical Fermi temperature, the star becomes unstable. Also, the temperature of the star varies significantly from core to the surface of the star. So, to debar all these restrictions and for a better undisputed M-R outline of PNS, we extended our calculations for fixed entropy EoS. The method of fixed entropy EoS is more appropriate to study the properties of PNS at finite temperature. By imposing the above stated conditions on the described Lagrangian, we can easily calculate the mass and radius of the static isotropic proto-neutron star for both kind of EoSs using the Tolman-Oppenheimer-Volkov (TOV) equations [85,86]. The TOV equations are Where E P N S and P P N S are the energy density and pressure defined in Eqs. (38) and ( produces the maximum mass M max for a given EoS is ρ max . The crust part of the neutron star also plays an important role in determining the properties of the proto-neutron star [87]. We have also used the crust EoS of the corresponding temperature for a complete contact of Mass-Radius profile [88]. The mass-radius profile of the proto-neutron star for fixed temperature EoS along with trapped neutrinos is plotted in Fig.10. The results from the precisely measured neutron stars masses, such as PSR J1614-2230 with mass M = 1.97 ± 0.04M [89] and PSR J0740+6620 with M = 2.15 +0. 10 −0.09 M [90] are shown in the horizontal bars in pink color. These observations suggest that the maximum mass predicted by any theoretical model should reach the limit ∼ 2.0M , and this condition is satisfied by all of the EoSs taken into consideration. Although G3 parameter set follow the observational constraint of PSR J0740+6620 for the entire assumptive temperature range. We observe that the inclusion of neutrino trapping flatten the M-R curve at the top for G3 and IUFSU parameter sets, which affect the radius of the proto-neutron star considerably. We also notice that the inclusion of temperature increases the pressure at a given baryon density which yields the increase in the mass-radius profile of the protoneutron star. The proto-neutron star has a little-bit large mass compared with that of the neutron star at zero temperature because the EoS is stiffer in the former case. The mass-radius profile of proto-neutron star using constant entropy EoS have also been studied and depicted in Fig. 11. By analysing the Figs. 10 and 11, we concluded that the maximum mass determined by the G3 parameter set strongly satisfy the constraints of the PSR J0740+6620 and PSR J1614-2230. We observe an increase in maximum mass and radius in comparison to the cold NS for both kind of EoS (fixed T and S). The maximum mass, radius and central tem-  Table 2. The maximum mass of the PNS for all the three parameter sets goes on increasing with increase in entropy of the system. Also, we observe that the massive PNSs have larger radius, so we concluded that more the entropy per baryon of the star, larger the radius of the system. In other words, we can say that the evolutionary process of the newly born star favours contraction in terms of size. It has also been reported that the maximum mass of a proto-neutron star is considerably affected by the earliest stage of its evolution i.e. ratheripe type proto-neutron star (∼ 1 sec after core bounce) have larger mass than the late type proto-neutron star [84]. So, the properties of the PNS depends strongly on how immaculately we gravitate the entropy per baryon of the system in its early stage. Also, to reconcile the results for maximum mass and radius determined by the two different approaches (i.e. constant entropy and constant temperature), we calculated the mass and radius of the PNS using constant temperature approach by fixing the temperature at T = T Ce for the corresponding parameter set. The results for the same have been presented in the Table 3. We observe that the fixed temperature approach provide larger mass and radius at T = T Ce of the corresponding entropy value. G3 parameter set demonstrate ∼ 2.5% and IU-FSU set shows ∼ 4.5% difference in the radius of the PNS calculated using constant entropy and temperature approach for the corresponding T = T Ce . At last, we concluded that the properties of a PNS can be decipher more appropriately by using the constant entropy perspective and as discussed in the previous sections also, since G3 parameter set respects all the experimental and observational constraints thoroughly, so it is one of the most compatible RMF parameter set.

Summary and conclusions
We have studied the consequences of finite temperature on the nuclear properties of SNM and the crucial section required for the cooling of remnants of supernovae explosion. We used well known NL3 and IU-FSU and the recently developed G3 parameter sets of the RMF model for a comparative study and concluded that NL3 with the possession of stiffest EoS does not provide the empirical values for most of the nuclear properties. However, IU-FSU satisfy most of the constraints on the NM properties but also neglect some of them. The variations of the binding energy and pressure with baryon density for different temperatures are qualitatively similar for all the three acquired parameter sets. However, G3 predicts the higher value of critical temperature for liquidgas phase transition, which is more proximate to the reported experimental values. We also observed a contrarious development in the value of K 0 for the defined parameter sets with increase of temperature. Both K and K sym being the second derivatives of the different forms of energy, this behaviour of K 0 may be influenced by the magnitude of K sym which is positive for NL3 and negative for the G3 parameter set. K sym value of the G3 parameter set lies in the range reported by the NICER and LIGO collaboration. This result shows that the G3 parameter set is more suitable to reproduce the appropriate form of the EoS and can be used to study the properties of neutron stars more accurately. We did not observe any significant variation in the Q sym parameter at saturation density for the G3 parameter set, while its value decreases for the NL3 and IU-FSU parameter sets. We emphasized on the cooling mechanism of high dense matter and studied the effects of all three parameter sets on the neutrino emissivity. We used the relativistic approach to derive the detailed expression of neutrino emissivity, which is more effective than the non-relativistic approach. We concluded some important remarks about the cooling mechanism of the newly born proto-neutron star on the basis of the outcomes of our detailed calculation. We observed that the magnitude of neutrino emission is directly proportional to the temperature of the neutron star. As the body cools down, the magnitude of the neutrino emission decreases. Also, the neutrino emissivity has maximum value around the saturation density, which indicates that the saturated matter cools more rapidly through direct Urca process. Another important aspect that we inspect is that the neutrino emissivity is higher for that parameter set which provides softer EoS. As we can see in the previous section that the stiff EoS measures higher mass NS, so we concluded that the lighter remnant of the supernovae explosion cools down more expeditiously through neutrino emissivity of the direct Urca process. Moreover, We look forward to see if any post-merger signal is explored observationally in near future by LIGO /VIRGO and NICER collaboration, which will help us to comprehend finite temperature NM physics more appropriately and extensively. We also take notice that the mass numerated by the G3 parameter set for whole temperature range fits in the envelope determined by the GW170817 experimental data. Finally, we adduced that G3 parameter set is more appropriate to study the properties of stellar objects through EoS.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .