Brown dwarfs in Eddington-inspired Born-Infeld and beyond Horndeski theories

We studied herein the mass and the radius of brown dwarfs predicted by beyond Horndeski (BH) and Eddington-inspired Born-Infeld (EiBI) gravity theories by numerically solving the modified non-relativistic hydrostatic equations of both theories. We used a recent compilation of brown dwarf masses and radii obtained from Ref. Bayliss et al. (Astrophys J 153:1, 2016) to constrain the free parameter of both theories. We obtain the range of the corresponding parameters with 1σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} and 5σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} confidence by using chi-squared analysis. Furthermore, the minimum chi-squared values can be reached for the cases of κ=0.17×102m5kg-1s-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa = 0.17 \times 10^2 ~\mathrm{m}^5\,\mathrm{kg}^{-1}\,\mathrm{s}^{-2}$$\end{document} and γ=-0.1207\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma = -0.1207$$\end{document} for EiBI and BH theories, respectively. The corresponding parameter values with the minimum chi-squared values are relatively small; therefore, they cannot significantly change the brown dwarf mass limits determined from the equivalence of nuclear and photosphere luminosities for the pp (hydrogen burning) and pp+pd (deuterium burning) reactions.


Introduction
Einstein's Theory of General Relativity (GR) is one of the most successful theories in physics. This theory has passed all precision tests in a medium-energy scale, and has a substantial conceptual foundation. However, this theory is still not perfect. In a low-energy scale, additional ingredients, such as dark matter and a nonzero cosmological constant (dark energy) should be introduced to fit with a variety of astrophysical observations of galaxies and clusters as well as with the cosmic acceleration of the universe. These ingredients raise two additional questions related to the cosmological constant and coincidence problems. In a high-energy scale, GR as a classical theory, also yields a singularity, which is a e-mail: anto.sulaksono@sci.ui.ac.id slightly unpleasant for some people. Therefore, until now, much interest has been observed in the study of alternative or modified theories of GR (Please see the reports of a recent review in modified gravity theories [2][3][4][5][6][7] and the references therein). After passing very stringent sets of some commonly accepted fundamental physical requirements, BH and EiBI theories can now be classified as plausible alternatives or modified gravity theories in the context of strong gravitation field tests (Please see details, in Ref. [2].). The EiBI and a broad class of BH theories of the G 3 type yield modifications of the Newtonian hydrostatic equilibrium equations. These modifications affect the properties of non-relativistic stars, such as white dwarfs (WDs), main sequences, and brown dwarfs (BDs). We need to note some recent works concerning the BH and EiBI theories and their applications in astronomical and cosmology objects. For BH theory, they can be found in Refs. [8][9][10][11][12][13][14][15][16][17][18] while those for EiBI theory, can be found in Refs. [19][20][21][22][23][24][25][26][27][28][29]). We also need to note that the recent review of the cosmological test of modified gravities has been discussed in Refs. [30,31], while the discussion about cosmological constraint on Horndeski gravity in the light of GW170817 can be found in Ref. [32].
BDs were theoretically predicted for the first time in 1963 by Kumar [33] and Hayashi and Nakano [34]. BDs are known as faint self-gravitating bodies. A BD generally has a sufficiently low mass, and its radius is almost independent to the mass. However, the BD elemental composition is substantially uniform. A BD's core is composed of hydrogen, helium, and other chemical compositions, and its body is mainly supported by a repulsive force from degenerate electrons with respect to gravitational collapse. However, detailed structure of the atmosphere and the surface luminosity of BDs is quite complex (Please see a recent review of the BD theory in [35][36][37] and the references therein.). A BD was detected for the first time more than 30 years later after its prediction. Since then hundreds of BDs have been observed (Please see details in Ref. [38] and the references therein.). Much progress has recently been made and reported from the observation investigations (Please see Ref. [1] and the references therein.). For example, the most precise masses and radii of BDs (LHS 6343A and LHS6343C) have been reported in [39,40] ( i.e., BD with a mass of 62.7 ± 2.4M J up and a corresponding radius of 0.833 ± 0.021R J up and BD with a mass of 62.1 ± 1.2M J up and a radius of 0.782 ± 0.013R J up (M = 10 3 M J up and R = 9.735R J up )). Note that, these results depended on the calculated mass with a precision of 1.9% and on the calculated radius with an accuracy of 1.4% [1,39,40]. One crucial challenge is to distinguish BDs and very low-mass stars (VLMs) from their maximum mass because from the BD models, the features are degenerate with their age, radius, and metallicity [41]. The maximum mass of BDs is known as the hydrogen-burning limit maximum mass (MMHB), where the corresponding mass is approximately ∼ 0.08M . The MMHB shows the limit where the collapse will not occur as far as the star's surface luminosity is still supplied by nuclear fusion. The precise value of the MMHB depends, but not significantly, on the corresponding elemental composition. A recent study using analytical models within GR theory [42] has shown an acceptable MMHB in the range of 0.064 − 0.087M . BD Kepler-503b has also recently reported to sit right at the hydrogen-burning limit mass, straddling the boundary between BDs and VLMs [43]. Another study [44] reported that objects that behave exactly like BDs are evidence, with their masses exceeding MMHB. These objects, which are termed "overmassive BDs," form a continuous sequence with traditional BDs in any property, such as mass, effective temperature, radius, and luminosity. Similar to the MMHB limit, a distinction between the planetary-mass objects below and the BDs above is also needed. The corresponding minumum mass limit of a BD is known as the deuterium-burning limit mass (MMDB). The MMDB value is generally ∼ 0.013M . The MMDB strongly depends on the elemental composition, such as helium abundance, the initial deuterium abundance, and model metallicity. Even though the MMDB value is ∼ 0.013M for most BDs, models ranging from zero metallicity to more than three times of solar metallicity have recently been reported to have a corresponding MMDB ranging from ∼ 0.011 to ∼ 0.016M [45]. Moreover, note that a recent study on deuterium burning in objects forming through the core accretion scenario has shown that the value of the MMDB (0.013M ) does not change much when none of the corresponding investigated parameters was able to change this mass limit by more than 0.0008M . Furthermore, the luminosity of hot and cold stars caused by hydrogen burning becomes comparable only after ∼ 200 Myr [46].
Studies on WDs within modified gravity theories, including the EiBI and BH theories have been done by many authors (please see Refs. [12,18,29,47,48] for examples).
Constraining the parameters of several modified gravity theories, such as the EiBI, BH, scalar-tensor-vector gravity, f(R), and forth-order gravity theories, WD masses and radii data from observations have also been used and reported in Refs. [12,18,29,48]. These studies, used chi-squared analysis to determine the allowed region of parameter space with its corresponding confidence. This analysis is very robust determining the acceptable parameter of the modified gravity theories from the observation data of WD properties. However, such study has not been done before for the case of BDs because of the limited masses and radii data. Therefore, it is interesting to know how robust the BD properties from the masses and radii observation data in extracting the acceptable parameter of modified gravity theories. Sakstein [10,16] recently studied the impact of modified gravity within BH theory on the MMHB. He found that the MMHB increases when the free parameter γ value of BH theory increases. He stated that the MMHB of BDs has the potential to provide an extremely stringent limit to constraint the γ value of BH theory.
The present study constrains the parameter of the EiBI and BH theories using the current masses and radii observation data of BDs [1]. We also use Chi-squared analysis to determine the allowed region of parameter space of both theories and the acceptable parameters with their confidences. The obtained free parameters are used to calculate the mass limits determined from the equivalence of nuclear and photon sphere luminosities for the pp (hydrogen burning) and pp+pd (deuterium burning) reactions predicted by the EiBI and HB theories and compared to the ones of GR. The results can be used to study the impact of both modified gravity theories to the BD mass limits from hydrogen and deuterium burning. Section 2, briefly discusses the BD equation of states (EOS), which is the modified hydrostatic equilibrium equation of both theories presented. We also discuss the chisquared analysis on BDs and the EOS profile of both theories. Section 3 presents the hydrogen and deuterium-burning mass. Finally, Sect. 4 provides the conclusions.

BDs in the EiBI and BH theories
We need to solve a modified hydrostatic equilibrium equation to calculate the mass and radius of non-relativistic stars, such as BDs, within a modified gravity theory. However, the EOS of stars as input for the hydrostatic equilibrium equation is still similar to that used for Newtonian. We discuss herein the EOS of BDs and the modification of the hydrostatic equilibrium equations within the EiBI and BH theories, chi-squared test of the mass and radius, and the corresponding EOS profiles of a BD.

Equation of states
The density of BDs is estimated between the range of 10 3 kg/m 3 and 10 6 kg/m 3 [48]. BDs preserve their lives from gravitational collapse mostly from the contribution of a repulsive Fermi pressure from degenerate electrons. The corresponding Fermi pressure in a non-relativistic limit is presented as follows: and and μ are energy density and chemical potential of electrons. K B and T are the Boltzmann constant and temperature, respectively. The P F integral equation can be analytically solved and compactly expressed using dilogaritmic Li 2 and trilogaritmic Li 3 functions. The result is as follows: We define the degeneracy parameter ψ for simplification. This parameter can link the central density and the central temperature of a BD and make the expression of the EOS in a polytropic form, then ψ is defined as: where μ F is the Fermi electron energy in the degeneracy limit, μ e = 1.143, 1 μ e = X + Y 2 ; X and Y are mass fractions from hydrogen and helium, respectively. m e and m H are the electron and hydrogen masses, respectively. Parameter ψ is a function of temperature and density. The BDs interior consisted of hydrogen and helium hence the total pressure has the combination effects of electrons and ions, where P ion = κρT μ 1 m H and μ 1 is the average molecular mass for the hydrogen and helium mixture. The explicit expression for the total pressure P = P F + P ion can be written as: with now α is 5μ e 2μ 1 . Equation (4) can now be recast into a polytropic form as: where Γ = 1 + 1 n with n = 3/2. K takes the following form: The term Υ is defined herein as: We use the polytropic expression of the EOS to convert the form of the hydrostatic equilibrium equation from a gradient of the pressure form into a gradient of the density form.

BH and EiBI gravity theories
BH theory is a modified gravity theory containing the screening mechanism if R << r V . r V is the Vainshtein radius and R is the radius of the corresponding astrophysical body. In other words r V defined as the transition between the screened and unscreened regimes inside astrophysical bodies (Please see the discussion of the Vainshtein mechanism in Ref. [14] and the references therein.). The action of BH theory is presented as [14]: π(r, t) and use the Newtonian form of Friedmann-Robertson-Walker (FRW) space-time metric.
with r 2 = δ i j x i x j and Φ(r, t) as well as Ψ (r, t) are time dependent gravitation potentials, the non-relativistic limit of the static gravitation potential Φ(r ) of BH theory derived from Eq. (8) yields a modified Poisson equation. Note that the BH dimensionless parameter is defined as γ = 2 Λ 4 with =φ 0 Λ 4 . The expression of the corresponding modified Poisson equation of BH theory is presented as [16]: Please see the detail derivation in Ref. [14]. We can extract the radial acceleration from the Poisson equation in Eq. (10) to obtain the Euler equation. From the Euler equation, we can obtain the hydrostatic equilibrium equation represented as: The EiBI theory of gravity belongs to one type of Born-Infeld gravity modification introduced by Banados and Ferreira [19], who used the auxiliary field approach. This theory can regularize the gravitational dynamics, leading to nonsingular cosmologies, and predict regular black hole spacetime without resorting to quantum gravity. The EiBI action is presented as [19]: where |.| denotes a determinant. Here κ is a free parameter of the EiBI theory. S m (g ab , χ m ) is the matter action, χ m denotes the matter field, g ab is a metric tensor and R ab is the symmetric part of the Ricci tensor, while λ is related to the cosmological constant, Λ = (λ−1) κ so that we can obtain the asymptotically flat solutions when λ = 1. The EiBI modified Poisson equation is obtained as follows from the non-relativistic limit of the EiBI field equation [19,20]: Similar to that of BH theory, we can easily arrive to the hydrostatic equilibrium equation: The gradients of the mass of both theories are the same in the non-relativistic limit: Together with Eq. (11) for BH theory or Eq. (14) for EiBI theory and using Eq. (5) as input, Eq. (15) can be numerically solved using the fourth-order Runge-Kutta method. We used herein M(≈ 0) ≈ 0, ρ(≈ 0) = ρ c and ρ(R) ≈ 0 as the the boundary conditions. The upper panel of Fig. 1 shows the mass-radius relation of the BDs for the current epoch predicted by BH theory, while the top panel of Fig. 2 depicts that for EiBI theory. The corresponding figures also illustrate those predicted by Newtonian and the observation data taken from Ref. [1] for comparison. The γ parameter can be positive or negative for BH theory. For γ > 0, the gravity correction term in the hydrostatic equation plays a role as an additional repulsive force. Meanwhile, for γ < 0, it plays a position as another additional attractive force. As a result, the mass and radius relation is modified because of the role of γ , where for the fixed mass value, the radius is increased or decreased by increasing or decreasing the γ value. This situation is also applicable to EiBI theory, where the role of the dimensionless parameter γ in BH is replaced by the role of parameter κ. However, the trend of the mass changes with respect to the radii of both theories because the functional form of the gravity correction of both theories is different. Note that for the mass-radius data taken from the observation [1], three points have rather a large error bar.

Chi-squared analysis and mass-radius relation
In this section, we determine whether or not the appropriate parameter values of the BH and EiBI theories which are the masses and the radii predictions, are compatible with the 12 observational masses and radii of BD data in Ref. [1] using chi-squared analysis. We selected one point from our theoretical mass and radius relation results of each theory with a variety of parameter values that have smaller deviation to the mass and radius observational data point of Ref. [1]. For each theory, we also collected all masses and radii having a small deviation together with the corresponding data from the observation. Subsequently, we used the following equation to calculate the chi-square (Δχ 2 i ) of each BD: where M th is the mass from the theoretical calculation point with a small deviation to M i , which is the mass of the corresponding observational of the ith BD. R th is the theoretical radius, and R i is the radius of observational ith BD. σ 2 M,i and σ 2 R,i are the corresponding mass and radius standard deviations, respectively. We then summed up the results of Eq. (16) for all 12 BDs using the following equation: To this end, we obtained the chi-squared result in every free parameter of the used modified gravity theories shown in the lower panels of Figs. 1 and 2 for the BH and EiBI theories, respectively. Figure 2 shows that for the EiBI model, the chi-squared limits of the corresponding parameter were in the range of −2.375 ≤ κρ A ≤ 1.266 or (−1.51 × 10 2 ≤ κ ≤ 0.81 × 10 2 ) m 5 kg −1 s −2 for the 1σ confidence level and −2.5 ≤ κρ A ≤ 1.822 or (−1.59 × 10 2 ≤ κ ≤ 1.16 × 10 2 ) m 5 kg −1 s −2 for the 5σ confidence level. The minimum of the chi-square was κρ A = 0.274 or κ = 0.17 × 10 2 m 5 kg −1 s −2 , which was quite close to that of Newtonian. Figure 1 also depicts that the chi-squared limits of the γ parameter for BH theory were in the range of −0.565 ≤ γ ≤ 0.234 for the 1σ confidence level and −0.6 ≤ γ ≤ 0.391 for the 5σ confidence level. This theory showed two minimums but the true minimum of the chisquare was for γ = −0.12072. This result was also quite close to that of Newtonian. However, note that the three data points taken from Baylisset al. [1] had quite a large error bar in radius. If in the future we have more data points with the values of the average error bar smaller than the one used herein, we will have a better prediction of γ and κ based on the BD properties. For comparison for EiBI theory, we need to note that the authors of Ref. [29] obtained from WDs mass and radius constraints, a quite close κ range to the ones of BDs i.e., (−16.0 × 10 2 < κ < 0.35 × 10 2 ) m 5 kg −1 s −2 for the 1σ confidence level. From neutron stars properties, it is reported by the authors of Ref. [21]that κ 10 −3 m 5 kg −1 s −2 and using the recent maximum mass constraints, it is reported a more restricted κ constraint i.e., (2.7 × 10 −4 < κ < 0.35 × 10 −4 ) m 5 kg −1 s −2 [24]. This leads to the fact that constraints from NS are tighter than those of WD and BD. While stellar equilibrium and solar constraints lead to |κ| 3 × 10 5 m 5 kg −1 s −2 [26]. |κ| 6 × 10 8 m 5 kg −1 s −2 constraint comes from primordial nucleosynthesis [23]. It is also reported recently that in EiBI theory, the speed of gravitation waves in matter deviates to c. Therefore, from the time delay in the arrival gravitation wave signals at Earth-based detectors, can be estimated that |κ| 10 11 m 5 kg −1 s −2 , while from time delay between the signals of GW170817 and GRB170817A in a background FRW universe |κ| 10 27 m 5 kg −1 s −2 [28]. Therefore, based on the variety predictions of free parameter of EiBI theory κ, we might suspect that the values of κ could depend on the compactness of the corresponding stellar and cosmological objects and compact objects provide tighter constraint than those of cosmological objects.
For the case of BH theory, several works [10,[12][13][14]16] have assessed the behavior of less compact stars using a nonrelativistic form of BH theory. They have found the corresponding γ constraint is −0.51 < γ < 0.027 at red-shift zero. γ = −0.51 lower limit comes from the consistency of the lowest mass of WD with Chandrasekhar mass [12], while γ = 0.027 upper limit comes from the minimum mass for hydrogen burning [10,16] (Please see also the discussion of this γ constraint range in Ref. [17]). The constraints from BDs with 1σ confidence is quite compatible with this γ constraint range. Furthermore, using the most recent, complete, and independent measurements of masses and radii of WDs as well as using a realistic EOS of WD, the authors of Ref. [18] has found tighter constraint than that of the upper limit from Hydrogen burning, i.e., γ < 0.14. Furthermore, the analysis of weak lensing and X-ray profiles of 58 galaxy clusters with an averaged red-shift of 0.33 have done by the authors of Ref. [17] obtain the constraint γ = −0.11 0.93 −0.67 . It is interesting to observe that the true minimum of the chisquare obtained in our result is quite compatible with this galaxy cluster constraint. Note that γ is related to the parameters appearing in the effective field theory (EFT) of dark energy. The EFT parameters characterize the linear cosmology of BH theory. Therefore, γ constrain deviation from GR on cosmological scale [17]. We also need to note that in the case of neutron stars, the authors of Ref. [8] have found that the configuration with γ −0.05 predict reasonable maximum masses M ∼ 3 M and radii R 14 km.

Brown dwarf profiles
For completeness, we also provide herein the contour plots of the EOS profiles (ρ-T-P plots) for a certain ρ c in the range of 10 4 −10 6 kg/m 3 [37] in Fig. 3 as predicted by the BH and EiBI theories using the γ and κ values with the minimum of the chi-square. The results of Newtonian were also given for comparison. The patterns of the results of the EiBI, and BH theories and Newtonian are quite similar because we used quite small values of γ for BH and κ for EiBI. The profiles will look more different if we used relatively large values of γ and κ.

Brown dwarf limiting masses
One must clearly determine the mass boundaries of WDs to determine the mass range where the WDs exist. The maximum and minimum masses of WDs can be used to differentiate BDs from other stars or planets. BDs maintain their shapes with respect to gravitational collapse using the pressure degeneracy of electrons and the role of any other ions. The core burning of hydrogen and helium significantly acts to determine the nuclear luminosity. The nuclear luminosity limit needs to be equal to the photosphere luminosity to prevent balance. The requirement is L N = L e . In this section, we want to know how further modifying the gravity will shift the limiting masses of BDs.

Nuclear burning luminosity
The primary reaction to differentiating BD and VLS is p + p → d + e + + ν e . The next important reaction is p + d → 3 He + γ . One can always parameterize the thermonuclear rates to power laws in T and ρ. This parameterization is valid only for a narrow range in T and ρ, but is with s = 6.31, u = 2.28, and c = pp for hydrogen burning (Please see details in [10]) and c = pp + pd for deuterium burning (e.g., discussed in [44]). The detailed process explains that the energy generation rates depend on T and ρ [37]. We use the following expression of the energy generation rates taken from Refs. [37,44]: where T 6 is defined as T /(10 6 K ). Deuterium is produced by the weak interaction in the p-p reaction. The two reactions mentioned in the earlier section are accountable for the deuterium chain production. The equilibrium value for deu- terium production is presented as [44]: where Q pd = 5.494 MeV and Q pp = 1.442−0.262 MeV = 1.18 MeV. The p-p reaction also shows that the average neutrino energy lost is 0.262 MeV [36].

Luminosity of hydrogen and deuterium burning
We can obtain the nuclear luminosity as follows: This becomes the following when we substitute the gradient mass into Eq. (22) with T /T c = (ρ/ρ c ) 2/3 . This equation will be solved together with other differential equations to obtain L N . Note that we use L N (r ≈ 0) ≈ 0 as an initial boundary [49].

Photosphere luminosity
When comparing the star radius, the surface can be really thin. The surface gravity g = G M/r 2 . We will use the approximation from d M/dr = 2M/r and d 2 M/dr 2 = 2M/r 2 [10] to find the modified surface gravity g e f f from the BH and EiBI gravity theories. Accordingly, g e f f in BH theory is given as follows: Meanwhile in the EiBI theory g e f f is obtain as: We used the ideal gas law P e = ρ e k B T e μm H , where the mean molecular mass μ = 0.593 and we can obtain the following equation using the polytropic EOS to find the effective temperature in the surface  with k R = 10 −2 cm 2 /g. With this effective temperature in hand we can calculate the photo sphere luminosity as follows: Figure 4 shows the nuclear and photo sphere luminosities as a function of the BD mass predicted by the BH and EiBI theories using κ and γ with the minimum chi-squared values. Those of GR with their corresponding hydrogenand deuterium-limiting mass boundaries were also given for comparison. Figure 4 further illustrates that for those of GR, our hydrogen-limiting mass result was compatible with the known standard MMHB value. However, our deuteriumburning limit mass (M (DB limit)) results were larger than those of the standard MMDB value. Therefore, the calculation of M (DB limit) using this procedure only roughly approximates the MMDB. A further refined calculation, such as a more realistic EOS, and a better approximation for L N and Le are needed to reproduce a correct MMDB value. By using κ and γ with the minimum chi-squared values, we found that the corresponding hydrogen-and deuteriumburning limits were not significantly changed. However, for EiBI theory, the range slightly shifted to a larger mass. In contrast, for BH, it slightly shifted to a smaller mass. Figure 5 shows the masses and the radii of the hydrogen-and deuterium-burning mass limits predicted by BH theory as a function of parameter γ . For that of hydrogen-burning, our mass and radius results were compatible to those obtained in Ref. [10]. the radii and the masses of hydrogen-and deuterium-burning limits increased with the increasing γ . Figure 6 shows those of EiBI theory. A similar trend was also evidenced from Fig. 6, wherein the radii of the hydrogen-and deuterium-burning mass limits increased with the increasing κ. The impact of the contribution of the corresponding modified gravity terms more significantly appeared in BH theory than in EiBI theory for the case of κ and γ with the minimum chi-squared values. Figures 5 and 6 also show the impacts of the maximum and minimum γ and κ values with the 1σ and 5σ confidences in the masses and the radii of the hydrogenand deuterium-burning mass limits. The maximum allowed parameter κ value (the one within 1σ confidence) of EiBI had a more significant effect on the masses of the hydrogenand deuterium-burning limits compared to the corresponding γ value of BH.

Conclusion
Using compact and less-compact stars to study the modified gravity theory is important. In this way we can check the possible deviation signature of the standard gravity theory and the dependency of the corresponding free parameter of modified gravity theories with the compactness of the stellar and cosmological objects. We studied herein the mass and the radius of BDs as a representation of lesscompact stars predicted by BH and EiBI gravity theories by numerically solving the modified hydrostatic equations predicted by both theories. We used the most recent compilation of BD masses and radii obtained from Ref. [1] to constrain the free parameter of the BH and EiBI theories using a minimum chi-squared analysis. We found that for EiBI theory, the chi-squared limit κ parameter is in the range of (−1.51 × 10 2 ≤ κ ≤ 0.81 × 10 2 ) m 5 kg −1 s −2 for the 1σ confidence level and in the range of (−1.59 × 10 2 ≤ κ ≤ 1.16 × 10 2 ) m 5 kg −1 s −2 for the 5σ confidence level. Meanwhile for BH theory, the chi-squared limit γ parameter is in the range of −0.565 ≤ γ ≤ 0.234 for the 1σ confidence level and in the range of −0.6 ≤ γ ≤ 0.391 for the 5σ confidence level, while κ = 0.17 × 10 2 m 5 kg −1 s −2 for the EiBI and γ =-0.1207 for the BH theories can yield minimum chi-squared values. The value of κ and γ parameters with minimal χ 2 were relatively small. Therefore, they cannot significantly provide effects to increase or decrease the MMHB and MMDB limits of a BD. We also obtained the hydrogen and approximate deuterium mass limits by calculating the BD nuclear and surface luminosity using c = pp for hydrogen burning and c = pp + pd for deuterium-burning reactions. Also, we have found that the EiBI free parameter κ constraint from BDs is tighter than that from cosmological objects, while the BH free parameter γ constraint from BDs is compatible with the one from cosmological objects.