Interior solutions of relativistic stars in the scale-dependent scenario

We study relativistic stars in the scale-dependent scenario, which is one of the approaches to quantum gravity, and where Newton's constant is promoted to a scale-dependent quantity. First, the generalized structure equations are derived here for the first time. Then they are integrated numerically assuming a linear equation-of-state in the simplest MIT bag model for quark matter. We compute the radius, the mass and the compactness of strange quarks stars, and we show that the energy conditions are fulfilled.


Introduction
Einstein's General Relativity (GR) [1] is widely accepted as a relativistic theory of gravitation, which is at the same time both beautiful and very successful [2,3]. Remarkably, a series of predictions have been confirmed observationally over the past 100 years or so. Apart from the classical tests and solar system tests [4], over the last four years or so two more predictions of GR, namely the existence of gravitational waves as well as the existence of black holes, have been confirmed as well thanks to the efforts of the LIGO/VIRGO collaborations [5].
Despite its success, it has been known for a long time now that GR is not compatible with quantum physics. It is well-known that the formulation of a consistent quantum theory of gravity is still an open task in modern theoretical physics. Although as of today several approaches to the problem do exist in the literature (for a partial list see e.g. [6][7][8][9][10][11][12][13][14] and references therein), there is one property in particular that all of those have in common. Namely, the basic quantities that enter into the action defining the model at hand, such as Newton's constant, electromagnetic coupling, the cosmological constant etc, become scale dependent (SD) quantities. This does not come as a surprise of course, since scale dependence at the level of the effective action is a generic feature of ordinary quantum field theory.
In black hole physics the impact of the SD scenario on properties of black holes has been studied over the a E-mail: grigorios.panotopoulos@tecnico.ulisboa.pt b E-mail: angel.rincon@pucv.cl c E-mail: ilidio.lopes@tecnico.ulisboa.pt last years, and it has been found that the scale dependence modifies the horizon, the thermodynamics as well as the quasinormal spectra of classical black hole backgrounds [15][16][17][18][19][20][21][22]. However, the astrophysical implications of the SD scenario should be investigated as well. In the present work we propose to obtain for the first time interior solutions of relativistic stars in the SD scenario, see however [23] for an alternative approach. In particular, here we shall focus on strange quark stars, which may be viewed as an alternative to the more conventional paradigm based on neutron stars. Although as of today it is a speculative class of compact objects, strange quarks stars cannot conclusively be ruled out yet. As a matter of fact, there are some claims in the literature that there are currently some observed compact objects with peculiar features (such as small radii for instance) that cannot be explained assuming the known hadronic equations-ofstate for neutron stars, see e.g. [24][25][26], and also Table 5 of [27] and references therein.
Our work is organized as follows: After this introduction, in the next section we briefly review the formalism of the SD scenario. In the third section we apply it to the case of relativistic stars, and we present and discuss our main numerical results. In the last section we close our work with some concluding remarks. We adopt the mostly positive metric signature, (−, +, +, +), and we work in geometrical units where the speed of light in vacuum as well as the classical Newton's constant are set to unity, c = 1 = G N .

Scale-dependent gravity
This section is devoted to summarize the main features about the formalism used along this manuscript, i.e., scaledependent gravity. The idea is largely inspired by the asymptotic safety program and related approaches as the well-known Renormalization group improvement method [28][29][30][31]. Taking them as inspiration, some authors introduced the now known scale-dependent gravity which has been systematically used in black holes physics [32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47], and recently, in cosmological problems [48,49]. The crucial point of scale-dependent gravity is to promote the classical parameter {A 0 , B 0 , (· · · ) 0 } to functions which depends on the energy scale k, namely {A k , B k , (· · · ) k }. Depending of the case, the arbitrary scale k is linked to the radial coordinate (assuming circular symmetry) or it is connected to the physical time. In what follows, we will consider three different contributions accounted into the following action: where the first term, S EH , is given by the well-known Einstein Hilbert action, the second term, S M , is associated with a perfect fluid distribution, and the third term, S SD , encodes the scale-dependent sector. In light of this, the scale-dependent coupling of the theory is the Newton's coupling G k only (which can be related with the gravitational coupling by κ k ≡ 8πG k ). In addition, there are two independent fields, which are the metric tensor g µν (x) and the scale field k(x). The equations of motion obtained from a variation of (1) with respect to g µν (x) to obtain: where the effective energy momentum tensor, T effec µν , is defined in such a way that it include both the usual matter fields T µν and the contribution of the G-varying part ∆t µν . Then, such a tensor is defined to be Also, notice that the two tensor defined in the right hand side are given by the expressions: and At this level, some comments are in order. Firstly notice that, in a quantum field theory, the corresponding renormalization scale k needs to be linked with the parameters of the physical system under study. For background solutions of the gap equations, it is not constant any more. Therefore, an arbitrarily non-constant k = k(x) implies that the set of equations of motion does not close consistently. The above means that the stress energy tensor is most likely not conserved for almost any choice of the functional dependence k = k(x). This pathology has been considerably investigated in the context of renormalization group improvement of black holes in asymptotic safety scenarios. To clarify the situation we should remark that the problem with the conservation law appears due one consistency equation is missing. The latter equation can be computed taking the variation of the corresponding action respect to the arbitrary scalar field k(r), namely which is usually interpreted as a variational scale setting procedure [31,50]. To guarantee the adequate conservation of the stress energy tensor, we should to combine (6) and the equations of motion. This, route, however, introduce a unavoidable problem, i.e., it is required first to compute the β-functions of the problem. It is well-known that the βfunctions are not unique, reason why use them introduce some uncertainties. A way to circumvent the aforementioned issue is to impose some supplementary condition.
In exterior solutions, the null energy condition (one of the four energy conditions demanded in general relativity) has been typically used as a supplementary equation. Thus, we promote the classical couplings to radial-dependent couplings (given that O(k(r)) → O(r)) and, supported by the NEC, we are able to solve the functions involved. Thus, this philosophy of assuring the consistency of the equations by imposing a null energy condition will also be applied for the first time in the following study on interior solutions of relativistic stars.

Hydrostatic equilibrium of relativistic stars
Here we briefly review relativistic stars in GR, and then we generalize the structure equations in the SD scenario.

Structure equations
For non-rotating objects we seek spherically symmetric solutions in Schwarzschild coordinates, (t, r, θ, φ), making as usual for the metric tensor the following ansatz where dΩ 2 = dθ 2 + sin 2 θ dφ 2 is the usual line element of the unit two-dimensional sphere. In addition, we introduce for convenience the mass function m(r), which is defined to be Assuming for matter a fluid characterized by a stressenergy tensor of the form with ρ being the energy density and p being the pressure, the tt and rr field equations yield respectively, where a prime denotes differentiation with respect to the radial coordinate r. Instead of the angular field equations one may employ the conservation of energy Therefore, the Tolman-Oppenheimer-Volkoff equations [51,52] read Therefore, the solution for ν(r) is given by where In the following step, we now wish to generalize the standard structure equations valid in GR in the SD approach to quantum gravity. To that end, we first recall that in this framework Newton's constant is promoted to a function of the radial coordinate, G(r), and that in the field equations the total (or effective) stress-energy tensor has two contributions, namely one from the ordinary matter, T µν , and another from the contribution of the G-varying part, ∆t µν where the G-varying part is given by [54] Therefore, it is straightforward to obtain the new structure equations, which are found to be (G(r)m(r)) = 4πG(r)r 2 ρ ef f (r) (22) ν (r) = G(r) m(r) + 4πr 3 p ef f (r) r 2 (1 − 2G(r)m(r)/r) and we spare the details for the last equation, as it is too long to be included here. We have checked that the new generalized structure equations are reduced to the usual TOV equations when Newton's constant is taken to be a constant, G (r) = 0 = G (r). Finally, there is an additional differential equation of second order for G(r), which is the following [54] 2 and which must be supplemented by two initial conditions at the centre of the star, G(0) = G c and G (0) = G 1 . In this framework the metric potential A(r) and the mass function m(r) are related via It should be pointed out that the exterior solution remains the same in the SD scenario assuming a vanishing cosmological constant. Therefore, the interior solution obtained assuming a r-varying gravitational constant will be matched to the usual Schwarzschild vacuum solution at the surface of the stars.

Equation-of-state
In order to close the system of equations, an equation-ofstate for quark matter must be assumed. Here we shall adopt the very popular and widely used in the literature "radiation plus constant" analytic function which is the simplest version of the MIT bag model [55,56]. The model is characterized by three parameters, namely the mass of the s quark, m s , the QCD coupling constant, α c , and the bag constant, B. In the simplest MIT bag model m s = 0 = α c , while B = 60 M eV /f m 3 [57,58]. More realistic and sophisticated EoSs have been developed and considered over the years in the literature. For instance, at asymptotically large densities color superconductivity effects [59,60] become important. Quark matter is in the color flavor locked (CFL) state [61,62], in which quarks form Cooper pairs of different color and flavor, and where all quarks have the same Fermi momentum and electrons cannot be present. Other possibilities include models that incorporate a chiral symmetry breaking [63,64], or models based on perturbative QCD studies [65,66], and others [67,68].

Numerical results
Our numerical results for the quantities of interest, namely p(r), ρ(r), m(r), ν(r), G(r), are summarized in Figures 1  and 2, and also in Table 1 where we show the mass, the  radius and the compactness of the stars as well as the initial conditions for G(r).
We have obtained five interior solutions assuming different central values for the gravitational coupling, and imposing different central values p c and ν c such that the following conditions are satisfied The first case (dashed black line) corresponds to GR, which is also shown for comparison reasons. In particular, the five solutions obtained in this work and shown in Figures 1 and 2 Fig. 1 (left panel) shows the scale-dependent gravitational coupling as a function of radial coordinate. Both positive and negative values of G(r) lead to realistic interior solutions. The horizontal line (dashed black line) corresponds to the case where G(r) = constant = 1, which is the GR case. All our numerical solutions tend to G(r = R) = 1 at the surface of the star, which was imposed right from the start. As it is expected, the deviation from GR is indeed quite small. Fig. 1 (middle panel) depicts the metric potential e 2ν(r) as a function of the radial coordinate. We notice that the inclusion of a r-varying gravitational coupling shifts the potential to higher values in comparison with the classical counterpart (dashed black line). Finally, in Fig. 1 (right)

m(r) in solar masses versus normalized radial coordinate.
Our results show that the mass function of the compact object is lower than in GR, irrespectively of the sign of the gravitational coupling. Fig. 2 (left) shows the energy density versus normalized radial coordinate for the five cases considered here. Once again, we observe that the scale-dependence inclusion induces a non-trivial deviation to the density, which is a physical quantity of the model. Similarly, the pressure is a monotonically decreasing function, lower than the GR prediction. In Table 1 we show the mass M , the radius R, the compactness C ≡ M/R as well as the initial conditions for the gravitational coupling, G c , G 1 . Clearly, the r-varying nature of G(r) leads to smaller and lighter strange quark stars, which are found to be less compact, irrespectively of the sign of G 1 .
Next, and before we finish, a check regarding the energy conditions should be made. The obtained solutions must be able to describe realistic astrophysical configurations. Therefore, we should investigate if the energy conditions are fulfilled or not. We require that [69][70][71] ρ ≥ 0 (28) We plot the (normalized) quantities ρ/B, p/B, 3p/B in Fig. 3 for GR (in black) and for one of the other four solutions (in red). We observe i) that they are positive and ii) that the energy density remains larger that both p and 3p throughout the star. Since all three quantities are positive, clearly the first, the second and the last energy conditions are fulfilled. Moreover, since ρ always remains higher than p and 3p, the other two energy conditions are fulfilled as well. We have obtained qualitatively similar behaviour for the rest of the solutions. We thus conclude that the interior solutions obtained in the present work are realistic solutions within the framework adopted here, and as such they are able to describe realistic astrophysical configurations.
As a final remark it should be stated here that in the present article we took a first step towards the investigation of spherically symmetric strange quark stars in the scale-dependent scenario assuming the simplest MIT bag model equation-of-state. Clearly, there is still a lot of work to be done. For instance, one may study i) more sophisticated equations-of-state for quark matter, ii) rotating stars, or iii) other compact objects, such as neutron stars or white dwarfs. We hope to be able to address some of those interesting issues in future works.

Conclusions
Summarizing our work, in the present article we have obtained interior solutions of relativistic stars in the scaledependence scenario. In particular, we have studied strange quark stars assuming for quark matter the simplest MIT bag model equation-of-state. First we presented the new structure equations describing the hydrostatic equilibrium of the stars. The new equations generalize the usual TOV equations of GR, and they boil down to those when Newton's constant is taken to be a constant, G (r) = 0 = G (r). Then we numerically integrated the structure equations, and we computed the radius and the mass of the stars for a varying Newton's constant, both increasing and decreasing throughout the objects. We also have shown that the energy conditions are fulfilled. Our numerical results indicate that the r-varying nature of the gravitational coupling in the SD scenario leads to less compact strange quark stars, which are found to be both lighter and smaller, irrespectively of the sign of G 1 .