Abstract
We analyze the stability of self-gravitating systems which dynamics is investigated using the collisionless Boltzmann equation, and the modified Poisson equation of Eddington-inspired Born–Infield gravity. These equations provide a description of the Jeans paradigm used to determine the critical scale above which such systems collapse. At equilibrium, the systems are described using the time-independent Maxwell–Boltzmann distribution function \(f_0(v)\). Considering small perturbations to this equilibrium state, we obtain a modified dispersion relation, and we find a new characteristic scale length. Our results indicate that the dynamics of self-gravitating astrophysical systems can be fully addressed in the Eddington-inspired Born–Infeld gravity. The latter modifies the Jeans instability in high densities environments, while its effects become negligible in star formation regions.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In General Relativity (GR), matter is minimally coupled with the metric and the Einstein–Hilbert Lagrangian, which is linear in the Ricci scalar, gives rise to second order field equations. These are able to explain the dynamics of the particles up to a solar system scale, but they fail at scales of galaxies and beyond. The dynamics of self-gravitating systems and the current period of accelerated expansion of the Universe cannot be explained by just baryonic matter. Thus, GR needs to incorporate two unknown components to explain the dynamics at both galactic/extragalactic and cosmological scales. Specifically, almost \(\sim \)68% of the total amount of the matter and energy in the Universe should be in form of the cosmological constant, or more in general of dark energy, while \(\sim \)26% should be in the form of invisible and exotic particles, named dark matter. Nevertheless their fundamental nature is still unknown [1,2,3,4]. The need to incorporate them has been interpreted as a breakdown of GR at astrophysical and cosmological scales, opening the door to alternative theories of gravity.
Generalizations of the gravitational action have been extensively explored to overcome the need of these two exotic components. On th one hand these are motivated by their capability to explain the dynamics of self-gravitating systems and the accelerated expansion of the Universe without resorting to dark matter and/or dark energy [5,6,7,8,9,10,11,12,13,14]. On the other hand, they are also motivated by the fact that GR is not the quantum theory of gravity needed to describe the space-time near the singularities, which, as is well known, seemingly cannot be avoided [15]. Although a quantum theory of gravity should be able to overcome such problems, there also exists the possibility to avoid singularities modifying the coupling between matter and gravity. In this context, Eddington-inspired Born–Infeld (EiBI) gravity has been recently proposed [16]. EiBI gravity is inspired by the Born–Infeld action for non-linear electrodynamics, with the Ricci tensor replacing the field tensor \(F_{\mu \nu }\). This structure was motived by some classes of string theories where the Born–Infeld electrodynamics arises as a low-energy effective theory [17, 18]. One of the most interesting features is that EiBI is equivalent to GR in the vacuum while it introduces modifications in dense matter environments, where GR is experimentally not well probed. EiBI is able to describe, with only a single extra parameter (\(\kappa \)), astrophysical objects such as the Sun [19] and the internal structure of compact objects [20,21,22,23,24,25,26], and the cosmological expansion of the Universe [22, 27,28,29,30,31] (for comprehensive reviews see [7, 12] and the references therein).
Briefly, the gravitational action of EiBI gravity takes the following form:
where \(R_{\mu \nu }\) is the symmetric part of the Ricci tensor, \(\phi _M\) represents the matter field, and \(\lambda \) is a constant. The latter is linked to the cosmological constant in such a way that one obtains asymptotically flat solutions setting \(\lambda =1\). Finally, the field equations are built varying the action as in the Palatini approach. As in other modified theories of gravity, the Palatini approach is not equivalent to a pure metric one. However, the latter contains ghosts that can be eliminated only adding extra terms in the gravitational action [32, 33]. The higher order curvature terms account for both non-linear matter coupling and for avoiding singularities. Such correction terms also appear in the non-relativistic limit, where EiBI gravity leads to a modified Poisson equation given by
where \(\varPhi (\mathbf {r},t)\) is the gravitational potential, and \(\rho (\mathbf {r},t)\) is the matter density. Let us note that by setting \(\kappa =0\), the previous equation immediately reduces to the standard Poisson equation \({{\mathbf {\nabla }}^2}\varPhi (\mathbf {r},t) = 4\pi G\rho (\mathbf {r},t)\). The tightest constraint in the literature on the EiBI parameter has been obtained comparing the electromagnetic and gravitational interactions inside atomic nuclei: \( |\kappa | < 10^{-3} \, \mathrm{kg^{-1} \, m^5 \, s^{-2}}\) [37].
In this paper, we analyze the kinetic theory of the Jean instability for self-gravitating systems in EiBI gravity. This mechanism constitutes, on the theoretical side, a remarkable instrument to retain/rule out modified theories of gravity at astrophysical level. In fact, although a self-gravitating system collapses under the gravitational force induced by the modified Poisson equation and gives rise to star formation, we expect that the effects of the EiBI gravity are totally negligible in star formation environments while they must show some deviations from GR in compact objects collapsing into Black Holes.
The paper is organized as follows: in Sect. 2 we compute the dispersion relation for an homogeneous self-gravitating system in EiBI gravity; in Sect. 3 we analyze the dispersion relation in low and high frequency regimes, and we study the unstable modes that led to the collapse of the structure; in Sect. 4 we give our conclusion and remarks.
2 Dispersion relation of a collisionless self-gravitating system
The standard approach to describing the collapse of a self-gravitating system, either a star formation regions such as interstellar clouds which physical conditions change from hot X-ray emitting plasma to cold molecular gas, or compact objects collapsing into Black Hole, is the Jeans instability [34]. The latter is usually described by a distribution function of the particles \(f(\mathbf {r},\mathbf {v},t)\) which is solution of the Boltzman–Vlasov system of equations:
where the mass density distribution reads
Here \(I_{\mathrm{coll}}\) is the collision term and, since we analyze a collisionless system, it can be neglected (\(I_{\mathrm{coll}}=0\)).
To study the effect of the EiBI gravity, Eqs. (3)–(5) must be linearized. Assuming that the unperturbed potential \(\varPhi _0\) is locally constant in the system, one can set \(\nabla \varPhi _0 =0\). This is generally known as Jeans swindle [34]. Thus, in small perturbations regime, the distribution function and the gravitational potential can be written as
where \(\varepsilon \ll 1\) for small perturbations. Thus, at the first order, the Boltzman–Vlasov system of equations can be recast in the Fourier space as follows:
and the dispersion relation reads
Assuming the local thermodynamical equilibrium, the background distribution function of the particles can be described using the Maxwell–Boltzman distribution. Therefore, \(f_0(\mathbf {v})\) is given by
where \(\rho _0\) is the matter density at equilibrium, and \(\sigma \) is the thermal dispersion velocity of the particles. Finally, inserting the Maxwell–Boltzman distribution in Eq. (10) the dispersion relation reads
The previous equation shows a singularity at \(\omega = kv\). Moreover, setting \(\kappa =0\), it reduces to the Newtonian dispersion relation. Thus, one can infer the limit for the collapse setting \(\omega =0\) and computing the maximum wavelength of the perturbations supported by the system, and above which the system collapses. Specifically, in the Newtonian case, setting \(\omega =0\) one obtains the so-called Jeans wavenumber
that can be used to define the Jeans mass as the mass enclosed in a sphere of radius \(\lambda _J=2\pi /k_J\), obtaining
On the one hand, perturbations having wavelength \(\lambda >\lambda _J\) are unstable, and they exponentially grow. On the other hand, if the wavelength is less than the Jeans limit then perturbations are strongly damped. In EiBI gravity, such limit as a result is modified by the additional term in the Poisson equation. Therefore, from Eq. (12) we obtain
which corresponds to a mass
Thus, the mass limit for the collapse of a self-gravitating system depends on the EiBI parameter. It can be higher or lower than the classical Jeans Mass and, as a consequence, EiBI gravity can favor or disfavor the gravitational collapse depending on the sign of \(\kappa \).
3 Analysis of the collisionless dispersion relation
Equation (12) can easily be rewritten as
where we have defined the following variables:
Moreover, in EiBI gravity there naturally arises a new wavelength
which is equal to the one found at cosmological scale and using fluid approach [22]. The EiBI wavelength allows us to rewrite the dispersion relation in a more compact form
Let us study the limit of high frequency perturbations \(\beta \gg 1\).
In this case, having no singularities, we can integrate Eq. (20) along the real axis (\(\omega = \omega _R + i\omega _I \approx \omega _R\)) obtaining
which is a quadratic equation for \(\omega _R^2\). Remembering the condition \(\beta \gg 1\), we find
which for \(\kappa =0\) reduces to the classical Newtonian relation. Thus, as for the Newtonian case, high frequency perturbations with \(k^2 >k_J(\frac{k^2_{\mathrm{EiBI}}}{3 k^2_{\mathrm{EiBI}} + k^2_{J}})\) are quickly damped by the self-gravitating system whenever the EiBI wavenumber satisfies the condition \(k^2_{\mathrm{EiBI}} > - k^2_{J}/2\). On the contrary, high frequency perturbations can be supported by the system as is shown in Fig. 1 (magenta line).
More interesting is to consider the case of low frequency perturbations \(\beta \ll 1\), which in the Newtonian case led to unstable modes. In such a limit, the integral in Eq. (20) can be recast into the form
Therefore, splitting \(\omega \) in its real and imaginary part, \(\omega =\omega _R + \imath \omega _I\), and setting \(\omega _R=0\), we find
which gives unstable modes when \(\omega _I >0\), that is, for
In other words, the system supports perturbations having a wavelength
while it shows a singularity for perturbations having \(\lambda = \lambda _{\mathrm{EiBI}}\). Perturbations having \(\lambda < \lambda ^*\) are quickly damped by the system staying stable. Next, while in Newtonian gravity perturbations having \(\lambda > \lambda _{J}\) we are capable to generate the collapse of the system, in EiBI gravity the limit for the collapse becoming lower or higher depending on the parameter \(\kappa \). This can be quickly understood looking at Eq. (16). Since the EiBI modification depends on the coupling between matter and gravity, which is larger as the density increases [35], its effects must be negligible in low density environments like star formation regions where the condition \(\lambda _{\mathrm{EiBI}}\ll \lambda _{J}\) must hold. This condition can be straightforwardly translated into a density threshold,
which means that EiBI gravity does not affect systems having density below \(\rho ^*_0\). On the one hand, the interstellar medium has temperatures ranging from 10 to 100 K, while matter density ranges from \(\rho _0 \sim 10^{-18}\) to \(\rho _0 \sim 10^{-16}\, \hbox {kg} \,\hbox {m}^{-3}\). Setting \( |\kappa | < 10^{-3} \, \mathrm{kg^{-1} \, m^5 \, s^{-2}}\) [37] it follows that the density threshold is \(\rho ^*_0 \approx 10^7\, \hbox {kg}\, \hbox {m}^{-3}\). As expected, \(\rho _0 < \rho ^*_0\) in a star formation region. Therefore, EiBI gravity can be neglected. On the other hand, high density and temperature systems such as a Hyper Massive Neutron Star (HMNS, [36]) are one of the most promising laboratories to probe EiBI gravity. Those systems arise from the merger of a neutron star binary, and they are expected to collapse to a rotating black hole having a massive accretion torus. The temperature of those systems lies in the range \([5,50]\times 10^{10}\) K with a particle number density \(\sim 10^{39}\,\hbox {cm}^{-3}\). The collapse of such a system is considered as one of the possible sources of short gamma ray bursts. For such systems, with \(\rho ^*_0 \approx 10^{15}\, \hbox {kg}\, \hbox {m}^{-3}\) and \(\rho _0 \approx 10^{18}\, \hbox {kg}\, \hbox {m}^{-3}\), EiBI gravity affects the kinetic instability producing deviations from the GR.
Figures 2 and 3 display the growth rate of the unstable roots of Eq. (24) as a function of the normalized wavenumber \(k/k_J\) and for different values of the EiBI parameter \(\kappa \). Specifically, Fig. 2 illustrates the growth rate for a few cases corresponding to negative values of \(\kappa \), while Fig. 3 is devoted to positive values of \(\kappa \). Figure 2 illustrates that the growth rate is larger for lower values of \(\kappa \), while as \(\kappa \) is closer to zero (which corresponds to the Newtonian solution) the system shows unstable modes for \(k<k^*\), and stable modes for \(k>k^*\). Nevertheless, the Jeans mass as a result is higher (see Eq. (16)); therefore the collapse of the HMNS is disfavored with respect to the Newtonian case.
More interesting is the growth rate for positive values of \(\kappa \) shown in Fig. 3. In this case the dispersion relation has a singularity at \(k=k_{\mathrm{EiBI}}\). The physical behavior around such a point can only be described by more complex physical theories where a singularity does not occur. Therefore, the study here presented shows also the limit of the Jeans instability in EiBI gravity, which cannot be blindly applied in the points around \(k = k_{\mathrm{EiBI}}\).
Notice that such dramatic changes in \(\omega _{I}\) correspond to a phase transition of the system. In thermodynamic systems, phase transitions occur when the free energy, or the partition function, has singularities for some choice of variables, such as the temperature. For example, in the phenomenon of superconductivity, certain materials, when cooled below a critical temperature \(T_{C}\), pass into the superconducting state characterized by zero electrical resistance and by the complete ejection of magnetic field lines from the interior of the superconductor. In other systems, other physical parameters play the role of the temperature. For example, quantum phase transitions can be obtained by varying the magnetic field or the pressure at zero temperature.
In our case, a significant physical parameter is the wavenumber k and the critical value of k, which marks the phase transition is \( k_{\mathrm{EiBI}}\). Indeed, waves having wavenumbers smaller than \(k_{\mathrm{EiBI}}\) and satisfying the condition \(k<k^*\) give rise to unstable modes that favor the collapse of the HMNS into a black hole. Meanwhile, waves having wavenumber higher than the EiBI wavenumber show only stable solutions that do not produce the collapse of the structure even for \(k_{\mathrm{EiBI}}<k<k^*\).
We further remark that the extra modes \(k \sim k_{\mathrm{EiBI}}\) do not exist in Newtonian gravity where the singularity in \(k = k_{\mathrm{EiBI}}\) is absent in such a model (see Fig. 3), and also that this mode does not depend on the approximation made to compute the dispersion relation in Eq. (24) since it is straightforward to highlight that such a singularity is also present in the general expression of the dispersion relation given in Eq. (10). Therefore, this sort of phase transition could be an indication that the general paradigm of the Jeans instability, which works fine in GR and in other modified gravity models [38], does not work anymore in EiBI gravity for \(k>k_{\mathrm{EiBI}}\).
4 Conclusions and remarks
We have investigated the impact of the recently proposed EiBI gravity theory on the gravitational collapse of a self-gravitating system. We have solved the collisionless Boltzmann equation together with the modified Poisson equation of EiBI gravity to study the kinetic instability of self-gravitating system, and we have computed the corresponding dispersion relation leading to a new gravitational scale length. We studied both the high and low frequency limits of the dispersion relation. In the high frequency limit, the self-gravitating system behaves as in Newtonian gravity, not supporting the propagation of the perturbations. In the low frequency limit EiBI gravity may introduce a modification to the Jeans instability. Although EiBI gravity modified the Jeans mass needed for the collapse, it does not affect the star formation because of the low density environments where it happens. Nevertheless, in higher density environments such as HMNS the higher order terms of EiBI gravity produce a departure from the Newtonian growth rate, which could be both an indication that the standard Jeans paradigm does not hold anymore, or that effective new modes are present that could serve as a mechanism to generate a black hole from massive stars.
References
Planck Collaboration, A&A 594, A13 (2016)
J.L. Feng, Ann. Rev. Astron. Astrophys. 48, 495 (2010)
M. Yu, Khlopov Mod. Phys. Lett. A 32, 1702001 (2017)
A. Capolupo, Adv. High Energy Phys. 2016, 8089142 (2016)
S. Capozziello, M. De Laurentis, Phys. Rep. 509, 167 (2011)
S. Capozziello, M. De Laurentis, Ann. Phys. 524, 545 (2012)
T. Clifton, P.G. Ferreira, A. Padilla, C. Skordis, Phys. Rep. 513, 1 (2012)
I. de Martino, M. De Laurentis, F. Atrio-Barandela, S. Capozziello, MNRAS 442(2), 921–928 (2014)
I. de Martino, M. De Laurentis, S. Capozziello, Universe 1, 123 (2015)
Y.-F. Cai, S. Capozziello, M. De Laurentis, E.N. Saridakis, Rep. Prog. Phys. 79(4), 106901 (2016)
I. de Martino, Phys. Rev. D 93, 124043 (2016)
J. Beltran Jimenez, L. Heisenberg, G.J. Olmo, D. Rubiera-Garcia (2017). arXiv:1704.03351
S. Nojiri, S.D. Odintsov, Phys. Rept. 505, 59–144 (2011). https://doi.org/10.1016/j.physrep.2011.04.001
S. Nojiri, S.D. Odintsov, V.K. Oikonomou, Phys. Rept. 692, 1–104 (2017). https://doi.org/10.1016/j.physrep.2017.06.001. arXiv:1705.11098
S.W. Hawking, G.F.R. Ellis, The Large Scale Structure of Space-Time (Cambridge University Press, Cambridge, 1973)
M. Banados, P.G. Ferreira, Phys. Rev. Lett. 105, 011101 (2010)
M. Born, L. Infeld, Proc. R. Soc. A 144, 425–451 (1934)
E. Fradkin, A.A. Tseytlin, Phys. Lett. B 163, 123 (1985)
J. Casanellas, P. Pani, I. Lopes, V. Cardoso, Astrophys. J. 745, 15 (2012)
P. Pani, V. Cardoso, T. Delsate, Phys. Rev. Lett. 107, 031101 (2011)
P. Pani, T. Delsate, V. Cardoso, Phys. Rev. D 85, 084020 (2012)
P.P. Avelino, Phys. Rev. D 85, 104053 (2012)
Y.-H. Sham, L.-M. Lin, P.T. Leung, Phys. Rev. D 86, 064015 (2012)
Y.-H. Sham, P.T. Leung, L.-M. Lin, Phys. Rev. D 87, 061503(R) (2013)
T. Harko, F.S.N. Lobo, M.K. Mak, S.V. Sushkov, Phys. Rev. D 88, 044032 (2013)
H. Sotami, Phys. Rev. D 89, 104005 (2014)
P.P. Avelino, R.Z. Ferreira, Phys. Rev. D 86, 041501 (2012)
A. De Felice, B. Gumjudpai, S. Jhingan, Phys. Rev. D 86, 043525 (2012)
M. Bouhmadi-López, I. Albarran, C.-Y. Chen, Universe 3(2), 36 (2017)
I. Albarran, M. Bouhmadi-López, C.-Y. Chen, P. Chen, Phys. Lett. B 772, 814–818 (2017)
M. Bouhmadi-López , C.-Y. Chen, P. Chen, Eur. Phys. J. C 75, 90 (2015)
S. Deser, G.W. Gibbons, Class. Quantum Gravity 15, L35 (1998)
D.N. Vollick, Phys. Rev. D 69, 064030 (2004)
J. Binney, S. Tremaine, Galactic Dynamics (Princeton University Press, Princeton, 1994)
T. Delsate, J. Steinhoff, Phys. Rev. Lett. 109, 021101 (2012)
L. Rezzolla, O. Zanotti, Relativistic Hydrodynamics (Oxford University Press, Oxford, 2013)
P.P. Avelino, JCAP 11, 022 (2012)
S. Capozziello, M. de Laurentis, I. de Martino, M. Formisano, S.D. Odintsov, Phys. Rev. D 85(4), 044022 (2012)
Acknowledgements
IDM acknowledges financial supports from University of the Basque Country UPV/EHU under the program “Convocatoria de contratación para la especialización de personal investigador doctor en la UPV/EHU 2015”, and from the Spanish Ministerio de Economía y Competitividad through the research project FIS2010-15492, and from the Basque Government through the research project IT-956-16. AC acknowledges partial financial support from MIUR and INFN. The authors acknowledges the COST Action CA1511 Cosmology and Astrophysics Network for Theoretical Advances and Training Actions (CANTATA).
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
De Martino, I., Capolupo, A. Kinetic theory of Jean instability in Eddington-inspired Born–Infeld gravity. Eur. Phys. J. C 77, 715 (2017). https://doi.org/10.1140/epjc/s10052-017-5300-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-017-5300-0