Expansion of magnetic neutron stars in an energy (in)dependent spacetime
Abstract
Regarding the strong magnetic field of neutron stars and high energy regime scenario which is based on high curvature region near the compact objects, one is motivated to study magnetic neutron stars in an energy dependent spacetime. In this paper, we show that such strong magnetic field and energy dependency of spacetime have considerable effects on the properties of neutron stars. We examine the variations of maximum mass and related radius, Schwarzschild radius, average density, gravitational redshift, Kretschmann scalar and Buchdahl theorem due to magnetic field and also energy dependency of metric. First, it will be shown that the maximum mass and radius of neutron stars are increasing function of magnetic field while average density, redshift, the strength of gravity and Kretschmann scalar are decreasing functions of it. These results are due to a repulsivelike force behavior for the magnetic field. Next, the effects of the gravity’s rainbow will be studied and it will be shown that by increasing the rainbow function, the neutron stars could enjoy an expansion in their structures. Then, we obtain a new relation for the upper mass limit of a static spherical neutron star with uniform density in gravity’s rainbow (Buchdahl limit) in which such upper limit is modified as . In addition, stability and energy conditions for the equation of state of neutron star matter are also investigated and a comparison with empirical results is done. It is notable that the numerical study in this paper is conducted by using the lowest order constrained variational (LOCV) approach in the presence of magnetic field employing AV18 potential.
I Introduction
Neutron stars with highdensity baryonic matter and in the presence strong magnetic fields are of considerable interests in astrophysics. Interior strong magnetic field of a neutron star can form from the compression of constant magnetic flux. This magnetic field can be as large as in the core of high density inhomogeneous gravitationally bound neutron stars (1). By comparing with the observational data, the magnetic field strength of neutron stars is obtained about (2). Besides, in the quark core of the high density neutron stars, the maximum field may reach up to about (1); (3). On the other hand, in order to study the compact objects such as neutron and quark stars, one has to consider the effects of their curvature in the context of general relativity (GR). The first hydrostatic equilibrium equation (HEE) for stars was obtained by Tolman, Oppenheimer, and Volkoff (TOV) (4); (5); (6). Employing TOV equation for investigating the properties of compact objects has been done by many authors (7); (8); (9); (10); (11); (12); (13); (14).
Strong magnetic fields have important influences on the structure of compact objects (15); (16); (17); (18); (19); (20); (21); (22); (23); (24). Using magnetohydrodynamic simulation in full GR, the stability of neutron stars with toroidal magnetic fields has been studied (17). It is found that nonrotating neutron stars are dynamically unstable for some value of magnetic field strength. Considering the equation of state (EoS) of the nuclear matter within a relativistic model subjected to a strong magnetic field, it has been shown that, the influence of magnetic field can increase the mass by more than for stars with hyperons interior (19). Employing an effective field theory model for the nuclear EoS and by considering two magnetic field geometries (tangled field and a forcefree field), Kamiab et al (25) have shown that these fields have specific effects on the mass of neutron stars, so that it leads to increasing of the maximum mass more that . The effects of a strong magnetic field on the nuclear EoS has been studied in (26). Cardall et al. (16) have investigated static neutron stars by considering polodal magnetic fields with a simple class of electric current distribution. They have used the nuclear EoS evaluated by Broderick et al. in Ref. (26), and found that the magnetic field increases the maximum mass of neutron stars, Astashenok et al. (27) studied neutron stars in the context of F(R) gravity. They have used a nuclear EoS which was similar to the one extracted by Broderick et al. (26), and by assuming the global effect of magnetic field pressure, they found that the magnetic field can also increase the maximum mass of neutron star. The effects of strong magnetic fields on the structure of neutron stars in a perturbative gravity and also a quantum hydrodynamics model have been investigated in (21). It is shown that there exists a set of parameters for the modified gravity and magnetic field strength in which even a soft EoS can accommodate a large maximum neutron star mass through the modified massradius relation. In the framework of GR and relativistic mean field theory, employing the coupled EinsteinMaxwellThomasFermi equations, a model for both static and uniformly rotating neutron stars has been developed (22). It is believed that the use of realistic parameters of rotating neutron stars obtained from numerical integration of the selfconsistent axisymmetric general relativistic equations of equilibrium leads to the values of magnetic field and radiation efficiency of pulsars very different from estimates based on fiducial parameters.
In solar system scale or other moderate gravitational regimes, GR is an effective theory with accurate results, whereas in context of strong gravity (near the compact objects), this theory has problems various in which we address some of them. From the cosmological point of view, the usual GR cannot explain the accelerating expansion of the universe (28); (29); (30). On the other hand, in order to construct a quantum theory of gravity, we have to evaluate the validity of theory in the ultraviolet (UV) regime. It was shown that GR is valid in infrared (IR) limit while in UV regime, it needs to improve for producing accurate results. Therefore, in order to include the UV regime, we have to modify GR. HoravaLifshitz gravity (31) and gravity’s rainbow (32) are two theories which include the UV regime by considering a modification of the usual energymomentum relation. Among these two theories, gravity’s rainbow has more interesting results. Garattini and Saridakis (33) have shown that considering a suitable choice of the rainbow functions, the HoravaLifshitz gravity can be related to the gravity’s rainbow.
Gravity’s rainbow can be originate from the generalization of double special relativity to the curved spacetime. The doubly special relativity (34) is the special relativity when we consider another upper limit constant to this theory. Indeed, in doubly special relativity there are two fundamental constants; the velocity of light (), and the Planck energy (). In doubly special relativity, a particle can not attain velocity and energy larger than and , respectively ( and , in which and are, respectively, the velocity and the energy of a test particle). In this theory, the metric specifying spacetime is constructed by employing the spectrum of energy that particles can obtain (35). Generalization of doubly special relativity in the presence of gravity leads to double GR or gravity’s rainbow (32). It is notable that, the gravity’s rainbow reduces to GR in the IR limit, in which the energy of particles has no direct effect on the curvature.
It was pointed out that although gravity’s rainbow modifies thermodynamical aspects of the black objects (36); (37); (38), the uncertainty principle is valid (39); (40). The effects of this gravity on information paradox posed by Hawking radiation and black holes evaporation have been studied in Ref. (41). Also, this theory has been investigated in context of black holes with different configurations (42); (43). Another interesting result of this gravity is related to early universe. Indeed by considering the energy dependent FriedmannRobertsonWalker metric, one can study the early universe without big bang singularity (44); (45); (46); (47). Therefore, it is important to investigate the structure of other compact objects such as neutron and quark stars in this theory.
In the present work, firstly, we want to study the effects of strong magnetic field on the structure of neutron star in Einstein gravity. Then, by generalizing Einstein gravity to gravity’s rainbow in the present of strong magnetic fields, we will simultaneously investigated the effects of these two factors (magnetic field and rainbow function) on the structure of neutron star matter.
Ii Structure properties of magnetized neutron stars
In this paper, we are going to investigate the properties of the magnetized neutron stars using the TOV equation and the EoS introduced in Appendix (see Fig. 7 in Appendix B). The TOV equation in Einstein gravity is given by (4); (5); (6)
(1) 
where , also and are, respectively, the pressure and density of the fluid which are measured by local observer. in the above equation is the gravitational constant.
Selecting central density under the boundary conditions , , we integrate the TOV equation outwards to radius , at which vanishes. This yields the radius and mass (see Ref. (48)). The effects of magnetic field on the gravitational mass and radius of neutron star are presented in table 1. Also, the effects of magnetic field on the diagrams related to masscentral mass density () and massradius () relations are presented in Fig. 1. We found that the maximum mass and radius of magnetized neutron stars increase when the magnetic field grows. In order to more investigation the effect of magnetic field on neutron star, we study other interesting properties of our solution.
One of interesting quantities of neutron star is related to the average density of the magnetized neutron star can be written as
(2) 
the obtained results of average density for different magnetic fields are interesting. According to table 1 one finds that the average density decreases when magnetic field increases. In other words, growing magnetic field, unlike gravitational attraction, leads to more expansions of neutron stars.
Another important quantity which is related to the compactness of a spherical object. It can be defined as the ratio of the Schwarzschild radius to the radius of object,
(3) 
which may be interpreted as the strength of gravity. Also the Schwarzschild radius is . Another quantity which shows the strength of gravity is related to the spacetime curvature. In Schwarzschild metric, the components of Ricci tensor () and the Ricci scalar () are zero outside the star and these quantities do not give us any information about the spacetime curvature. Therefore, we use another quantity in order to more investigation of the curvature of spacetime. The quantity which can help us to understand the curvature of spacetime is the Riemann tensor. According to this fact the Riemann tensor may have more components and for simplicity we can study the Kretschmann scalar for measuring of the curvature in vacuum. Therefore, the curvature at the surface of a neutron star in Einstein gravity is given as (49); (50)
(4) 
The numerical results confirm that by increasing the magnetic field, the strength of gravity decreases and this may increase the radius of these stars.
Another known parameter of neutron stars is the gravitational redshift which can be written in the following form
(5) 
as one can see, increasing , leads to decreasing the redshift.
The dynamical stability of the stellar model against the infinitesimal radial adiabatic perturbation was introduced by Chandrasekhar (51). This stability condition was developed and applied to astrophysical cases by many authors (52); (53); (54); (55). The dynamical stability condition is satisfied when the adiabatic index, , is larger than , i.e, , everywhere within the isotropic star. The adiabatic index is defined in the following form
(6) 
In order to investigate the dynamical stability of magnetized neutron star, we plot the adiabatic index versus the radius in Fig. 2. As one can see, these stars enjoy interior dynamical stability.
Here, we want to investigate the upper mass limit of a static spherical neutron star with uniform density in GR, socalled Buchdahl theorem (56). The GR compactness limit is given by (56)
(7) 
in which the upper mass limit is . Our numerical results confirm that the obtained mass of magnetic neutron stars in GR respects Eq. (7).
Studying Fig. 7 (Appendix B), shows that the core of neutron star is more rigid when the magnetic field increases, whiles the radius of magnetic neutron star is an increasing function of magnetic field, but the average density, compactness, redshift and Kretschmann scalar are decreasing functions of it. Therefore, by increasing the magnetic field, although core of the neutron star remains compact, the outer layers expand. As a results, the magnetic field has repulsive property.
Theoretically, it has been shown that the magnetic field contributes directly to the pressure in the Einstein equations in which pressure fits with square magnetic field () , see Refs. (24); (27); (57); (58) for more details. This shows that when magnetic field is increased, associated magnetic pressure, due to its repulsive nature, increases the radius of neutron star. In other words, by adding the magnetic field to our system, its pressure increases. Such increasing in the pressure leads to increasing the repulsive force and as a result the system will be unstable. In order to have a stable case, the mass of the neutron star must increase. Therefore, it is a direct effect that when magnetic field is increased, the total mass of system is increased, too.
In the next section, we want to evaluate the effects of gravity’s rainbow on the structure of neutron stars.
Iii Structure properties of neutron stars in gravity’s rainbow
Here we introduce the spherical symmetric metric in gravity’s rainbow as
(8) 
in which and are rainbow functions. In order to study the structure of neutron stars in gravity’s rainbow, firstly we should obtain the HEE in this gravity. The HEE in gravity’s rainbow by considering the metric (8) and the following field equation
(9) 
where and are the Einstein and energymomentum tensors, respectively, is given by (59)
(10) 
where
(11) 
Here, we consider an upper limit for energies which (, which and are particle and Planck energies, respectively. See Refs. (32); (33); (34); (35); (36); (37); (38) for more details). For , the equation (10) yields to the usual TOV equation (4); (5); (6) (see Eq. (1)).
It is notable that there are known three cases for the rainbow function. The first one is related to the hard spectra from gammaray bursts (60), with the following form
(12) 
The second one is motivated by studies conducted in loop quantum gravity and noncommutative geometry (61); (62) as
(13) 
The third one of energy function is due to the consideration of constancy of the velocity of light (63)
(14) 
where in the above models , and are constants which can be adjusted by experiment. Here, we are going to investigate more properties of the neutron stars by using the obtained HEE in gravity’s rainbow. As it was shown, the obtained HEE for gravity’s rainbow only depends to rainbow function of . This shows that, the first model can not be useful for investigating the structure of neutron stars, because in this case the extracted HEE will yield Einstein’s HEE (TOV equation) and we do not see any effects of rainbow functions. We consider the third model for three reasons; i) more effectiveness; this model (third model) affects temporal and spatial coordinates, simultaneously, on the contrary, the effects of the first and second models are only on temporal coordinate and spatial coordinates, respectively. ii) low number of free parameters; there are two free parameters for the second model, whereas the third model has one free parameter. iii) phenomenologically logical; third set of rainbow functions comes from the constancy of light velocity, which we would like to respect it. The results are collected in table 2. As one can see, by increasing rainbow function, the maximum mass and the corresponding radius of neutron star increase. This emphasizes the contributions of the gravity’s rainbow on properties of neutron stars.
The average density of neutron star in this gravity can be written as
(15) 
Table 2 shows that the average density decreases when increases. The Schwarzschild radius in this gravity is given as (59)
(16) 
the Schwarzschild radius increases when increases.
We obtain the compactness of a spherical object (), by using the obtained Schwarzschild radius (Eq. (16)) and Eqs. (3) and (11), as
(17) 
The gravitational redshift in gravity’s rainbow can be written in the following form (59)
(18) 
The above equations show that the compactness of a spherical object and the gravitational redshift depend on both the rainbow function () and the radius of neutron star (). On the other hand, numerical calculations showed that increasing the rainbow function leads to increasing the radius of neutron star, too. So, based on numerical calculations, we have concluded that changing the rainbow function does not have any significant role in and up to our numeric precision (see table. 2). Also, it was shown that, these stars have the dynamical stability (see Figs. 5 and 6 of Ref. (59), for more details).
Here we want to obtain the curvature at the surface of a neutron star in gravity’s rainbow. Considering the metric (8) with field equation (9), we can obtain the Kretschmann scalar as
(19) 
It is notable that when , the Kretschmann scalar of Einstein gravity (Eq. 4) will be recovered. The results show that, by increasing rainbow function which leads to increasing of radius and mass of neutron stars, but the strength of gravity decreases.
Another important quantity which is interesting for evaluation is related to the upper mass limit of a static spherical neutron star with uniform density in gravity’s rainbow. For this purpose, we have to obtain Buchdahl theorem in an energy dependent spacetime. Solving the integral given in Eq. (11), we obtain
(20) 
Integrating from a central pressure , we have
(22) 
in which the integral on the righthand side of the above equation is
(23) 
which we have used of in the above equation. On the other hand, the lefthand side of Eq. (22) is
(24) 
Combining Eqs. (23) and (24), we can write
(25) 
so we have
(26) 
in which by using Eq. (20), we can rewrite it as
(27) 
In order to obtain an expression for finding how the radius of star depends on its mass and density, we use the fact that the pressure of a star on its surface vanishes, i.e., for (). So we can obtain
(28) 
Introducing a new dimensionless variable as , we can rewrite Eq. (28) with the following form
(29) 
where . According to this fact the variable is functions of and in which both them can never be negative, hence can not be negative and also due to the negative sign of the , Eq. (29) reaches its maximum value when reaches its minimum value. In order to find the minimum of the , we find the derivative as . This is always positive. Thus increases with increasing , and the minimum of is located at . Inserting in Eq. (29), we have
(30) 
The equation (30) gives the value for the maximum mass of a star of given radius in gravity’s rainbow. As a results of our calculation for finding the upper mass limit for a static spherical neutron star with uniform density in gravity’s rainbow, the obtained mass of neutron stars we must be in the range
(31) 
It is notable that, for , this limit reduces to Einstein case (). Our results of mass of magnetic neutron stars in gravity’s rainbow show that these masses are in the obtained range of Eq. (31).
The maximum mass and radius of neutron star are increasing functions of the rainbow function while the average density and the Kretschmann scalar were a decreasing function of it. These properties indicate that increasing the rainbow function results into expansion of the neuron star.
Iv Structure properties of magnetized neutron stars in gravity’s rainbow
Considering the EoS introduced in Appendix A and B, we want to calculate the properties of magnetic neutron star in gravity’s rainbow in this section. The effects of rainbow function and magnetic field on the diagrams related to and relations are presented in Figs. 3 and 4. As one can see, by increasing rainbow function and magnetic field, the maximum mass and radius of magnetized neutron stars increase. It is notable that for , the maximum mass of these stars reduces to the results obtained in the Einstein gravity (see table 1), as expected. By increasing rainbow function () and magnetic field, the maximum mass of these stars will be larger than the case of the Einstein gravity (see table 3 and Figs. 3 and 4). This emphasizes the contributions of the gravity’s rainbow and magnetic field on properties of neutron stars. The effects of magnetic fields on the gravitational masses of neutron stars and MR relation for a fixed amount of the rainbow function are presented in Fig. 5. The maximum mass and radius of magnetized neutron star for each value of the rainbow function increase when the magnetic field grows (see table 3). Our results cover the mass measurement of massive magnetized neutron stars. For example; about for Vela X1 (64), PSR J16142230 (65) about , PSR J0348+0432 (66) about , 4U 1700377 (67) about , and J17482021B (68) about .
Previously, Rhoades and Ruffini have shown that by employing the principles of causality and Le Chantelier in Einstein gravity, the mass of neutron star can not be larger than (69). Here, it was shown that the principles of causality and Le Chantelier are valid for the EoS in this paper (see Appendix B). In addition, we find that depending on choices of rainbow function, the mass of the neutron star could be larger than .
The results for the Schwarzschild radius show that for each value of the magnetic field with stronger gravity’s rainbow, the Schwarzschild radius has higher value. In addition, with a fixed amount of rainbow function, the Schwarzschild radius grows when the magnetic field increases. Moreover, the average density of the magnetized neutron star shows that the average density decreases when both and increase.
The compactness of this star confirms when increases, the strength of gravity is the same, but by increasing the magnetic field, the strength of gravity decreases and this may increase the radius of this star.
Using Eq. (19), we find that for a fixed value of the magnetic field by increasing rainbow function, the strength of gravity of a neutron star has smaller values. In addition, with a fixed amount of rainbow function, the Kretschmann scalar decreases when the magnetic field increases.
It is notable that, by increasing leads to decreasing the gravitational redshift. On the other hand, considering a fixed value of magnetic field, the rainbow function does not affect the gravitational redshift. Our results show that, the maximum redshift for this star is about , which is lower than the upper bound on the surface redshift for a subluminal equation of state, i.e (70).
In order to investigate the dynamical stability of magnetized neutron star, we plot the adiabatic index versus the radius in Fig. 6. As one can see, these stars have the dynamical stability inside.
As previously mentioned, the pressure can be written as . In other words, by adding the magnetic pressure () to the pressure of the fluid (), the pressure power increases. On the other hand, in order to have a stable magnetic neutron star, the gravitational force has to increase and this leads to increasing of mass of magnetic neutron star.
V Comparison between theory and observations
One of the interesting features of our calculations is related to comparison of the theory and its predictions with the empirical evidences. For this purpose, we compare our results with the observational data. We have presented this comparison in table 4. It is clear that the results derived in this theory are in agreement with the results obtained through observations. It is notable that for nonmagnetized neutron stars, the radius obtained through gravity’s rainbow are smaller than the observed values. However, the radius of the magnetized neutron stars are in the range of the observed values. This is due to the fact that we have considered the magnetic field of neutron stars in our calculations (see table 4 for more details).

It is worthwhile to mention that obtaining radius with the observational measurement has specific difficulties. However, with using this theory, we can obtain the radius of neutron stars. The results of the theoretical calculations are presented in table 4. For observational cases, the theory under consideration with strong magnetic field predicts a set of radii for some cases such as J17482021B and 4U 1700377, which are about (km) and (km), respectively.
Vi Conclusions
The paper at hand studied the structure of magnetized neutron stars through the use of LOCV method and potential employing the hydrostatic equilibrium equation in Einstein gravity. It was shown that while the maximum mass and radius of neutron star are increasing functions of the magnetic field, the average density, compactness, redshift and Kretschmann scalar are decreasing functions of it. This results into expansion of the neutron star in case of increasing power of the magnetic field. In some senses, the magnetic field may have repulsive force property.
Next, the effects of the gravity’s rainbow were explored. It was shown that the compactness and the gravitational redshift are not sensitive to variation of the rainbow function, the maximum mass and radius of neutron star were increasing functions of it and the average density and Kretschmann scalar were decreasing functions of it. Therefore, increasing the rainbow function leads to expansion of the neutron stars. Furthermore, it was also pointed out that it is possible to have magnetized neutron stars with mass larger than in gravity’s rainbow.
Finally, a comparison with the observational data was done and it was shown that our theoretical results are in agreement with the empirical evidences of neutron stars.
Briefly, we obtained the quite interesting results for magnetic neutron star in an energy dependent spacetime such as: I) EoS derived from microscopic calculations satisfied Le Chatelier’s principle and also both energy and stability conditions, simultaneously (see Appendix B). II) The maximum mass and radius of neutron star were increasing functions of magnetic field and rainbow function. III) The obtained maximum mass of magnetic neutron star can be more than , due to considering an energy dependent spacetime. IV) Magnetic neutron star in gravity’s rainbow was dynamically stable. V) Obtained results of magnetic neutron stars in gravity’s rainbow are consistent with all the measured masses of pulsars and neutron stars such as; Vela X1 (about ) (64), PSR J16142230 (about ) (65), PSR J0348+0432 (about ) (66), 4U 1700377 (about ) (67) and J17482021B (about ) (68). VI) We have shown the magnetic field and rainbow function may have repulsive force property. VII) We obtained the Kretschmann scalar in gravity’s rainbow as a characteristic of gravity strength and evaluated the its behavior. VIII) We have extracted a new relation for the maximum mass of a star of given radius in gravity’s rainbow as .
It will be worthwhile to study the effects of magnetic field on properties of other compact objects such as quark star and white dwarf in the context of gravity’s rainbow. In this paper we considered a neutron star in static spherical symmetric spacetime and studied simultaneously the effects of magnetic field and rainbow function on its properties. Generalization of static compact objects to anisotropic (71); (72); (73); (74); (75), rotating (76); (77); (78); (79); (80); (81); (82), rapidly rotating (83); (84); (85); (86); (87) of compact objects and investigating the effect of magnetic field in gravity’s rainbow can be interesting topics. We leave these issues for the future works.
Acknowledgements.
We wish to thank Shiraz University Research Council. This work has been supported financially by the Research Institute for Astronomy and Astrophysics of Maragha, Iran. Appendix: A LOCV formalism for magnetized neutron matter We suppose a pure homogeneous system composed of spinup and spindown neutrons with the number densities of and , respectively. We denote the spin polarization parameter by and the total density of system by . We use LOCV method to calculate the energy of our system (88); (89). Although in order to obtain accurate results we should regard temperature of our system, here we consider zero temperature for simplicity. A trial manybody wave function of the form is considered where is the uncorrelated groundstate wave function of independent neutrons, and is a proper body correlation function. Applying Jastrow approximation ?, is replaced by , where is a symmetrizing operator. In addition, a cluster expansion of the energy functional up to the twobody term, , is considered. Taking the uniform magnetic field along the direction, , the spin up and down particles correspond to parallel and antiparallel spins with respect to the magnetic field. The energy per particle up to the twobody term iswhere and are onebody and twobody energy terms, respectively. The operator is nuclear potential which has been given in Ref. (88), and is the neutron magnetic moment (in units of the nuclear magneton). In the next step, we minimize the twobody energy with respect to the variations in the trial manybody wave function subject to a normalization constraint presented in Ref. (88). The result of minimization of the twobody cluster energy is a set of differential equations (88). Solving the differential equations leads to the twobody energy of the system.
Appendix: B EoS of magnetized neutron star matter
From the energy of magnetized neutron matter, we have to consider a suitable EoS. Here we evaluate the pressure using the following relation
Our results for the EoS of magnetized neutron matter for different values of the magnetic field have been presented in Fig. 7. It confirms that the EoS of magnetized neutron matter becomes stiffer at higher magnetic fields which is due to the inclusion of neutron anomalous magnetic moments.
In order to do more investigations of the EoS introduced in this paper, we analyze both energy and stability conditions. We examine the energy conditions, such as null energy condition (NEC), weak energy condition (WEC), strong energy condition (SEC), and also dominant energy condition (DEC) at the center of neutron stars. The requirements of each energy conditions can be summarize as
(32)  
(33)  
(34)  
(35) 
in which and are, respectively, the density and pressure at the center () of a magnetized neutron star. Using Fig. 7 and the above conditions (3235), we obtain the results related to the energy conditions in table 5.
Our results show that all energy conditions are satisfied (see table 5). So, this EoS of magnetized neutron star matter is suitable to study the properties of magnetized neutron stars. In order to examine the stability of EoS for a physically acceptable model, one expects that the velocity of sound, , be less than the light’s velocity, , (90); (91). Considering the relation, , we have presented the sound speed versus density in Fig. 8. It is evident that for all magnetic fields, the EoS of magnetized neutron star matter satisfies the inequality (see Fig. 8). Furthermore Le Chatelier’s principle is another stability condition that should be fulfilled. The establishment of this condition indicates that both as a whole and also with respect to the nonequilibrium elementary regions with spontaneous contraction or expansion, the object is stable (69); (92). It is evident from Fig. 7 that Le Chatelier’s principle is satisfied. Therefore, the EoS of magnetized neutron star matter satisfied both energy and stability conditions.
Footnotes
 email address: behzad.eslampanah@gmail.com
 email address: ghbordbar@shirazu.ac.ir
 email address: hendi@shirazu.ac.ir
 email address: ruffini@icra.it
 email address: zrezaei@shirazu.ac.ir
 email address: rahim.moradi@icranet.org
References
 E. J. Ferrer, V. de la Incera, J. P. Keith, I. Portillo and P. L. Springsteen, Phys. Rev. C 82, 065802 (2010).
 Y. F. Yuan and J. L. Zhang, Astron. Astrophys. 335, 969 (1998).
 T. Tatsumi, Phys. Lett. B 489, 280 (2000).
 R. C. Tolman, Proc. Nat. Acad. Sc. 20, 169 (1934).
 R. C. Tolman, Phys. Rev. 55, 364 (1939).
 J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 374 (1939).
 R. R. Silbar and S. Reddy, Am. J. Phys. 72, 892 (2004).
 G. Narain, J. SchaffnerBielich and I. N. Mishustin, Phys. Rev. D 74, 063003 (2006).
 G. H. Bordbar, M. Bigdeli and T. Yazdizade, Int. J. Mod. Phys. A 21, 5991 (2006).
 P. Boonserm, M. Visser and S. Weinfurtner, Phys. Rev. D 76, 044024 (2007).
 X. Li, F. Wang and K. S. Cheng, JCAP 10, 031 (2012).
 A. M. Oliveira, H. E. S. Velten, J. C. Fabris and I. G. Salako, Eur. Phys. J. C 74, 3170 (2014).
 X. T. He, F. J. Fattoyev, B. A. Li and W. G. Newton, Phys. Rev. C 91, 015810 (2015).
 N. Yunes and M. Visser, Int. J. Mod Phys. A 18, 3433 (2003).
 J. S. Heyl and L. Hernquist, Mon. Not. R. Astron. Soc. 324, 292 (2001).
 C. Y. Cardall, M. Prakash and J. M. Lattimer, Astrophys. J. 554, 322 (2001).
 K. Kiuchi, M. Shibata and S. Yoshida, Phys. Rev. D 78, 024029 (2008).
 R. Ciolfi, V. Ferrari and L. Gualtieri, Mon. Not. R. Astron. Soc. 406, 2540 (2010).
 L. L. Lopes and D. P. Menezes, Braz. J. Phys. 42, 428 (2012).
 C. Chirenti and J. Skakala, Phys. Rev. D 88, 104018 (2013).
 M. K. Cheoun, C. Deliduman, C. Güngör, V. Keles, C. Y. Ryu, T. Kajino and G. J. Mathews, JCAP 10, 021 (2013).
 R. Belvedere, J. A. Rueda and R. Ruffini, Astrophys. J. 799, 23 (2015).
 D. Neilsen, S. L. Liebling, M. Anderson, L. Lehner, E. O’Connor and C. Palenzuela, Phys. Rev. D 89, 104029 (2014).
 L. L. Lopes and D. P. Menezes, JCAP 08, 002 (2015).
 F. Kamiab, A. E. Broderick and N. Afshordi, arXiv:1503.03898
 A. E. Broderick, M. Prakash and J. M. Lattimer, Astrophys. J. 537, 351 (2000).
 A. V. Astashenok, S. Capozziello and S. D. Odinstov, Astrophys. Space Sci. 355, 333 (2015).
 S. Perlmutter et al., Astrophys. J. 517, 565 (1999).
 S. Perlmutter, M. S. Turner and M. White, Phys. Rev. Lett. 83, 670 (1999).
 A. G. Riess et al., Astrophys. J. 607, 665 (2004).
 P. Horava, Phys. Rev. Lett. 102, 161301 (2009).
 J. Magueijo and L. Smolin, Class. Quantum Grav. 21, 1725 (2004).
 R. Garattini and E. N. Saridakis, Eur. Phys. J. C 75, 343 (2015).
 G. AmelinoCamelia, Nature. 418, 34 (2002).
 J. J. Peng and S. Q. Wu, Gen. Relativ. Gravit. 40, 2619 (2008).
 A. F. Ali, Phys. Rev. D 89, 104040 (2014).
 S. H. Hendi, M. Faizal, B. Eslam Panah and S. Panahiyan, Eur. Phys. J. C 76, 296 (2016).
 S. H. Hendi, S. Panahiyan, B. Eslam Panah and M. Momennia, Eur. Phys. J. C 76, 150 (2016).
 Y. Ling, X. Li and H. Zhang, Mod. Phys. Lett. A 22, 2749 (2007).
 H. Li, Y. Ling and X. Han, Class. Quantum Grav. 26, 065004 (2009).
 A. F. Ali, M. Faizal and B. Majumder, EPL 109, 20001 (2015).
 P. Galan and G. A. Mena Marugan, Phys. Rev. D 74, 044035 (2006).
 Y. Gim and W. Kim, JCAP 05, 002 (2015).
 A. Awad, A. F. Ali and B. Majumder, JCAP 10, 052 (2013).
 S. H. Hendi, M. Momennia, B. Eslam Panah and M. Faizal, Astrophys. J. 827,153 (2016).
 Y. Ling, JCAP 08, 017 (2007).
 S. H. Hendi, M. Momennia, B. Eslam Panah and S. Panahiyan, Phys. Dark Universe 16, 26 (2017).
 S. Shapiro and S. Teukolsky, Black Holes, White Dwarfs and Neutron Stars. Wiley, New York (1983).
 D. Psaltis, Living Rev. Relativ. 11, 9 (2008).
 K. Y. Eksi, C. Gungor and M. M. Turkoglu, Phys. Rev. D 89, 063003 (2014).
 S. Chandrasekhar, Astrophys. J. 140, 417 (1964).
 J. M. Bardeen, K. S. Thonre and D. W. Meltzer, Astrophys. J. 145, 505 (1966).
 H. Kuntsem, Mon. Not. R. Astron. Soc. 232, 163 (1988).
 M. K. Mak and T. Harko, Eur. Phys. J. C 73, 2585 (2013).
 M. Kalam, S. M. Hossein and S. Molla, arXiv:1510.07015
 H. A. Buchdahl, Phys. Rev. 116, 1027 (1959).
 R. Moradi, C. Stahl, J. Firouzjaee and S. S. Xue, arXiv:1705.04168
 S. Ray, A. L. Espindola, M. Malheiro, J. P. S. Lemos and V. T. Zanchin, Phys. Rev. D 68, 084004 (2003).
 S. H. Hendi, G. H. Bordbar, B. Eslam Panah and S. Panahiyan, JCAP 09, 013 (2016).
 G. AmelinoCamelia, J. R. Ellis, N. Mavromatos, D. V. Nanopoulos and S. Sarkar, Nature. 393, 763 (1998).
 U. Jacob, F. Mercati, G. AmelinoCamelia and T. Piran, Phys. Rev. D 82, 084021 (2010).
 G. AmelinoCamelia, Living Rev. Relativ. 16, 5 (2013).
 J. Magueijo and L. Smolin, Phys. Rev. Lett. 88, 190403 (2002).
 M. L. Rawls et al., Astrophys. J. 730, 25 (2011).
 P. Demorest, T. Pennucci, S. Ransom, M. Roberts and J. Hessels, Nature. 467, 1081 (2010).
 J. Antoniadis et al. Science. 340, 6131 (2013).
 J. S. Clark et al., Astron. Astrophys. 392, 909 (2002).
 P. C. C. Freire et al., Astrophys. J. 675, 670 (2008).
 Jr. C. E. Rhoades and R. Ruffini, Phys. Rev. Lett. 32, 324 (1974).
 P. Haensel, J. P. Lasota and J. L. Zdunik, Astron. Astrophys. 344, 151 (1999).
 C. G. Boehmer and T. Harko, Class. Quantum Grav. 23, 6479 (2006).
 B. C. Paul and R. Deb, Astrophys. Space Sci. 354, 421 (2014).
 S. K. Maurya, Y. K. Gupta, S. Ray and B. Dayanandan, Eur. Phys. J. C 75, 225 (2015).
 D. K. Matondo and S. D. Maharaj, Astrophys. Space Sci. 361, 221 (2016).
 B. S. Ratanpal, V. O. Thomas and D. M. Pandya, Astrophys. Space Sci. 361, 65 (2016).
 N. Andersson and G. L. Comer, Class. Quantum Grav. 18, 969 (2001).
 Z. B. Etienne, Y. T. Liu and S. L. Shapiro, Phys. Rev. D 74, 044030 (2006).
 P. Pani and E. Berti, Phys. Rev. D 90, 024025 (2014).
 R. F. P. Mendes, G. E. A. Matsas and D. A. T. Vanzella, Phys. Rev. D 90, 044053 (2014).
 K. V. Staykov, D. D. Doneva, S. S. Yazadjiev and K. D. Kokkotas, JCAP 10, 006 (2014).
 A. Cisterna, T. Delsate, L. Ducobu and M. Rinaldi, Phys. Rev. D 93, 084046 (2016).
 J. G. Coelho, D. L. Cáceres, R. C. R. de Lima, M. Malheiro, J. A. Rueda and R. Ruffini, Astron. Astrophys. 599, A87 (2017).
 D. Lai and S. L. Shapiro, Astrophys. J. 442, 259 (1995).
 S. Yoshida, S. Karino, S. Yoshida and Y. Eriguchi, Mon. Not. R. Astron. Soc. 316, L1 (2000).
 D. D. Doneva, S. S. Yazadjiev, N. Stergioulas and K. D. Kokkotas, Phys. Rev. D 88, 084060 (2013).
 B. Kleihaus, J. Kunz, S. Mojica and M. Zagermann, Phys. Rev. D 93, 064077 (2016).
 F. Cipolletta, C. Cherubini, S. Filippi, J. A. Rueda and R. Ruffini, Phys. Rev. D 96, 024046 (2017).
 G. H. Bordbar, Z. Rezaei and A. Montakhab, Phys. Rev. C 83, 044310 (2011).
 Z. Rezaei and G. H. Bordbar, Eur. Phys. J. A 52, 132 (2016).
 L. Herrera, Phys. Lett. A 165, 206 (1992).
 H. Abreu, H. Hernandez and L. A. Nunes, Class. Quantum Grav. 24, 4631 (2007).
 N. K. Glendenning, Phys. Rev. Lett. 85, 1150 (2000).