## Investigating the role of nuclear parameters on oscillation modes in hot Neutron Stars

Nilaksha Barman , ∗ Bikram Keshari Pradhan , † and Debarati Chatterjee ‡

Inter-University Centre for Astronomy and Astrophysics, Pune University Campus, Pune 411007, India

Recent studies have revealed that certain nuclear parameters are more dominant than others in governing global neutron star properties, such as its structure or oscillation mode characteristics. Although neutron stars can in general assumed to be cold, in astrophysical scenarios such as newly born neutron stars or remnants of binary neutron star mergers, finite temperature effects play a non-negligible role. In this work, we perform a consistent and systematic investigation of the role of nuclear parameters and thermal effects on neutron star properties and fluid oscillation modes within a full general relativistic scheme. We impose constraints on the parameter space of the relativistic mean field model using state-of-the-art information from terrestrial experiments and multi-messenger astrophysical data. We find effective nucleon mass to be the most important nuclear parameter controlling astrophysical observables of hot neutron stars, similar to the cold beta equilibrated matter. However, we conclude that the interplay among saturation properties and astrophysical observables depends not only on the thermal configurations considered but also on the constraints imposed. We also investigated the role of nuclear saturation parameters on some universal relations for hot NSs which are important in gravitational wave asteroseismology. Our investigation confirmed that these relations are mostly insensitive to nuclear saturation properties and mainly affected by variation of charge fraction in the star.

## I. INTRODUCTION

Neutron stars (NSs) are the densest known compact objects in the observable universe. They provide us unique opportunity to study matter under extreme conditions and hence sometimes referred to as cosmic laboratories. They are formed at the end stage of stellar evolution of massive stars known as core-collapse supernovae or type-II supernovae. Hot newly born neutron stars undergo cooling via neutrino emission and mature into cold degenerate compact objects. Cold NSs are primarily composed of neutrons, protons, electrons and muons (free neutrinos escape at low temperature &lt; 1 MeV). However, in the high density environment of the NS core other exotic particle degrees of freedom such as hyperons, quarks, kaon condensates etc. might appear in addition to nucleons via weak interactions and can be in equilibrium with the mixture [1-5].

saturation properties. Chiral effective field theory ( χ EFT) [13] is a microscopic method which uses chiral expansion for a reliable description of the sub-saturation region of pure neutron matter. χ EFT results are being used widely in phenomenological or parametric models as a robust method to constrain the parameter space of dense matter EoS at low densities e.g. [14-17].

The composition of the cold NS interior can be described by Equation of State (EoS) models that relate the pressure and the density [6]. The EoS connects the microscopic properties of the NS to its global macroscopic structure, such as its mass or its radius. However, the ultradense nuclear matter in the NS core requires solving the complicated many-body problem. There are different approaches e.g. microscopic models [7-10] which calculate the nucleon-nucleon interactions abinitio, while phenomenological models [11, 12] adjust the nucleonic parameters to satisfy observed nuclear

∗ nilaksha.barman@iucaa.in

† bikramp@iucaa.in

‡ debarati@iucaa.in

Being highly relativistic and compact, neutron stars can emit gravitational waves (GWs) via both binary neutron star (BNS) mergers as well as in isolated systems. Any non-axisymmetric perturbations, e.g. quasi-normal modes, deformations such as mountains or magnetic deformations, can lead to time varying mass quadrupole moment in isolated neutron stars and may manifest themselves as gravitational waves. The frequencies of some of the oscillation modes lie in the high frequency band of GWs (e.g. f -mode frequencies lie in 1-2 KHz) and provide interesting prospects of detection by current and future generation GW detectors. The first detection of BNS (GW170817) via GW and electromagnetic (EM) signals opened up a new window in the multimessenger astronomy [18-20]. In this context previous works have shown that spin and eccentricity effects lead to excitation of f -modes during the inspiral phase of BNS mergers [21, 22]. With LIGO O4 run and subsequent future GW observations many more BNS merger events are expected to be detected.

Any realistic model of neutron stars should satisfy astrophysical constraints from EM and GW observations. The maximum predicted theoretical NS mass should be able to explain the highest observed NS mass of 2 08 . +0 07 . -0 07 . M ⊙ [23] or 2 078 . +0 067 . -0 066 . M ⊙ [24]. EM observations from the recently launched NICER put constraints on masses and radii two pulsars namely PSR J0030-0451

[25, 26] and PSR J0740-6620 [23, 24]. Analysis of GW170817 revealed that the radius of a 1.4M ⊙ neutron star is 13 4 . ± 0 2 km [27] and the estimates for effective . tidal deformability for the binary system is 300 +420 -230 [20]. These multi-messenger observations are being used to constrain the high and intermediate density region of nuclear EoSs [14-17, 28].

Astrophysical events such as proto-neutron stars (PNSs) [29-32] born out of supernovae (SN) and binary neutron star (BNS) mergers [33-38] exceed temperatures of tens of MeV and the cold matter approximations are no longer valid. It is important to take into account thermal effects as they play a significant role in supernova dynamics, oscillations in PNS as well as in determining stability of the post-merger remnants, which along with the composition can change the threshold mass limit for prompt collapse to a black hole [39]. The EoS of a finite temperature system can be uniquely characterized by three quantities, temperature T , baryon number density n b and the charge fraction of strongly interacting particles Y Q . In order to be applicable in numerical simulations, a wide range and finely spaced sets of parameters should be considered, along with different possible composition in the neutron star interior. The EoS of hot homogeneous matter spans 0-100 MeV in temperature, 0.01-0.60 in charge fraction and 0.5-6.0 times in nuclear saturation densities in BNS mergers and supernovae [6]. In the absence of EoS models spanning the required parameter ranges and composition, the so-called Γ-law (supplementing the cold EoS with ideal gas thermal contributions) has often been employed to incorporate finite temperature effects in simulations. However, this method of constructing EoS has been shown to be inappropriate particularly for matter consisting of non-nucleonic degrees of freedom [40-43]. It has been described in a recent work in quark-hadron phase transition [43] that Γ value should be different in hadronic (Γ ∼ 1 75) and quark (Γ . ∼ 4 3) / phases. The Γ-law also does not take into account the effect of charge fraction across all densities. Therefore it is crucial to incorporate composition and thermal EoS effects consistently into simulations.

Currently, there exist only a limited number of available finite temperature EoSs which are widely applied in simulations e.g. Lattimer Swesty (LS220) [44] and Shen et al. (STOS) [45], both for nucleonic matter and based on single nucleus approximation (SNA). The more recent EoS model of Hempel and Schaffner-Bielich (HS) [46] is based on nuclear statistical equilibrium (NSE) which incorporates realistic nuclear masses and goes beyond single-nucleus approximation. The high density core can be modelled using the RMF formalism and is used in constructing non-linear (NL) RMF EoSs such as SFHo/SFHx [47], HS(IUF) [48], QMC-RMF [49]. The RMF EoSs based on linear density dependent (DD) models use the formalism of Typel and Wolter [12].

They include nucleonic EoSs such as HS(DD2) [46] as well as models with exotic matter such as hyperons and ∆-baryons e.g. BHB(DD2L ) [50], OMHN(DD2Y) [51], ϕ R(DD2Y∆) [40]. Other finite temperature EoSs include hybrid Chiral Mean field (CMF) model [52, 53].

The existing finite temperature EoS models, which consider a fixed parametrization, correspond to certain values of nuclear saturation properties but they do not span the entire range of uncertainty allowed by current terrestrial experiments. However recent studies have shown that nuclear saturation properties are correlated with astrophysical observables such as radius and tidal deformability in cold neutron stars [14, 15, 17, 54-56]. These studies provide an interesting prospect to investigate the role of nuclear saturation effects on NS observable properties at finite temperature. There have been a few recent works in this direction. An EoS formalism within the covariant density functional approach with density dependent couplings employing the posteriors of a Bayesian inference (DDB*) formalism [17] has been extended to finite temperature by Raduta et al. [57]. In another approach, the effect of slope of symmetry energy L sym and skewness Q sat on finite temperature EoS has been studied within the covariant density functional scheme for fixed DD-RMF parametrizations in [58, 59].

Significant efforts have been dedicated for advancing our understanding of GW asteroseismology for PNSs, through the study of GWs emitted by PNS oscillation modes [60-62]. Both simulations [62-66] and perturbative approaches [56, 67] have been considered in studying the PNS oscillations. In the simulation by Torres et al. [62], it was shown that most of the GW energy is stored in the lowest-order eigenmodes, in particular 2 g 1 and 2 f . The same group in their previous study in Cowling approximation showed that higher order p i -modes upto 5 nodes and g 1 -modes appear in GW spectrogram [68]. PNS simulations by Sotani et al. [64, 67, 69] also indicate that g 1 and f -modes appear to dominate GW spectra in the earlier and later stages of PNS evolution respectively. Further analyses involve assessing the detectability of GWs from PNSs and inferring PNS properties from the GW data [70, 71].

Cold NSs are known to follow certain EoS independent universality such as C -Love, ¯ I -Λ -¯ , Q f -Love [54, 7276] among others. Much effort has been dedicated into deriving universal relations (URs) that link the oscillation mode properties of PNSs with their physical characteristics, such as mass and radius, which are crucial for GW asteroseismology [61, 63, 66, 77] of protoneutron stars. Raduta et al. [78] obtained ¯ I -Λ -¯ Q -C relations for hot NSs using a DD-RMF parameterization. This was extended by Khadkikar et al. [79] to explore URs between static and Kepler rotation properties such as mass, radius and frequency etc. and fit relations were

provided for ¯ I -C and ¯ Q -C for different thermal configurations with various tabulated EoSs. Similar UR study for nucleonic and hybrid hot stars were performed by Largani et al. [80] where they provided fits for normalized moment of inertia with compactness at Keplerian limit for fixed values of entropy per baryon and lepton fraction. Laskos-Patkos et al. [81] studied deviations of hot neutron stars from cold universality for C -Love relations.

In a series of recent works [14-16], we developed a formalism to constrain the EoS of cold dense NS matter using χ EFT at low densities, heavy-ion data at intermediate densities and multi-messenger astrophysical data at high densities, within a non-linear relativistic mean field model (NL-RMF) framework. These works also investigated the correlation of nuclear empirical parameters with NS macroscopic properties (radius R , tidal deformability Λ) and found the effective nucleon mass m ∗ to be strongly correlated to these observables. Further, recent investigations within the same formalism also showed that m ∗ is the dominant parameter related to oscillation modes such as fundamental ( f -) mode in NSs [54, 82, 83]. However, it is not clear whether such correlations will hold in finite temperature scenarios, such as oscillations in newly born neutron stars or binary neutron star mergers. A few interesting studies recently probed thermal effects on the NS EoS imposing constraints from χ EFT and NS maximum mass within the NL-RMF formalism (see e.g. [49]). In another recent work [56], thermal effects were studied by considering a few density-dependent (DD) RMF tabulated EoSs from CompOSE [84]. The effect of finite temperature on f -and first overtone pressure ( p 1 )-modes was studied within Cowling approximation for the chosen EoSs. However, the correlations among nuclear and NS properties were investigated only for cold NS matter employing the DDB* model. From these studies, the correlations of the nuclear saturation properties on NS global properties and general relativistic mode characteristics at finite temperature are not directly evident.

The aim of this work is to perform a systematic and consistent investigation of the role of nuclear saturation parameters as well as thermal properties i.e., entropy per baryon, S/A and out-ofβ equilibrium condition (charge/electron fraction, Y Q ) on the EoS of hot and dense matter. For this we extend our recently developed formalism of imposing constraints at different densities on the parameter space in the RMF scheme for cold matter EoS to incorporate thermal effects. The finite temperature EoS is then used to derive macroscopic properties of hot neutron stars. The EoS framework is further applied for the study of oscillation modes within a fully general relativistic treatment. With the posteriors, we search for possible correlations between nuclear saturation parameters as well as thermal properties with astrophysical observables. We also investigate the effect of nuclear saturation and thermal properties on the Universal relations. Note that we have not included the effects of neutrino trapping in our study.

The paper is structured as follows: In section II we provide a description of the formalism involved. In subsection II A the RMF model and its extension to include thermal effects have been given. We also derive the expressions of various thermodynamic quantities of interest. In the next subsection II B we determine macroscopic NS properties, while in the subsection II C the fully General Relativistic (GR) description of calculating the oscillation modes are given. In section III we discuss results for thermal effects on EoSs consistent with χ EFT, astrophysical and heavy ion collision constraints as shown in subsection III A and their global properties in subsection III B. We then discuss results for oscillation mode frequencies and damping timescales in III C. We also present the results of the investigation of correlations of the nuclear properties with observable NS properties. In III D we show deviation of universality of oscillation modes at finite temperature compared to cold β equilibrated matter. Then we summarise our results in section IV. Throughout the paper we use the natural units ( c = ℏ = k B = G = 1).

## II. FORMALISM

## A. Equation of State

The EoS of the NS core is modelled using a non-linear relativistic mean field (NL-RMF) formalism [11, 85]. The nucleon interactions occur via exchange of σ , ω and ⃗ ρ mesons. σ is a scalar meson that facilitates long range attraction while ω is a vector meson for nucleon-nucleon repulsion. ⃗ ρ vector meson is responsible for isospin asymmetry and to maintain charge neutrality. The non-linear terms are the scalar self-interaction term U σ given by U σ = 1 3 bm N ( g σN σ ) 3 + 1 4 c g ( σN σ ) 4 , and Λ ω which mimics the ω -ρ interaction. The Lagrangian density for the homogenous nuclear matter is given below:

$$\text{is and} \quad & \homogenous nuclear matter is given below: \\ \text{ic and} \quad & \mathcal { L } = \sum _ { i \in B } \bar { \psi } _ { i } ( i \gamma _ { \mu } \partial ^ { \mu } - m _ { i } + g _ { \sigma i } \sigma - g _ { \omega i } \gamma ^ { \mu } \omega _ { \mu } - g _ { \rho i } \gamma ^ { \mu } \rho _ { \mu } ^ { \mu } \tau ^ { a } ) \psi _ { i } \\ \text{tropy} \\ \text{dition} \\ \text{t and} \quad & \quad + \frac { 1 } { 2 } \partial ^ { \mu } \sigma \partial _ { \mu } \sigma - \frac { 1 } { 2 } m _ { \sigma } ^ { 2 } \sigma ^ { 2 } - \frac { 1 } { 4 } \omega _ { \mu \nu } \omega ^ { \mu \nu } + \frac { 1 } { 2 } m _ { \omega } ^ { 2 } \omega ^ { 2 } \\ \text{loped} \\ \text{sities} \\ \quad & \quad - \frac { 1 } { 2 } \vec { \rho } _ { \mu \nu } \vec { \rho } ^ { \mu \nu } + \frac { 1 } { 2 } m _ { \rho } ^ { 2 } \vec { \rho } ^ { 2 } + \Lambda _ { \omega } ( g _ { \rho N } ^ { 2 } \vec { \rho } _ { \mu } \cdot \vec { \rho } ^ { \mu } ) ( g _ { \omega N } ^ { 2 } \omega _ { \mu } \omega ^ { \mu } ) \\ \text{r cold} \\ \text{finite} \\ \text{sopic} \\ \text{ework} \\ \text{over} \\ \quad & \quad \quad \quad$$

where, ψ i are the Dirac spinor for nucleon(s). The isoscalar couplings ( g σB , g ωB , b, c ) are fitted to reproduce isoscalar saturation properties (saturation density n sat , energy per particle E sat , compressibility K sat , effective nucleon mass m ∗ ) while isovector coupling constants

( g ρB , Λ ) are fitted ω to reproduce isovector saturation properties (symmetry energy J sym and its slope L sym ). The leptonic species present are non-interacting electrons ( e -) obeying Fermi-Dirac (FD) statistics; muons ( µ -) are not considered as they are expected to appear after the star reaches β equilibrium. Photons ( γ ) are bosonic non-interacting degrees of freedom described by A µ and F µν = ∂ µ A ν -∂ ν A µ . Under RMF approximation the meson mean field equations become,

$$m _ { \sigma } ^ { 2 } \bar { \sigma } = \sum _ { i \in B } g _ { \sigma i } n _ { i } ^ { s } - \frac { \partial U } { \partial \bar { \sigma } }$$

$$m _ { \omega } ^ { 2 } \bar { \omega } _ { 0 } = \sum _ { i \in B } g _ { \omega i } n _ { i } - 2 \Lambda _ { \omega } g _ { \rho } ^ { 2 } g _ { \omega } ^ { 2 } \bar { \rho } _ { 0 3 } ^ { 2 } \bar { \omega } _ { 0 } \quad \quad ( 3 ) \quad \text{when} \quad \dots$$

$$m _ { \rho } ^ { 2 } \bar { \rho } _ { 0 3 } = \sum _ { i \in B } g _ { \rho i } t _ { 3 i } n _ { i } - 2 \Lambda _ { \omega } g _ { \rho } ^ { 2 } g _ { \omega } ^ { 2 } \bar { \rho } _ { 0 3 } \bar { \omega } _ { 0 } ^ { 2 } \, \quad ( 4 ) \quad \text{thi} \quad \\ \text{stitu}$$

where ¯, σ ¯ , ω 0 ¯ ρ 03 are the mean meson fields and t 3 i is the third component of the isospin of the i -th baryonic species.

Including thermal effects, the number density and scalar density of the i -th baryonic species respectively are [40, 57]:

$$\intertext { a v e } \lfloor \omega, \, \vartheta \rfloor. & & \text{the c} \\ n _ { i } & = \frac { 2 J _ { i } + 1 } { 2 \pi ^ { 2 } } \int d k k ^ { 2 } \left [ f _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) & & \epsilon \\ - & f _ { F D } \left ( \frac { E _ { i } ( k ) + \mu _ { i } ^ { * } } { T } \right ) \right ] \,, & & ( 5 ) \\ n _ { i } ^ { s } & = \frac { 2 J _ { i } + 1 } { 2 \pi ^ { 2 } } \int d k \frac { k ^ { 2 } m _ { i } ^ { * } } { E _ { i } ( k ) } \left [ f _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) & & \\ + & f _ { F D } \left ( \frac { E _ { i } ( k ) + \mu _ { i } ^ { * } } { T } \right ) \right ] \,. & & ( 6 ) \\ \intertext { In a v e } \, \mathbb { s o n d } \, \mathcal { K } \, \mathbb { m } ^ { * } - \mathbb { m }. - \alpha \, \mathcal { \bar { \sigma } } \, \mathfrak { i t h o d o f o t w o n o n }$$

In eqs. 5 and 6, m ∗ i = m i -g σi ¯ is the effective baryon σ mass and E k i ( ) = √ k 2 + m ∗ i 2 , the effective energy dispersion of the i -th baryonic species. J i is the spin of the i -th baryonic species. f FD ( x ) = 1 (1 + / e x ) is the FD distribution function, T is the equilibrium temperature of the mixture and µ ∗ i is the effective chemical potential of the i -th species given by,

$$\mu _ { i } ^ { * } = \mu _ { i } - g _ { \omega i } \bar { \omega } _ { 0 } - g _ { \rho i } t _ { 3 i } \bar { \rho } _ { 0 3 } \. \quad \quad ( 7 ) \quad \text{em} \\ \text{task}$$

where µ i is the actual chemical potential of the corresponding species and can be written in terms of baryon ( µ B ) and charge ( µ Q ) chemical potentials. If B i and Q i are baryon and electric charges of the i -th species then, µ i = B µ i B + Q µ i Q . Both particle and anti-particle species contributes to entropy, as do leptons/anti-leptons and radiation. The effective chemical potentials for antiparticles are negative of corresponding particle effective chemical potentials. Chemical potential for electrons is µ e = µ L -µ Q , where µ L is the lepton chemical potential.

The total entropy density ( s ) is given by,

$$\begin{array} { c } \text{ } The total entropy density ( s ) \text{ is given by}, \\ \text{ } s = - \sum _ { i \in B } \frac { 2 J _ { i } + 1 } { 2 \pi ^ { 2 } } \int d k k ^ { 2 } \left [ f _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) \\ \text{ } \ln f _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) + \bar { f } _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) \\ \text{ } \ln \bar { f } _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) \right ] \\ \text{ } \text{ } \text{ } + a n t i - b a r y o n \, c o n t r i b u t i o n \, ( \mu _ { i } ^ { * } \to - \mu _ { i } ^ { * } ) \\ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text { } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
 \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{
 \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\}} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{	\text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ }
 \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
\text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ }; \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{]} \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{


 \end{array}$$

where, ¯ f FD ( x ) = 1 -f FD ( x ) is the Fermi-Dirac distribution function for vacancy of state. The electronic contribution to the entropy density ( s e -) is found by substituting electron chemical potential ( µ e ) in place of effective chemical potential ( µ ∗ i ) in eq. 8, while the radiation ( γ ) contribution is s γ = 4 π T 2 3 / 45. The entropy per baryon is then given by S/A = s/n B . Also radiation energy density, ϵ γ = π T 2 4 / 15. From the Lagrangian described in eq. 1, the energy density can be evaluated from the energy momentum tensor T µν in the mean field approximation, and the pressure can be evaluated using the Gibbs-Duhem relation, as given below:

$$\begin{array} { l } \text{the Gibbs-Duhem relation, as given below:} \\ \\ \epsilon = < T ^ { 0 0 } > \\ = \frac { 1 } { 2 } m _ { \sigma } ^ { 2 } \bar { \sigma } ^ { 2 } + \frac { 1 } { 2 } m _ { \omega } ^ { 2 } + \frac { 1 } { 2 } m _ { \rho } ^ { 2 } \bar { \rho } _ { 0 3 } ^ { 2 } + U _ { \sigma } \\ + 3 \Lambda _ { \omega } ( g _ { \rho N } g _ { \omega N } \bar { \rho } _ { 0 3 }. \bar { \omega } _ { 0 } ) ^ { 2 } \\ + \sum _ { i } \frac { 2 J _ { i } + 1 } { 2 \pi ^ { 2 } } \int d k k ^ { 2 } E _ { i } ( k ) \left [ f _ { F D } \left ( \frac { E _ { i } ( k ) - \mu _ { i } ^ { * } } { T } \right ) \\ + f _ { F D } \left ( \frac { E _ { i } ( k ) + \mu _ { i } ^ { * } } { T } \right ) \right ] \\ + \epsilon _ { e ^ { - } } + \epsilon _ { e ^ { + } } \\ i s _ { i } = 0 \\ p = \sum _ { i = B, e } \mu _ { i } n _ { i } + T s - \epsilon \. \quad & ( 1 0 ) \\ \text{re} \\ \text{ial} \\ \text{The electron contribution to the energy density is again} \\ \text{found by replacing effective chemical potential (} u ^ { * } ) \text{ with} \end{array}$$

The electron contribution to the energy density is again found by replacing effective chemical potential ( µ ∗ i ) with electron chemical potential ( µ e ). For describing the low density part of the EoS the HS(DD2) model is used [46], taken from CompOSE database [84]. For each thermal case, the crust and core are matched by keeping S/A and Y Q same for crust [HS(DD2)] and core (NL-RMF) ensuring thermodynamic consistency.

The direct consequences of finite temperature are effects such as out-ofβ equilibrium, neutrino trapping among others. These effects can manifest themselves in the hot EoS and consequently in the macroscopic properties. At postbounce evolution in PNSs, the thermal properties control the EoS until the neutron

star becomes cold and degenerate. These thermal properties are entropy per baryon ( S/A ) and electric charge (out-ofβ equilibrium) and/or lepton fractions (in case of neutrino trapping). During PNS evolution and also BNS merger scenarios, the system goes through different stages of S/A and Y Q . We discuss two such thermal configurations relevant for PNS evolution in section III and show the effect of variation of nuclear properties on these configurations subsequently. It has been shown in previous works that this leads to deviations from universal relations (URs) usually observed in cold neutron stars [51, 78, 86].

To verify expected behaviour of thermal properties as given in existing literature, we first use a fixed nuclear parameterization (from table I) and vary only the thermal properties of the star. The core temperature profiles corresponding to different thermodynamic conditions for nucleonic matter N S/A,Y ( Q ) i.e., different values of entropy per baryon S/A and charge fraction Y Q , are shown in fig. 1. It can be seen that higher entropy leads to larger temperatures and higher proton fraction results in lower temperatures at high densities.

FIG. 1: Temperature profiles in the NS core, for thermodynamic conditions N S/A ( = 1 , Y Q = 0 2), . N S/A ( = 1 , Y Q = 0 4), . N S/A ( = 2 , Y Q = 0 2) and . N S/A ( = 2 , Y Q = 0 4). . The saturation properties have been fixed to n sat = 0 15 fm . -3 , E sat = -16 0 MeV, . K sat = 240 MeV, J sym = 32 MeV, L sym = 60 MeV, m /m ∗ = 0 65. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-d0097c90-7eb3-4217-80d5-04d4d913bc66.jpg)

## B. Macroscopic and Tidal properties

We assume a perfect fluid configuration which is static and spherically symmetric. The space time geometry for such a configuration is given by Schwarzschild interior metric

$$d s ^ { 2 } = - e ^ { 2 \Phi ( r ) } d t ^ { 2 } + e ^ { 2 \lambda ( r ) } d r ^ { 2 } + r ^ { 2 } ( d \theta ^ { 2 } + s i n ^ { 2 } \theta d \phi ^ { 2 } ) \. \ ( 1 1 ) \quad \text{ample}$$

consisting of metric functions Φ( r ) and λ r ( ) which depend on the underlying equation of state [ p = p ϵ ( )]. The macroscopic configuration of the hot neutron star is determined by the hydrostatic equilibrium condition and is obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equations:

$$\frac { d m ( r ) } { d r } = 4 \pi r ^ { 2 } \epsilon ( r ) \, \quad \ \ \ \ \ \ \ \ \ \ ( 1 2 )$$

$$\frac { \overset { \omega } { d p } ( r ) } { d r } = - [ p ( r ) + \epsilon ( r ) ] \frac { [ m ( r ) + 4 \pi r ^ { 3 } p ( r ) ] } { \left [ r ^ { 2 } \left ( 1 - \frac { 2 m ( r ) } { r } \right ) \right ] } \, \quad ( 1 3 ) \\ \dots \, \quad 1$$

$$e ^ { 2 \lambda ( r ) } = \frac { 1 } { \left [ 1 - \frac { 2 m ( r ) } { r } \right ] } \, \quad \ \quad \ \ \ \ \ \ ( 1 4 ) \\ d \theta ( r ) \ \ 1 \quad d n ( r )$$

$$\frac { d \phi ( r ) } { d r } = - \frac { \mu } { [ \epsilon ( r ) + p ( r ) ] } \frac { d p ( r ) } { d r } \. \quad \quad ( 1 5 )$$

Here m r ( ) is the mass enclosed within a spherical shell of radius r and therefore if R is the radius of the star then m r ( = R ) = M , the NS mass. The tidal properties of the star are encoded in tidal love number k 2 [73] and it is related to dimensionless tidal deformability Λ through stellar compactness C = M/R .

$$\Lambda = \frac { 2 } { 3 } \frac { k _ { 2 } } { C ^ { 5 } } \.$$

We display the NS mass-radius relations in fig. 2 corresponding to the thermodynamic conditions mentioned in the previous section. The case N T ( = 0) stands for β equilibrated nucleonic matter.

FIG. 2: Mass-Radius relations for the different thermal configurations in fig. 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-3ec98f29-262d-44c7-9d1d-bac5d5cf6d3c.jpg)

It can be observed from fig. 2 that the increase of charge fraction leads to a decrease in maximum mass supported by the hot neutron star. Both the increase in S/A and Y Q results in higher radii for the same mass configurations in intermediate mass ranges. As an example we have also shown the isothermal configuration

with T = 20 MeV and Y Q = 0 2. . For this configuration low density part is taken again from HS(DD2) EoS with T = 20 MeV and Y Q = 0 2 from CompOSE [84] to ensure . thermodynamic consistency. This configuration has the flattest mass-radius relations among all the cases.

## C. Oscillation modes

In general relativity (GR), NS oscillation modes are analyzed by considering perturbations in both the fluid and the background metric. These modes are classified into polar and axial types depending on their transformation under parity [87]. Our focus here is on polar (which are strongly coupled to the fluid) and non-radial (dominant in the GW emission) oscillations. The theory of non-radial oscillations of non-rotating NSs, along with the mathematical formalism in the context of GR, was developed in the pioneering work of Thorne [87]. In GR, these modes are quasi-normal, characterized by complex eigenfrequencies ( ω ). The real part of ω corresponds to the mode frequency ( ν ), while the imaginary part is related to the gravitational damping time ( τ ).

A simplified approach known as the relativistic Cowling approximation, which neglects metric perturbations, is frequently used to solve for the oscillation modes, providing information solely on the oscillation frequency. While GR offers a comprehensive understanding of oscillation modes and GW damping times, it is more complex and computationally involved. Meanwhile, the simplified relativistic Cowling approximation has been used to draw many qualitative conclusions and study PNS oscillations. The accuracy of the oscillation mode solutions in the relativistic Cowling approximation of the PNS has been demonstrated in [67]. Similar to the cold NSs, the relativistic Cowling approximation can overestimate the f -mode frequency by ∼ 20 % compared to that obtained in the GR formalism. Further investigations involve searching for the dependence of the mode properties on the outer boundary conditions [88] and numerical dimensions [64]. Additionally, in the work of [62], the postbounce oscillation modes of PNSs have also been studied in general relativity, including spacetime perturbations.

Several approaches have been developed in the past few decades to solve for the oscillation modes within the linearized theory of GR, including the direct integration method [89, 90], WKB approximation, the CF method by Chandrasekhar and Ferrari [91] and the method of continued fraction [92, 93]. The differential equations for the fluid and metric variables for solving the nonradial modes can be found in [54, 67]. We adopt the direct integration method developed by Lindblom and Detweiler [90] to solve for the non-radial QNM modes. In short, the coupled perturbation equations for perturbed metric and fluid variables are integrated within the NS interior with appropriate boundary conditions. We define the PNS boundary as the radius obtained from the

TOVsolutions. Further, the fluid variables are set to zero outside the star, and the Zerilli wave equation [94] is integrated far away from the star. Further, a search is carried out for the complex mode frequency ( ω = 2 πν + ) corre1 τ sponding to only the outgoing wave solution to the Zerilli equation at infinity. For finding the mode characteristics, we use the numerical methods developed in our previous work [54].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-32d4e48c-5482-4dbf-bce4-983e1eef0895.jpg)

/circledot

(a) f -mode frequency v/s mass

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-d9eb111b-899c-49fb-8194-aaac718337b4.jpg)

/circledot

(b) p 1 -mode frequency v/s mass

FIG. 3: f ( p 1 )-mode frequencies as function(s) of mass. The nuclear saturation properties are the same as in fig. 1.

We display f - and p 1 -mode frequencies of some specific thermal cases in Fig. 3a and 3b respectively calculated in general relativistic treatment. The p 1 -mode frequencies are highest for the cold configuration and they are found to decrease with an increase in entropy per baryon and charge fraction. The isothermal configuration for

T = 20 MeV has the lowest p 1 -mode frequencies. A crossing appears between N (1 0 4) and , . N (2 0 2) config-, . urations for the p 1 -mode at high mass. This behavior is also seen in the full parameter space allowed by various constraints (discussed in subsection III C). On the other hand, f -mode frequencies for all cases have some overlapping regions at intermediate and high mass ranges. From fig. 3, it can be concluded that both the f - and p 1 -mode frequencies decrease in thermally configured stars compared to cold NSs. Notably, the thermal effects have a more pronounced impact on the p 1 -mode frequencies than on the f -mode frequencies. The reduction in mode frequencies with thermal effects aligns with the findings of previous studies [56, 95, 96]. This behavior can be attributed to the fact that both f -mode and p -mode frequencies scale inversely with a power of the stellar radius (see [61, 97]). As thermally configured stars of a given mass exhibit lower compactness than their cold counterparts, the increase in radius leads to a reduction in mode frequencies. Thus, the observed variations in frequency are a direct consequence of the underlying differences in compactness between hot and cold configurations (see fig. 2).

## D. Constraining the model parameters

The parameter space of the EoS considered should be consistent with microscopic calculations, state-of-the-art astrophysical observations and terrestrial nuclear experiments. In previous works by some of the authors of this paper [14, 15], a Bayesian-like cut-off scheme described in this section has been used to restrict the parameter space of EoS of cold nucleonic and hyperonic NS matter by imposing the above constraints. In Shirke et al. (2023) [16], this formalism was extended to cold hybrid stars with phase transition using constraints from perturbative QCD (p-QCD). In this work, we have applied a similar scheme to extend the study to finite temperature, as described below.

- 1. Chiral Effective Field Theory ( χ EFT): The χ EFT takes into account nuclear many body interactions using order-by-order expansion in terms of contact interactions and long-range pion exchange interactions and low momentum expansion of nuclear forces related to QCD symmetries are applied. These are based on microscopic calculations from the works of Drischler et al. [13, 98]. These constraints can been used in the low density range 0 07 . -0 20 . fm -3 for cold Pure Neutron Matter (PNM). It should be noted that χ EFT calculations have recently been extended to finite temperature asymmetric nuclear systems [99, 100]. These have also been applied recently to impose constraints on finite temperature EoSs [49]. However as these calculations deal with isothermal (fixed temperature) systems, these results have not been implemented
- in this work as it focuses on isentropic (fixed entropy per baryon) cases.
- 2. Astrophysical Observations (Astro): These are appropriate for the high density part of Asymmetric Nuclear Matter (ANM) and β equilibrated EoS. We impose astrophysical constraints in which the maximum mass supported by a cold neutron star has a lower limit of 2.01 M ⊙ , compatible with the current data. GW170817 observation [18] puts an upper limit on effective tidal deformability ( ˜ ) Λ to be 720 [27] and radius of a 1.4 M ⊙ neutron star to be 13.5 km [101, 102].
- 3. Heavy Ion Collisions (HIC): These constraints have been derived from heavy-ion experiments at GSI in Germany, namely KaoS [103], FOPI [104] and ASY-EOS [105], which leads to constraints at intermediate densities in the range 1-3 n sat . The KaoS experiment can be used to put upper limits on nucleon potential ( U N ) for ANM in the density range 2-3 n sat . The FOPI experiment can be used to constrain binding energy per nucleon ( E/A ) of Symmetric Nuclear Matter (SNM) in the density range 1.4-2 n sat , while ASY-EOS results put bounds on the Symmetry Energy of ANM in the suprasaturation density range 1.1-2 n sat .

TABLE I: Nuclear saturation parameters considered in this work [15]: The top row is the fixed set of the nuclear parameters used in figures 1, 2 and 3 and the bottom row shows the range of varied parameters used in the rest of this study.

|        | n sat (fm - 3 )   | E sat (MeV)   | K sat (MeV)   | J sym (MeV)   | L sym (MeV)   | m ∗ /m   |
|--------|-------------------|---------------|---------------|---------------|---------------|----------|
| Fixed  | 0.15              | -16.0         | 240           | 32            | 60            | .65      |
| Varied | 0.14 -            | -16.2 -       | 200 -         | 28 -          | 40 -          | 0.55 -   |
| Varied | 0.17              | -15.8         | 300           | 34            | 70            | 0.75     |

We provide the uncertainty ranges in nuclear saturation parameters from state-of-the-art experimental data in Table I As in [14, 15], a hard cut-off filter scheme has been used, in which nuclear saturation properties are varied within their allowed uncertainty range to generate a flat prior and only the results which are within the permitted uncertainty region of the above mentioned constraints are accepted. Imposing these constraints reduces the parameter space of nuclear saturation properties and in turn also of the couplings used in constructing the NL-RMF Lagrangian 1.

## III. RESULTS

The composition as well as thermodynamical properties of the neutron star change with the evolving temperature during the PNS formation or neutron star merger. In PNS early evolution the star goes through different

stages in which the core evolves from a low temperature (low entropy state) and gets heated up (high entropy) with a decrease in charge fraction [29, 56, 106]. It then starts cooling and eventually becomes cold and attains β equilibrium. In case of BNS mergers a recent simulation has shown that the entropy per baryon can reach upto values of 3 [34]. To investigate the effects of temperature and out-ofβ equilibrium, in this work we consider two thermal configurations, N S/A ( = 1 , Y Q = 0 4) and . N S/A ( = 2 , Y Q = 0 2) representing two different stages . of evolution of a hot PNS. The former state corresponds to a hot lepton/electron rich configuration while the latter corresponds to a hotter core with deleptonized configuration [56]. We also show comparison with cold β equilibrated state in our study.

## A. Equation of State

We investigate the effect of non-zero temperature and out-ofβ equilibrium condition on the macroscopic NS properties. The results after imposing constraints from χ EFT, Astro and HIC constraints are shown for the thermal case N (2 0 2) , . in fig. 4. The bands of EoSs corresponds to different sets of nuclear saturation properties consistent with filters mentioned above. The imposition of Astro and HIC constraints eliminates certain stiffer EoSs from the parameter set as evident from the figure.

FIG. 4: EoS bands for N (2 0 2) generated after imposing , . χ EFT, Astro and HIC constraints on nuclear saturation properties.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-d174f8ab-0309-45e9-add9-669539db03cb.jpg)

The temperature profiles in the core of these stars using the χ EFT+Astro constraints are shown in fig. 5. For low density ( n B &lt; 0 1 fm . -3 ) there is no significant variation in the temperature profiles for the thermal cases. However, saturation properties do induce appreciable changes in temperature at intermediate densities. For N (1 , 0 4) . the variation is limited to about n B &lt; 0 9 . fm -3 while that of N (2 0 2) extends to the entire density range con-, .

sidered. fig. 6 shows temperature profiles after inclusion of HIC filter and we see that temperature spreads across intermediate and high densities have reduced. However qualitative difference of spreads between N (1 0 4) , . and N (2 0 2) across densities still persists. , .

FIG. 5: Core temperature profiles obtained on imposing χ EFT and Astro constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-3e340a8f-c5a8-4f6e-85cc-ca78dc8761a1.jpg)

FIG. 6: Core temperature profiles obtained on imposing χ EFT + Astro + HIC constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-d038fe07-941b-4dc0-82fd-0d61b02fcdec.jpg)

## B. Mass-Radius Relations

The Mass-Radius relations for N (2 , 0 2) are shown in . fig. 7. The filter imposed by Astro constraints on stiffer EoSs translates to lower values of maximum masses and radii supported by the configurations. We will see later on in subsection III E that the maximum mass is strongly anti-correlated with effective nucleon mass ( m /m ∗ ) and

therefore this parameter controls the stiffness of the EoS. It also strongly effects other macroscopic properties such as radii and tidal deformability.

FIG. 7: Mass-Radius posterior plots for N (2 0 2). , . Constraints are same as in fig. 4. We observe that maximum mass supported after application of different filters decreases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-9d0dc9e7-79a1-40e1-a2ad-a4626374c2c2.jpg)

FIG. 8: Mass-Radius configuration for both the thermal cases along with N T ( = 0) β equilibrated matter. The constraints used are χ EFT + Astro. The solid lines represent the stiffest and softest configurations for each case.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-ed08fa5c-55e3-4dbb-9ce2-285066cb3571.jpg)

We also display these results for the thermal case N (1 0 4) and their comparison with cold , . β equilibrated matter in fig. 8. One can see that the cold β equilibrated configuration is the most compact one and the highest maximum mass it supports is almost identical to N (2 0 2). , . The N (2 , 0 2) configuration is more compact . for higher masses while for N (1 0 4) this is the case for , . low and intermediate mass configurations. After impo- sition of HIC constraints (fig. 9), the maximum masses supported by each thermal configuration further reduce. Most notably M max &gt; 2 5 . M ⊙ configurations are eliminated for all three cases. We discuss in subsection III E how the imposition of these constraints affect masses, radii and tidal deformability.

FIG. 9: Same as fig. 8. The constraints used are χ EFT + Astro + HIC.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-7617575f-a21f-4428-90ed-ec646a73fe24.jpg)

## C. f -mode frequency and damping time

As discussed in subsection II C, after imposing the constraints mentioned above, the f -mode frequencies and damping times are calculated. The results for the case N (2 0 2) are shown in fig. 10 and fig. 11. , .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-61deede3-6887-4dac-9588-988201e704c1.jpg)

/circledot

FIG. 10: f -mode frequency v/s mass posterior plots for N (2 0 2). , . Constraints are same as in fig. 4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-58dd0c02-0ef9-4aff-87f8-bda534677d89.jpg)

/circledot

FIG. 11: f -mode damping time v/s mass posterior plots for N (2 , 0 2). . Constraints are same as in fig. 4.

Next we show the combined f -mode spectrum for all the three cases with χ EFT + Astro constraints in fig. 12. From the figure we can see that f -mode frequencies are lower for low mass ( M ≲ 1 1 M . ⊙ ) NSs for both the thermal cases in comparison with cold and β equilibrated stars. For higher masses all the three configurations have overlapping regions similar to their mass-radius relations (fig. 8). The effective nucleon mass again plays a vital role in determining f -mode frequencies and damping times as have discussed in subsection III E. n sat also plays some role in f -mode oscillations of intermediate mass configurations (see subsection III E).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-33b0863d-960a-439f-8e04-b11b9bb11c1d.jpg)

/circledot

FIG. 12: f -mode frequencies with χ EFT + Astro constraints.

Fig. 13 shows damping timescales for the three cases. Cold neutron stars have the fastest damping for low mass (M ≲ 1 1M . ⊙ ) cases. N (2 0 2) , . configuration have the slowest damping in low and intermediate mass regions.

For higher masses we again have overlapping regions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-58cd79be-b9ac-4395-893b-77ca07d3cbdb.jpg)

/circledot

FIG. 13: Damping of f -modes from χ EFT + Astro constraints. We see overlapping regions in high mass configurations.

The f -mode frequencies and damping times obtained after applying HIC constraints are shown in figs. 14 and 15. These constraints eliminate certain regions of figs. 12 and 13, however the qualitative behaviour does not change much. We see that the frequencies and damping times are inversely correlated for the thermal cases, a fact which is already known for cold NSs. The correlations of f -mode frequencies and damping times with nuclear saturation properties are discussed in III E.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-b6ee82af-e1bb-49b7-91e5-a919cf738b31.jpg)

/circledot

FIG. 14: f -mode frequencies with χ EFT + Astro + HIC constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-e8d7048c-c573-4452-9660-9c61194d7cf7.jpg)

/circledot

FIG. 15: Damping of f -modes from χ EFT + Astro + HIC constraints.

## D. Universal relations

In this section, we discuss the effects of nuclear saturation properties on the Universal Relations (URs) of hot NSs. We focus on the two well-studied URs of cold neutron stars, i.e., Love-Compactness (Λ -C ) and f -mode frequency-Love ( f -Λ) relations. In addition to the two thermal cases considered earlier, we also show the variation of these URs with a wider range of thermal conditions.

## 1. C -Love Relations

The EoS independent UR between dimensionless tidal deformability (Λ) and compactness ( C ) was studied in Maselli et al., where they provide the fitting coefficients in equation (4) of their paper [72]. We display ln(Λ) -C plots from EoSs generated using constraints from χ EFT + Astro in fig. 16. We see that like cold β equilibrated stars ( N T ( = 0)), the plots for two thermal cases also fall in their respective curves which are distinct from that of N T ( = 0) curve, indicating that these configurations do not follow the same URs. The deviations from the cold case are more for higher values of Λ. We also see that compactness corresponding to a particular Λ is highest for N T ( = 0) followed by N (2 0 2) and , . N (1 0 4) indicat-, . ing the order in compactness in the three configurations. The C -ln(Λ) polynomial fitting can be given as -

$$C = \sum _ { i = 0 } ^ { 2 } a _ { i } [ \ln ( \Lambda ) ] ^ { i } \. \quad \quad ( 1 7 ) \quad \frac { \rfloor } { \tilde { \ c h i } }$$

FIG. 16: Λ v/s C universality for the three cases.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-43e668d9-baee-4ed1-bd64-4a7d231487af.jpg)

TABLE II: Fit relations for C -ln(Λ) UR for the two thermal configurations

|               | a 0              | a 1                | a 2              |
|---------------|------------------|--------------------|------------------|
| N (1 , 0 . 4) | 3 . 556 × 10 - 1 | - 4 . 209 × 10 - 2 | 1 . 173 × 10 - 3 |
| N (2 , 0 . 2) | 3 . 598 10 - 1   | 4 . 070 10 - 2     | 1 . 182 10 - 2   |

×

-

×

×

To further investigate how the thermal effects break this universality, we separately vary entropy per baryon ( S/A ) and charge fractions ( Y Q ) keeping nuclear saturation properties fixed. First we show the variation on changing S/A from 1 -2 by fixing Y Q = 0 4, and in . β equilibrium. Next we fix S/A = 1 and vary Y Q in the range 0 01 . -0 40. .

FIG. 17: Effect of S/A on Λ -C universality. Nuclear saturation parameters are fixed to same values as in fig. 1. S/A has been varied between 1 -2 in steps of 0.01.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-9b7f7e94-a6a7-426a-8220-39520a27262b.jpg)

It is clear that Λ -C relation depends not only on charge fraction but also on whether the matter is in β

equilibrium when we vary S/A . The spread in universality in fig. 17 is much less than that of fig. 18. From the results in fig. 17 one can obtain universal fits similar to the previous works [72, 75]. We provide the universal fits for the two thermal configurations in Table II. On the other hand, we see from fig. 18 that Y Q , which is a measure of out-ofβ equilibrium, is mainly responsible for breaking this universality and there is no possibility to obtain a fit relation.

FIG. 18: Effect of Y Q on Λ -C universality. Nuclear saturation properties are fixed as in fig. 1. Y Q has been varied between 0.01-0.40 in steps of 0.01.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-6ce8b411-d2d4-494e-9970-7e3b4d8a855a.jpg)

## 2. f -Love Relations

It is known that for cold NSs f -mode frequencies scale with certain combinations of astrophysical observables, such as M R , and Λ. It was shown that f -mode frequency has an approximately linear relation with square root of average density ( √ M/R 3 ) while the normalized damping time ( τ f R /M 4 3 ) scales with stellar compactness [107]. These fit relations are found in the works of Andersson and Kokkotas [107], Benhar and Ferrari [108], Doneva et al. [109], Pradhan and Chatterjee [83].

However, it turns out that these fits are composition dependent and have large deviations depending on the EoSs considered. Pradhan et al. [54] also provided other improved universal fits which correlate mass-scaled frequencies ( Mω ) (real and imaginary) with logarithm of dimensionless tidal deformability (ln(Λ)). This polynomial fitting can be given as -

$$M \omega = \sum _ { i = 0 } ^ { 5 } \alpha _ { i } [ \ln ( \Lambda ) ] ^ { i } \. \quad \quad ( 1 8 )$$

Here we examine this f -Love universality in the context of hot NSs. The f -Love relations after imposing constraints from χ EFT and Astro are shown in figs. 19a and 19b. Unlike fig. 16, we observe negligible spreads for both thermal cases. Small deviations from the universality appear for NSs with Λ &gt; e 8 ∼ 3000 or with masses close to or less than solar mass. After the postbounce stage, a PNS becomes more and more compact as it evolves, consequently Λ goes from high to low and the star approaches a cold β equilibrated state. The f -Love fit coefficients for the two thermal configurations are given in Table III.

These findings are reasonable, as at higher Λ the star can be assumed to be in the N (1 , 0 4) representative state . and this configuration shows the highest deviation from cold f -Love universality. The N (2 0 2) , . state is representative of intermediate Λ and it has close resemblance with the cold β equilibrated configuration. The close resemblance these URs have with N T ( = 0) promotes us to study them across ranges of S/A and Y Q . The conclusions are somewhat different to those for C -Love relations (see Appendix sections A and B).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-7d21d5b8-2261-4d54-b00f-eb5fd9fa0758.jpg)

(a) Real(

Mω v/s

)

ln(Λ)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-ad061d88-2b8f-44ba-ac09-7d76375ffa43.jpg)

(b) Im(

Mω v/s

)

ln(Λ)

FIG. 19: f -Love universality for the three cases. The relations closely resembles N T ( = 0) state and deviations are seen only for high Λ or consequently less compact and less massive configurations.

TABLE III: Fit relations for f -Love UR for the two thermal configurations

|               | α 0                                   | α 1                                     | α 2                                     | α 3                                   | α 4                                     | α 5                                   |
|---------------|---------------------------------------|-----------------------------------------|-----------------------------------------|---------------------------------------|-----------------------------------------|---------------------------------------|
| N (1 , 0 . 4) | 1 . 812 × 10 - 1 + i 2 . 237 × 10 - 5 | - 5 . 435 × 10 - 3 + i 5 . 599 × 10 - 5 | - 4 . 999 × 10 - 3 - i 1 . 863 × 10 - 5 | 7 . 106 × 10 - 4 + i 2 . 142 × 10 - 6 | - 3 . 798 × 10 - 5 - i 1 . 004 × 10 - 7 | 7 . 034 × 10 - 7 + i 1 . 452 × 10 - 9 |
| N (2 , 0 . 2) | 1 . 809 × 10 - 1 + i 1 . 411 10 - 5   | - 5 . 110 × 10 - 3 + i 6 . 519 10 - 5   | - 5 . 141 × 10 - 3 i 2 . 239 10 - 5     | 7 . 405 × 10 - 4 + i 2 . 849 10 - 6   | - 4 . 114 × 10 - 5 i 1 . 619 10 - 7     | 8 . 486 × 10 - 7 + i 3 . 466 10 - 9   |

×

×

-

×

×

-

×

×

## E. Correlation studies

From the correlation studies in the case of cold β equilibrated matter, it was found that the most important nuclear saturation parameter which affects the NS observable properties is the nucleon effective mass m ∗ and to a certain extent the nuclear saturation density n sat [14, 54]. In this section we will probe possible correlations of the nuclear parameters with NS observables and also among each for the thermal cases. First we study the effect of imposing the χ EFT + Astro constraints on the correlations. Next, we also include the HIC constraints, and investigate how the correlations get altered. As the HIC results are known to be dependent on the transport models, we have treated them separately with caution (see [14] for a more detailed discussion).

## 1. χ EFT + Astro

The correlations among the nuclear saturation parameters and the NS observable properties on imposing χ EFT + Astro constraints are shown for the two thermal cases in figs. 20 and 21. We describe the findings in detail below.

## 1. Correlations among nuclear saturation properties:

- · Effective nucleon mass m ∗ shows moderate anti-correlation with n sat .
- · Relatively strong correlation is observed between J sym and L sym . This is consistent with previous studies [14].
- · E sat and K sat do not show correlations with other saturation properties.

## 2. Correlations of nuclear saturation properties with NS observables:

- · For both thermal cases, n sat shows moderate correlations/anti-correlations with properties ( R , Λ, ν f , τ f ) of canonical 1 4M . ⊙ hot NSs.
- · On the other hand n sat shows weak correlations/anti-correlations with properties ( R , Λ, ν f , τ f ) of canonical 2M ⊙ hot NSs and it does not have any effect on maximum mass supported by the corresponding EoS.
- · Effective nucleon mass m ∗ has strong anticorrelation with M max , it also shows relatively stronger correlations/anti-correlations with R 2M ⊙ , Λ 2M ⊙ , ν f, 2M ⊙ and τ f, 2M ⊙ for both cases.
- · However for 1 4M . ⊙ hot NSs, m /m ∗ has weak correlations/anti-correlations with the observables for the case N (1 , 0 4) . but moderate correlations/anti-correlations for N (2 0 2) ob-, . servables. Previous study with cold NSs have also shown moderate correlations of m /m ∗ with 1 4M . ⊙ observables [14]. We can conclude that during the early phase of hot NS with low entropy and high lepton fraction, m /m ∗ have lesser effect on properties of intermediate mass configurations.

## 3. Correlations among NS observables:

- · Maximum mass is strongly correlated/anticorrelated with properties of canonical 2M ⊙ observables for both cases.
- · With observables of 1 4M . ⊙ , these correlations are moderate for the state N (1 , 0 4) and . strong for the case N (2 , 0 2). .
- · Correlations among radii, Λ, f -mode frequencies and damping times are similar in both cases and they are similar to that of cold NSs [54].

## 2. χ EFT + Astro + HIC

The correlations on imposing χ EFT + Astro + HIC constraints are shown for the two thermal cases in figs. 22 and 23. We describe the key results and differences from χ EFT + Astro filters below.

## 1. Correlations among nuclear saturation properties:

- · Correlation between J sym and L sym gets weaker.
- · Anti-correlation between m /m ∗ and n sat becomes relatively stronger.
- · K sat and L sym now have moderate correlations with m /m ∗ .

## 2. Correlations of nuclear saturation properties with NS observables:

- · m /m ∗ anti-correlation with maximum mass remains strong for both cases but they are relatively weaker than those of correlations without HIC constraints.
- · Masses, radii, f -mode frequencies and damping times of 1 4M . ⊙ PNSs now have poor correlation with effective mass.
- · Correlations/anti-correlations of effective mass with properties ( R , Λ, ν f , τ f ) of canon-
- ical 2M ⊙ PNSs become moderate for both cases.
- · Correlations/anti-correlations of n sat with observables ( R , Λ, ν f , τ f ) of canonical 1 4M . ⊙ PNSs become stronger for both cases. However, its correlation with canonical 2M ⊙ hot NS observables remain weak.

## 3. Correlations among NS observables:

- · Correlations/anti-correlations of M max with observables of 1 4M . ⊙ hot NSs become weaker for both cases.

FIG. 20: Correlations among saturation properties and observables for the state N (1 , 0 4) for . χ EFT + Astro constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-74636288-48b1-47a0-ab9b-880b7794d75b.jpg)

FIG. 21: Correlations among saturation properties and observables for the state N (2 , 0 2) for . χ EFT + Astro constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-33ea6baf-6345-4c2c-8d51-d06383f23096.jpg)

FIG. 22: Correlations among saturation properties and observables for the state N (1 , 0 4) for . χ EFT + Astro + HIC constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-17899fec-c2f7-41cd-8a0a-09f34890ee80.jpg)

FIG. 23: Correlations among saturation properties and observables for the state N (2 , 0 2) for . χ EFT + Astro + HIC constraints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-6088e922-1c0f-4d26-bb86-c7ccec6a3c1e.jpg)

## IV. DISCUSSIONS

Summary of the results: In this work, we develop a systematic and consistent formalism for nucleonic equations of state in neutron stars at finite temperature, by extending a cold matter EoS whose parameter space of allowed uncertainties of the model is restricted by demanding compatibility with current nuclear experimental data, multi-messenger astrophysical observations as well as heavy-ion physics. The developed framework was then applied to the calculation of polar f -modes and p 1 -modes within a general relativistic scheme. The full ranges of allowed nuclear saturation properties for nucleonic matter used in our work allowed for the extensive study of correlations of the nuclear and thermal properties on macroscopic and oscillation properties of hot NSs. We specifically chose a low entropy high charge fraction state and a high entropy deleptonized low charge fraction state and compared it with cold β equilibrated NSs.

On including thermal contributions to the EoS, we found that increasing S/A leads to higher core temperatures while Y Q has an opposite effect and deleptonized configurations are more compact. Increasing S/A results in larger radii and flatter mass-radius relations. These observations are consistent with previous studies [56, 78]. The f -and p 1 -mode frequencies calculated in a full general relativistic scheme are found to be the highest for the cold configuration and decrease with an increase in entropy per baryon and charge fraction.

On imposing χ EFT, Astro and HIC constraints, we found that nuclear saturation parameters affect temperature profiles at intermediate and high densities. There are overlapping regions in the posteriors for M -R f , -mode frequency and damping timescales for high masses among N T ( = 0) , N (1 0 4) , . , N (2 , 0 2) . states. From the correlation study of hot NSs, we found on imposing χ EFT and Astro constraints that the most important saturation properties having effect on astrophysical observables, such as maximum mass, radii, tidal deformation as well as f -mode characteristics, are found to be m /m ∗ followed by n sat . However their correlations are qualitatively different depending on the thermal state after imposing χ EFT + Astro constraints and N (2 0 2) state seems to resemble closely to that of , . cold NS results [14]. After applying HIC constraints, we found that correlations of m ∗ with NS observables are significantly weaker. These findings are qualitatively similar to the results in cold β equilibrated matter case [14, 54]. Hence, the interplay among saturation properties and astrophysical observables depends not only on the thermal configurations considered but also on the constraints imposed. It should be noted that previous study of cold β equilibrated NSs by Pradhan et al. [54] found strong correlations/anti-correlations of effective nucleon mass with all astrophysical observables but did not include constraints from χ EFT and HIC.

Using the posteriors from the constraint filters, we further investigated the role of nuclear saturation parameters on the universal relations for hot NSs, in particular Love-compactness and f -Love relations. We concluded that, hot NS configurations at higher Λ values show significant deviations in C -Love and f -Love relations compared to the cold β equilibrated case, although these variations are more prominent in the former case. Further investigation confirmed that these URs are mostly insensitive to nuclear saturation properties and is mainly affected by charge fraction variation, which is a measure of out-ofβ equilibrium condition.

Comparison with other works: There exist other EoS models in the RMF scheme that have been extended to finite temperature and tested against current constraints. In the extension of the FSU2H parametrization based on the RMF scheme to finite temperature [41], nucleonic and hyperonic matter have been considered. The resulting FSU2H* model is compatible with astrophysical constraints (maximum mass of 2 M ⊙ and radii below 13 km) as well as saturation properties of nuclear matter and finite nuclei, and pressure constraint from heavy-ion collisions. Alford et al. [49] also constructed QMC-RMF EoS models for finite temperature NSs, informed by χ EFT and astrophysical constraints. These works demonstrate the strong influence of thermal corrections on the EoS and thermodynamic properties. However, they only consider fixed RMF parametrizations and therefore the role of the nuclear parameters in governing NS observable properties is not evident from such studies, but this is possible with a systematic variation of nuclear parameters as done in our work. Further, as these calculations are performed for isothermal and not isentropic conditions, it is not easy to estimate their effects on astrophysical observables.

It is interesting to note that, several numerical simulations of CCSN, PNS as well as BNS mergers, employing different EoS models, have also reported the influence of the nucleon effective mass [35, 110-112]. Using a selection of purely nucleonic EoS models from CompOSE (including Skyrme, covariant density functional models, CMF) in agreement with constraints from astrophysical observations and nuclear data, Raduta et al. [113], reviewed thermal properties relevant for PNS and BNS mergers and concluded that they are strongly influenced by the effective nucleon mass. This was further corroborated by using a covariant density functional theory with density dependent couplings employing the posteriors of a Bayesian inference (DDB*), which encompasses a large number of EoS models [57]. These EoSs were also employed in Thapa et al. [56] to calculate the f -and p 1 -mode frequencies within Cowling approximation. Although the thermal and composition effects on the oscillation modes were studied for the selection of EoSs

from CompOSE, the correlations were calculated within the cold EoS approximation, and therefore did not provide any information about the thermal conditions.

A few recent works by Ghosh et al. [114] and Kumar et al. [115] studied oscillations within Cowling approximation including neutrino contributions for hot nucleonic and hybrid neutron stars. It should be noted that Burgio et al. [95] used a model where core and envelope entropy per baryon are different, and the former rises and then falls and the latter falls along with deleptonization during its evolution, to study g 1 -, f -and p 1 -modes using a BHF model extended to finite temperature. In our work, the entropy per baryon has been kept fixed throughout the star in our study.

In the recent work by Tsiopelas et al. [58], which employed a density dependent covariant density functional model, the effect of L sym was found to be diminishing with increasing temperature. In our work using NLRMF formalism we found that after applying χ EFT + Astro filters, L sym has no impact on macroscopic or oscillation properties for N (1 0 4) configuration. , . On the other hand for N (2 , 0 2) configuration, these correlations . are qualitatively similar to that of cold NS results [14].

Previous studies have indicated deviations of ¯ I -¯ Q relations in hot environment [51, 78, 116, 117]. A few studies [78, 79] showed that universality is restored if finite T configurations for different EoS with same S/A and electron fraction are considered. PNS simulations by Sotani et al. [97] explored the prospect of measurement PNS mass and radius using simultaneous measurements of f -(or p 1 -) and w 1 -mode frequencies by relating those to average mass density and stellar compactness respectively. Torres-Forn´ et al. [63] have demonstrated e how g 1 and f -mode frequencies scale with surface gravity in their PNS simulations. In our article, we have shown that C -Love and f -Love URs are mostly not sensitive to nuclear saturation properties but they change depending on the thermal state. In particular charge fraction significantly affects these deviations and they occur at high Λ-values (or low compactness). A recent simulation [118] has also shown that due to thermal effects deviations from URs i.e., f -√ M/R 3 , ¯ I -Λ -¯ , Q f -Love occur during PNS evolution at earlier stages where the configuration has large tidal deformability. It should be noted that in the case of proto-neutron stars, in addition to thermal effects one must also take into account the effect of strong magnetic fields which are known to break the universality of ¯ I -Love -¯ relations [119]. Q

Limitations and future work: It is worth mentioning that, we define the modes following the classification by Cowling [120]. However, a recent investigation done in [65], considering the nucleonic EoS models, demonstrates that in the initial phase (within 0.4 s) of the PNS evolution, the Cowling classification scheme does not hold but after 0.4 sec the classification of oscillation modes with different approaches are in very good agreement. We further note that neutrinos have not been included in this work. Neutrinos may not remain free streaming for the entire temperature regime considered in this study, however we restrict ourselves to this approximation for this work and defer calculations including neutrino trapping to our forthcoming studies.

This study can be extended to include non-nucleonic degrees of freedom such as hyperons, delta baryons and quark matter in hot EoS that satisfy microscopic, experimental and astrophysical constraints. These new degrees of freedom are known to soften the EoSs and might also lead to decrease in core temperature due to entropy being contributed by these new species.

It would also be interesting to include the effect of neutrinos as they are expected to be in equilibrium with the thermal/isentropic bath at high temperatures. Previous studies on BNS mergers [121, 122] have shown that in the temperature regime T ∼ 1 -5 MeV, nuclear matter is transparent to neutrinos due to their long mean free path. At high temperatures (beyond ∼ 5 MeV) neutrinos may be trapped, depending on the neutrino mean free path at a given density. As a first step, in this work we developed the formalism for the ν -transparent case (valid for temperatures upto 5 MeV) and plan to extend this to include trapped neutrinos in a future work.

## ACKNOWLEDGMENTS

The authors are grateful to Suprovo Ghosh for providing the posterior data that was used in this work. D.C. would like to thank Adriana R. Raduta and Mikhail V. Beznogov for their warm hospitality at IFIN-HH, Bucharest which led to several interesting discussions related to this work. D.C. would also like to thank Prof. J. Schaffner-Bielich and his students at ITP Frankfurt for insightful discussions on finite temperature models of neutron stars. D.C. is grateful to Prof Mark Alford, Alexander Haber and their group at Washington University in St. Louis for their hospitality and fruitful exchanges. The authors also thank Prof Luciano Rezzolla for pointing them to his relevant works. N.B., B.K.P. and D.C. acknowledge the usage of the IUCAA HPC computing facility for numerical calculations of oscillation modes.

## V. DATA AVAILABILITY

The posterior data of both the nuclear saturation properties and the coupling constants obtained after imposing χ EFT + Astro and χ EFT + Astro + HIC constraints which support this work can be found at [123]. Alterna-

tively the posterior sets of only nuclear saturation properties are available in the arXiv ancillary files of previously published work of Ghosh et al. [14].

## Appendix A: Effect of entropy per baryon on f -Love relations

We fix the nuclear saturation properties as in fig. 1. Figs. 24a and 24b show effect of changing S/A between 1 -2 in Y Q = 0 4 and . β equilibrated states. From fig. 24a, we see that Y Q = 0 4 case shows deviations from . β equilibtrated hot NSs for higher Λ and they are overlapping otherwise. This is in contrast to Λ -C URs shown in fig. 17 where we did not get any overlapping region. The relation between Im( Mω ) and ln(Λ) are shown in fig. 24b also establishes deviations at high Λ values. There are also qualitatively negligible deviations for both real and imaginary parts across the entropy per baryon range considered as compared to Λ -C URs.

FIG. 24: Effect of entropy per baryon on f -Love universality.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-14f0411a-c6ac-4f1a-a9e7-4ac8931d53da.jpg)

These deviations shown between states with β equi- libration and Y Q = 0 4 clearly shows that at higher Λ . universality is broken by out-ofβ equilibrium effects. The spreads however in figs. 24a and 24b are qualitatively much lower than that of Λ -C relations shown in fig. 18 where deviations are observed across the entire range of dimensionless tidal deformability.

## Appendix B: Effect of charge fraction on f -Love relations

Here we fix entropy per baryon at S/A = 1 and vary Y Q from 0.01 - 0.40. We see in figs. 25a and 25b that this has no effect on both real and imaginary parts of f -Love relations for low Λ values and at higher values there are small spreads. These deviations appear for Λ &gt; e 8 ∼ 3000. The conclusion is that charge fraction is responsible for deviations in the universality like C -Love relations albeit to a much lesser degree and only at higher Λ.

FIG. 25: Effect of charge fraction on f -Love universality.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00739v2-figure-0af0c4a3-647c-489a-a888-a2e4a273ba44.jpg)

- [1] N. K. Glendenning, Compact Stars (Springer New York, NY, 1997).
- [2] J. Schaffner-Bielich, Compact Star Physics (Cambridge University Press, 2020).
- [3] J. M. Lattimer and M. Prakash, The physics of neutron stars, Science 304 , 536 (2004).
- [4] J. M. Lattimer, The nuclear equation of state and neutron star masses, Annual Review of Nuclear and Particle Science 62 , 485 (2012).
- [5] D. Chatterjee and I. Vida˜a, Do hyperons exist in the n interior of neutron stars?, The European physical journal A 52 , 1 (2016).
- [6] M. Oertel, M. Hempel, T. Kl¨hn, and S. Typel, Equaa tions of state for supernovae and compact stars, Rev. Mod. Phys. 89 , 015007 (2017).
- [7] G. Taranto, M. Baldo, and G. F. Burgio, Selecting microscopic equations of state, Phys. Rev. C 87 , 045803 (2013).
- [8] H. Tong, C. Wang, and S. Wang, Nuclear matter and neutron stars from relativistic brueckner-hartree-fock theory, The Astrophysical Journal 930 , 137 (2022).
- [9] H.-J. Schulze, Brueckner-hartree-fock approach to neutron star structure, Progress of Theoretical Physics Supplement 186 , 57 (2010).
- [10] J. H. C. Drischler and C. Wellenhofer, Chiral effective field theory and the high-density nuclear equation of state, Annual Review of Nuclear and Particle Science 71 , 403 (2021).
- [11] N. Hornick, L. Tolos, A. Zacchi, J.-E. Christian, and J. Schaffner-Bielich, Relativistic parameterizations of neutron matter and implications for neutron stars, Physical Review C 98 , 065804 (2018).
- [12] S. Typel and H. Wolter, Relativistic mean field calculations with density-dependent meson-nucleon coupling, Nuclear Physics A 656 , 331 (1999).
- [13] C. Drischler, A. Carbone, K. Hebeler, and A. Schwenk, Neutron matter from chiral two- and three-nucleon calculations up to n LO, Q. J. 3 Mech. Appl. Math. 94 , 054307 (2016).
- [14] S. Ghosh, D. Chatterjee, and J. Schaffner-Bielich, Imposing multi-physics constraints at different densities on the neutron star equation of state, The European physical journal A 58 , 1 (2022).
- [15] S. Ghosh, B. K. Pradhan, D. Chatterjee, and J. Schaffner-Bielich, Multi-physics constraints at different densities to probe nuclear symmetry energy in hyperonic neutron stars, Frontiers in Astronomy and Space Sciences 9 , 864294 (2022).
- [16] S. Shirke, S. Ghosh, and D. Chatterjee, Constraining the equation of state of hybrid stars using recent information from multidisciplinary physics, The Astrophysical Journal 944 , 7 (2023).
- [17] M. V. Beznogov and A. R. Raduta, Bayesian inference of the dense matter equation of state built upon covariant density functionals, Phys. Rev. C 107 , 045803 (2023).
- [18] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Gw170817: Observation of gravitational waves from a binary neutron star inspiral, Phys. Rev. Lett. 119 , 161101 (2017).
- [19] B. P. Abbott et al. (The LIGO Scientific Collaboration and the Virgo Collaboration), Gw170817: Measure-
- ments of neutron star radii and equation of state, Phys. Rev. Lett. 121 , 161101 (2018).
- [20] B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Properties of the binary neutron star merger gw170817, Phys. Rev. X 9 , 011001 (2019).
- [21] C. Chirenti, R. Gold, and M. C. Miller, Gravitational waves from f-modes excited by the inspiral of highly eccentric neutron star binaries, The Astrophysical Journal 837 , 67 (2017).
- [22] J. Steinhoff, T. Hinderer, T. Dietrich, and F. Foucart, Spin effects on neutron star fundamental-mode dynamical tides: Phenomenology and comparison to numerical simulations, Phys. Rev. Res. 3 , 033129 (2021).
- [23] M. C. Miller et al., The radius of psr j0740+6620 from nicer and xmm-newton data, The Astrophysical Journal Letters 918 , L28 (2021).
- [24] T. E. Riley et al., A nicer view of the massive pulsar psr j0740+6620 informed by radio timing and xmm-newton spectroscopy, The Astrophysical Journal Letters 918 , L27 (2021).
- [25] M. C. Miller et al., Psr j0030+0451 mass and radius from nicer data and implications for the properties of neutron star matter, The Astrophysical Journal Letters 887 , L24 (2019).
- [26] T. E. Riley et al., A nicer view of psr j0030+0451: Millisecond pulsar parameter estimation, The Astrophysical Journal Letters 887 , L21 (2019).
- [27] H. Tong, P. Zhao, and J. Meng, Symmetry energy at supra-saturation densities via the gravitational waves from gw170817, Phys. Rev. C 101 , 035802 (2020).
- [28] S. De, D. Finstad, J. M. Lattimer, D. A. Brown, E. Berger, and C. M. Biwer, Tidal deformabilities and radii of neutron stars from the observation of gw170817, Phys. Rev. Lett. 121 , 091102 (2018).
- [29] J. A. Pons, S. Reddy, M. Prakash, J. M. Lattimer, and J. A. Miralles, Evolution of proto-neutron stars, The Astrophysical Journal 513 , 780 (1999).
- [30] E. J. Lentz, S. W. Bruenn, W. R. Hix, A. Mezzacappa, O. E. B. Messer, E. Endeve, J. M. Blondin, J. A. Harris, P. Marronetti, and K. N. Yakunin, Three-dimensional core-collapse supernova simulated using a 15 m ⊙ progenitor, The Astrophysical Journal Letters 807 , L31 (2015).
- [31] H.-T. Janka, K. Langanke, A. Marek, G. Mart´ ınezPinedo, and B. M¨ller, Theory of core-collapse superu novae, Physics Reports 442 , 38 (2007), the Hans Bethe Centennial Volume 1906-2006.
- [32] A. Burrows and D. Vartanyan, Core-collapse supernova explosion theory, Nature 589 , 29 (2020).
- [33] L. Baiotti and L. Rezzolla, Binary neutron star mergers: a review of einstein's richest laboratory, Reports on Progress in Physics 80 , 096901 (2017).
- [34] E. R. Most, A. Motornenko, J. Steinheimer, V. Dexheimer, M. Hanauske, L. Rezzolla, and H. Stoecker, Probing neutron-star matter in the lab: Similarities and differences between binary mergers and heavy-ion collisions, Phys. Rev. D 107 , 043034 (2023).
- [35] J. Fields, A. Prakash, M. Breschi, D. Radice, S. Bernuzzi, and A. da Silva Schneider, Thermal effects in binary neutron star mergers, A. P. J. L. 952 , L36 (2023).

- [36] D. Radice, S. Bernuzzi, and A. Perego, The dynamics of binary neutron star mergers and gw170817, Annual Review of Nuclear and Particle Science 70 , 95 (2020).
- [37] N. Sarin and P. Lasky, The evolution of binary neutron star post-merger remnants: a review, General Relativity and Gravitation 53 , 59 (2021).
- [38] L. Baiotti, Gravitational waves from neutron star mergers and their relation to the nuclear equation of state, Progress in Particle and Nuclear Physics 109 , 103714 (2019).
- [39] K. P. Nunna, S. Banik, and D. Chatterjee, Signatures of strangeness in neutron star merger remnants, The Astrophysical Journal 896 , 109 (2020).
- [40] A. R. Raduta, Equations of state for hot neutron starsii. the role of exotic particle degrees of freedom, Eur. Phys. J. A 58 , 1 (2022).
- [41] H. Kochankovski, A. Ramos, and L. Tolos, Equation of state for hot hyperonic neutron star matter, Monthly Notices of the Royal Astronomical Society 517 , 507 (2022).
- [42] S. Blacker, H. Kochankovski, A. Bauswein, A. Ramos, and L. Tolos, Thermal behavior as indicator for hyperons in binary neutron star merger remnants, Phys. Rev. D 109 , 043015 (2024).
- [43] S. Blacker, A. Bauswein, and S. Typel, Exploring thermal effects of the hadron-quark matter transition in neutron star mergers, Phys. Rev. D 108 , 063032 (2023).
- [44] J. M. Lattimer and D. F. Swesty, A generalized equation of state for hot, dense matter, Nuclear Phsyics A 535 , 331 (1991).
- [45] H. Shen, H. Toki, K. Oyamatsu, and K. Sumiyoshi, Relativistic equation of state of nuclear matter for supernova and neutron star, Nuclear Phsyics A 637 , 435 (1998).
- [46] M. Hempel, T. Fischer, J. Schaffner-Bielich, and M. Liebend¨rfer, New Equations of State in Simulations o of Core-collapse Supernovae, The Astrophysical Journal 748 , 70 (2012).
- [47] A. W. Steiner, M. Hempel, and T. Fischer, Core-collapse supernova equations of state based on neutron star observations, The Astrophysical Journal 774 , 17 (2013).
- [48] F. J. Fattoyev, C. J. Horowitz, J. Piekarewicz, and G. Shen, Relativistic effective interaction for nuclei, giant resonances, and neutron stars, Phys. Rev. C 82 , 055803 (2010).
- [49] M. G. Alford, L. Brodie, A. Haber, and I. Tews, Tabulated equations of state from models informed by chiral effective field theory, Physica Scripta 98 , 125302 (2023).
- [50] S. Banik, M. Hempel, and D. Bandyopadhyay, New hyperon equations of state for supernovae and neutron stars in density-dependent hadron field theory, The Astrophysical Journal Supplement Series 214 , 22 (2014).
- [51] M. Marques, M. Oertel, M. Hempel, and J. Novak, New temperature dependent hyperonic equation of state: Application to rotating neutron star models and i -q relations, Phys. Rev. C 96 , 045806 (2017).
- [52] P. Papazoglou, D. Zschiesche, S. Schramm, J. SchaffnerBielich, H. St¨cker, and W. Greiner, Nuclei in a chiral o su(3) model, Phys. Rev. C 59 , 411 (1999).
- [53] V. Dexheimer and S. Schramm, Proto-neutron and neutron stars in a chiral su(3) model, The Astrophysical Journal 683 , 943 (2008).
- [54] B. K. Pradhan, D. Chatterjee, M. Lanoye, and P. Jaikumar, General relativistic treatment of f -mode oscilla-
- tions of hyperonic stars, Physical Review C 106 , 015805 (2022).
- [55] B. K. Pradhan, D. Chatterjee, R. Gandhi, and J. Schaffner-Bielich, Role of vector self-interaction in neutron star properties, Nuclear Physics A 1030 , 122578 (2023).
- [56] V. B. Thapa, M. V. Beznogov, A. R. Raduta, and P. Thakur, Frequencies of f - and p -oscillation modes in cold and hot compact stars, Phys. Rev. D 107 , 103054 (2023).
- [57] A. R. Raduta, M. V. Beznogov, and M. Oertel, Bayesian inference of thermal effects in dense matter within the covariant density functional theory, Physics Letters B 853 , 138696 (2024).
- [58] S. Tsiopelas, A. Sedrakian, and M. Oertel, Finitetemperature equations of state of compact stars with hyperons: three-dimensional tables (2024), arXiv:2406.00484 [nucl-th].
- [59] J.-J. Li and A. Sedrakian, New covariant density functionals of nuclear matter for compact star simulations, A. P. J. 957 , 41 (2023).
- [60] H. Sotani and T. Takiwaki, Gravitational wave asteroseismology with protoneutron stars, Phys. Rev. D 94 , 044043 (2016), arXiv:1608.01048 [astro-ph.HE].
- [61] H. Sotani, Gravitational waves from protoneutron stars and nuclear equation of state, Astron. Nachr. 340 , 217 (2019).
- [62] A. Torres-Forn´, e P. Cerd´ a-Dur´n, a A. Passamonti, M. Obergaulinger, and J. A. Font, Towards asteroseismology of core-collapse supernovae with gravitational wave observations - II. Inclusion of space-time perturbations, Mon. Not. Roy. Astron. Soc. 482 , 3967 (2019).
- [63] A. Torres-Forn´, e P. Cerd´ a-Dur´n, a M. Obergaulinger, B. M¨ uller, and J. A. Font, Universal Relations for Gravitational-Wave Asteroseismology of Protoneutron Stars, Phys. Rev. Lett. 123 , 051102 (2019), [Erratum: Phys.Rev.Lett. 127, 239901 (2021)], arXiv:1902.10048 [gr-qc].
- [64] H. Sotani and T. Takiwaki, Dimension dependence of numerical simulations on gravitational waves from protoneutron stars, Phys. Rev. D 102 , 023028 (2020), arXiv:2004.09871 [astro-ph.HE].
- [65] M. C. Rodriguez, I. F. Ranea-Sandoval, C. Chirenti, and D. Radice, Three approaches for the classification of protoneutron star oscillation modes, Mon. Not. Roy. Astron. Soc. 523 , 2236 (2023), arXiv:2304.00033 [astroph.HE].
- [66] H. Sotani, B. M¨ller, and T. Takiwaki, Universality in u supernova gravitational waves with proto-neutron star properties (2024), arXiv:2405.09030 [astro-ph.HE].
- [67] H. Sotani and T. Takiwaki, Accuracy of the relativistic cowling approximation in protoneutron star asteroseismology, Phys. Rev. D 102 , 063025 (2020).
- [68] A. Torres-Forn´, e P. Cerd´ a-Dur´n, a A. Passamonti, and J. A. Font, Towards asteroseismology of corecollapse supernovae with gravitational-wave observations -i. cowling approximation, Monthly Notices of the Royal Astronomical Society 474 , 5272 (2017), https://academic.oup.com/mnras/articlepdf/474/4/5272/23231251/stx3067.pdf.
- [69] H. Sotani and T. Takiwaki, Avoided crossing in gravitational wave spectra from protoneutron star, Mon. Not. Roy. Astron. Soc. 498 , 3503 (2020), arXiv:2008.00419 [astro-ph.HE].

- [70] T. Bruel, M.-A. Bizouard, M. Obergaulinger, P. Maturana-Russel, A. Torres-Forn´, P. Cerd´-Dur´n, e a a N. Christensen, J. A. Font, and R. Meyer, Inference of protoneutron star properties in core-collapse supernovae from a gravitational-wave detector network, Phys. Rev. D 107 , 083029 (2023), arXiv:2301.10019 [astro-ph.HE].
- [71] M.-A. Bizouard, P. Maturana-Russel, A. Torres-Forn´, e M. Obergaulinger, P. Cerd´ a-Dur´n, a N. Christensen, J. A. Font, and R. Meyer, Inference of protoneutron star properties from gravitational-wave data in corecollapse supernovae, Phys. Rev. D 103 , 063006 (2021), arXiv:2012.00846 [gr-qc].
- [72] A. Maselli, V. Cardoso, V. Ferrari, L. Gualtieri, and P. Pani, Equation-of-state-independent relations in neutron stars, Phys. Rev. D 88 , 023007 (2013).
- [73] K. Yagi and N. Yunes, I-love-q relations in neutron stars and their applications to astrophysics, gravitational waves, and fundamental physics, Phys. Rev. D 88 , 023009 (2013).
- [74] B. K. Pradhan, D. Chatterjee, and D. E. AlvarezCastillo, Probing hadron-quark phase transition in twin stars using f-modes, Monthly Notices of the Royal Astronomical Society 531 , 4640 (2024), https://academic.oup.com/mnras/articlepdf/531/4/4640/58296587/stae1425.pdf.
- [75] B. K. Pradhan, A. Vijaykumar, and D. Chatterjee, Impact of updated multipole love numbers and f -love universal relations in the context of binary neutron stars, Phys. Rev. D 107 , 023010 (2023).
- [76] D.-H. Wen, B.-A. Li, H.-Y. Chen, and N.-B. Zhang, Gw170817 implications on the frequency and damping time of f -mode oscillations of neutron stars, Phys. Rev. C 99 , 045806 (2019).
- [77] B.-A. Li, P. G. Krastev, D.-H. Wen, and N.-B. Zhang, Towards Understanding Astrophysical Effects of Nuclear Symmetry Energy, Eur. Phys. J. A 55 , 117 (2019), arXiv:1905.13175 [nucl-th].
- [78] A. R. Raduta, M. Oertel, and A. Sedrakian, Protoneutron stars with heavy baryons and universal relations, Monthly Notices of the Royal Astronomical Society 499 , 914 (2020).
- [79] S. Khadkikar, A. R. Raduta, M. Oertel, and A. Sedrakian, Maximum mass of compact stars from gravitational wave events with finite-temperature equations of state, Phys. Rev. C 103 , 055811 (2021).
- [80] N. Khosravi Largani, T. Fischer, A. Sedrakian, M. Cierniak, D. E. Alvarez-Castillo, and D. B. Blaschke, Universal relations for rapidly rotating cold and hot hybrid stars, Monthly Notices of the Royal Astronomical Society 515 , 3539 (2022).
- [81] P. Laskos-Patkos, P. S. Koliogiannis, A. Kanakis-Pegios, and C. C. Moustakidis, Thermodynamics of hot neutron stars and universal relations, Universe 8 , 395 (2022).
- [82] S. Jaiswal and D. Chatterjee, Constraining dense matter physics using f-mode oscillations in neutron stars, Physics 3 , 302 (2021).
- [83] B. K. Pradhan and D. Chatterjee, Effect of hyperons on f -mode oscillations in neutron stars, Physical Review C 103 , 035810 (2021).
- [84] S. Typel et al., CompOSE (2022).
- [85] W.-C. Chen and J. Piekarewicz, Building relativistic mean field models for finite nuclei and neutron stars, Phys. Rev. C 90 , 044305 (2014).
- [86] P. Laskos-Patkos, P. S. Koliogiannis, A. Kanakis-Pegios, and C. C. Moustakidis, Universal relations and finite temperature neutron stars, HNPS Adv. Nucl. Phys. 29 , 94 (2023).
- [87] K. S. Thorne and A. Campolattaro, Non-Radial Pulsation of General-Relativistic Stellar Models. I. Analytic Analysis for L &gt; = 2, Astro Phys. Journal 149 , 591 (1967).
- [88] H. Sotani, T. Kuroda, T. Takiwaki, and K. Kotake, Dependence of the outer boundary condition on protoneutron star asteroseismology with gravitationalwave signatures, Phys. Rev. D 99 , 123024 (2019), arXiv:1906.04354 [astro-ph.HE].
- [89] L. Lindblom and S. L. Detweiler, The quadrupole oscillations of neutron stars., Astrophys. J. 53 , 73 (1983).
- [90] S. Detweiler and L. Lindblom, On the nonradial pulsations of general relativistic stellar models, Astrophys. J. 292 , 12 (1985).
- [91] S. Chandrasekhar and V. Ferrari, On the non-radial oscillations of a star, Proc. Roy. Soc. Lond. A 432 , 247 (1991).
- [92] H. Sotani, K. Tominaga, and K.-i. Maeda, Density discontinuity of a neutron star and gravitational waves, Phys. Rev. D 65 , 024010 (2001).
- [93] M. Leins, H. P. Nollert, and M. H. Soffel, Nonradial oscillations of neutron stars: A New branch of strongly damped normal modes, Phys. Rev. D 48 , 3467 (1993).
- [94] F. J. Zerilli, Effective potential for even-parity reggewheeler gravitational perturbation equations, Phys. Rev. Lett. 24 , 737 (1970).
- [95] G. F. Burgio, V. Ferrari, L. Gualtieri, and H.-J. Schulze, Oscillations of hot, young neutron stars: Gravitational wave frequencies and damping times, Phys. Rev. D 84 , 044017 (2011).
- [96] V. Ferrari, G. Miniutti, and J. A. Pons, Gravitational waves from newly born, hot neutron stars, Monthly Notices of the Royal Astronomical Society 342 , 629 (2003), https://academic.oup.com/mnras/articlepdf/342/2/629/3433686/342-2-629.pdf.
- [97] H. Sotani, T. Kuroda, T. Takiwaki, and K. Kotake, Probing mass-radius relation of protoneutron stars from gravitational-wave asteroseismology, Phys. Rev. D 96 , 063005 (2017).
- [98] C. Drischler, K. Hebeler, and A. Schwenk, Chiral interactions up to next-to-next-to-next-to-leading order and nuclear saturation, Phys. Rev. Lett. 122 , 042501 (2019).
- [99] J. Keller, C. Wellenhofer, K. Hebeler, and A. Schwenk, Neutron matter at finite temperature based on chiral effective field theory interactions, Phys. Rev. C 103 , 055806 (2021).
- [100] J. Keller, K. Hebeler, and A. Schwenk, Nuclear equation of state for arbitrary proton fraction and temperature based on chiral effective field theory and a gaussian process emulator, Phys. Rev. Lett. 130 , 072701 (2023).
- [101] E. R. Most, L. R. Weih, L. Rezzolla, and J. SchaffnerBielich, New constraints on radii and tidal deformabilities of neutron stars from gw170817, Phys. Rev. Lett. 120 , 261103 (2018).
- [102] E. Annala, T. Gorda, A. Kurkela, and A. Vuorinen, Gravitational-wave constraints on the neutronstar-matter equation of state, Phys. Rev. Lett. 120 , 172703 (2018).
- [103] C. Hartnack, H. Oeschler, and J. Aichelin, Hadronic matter is soft, Phys. Rev. Lett. 96 , 012302 (2006).

- [104] A. Le F` evre, Y. Leifels, W. Reisdorf, J. Aichelin, and C. Hartnack, Constraining the nuclear matter equation of state around twice saturation density, Nuclear Physics A 945 , 112 (2016).
- [105] P. Russotto et al., Results of the asy-eos experiment at gsi: The symmetry energy at suprasaturation density, Phys. Rev. C 94 , 034608 (2016).
- [106] M. Prakash, I. Bombaci, M. Prakash, P. J. Ellis, J. M. Lattimer, and R. Knorren, Composition and structure of protoneutron stars, Physics Reports 280 , 1 (1997).
- [107] K. D. K. Nils Andersson, Towards gravitational wave asteroseismology, Monthly Notices of the Royal Astronomical Society 299 , 1059-1068 (1998).
- [108] O. Benhar, V. Ferrari, and L. Gualtieri, Gravitational wave asteroseismology reexamined, Phys. Rev. D 70 , 124015 (2004).
- [109] D. D. Doneva, E. Gaertig, K. D. Kokkotas, and C. Kr¨ger, Gravitational wave asteroseismology of fast u rotating neutron stars with realistic equations of state, Phys. Rev. D 88 , 044052 (2013).
- [110] C. A. Raithel and V. Paschalidis, Influence of stellar compactness on finite-temperature effects in neutron star merger simulations, Phys. Rev. D 108 , 083029 (2023).
- [111] O. E. Andersen, S. Zha, A. da Silva Schneider, A. Betranhandy, S. M. Couch, and E. P. O'Connor, Equationof-state dependence of gravitational waves in corecollapse supernovae, The Astrophysical Journal 923 , 201 (2021).
- [112] A. S. Schneider, C. Constantinou, B. Muccioli, and M. Prakash, Akmal-pandharipande-ravenhall equation of state for simulations of supernovae, neutron stars, and binary mergers, Phys. Rev. C 100 , 025803 (2019).
- [113] A. Raduta, F. Nacu, and M. Oertel, Equations of state for hot neutron stars, The European Physical Journal A 57 , 1 (2021).
- [114] S. Ghosh, S. Shaikh, P. J. Kalita, P. Routaray, B. Kumar, and B. K. Agrawal, Exploring the macroscopic properties and nonradial oscillations of proto-neutron stars: Effects of temperature, entropy, and lepton fraction (2024), arXiv:2307.06892 [nucl-th].
- [115] A. Kumar, P. Thakur, and M. Sinha, Non-radial oscillations in newly born compact star considering effects of phase transition, Monthly Notices of the Royal Astronomical Society 530 , 501-513 (2024).
- [116] G. Martinon, A. Maselli, L. Gualtieri, and V. Ferrari, Rotating protoneutron stars: Spin evolution, maximum mass, and i-love-q relations, Phys. Rev. D 90 , 064026 (2014).
- [117] S. S. Lenka, P. Char, and S. Banik, New covariant density functionals of nuclear matter for compact star simulations, J. Phys. G: Nucl. Part. Phys. 46 , 105201 (2019).
- [118] V. Guedes, S. Y. Lau, C. Chirenti, and K. Yagi, The broadening of universal relations at the birth and death of a neutron star (2024), arXiv:2402.10868 [astroph.HE].
- [119] B. Haskell, R. Ciolfi, F. Pannarale, and L. Rezzolla, On the universality of I-Love-Q relations in magnetized neutron stars, Monthly Notices of the Royal Astronomical Society: Letters 438 , L71 (2013).
- [120] T. G. Cowling, The Non-radial Oscillations of Polytropic Stars, Monthly Notices of the Royal Astronomical Society 101 , 367 (1941), https://academic.oup.com/mnras/articlepdf/101/8/367/8071901/mnras101-0367.pdf.
- [121] M. G. Alford and S. P. Harris, Beta equilibrium in neutron-star mergers, Physical Review C 98 , 10.1103/physrevc.98.065806 (2018).
- [122] M. G. Alford, A. Haber, S. P. Harris, and Z. Zhang, Beta equilibrium under neutron star merger conditions, Universe 7 , 399 (2021).
- [123] N. Barman, Repository: coupling constants for nucleonic rmf, (2024).
Nilaksha Barman Bikram Keshari Pradhan Debarati Chatterjee
Investigating the role of nuclear parameters on oscillation modes in hot Neutron Stars
Recent studies have revealed that certain nuclear parameters are more
dominant than others in governing global neutron star properties, such as its
structure or oscillation mode characteristics. Although neutron stars can in
general assumed to be cold, in astrophysical scenarios such as newly born
neutron stars or remnants of binary neutron star mergers, finite temperature
effects play a non-negligible role. In this work, we perform a consistent and
systematic investigation of the role of nuclear parameters and thermal effects
on neutron star properties and fluid oscillation modes within a full general
relativistic scheme. We impose constraints on the parameter space of the
relativistic mean field model using state-of-the-art information from
terrestrial experiments and multi-messenger astrophysical data. We find
effective nucleon mass to be the most important nuclear parameter controlling
astrophysical observables of hot neutron stars, similar to the cold beta
equilibrated matter. However, we conclude that the interplay among saturation
properties and astrophysical observables depends not only on the thermal
configurations considered but also on the constraints imposed. We also
investigated the role of nuclear saturation parameters on some universal
relations for hot NSs which are important in gravitational wave
asteroseismology. Our investigation confirmed that these relations are mostly
insensitive to nuclear saturation properties and mainly affected by variation
of charge fraction in the star.