Phase transition and stiffer core fluid in neutron stars: effects on stellar configurations, dynamical stability, and tidal deformability

In this work, we investigate the influence of the phase transition and a stiffer fluid in neutron stars’ cores on the static equilibrium configuration, dynamical stability, and tidal deformability. For this aim, it is taken into account that the fluid in the core and the envelope follow the relativistic polytropic equation of state. We find that the phase transition and a stiffer fluid in the core will reflect in the total mass, radius, speed of sound, core radius, radial stability with a slow and rapid conversion at the interface, and tidal deformability. We also investigate the dimensionless tidal deformability Λ1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda _1$$\end{document} and Λ2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda _2$$\end{document} for a binary neutron stars system with chirp mass equal to GW170817. Finally, we contrast our results with observational data to show the role that phase transition and a stiffer core fluid could play in the study of neutron stars.


Introduction
The direct multimessenger detection from binary black holes merger carried out by the LIGO-Virgo scientific network [1][2][3][4][5] has marked the starting of the era of Gravitational Waves (GWs) astronomy.The detection of GWs has opened a new window to explore the cosmos and supplied some astrophysics and fundamental physics implications (check, e.g., [6][7][8]).Another important event of GWs comes from a merger of a pair of neutron stars (NSs) [9], known as event GW170817, which was also reported by the LIGO-Virgo scientific network.This new signal opened the GWs multimessenger astronomy, being the first detection with electromagnetic counterpart [10], and providing a set of valuable information about the properties of NSs and their equation of state (EOS).After the first detection of GWs from the NSs binary system, many important efforts have been realized to constrain, e.g., their radii and EOS [11][12][13][14][15]. Addia e-mail: jose.arbanil@upn.petional constraints are feasible because of the implication of tidal deformation [16][17][18].
It is known that the matter density that makes up NSs reach densities up to a few times the nuclear saturation density, however, until these days, detailed information about the characteristics and nature of their deep interiors is still lacking.Future multi-messenger signatures hold the promise of identifying the specific internal aspect of NSs.Theoretically, asteroseismology is widely employed to analyze the internal structure of compact stars -the name used for white dwarfs, neutron stars, hybrid stars, or strange quark starsto investigate the thermodynamic properties inside these objects.Through this diagnostic technique, analyzing the frequency modes can obtain a solid way to learn more about the physics inside compact stars.For example, if inside these stars a single component fluid is present [19][20][21][22][23][24][25][26][27][28][29][30][31][32] or the existence of a phase transition between layers with different mechanical properties [33][34][35][36][37][38][39][40][41][42][43][44][45][46].
As regards the radial perturbations of compact stars with a sharp interface, the set of equations must be solved by taking into account the additional boundary conditions at the phase-splitting interface.Around this point, there are two types of physical behavior due to radial perturbations: the slow and rapid phase conversion [35,36].In the case of slow conversion, there is no change of matter over the pulsating interface.On the contrary, the rapid conversion case involves a flow of mass from one phase to the other, and vice-versa, through the moving phase boundary.In recent years, the im-pact of the phase transition on the radial oscillations of compact stars has been reported in different articles.In the case of slow transition, for example, the authors concentrate on investigating the effects of the core formation [34], a sharp phase transition [37,38], the mixed-phase [39][40][41][42][43], and electric charge [44].On the other hand, among those reported considering both phase transitions in compact stars, we find: the ones that analyze the fast and slow conversion in the context of general relativity [36], how these are affected in the face of a magnetic field [45], and their influences on nonradial oscillations [35].In the rapid phase transition case, unlike the slow phase transition, in a sequence of equilibrium configuration, the maximum mass peak marks the beginning of radial instability.Indeed, in slow transitions, after this turning point, it is possible to find additional stable equilibrium configurations.Therefore, in a sequence of equilibrium configurations with increasing central energy density, some stars with the same mass but different radii are obtained.These stars are known as twin stars.
In the aforementioned articles, different models of equations of state are studied from the perspective of observational deformability data from the event GW170917.Some of these works study this phenomenon against the possibility of the existence of phase transitions inside the compact star, some taking the aspect of an analysis of the stability of these stars, and calculating radial oscillations.However, these same articles often make these calculations assuming only slow transitions, which allow the appearance of stable regions, in the mass-radius diagram, after the maximum mass.In this work, we present a detailed study of the influence of the phase transition and a stiffer fluid in NSs core on the equilibrium configuration, radial stability, and tidal deformability.In this sense, we analyze how the radius, mass, speed of sound, core radius, radial frequency of oscillation, and tidal deformation change when a phase transition and stiffer fluid in NSs core are considered.In the analysis of the radial stability of NSs, we will focus on the slow and rapid phase conversions.We also contrast our results with observational data to see the role that phase transition and a stiffer core fluid could play in the study of NSs.
The present article is arranged as follows: Section 2 presents the equilibrium equations and radial stability equations; moreover, this section is also devoted to presenting the junction conditions at the interface of the two-phase, which are required to investigate the slow and rapid phase conversions.Section 3 presents the EOSs employed for NSs, as well as the numerical method used to solve the complete set of equations required to investigate the equilibrium and radial stability.In Section 4 we show the numerical results for equilibrium configurations, radial stability, and tidal deformability of NSs with two-phase.Finally, we conclude in Section 5. Throughout the paper, we work with geometric units, i.e., c = 1 = G, and the metric signature +2.

Equilibrium equations
We take into account that the unperturbed neutron star is made up of layers of effective perfect fluids, whose energymomentum tensors can be expressed as with ρ, p, and u µ representing respectively the energy density, the fluid pressure, and the four-velocity.
To analyze the effect of phase transition in the dense matter on the equilibrium and radial stability of neutron stars, we set the space-time metric, in Schwarzschild coordinates, as ( The potential metric functions ν = ν(r) and β = β (r) depend on the radial coordinate r only.
For the energy-momentum tensor (Eq.( 1)) and line element (Eq.( 2)) adopted, with the potential metric e −β = (1 − 2m/r), we derive the set of stellar structure equations where the parameter m represents the mass inside the sphere of radius r.Eq. ( 4) is known as the hydrostatic equilibrium equation for a spherically symmetric static astrophysical object, also called as Tolman-Oppenheimer-Volkoff equation [47,48].
The stellar structure equations (3)-( 5) are integrated from the center toward the star's surface.At the center (r = 0) the integration starts The surface of the star (r = R) is determined by At this point, the interior solution connects smoothly with the exterior Schwarzschild vacuum solution.This indicates that at the star's surface the interior and exterior potential metrics are related through of the form: with M being the total mass of the star.

Radial oscillations equations
The radial pulsation equation is obtained by Chandrasekhar [49] making perturbation in the fluid and space-time variables.The perturbed quantities are placed into Einstein's field equation and in the linearized form of the conservation of stress-energy tensor.
The solution of the radial oscillation equation provides information about the eigenfrequency of oscillations ω.Intending to set this equation in a more appropriate form for numerical integration, we place it into two first-order equations for the variables ∆ r/r and ∆ p [24]; with ∆ r and ∆ p representing respectively the relative radial displacement and Lagrangian perturbations of pressure.Thus, the system of equations is, ξ = ∆ r/r: where dρ .The variables ξ and ∆ p have a time dependence of the form e iωt , with ω being the eigenfrequency.
To solve the differential equations ( 9) and ( 10), boundary conditions in the center and on the star's surface are required.Moreover, to find regular solutions in the center of the star, the second term of the right-hand side of Eq. ( 9) must vanish in r → 0. In this way, it is considered At this point, for normalized eigenfunctions, we regard ξ (r = 0) = 1.On the other hand, as established above, the surface of the star is determined when p(R) = 0.It implies 2.3 Tidal deformability Tidal effects are very common in the context of NSs binary systems.In fact, the gravitational field generated by one star in a binary system can result in deformation in its companion.The parameter of tidal deformability is the measure of the deformation in compact stars due to an external field.
From a mathematical point of view, this parameter can be expressed in terms of the fraction, where Q i j is the quadrupole moment perturbed by an external tidal field ε i j [16][17][18].The tidal deformability parameter λ 1 is connected with the Love number k 2 through the relation k 2 = 3 2 λ 1 R −5 .Moreover, the dimensionless tidal deformability Λ can be written in terms of the Love number k 2 as with C = M/R being the compactness parameter.k 2 can be expressed in terms of parameter y R ≡ y(r = R) as follows The parameter y(r) is calculated along the whole of the star -from the center to the surface of the star-integrating the equation together with the set of equations ( 3)-( 5), considering at the center of the star y(r = 0) = 2.The functions F = F(r) and Q = Q(r) are represented by the relations:

Junction conditions at the interface
In the last years, compact stars with two different phases have been considered a real possibility.However, there are still a bunch of open questions such as, for instance, the density at which a hadron-quark phase transition occurs and some discussions about the kind of phase transition depending on the surface tension between the phases [35,36,45,50,51].In our case, a first-order phase transition is considered, which results in the presence of a finite energy density discontinuity.Moreover, some approaches are regarded to investigate the radial stability and deformability of stars.

Radial oscillations
The phase transitions can be classified as slow or rapid depending of the time scale of the reaction of the matter in the neighborhood of the hadron-quark interface [36,52].A scenario of slow phase transition appears when the rate of reaction transforming one phase into another is much greater than those of the radial perturbations.In such circumstances, there is no flow of matter across the surface splitting the two phases.Such a condition implies that ξ must always be continuous at the interface across the interface, i.e., Moreover, this also leads to a continuity pressure at the interface of the two phases.At this point, it assures that In the case of a rapid phase transition, the rate of reaction transforming one phase into another is much lower than those of the radial perturbations.In this scenario, a change of mass flow through the interface occurs.The diminution of mass on one side should be equal to the increase of mass on the other side.This condition, together with the demand for the continuity of pressure, lead to where the prime stands the operation with respect to the radial derivative, and

Tidal deformability
At the interface of the two layers, we can clearly note that exists a singularity in Eq. ( 18) due to the speed of sound (d p/dρ) since, at this point, we have the same value of the fluid pressure for two different energy densities.In Ref. [17], the authors discussed this problem in the context of the surface vacuum discontinuity for incompressible stars.After in Refs [53][54][55], this approach was extended for the case of first-order transitions inside hybrid-stars, where authors concluded that at this point the function y(r) must follow the next condition: with r tr and ε being the radial position where phase-transition occurs inside of the star and ε an infinitesimal parameter, respectively.Obviously, the region r < r tr represents the core of the star, and the region r > r tr depicts the envelope of the star.
3 Equation of state and numerical method

Equation of state
To describe the matter that makes up the compact object, in the two-phase configurations, the relativistic polytropic equation of state [56] is adopted.Then, the energy density and fluid pressure of each phase are respectively connected through the relations: These two relations bear parameters from the inner and outer regions denoted by the sub-indexes "inn" and "out", respectively; namely, K inn and K out are the polytropic constants, Γ inn and Γ out being the polytropic exponents, and p dis inn and p dis out represent the phase transition pressure.Following [27], we set the inner polytropic constant value: Since the fluid pressure should be continuous along the star, at the phase transition point, where the inner and outer phase transition pressures meet the condition p dis out = p dis inn , the outer polytropic constant takes the form: At this point, the inner and outer phase transition energy densities are related by: where λ is known as the density jump parameter.Some examples of the EOS with a sharp density jump are presented in Fig. 1.
In the next section, we consider the value of Γ inn = 2.4, Γ inn = 2.6 and Γ out = 2.4 and the phase transition parameter λ between the values 0.5 and 1.0, in order to analyze both the role of a stiffer core fluid and the effects of a more abrupt phase transition in the star equilibrium configuration, respectively.The values of Γ inn ,Γ out and λ chosen allow us to obtain comparable results with the data found through observation, which are reported by the LIGO-Virgo network in [57] and by NICER in [58][59][60][61][62][63][64].

Numerical method
The effects of the phase transition on the equilibrium configurations and tidal deformations are investigated through the numerical solution of the system of equations, boundary conditions, and junction conditions established in the section 2, for each ρ dis inn , λ , Γ inn , and Γ out .This system of equations is integrated from the center toward the star's surface.
The analysis of the radial stability starts by solving the stellar structure equations by using the Runge-Kutta fourthorder method in order to determine the radial pulsation coefficients.After, we begin at the star's core with the solution of Eqs. ( 9)- (10) for a test value of ω 2 .These equations are numerically integrated outwards until the interface is found, where the junction conditions are employed (for the slow case Eqs. ( 19) and ( 20) and for the rapid case Eqs. ( 21) and ( 22)) to find the values of ξ and ∆ p at the other side of the interface.Then, the numerical integration continues towards the surface of the stars attempting to reach the conditions (7) and (12).Whether, after each integration, the equality (12) is not fulfilled, ω 2 is corrected until satisfying this equality in the next integration.The parameters ω 2 which satisfy the oscillation are called eigenvalues of the radial pulsation equation and ω of eigenfrequencies (review [36]).
We are interested in analyzing the radial stability of stars, we only analyze the lowest eigenvalue, i.e., ω 2 0 .When ω 2 0 > 0, the star is stable against small radial perturbations.ω 0 is known as the eigenfrequency of the fundamental mode.In the panel, it is noted the diminution of the total mass with the density jump parameter.This is associated with the fact that the fluid pressure decays abruptly due to the presence of a phase transition, being this declines greater for lower λ (see Fig. 1).The change of the mass with Γ inn is also observed in Fig. 2. For a fixed ρ c and λ , for a greater Γ inn , a larger total mass is derived.This point can be understood since for larger interior polytropic exponents larger central pressures are obtained.Thus, larger central pressure supports more mass against gravitational collapse.In addition, it is important to say that, for larger Γ inn , neutron stars with more compact cores are obtained (review, e.g., [65,66]).

Equilibrium configuration of neutron stars
Fig. 3 shows the mass as a function of the total radius for two values of Γ inn , for ρ dis inn = 8 × 10 17 [kg/m 3 ], few values of λ , and Γ out = 2.4.As in the case of Fig. 2, in this figure is considered ρ c between 10 18 and 10 20 [kg/m 3 ].The panel shows an increment of M/M with the diminution of R until to find a M max /M .After this point, the curves turn anticlockwise, to M(R) starts to decrease with R until to reach  R min .From here on, the mass decays with the increment of the total radius.In Fig. 3, for some range of central energy density, we also find an increment of the total radius with the diminution of the density jump parameter.This is due to the fact that in these compact stars, the pressure of the fluid decays slower with the increase of the radial coordinate, thus obtaining larger radii.In addition, for some range of ρ c and λ , we can also observe decrements of the total radius with the increment of Γ inn .Despite having an increase of the central pressure with Γ inn , it decays faster with the growth of the radial coordinate.In this way, compact objects have a smaller total radius.Finally, in Fig. 3, we also see that the radius of the stars with maximum mass decreases with the jump in phase transition density.This is due to the fact that the central pressure of the star decreases with λ , thus, the pressure decays faster with the radial coordinate.
The total mass versus the speed of sound is plotted in Fig. 4, for ρ dis inn = 8 × 10 17 [kg/m 3 ], some values of λ , and Γ out = 2.4.In this figure, we only present equilibrium configurations with the speed of sound lower than the speed of light c 2 s = 1.0.As can be seen, in particular, at stars that present first-order phase transition with low central energy densities (stars with low total masses), the speed of sound never exceeds the conformal limit c 2 s = 1/3.However, in stars with larger total masses, the speed of sound exceeds the conformal limit value but is far to attain the speed of light.These results are in concordance with those reported in the article [32].

Radial stability of neutron stars
Fig. 5 shows the behavior of the slow and rapid eigenfrequency of the fundamental oscillations as a function of the total mass and against the central energy density for some different values of λ , ρ dis inn = 8×10 17 [kg/m 3 ], Γ out = 2.4, and Γ inn = 2.4 on the left panels, and Γ inn = 2.6 on the right panels.From the figure, for λ = 1, it can be noted that the maximum total mass is found at the zero eigenfrequencies of oscillation.This case represents the usual study of radial oscillation of neutron stars in the absence of phase transition.In turn, for λ < 1, as well as in [36], we note that the total mass at the null eigenfrequency of oscillation depends on the type of the phase transition.At this point, when Γ inn = Γ out = 2.4, the difference of the mass attained in the rapid and slow case is almost 0.246% for λ = 0.9, 0.470% for λ = 0.8, 0.397% for λ = 0.7, 0.299% for λ = 0.6, and 0.164% for λ = 0.5.In the case Γ inn = 2.6 and Γ out = 2.4, the difference of the mass is around 0.197% for λ = 0.9, 0.358% for λ = 0.8, 0.330% for λ = 0.7, 0.161% for λ = 0.6, and 0.0371% for λ = 0.5 (review Table 1).
In Fig. 5, in the slow case, the total mass at the zero eigenfrequencies of oscillation is derived at ρ c larger than the one employed to obtain the maximum mass value; i.e., twins stars are derived, stable stars with the same total mass but with both different central energy densities and total radii.However, in the rapid case, the maximum mass and the null eigenfrequency of oscillation are obtained by using the same value of ρ c .It indicates that in a sequence of static equilibrium configurations, the maximum mass point marks the beginning of the instability against small radial perturbations.This characteristic in each phase transition could be useful to differentiate them.Table 1: The central energy density ρ c , the total mass M/M , the total radius R, and the core radius R core where the null eigenfrequency of oscillation of the slow case ω 0,s and rapid case ω 0,r are derived.These parameter are found for ρ dis inn = 8 × 10 17 [kg/m 3 ], Γ out = 2.4 and different values of Γ inn and λ .
Table 1 presents the central energy densities and total masses where the zero eigenfrequencies of oscillations for the slow and the rapid phase transition.These parameters are derived for ρ dis inn = 8 × 10 17 [kg/m 3 ], Γ out = 2.4 and some values of Γ inn and λ .In the table, for a fixed Γ inn , at the null eigenfrequency of oscillation for the slow and rapid case, we note that the total mass, the total radius, and the core radius decrease with the density jump parameter.This could be understood since p c decays with λ , in this way, the total pressure diminishes faster with the growth of the radial coordinate.On the other hand, when Γ inn is increased, i.e., when a stiffer core fluid is considered, stars with a core radius R core closer to the total radius R are found.

Tidal deformability in the light of GW170817
Tidal deformability against the total mass of stable NS is plotted at the top panel of Fig. 6 for ρ dis inn = 8 × 10 17 [kg/m 3 ], Γ out = 2.4, different values of λ and Γ inn .On the panel, it is also presented the tidal deformability constrained by the event GW170817 for a star of 1.4M to be 70 ≤ Λ 1.4M ≤ 580 [57], at 90% confidence level, for low-spin priors.
For a fixed λ and masses range, at the top of Fig. 6, we note an increment of the tidal deformability when a stiffer fluid is considered in the core.On the other hand, by setting the parameter Γ inn for different values of the phase transition parameter λ , we can notice that the importance of the parameter Λ changes in a more relevant way.The lower the value of λ , the less the pressure of transition (see Fig. 1) and the smaller the value of the deformability for the same mass.From these results, the effect of phase transition and stiffer fluid in the core are noticeable in the tidal deformability.Nonetheless, between these two factors, the phase transition parameter affects the NS properties more.
At the bottom panel of Fig. 6, the curves Λ 1 − Λ 2 are presented for a binary NS system with chirp mass equal to GW170817.Since the total mass is associated with the dimensionless tidal deformability (top panel of Fig. 6), the curves Λ 1 − Λ 2 are obtained once chosen a value of m 1 and calculating m 2 for the fixed value of the chirp mass M = 1.188M [9], which is theoretically calculated by the relation: The values considered for m 1 and m 2 run from 1.36M ≤ m 1 ≤ 1.60M and 1.17M ≤ m 2 ≤ 1.36M , respectively.At the bottom panel of Fig. 6, we investigate the effects of the phase transition and a stiffer fluid in the core (for ρ dis inn = 8×10 17 [kg/m 3 ]) of two neutron stars in a binary system.In this case, we note that both the phase transition and a stiffer equation of state could play an important role in the detection of these compact objects.From the results, we note that there is an interval for the density jump parameter λ , 0.7 ≤ λ ≤ 0.8, which plays the compact stars inside of 50% and 90% regions.The compact stars with λ ≥ 0.9 are Fig.6: Top: Dimensionless tidal deformability Λ as a function of the total mass in Solar masses.The vertical dash-dotdot line marks the tidal deformability of event GW170817 estimated in [57].Bottom: The dimensionless tidal deformability Λ 1 and Λ 2 for a binary NS system with masses m 1 and m 2 and the same chirp mass as the event GW170817 [9].We only take into account the combination with m 1 > m 2 .The diagonal short dot line denotes the Λ 1 = Λ 2 limit.The solid top yellow line indicates the 90% credibility level and the bottom yellow solid line is the 50% level established by LIGO-Virgo scientific network in the low-spin prior scenario.Only stable equilibrium configurations with slow conversions at the interface are shown.outside the 90% region and do not appear in the panel.This is realized with the aim that the curves shown inside 50% and 90% can be seen clearly.Furthermore, we observe that for smaller values of λ the curve change to smaller values of dimensionless deformability.On the other hand, it can be noted that a stiffer fluid in the core produces larger values of deformability.
4.4 Change of stellar physical parameters with phase transition energy density and its comparison with observational data In Fig. 7, the mass-radius curves are compared with the observational data considering Γ inn = Γ out = 2.4, some values of λ and two ρ dis inn .On the top and bottom panel are employed ρ dis inn = 7.0×10 17 [kg/m 3 ] and ρ dis inn = 9.0×10 17 [kg/m 3 ], respectively.In this figure is used 10 18 ≤ ρ c ≤ 10 20 [kg/m 3 ].The observation data correspond to the NICER constraints obtained from the pulsars PSR J0030 + 0451 [58,59] and PSR J0740 + 6620 [60,61].The corresponding bands of the pulsars PSR J0740 + 6620 [62], PSR J0348 + 0432 [63] and PSR J1614 + 2230 [64] are also presented.From the figure, we observe that the change of ρ dis inn affects the stellar struc-ture configuration; being this change more noticeable in the range of low central energy densities.In this interval, for larger ρ dis inn , greater total mass and smaller total radius are found.From these results, we note that the change of internal phase transition energy density allows us to obtain some results more accurate and closer to empirical evidence of the neutron stars PSR J0030 + 0451.In addition, we see that with the increment of ρ dis inn , grow the possibility of having equilibrium solutions with sharper phase transition (smaller density jump parameter λ ) within PSR J0030+0451.On the other hand, from Figs. 3 and 7, we observe that in the range of larger total mass, a stiffer fluid in the core could help to reach empirical evidence of neutron stars PSR J0740+6620, PSR J0348 + 0432, and PSR J1614 + 2230.
In Fig. 8, the top and bottom panels, respectively, present the tidal deformability against the total mass and Λ 1 -Λ 2 curves for a binary NS system with chirp mass equal to GW170817 considering the relation (29) where m 1 and m 2 goes from 1.36M ≤ m 1 ≤ 1.60M and 1.17M ≤ m 2 ≤ 1.36M .The inner phase transition energy density considered on the left and right panels are 7.0×10 17 [kg/m 3 ] and 9.0×10 17 [kg/m 3 ], respectively.In all panels of Fig. 8, only stable equilibrium configurations for the rapid transition case, employing Γ inn = Γ out = 2.4, are considered.As stated above, the observation data correspond to the event GW170817.
In Fig. 8, on the top panels, it can be seen that all curves tidal deformability (Λ )-total mass (M/M ) decay with the increment of ρ dis inn .From this, we understand that this phenomenon allows us to have equilibrium configurations with a lower density jump parameter λ (a sharper phase transition) within the observational data of the GW170817 event.
From the bottom panels, we note that some equilibrium solutions located outside the range of observational data derived from the GW170817 event fall within this interval when we increase the phase transition energy density ρ dis inn .In addition, from Figs. 8 and 6, we can also say that having a stiffer fluid in neutron stars' cores helps us to have static equilibrium configurations within the range of observational data.

Conclusions
In this work, we investigated the influence of the phase transition on the equilibrium, radial stability, and tidal deformability of NSs with a stiffer fluid in the core.In the core and the envelope of the star, the relativistic polytropic equation of state is considered.The spherical equilibrium configurations are connected smoothly with the Schwarzschild exterior spacetime.We examined the change of the mass, radius, speed of sound, core radius, the eigenfrequency of the fundamental mode of the star with a slow and rapid phase conversion at the interface, and tidal deformability for different density jump parameters λ , phase transition energy Fig. 8: Top: Dimensionless tidal deformability Λ against the total mass in Sun masses.The vertical line of the dash-dotdot plots the tidal deformability of the GW170817 event measured at [57].Bottom: The dimensionless tidal deformability Λ 1 and Λ 2 for a binary NS system with masses m 1 and m 2 and a chirp mass of 1.4M considering the combination with m 1 > m 2 .The diagonal short dot line marks Λ 1 = Λ 2 limit.The solid top and bottom yellow lines denote respectively 90% and 50% level established by LIGO-Virgo scientific network in the low-spin prior scenario.On the left and right panels are used ρ dis inn = 7.0 × 10 17 [kg/m 3 ] and ρ dis inn = 9.0 × 10 17 [kg/m 3 ], respectively.Only stable equilibrium configurations with rapid conversions at the interface considering Γ inn = Γ out = 2.4 are presented.densities ρ dis inn , interior polytropic exponents Γ inn , and exterior polytropic exponent Γ out = 2.4.
As well as in the study of NSs developed in [67], which employs a non-relativistic polytropic equation of state, we note that some aspects of the static equilibrium configurations -such as the mass and radius-are affected by the phase transition, stiffer fluid in the core (which change with Γ inn ), and phase transition energy density.
For the values Γ inn and λ employed, in the slow case, the zero eigenfrequencies of the fundamental mode are attained beyond the maximum mass points and, in the rapid case, the maximum masses points mark the beginning of the radial instability thus indicating that the regions constituted by stable and unstable stars can be recognized by the conditions dM/dρ c > 0 and dM/dρ c < 0, respectively.These results are in concordance with those one reported in the works [35,36,45].
The change of the tidal deformability for a NS (Λ ) and a binary NS system (Λ 1 and Λ 2 ), with equal chirp mass as GW170817 event, as a function of Γ inn and λ , has been analyzed.We obtained a dependence of the dimensionless tidal deformability with these two factors in the aforementioned frameworks.For NS configurations, for some interval of masses, we noted that Λ grows and decreases with the increment of Γ inn and diminution of λ .In turn, for a binary NS scenario, we showed that the phase transition and stiffer fluid in the NS core could also play an important role in the detection of NSs.These results are in agreement with those ones published in [38].
We also investigated the dependence of some physical parameters of NSs with the phase transition energy density.
At the interval of low central energy densities, we found that ρ dis inn can also be important in the study of NSs since their physical parameters could be significantly affected by the value of phase transition energy density.
Finally, we noted that a change in jump density parameter λ , phase transition energy density ρ dis inn , and stiffer core fluid Γ inn could lead to the possibility that some equations of state that are outside of the observational data, can be inside this framework for accurate values of λ , ρ dis inn and Γ inn .

Fig. 5 :
Fig. 5: The slow (s) and rapid (r) eigenfrequency of the fundamental mode ω 0 as a function of the total mass M/M and against the central energy density for ρ dis inn = 8 × 10 17 [kg/m 3 ], Γ out = 2.4, four values of λ , and two values of Γ inn .On the left panels, it is used Γ inn = 2.4 and on the right panels, it is employed Γ inn = 2.6.
O a n d V i r g o ( 5 0 % ) L I G O a n d V i r g o ( 9 0 % )