Neutron Stars : A Comparative Study
Abstract
The inner structure of neutron star is considered from theoretical point of view and is compared with the observed data. We have proposed a form of an equation of state relating pressure with matter density which indicates the stiff equation of state of neutron stars. From our study we have calculated mass(M), compactness(u) and surface redshift() for the neutron stars namely PSR J16142230, PSR J1903+327, Cen X3, SMC X1, Vela X1, Her X1 and compared with the recent observational data. We have also indicated the possible radii of the different stars which needs further study. Finally we have examined the stability for such type of theoretical structure.
pacs:
24.10.Jv, 04.40.Dg, 26.60.Kp, 97.10.Nf, 97.10.PgI Introduction
Study of compact objects like star takes much attention to the astrophysicists for the last few Decades. It plays the crucial role like a bridge among astrophysics,
nuclear physics and particle physics. One of the possible existing compact objects in our universe is neutron star (or strange star). Neutron stars are the most
acceptable and known compact objects for the study of dense matter physics. In general, neutron stars are composed
almost entirely of neutrons, while strange stars are made from strange quark matter(SQM). If we consider strange quark matter is stable then strange stars may be formed
during supernova explosions. As a consequence of this neutron stars can be converted to strange stars by a number of different mechanisms such as:
a) pressureinduced transformation to udsquark matter
via udquark matter, b) sparking by highenergy neutrinos, c) triggering due to the intrusion of a quark nugget. All of these possibilities were described
by C. Alcock, E. Farhi, and A. Olinto(C. Alcock et al., 1986). Scientific community beleive that a much more rapid cooling of SQM in strange stars takes place than neutron stars
due to neutrino emitting weak interactions involving the quarks(C. Alcock et al., 1986). Therefore, a strange star was known as much colder than a neutron star of similar age.
The most massive compact stars are entirely made of quarks, i.e. they are quark stars(A. Drago, 2014; P. Haensel, 1986).
Famous scientist Edward Witten(E. Witten, 1984) concluded that “If quark matter is stable,
it is probably necessary to assume that ordinary neutron stars are really quark stars”.
Recently observed neutron star or strange star have been detected in the Milky Way and only very few of them have been discovered in globular clusters. Some of the possible
existing compact stars(neutron star) are PSR J16142230, PSR J1903+327, Cen X3, SMC X1, Vela X1,
Her X1(P B Demorest et al., 2010; P. C. C. Freire et al., 2011; M L Rawls et al., 2011; M. J. Coe et al., 2013). Many researcher
studied(Rahaman, 2012; Kalam, 2012a; Hossein, 2012; Rahaman, 2012; Kalam, 2012b, 2013, 2014a, 2014b, 2014c; Hossein, 2014; Lobo, 2006; Bronnikov, 2006; Egeland, 2007; Dymnikova, 2002)
compact stars in various directions. Scientists used different techniques such as computational, observational or theoretical analysis to study astrophysical objects.
Since the uncertainty in the behavior of matter inside the compact stars ( ”normal matter for neutron stars” or
“strange quark matter for strange stars”), their physical
structure and some physical property can be obtained by applicable analytic solutions of Einsteinâs
gravitational field equations. The Relativistic stellar model have been studied by Karl Schwarzschild in 1916(K. Schwarzschild, 1916), occasionally the first solution of
Einsteinâs field equation for the
interior of a compact object which are in hydrostatic equilibrium. In this regard, an important
investigation was done by Tolman (R. C. Tolman, 1939) in 1939. In his paper, he proposed eight evident analytical solutions of the field equations.
Clifford E. Rhoades, Jr. and Remo Ruffini(Clifford E. Rhoades, 1974) showed that maximum possible mass of a neutron star can not be larger
than 3.2.
Buchdahl(H. A. Buchdahl, 1959) also contributed an important study of
the fluid spheres by obtaining the famous bound to the
mass(M)  radius(R) ratio for stable general relativistic spheres, which is .
Inspite of considerable progress in current years, still there are some important features to be addressed about neutron stars. Observations of masses and radii
confirms the theoretical model for white dwarf stars. However, due to the uncertainities in the equation of state of neutron stars, observations of masses and radii are
used to test theories of neutron stars.
Motivated by the above facts we are introducing a theoretical model of neutron stars which indicates the possible equation of state, determine the masses and radii of various
neutron stars. The solution presented in this article satisfies all the energy conditions, TOVequation, and Buchdahl massradius relation(H. A. Buchdahl, 1959)
which is required for the stability of the stars.
It also satisfies other stability condition like speed of sound () inside the star and adiabetic index () in radial adiabatic perturbations.
In this article we are able to give an equation of state(EOS) relating pressure and matter density of
the star which is also a support for stiff equation of state for neutron star.
Interestingly, the star masses calculated from our model is well matched with the observed stars masses, which confirm the validity of our model for real case.
In Sec II, we have discuss the interior spacetime and behavior of the star. In Sec. III, we have studied some special features of the star namely, TOV equation,
Energy conditions, Stability, Massradius relation, Compactness and Surface redshift in different subsections. In Sec. IV, The article concluded
with a short discussion with numerical data and concluding remarks. Before we start it is worthy to mention that, we use geometric units throughout the article.
Ii Interior Spacetime of the star
In our stellar model, we consider a static and spherically symmetric matter distribution whose interior space time as:
(1) 
According to H. Heintzmann(H. Heintzmann, 1969),
(2) 
and
(3) 
where A (dimensionless), C(dimensionless) and (length) are constants.
The matter within the star is to be considered as perfect fluid (locally) and consequently the energy momentum tensor for the fluid distribution can be
written as
where is the four velocity, is the matter density and is the fluid pressure.
Solving the Einstein’s field equation we get,
and
where and are the central density and central pressure of the star respectively.
We observe that, pressure and density are maximum at the centre and decreases monotonically towards the boundary(FIG. 1 and FIG. 2). Therefore, they are well behaved in the interior of the stellar structure.
Iii Some special features
iii.1 TOV equation
The modified TolmanOppenheimerVolkoff (TOV) equation describes the equilibrium condition for the compact star subject to the gravitational() and hydrostatic() forces,
(4) 
where,
(5)  
(6) 
Therefore, the static equilibrium configurations do exist in the presence of gravitational and hydrostatic forces (FIG 3).
iii.2 Energy conditions
We verify whether all the energy conditions, like, null energy condition(NEC), weak energy condition(WEC), strong energy condition(SEC)
and dominant energy condition(DEC) are satisfied at the centre () of the star or not. For this purpose the following equations to be satisfied:
(km)  (km)  + (km)  3+ (km) 

0.000543987  0.000296278  0.000840265  0.001432821 
(i) NEC: ,
(ii) WEC: , ,
(iii) SEC: , ,
(iv) DEC: .
From FIG. 1 and FIG. 2 we observe that all the energy conditions are satisfied (see Table:I).
iii.3 Stability
For a physically acceptable model, one expects that the velocity of
sound should be less than light velocity i.e. within the range (Herrera, 1992; Abreu et al., 2007).
We plot the sound speed in FIG.4 and observe that it satisfy the inequality everywhere within the stellar object. This indicates the stability of the model.
We also investigate the dynamical stability of the stellar model against the infinitesimal radial adiabatic perturbation which was introduced by
S. Chandrasekhar(S. Chandrasekhar, 1964). Later this stability condition was devoloped and applied to astrophysical cases by
J. M. Bardeen, K. S. Thorne, and D. W. Meltzer(J. M. Bardeen et al., 1966), H. Knutsen(H. Knutsen, 1988), M. K. Mak and T. Harko(T. Harko, 2013), gradually.
The adiabatic index() is defined as
(7) 
The above equation can be written as (assume c=1)
(8) 
Since should be everywhere within the isotropic stable star, we plot the adiabetic index for our model (FIG. 5) and
observe that everywhere within the star. Therefore this stellar model is stable against the radial adiabatic infinitesimal perturbations.
iii.4 MassRadius relation and Surface redshift
According to Buchdahl (H. A. Buchdahl, 1959), for a static spherically symmetric perfect fluid sphere, maximum allowable massradius ratio should be . Mak.(Mak et al., 2001) also proposed more simplified expression. In our model, the gravitational mass (M) in terms of the energy density () can be expressed as
(9) 
where is the radius of the star.
The compactness(u) of the star should be
(10) 
The nature of the Mass and Compactness of the star are shown in FIG. 6 and FIG. 7.
The surface redshift () corresponding to the above compactness () can be written as
(11) 
where
(12) 
Therefore, from FIG. 8, the maximum surface redshift for the isotropic neutron stars of different radii can be obtain. Mass, compactness, surface redshift of the neutron stars are evaluated from the above equations ( 9, 10, 12) and a comparative analysis is done in Table II.
Star  Radius(in km)  Observed Mass()  Mass from Model()  Redshift  Compactness 
PSR J16142230  10.3  1.97 0.04  1.86932  0.252915  0.2677 
Vela X1  9.99  1.77 0.08  1.73007  0.23689  0.2554 
PSR J1903+327  9.82  1.667 0.021  1.65602  0.228382  0.2487 
Cen X3  9.51  1.49 0.08  1.52525  0.213363  0.2366 
SMC X1  9.13  1.04 0.09  1.37254  0.195798  0.2217 
Her X1  7.7  1.073 0.36  0.874387  0.137476  0.1675 
Iv Discussion and Concluding remarks
For a binary system, Jacoby et al.Jacoby et al. (2005) and Verbiest et al. Verbiest et al. (2008) used the detection of Shapiro delay to measure the masses of
both the neutron star and its binary companion. Using the same approach, Demorest et al.Demorest et al. (2010) has performed radio timing observations for the binary
millisecond pulsar PSR J16142230. The measured mass for the above pulsar in Demorest et al. (2010) is up to the highest precession reported yet.
With the help of Arecibo and Green Bank radio timing observations and considering the relativistic Shapiro delay, Freire et alP. C. C. Freire et al. (2011) obtained new constraints
on the mass of the pulsar and its companion and determine the accurate mass for PSR J1903+327 as through a detailed analysis. Rowls et al.M L Rawls et al. (2011)
have determine the mass of neutron stars such as (Vela X1, Cen X3, SMC X1, Her X1) in eclipsing Xray pulsar binaries. Their measured values are for Vela X1,
for Cen X3, for SMC X1 and for Her X1. We will restrict our discussions to these six stars only
though there are similar data for other stars also.
In this article, we have investigated the physical behavior of the neutron
stars by considering it as isotropic pressure in nature and the spacetime
of it to be described by metric Heint IIa(H. Heintzmann, 1969).
The observed masses and proposed radii for the neutron stars PSR J16142230, PSR J1903+327, Cen X3, SMC X1, Vela X1, Her X1 are used in this paper for an
effective model of neutron stars.
We have obtained quite interesting results from this model, which are as follows:

The solutions are regular at the origin.

Density and pressure variation at the interior of the neutron stars are well behaved(positive definite at the origin)[FIG.1 and FIG.2].

Pressure reduced to zero at the boundary but density shows some finite value which is physically acceptable.

In our model both the pressure and energy density are monotonically decreasing to the boundary.

It also satisfies TOV equation and energy conditions.

Our isotropic neutron star model satisfy Herrera’s stability condition (Herrera, 1992).

From the massradius relation (9, 10), all desired interior features of a neutron stars can be evaluated.

Our model satisfies Buchdahl massradius relation()(H. A. Buchdahl, 1959).

According to our model the surface redshift of the neutron stars is found within standard value i.e. which is satisfactory(Haensel et al., 2000).

This model is stable with respect to infinitesimal radial perturbations.

From FIG. 6 it is clear that our derived mass function is very much acceptable as the observed masses of various neutron stars are lying on the graph.

We also give two possible equations of state(EOS) of matter at the interior of the star, one is polynomial(FIG. 9) and another is exponential(FIG. 10). But we see that best fitted equation of state(EOS) is exponential.

We also estimated from our model that the equation of state of a neutron star would be like where are constant.

According to our model, equation of state (FIG 11) for a neutron star should not be soft equation of state rather it would be a stiff equation of state which is also evidently proved in a result of F. zel F. Ozel (2006).
Therefore our model is justified for the neutron star PSR J16142230, PSR J1903+327, Cen X3, SMC X1, Vela X1 and Her X1. Hence we can conclude with the hope that it may be helpful to analyse the other neutron stars also.
Acknowledgment
MK and SMH would like to thank IUCAA, Pune India and IMSc, Chennai, India for providing research facilities under Visiting Associateship.
References
 C. Alcock, E. Farhi, and A. Olinto, Astrophys. J. 310, 261 (1986).
 A. Drago, A. Lavagno, G. Pagliara, Phys.Rev. D89,043014 (2014).
 P. Haensel, J. Zdunik, R. Schaeffer, Astron.Astrophys. 160, 121 (1986)
 E. Witten, Phys. Rev. D 30, 272 (1984).
 P B Demorest et al, Nature, 467, 1081 (2010);
 P. C. C. Freire et al, MNRAS, 412, 2763 (2011) ;
 M L Rawls et al, ApJ, 730, 25(2011) ;
 M. J. Coe et al, MNRAS, 433, 746(2013) ;
 F. Rahaman et al., Gen. Relativ. Gravit. 44, 107 (2012).
 F. Rahaman et al., Eur. Phys. J. C 72, 2071 (2012).
 M. Kalam et al., Eur. Phys. J. C 72, 2248 (2012).
 Sk M. Hossein et al, Int. J. Mod. Phys. D 21, 1250088 (2012).
 M. Kalam et al., Int.J.Theor.Phys. 52, 33193328 (2013).
 M. Kalam et al., Eur. Phys. J. C 73, 2409 (2013).
 M. Kalam et al., Eur. Phys. J. C 74, 2971 (2014).
 M. Kalam et al., Astrophys. Space Sci. 349, 865 (2014)
 M. Kalam et al., arXiv: 1410.0199 [grqc].
 Sk M. Hossein et al., arXiv: 1408.2412 [grqc].
 F. Lobo, Class. Quantum. Grav. 23, 1525 (2006).
 K. Bronnikov and J.C. Fabris, Phys. Rev. Lett. 96, 251101 (2006).
 E. Egeland, Compact Star, Trondheim, Norway (2007).
 I. Dymnikova, Class. Quantum. Gravit. 19, 725 (2002).
 K. Schwarzschild, Sitzungsberichte der Koniglich Preussischen Akademie der Wissenschaften zu Berlin, Phys.Math. Klasse 23, 189 (1916).
 R. C. Tolman, Phys. Rev. 55, 364 (1939).
 Clifford E. Rhoades, Jr. and Remo Ruffini, Phys. Rev. Lett. 32, 324 (1974).
 H. A. Buchdahl, Phys. Rev. 116, 1027 (1959).
 H. Heintzmann, Z. Phys., 228, 489 (1969).
 Herrera L., Phys. Lett. A, 165, 206 (1992)
 Abreu H., Hernandez H. and Nunez L. A., Class. Quantum. Grav., 24, 4631 (2007)
 S. Chandrasekhar, Astrophys. J. , 140, 417 (1964)
 J. M. Bardeen, K. S. Thorne, and D. W. Meltzer, Astrophys. J. , 145, 505 (1966)
 H. Knutsen, MNRAS , 232, 163 (1988)
 M. K. Mak and T. Harko , Eur. Phys. J. C 73, 2585 (2013)
 Mak M. K., Dobson P. N. and Harko T., Europhys. Lett., 55, 310 (2001)
 B. A. Jacoby, A. Hotan, M. Bailes, S. Ord, S. R. Kulkarni, Astrophys. J. , 629, L113 (2005)
 J. P. W. Verbiest, M. Bailes, W. van Straten, G.B.Hobbs, R.T. Edwaris, R.N. Manchester, N.D.R. Bhat, J.M. Sarkissian, B.A. Jacoby, S.R. Kulkarni, Astrophys. J. , 676, 675 (2008)
 P.B. Demorest, T. Pennucci, S. M. ransom, M.S.E. Roberts, W.T. Hessels, Nature, 467, 1081 (2010)
 P. Haensel, J. P. Lasopa, J.L. Zdunik, Nucl. Phys. Proc. Suppl. 80, 1110 (2000)
 F. zel, Nature, 441, 1115 (2006)