The Palatini star: exact solutions of the modified Lane–Emden equation

Two exact solutions for n=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=0$$\end{document} and n=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=1$$\end{document} of the Palatini-modified Lane–Emden equation are found. We have employed these solutions to describe a Palatini–Newtonian neutron or low mass brown dwarf star and compared the result with the pure Newtonian counterpart. It turned out that for the negative parameter of the Starobinsky model the star is heavier and larger.


Introduction
Although there is no doubt of the beauty and the overall correctness of the General Relativity [1][2][3], the theory suffers from a number of shortcomings [4][5][6][7][8] which could be addressed by extensions, provided by the Extended Theories of Gravity (ETGs), of the Einstein's proposal. Apart from the cosmological arguments on the need of looking for modified theories of gravity [9][10][11][12][13][14][15][16], there are also indications coming from astrophysics strengthening reasons for searching for the alternatives. In the light of the very recent observation of the pulsar PSR J2215+5135 with the mass M 1 = 2.27 +0. 17 −0.15 M for the neutron star [17] putting in trouble exotic forms of matter (hyperons or deconfined quarks) and the previous discoveries of the massive neutron stars as 2.01±0.04M (PSR J0348+0432) [18], and 1.97 ± 0.04M (PSR J1614-2230) [19], the ETGs [20,21] allow to exceed the maximal mass of 2M demonstrated by GR without introducing exotic components.
Let us just briefly recall the Palatini formalism which assumes that the spacetime geometry is described by two independent structures: a class of Lorentzian metrics which are conformally related to each other, and a connection which later on, as a dynamical feature of the action, will turn out to be a Levi-Civita connection with respect to a metric conformally related to the metric g. Thus, we consider a standard action of the f (R) gravity of the form whereR =R μν g μν is the Palatini-Ricci scalar constructed using two independent objects, the connection˜ and the metric g. The variation with respect to the metric g and the independent connection˜ is given, respectively, by the formulas where T μν is the energy momentum tensor and f (R) = d f/dR. Equations (3) imply that the connection∇ β is the Levi-Civita connection of the conformal metric h (up to a projection, for the detailed discussion see e.g., [62]): Let us observe that Palatini gravity is equivalent to GR when one considers the linear function f of the form f (R) = R − 2 , and the independent connection turns out to be Levi-Civita connection of the metric g.
We are going to use here the following interesting fact: the equations in question can be rewritten [25] in the terms of the conformal metric h μν [63,64] and the scalar field defined as = f (R): where the bar quantities are defined as follows:Ū ( ) = (R − f (R))/ 2 , and the appropriate energy momentum tensor readsT μν = −1 T μν . One should also bear in mind thatR μν =R μν ,R = h μνR μν = −1R and h μνR = g μνR . The above equations are field equations of the scalartensor action for the metric h μν and (non-dynamical) scalar field .
We will now focus on the spherical-symmetric spacetime whose g-metric is given by ds 2 = −B(r )dt 2 + A(r )dr 2 + r 2 dθ 2 + r 2 sin 2 θ dϕ 2 (6) while the conformal metric is given by (4). Solving the equations in the Einstein frame, using the properties and the interpretation of the Palatini gravity it was shown [54,65] that the generalized Tolman-Oppenheimer-Volkoff (TOV) equations read where the generalized energy density Q and the pressure are defined as Recall thatŪ and depend on the choice of the model we are interested in. Let us notice that the r -coordinate as well asÃ are objects from the Einstein frame thus one applies the conformal transformation: In what follows we will focus on the Starobinsky model f (R) =R + βR 2 [11] considered in [66] in the context of the Lane-Emden equation. The generalized versions of Lane-Emden equation were already studied in [66][67][68][69][70][71]74,75], and some of them were obtained from the generalized TOV equations [65,[76][77][78][79][80][81][82][83].
We follow the Weinberg convention [84], and hence the signature of the metric is Given the prominent role played by exact solutions and exactly solvable models in modern physics, see e.g., [85,86] and references therein, it is natural to look for exact solutions of the generalized Lane-Emden equation introduced in [66] in the context of Palatini gravity. In the next section we present such solutions for two cases: that of an incompressible star (n = 0), when the solution coincides with a well-known one, and that of a neutron star (n = 1), when the solution we found is, to the best of our knowledge, new. In Sect. 3 the physical implications of this new solution are explored using both analytical and numerical techniques, and Sect. 4 contains discussion.

Exact solutions of the Lane-Emden equation
The standard Lane-Emden equation coming from the Newtonian limit of the GR equilibrium equations (TOV equations) and polytropic equation of state p = Kρ γ , where K and γ are polytropic parameters, has the following form: here n = 1 γ −1 . Recall the relation among dimensionless quantities appearing in (11) and the physical ones: where p c and ρ c are central pressure and energy density, respectively. Equation (11) is known to possess three exact solutions for n = {0, 1, 5} and, generally speaking, it proved pretty much impossible to find exact solutions of the generalized Lane-Emden equations considered in [67][68][69][70][71]74,75]. Thus, for other values of n or in the case of generalized Lane-Emden equations, the approximate and/or numerical methods were applied in order to examine the properties of stars.
As shown in [66], the generalized Lane-Emden equation for the Starobinsky Lagrangian f (R) =R + βR 2 in the Einstein frame is given by 1 ξ where α = κc 2 βρ c , with ρ c being a central density. The extra term including α occurring in (14) comes from the potential U appearing in (5a) and definitions (9) whose form depends on the Lagrangian. The details of the derivation can be found in [66]. Under the conformal transformationξ 2 = ξ 2 , where = 1 + 2αθ n , the above equation takes the form An approximation of (15) was analyzed from the perspective of finding the numerical solutions in [66]. Notice that for α = 0 we have = 1 and (15) boils down to (11). Now note that for n = 0 equation (14) becomes linear, and its general solution is readily found to read where C i are arbitrary constants. The transformationξ 2 = ξ 2 in this case amounts to a rescalingξ = (1 + 2α) 1/2 ξ , so the general solution of (15) for n = 0 reads where we have rescaled C 1 for the sake of convenience. If C 1 = 0 then this solution has singularity at ξ = 0, which is unphysical, so we should assume that C 1 = 0. Upon further imposing the boundary conditions in the Jordan frame, θ(0) = 1, θ (0) = 0, we find that (18) becomes Notice that it has the same form as the solution of the standard Lane-Emden equation (11) for n = 0 describing incompressible stars, and for this reason it has been discussed in many textbooks (see e.g., [84]).
Our key observation now is that for n = 1 the generalized Lane-Emden equation (15) has an exact solution which appears to be new and can be employed to describe a Newtonian neutron star.

Newtonian neutron star in Palatini gravity
As we have already mentioned, the solutions of the Lane-Emden equation have to satisfy the boundary conditions, namely θ(0) = 1 and θ (0) = 0. The solution (20) for n = 1 turns out to satisfy them only if α = 3 4 . This value lies in the allowed range α > −1/2. Indeed, because of the boundary condition θ(0) = 1 the conformal factor for ξ = 0 takes the value 1 + 2α. Now, for α = −1/2 this value becomes zero which is not allowed since the conformal factor may not change sign, cf.
[66] for details. Moreover, the numerical analysis in [66] also indicates that the values of α smaller than −1/2 are unphysical.
This means that when we assume that the value of the central density of an average neutron star is ρ c ∼ 8×10 17 kg m 3 , the Starobinsky parameter is found to be β ∼ −5.02978 × 10 6 m 2 .
Let us recall that α = κc 2 βρ c and because of the Weinberg's convention (κ = −8π G/c 4 ) that we follow throughout the paper, κ introduces a minus to the definition of the parameter α. Therefore, for positive α we deal with a negative value of the Starobinsky parameter β. In order for the solution (20) to describe a polytropic (with n = 1) star, one obtains α = 3/4 which immediately gives values of the parameter β, which however will depend on a type of stars. Thus, the denser star, the larger absolute value of the Starobinsky parameter. This example also demonstrates how the theories of Palatini family are sensitive to the local energy density distribution [60]. Now we are able to find the star's mass, radius, and central density. Recall that by introducing the dimensionless quantities [84,87]: we can write down the star's mass, the mass-radius relation, central density, and temperature, in the form where the last equation for the temperature profile is given under the assumption that the gas is ideal with the equation of state T = k B ρT μ , with k B and μ being the Boltzmann constant and mean molecular weight, respectively.
However, let us notice that the above quantities were defined for the General Relativity setting. Since we wish to apply them to the Palatini gravity, we must keep in mind that the above expressions are the ones in Einstein frame. Thus, the quantities ω n and δ n should be rewritten, see Eqs. (31) and (32) below; notice that the mass (8) is still written in the Einstein frame because the r -coordinate comes from the conformal metric h. This can be directly verified by applying the generalized Lane-Emden equation (15) to (8) while taking into account that the conformal transformation in the case of Starobinsky model preserves the polytropic equation of state for small values of p [61]. Thus, in the case of Palatini gravity we should have written Then, applying the conformal transformation relationξ 2 = ξ 2 , we find that Now we are able to calculate the physical values of the star's mass. Moreover, since the solution for n = 1 in the case of the standard Lane-Emden equation is known, θ N = sin ξ ξ , we can compare Newtonian neutron stars in the two models. Here and below we will denote the solution and values obtained from the standard Lane-Emden equation by the subscript N .
We plot the solutions in Fig. 1. Notice that Fig. 1 also represents the dimensionless temperature profile as T ∼ θ .
The neutron star's radius ξ R is defined by the first zero of the solution, θ(ξ R ) = 0, and hence we find that ξ R = √ 15 while ω 1 ∼ 3.87, and δ 1 = 5. Therefore, we can immediately compare the masses and radii of the stars where we have used (ξ R ) N = π and ω N = π [84]. Let us notice a very important result here: Since one describes the low mass brown dwarfs, that is, the stars with the masses M 4 × 10 −3 M , by the polytropic models 2α(10+ξ 2 ) , with α = 3/4. The solution for n = 0, which is the same for both cases, is given by the dotted curve with n = 1, we may easily compare the result coming from Palatini gravity with the one from GR. From (26) we see that the radius for n = 1 is independent of the star's mass whose value in the case of GR is ≈ 0.1R being close to the observed value [71]. However, in the modified gravity, because the polytropic constant K is independent of the theory of gravity one deals with and thus immediately we have in our case that the radius of a low mass dwarf star is ≈ 0.123R . This can be compared to the measurements of the mass-radius relation for such objects as soon as for example GAIA survays will release the data [72,73].

Conclusions
We have found two exact solutions for n = 0 and n = 1 of the modified Lane-Emden equation (15) using the Einstein-Jordan frame correspondence [88][89][90][91], keeping in mind that physical properties of the star such as mass and radius are the ones appearing in Jordan frame, which we demonstrated clearly in the discussion preceding the formulas (29) and (30). We have not examined the result for n = 0 since it would describe an incompressible star. We should notice that the latter solution has exactly the same form as the solution of the standard Lane-Emden equation. However, the case of the Newtonian neutron star (the solution for n = 1), although being just a toy model, shows that one deals with the stellar objects bigger than the stars of the pure GR models. Let us also stress that in this work we have used the exact formula for obtaining the star's mass, that is, Eqs. (25) and (31), in contrast to the previous work [66] where an approximate formula for mass was used. Moreover, with the result on the radius of a low mass dwarf star which is also modelled by the polytropic EoS with n = 1 we have shown that the theory can be probed thanks to the observation of such objects. The modification introduced by Palatini gravity changes the value of their radii with respect to the GR one and hence the theories could be distinguished from one another. We expect that the examination of the full relativistic theory, described by the modified TOV Eqs. (7) and (8), will also provide a similar conclusion, that is, that the Palatini star is larger and heavier for the negative Starobinsky parameter β (positive α parameter), without a need of introducing exotic matter in their internal structure. It should be clarified here that the result was obtained for the negative Starobinsky parameter β: the parameter α appearing in the modified Lane-Emden equation consists of the constant κ, which had been defined to be negative.
Moreover, requiring that the solution (20) has to satisfy the boundary conditions accompanying the Lane-Emden equation, we obtained the exact value of the parameter α = 3 4 which led to the fact that the Starobinsky parameter is negative and of the order 10 6 m 2 under the assumption of the central energy density ρ c ∼ 8 × 10 17 kg m 3 . Let us also mention a discussion in [90] that in the Palatini theories (for example, f (R) gravity and Born-Infeld inspired gravity, cf. e.g., [98][99][100][101]) the energy density must be continuous and differentiable function in order to avoid divergences in the field equations. This comes from the fact that one obtainsR =R(ρ) which appears in the conformal function taking part in the conformal transformation which includes the derivatives of ρ.
Therefore, we would like to conclude that the Palatini gravity is a viable alternative to General Relativity since it passes Solar System tests (vacuum solutions are equivalent to General Relativity with cosmological constant) [92][93][94], introduces inflation preserving late time accelerated expansion [63,64,96,97], gives clues on the Dark Matter problem [102], and provides conditions for existence of stable relativistic stars [54] that are similar to General Relativity.
Acknowledgements AW acknowledges financial support from FAPES (Brazil). The research of AS was supported in part by the Grant Agency of the Czech Republic (GAČR) under grant P201/12/G028, under RVO funding for IČ47813059, and the INTER-EXCELLENCE project No. LTI17018 that supports the collaboration between the Silesian University in Opava and the Astronomical Institute in Prague. AS is pleased to thank the Institute for Theoretical Physics of the Wrocław University and, in particular, Prof. Andrzej Borowiec, for the warm hospitality extended to him in the course of his visits to Wrocław where the present research project was initiated. Furthermore, the authors warmly thank Andrzej Borowiec, Gonzalo Olmo, Diego Rubiera-Garcia, and Her-mano Velten for stimulating discussions and comments. We also thank the anonymous referee for useful suggestions.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This paper is of theoretical nature; no experimental datasets were used.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .