Arxiv:1912.00685V2 [Gr-Qc] 30 Mar 2020

USTC-ICTS-19-29 Relativistic stars in mass-varying massive gravity Xue Sun1 and Shuang-Yong Zhou2 1School of Physical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, China 2Interdisciplinary Center for Theoretical Study, University of Science and Technology of China, Hefei, Anhui 230026, China and Peng Huanwu Center for Fundamental Theory, Hefei, Anhui 230026, China (Dated: March 31, 2020) Mass-varying massive gravity allows the graviton mass to vary according to different environments. We investigate neutron star and white dwarf solutions in this theory and find that the graviton mass can become very large near the compact stars and settle down quickly to small cosmological values away the stars, similar to that of black holes in the theory. It is found that there exists a tower of compact star solutions where the graviton mass decreases radially to zero non-trivially. We compute the massive graviton effects on the mass-radius relations of the compact stars, and also compare the relative strengths between neutron stars and white dwarfs in constraining the parameter space of mass-varying massive gravity. I. INTRODUCTION ton mass near the horizon, there is an extra potential barrier to the left of the photosphere barrier in the mod- ified Regge-Wheeler-Zerilli equation, and this can lead to Recent advances in astrophysics, particularly the ar- gravitational wave echoes [18, 19] in the late time ring- rival of gravitational wave astronomy [1, 2], have pro- down waveform when the black hole is perturbed [16]. vided fresh new opportunities to test gravity in the strong field regime with compact astronomical objects. Whilst The detection of gravitational waves from neutron star in Einstein gravity the graviton is massless, an interest- mergers [2] has added a new dimension to our ability to ing class of alternative theories of gravity is to let the unravel the nature of strong gravity and the structure of graviton become massive [3{6]. See [7{9] for a recent relativistic compact stars. Neutron star solutions have review on recent developments in constructing massive been previously investigated in the original dRGT model gravity and bi-gravity models [3{6] that are free of the [20] and in bi-gravity [21] with a graviton mass that is Boulware-Deser ghost [10] and their applications. Typi- much greater than typical cosmological scales. Relativis- cally, the graviton mass is set to be close to the current tic star solutions in dRGT-like models with a singular Hubble scale [11], as one of the main phenomenological reference metric have also been studied [22, 23]. applications of these models is to explain the late time In this paper, we will further investigate the compact cosmic acceleration [12, 13]. For such a small graviton star solutions in mass-varying massive gravity. We will mass, or such a long Compton wavelength, simple di- compute the neutron star and white dwarf solutions, us- mension analysis suggests that these models can be best ing the APR [24] and SLy [25] equations of state for neu- constrained by large scale or cosmological observations, tron stars and the Chandrasekhar equation of state for although sometimes the solar system gravity tests may white dwarfs, and analyze the behavior of the graviton achieve comparable bounds, thanks to their superb accu- mass near compact stars. We will study how the effects of racy [11]. a varying graviton mass affect the mass-radius diagram However, the graviton mass does not have to be uni- of compact stars and compare the relative strengths of formly small across the whole spacetime. In mass- neutron stars and white dwarfs in constraining the pa- varying massive gravity [14, 15], the original de Rham- rameter space of mass-varying massive gravity. Gabadadze-Tolley (dRGT) model [4] is augmented with The paper is organized as follows: We introduce mass- an extra environmental scalar field and the graviton mass varying massive gravity in Section II, and in Section III, arXiv:1912.00685v2 [gr-qc] 30 Mar 2020 now depends on the environmental field, which can take we reduce the field equations into a system of ordinary different values in various different astrophysical environ- differential equations for compact stars and describe our ments. In [16], it is shown that the graviton mass can be- numerical setup; In Section IV, we set up the equations come extremely large near the event horizon of the hairy of state for neutron stars, compute their solutions and black hole in mass-varying massive gravity [17], while set- investigate how the mass-radius diagrams change with ting down to the small cosmological value away from the the model parameters; In Section V, we investigate the horizon, thus satisfying all the current tests of gravity and white dwarf solutions and compare the relative strengths yet still giving rise to interesting deviations from general of neutron stars and white dwarfs in constraining the relativity that are testable in the current and upcoming parameters in mass-varying massive gravity; We conclude experiments. An important feature of those hairy black in Section VII. We shall use geometric units which set holes is that because of the rapid increase of the gravi- c = G = 1 and express dimensionful quantities in terms 2 of kilometers. α and β are related to α3 and α4 by α = 1 + α3; β = α3 + α4: (8) II. MASS-VARYING MASSIVE GRAVITY Mass-varying massive gravity (MVMG) [14] is ex- III. SETUP tended dRGT massive gravity [4] where the mass of the graviton can vary according to the value of a scalar field We consider a static and spherically symmetric Ansatz in different environments. Its action is given by for the scalar field σ = σ(r) and for the metric 1 Z p R ds2 = −a(r)dt2 + 2b(r)drdt + e(r)dr2 + d(r)dΩ2; (9) S = d4x −g + V (σ)U(K) 8π 2 2 2 2 2 2 dsη = −dt + dr + r dΩ : (10) 1 − gµν @ σ@ σ − W (σ) + S [g ] ; (1) 2 µ ν m µν Plugging these into the Einstein tensor, the nonzero components are where Sm is the action of the conventional matter which only couples to the dynamic metric g and the dRGT 2a0cdd0 + 4acdd00 − 4c2d − 2ac0dd0 − acd02 µν Gt = ; (11) graviton potential is given by t 4c2d2 0 0 02 00 µ ν µ ν ρ µ ν ρ σ t bc dd + bcd − 2bcdd U(K) = K K + α3K K K + α4K K K K ; (2) G = ; (12) [µ ν] [µ ν ρ] [µ ν ρ σ] r 2c2d2 0 0 02 µ µ p −1 µ −1 µν −4cd + 2a dd + ad with Kν ≡ δν − g η , g = (g ) being the in- r ν Gr = 2 ; (13) verse metric and η = (ηµν ) being the Minkowski metric. 4cd 1 The anti-symmetrization is defined with unit weight, e.g., Gθ = − a0c0d2 + 2ca0dd0 − ac0dd0 − acd02 µ ν µ ν µ ν θ 4c2d2 K[µ K ν] = (KµKν − Kν Kµ)=2. There are of course many 00 2 00 choices for V (σ) and W (σ), which are functions of the + 2a cd + 2acdd ; (14) environment scalar field σ. We consider a simple model 0 where these potentials are [16] where denotes a derivative with respective to the radius φ θ r and by symmetry we have Gφ = Gθ. The matter energy 1 V (σ) = m2 + m2σ4;W (σ) = m2 σ2 + λ σ6; (3) momentum tensor is taken to be a perfect fluid 0 2 σ σ Tµν = (ρ(r) + p(r))uµuν + p(r)gµν ; (15) where m0 is chosen to be at most of the order of the Hubble scale, potentially generating the late time cosmic where ρ is the matter energy density, p is the matter pres- acceleration. For star solutions, which are well within sure and uµ is the 4-velocity of the perfect fluid, which, in 2 the Hubble horizon, we can safely neglect the m0 term. the static case, is determined completely by the normal- 6 The σ term, where λσ is small, is added for the stability µ ν ization gµν u u = −1. Then the non-trivial components of the vacuum and can be neglected for local star solu- of the matter energy momentum tensor are tions [16]. So apart from α3 and α4, the essential theory parameters are m and m , the dependence of which will (p + ρ)b σ T t = −ρ, T t = ;T r = T θ = T φ = p: (16) be explored for the compact stars in the next sections. t r a r θ φ The equations of motion for this model are given by On the other hand, the nonzero components of the (σ) Gµν = 8πTµν + Tµν + V (σ)Xµν ; (4) energy-momentum tensor generated by the scalar σ are p µv p @µ −gg @vσ = −g (Wσ − VσU) ; (5) 02 (σ)t (σ)θ (σ)φ aσ T t = T θ = T φ = −W (σ) − ; (17) where Vσ = dV=dσ, Wσ = dW=dσ, Tµν is the energy 2c momentum tensor from S , m and the nonzero components of the effective energy- 1 momentum tensor from the massive graviton are T (σ) = @ σ@ σ − g gργ @ σ@ σ + W (σ) (6) µν µ ν µν 2 ρ γ b2n2k − k K Xt = 1 2 3 + k (αk + 2) ; (18) and t b2n2 − 1 3 3 bK X = − g Kρ − Kµ1 g t r 3 µν µρ ν µ1 µν Xr = −Xt = ; (19) c ρ σ µ1 ρ µ1 µ2 2 2 + α gµρKσKν − Kµ gµρKν + K K gµν b n k − k K 1 [µ1 µ2] Xr = 2 1 3 + k (αk + 2) ; (20) r 2 2 3 3 ρ σ τ µ1 ρ σ b n − 1 − β gµρKσKτ Kν − Kµ gµρKσKν 1 Xθ = Xφ = α (k k + k k + k k ) + βk k k θ φ 1 2 1 3 2 3 1 2 3 µ1 µ2 ρ µ1 µ2 µ3 + K K gµρK − K K K gµν : (7) [µ1 µ2] ν [µ1 µ2 µ3] + k1 + k2 + k3; (21) 3 2 where K3 = βk3 + 2αk3 + 1 and k1, k2 and k3 are the µ 1.0 eigenvalues of matrix Kν , explicitly given by 0.9 s k3+ 2 0.8 k1 = 1 − ; (22) p 2 k3- a + e + (a + e) − 4c 0.7 s 2 A 0.6 k2 = 1 − ; (23) p 2 a + e − (a + e) − 4c 0.5 r r2 0.4 k3 = 1 − ; (24) d 0.3 with n = 2=(p(a + e)2 − 4c + a − e) and c = b2 + ae.

Load more