Neutron Stars in Mimetic Gravity

In this paper, a modified version of the hydrostatic equilibrium equation based on the mimetic gravity in the presence of perfect fluid is revisited. By using the different known equation of states, the structural properties of neutron stars are investigated in general relativity and mimetic gravity. Comparing the obtained results, we show that, unlike general relativity, we can find the appropriate equation of states that support observational data in the context of mimetic gravity. We also find that the results of relativistic mean-field-based models of the equation of states are in better agreement with observational data than non-relativistic models.


I. INTRODUCTION
General Relativity (GR) and its black hole solutions are interesting topics in gravitating systems. Despite having acceptable information from the region beyond the event horizon of black holes, the interior solution of horizonless massive objects is one of the open questions in physics. Indeed, the interior properties of stars and their time evolution are of interest during recent decades. The Tolman-Oppenheimer-Volkov (TOV) equation [1,2] is one of the primitive attempts to describe the interior properties of spherical static perfect fluid objects. However, for solving the field equations and find the physical properties of a typical star, we have to consider an equation of state (EoS) explaining the relationship between two physical quantities; pressure and density.
We should emphasize that the interior solutions are based on three components: i) metric ansatz, ii) gravitational theory with an appropriate energy-momentum tensor such as GR with a perfect fluid which leads to the equation of hydrostatic equilibrium known as TOV, iii) EoS such as polytropic relation between pressure and density. Since the theoretical results of the TOV equation arisen from GR are not consistent with the observational evidence, different attempts have been considered to modified one (some) of the mentioned components. All these attempts are motivated by obtaining a better agreement between theory and observation, or testing alternative theories of gravitation. In this regard, different models of gravity are considered to investigate the compact objects and their modified TOV equations. For example neutron stars in an energy dependent spacetime are studied in Ref. [3]. Modified TOV equation in vectortensor-Horndeski theory and dilaton gravity are reported in Refs. [4] and [5], respectively. Moreover, the neutron star structure is investigated in the context of F (R), F (G), F (T ), F (R, T ) and F (R, L) theories of gravity [6][7][8][9][10][11][12][13][14][15][16][17][18][19]. Furthermore, neutron stars in massive gravity and magnetized neutron stars in gravity's rainbow are analyzed in [20] and [21], respectively. The structure of neutron stars in massive scalar-Gauss-Bonnet gravity is studied in Refs. [22][23][24]. Besides, neutron stars in the context of Eddington-inspired Born-Infeld theory of gravity are extracted in Refs. [25]. Considering the scalar-tensor theories of gravity, some interesting properties of neutron stars are studied in Refs. [26][27][28][29][30][31]. In addition, neutron star structure in Hořava and Hořava-Lifshitz theories of gravity are evaluated in Refs. [32] and [33], respectively. See Refs.  for additional information about neutron stars in modified gravity theories.
In this paper, we regard a modified TOV equation comes from the mimetic gravity in the presence of perfect fluid. As a matter of fact, mimetic gravity is a Weyl-symmetric extension of GR in which a dust-like perfect fluid can mimic cold dark matter from the cosmological point of view [62]. So, one of the main motivations of considering mimetic gravity is providing an interesting geometrical description for challenging topics such as the late-time acceleration [63], inflation [64] and dark matter [62]. Indeed, considering the mimetic scalar field in gravitating systems, one finds a dynamical longitudinal degree of freedom in addition to two traditional transverse degrees of freedom of GR. It is shown that such an additional degree of freedom could play the role of mimetic dark matter. The main idea of mimetic gravity backs to the interesting work of Chamseddine and Mukhanov [62], that proposed isolation of the conformal degree of freedom of gravity by introducing a parameterized physical metric g µν in terms of an auxiliary metricg µν and a (mimetic) scalar field φ, as follows g µν = −g µνg γδ ∂ γ φ∂ δ φ, (1.1) where confirms that the physical metric is invariant under conformal transformations of the auxiliary metric. Besides, taking a consistency condition into account, it is easy to show that the mimetic scalar field should satisfy the following constraint g γδ ∂ γ φ∂ δ φ = −1.
(1. 2) In this work, we follow the conventional notation of mimetic gravity in which this constraint can be appeared at the level of action formalism by including a Lagrange multiplier. Following Ref. [62], it is worth mentioning that the conformal degree of freedom becomes dynamical even in the absence of matter field, and therefore, the mimetic gravity may admit non-trivial solutions.
Motivated by what was mentioned above, the mimetic gravity has arisen a lot of attention in the past few years. From the cosmological point of view, it is a powerful theory to explain the flat rotation curves of spiral galaxies [65], gravitational wave [66,67] and the cosmological singularity [68]. This theory can resolve the singularity in the center of black holes [69] and introduce new black hole solutions from the gravitational viewpoint. Thus, it will be interesting to investigate its effect on astrophysical objects such as neutron stars and study their structures.
Without exaggeration, neutron stars are proper laboratories to examine fundamental physics at high energy regime; from strong gravitational field to considerable nuclear densities. Actually, neutron stars help us to investigate the high density low temperature regime of matter field that is complementary to terrestrial laboratories. The density of neutron stars is distributed from few g/cm 3 to more than 10 15 g/cm 3 at their surface and the center, respectively. The different layers of a neutron star can be categorized into the atmosphere, the outer and inner crusts, and also outer and inner cores. The outer crust can be described by comparing it with experimental data of atomic nuclei. Nevertheless, uncertainty in the EoS of neutron stars for regions with densities above nuclear matter density is still a substantial challenge since analyzing the astrophysical observations depends properly on the neutron star structure. So, there are different proposed EoSs that can describe the inner neutron stars with different structures such as the nucleon, hybrid, and the strange quark. The main theoretical techniques for determining EoSs are based on variational method, relativistic mean-field (RMF) models, effective interactions, perturbation expansion, effective energy-density functionals, and so on [70].
In this work, we consider the nuclear and hadronic stars which are in the range of the gravitational wave event GW170817 [71] and they have the maximum mass M max > 2M . For pure nuclear matter, we consider different EoSs: SLy4 [72] of the Lyon group that is based on energy density functional, WFF1 [73] that is obtained by the variational method. We also include BSK21 [74] and MPA1 [75] constructed from the generalized Skyrme nuclear interaction arising from the Argonne V 18 potential plus a microscopic nucleonic three-body force and Brueckner-Hartree-Fock theory, respectively. Besides, we choose the density-dependent relativistic mean-field (RMF) HS (DD2) EoS [76,77] and the recently proposed FSU2H [78][79][80]. For the hadronic stars, we use various hadronic EoSs: BHBΛΦ [81] and SFHoY [76,82,83].
In this paper, we consider the different observational evidence of mass-radius relation of neutron stars, such as GW 170817, PSR J0740 + 6620, PSR J2215 + 5135, NICER data on PSR J0030 + 0451 and GW190814. In fact, taking different observational data into account, we want to demonstrate how the mimetic gravity can be adapted to the appropriate EoS to obtain theoretical results consistent with observational evidence.
The outline of this paper is as follows. In the next section, we consider the spherically symmetric metric and obtain the required equations for solving, numerically, the generalized TOV in mimetic gravity. In section III, we introduce a set of nuclear matter EoSs and three potentials to drive M − R diagrams and discuss the properties of neutron stars. Finally, we finish our paper with conclusions.

II. HYDROSTATIC EQUILIBRIUM EQUATION IN MIMETIC GRAVITY
The starting point is the action of mimetic gravity in four dimension in the presence of a matter field as where R and φ are, respectively, the Ricci scalar and the mimetic scalar field, λ is the Lagrange multiplier while the constant = ±1 can be fixed by spacelike or timelike nature of ∂ µ φ. Besides, V (φ) is the potential related to the mimetic scalar field and I m denotes the action of matter field which we consider it as a perfect fluid. It is straightforward to vary the action (2.1) with respect to the metric tensor g µν and the mimetic scalar field φ to obtain the equations of motion as where κ = 8πG c 4 and G is the four dimensional gravitational constant. In addition, G µν is the Einstein tensor and T µν = −2 √ −g δIm δg µν is the energy-momentum tensor with the following explicit form where ρ and P are, respectively, the density and pressure of the perfect fluid measuring by a local observer, and u ν is the fluid four-velocity (u ν u ν = c 2 ). Moreover, one can vary the above action with respect to the Lagrange multiplier to obtain g µν ∂ µ φ∂ ν φ = . (2.5) Taking the trace of Eq. (2.2) and inserting Eq. (2.5), one finds Now, we can insert λ (Eq. (2.6)) into the field equations to achieve In order to find the equation of hydrostatic equilibrium for compact stars, we have to regard a suitable line element. Here, we assume a spherical symmetric spacetime with the following ansatz where we should analyze the properties of the radial metric functions f (r) and g(r). Using the metric introduced in Eq. (2.9), we can obtain the components of energy-momentum tensor (2.4) as Taking into account Eq. (2.5) with the metric ansatz (2.9), the mimetic scalar field can be given by the spatial metric function with the following form Inserting the metric ansatz (2.9) with Eqs. (2.10) and (2.11) into the gravitating field equation, we can simplify the components of Eq. (2.7) where f , g, ρ and P are functions of r and also the prime and double prime are, respectively, the first and second derivatives with respect to r. Taking Eq. (2.11) into account, we can remove φ in Eq. (2.13) as follows Considering Eq. (2.12) as a differential equation, we can obtain g(r) as a function of from ρ and V (φ) Now, we derive, algebraically, g (r) from Eq. (2.12) and regard g(r) from Eq. (2.16), and insert them into Eq. (2.15). After some manipulations, we can find that Eq. (2.15) can be rewritten as where M = 4πr 2 ρ dr. Here, we can define a new variable w = df dr to change Eq. (2.17) into the following first-order differential equation Now, we are in a position to talk about the equation of hydrostatic equilibrium by use of the conservation of energy-momentum tensor, ∇ µ T µν = 0, in which for ν = r, we can write Now, we have Eqs. (2.11), (2.12), (2.18), (2.19) and we should solve them numerically for a given EoS (P (ρ)).

III. EOS, NEUTRON STAR STRUCTURE AND OBSERVATIONAL DATA
There is a theoretical challenge to applying a proper EoS for the description of the internal structure of neutron stars. Here, we consider different known models of EoS in the context of GR and mimetic gravity to obtain an appropriate conclusion. So far, a large number of the EoS with various motivations have been introduced, and in this paper, we examine eight relativistic/non-relativistic EoSs that seem more realistic and have already received more attention.
The first non-relativistic EoS is SLy4 has been derived from energy density functional methods with an effective nucleon-nucleon interaction [72]. The second model of non-relativistic EoSs refers to BSK21 that is constructed from the generalized Skyrme nuclear interaction arising from the Argonne V 18 potential plus a microscopic nucleonic three-body force [74]. The third non-relativistic model is defined according to a variational method known as WFF1 [73]. Now, we point out the relativistic category of EoSs that we considered in this paper. The first model is devoted to MPA1 arising from the relativistic Dirac-Brueckner-Hartree-Fock approach with the π-ρ mesons exchange [75]. We also choose the HS (DD2) EoS [76,77] as the second model of this category that is based on the density-dependent relativistic mean-field (RMF) interactions of Typel et al. [84]. The third model is known as FSU2H EoS. The nucleonic FSU2 model of [85] is modified in the context of relativistic mean-field theory to describe both the nucleon and hyperon interactions, and so FSU2H EoS is created [78][79][80].
We also use two hadronic EoSs, BHBΛΦ [81] and SFHoY [76,82,83], as other realistic models. The BHBΛΦ EoS additionally includes Λ hyperons and hyperon-hyperon interactions mediated by Φ mesons. The Λ hyperon makes the EoS stiffer resulting in 2.1M maximum mass neutron star corresponding to a radius 11.58 km. Also, SFHoY is based on the statistical model with an excluded volume and interactions of Hempel and Schaffner-Bielich(HS) [76] with RMF interactions SFHo [83]. Although these EoSs are obtained with RMF interactions, we may exclude them from the full relativistic category since they are formulated in densities below saturation density and low temperatures. Before proceeding, it is notable that we examine the mentioned models of EoSs with different accepted functions of scalar potentials V (φ).
It is worth mentioning that in this paper, for nucleonic EoSs, we use from the parameterized piecewise-polytrope representation applied in [86][87][88], which we now describe in brief.
A piecewise polytropic EoS with four segment pieces satisfy the following polytropic relation, where P is the fluid pressure, and K i and Γ i denote, respectively, the polytropic constant and the adiabatic index. The first polytrope piece shows the crust EoS that is determined by K 0 = 35939 × 10 13 [cgs] and Γ 0 = 13572 [86]. Furthermore, two dividing high densities between core pieces are chosen as ρ 1 = 10 14.7 g/cm 3 and ρ 2 = 10 15 g/cm 3 . In this method, for a chosen of four free parameters P 1 , Γ 1 , Γ 2 , Γ 3 (where P 1 = P (ρ 1 ) and Γ 1 , Γ 2 and Γ 3 are the adiabatic indices) and by considering the continuity of the pressure, we can fix the polytropic constants K 1 , K 2 , and K 3 in the following form Considering the mentioned points and following Ref. [86], the parameterized EoS is completely determined (For the details see Ref. [86]). Moreover, taking the GR limit into account, one can find that all of the chosen EoSs are in the range of the gravitational wave event GW170817 [71] and they produce a neutron star with maximum mass M max > 2m . However, these EoSs did not support some of the observational evidence (see Fig. 1 and table I for more details).
Before comparing the results of neutron star properties in GR and mimetic gravity, we should note that in Figs. 1-5, we consider 1σ(%63) and 2σ(%95) mass-radius constraints from the gravitational wave event GW170817 by the blue and orange clouded regions corresponding to the heavier and lighter neutron stars, respectively [71]. In these figures, the purple region also shows the 1σ and 2σ confidential levels of the massive pulsar PSR J0740+6620 with mass 2.14 +. 20 −0.18 M [89]. Also, we illustrate the mass range of the pulsar PSR J2215+5135 [90], with a mass 2.27M , by light green region in the mentioned figures. Besides, regarding these figures, the black error bars and the red band show the NICER (Neutron Star Interior Composition Explorer) mass-radius measurement on PSR J0030+0451 [91,92] and GW190814 event [93], respectively.
Since we have the results of different EoSs in GR, here, we use the modified theory of gravity (mimetic gravity) to calculate the mass-radius relation for each EoS to find matching with the observational data. Strictly speaking, we investigate the mass-radius relation of neutron stars in mimetic gravity for three scalar potentials; e Kφ 2 , and V 3 (φ) = α3φ 2 1+e Kφ 2 . Therefore, we have the free parameters α i (i = 1..3) and φ(0) that can affect more or less the maximum mass and radius of the neutron stars (Note: φ(0) = φ(r)| r=0 is the initial value of mimetic scalar field at the center of the neutron star which should be fixed for numerical analysis).
The first step is devoted to considering one of the well-known EoSs, SLy4 EoS. It is easy to show that in the GR limit, this EoS can reach 2M and it is within the GW170817 region. However, the mass-radius relation of this EoS could not support one of the ranges of J0030+0451 error bars. For SLy4 EoS, we plot the M −R relation for potentials Fig. 2 with φ(0) = 1 (left panels) and φ(0) = 0.8 (right panels). We observe for φ(0) = 1, the maximum mass generally increases by increasing the parameter α i . Recall that the sign of α i is negative, and hereafter, we discuss the absolute value of α i for the sake of simplicity. The increase of the mass at V 3 (φ) is greater than V 1 (φ) and V 2 (φ). In order to have a easier comparison, we calculate the maximum mass and radius for three     different values for α i and two values for φ(0) in table II. Considering the reported maximum mass, one can find that the larger one belongs to potential V 3 (φ). Furthermore, according to Fig. 2, we can find that by increasing of |α i |, there is a small increment in the radius corresponding to potential V 2 (φ) and a bigger one corresponding to V 1 (φ) and V 3 (φ). Therefore, considering the potentials V 1 (φ) and V 3 (φ) for α i = −0.05, −0.10 (i = 1, 3) , not only the curves is still within the GW 170817 region, but also it can reach both the PSR J0030 + 0451 error bars. We also consider the case φ(0) = 0.8 in right panels of Fig. 2. Obviously, the mass-radius yields a similar behavior as in the previous case φ(0) = 1, i.e., as we increase the value of parameter α i , an enhancement is observed in the maximum mass. However, by comparing the left and right panels of Fig. 2, one finds that reducing φ(0) leads to increasing the maximum mass. Such a result can be confirmed by comparison of maximum masses for φ(0) = 1 and φ(0) = 0.8 in table II more clearly. It is observed that for potential V 3 (φ) (|α 3 |= 0.1), the maximum mass reaches 2.14M , so in addition to the fact that it still lies within two GW 170817 and NICER regions, it also reaches the PSR J2215 + 5135.
In order to investigate the effect of MPA1 EoS on the neutron star structure, we plot the corresponding mass-radius relation in Fig. 3. By increasing parameter α i or reducing φ(0), the mass-radius relation yields a similar behavior such as previous case (the SLy4 case). Also, similar to the SLy4 case, we observe a greater maximum mass for potential V 3 (φ) comparing to the potentials V 2 (φ) and V 1 (φ) (see table III for more details). It is notable that the curves for all three potentials can support all the observation regions. However, some of the curves for parts (a) and (b) of Fig.  3 is discredited, since it lies out of the LIGO-VIRGO region. As a final result, it is worth mentioning that to support all observational evidence within this EoS, the best value of parameter α 1 is about −0.05 for potential V 1 (φ). Following the previous calculations, we investigate the mass-radius behaviour of neutron stars for the BSK21 EoS in Fig. 4 and table III. It is worth mentioning that similar to the case of MPA1 EoS, there are some curves for parts (a) and (b) of Fig. 4 that are out of LIGO-VIRGO data, and therefore, parameter α i can be considered around −0.05 which is similar to MPA1 EoS. We desire to investigate the mass-radius relation for WFF1 EoS that is reflected in Fig. 5. The general behavior of the mass-radius relation is similar to those of the previous EoSs, i.e, increasing |α i | leads to increasing maximum mass (see table II). It is also observed that for GR limit this EoS could not support the pulsar J2215+5135 and its radius is small (see Fig. 1). So it could not reach NICER observations. Taking Fig. 5 into account, we find that in mimetic gravity for all three potentials, for large α i , the maximum mass is within the pulsar J2215+5135 region. Besides, for potential V 1 (φ) and by considering φ(0) = 1, the mass-radius plot supports the first NICER error bar, III: Mmax and corresponding radius for relativistic EoSs with K = 0.5, φ(0) = 1 (φ(0) = 0.8) and various αi(i = 1.. 3) .       while by changing φ(0) to 0.8, we find that both of NICER error bars are supported by mass-radius diagram. Also, we investigate the DD2 EoS in Fig. 6 in mimetic gravity. It is seen that this EoS can reach the GW190814 event for all potentials. However, for V 1 (φ) model with α 1 = −0.1, we observe that radii are large and the mass-radius diagram cannot support the GW 170817 region. As for the final nucleonic EoS, we consider the FSU2H EoS in Fig. 7. It can be seen that for V 3 (φ), by increasing α, the maximum mass lies in the GW190814 observations while it could not reach this region in GR limit. Furthermore, for all of potentials, the radii of this hadronic EoS increase and they are not in the GW170817 region anymore.
In Fig. 8 and 9, the mass-radius curves are exhibited for hadronic EoSs BHBΛΦ and SFHoY. We calculate maximum mass and radius for different values of α and potentials in table IV. According to Fig. 1, in the GR limit, the maximum mass of star for BHBλΦ could not reach the PSR J2215+5135. However, according to Fig. 8 we observe that the maximum mass is reached in this region for all of V 1 (φ), V 2 (φ) and V 3 (φ) potentials. But the radius for V 1 (φ) is expanded and it is out of GW170817 event. Nevertheless, for the V 2 (φ) and V 3 (φ) potentials, the curves still lie in the GW170817 region.
The second hadronic EoS is SFHoY EoS. In Fig 9, it can be observed this EoS supports the first or even second NICER error bars in mimetic gravity by increasing |α|. But in GR limit, it does not lie in this region. In addition, for all of choice of α parameter and all potentials V 1 (φ), V 2 (φ) and V 3 (φ), the curves support the GW170817 event.
In order to compare the effectiveness of different EoSs, we collect the results of GR and mimetic gravity in table V. In general, we find that mimetic gravity can improve the results and it can support more observational evidence than GR case. Besides, it is evident that MPA1, DD2 and FSU2H EoSs in the context of mimetic gravity can support the mentioned observational data as indicated in table V. According to this interesting result, we focus on these three EoSs to drive other physical properties of neutron stars such as Schwarzschild radius, redshift and compactness.
Schwarzschild radius: considering the obtained function g(r) in Eq. (2.16), and using the horizon radius constraint ( g (r)| r=R = 0), we can extract the Schwarzschild radius (R Sch ). It is clear that by applying the mimetic scalar field (Eq. (2.11)), and using g (r)| r=R = 0, the mentioned scalar potentials become zero on surface of neutron star, i.e., V (φ)| r=R = 0. Therefore, we obtain the Schwarzschild radius for the mimetic gravity (which is the same GR gravity) as As one can see in the equation (3.3) and tables VI-VIII, by increasing maximum mass of neutron star, the Schwarzschild radius increases. Gravitational redshift and compactness: other important quantities in which we intend to extract them, are related to the gravitational redshift and compactness of neutron stars in the mimetic gravity. These quantities may give us information about the strength of gravity of neutron star.
We use from definition of the gravitational redshift and this fact that the mentioned scalar potentials are zero on surface of neutron star, we extract the gravitational redshift as where R is radius of neutron star. The compactness of a spherical compact object may be defined by the ratio of the mass to radius of that compact Now, we are in a position to investigate the gravitational redshift and compactness of the obtained neutron stars of MPA1 EoS in the mimetic gravity. Considering Eqs. (3.4), (3.5), and the obtained results in table III, we find an interesting property of neutron stars. The gravitational redshift and compactness of neutron star in the mimetic gravity are dependent on the scalar potential models. In other words, there are two different manners for these quantities. They decrease by increasing the maximum mass of neutron star in V 1 (φ) model. However, they increase by increasing the maximum mass for V 2 (φ) and V 3 (φ) models (see table VI, for more details). In addition, there are the same behavior for DD2 and FSU2H EoSs (see tables VII-VIII, for more details). Therefore the gravitational redshift and compactness of neutron stars in the mimetic gravity are completely depending on the scalar potential models.

IV. CONCLUSION
The study of neutron stars, one of the possible endpoints of stellar evolution, is an interesting hot topic in physics communities since they provide a natural laboratory for the investigation of high density and pressure medium. The neutron star structure can be considered from theoretical and observational perspectives. In recent years, different properties of diverse neutron stars are accessible from the observational evidence, such as GW170817, PSR J0740+6620, PSR J2215+5135, NICER data on PSR J0030+0451 and GW190814. However, from theoretical point of view, different EoSs in the context of GR and other modified theories of gravitation could not support all observational data, simultaneously. Therefore, to describe the neutron stars, research on suitable EoS and appropriate theory

GR
GW170817 GW190814 PSR J0740+6620 PSR J2215+5135 PSR J0030+04511 PSR J0030+04511 (the first error bar) (the second error bar) SLy4 Mimetic gravity GW170817 GW190814 PSR J0740+6620 PSR J2215+5135 PSR J0030+0451 PSR J0030+0451 (the first error bar) (the second error bar)          of gravitation, consistent with astrophysical observations, is ongoing with great motivation. VIII: Structure properties of neutron star in mimetic gravity for FSU2H EoS with φ(0) = 1 (φ(0) = 0.8) and K = 0.5.  This paper is devoted to study the structure of neutron stars through the use of different state equations in mimetic gravity. Various models of EoS are examined in the context of GR and mimetic gravity. Differentiating between the models of EoS indicated by the measurements of the radius and mass of neutron stars in a theory of gravitation. It was shown that considering a few suitable EoSs, the maximum mass and related radius in mimetic gravity can support the mentioned observational data.
In this paper, we examined eight relativistic/non-relativistic EoSs that seem more realistic and have already received more attention in literature. For the non-relativistic category of EoSs, we considered SLy4, BSK21 and WFF1 models. While for the relativistic set of EoSs, we regarded MPA1, HS (DD2) and FSU2H cases. We also used two hadronic EoSs known as BHBΛΦ and SFHoY as other realistic models. According to the obtained results reflected in different figures and tables, we have found that, in general, the results of mimetic gravity are more consistent with observational data than the GR. Besides, we have shown that the results of mimetic gravity are more or less depending on the functional form of potential and its free parameters. Nonetheless, we have found that the results of full relativistic mean-field-based models of EoSs are in better agreement with observational data than non-relativistic models.
For the three non-relativistic models of EoSs, SLy4, WFF1 and BSK21, we have shown that the mentioned state equations can be improved when we generalize the GR to mimetic gravity. Strictly speaking, regardless of GW190814, other observational evidence can be supported by the results of theoretical counterparts in mimetic gravity for SLy4 and WFF1 models. We should note that the same result can be obtained, equivalently, for BSK21 in GR and mimetic gravity. After that, a relativistic model of EoS can be regarded to examine the possible support of GW190814 data. We have taken the relativistic EoSs MPA1, DD2 and FSU2H into account in GR and mimetic gravity, and found that although such models could not solve our problem in GR, they could remove the inconsistency between the GW190814 data and the theoretical results of previous models in mimetic gravity. In other words, considering these three EoSs are excellent candidates for investigation of neutron star structures in mimetic gravity.
According to the presented tables, one can find that for numerical analysis, we have considered the effects of mimetic gravity generalization as a small changes with respect to GR. In other words, we have investigated small values of potential since we did not have the allowed region of free parameters. So, it is interesting to find permitted region of free parameters of mimetic gravity (its potential) based on observational cosmology or fundamental conceptions of high energy physics, and study the effects of mimetic gravity when it is far from the GR. Besides, taking the modification of the mass-radius relation in mimetic gravity into account, it will be interesting to study the thermal evolution of neutron stars. In other words, it is attractive to investigate the cooling process of isolated neutron stars to compare mimetic gravity and GR. So, one can study the possibility of the direct Urca mechanism for rapid cooling or modern theory of cooling based on the nucleon superfluidity depending on the chosen EoS. Keeping the results of this paper and the mentioned suggestions in mind, we can analyze the behavior of hybrid neutron stars or quark stars as well as other compact objects in mimetic gravity. All these points are under considerations.
As a final point, Astashenok et al. in Ref. [94] have studied the structure of neutron stars in mimetic gravity. They have obtained some interesting results about the structure of massive neutron stars. But, they had not compared their results with observational data. Nonetheless, observational data such as GW170817, PSR J0740+6620, PSR J2215+5135, NICER data on PSR J0030+0451, and GW190814 give us essential information on massive compact objects e.g., neutron stars. In this work, we extended the work of Astashenok et al. by comparing observational data with mimetic theory of gravity. In other words, to obtain a good agreement between theory and observational data, we evaluated the structure of neutron stars in GR and mimetic gravity.