A generalized model for compact stars

By virtue of the maximum entropy principle, we get an Euler–Lagrange equation which is a highly nonlinear differential equation containing the mass function and its derivatives. Solving the equation by a homotopy perturbation method we derive a generalized expression for the mass which is a polynomial function of the radial distance. Using the mass function we find a partially stable configuration and its characteristics. We show that different physical features of the known compact stars, viz. HerX-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Her~X-1$$\end{document}, RXJ1856-37\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$RXJ~1856-37$$\end{document}, SAX J (SS1), SAX J (SS2), and PSRJ1614-2230\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$PSR~J~1614-2230$$\end{document}, can be explained by the present model.


Introduction
To find an analytical solution of nonlinear equations has always been a great challenge in relativistic astrophysics as well as cosmology. Considering different methods and approximations mathematicians and relativists attempted to find the suitable interior solutions of physical interest. It is in this aspect seen that the homotopy perturbation method (HPM) is a powerful yet simple tool to solve several highly nonlinear equations with the least number of assumptions as can be seen in the literature [1][2][3][4][5][6][7][8][9][10][11][12][13].
It is interesting to note that in one of his earlier works He [4] proposed a coupling method of a homotopy technique and a perturbation technique to solve a nonlinear equation. In contrast to the traditional perturbation methods, the proposed method by He did not require any small perturbation in the equation to obtain an effective and simple solution. Thereafter several researchers have employed the HPM in the diverse fields of mathematics and physics [14][15][16][17][18].
For special mention in the field of astrophysics, it is observed that by using the MIT bag model, Rahaman et a e-mail: azizmail2012@gmail.com b e-mail: saibal@associates.iucaa.in c e-mail: rahaman@associates.iucaa.in al. [19] obtained a deterministic model for a strange star where they considered a mass polynomial and analyzed critically all the physical properties from the model. However, their model was unable to explain the physical properties up to 6 km from the center of the spherical system.
Deb et al. [20] attempted to find a new class of solutions relating to a radiation model by using the HPM. The spherical distribution was considered to be made of a perfect fluid where radiation is along the radially outward direction of the system. Here we have employed the HPM as a new tool for astrophysical systems which via the mass polynomial facilitates solving of the field equations. A set of interior solutions has been found on the basis of the equation of state in the form p = 1 3 ρ, where p and ρ are the radial pressure and density of the matter distribution, with the requirement that the radiation model is consistent with the Schwarzschild metric. It is observed that the set of solutions thus developed provides a new metric to explain strange stars, which, however, seems to be not free from a singularity.
In another work Rahaman et al. [21] employed the HPM for a spherically symmetric system of a radiating star which was shown to suffer from an instability problem. It is observed from the mass function that it was possible to explain the features of an E0 type brown dwarf star with 0.065 M (i.e. 45227 km) whereas the stability of the brown dwarf could not be explained which suffers from an instability within the region of the core radius 41800 km. However, in the present investigation we propose a general mass function and use it for different cases.
Therefore the motivation of the present work is to provide a general model for compact stars by employing a polynomial parameter n and to find the different features of compact stars with this mass function. We derive expressions for the density, pressure, redshift etc. and have shown that several known compact stars, e.g. Her X − 1, R X J 1856 − 37, S AX J (SS1), S AX J (SS2), and P S R J 1614 − 2230, can be well explained by this model. The plan of the investigation is as follows: in Sect. 2 we have calculated the mass of the system by using the maximum entropy principle (MEP), and thereafter by using the HPM we have obtained a suitable expression for the mass. We have reformatted Einstein's field equations by using the mass function and obtain solutions for the system in Sect. 3. Also in Sect. 3 we have explored different physical features, viz. the density and mass, pressure and anisotropy, stability, energy conditions, TOV equations, compactness, redshift, etc. with elaborate discussion and a comparative study of the available data from observation with that of the predicted values. In Sect. 4 we have made some concluding remarks regarding different aspects of the present model.

Maximum entropy principle
Let us consider the metric of a system which has spherical symmetry as follows where m(r ) is the mass distribution of the spherically symmetric body and g tt (r ) is the time-time component of the metric. These are functions of the radial coordinate r only. Also we consider a linear equation of state for the system with thermal radiation given by where p r and ρ are the radial pressure and density of the matter distribution. Compatible with the spherically symmetry, we assume the general energy momentum tensor to be with u μ u μ = −η μ η μ = 1.
The energy density and entropy density take the forms for locally measured temperature T [22] as follows: where b is a constant of order unity (in Planck units G = c = h = k = 1), on the assumption that the number of species of radiation is of order unity [23]. For the matter distribution up to radius r ≤ R, the total entropy is given by which is simplified as where the Lagrangian L is given by Here, s = α(ρ) 3 4 , α = 4 3 b 1 4 , and the Einstein equations yield which can be reduced to We see that the equation contains the mass, its derivatives, and products of them. So it will be very difficult to get an exact solution for this highly nonlinear differential equation. Therefore we use a homotopy perturbation method to solve the equation and thus we get an approximate solution for the mass.

Homotopy perturbation method
In the HPM we consider the linear and nonlinear parts of the differential equation separately. Then we build the homotopy structure [4] perturbed with a parameter, say , as follows: For the present problem we take the linear part L(m) as dm dr and the rest of the terms as the nonlinear part. So the above equation becomes where m 0 is the initial solution of the linear part of the equation. Now we assume the general solution By varying the value of from 0 to 1 we shall get the desired solution of the nonlinear equation (10).
Since we are looking for a more general expression for the mass we assume the initial solution as m 0 = ar n where we choose 'a' as model parameter, r as radial distance; n is an index assumed to have any value.
Substituting this in Eq. (12) and equating coefficients of the different orders , we get the solution for the mass up to a second order correction as Again, adding the correction terms to the initial solution we get the general expression for the mass as where A and B are the constants of integration.
Here we get the mass as a polynomial function of the radial distance with different powers of n. Now this equation can be utilized to get a different expression of the mass for different values of n. It is the mass function here which actually governs different features of the stars. So, primarily the goal of the paper is to develop a generalized mass function using the HPM. At a first glance, we observe that the mass function becomes undefined for n = 1 and some fractional values of n, e.g. 1/2, 2/3. The values of the n should be such that the resultant mass function be a well-defined, positive, and increasing function of the radial distance. For this purpose we need to check the physical validity of the mass using the real values (integer or non-integer) of n. In the present model we therefore consider only those integer values of n which show physical validity of the mass and other basic features of the stars.
We begin our investigation by taking integer values of n as follows: (i) For n = 0, we find that mass is independent of the radial distance i.e. it is constant in nature. But variation of the mass of the star is mainly responsible to all the interstellar phenomena. (ii) For n = 1, the mass function is undefined.
(iii) For n = 2, mass is a well-defined function of the radial distance.
Hence the mass function (15), under the above restrictive values of n, takes the form Now, we use the boundary conditions with R being the radius of the sphere. Thus, we get and Substituting A and B in Eq. (16), we get To explore different features of a star we use this mass function in Sect. 3.

Einstein's field equations and the physical parameters of a compact star
The Einstein field equations for the matter distribution given in Eq. (3) can be written as 2m 2m From Eq. (22) we get the density Using Eq. (18) we get an expression for the density: Again, by using Eq. (2) we solve Eq. (23) for g tt as follows: where K is a constant.
It can be noted that the exponential term cannot be evaluated analytically. So we have to make approximations to evaluate it. Also to find the constant K we need boundary conditions. Actually we do not need to calculate the proper form of g tt in terms of the radial distance as we are only interested in finding different features such as density, pressure, etc. By intuition we see that these features are derivable from the Einstein field equations and the mass function. This is because the Einstein field equations contain ratios of g tt , g tt , g tt .
The total mass of the system is given by We consider 2m(r ) r < 1 and simplifying the above integration we get the expression for the total mass as

Mass and density
We begin our calculations for the group of stars of E, M, K , G type etc. Calculations for these stars show that the nature of the tangential pressure is negative in nature which is against the regularity condition for a perfect fluid. Also it is observed that these stars suffer from instability. While investigating a compact star we find that there are some specific values of R and 'a' for which our model could describe them. For example, we try to explain features of a compact star, namely Her X − 1, of radius 6.7 km and mass of 0.9824 solar mass. We see that for the model parameter R = 10 km our model could describe all the features of the compact star. For this purpose we find a suitable value of the parameter 'a', which amounts to −0.01860326472 km −1 . The variations of the mass and the density are shown in Fig. 1. We see that the mass and density are positive and they are decreasing functions.

Pressure and anisotropy
Using Eqs. (2) and (26) we get an expression for the radial pressure, From Eq. (24) we obtain the expression for the tangential pressure  Fig. 2 we see that the radial pressure of the system are positive and decreasing functions of the radial distance whereas the tangential pressure switches sign at r = 6.788 km. Now for a perfect fluid, positivity of the pressure and the decreasing nature of the pressure are the criteria to hold good. Therefore we require the maximum allowed radius for the star to be 6.788 km. Also we find that the total mass of star for this radius is 1.02 solar mass. These data almost match with the observed data for Her X −1. In the following sections we shall show the variations of different features of the system by graphs.
The anisotropy for the perfect fluid is = ( p t − p r ), plotted against r in Fig. 2. We see that is negative for any r . So, the anisotropic force is inward in this self-gravitating system.

Stability condition
We shall use Herrera's approach [24] [usually known as the concept of cracking (or overturning)] to check potential instability or stability of our configuration. According to Herrera's theorem, the region for which v 2 st − v 2 sr < 0 is a stable region. We observe from Fig. 3 that v 2 st − v 2 sr < 0 in the region for which r < 6.788 km. This implies that our configuration is potentially stable within this region. The radial and tangential sound speeds are defined as v 2 si = d p i dρ (Fig. 3).

Energy conditions
We shall now check the energy conditions for the present model for each case: We see from Fig. 4 that different energy conditions in our model are satisfied throughout the region up to approximately 8.5 km.

TOV equation
The relativistic Tolman-Oppenheimer-Volkoff (TOV) equation can be written in the following form: where M G = M G (r ) is the effective gravitational mass within the sphere up to radius r , given by The above equation indicates that the equilibrium situation is achieved for fluid elements subject to gravitational and hydrostatic forces, and forces due to anisotropy. Thus the above equation assumes the form where The profiles of F g , F h , and F a for the sources are shown in Fig. 4. It is to be noted that < 0 implies an attractive force, which means a hydrostatic force produces a repulsive force and essentially this stabilizes the system from collapsing due to the attractive forces F a and F g . These profiles of three forces indicate that our model provides a stable configuration.

(40)
The variations of the compactness and the redshift are shown in Fig. 5.

Minimum mass-radius relation
Buchdahl's condition is m(R) R < 4 9 < 0.44. Table 2 shows that m(D) D = 0.22 for the compact star Her X − 1. So Buchdahl's condition is satisfied in our model. Also the value of the surface redshift for this star is 0.19. We show the mass-radius ratios for some other stars also in Tables 1 and 2.
In Tables 1 and 2 the following symbols are used: R obs = observed radius, D pre = predicted radius, M R = observed mass, M D = predicted mass, M(D)/D = predicted massradius ratio, and Z (D) = predicted surface redshift. Here the data are considered from Ref. [25].

Conclusion
First in our concluding remarks and discussions we would like to make some notes on observations regarding the radial pressure, which is expected to be positive at the boundary r = 10 km. However, the tangential pressure in the present model actually becomes zero at r = 6.788 km. Also we observe that the stability condition by the Herrera approach [24] is well satisfied within the radius r = 6.788 km. Therefore, with reference to these observations, we can think of the formation of our stellar model of two phases. In phase I there is a shell of 2.3 km forming the outer region of the star whereas the phase II consists of a core of radius 6.7 km forming the inner region of the star. Now we observe from the related graphs that all the features of the star are physically valid within the core but some of the features do not show physical validity beyond the core i.e. within the outer shell. It is interesting to note that we find a well-known star, Her X − 1 of radius 6.7 km, which can be explained by our model very well if we take only the core part of the star into our physical consideration. For the present model of radiating compact stars as described above we observe the following interesting and salient features: (1) We obtain a general mass function by solving a highly nonlinear equation of the mass by the HPM. In the HPM the homotopy structure consists of a linear and a nonlinear functional part. The choices of these functional operators are such that the trivial problem (the linear part) is continuously deformed to the original problem with variation of the embedding parameter ∈ from 0 to 1. In the present study, however, we use a particular homotopy structure to solve the nonlinear differential equation for the mass function so that all results obtained in this paper are valid subject to an approximation given by this particular homotopy. (2) For n = 0 the mass has no dependence on the radial distance. Again for n = 1 the mass function becomes indeterminate. So we start our investigation with n = 2.
(3) We find that the tangential pressure becomes negative throughout the whole region for the different groups of stars (E, M, K , G type). However, in the literature we note that this type of special solutions could exist [26]. (4) The features of the compact star of radius 6.788 km are well explained. We figure out that the mass of the star is 1.02 solar mass. This matches the observed compact star Her X − 1 of radius 6.7 km and a mass of 0.9824 solar mass. Some other compact stars are also shown in Tables  1 and 2. (5) For n = 3, as a particular case of the present paper, we have already shown that the model supports a brown dwarf of E0 type rather than compact stars. Also we have shown there that the brown dwarf is partially stable [21].
Regarding the MEP we have made some specific observations: (1) Based on the work by Rahaman et al. [21] we note that there are several papers available in the literature where one may find the application of the MEP in the anisotropic case. In this connection a comment by Richstone and Tremaine [27] seems convincing: "It is sometimes argued that in maximum-entropy models the velocity-dispersion tensor must be isotropic and hence the distribution function must be a function of energy alone. In fact this is not correct when we seek the maximum-entropy model subject to the constraint of a measured density profile. In this case, might gradients in the dispersion cause anisotropy in the dispersion tensor of a maximum-entropy distribution function" Thus, as a future plan we would like to consider the anisotropic case under the MEP as has been performed by several researchers under different requirements [28][29][30]. However, applications of the HPM can be found in the field of astrophysics in different contexts paving the way to a new research arena [31,32]. (2) Rahaman et al. [21] argued that instability may be an inherent property of any radiating compact stars under the HPM and the MEP as observed by them with brown dwarf stars. Instability has also been observed in the extension of the work of Bhar et al. [33], related to highly compact stars such as neutron stars and strange stars. In the present model we also observe the same feature of instability and may conclude that the applicability of the HPM and the MEP to radiating and highly compact stars may have some intrinsic problem to provide the entire physical validity of the models.