## Universal relationships for neutron stars from perturbative approach

Debasis Atta 1 † , Vinay Singh 2 § and D. N. Basu 3 §

† Department of Higher Education,

Government of West Bengal, ∗ (Dated: August 2, 2024)

Bikash Bhavan, Saltlake, Kolkata 700091, India and

§ Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700064, India

The universal relationships for compact stars have been investigated employing perturbative approach using canonical (APR) and Brussels-Montreal Skyrme (BSk22, BSk24, BSk26) equations of state describing hadronic matter of neutron stars. The neutron star matter has been considered to be β -equilibrated neutron-proton-electron-muon matter at the core with a rigid crust. The multipole moments of a slowly rotating neutron star characterize its external gravitational field. These variables are dependent on the interior structure of the neutron star described by the equation of state of the neutron star matter. The properties of neutron stars, such as the mass, the radius, the dimensionless moment of inertia, the compactness, the Love number, the dimensionless tidal deformability and the dimensionless quadrupole moment have been calculated and relations among these quantities have been explored. It is found that most of these relations do not depend sensitively on the details of the internal structure of neutron stars. Such universality implies that the measurement of a single quantity appearing in a universal relation would automatically provide information about the others, notwithstanding the fact that those may not be accessible observationally. These can be used to estimate deformability of compact stars through moment of inertia measurements, to quantify spin in binary inspirals by breaking degeneracies in the detection of gravitational waves and test General Relativity in a way that is independent of nuclear structure.

Keywords : Neutron Star; Love number; Tidal deformability; Quadrupole moment.

PACS numbers: 21.65.-f, 26.60.-c, 04.30.-w, 26.60.Dd, 26.60.Gj, 97.60.Jd, 04.40.Dg, 21.30.Fe

## I. INTRODUCTION

The compact stars provide the most conspicuous natural laboratories for exploring the matter constituting these objects [1-4] at the extremes of density, pressure and isospin asymmetry. Gravitational waves (GWs) detected from recently found pulsars, like GW170817, have provided fresh insights into the internal structure of these stars, primarily regarding composition [5-7]. While the nuclear matter equation of state (EoS) is well known up to nuclear saturation density, characterizing matter at extremely high densities realized in the interior of these stars is a challenge. The nature of the matter at these densities being currently unknown, makes it difficult to create stellar models that would account for the observed findings. Probable candidates include weird quark stars made up of deconfined quarks, hybrid stars made up of deconfined quark cores with hadronic outer shells and pure neutron stars (NSs) made up of hadrons. In the present work, universal relationships for Nss have been investigated. The findings are extremely significant, particularly for the study of dense matter physics [8-11].

from observations was presented in Ref. [16]. Some of the recent works dedicated to the study of the NS interior and possible existence of twin stars can be found in Refs.[17-61].

The first time, a few instances of EoSs that fulfill the maximum mass constraint arisen from the existence of 2 M ⊙ pulsars [12, 13] were presented and were subsequently elaborated by Benic et al. [14, 15]. A systematic Bayesian analysis of the NS EoS with constraints deduced

∗ E-mail 1: debasisa906@gmail.com; E-mail 3: vsingh@vecc.gov.in; E-mail 3: dnb@vecc.gov.in

At a frequency of f = 709 Hz, pulsar PSR J09520607 was initially identified by Bassa et al. [62] as the fastest known spinning pulsar in the disk of Milky Way. Recently, Romani et al. [63] have discovered that the mass of PSR J0952-0607 is M = 2 35 . ± 0 17 . M ⊙ which is the precisely measured largest mass found till date. The aforementioned findings led to investigations into the limits of the evolution of dense nuclear matter and may even reexamine many theoretical predictions about fundamental characteristics of NSs [64] like the presence of twin stars. The twin stars, however, might not exist if there is a supermassive NS [23]. A methodical and meticulous examination ought to be contemplated, as the detection of these supermassive stars is crucial in verifying the reliability of the associated theoretical predictions. Conversely, recent discoveries about the radius and tidal deformability of stars with mass close to 1 4 . M ⊙ should be taken into account when combining theoretical models for the maximum mass of NSs. A nuclear model is often evaluated based on its ability to accurately estimate the observable radial and tidal deformability values while predicting the maximum achievable mass. Numerous studies have examined the rotating NSs with unusual degrees of freedom in their cores [38, 55, 58, 59].

In the present study, we employ the constraints on mass and radius of the recently observed compact objects, namely, the supernova remnant HESS J1731-

FIG. 1. The mass-radius plots for APR, BSk22, BSk24 and BSk26 EoSs. The shaded regions represent the HESS J1731347 remnant [65], the GW170817 event [66], PSR J1614-2230 [67], PSR J0348+0432 [68], PSR J0740+6620 [69], and PSR J0952-0607 [70] pulsar observations for the possible maximum mass.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-0e5a75d4-65a7-4f66-bc0e-75a3e38dc8be.jpg)

347 [65], the GW170817 event [66], PSR J1614-2230 [67], PSR J0348+0432 [68], PSR J0740+6620 [69] and PSR J0952-0607 [70] pulsars. The the supernova remnant HESS has M = 0 77 . +0 20 . -0 17 . M ⊙ and R = 10 4 . +0 86 . -0 78 . km, while the observations were made with the help of X-ray spectrum modeling and a reliable distance estimate from Gaia observations [65]. The properties of neutron stars, including mass ( M ), radius ( R ), compactness ( C ), Love number ( k 2 ), dimensionless tidal deformability (Λ), dimensionless quadrupole moment( Q ) and the dimensionless moment of inertia ( I ), have been calculated. It has been explored how these quantities are universally related.

The paper is arranged in the following manner. In Sec. II we provide basic formalism leading to the tidal deformability while rate of rotation of an inertial frame and quadrupole deformation are described briefly in Sec. III. In Sec. IV results of the present calculations are presented and discussed. Finally, summary of the present theoretical endeavor along with concluding remarks are furnished in Sec. V. In the present work we use units G = c = 1 throughout.

## II. TIDAL DEFORMABILITY

## A. Tolman-Oppenheimer-Volkoff equations

The hydrostatic equilibrium of compact star matter is well determined by a system of differential equations

FIG. 2. Fitting curve (black dashed), given in Eq.(25), and numerical results (continuous lines) of the universal relation for dimensionless quadrupole moment ( Q ) vs mass ( M ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-e8ef34b1-2f4c-41bf-9dd4-c236889b0eff.jpg)

known as the Tolman-Oppenheimer-Volkoff (TOV) equations [71, 72]. The structure of a slowly rotating, spherical star can, therefore, be well described by the TOV equations

$$& \frac { d M ( r ) } { d r } = 4 \pi r ^ { 2 } \varepsilon ( r ), & ( 1 ) \\ & \frac { d \nu ( r ) } { d r } = 2 \frac { [ M ( r ) + 4 \pi r ^ { 3 } P ( r ) ] } { r [ r - 2 M ( r ) ] }, \\ & \frac { d P ( r ) } { d r } = - \frac { [ P ( r ) + \varepsilon ( r ) ] } { 2 } \frac { d \nu ( r ) } { d r }. \\ ^ { \prime } \, \cdots \quad, \, \cdots \quad, \quad.. \quad.$$

where ε r ( ) and P r ( ) are the energy density and pressure distribution inside the compact star. Since all functions considered here depend on radial co-ordinate r , this dependency will be omitted henceforth. The quantities, R and M R ( ), hereafter will indicate, respectively, the radius and the mass of the compact star, obtained by solving Eqs.(1) for a designated EoS with the boundary condition that P r ( ) → 0 on the surface of the compact star at r = R .

## B. Love number and tidal deformability

The GW emission during the late phase of inspiraling of a binary NS system ahead of merger [73-75] is an important source for the GW detectors. In the inspiral phase the tidal effects can be detected [74]. In recent years, observations of GWs resulting from sources like

FIG. 3. Fitting curve (black dashed), given in Eq.(25), and numerical results (continuous lines) of the universal relation for Love number ( k 2 ) vs mass ( M ) for APR, BSk22, BSk24 and BSk26 EoSs. The shaded regions represent the HESS J1731-347 remnant [65], the GW170817 event [66], PSR J1614-2230 [67], PSR J0348+0432 [68], PSR J0740+6620 [69], and PSR J0952-0607 [70] pulsar observations for the possible maximum mass.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-cedccfd9-9e7b-4402-ad87-e74223361f8f.jpg)

the merger of black hole-NS and NS-NS binary systems induce measurements of various properties of NSs.

The quadrupole moment Q ij of compact star and the external tidal field E ij are defined as coefficients in an asymptotic expansion of the total metric at large distances r from the star. The tidal Love number k 2 depends upon EoS and characterizes the response of NS to the tidal field E ij [74]. This relation can be expressed as

$$Q _ { i j } = - \frac { 2 } { 3 } k _ { 2 } R ^ { 5 } \mathcal { E } _ { i j } = - \lambda \mathcal { E } _ { i j } \quad \quad ( 2 ) \quad \text{as } \mathcal { \ } \\ \text{and}$$

where R is radius of NS and λ = 2 k R / 2 5 3 is the tidal deformability. The l = 2 tidal Love number is given by [74, 75]

$$k _ { 2 } \, = \, \frac { 8 C ^ { 5 } } { 5 } ( 1 - 2 C ) ^ { 2 } [ 2 + 2 C ( y - 1 ) - y ] & \frac { d r } { d r } \ \left. \\ & \times \left \{ 2 C [ 4 ( 1 + y ) C ^ { 4 } + ( 6 y - 4 ) C ^ { 3 } + ( 2 6 - 2 2 y ) C ^ { 2 } \\ & + 3 ( 5 y - 8 ) C - 3 y + 6 ] + 3 ( 1 - 2 C ) ^ { 2 } \\ & \times [ 2 C ( y - 1 ) - y + 2 ] \ln ( 1 - 2 C ) \right \} ^ { - 1 }, \\ \text{where the dimensionless quantity $C$ = $M(R)/R$ is the solved}$$

where the dimensionless quantity C = M R /R ( ) is the compactness of the star and the dimensionless quantity y is defined as

FIG. 4. Fitting curve (black dashed), given in Eq.(25), and numerical results (continuous lines) of the universal relation for dimensionless tidal deformability (Λ) vs mass ( M ) for APR, BSk22, BSk24 and BSk26 EoSs. The shaded regions represent the HESS J1731-347 remnant [65], the GW170817 event [66], PSR J1614-2230 [67], PSR J0348+0432 [68], PSR J0740+6620 [69], and PSR J0952-0607 [70] pulsar observations for the possible maximum mass.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-f0083b37-c70c-4f57-b041-c3b75a90897d.jpg)

$$y = \frac { R \, \beta ( R ) } { H ( R ) }$$

for the internal solution determined by numerically solving the second-order differential equation for H r ( ). The second-order differential equation for H can be separated into a first-order system of ordinary differential equations in terms of the usual TOV quantities M p r , ( ), and ε p ( ), as well as the additional functions H r ( ), β r ( ) = dH/dr , and the EoS function f ( p ) = dε/dp :

$$\dots & \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dia \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dots \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dia \quad \dia \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dia \quad \d� \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \digma \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \d� \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dia \quad \digma \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \digma \quad \d� \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dag \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \d } \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \d  \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dangle \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \d | \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \dagger \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \d] \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \delta \quad \delta \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \quad \delta \$$

The above equations combined with TOV Eqs.(1) is then solved simultaneously. The system is integrated outward starting just outside the center using the expansions H r ( ) = a r 0 2 and β r ( ) = 2 a r 0 as r → 0. The constant

FIG. 5. Plots of moment of inertia I 45 in units of 10 45 g cm 2 vs M for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-9dbd4054-63c5-44d2-ba6d-60ece6bc733d.jpg)

FIG. 6. Plots of compactness ( C ) vs dimensionless moment of inertia ( I ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-6159b202-dada-4a90-8b05-ee16b4ae75fc.jpg)

a 0 determines how much the star is deformed and can be chosen arbitrarily as it cancels in the expression for the Love number. The boundary conditions that determine the unique choice of this solution follow from matching the interior and exterior solutions and their first derivatives at the boundary of the star, where r = R . The dimensionless tidal deformability related to the Love number can now be given by

FIG. 7. Plots of dimensionless moment of inertia ( I ) vs mass ( M ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-e8a356d7-266f-4add-80f5-1f343af2e127.jpg)

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

where C is the compactness of the star defined previously. One of the primary quantities that can be accurately determined by the detection of the associated GWs is the tidal deformability, or Λ.

## III. RATE ROTATION OF INERTIAL FRAME

The compact star has been assumed to be a slowly rotating object having a uniform angular velocity Ω [76, 77]. To the first order in Ω the shape of the star remains spherical and there is only one function, ω , to be determined, which is the angular velocity of the local inertial frame [76, 78, 79]. By introducing the function ϖ = Ω -ω , the angular velocity of a point in the star measured with respect to the angular velocity of the local inertial frame, it can be shown that it satisfies the following equation [78, 80]

$$\frac { 1 } { r ^ { 4 } } \frac { d } { d r } ( r ^ { 4 } j \frac { d \varpi } { d r } ) + \frac { 4 } { r } \frac { d j } { d r } \varpi = 0,$$

where j ( r ) = e -ν/ 2 √ 1 -2 M r which → 1 at r → R . The Eq.(8) can be solved subject to the boundary conditions that ϖ r ( ) is regular as r → 0 and ϖ r ( ) → Ω as r →∞ . The moment of inertia of the star can then be calculated dividing the total angular momentum J by Ω as

$$\mathcal { I } = \frac { 1 } { 6 \Omega } R ^ { 4 } \frac { d \varpi } { d r } \Big | _ { r = R }.$$

FIG. 8. Fitting curve (black dashed), given in Eq.(25), and numerical results (points) of the universal relation for dimensionless moment of inertia ( I ) vs dimensionless tidal deformability (Λ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-323e4768-8028-4a75-aa82-de7724e93b94.jpg)

FIG. 9. Relative fractional errors between the fitting curve and numerical results of Fig.-8.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-03ccccc8-d96e-40f3-8319-d7357c291d64.jpg)

In the dimensionless form moment of inertia I of the star is given by I = I /M 3 .

The Eq.(8) can be reduced to two first order equations by defining two auxiliary functions of r

$$\chi = j \varpi, \quad \ \ u = r ^ { 4 } j \frac { d \varpi } { d r } \quad \ \ ( 1 0 ) \quad \ \ t h e$$

FIG. 10. Fitting curve (black dashed), given in Eq.(25), and numerical results (points) of the universal relation for compactness ( C ) vs dimensionless tidal deformability (Λ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-c511e87f-72a4-492a-9c54-6a3784585158.jpg)

FIG. 11. Relative fractional errors between the fitting curve and numerical results of Fig.-10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-cb19d892-b59d-4539-ab96-2ee060e347a3.jpg)

that satisfy

$$r \leq R \quad & \frac { d \chi } { d r } = \frac { u } { r ^ { 4 } } - \frac { 4 \pi r ^ { 2 } ( \varepsilon + P ) \chi } { r - 2 M } \\ & \frac { d u } { d r } = \frac { 1 6 \pi r ^ { 5 } ( \varepsilon + P ) \chi } { r - 2 M }. \quad & ( 1 1 ) \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

The above equations have to be integrated outward from the center of the star r = 0 to its surface at r = R

FIG. 12. Fitting curve (black dashed), given in Eq.(25), and numerical results (points) of the universal relation for dimensionless moment of inertia ( I ) vs dimensionless quadrupole moment ( Q ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-ec2edbf5-ba7e-4472-b000-711fda416c1a.jpg)

FIG. 13. Relative fractional errors between the fitting curve and numerical results of Fig.-12.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-4f3ab84f-0a11-432a-80c9-bad7f6668c20.jpg)

imposing ensuing boundary (initial) conditions that

$$\chi ( 0 ) = j ( 0 ) \varpi _ { c }, \quad u$$

These initial conditions result from the behavior of ϖ near r → 0

$$\varpi \sim \varpi _ { c } ( 1 + \varpi _$$

FIG. 14. Fitting curve (black dashed), given in Eq.(25), and numerical results (points) of the universal relation for dimensionless quadrupole moment ( Q ) vs dimensionless tidal deformability (Λ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-1a159e3c-7afb-4c6d-a96d-7b83094ee64a.jpg)

FIG. 15. Relative fractional errors between the fitting curve and numerical results of Fig.-14.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-2f0c1cdd-e22e-4c29-a3e9-29c404d55ab2.jpg)

where ϖ c is a constant. The constant ϖ c can be chosen arbitrarily. When one reaches the surface and only then one can determine the angular velocity Ω and angular momentum J , corresponding to ϖ c . It is worthwhile to mention here that from Eq.(10) and Eq.(13) it follows that the entity u goes to zero more rapidly than r 4 .

For r ≥ R P , = 0, ε = 0, M ≡ M R ( ), j ( R ) = 1 and the solution of Eqs.(11) can be given by

$$r \geq R \quad \chi ( r ) = \Omega - \frac { 2 J$$

where J is a constant of motion that, to first order in Ω, reflects the star's angular momentum. The requirement of continuity of the exterior and interior solutions demands that they be matched at R , which yields the values of Ω and J . This indicates that J = u R / ( ) 6 and Ω = χ R ( ) + 2 J/R 3 . If a different value of Ω is desired, one rescales the function ϖ r ( ) through ϖ c to obtain it.

## A. Quadrupole deformation

The equations for h 2 ( r ) and v 2 ( r ) are given by

FIG. 16. Fitting curve (black dashed), given in Eq.(25), and numerical results (points) of the universal relation for compactness ( C ) vs dimensionless quadrupole moment ( Q ) for APR, BSk22, BSk24 and BSk26 EoSs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-14208c8d-ae52-44a9-9934-a78d6e0e2173.jpg)

$$\begin{array} {$$

which are responsible for the quadrupole ( l = 2) deformation of a compact star. In order to solve these equations for r ≤ R , it is important to note that a regular solution of Eqs.(15) at r → 0 must behave as

$$r \to 0 \quad h _ { 2 } ( r ) \sim A r ^ { 2 }$$

where the equation that describes how the constants A and B are connected is given by

$$B + 2 \pi \left [ \frac { 1 } { 3 } \varepsilon ( 0 ) + P ( 0 ) \right ] A = \frac { 2 } { 3 } \pi \left [ \varepsilon ( 0 ) + P ( 0 ) \right ] \left ( j ( 0 ) \varpi _ { c } \right ) ^ { 2 }. \\ \pi _ { c } \pi _ { \dots } \dots \dots \pi _ { \dots } \varepsilon ( \pi _ { \dots } \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \circ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots.$$

Following Refs.[78, 80], the general solution of Eqs.(15) for r ≤ R can be written as

$$h _ { 2 } ( r ) & = h _ { 2 } ^ { P } + C _ { 2 } \ h _ { 2 } ^ { H } \quad \quad ( 1 8 ) \\ v _ { 2 } ( r ) & = v _ { 2 } ^ { P } + C _ { 2 } \ v _ { 2 } ^ { H }.$$

The constant C 2 remains to be determined while ( h P 2 , v P 2 ) are a particular solution which can be found by integrating Eqs.(15) with the initial conditions, say, A = 1 and B given by Eq.(17) and ( h H 2 , v H 2 ) are the solution of the homogeneous equations which can be obtained by setting χ = 0 and u = 0 in Eqs.(15) and integrating numerically with the following boundary conditions

FIG. 17. Relative fractional errors between the fitting curve and numerical results of Fig.-16.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-f6354782-0407-4227-9caf-c168a9a379ab.jpg)

$$r \to 0 \quad & h _ { 2 } ^ { H } ( r ) \sim r ^ { 2 } \quad & ( 1 9 ) \\ & v _ { 2 } ^ { H } ( r ) \sim - 2 \pi \ \left [ \frac { 1 } { 3 } \varepsilon ( 0 ) + P ( 0 ) \right ] r ^ { 4 }.$$

For r ≥ R , analytic solutions for Eqs.(15) in terms of the associated Legendre functions of second kind can be deduced [76, 78] which are

FIG. 18. Plots of chirp mass versus effective tidal deformability (Λ ∗ ) for binary mass ratio q = 0 7. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-3a73719d-c30a-4abd-a753-19f1706a80ce.jpg)

FIG. 19. Plots of chirp mass versus effective tidal deformability (Λ ∗ ) for binary mass ratio q = 0 9. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-79477660-d2d9-4255-a991-abaf77faaf7a.jpg)

$$r \geq R \quad h _ { 2 } ( r ) & = J ^ { 2 } ( \frac { 1 } { M ( R ) r ^ { 3 } } + \frac { 1 } { r ^ { 4 } } ) + K \ Q _ { 2 } ^ { 2 } ( \xi ) \quad ( 2 0 ) \\ v _ { 2 } ( r ) & = - \frac { J ^ { 2 } } { r ^ { 4 } } + K \ \frac { 2 M ( R ) } { [ r ( r - 2 M ( R ) ) ] ^ { 1 / 2 } } Q _ { 2 } ^ { 1 } ( \xi )$$

where K is a constant and

$$Q _ { 2 } ^ { 2 } ( \xi ) = [ \frac { 3 } { 2 } ( \xi ^ { 2 } - 1 ) \log ( \frac { \xi + 1 } { \xi - 1 } ) - \frac { 3 \xi ^ { 3 } - 5 \xi } { \xi ^ { 2 } - 1 } ] \quad ( 2 1 ) \quad \overset { \cdot \cdot } { \text{vis} }$$

FIG. 20. Plots of chirp mass versus effective tidal deformability (Λ ∗ ) for binary mass ratio q = 0 99. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-f3d8d800-88e6-42b4-b2ba-fc479c93cea2.jpg)

FIG. 21. Plots of chirp mass versus effective tidal deformability (Λ ∗ ) for binary mass ratio q = 1 0. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00646v1-figure-969cbc72-b90d-4918-ac63-0c87161267eb.jpg)

$$Q _ { 1 } ^ { 2 } ( \xi ) = ( \xi ^ { 2 } - 1 ) ^ { 1 / 2 } [ \frac { 3 \xi ^ { 2 } - 2 } { \xi ^ { 2 } - 1 } - \frac { 3 } { 2 } \xi \log ( \frac { \xi + 1 } { \xi - 1 } ) ].$$

with

$$\xi = \frac { r } { M ( R ) } - 1. \text{ \quad \ \ } ( 2 2 )$$

The constants C 2 and K can now be evaluated by imposing the condition that the solutions (18) and the solutions (20) are continuous at r = R and hence equal.

Therefore, the functions h 2 and v 2 are, thus, completely determined.

## B. Quadrupole moment

Once the constants C 2 and K have been determined by imposing the continuity condition that the solutions (18) and the solutions (20) be matched and equated at r = R , the quadrupole moment Q can be determined in terms of K as

$$\mathcal { Q } = \frac { J ^ { 2 } } { M } + \frac { 8 } { 5 } K M ^ { 3 } \quad \quad ( 2 3 ) \quad \text{ w}$$

which can be expressed in the dimensionless form as

$$Q = \frac { \mathcal { Q } } { M ^ { 3 } ( J / M ^ { 2 } ) ^ { 2 } } = 1 + \frac { 8 } { 5 } K \frac { M ^ { 4 } } { J ^ { 2 } } \, \quad \ \ ( 2 4 )$$

## IV. THEORETICAL CALCULATIONS AND RESULTS

In this work, EoSs obtained from the effective nucleon-nucleon (NN) interaction described by APR [81] and Brussels-Montreal Skyrme effective interaction with BSk22, BSk24 and BSk26 parameter sets [82] have been used. The β -equilibrated npeµ neutron star matter has been used for subsequent calculations. The calculations for masses and radii have been performed using FMT [83], BPS [84] and BBP [85] the EoSs up to the number density of 0.0582 fm -3 for the outer crust and the β -equilibrated NS matter beyond covering the inner crustal and core regions of a compact star. The location of the inner edge of the NS crust, the core-crust transition density and pressure were determined [86] by dynamical method. The observations of the binary millisecond pulsar J16142230 suggest that its mass lies in the range 1 97 . ± 0 04 . M ⊙ [67, 87]. The measurements of radio timing for pulsar PSR J0348+0432 and its companion (white dwarf) have confirmed the mass of the pulsar to be in the range of 1.97 -2.18 M ⊙ [68]. Very recently, the studies for PSR J0740+6620 [88] and for PSR J0952-0607 [70] find masses of 2.08 ± 0.07 M ⊙ and 2.35 ± 0.17 M ⊙ , respectively. Current observations of PSR J0740+6620 also suggest its mass to be 2.072 +0 067 . -0 066 . M ⊙ [89, 90]. The mass versus radius of NSs obtained using APR, BSk22, BSk24 and BSk26 EoSs have been plotted in Fig.-1. The maximum masses obtained using APR, BSk22, BSk24 and BSk26 are, repectively, 2.19 M ⊙ , 2.27 M ⊙ , 2.28 M ⊙ and 2.18 M ⊙ [86] which provide reasonably good estimates for latest observations of NS masses. The shaded regions represent the HESS J1731-347 remnant [65], the GW170817 event [66], PSR J1614-2230 [67], PSR J0348+0432 [68], PSR J0740+6620 [69], and PSR J0952-0607 [70] pulsar observations for the possible maximum mass.

## A. Universal relationships

The universal relations among M I , , Λ, C and Q have been investigated for APR, BSk22, BSk24 and BSk26 EoSs. It has been observed that these relations hold universally for NS sequences, essentially independently of their EoS. Such relations can be numerically fitted with a polynomial given by

$$Y _ { f i t } = \sum _ { n = 0 } ^ { 4 } a _ { n } X ^ { n } \, \text{ \quad \ \ } \text{(2.5)}$$

where the variables and fitted coefficients have been summarized in Table-I. The relative fractional error ∆ Y between the fitted values Y fit and calculated numerical results Y has been defined as

$$\Delta Y = \frac { | Y - Y _ { f i t } | } { Y _ { f i t } }$$

which provides the fractional difference between the numerical results and the corresponding analytic fits.

The variation of dimensionless quadrupole moment ( Q ) with NS mass ( M ) obtained for APR, BSk22, BSk24 and BSk26 EoSs are displayed in Fig.-2. The fitting curve (black dashed), given by Eq.(25), for the universal relation is also shown. In Figs.-3 and 4, respectively, variations of Love number ( k 2 ) and the dimensionless tidal deformability (Λ) with mass ( M ) are displayed for APR, BSk22, BSk24 and BSk26 EoSs. The fitting curves (black dashed), given by Eq.(25), for the universal relation for these are also shown. The shaded regions represent the HESS J1731-347 remnant [65], the GW170817 event [66], PSR J1614-2230 [67], PSR J0348+0432 [68], PSR J0740+6620 [69], and PSR J0952-0607 [70] pulsar observations for the possible maximum mass.

The moment of inertia I 45 in units of 10 45 g cm 2 versus mass M for APR, BSk22, BSk24 and BSk26 EoSs are plotted in Fig.-5. In Fig.-6 the plots of compactness ( C ) against dimensionless moment of inertia ( I ) for APR, BSk22, BSk24 and BSk26 EoSs corresponding to maximum masses are shown. In Fig.-7 the variation of dimensionless moment of inertia ( I ) with mass ( M ) for APR, BSk22, BSk24 and BSk26 EoSs are shown.

In Figs.-8,9,10,11,12,13,14,15,16,17, respectively, the plots of the dimensionless moment of inertia ( I ) versus dimensionless tidal deformability (Λ), compactness ( C ) versus the dimensionless tidal deformability (Λ), the dimensionless moment of inertia ( I ) versus dimensionless quadrupole moment ( Q ), dimensionless quadrupole moment ( Q ) versus the dimensionless tidal deformability (Λ) and compactness ( C ) versus dimensionless quadrupole moment ( Q ) are shown for APR, BSk22, BSk24 and BSk26 EoSs corresponding to maximum masses along with respective relative fractional errors between the fitting curve and results of numerical calculations. The fitting curves (black dashed), given by Eq.(25), for the

TABLE I. The variables and fitted coefficients.

| Y    | X    |       a 0 |        a 1 |         a 2 |          a 3 |          a 4 |
|------|------|-----------|------------|-------------|--------------|--------------|
| ln Q | ln M | 2.32466   | -1.66851   | -0.70994    | -0.15927     |  0.0510313   |
| k 2  | ln M | 0.0970581 | -0.0214968 | -0.0773239  | -0.0224548   | -0.000390682 |
| lnΛ  | ln M | 8.03693   | -5.39839   | -2.3464     | -1.06376     |  0.00729592  |
| ln I | lnΛ  | 1.58383   |  0.113245  |  0.0183125  | -0.00078812  |  3.31898e-05 |
| C    | lnΛ  | 0.360812  | -0.0357859 |  0.00031568 |  6.15497e-05 | -1.62488e-06 |
| ln I | ln Q | 1.33849   |  1.20787   | -0.479594   |  0.174713    | -0.0139314   |
| ln Q | lnΛ  | 0.184514  |  0.107931  |  0.043545   | -0.00361704  |  9.42312e-05 |
| C    | ln Q | 0.354018  | -0.117589  |  0.00197844 |  0.00261095  | -0.000217161 |

universal relation for these are also shown. It may be discerned that the universal relations strongly emerge except in cases of variations of Q k , 2 and Λ with NS mass M where the universal relations do not crop up so strongly.

## B. Chirp Mass

The chirp mass M c of a binary NSs system, having constituent masses M 1 and M 2 , is an reliably measurable quantity by the ground-based detectors [7] for GWs. It is defined in terms of constituent masses as

$$M _ { c } = \frac { ( \mathcal { M } _ { 1 } \mathcal { M } _ { 2 } ) ^ { 3 / 5 } } { ( \mathcal { M } _ { 1 } + \mathcal { M } _ { 2 } ) ^ { 1 / 5 } } = \mathcal { M } _ { 1 } \frac { q ^ { 3 / 5 } } { ( 1 + q ) ^ { 1 / 5 } }, \quad ( 2 7 ) \quad \underset { \text{nary} } { \text{serv} }$$

where by convention M ≥M 1 2 has been chosen so that the binary mass ratio q = M M 2 / 1 lies within the range 0 &lt; q ≤ 1.

( R ), compactness ( C ), Love number ( k 2 ), dimensionless tidal deformability (Λ), dimensionless quadrupole moment( Q ) and the dimensionless moment of inertia ( I ) within a perturbative approach using canonical (APR) and Brussels-Montreal Skyrme (BSk22, BSk24, BSk26) EoSs describing hadronic matter of neutron stars. Various relations among these calculated quantities have been studied. Most of these relations have been found not to depend sensitively on the structure of neutron star interior implying universality. According to this concept of universality, measurements of one quantity that occurs in a universal relation inevitably provide information about the others, even though those others might not be observable. These can be used to test General Relativity independently of nuclear structure, quantify spin in binary inspirals by breaking degeneracies in gravitational wave detection and assess the deformability of compact stars by moment of inertia measurements.

Another reliably measurable quantity is the dimensionless effective tidal deformability Λ ∗ which is defined as [7]

$$\Lambda ^ { * } & = \frac { 1 6 } { 1 3 } \frac { ( \mathcal { M } _ { 1 } + 1 2 \mathcal { M } _ { 2 } ) \mathcal { M } _ { 1 } ^ { 4 } \Lambda _ { 1 } + ( \mathcal { M } _ { 2 } + 1 2 \mathcal { M } _ { 1 } ) \mathcal { M } _ { 2 } ^ { 4 } \Lambda _ { 2 } } { ( \mathcal { M } _ { 1 } + \mathcal { M } _ { 2 } ) ^ { 5 } }, \quad \text{$\text{$\text{$\text{$eval$rel}$} \\ \text{$\text{$\text{$\text{$speed$} } \times \text{$\text{$\text{$\text{$speed$} } \times \text{$\text{$\text{$speed$} } \times \text{$\text{$\text{$speed$} } \times \text{$\text{$\text{$speed$} } \times \text{$\text{$\text{$speed$} } \times \text{$\text{$\text{$speed$} } \times \text{$\text{$speed$} } \times \text{$\text{$speed$} } \times \text{$\text{$speed$} } \times \text{$\text{$speed$} } \times \text{$\text{$speed$} }.$$

such that Λ ∗ = Λ 1 = Λ 2 when M 1 = M 2 where Λ 1 and Λ 2 are the dimensionless tidal deformabilities [7] for the corresponding masses. The variations of chirp mass (M ) as c a function of effective tidal deformability (Λ ) ∗ have been shown for binary mass ratios of q = 0 7 0 9 0 99 1 0, respectively, in Figs.-18,19,20,21. . , . , . , .

## V. SUMMARY AND CONCLUSIONS

Our objective in the present work is to explore the properties of neutron stars such as mass ( M ), radius

Exciting applications in fundamental physics, GW theory and astrophysics can be made possible by the universal relations. Numerous follow-up investigations can be possible by the examining the universal relations provided here. One could find out, for instance, if these relations hold for compact stars that rotate rapidly, have significant internal magnetic fields and anisotropic pressure. Perhaps the most severe approximation is the slowrotation approximation which is appropriate only for NSs with periods higher than or equal to one millisecond. It may be anticipated that the EoS universality identified here would withstand even for rapid rotation, but with distinct universal relations for stars with different spin periods. But this variation should be limited within 10% [80, 91-95]. To allow for fast spinning NSs, one possible expansion of this work would be to improve the universal relations.

## ACKNOWLEDGMENTS

D. N. Basu acknowledges support from Science and Engineering Research Board, Department of Science and

- [1] N. K. Glendenning, Compact stars -Nuclear physics, Particle physics and General Relativity, Springer (1997).
- [2] P. Haensel, A. Y. Potekhin, and D. G. Yakovlev, Neutron Stars 1: Equation of State and Structure (SpringerVerlag, New York, 2007).
- [3] J. Schaffner-Bielich, Compact Star Physics (Cambridge University Press, Cambridge, England,2020).
- [4] ] F. Weber, Pulsars as Astrophysical Laboratories for Nuclear and Particle Physics (Institute of Physics, Bristol, U.K., 1999).
- [5] B.P. Abbott et al. , Phys. Rev. Lett. 119 , 161101 (2017).
- [6] B. P. Abbott et al. , Phys. Rev. Lett. 121 , 161101 (2018).
- [7] B. P. Abbott et al. , Phys. Rev. X 9 , 011001 (2019).
- [8] G. Baym, T. Hatsuda, T. Kojo, P.D. Powell, Y. Song, and T. Takatsuka, Rep. Prog. Phys. 81 056902 (2018).
- [9] M.G. Alford, S. Han, and M. Prakash, Phys. Rev. D 88 , 083013 (2013).
- [10] J.E. Christian, A. Zacchi, and J. Schaffner-Bielich, Eur. Phys. J. A 54 , 28 (2018).
- [11] G. Montana, L. Tolos, M. Hanauske, L. Rezzolla, Phys. Rev. D 99 , 103009 (2019).
- [12] D. Blaschke, D.E. Alvarez-Castillo, and S. Benic, PoS CPOD2013 063 (2013); [arXiv:1310.3803].
- [13] D.E. Alvarez-Castillo and D. Blaschke, arXiv:1304.7758.
- [14] S. Benic, D. Blaschke, D.E. Alvarez-Castillo, T. Fischer, and S. Typel, Astron. Astrophys. 577 , A40 (2015).
- [15] D.E. Alvarez-Castillo and D. Blaschke, Phys. Part. Nucl. 46 , 846 (2015).
- [16] D.E. Alvarez-Castillo, A. Ayriyan, S. Benic, D. Blaschke. H. Grigorian, and S. Typel, Eur. Phys. J. A 52 , 69(2016).
- [17] A. Ayriyan, N.-U. Bastian, D. Blaschke, H. Grigorian, K. Maslov, and D. N. Voskresensky, Phys. Rev. C 97 , 045802 (2018).
- [18] K. Maslov, N. Yasutake, D. Blaschke, A. Ayriyan, H. Grigorian, and T. Maruyama, Phys. Rev. C 100 , 025802 (2019).
- [19] D. E. Alvarez-Castillo and D. B. Blaschke, Phys. Rev. C 96 , 045809, (2017)
- [20] M. Bejger, D. Blaschke, P. Haensel, J.L. Zdunik, and M. Fortin, Astron. Astrophys. 600 , A39 (2017).
- [21] A. Bhattacharyya, I. N. Mishustin, W. Greiner, J. Phys. G 37 , 025201 (2010).
- [22] J.E. Christian, A. Zacchi, and J. Schaffner-Bielich, Phys. Rev. D 99 , 023009 (2019).
- [23] J.E. Christian, and J. Schaffner-Bielich, Phys. Rev. D 103 , 063042. (2021)
- [24] J.E. Christian, and J. Schaffner-Bielich, Astrophys. J. 935 , 122 (2022).
- [25] S. Han and A.W. Steiner, Phys. Rev. D 99 , 083014 (2019).
- [26] P. L. Espino and V. Paschalidis, Phys. Rev. D 105 , 043014 (2022).
- [27] H. Tan, T. Dore, V. Dexheimer, J. Noronha-Hostler, and N. Yunes, Phys. Rev. D 105 , 023018 (2022).
- [28] J. J. Li, A. Sedrakian, and M. Alford, Phys. Rev. D 104 ,

Technology, Government of India, through Grant No. CRG/2021/007333.

L121302 (2021).

- [29] Z. Sharifi, M. Bigdeli, and D. Alvarez-Castillo, Phys. Rev. D 103 , 103011 (2021).
- [30] E. Benitez, J. Weller, V. Guedes, C. Chirenti, and M. Coleman Miller, Phys. Rev. D 103 , 023007 (2021).
- [31] J.J. Li, A. Sedrakian, and M. Alford, Phys. Rev. D 101 , 063022 (2020).
- [32] D.E. Alvarez-Castillo, D.B. Blaschke, A.G. Grunfeld, and V.P. Pagura, Phys. Rev. D 99 , 063010 (2019).
- [33] V. Paschalidis, K. Yagi, D. Alvarez-Castillo, D.B. Blaschke, and A. Sedrakian, Phys. Rev. D 97 , 084038 (2018).
- [34] I.F. Ranea-Sandoval, M.G. Orsaria, S. Han, F. Weber, and W.M. Spinella, Phys. Rev. C 96 , 065807 (2017).
- [35] M. Alford and A. Sedrakian, Phys. Rev. Lett. 119 , 161104 (2017).
- [36] A. Zacchi, L. Tolos, and J. Schaffner-Bielich, Phys. Rev. D 95 , 103008 (2017).
- [37] A. Zacchi, M. Hanauske, and J. Schaffner-Bielich, Phys. Rev. D 93 , 065011 (2016).
- [38] A. Bhattacharyya, S. K. Ghosh, M. Hanauske, and S. Raha, Phys. Rev. C 71 , 048801 (2005).
- [39] D. Sen, N. Alam, and G. Chaudhuri, Phys. Rev. D 106 , 083008 (2022).
- [40] T. Minamikawa, T. Kojo, and M. Harada, Phys. Rev. C 103 , 045205 (2021).
- [41] R. De Pietri, A. Drago, A. Feo, G. Pagliara, M. Pasquali, S. Traversi, and G. Wiktorowicz, Astrophys. J. 881 , 122 (2019).
- [42] J.L. Zdunik and P. Haensel, Astron. Astrophys. 551 , A61 (2013).
- [43] S. Han and M. Prakash, Astrophys. J. 899 , 164 (2020).
- [44] A. Li, Z.Y. Zhu, E.P. Zhou, J.M. Dong, J.N. Hu, C.J. Xia, Journ. High Energ. Astroph. 28 , 19 (2020).
- [45] N.K. Largani, T. Fischer, A. Sedrakian, M. Cierniak, D.E. Alvarez-Castillo, and D. B. Blaschke, Mon. Not. Roy. Astron. Soc. 515 , 3539 (2022).
- [46] O. Ivanytskyi and D. Blaschke, Phys. Rev. D 105 , 114042 (2022).
- [47] G.A. Contrera, D. Blaschke, J.P. Carlomagno, A.G. Grunfeld, and S. Liebing, Phys. Rev. C 105 , 045808 (2022).
- [48] S. Schramm, V. Dexheimer, and R. Negreiros, Eur. Phys. J. A 52 , 14 (2016).
- [49] G.F. Burgio, A. Drago, G. Pagliara, H.J. Schulze, and J.B. Wei, Astrophys. J. 860 , 139 (2018).
- [50] M. Sieniawska, W. Turczanski, M. Bejger, and J.L. Zdunik, Astron. Astrophys. 622 , A174 (2019).
- [51] E.R. Most, L.R. Weih, L. Rezzolla, and J. SchaffnerBielich, Phys. Rev. Lett. 120 , 261103 (2018).
- [52] R. Nandi and P. Char, Astrophys. J. 857 , 12 (2018).
- [53] T. Deloudis, P. Koliogiannis, Ch. Moustakidis, EPJ Web Conf. 252 , 06001 (2021).
- [54] Q.W. Wang, C. Shi, Y. Yan, and H.S. Zong, Nucl. Phys. A 1025 , 122489 (2022).

- [55] S. Banik, M. Hanauske, D. Bandyopadhyay, and W. Greiner, Phys. Rev. D 70 , 123004 (2004).
- [56] S. Banik and D. Bandyopadhyay, Phys. Rev. D 64 , 055805 (2001).
- [57] S. Banik and D. Bandyopadhyay, Phys. Rev. D 67 , 123003 (2003).
- [58] P. Haensel, M. Bejger, M. Fortin, and L. Zdunik, Eur. Phys. J. A 52 , 59 (2016).
- [59] G. Bozzola, P.L. Espino, C.D. Lewin, and V. Paschalidis, Europ. Phys. J. A 55 , 149 (2019).
- [60] J.P. Pereira, M. Bejger, N. Andersson, and F. Gittins, Astrophys. J. 895 , 28 (2020).
- [61] S. Han, M.A.A. Mamun, S. Lalit, C. Constantinou, M. Prakash, Phys. Rev. D 100 , 103022 (2019).
- [62] C.G. Bassa, Z. Pleunis, J.W.T. Hessels, et al. , Astrophys. J. Lett. 846 , L20 (2017).
- [63] R.G. Romani, D. Kandel, A.V. Filippenko, T.G. Brink, W. Zheng, Astrophys. J. Lett. 934 , L17 (2022).
- [64] C. Ecker and L. Rezzolla, arXiv:2209.08101.
- [65] V. Doroshenko, V. Suleimanov, G. P¨hlhofer and Andrea u Santangelo, Nat. Astron. 6 , 1444 (2022).
- [66] B. P. Abbott et al., Phys. Rev. X 9 , 011001 (2019).
- [67] Z. Arzoumanian, A. Brazier, S. Burke-Spolaor et al., Astrophys. J. Suppl. Ser. 235 , 37 (2018).
- [68] J. Antoniadis, P. Freire, N. Wex et al., Science 340 , 448 (2013).
- [69] H. Cromartie, E. Fonseca, S. Ransom et al., Nat. Astron. 4 , 72 (2020).
- [70] R. G. Romani, D. Kandel, A. V. Filippenko, T. G. Brink, and W. Zheng, Astrophys. J. Lett. 934 , L17 (2022).
- [71] R. C. Tolman, Phys. Rev. 55 , 364 (1939).
- [72] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55 , 374 (1939).
- [73] S. Postnikov, M. Prakash, and J. M. Lattimer, Phys. Rev. D 82 , 024016 (2010).
- [74] E. E. Flanagan and T. Hinderer, Phys. Rev. D 77 , 021502(R) (2008).
- [75] T. Hinderer, Astrophys. J. 677 , 1216 (2008).
- [76] J. B. Hartle, Astrophys. J. 150 , 1005 (1967).
- [77] W D Arnett and R L Bowers, Astrophys. J. Suppl. 33 , 415 (1977).
- [78] J. B. Hartle and K. S. Thorne, Astrophys. J. 153 , 807 (1968).
- [79] J. B. Hartle, Astrophys. Space Sci. 24 , 385 (1973).
- [80] Omar Benhar, Valeria Ferrari, Leonardo Gualtieri and Stefania Marassi1, Phys. Rev. D 72 , 044028 (2005).
- [81] A. Akmal, V.R. Pandharipande and D.G. Ravenhall, Phys. Rev. C 58 , 1804 (1998).
- [82] S. Goriely, N. Chamel, and J. M. Pearson, Phys. Rev. C 88 , 024308 (2013).
- [83] R P Feynman, N Metropolis and E Teller, Phys. Rev. 75 , 1561 (1949).
- [84] GBaym, C Pethick and P Sutherland, Astrophys. J. 170 , 299 (1971).
- [85] GBaym, H A Bethe and C J Pethick, Nucl. Phys. A 175 , 225 (1971).
- [86] Joydev Lahiri, Debasis Atta and D. N. Basu, Pramana J. Phys. Rev. 98 , 17 (2024).
- [87] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts and J. W. T. Hessels, Nature 467 , 1081 (2010).
- [88] E. Fonseca et al., Astrophys. J. Lett. 915 , L12 (2021).
- [89] Isaac Legred, Katerina Chatziioannou, Reed Essick, Sophia Han and Philippe Landrya, Phys. Rev. D 104 , 063003 (2021).
- [90] Thomas E. Riley et al., Astrophys. J. Lett. 918 , L27 (2021).
- [91] E. Berti and N. Stergioulas, Mon. Not. R. Astron. Soc. 350, 1416 (2004).
- [92] E. Berti, F. White, A. Maniopoulou, and M. Bruni, Mon. Not. R. Astron. Soc. 358 , 923 (2005).
- [93] G. Pappas and T. A. Apostolatos, Phys. Rev. Lett. 108 , 231104 (2012).
- [94] Kent Yagi and Nicola´ s Yunes, Phys. Rev. D 88 , 023009 (2013).
- [95] Kent Yagi and Nicola´ s Yunes, Phys. Rev. D 91 , 123008 (2015).
Debasis Atta Vinay Singh D. N. Basu
Universal relationships for neutron stars from perturbative approach
The universal relationships for compact stars have been investigated
employing perturbative approach using canonical (APR) and Brussels-Montreal
Skyrme (BSk22, BSk24, BSk26) equations of state describing hadronic matter of
neutron stars. The neutron star matter has been considered to be
$\beta$-equilibrated neutron-proton-electron-muon matter at the core with a
rigid crust. The multipole moments of a slowly rotating neutron star
characterize its external gravitational field. These variables are dependent on
the interior structure of the neutron star described by the equation of state
of the neutron star matter. The properties of neutron stars, such as the mass,
the radius, the dimensionless moment of inertia, the compactness, the Love
number, the dimensionless tidal deformability and the dimensionless quadrupole
moment have been calculated and relations among these quantities have been
explored. It is found that most of these relations do not depend sensitively on
the details of the internal structure of neutron stars. Such universality
implies that the measurement of a single quantity appearing in a universal
relation would automatically provide information about the others,
notwithstanding the fact that those may not be accessible observationally.
These can be used to estimate deformability of compact stars through moment of
inertia measurements, to quantify spin in binary inspirals by breaking
degeneracies in the detection of gravitational waves and test General
Relativity in a way that is independent of nuclear structure.