Role of Anisotropy on the Tidal Deformability of Compact Stellar Objects †

: In this paper, we introduce a framework to study the tidal deformation of relativistic anisotropic compact stars. Anisotropic stresses are ubiquitous in nature and widely used in modelling compact stellar objects. Tidal deformability of astrophysical compact objects is a natural effect of gravity, such as one produced by a companion in a binary system. In general relativity, the existence of this measurable effect of gravity can be quantiﬁed by their tidal Love numbers (TLN), which characterize the deformability of a neutron star (NS) from sphericity. The tidal deformability or polarizability parameter of an NS depends on its complex internal structure, and hence, the nature of the compact object can be studied by measuring the TLN. We choose a particular solution, which is the anisotropic generalization of the Tolman IV model, as the interior of the compact stellar object. The physical acceptability of the model has been shown graphically by considering the pulsar 4 U 1608-52 with their current estimated mass and radius. By computing the quadrupole moment, we found that the TLN is dependent on anisotropy of the compact object. We graphically analyze the variation of the TLN against anisotropy for different compact objects with a compactness factor. The numerical value of TLN is given for different compact objects for physically acceptable values of the anisotropic parameter.


Introduction
Compact objects are extremely dense astro-physical objects that provide strong gravity and high density, allowing for the study of fundamental physics related to nuclear matter properties. In general, compact objects exist with their binary companion is a natural setup. In this binary setup, a compact star is assumed to be immersed in the tidal field of its companion and hence produce tidal deformation. The possibility of the generation of the gravitational waves during extremely fast rotational motion of a binary system was predicted. Recently, a gravitational wave has been detected by advanced astronomical observations of LIGO and Virgo collaborations from the binary neutron star merger event GW170817 [1]. A neutron star placed in a perturbing external gravitational field is deformed and induces a quadrupole moment, which affects the binding energy of the system and increases the rate of emission of gravitational waves during the late stage inspirals. The recent observational data based on the measurement of the tidal deformability impose a stringent constraint on the allowed equation of state (EOS). In particular, the theoretical prediction of the mass and radius of a neutron star (NS) mostly depends on the nature of the nuclear EOS at supra-nuclear densities. In this context, tidal deformability can be used to study their interiors. The EOS of neutron stars involves their microscopic properties and uniquely determines the macroscopic properties, such as the maximum allowed NS mass, radius and tidal effects. The tidal response is the astrophysical constraints that can be employed as probes of NS properties. It is the astrophysically observable macroscopic property of NS that can be defined as the ratio of the induced multipole moment of a star over the induced tidal field from its companion. The tidal Love number (TLN), which is the ratio of the induced quadrupole moment to the perturbing tidal gravitational field, can be expressed by a relatively simple analytical formula.
On the other hand, the existence of pressure anisotropy, the difference of radial and transverse pressures, is ubiquitous in a compact star. The source of pressure anisotropy in a compact star can be due to various reasons, e.g., pion and kaon condensates [2,3], high density, existence of a solid core or type 3A superfluid [4,5], strong magnetic fields [6], a mixture of a perfect and a null fluid, viscosity, and phase transition [7], etc. There are several works available in the literature where incorporating anisotropy into the matter distribution of compact objects in the background of general relativity (GR) has been addressed for various issues of the compact structures [8][9][10][11][12][13][14].
In this paper, we assume a known solution, which is an anisotropic generalization of the Tolman IV model, to describe the compact star. We have calculated the TLN, which actually measures the tidal deformability of the compact object induced by the external field.

Physical Features and Tidal Love Number
The tidal distortion of NSs in a binary system connects the EOS, describing the nature of the matter composition star with that of the gravitational wave emission during the inspiral [15]. We consider a static spherically symmetric star, immersed in an external quadrupolar tidal field E ij [16,17] arising due to its binary companion. The star in response to the tidal field develops quadrupole moments Q ij , which can be related to the linear order external tidal field E ij as [17].
where Λ is the tidal deformability of the NS and is related to the l = 2 dimensionless TLN k 2 as [17] The background geometry of spacetime of a spherical static star can be written as For the spherically static metric (3), the stress-energy tensor is given as where u ξ u ξ = −1, η ξ η ξ = 1 and η ξ u ξ = 0.
We choose a particular model, which is an anisotropic generalization of the Tolman IV Model [18], given as It is interesting to note that for α = 0, this solution reduces to the well-known Tolman IV solution [19].
For the line element (3), the independent set of the Einstein field equations are then obtained as 8π where primes ( ) denote differentiation with respect to r. In the field Equations (7)- (9), we have assumed G = 1 = c. The system of equations determines the behaviour of the gravitational field of an anisotropic imperfect fluid sphere.
For the assumed solutions we have where we define ∆ = 8π(p t − p r ) as the measure of anisotropy of the spherical system. The exterior Schwarzschild metric is across the boundary of the star r = R, where M is the total mass of the sphere. Making use of the junction conditions, the constants A, B, C are determined as Now the background metric (0) g µν (x ν ) under the influence of an external tidal field with a small perturbation h µν (x ν ) gets modified as For the linearized metric perturbation h µν , using the method as in [20,21], without loss of generality, we restrict ourselves to static l = 2, m = 0 even parity perturbation. The perturbed metric, with the assumption that the tidal deformation will be axis symmetric around the line connecting the two stars, which we take as the axis of spherical harmonic decomposition, becomes h µν = diag H 0 (r)e 2ν , H 2 (r)e 2λ , r 2 K(r), r 2 sin 2 θK(r) Y 2m (θ, φ).
The vacuum exterior of the star is of Schwarzschild type, so that by setting ρ = 0, p r = 0, p t = 0 and e 2λ = 1/(1 − 2M/r), the master Equation (27) becomes The solution to this second order differential Equation (30) is where c 1 and c 2 are integration constants. In order to obtain the expressions for these constants, we perform a series expansion of Equation (31) Now, at large r, the metric coefficient g tt is given by [17]: where n i = x i /r. Matching the asymptotic solution from Equation (32) with the expansion from Equation (33) and using the Equation (1), we have Using Equations (34), (31) and (2), we obtain the expression for TLN k 2 as follows: where Here, the compactness C = M R and y depend on r, H and its derivatives evaluated at R with To calculate numerically the value of k 2 for a particular NS [22], one needs to modify the master Equation (27) using Equation (37) as ry + y 2 + (rR − 1)y + r 2 S = 0. (38)

Results
The recent data available from the pulsar 4U 1608 − 52 (for the star to be composed of an anisotropic fluid distribution with α = 0.5), whose estimated mass and radius are M = 1.57 M and R = 9.8 km, respectively [23,24], are used to find the constants are calculated as A = 0.53953, B = 0.291097, C = 0.008452. We set a = 1. Making use of these values, we show the nature of all the physically meaningful quantities graphically in Figure 1. Using the initial condition y(0) = 2 and all the mentioned equations, Equation (38) for a particular NS, can be numerically calculated. Having the numerical value of y, for a particular NS, from Equation (35), the TLN k 2 can be obtained numerically.

Conclusions
The plots clearly show that all the quantities comply with the requirements of a physically viable realistic star. In particular, the figures highlight the effect of anisotropy on the gross physical behaviour of the compact star. In Figure 2, the TLN k 2 is plotted against α for different compact objects with compactness C. From this panel of figures, we note that k 2 decreases monotonically with increasing α. This is a much expected physical property of a compact object with anisotropy.