Charged Quark Stars and Extreme Compact Objects in Regularized 4D Einstein-Gauss-Bonnet Gravity (2024)

Michael Gammonmgammon@uwaterloo.caDepartment of Physics and Astronomy, University of Waterloo, Waterloo, Ontario, Canada, N2L 3G1 Robert B. Mannrbmann@uwaterloo.caDepartment of Physics and Astronomy, University of Waterloo, Waterloo, Ontario, Canada, N2L 3G1 Sarah Rourkesarah.a.rourke@mail.mcgill.caDepartment of Physics, McGill University, Montreal, Quebec, Canada

Abstract

Since the derivation of a well-defined $D\rightarrow 4$ limit for 4dimensional Einstein Gauss-Bonnet (4DEGB) gravity coupled to a scalar field, there has been interest in testing it as an alternative to Einstein’s general theory of relativity. Using the Tolman-Oppenheimer-Volkoff (TOV) equations modified for charge and 4DEGB gravity, we model the stellar structure of charged, non-interacting quark stars. We find that increasing the Gauss-Bonnet coupling constant $\alpha$ or the charge $Q$ both tend to increase the mass-radius profiles of quark stars described by this theory, allowing a given central pressure to support larger quark stars in general. We also derive a generalization of the Buchdahl bound for charged stars in 4DEGB gravity. As in the uncharged case, we find that quark stars can exist below the general relativistic Buchdahl bound (BB) and Schwarzschild radius $R=2M$ , due to the lack of a mass gap between black holes and compact stars in the 4DEGB theory. Even for $\alpha$ well within current observational constraints, we find that quark star solutions in this theory can describe Extreme Compact Charged Objects (ECCOs), objects whose radii are smaller than what is allowed by general relativity.

I Introduction

Modified theories of gravity continue to attract attention despite the empirical success of general relativity (GR). These theories are motivated by a variety of problems, including addressing issues in modern cosmology bueno2016 ; sotiriou_2010 ; nojiribook ; clifton_2012 , quantizing gravity ahmed2017 ; stelle1977 , eliminating singularities brandenberger1992nonsingular ; BRANDENBERGER_1993 ; brandenberger1995implementing , and, perhaps most importantly, finding viable phenomenological competitors against which GR can be tested in the most stringent manner possible.

Higher curvature theories (or HCTs) are amongst the most popular modifications. An HCT modifies the assumed linear relationship in GR between the curvature and the stress-energy, replacing the former with an arbitrary sum of powers of the curvature tensor (appropriately contracted to two indices). Such modificationsprovide us with a foil to further challenge the empirical success of GR, while also making new testable predictions.

Lovelock theories lovelock1971 have long been at the forefront of this search, since they possess the distinctive feature of having 2nd order differential equations of motion. The physical significance of such theories has been unclear, however, since their higher order terms yield non-trivial contributions to the equations of motion only in more than four spacetime dimensions ( $D>4$ ).

Recently this restriction was circumvented hennigar_2020_on ; Fernandes:2020nbq in the quadratic case, or what is better known as “Einstein-Gauss-Bonnet” gravity . The Gauss-Bonnet (GB) contribution to the gravitational action is

S_{D}^{GB}=\alpha\int d^{D}x\sqrt{-g}\left[R^{\mu\nu\rho\tau}R_{\mu\nu\rho\tau%}-4R^{\mu\nu}R_{\mu\nu}+R^{2}\right]\equiv\alpha\int d^{D}x\sqrt{-g}\mathcal{G}

(1)

(where $R_{\mu\nu\rho\tau}$ is the Riemann curvature tensor), which becomes the integral of a total derivative in $D=4$ , and thus cannot contribute to a system’s gravitational dynamics in less than five dimensions. For this reason it is often referred to as a “topological term” having no relevance to physical problems. Indeed, the Lovelock theorem lovelock1971 ensures that a $D=4$ dimensional metric theory of gravity must incorporate additional fields in order to have second order equations of motion and diffeomorphism invariance.

A few years ago it was noted glavan2020 that several exact solutions to $D$ -dimensional Einstein-Gauss Bonnet gravity have a sensible limit under the rescaling

\lim_{D\to 4}(D-4)\alpha\rightarrow\alpha,

(2)

of the Gauss-Bonnet coupling constant. Using this approach a variety of 4-dimensional metrics can be obtained, including cosmological glavan2020 ; li2020 ; kobayashi2020 , spherical black hole glavan2020 ; kumar2022 ; fernandes2020 ; kumar2020 ; kumar2022_2 , collapsing malafarina2020 , star-like doneva2021 ; charmousis2022 , and radiating ghosh2020 metrics, each carrying imprints of the quadratic curvature effects of their $D>4$ counterparts. However a number of objections to this approach were subsequently raised gurses2020 ; ai2020 ; shu2020 , based on the fact thatthe existence of a limiting solution does not imply the existence of a well-defined 4D theory whose field equations have that solution. This shortcoming was quickly addressed when it was shown that the $D\to 4$ limit in (2) can be taken in the gravitational action hennigar_2020_on ; Fernandes:2020nbq , generalizing an earlier procedure employed in obtaining the $D\to 2$ limit of GR Mann:1992ar . One can also compactify $D$ -dimensional Gauss-Bonnet gravity on a $(D-4)$ -dimensional maximally symmetric space and then use (2) to obtain a $D=4$ HCT Lu:2020iav . This approach yields the same result (up to trivial field redefinitions), in addition to terms depending on the curvature of the maximally symmetric $(D-4)$ -dimensional space. Taking this to vanish yields

\displaystyle S_{4}^{GB}=

\displaystyle\alpha\int d^{4}x\sqrt{-g}\left[\phi\mathcal{G}+4G_{\mu\nu}\nabla%^{\mu}\phi\nabla^{\nu}\phi-4(\nabla\phi)^{2}\square\phi+2(\nabla\phi)^{4}\right]

(3)

where we see that an additional scalar field $\phi$ appears. Surprisingly, the spherically symmetric black hole solutions to the field equations match those from the naïve $D\rightarrow 4$ limit of $D>4$ solutions glavan2020 . The resultant 4D scalar-tensor theory is a particular type of Horndeski theory horndeski1974 , and solutions to its equations of motion can be obtained without ever referencing a higher dimensional spacetime hennigar_2020_on .

We are interested here in what is called 4D Einstein-Gauss-Bonnet gravity (4DEGB), whose action is given by (3) plus the Einstein-Hilbert term:

	$\displaystyle S$	$\displaystyle=S^{GR}+S_{4}^{GB}$
		$\displaystyle=\frac{1}{8\pi G}\int d^{4}x\sqrt{-g}\left[R+\alpha\left\{\phi%\mathcal{G}+4G_{\mu\nu}\nabla^{\mu}\phi\nabla^{\nu}\phi-4(\nabla\phi)^{2}%\square\phi+2(\nabla\phi)^{4}\right\}\right]$		(4)

which has been shown to be an interesting phenomenological competitor to GR Clifton:2020xhc ; charmousis2022 ; Zanoletti:2023ori .Despite much exploration Fernandes:2022zrq , the role played by thesehigher curvature terms in real gravitational dynamics is still not fully understood. One important arena for testing such theories against standard general relativity is via observations of compact astrophysical objects like neutron stars. The correct theory should be able to accurately describe recent gravitational wave observations of astrophysical objects existing in the mass gap between the heaviest compact stars and the lightest black holes.

Modern observational astrophysics is rich in findings of compact objects and as such our understanding of highly dense gravitational objects is rapidly advancing. However, there is as of yet no strong consensus on their underlying physics.A number of such objects have been recently observed that are inconsistent with standard GR and a simple neutron star equation of state. It was recently shown zhang_2021_unified ; zhang_2021_stellar that in standard general relativity, the secondary component of the merger GW190814 could feasibly be a quark star with an interacting equation of state governed by a single parameter $\lambda$ . This parameterization of strong interaction effects was inspired by another recent theoretical study showing that non-strange quark matter (QM) could feasibly be the ground state of baryonic matter at sufficient density and temperature holdom_2018_quark . Similar analyses with a different equation of state and/or QM phase miao_2021_bayesian ; lopes_2022_onthe ; oikonomou_2023_colourflavour found similarly promising results.This same object was subsequently shown to be well described as a slowly-rotating neutron star in the 4DEGB theory without resorting to exotic quark matter EOSs, while also demonstrating that the equilibrium sequence of neutron stars asymptotically matches the black hole limit, thus closing the mass gap between NS/black holes of the same radius charmousis2022 . More recently, some groups have also been interested in modelling ECOs (extreme compact objects) zhang_2023_rescaling ; gammon_2024 ; Mann:2021mnc ; Conklin:2021cbc as well as unusually light compact stars horvath_2023_alight ; oikonomou_2023_colourflavour (like that in the gamma ray remnant J1731-347, which is inconsistent with minimum mass calculations of neutron stars generated by iron cores) as quark stars to explain their unusual properties.

To further illuminate the range of possibilities, we consider in this paper charged quark star solutions to the 4DEGB theory, as the confining nature of the strong interaction make quark stars one of the few remaining candidates for charged stellar objects. Although charged quark star solutions have previously been considered in the context of 4DEGB pretel_2022 , the upper limit of the coupling $\alpha$ was taken to be much smaller than that allowed by current observational constraints Fernandes:2022zrq ; Clifton:2020xhc , and the relation with the charged 4DEGB Buchdahl bound was not discussed (as the bound was not derived). In considering values of $\alpha$ closer to presently allowed bounds, we obtain a number of interesting novel results.

Our most intriguing result (extending similar results found in gammon_2024 ) is that charged quark stars in 4DEGB can be Extreme Compact Charged Objects (ECCOs), objects whose radii are smaller than that allowed by the Buchdahl bound and Schwarzschild radius in GR. Indeed, there exist uncharged quark stars in 4DEGB whose radii are smaller than that of a corresponding black hole of the same mass in GR¹¹1This phenomenon was also shown to be present for neutron starscharmousis2022 , though the relationship with the Buchdahl bound was not noted.. Observations of these latter objects, apart from indicating a new class of astrophysical phenomena Mann:2021mnc , would provide strong evidence for 4DEGB as a physical theory. These ECOs respect a generalization of the Buchdahl bound, whose small-radius limit is that of the horizon radius of the corresponding minimum mass black hole. ECCOs likewise respect a generalization of the charged Buchdahl bound as we shall demonstrate.

We also find in general that for a given central pressure, charged quark stars in 4DEGB have a larger mass and radius than their GR counterparts with the same pressure. This can be attributed to the ‘less attractive’ nature of gravity in 4DEGB: for $\alpha>0$ the gradient of the effective potential yields a weaker force than pure GR. We consequently find a larger maximal mass for a given value of $Q$ in 4DEGB than in general relativity. Similarly, increasing $Q$ also weakens the overall force gradient compared to the uncharged case, and leads to larger quark stars relative to their uncharged counterparts.

The outline of our paper is as follows: In section II we introduce the basic theory underlying 4DEGB gravity, as well as the non-interacting quark matter equation of state that we make use of. Following this, a electromagnetic perfect fluid stress-energy tensor is employed to derive the charged 4DEGB TOV equations, and current observational constraints on the coupling constant are briefly discussed. We use section III to derive a generalization of Buchdahl’s bound for charged stars in 4DEGB gravity to be plotted alongside the solutions to the modified TOV equations. Section IV outlines the results of our numerical calculations. We conclude our results with a brief analysis of the stability of charged 4DEGB quark stars. Section VI summarizes our key findings and suggests topics for future study.

II Theory

II.1 4D Einstein-Gauss-Bonnet Gravity

To model the structure of a quark star in 4DEGB gravity, we add tothe action (4)a perfect-fluid term corresponding to the energy of the quark matter and

S_{\mathrm{Max}}=-\frac{1}{8\pi G}\int\sqrt{-g}\left({\frac{1}{4}}F_{\mu\nu}F^%{\mu\nu}+A_{\mu}J^{\mu}\right)

(5)

where $J^{\mu}$ is the electromagnetic current of the charged quark matter.The field equations of 4DEGB are obtained from a straightforward variational principle applied to this action.

Variation with respect to the scalar $\phi$ yields

$\displaystyle\mathcal{E}_{\phi}=$	$\displaystyle-\mathcal{G}+8G^{\mu\nu}\nabla_{\nu}\nabla_{\mu}\phi+8R^{\mu\nu}%\nabla_{\mu}\phi\nabla_{\nu}\phi-8(\square\phi)^{2}+8(\nabla\phi)^{2}\square%\phi+16\nabla^{a}\phi\nabla^{\nu}\phi\nabla_{\nu}\nabla_{\mu}\phi$	(6)
	$\displaystyle\qquad+8\nabla_{\nu}\nabla_{\mu}\phi\nabla^{\nu}\nabla^{\mu}\phi$
$\displaystyle=$	$\displaystyle\;0$

while the variation with respect to the metric gives

$\displaystyle\mathcal{E}_{\mu\nu}$	$\displaystyle=\Lambda g_{\mu\nu}+G_{\mu\nu}+\alpha\left[\phi H_{\mu\nu}-2R%\left[\left(\nabla_{\mu}\phi\right)\left(\nabla_{\nu}\phi\right)+\nabla_{\nu}%\nabla_{\mu}\phi\right]+8R_{(\mu}^{\sigma}\nabla_{\nu)}\nabla_{\sigma}\phi+8R_%{(\mu}^{\sigma}\left(\nabla_{\nu)}\phi\right)\left(\nabla_{\sigma}\phi\right)\right.$	(7)
	$\displaystyle-2G_{\mu\nu}\left[(\nabla\phi)^{2}+2\square\phi\right]-4\left[%\left(\nabla_{\mu}\phi\right)\left(\nabla_{\nu}\phi\right)+\nabla_{\nu}\nabla_%{\mu}\phi\right]\square\phi-\left[g_{\mu\nu}(\nabla\phi)^{2}-4\left(\nabla_{%\mu}\phi\right)\left(\nabla_{\nu}\phi\right)\right](\nabla\phi)^{2}$
	$\displaystyle+8\left(\nabla_{(\mu}\phi\right)\left(\nabla_{\nu)}\nabla_{\sigma%}\phi\right)\nabla^{\sigma}\phi-4g_{\mu\nu}R^{\sigma\rho}\left[\nabla_{\sigma}%\nabla_{\rho}\phi+\left(\nabla_{\sigma}\phi\right)\left(\nabla_{\rho}\phi%\right)\right]+2g_{\mu\nu}(\square\phi)^{2}$
	$\displaystyle-4g_{\mu\nu}\left(\nabla^{\sigma}\phi\right)\left(\nabla^{\rho}%\phi\right)\left(\nabla_{\sigma}\nabla_{\rho}\phi\right)+4\left(\nabla_{\sigma%}\nabla_{\nu}\phi\right)\left(\nabla^{\sigma}\nabla_{\mu}\phi\right)$
	$\displaystyle\left.-2g_{\mu\nu}\left(\nabla_{\sigma}\nabla_{\rho}\phi\right)%\left(\nabla^{\sigma}\nabla^{\rho}\phi\right)+4R_{\mu\nu\sigma\rho}\left[\left%(\nabla^{\sigma}\phi\right)\left(\nabla^{\rho}\phi\right)+\nabla^{\rho}\nabla^%{\sigma}\phi\right]\right]$
	$\displaystyle=\;T_{\mu\nu}$

where $T_{\mu\nu}$ is the stress-energy tensor of the charged quark matter and the electromagnetic field, and

\displaystyle H_{\mu\nu}=2\Big{[}RR_{\mu\nu}-2R_{\mu\alpha\nu\beta}R^{\alpha%\beta}+R_{\mu\alpha\beta\sigma}R_{\nu}^{\alpha\beta\sigma}-2R_{\mu\alpha}R_{%\nu}^{\alpha}-\frac{1}{4}g_{\mu\nu}\mathcal{G}\Big{]}

(8)

is the Gauss-Bonnet tensor. These field equations satisfy the following relationship

g^{\mu\nu}T_{\mu\nu}=g^{\mu\nu}\mathcal{E}_{\mu\nu}+\frac{\alpha}{2}\mathcal{E%}_{\phi}=4\Lambda-R-\frac{\alpha}{2}\mathcal{G}

(9)

which can act as a useful consistency check to see whether prior solutions generatedvia the solution-limit method glavan2020 are even possible solutions to the theory.For example, using (9)it is easy to verify that the rotating metrics generated from a Newman-Janis algorithm kumar2020 ; wei2020 are not solutions to the field equations of the scalar-tensor 4DEGB theory.

II.2 Charged 4DEGB TOV Equations

Throughout the following we assume a non-interacting quark equation of state

P(r)=\frac{1}{3}(\rho-4B_{\mathrm{eff}})

(10)

where $B_{\mathrm{eff}}$ is the MIT bag constant (for which we use a benchmark value of $60\;\mathrm{MeV}/\mathrm{fm}^{3}$ zhang_2021_stellar ).

The standard Tolman-Oppenheimer-Volkoff (TOV) equations for stellar structure are well-known in GR. Variation of the action with respect to the gauge potential $A_{\mu}$ , together with (6) and(7) yield the TOV equations for the charged 4DEGB quark star. To obtain these we begin with a static, spherically symmetric metric ansatz in natural units ( $G=c=1$ ):

ds^{2}=-e^{\Phi(r)}c^{2}dt^{2}+e^{\Lambda(r)}dr^{2}+r^{2}d\Omega^{2}.

(11)

As usual hennigar_2020_on ; gammon_2024 , so long as $e^{\Phi(r)}=e^{-\Lambda(r)}$ outside the star, the combination $\mathcal{E}_{0}^{0}-\mathcal{E}_{1}^{1}$ of the field equations can be used to derive the following equation for the scalar field:

II.3 Observational Constraints on the 4DEGB Coupling Constant

An investigation of the observational constraints on the coupling $\alpha$ yieldedClifton:2020xhc ; Fernandes:2022zrq

-10^{-30}\;\textrm{m}^{2}<\alpha<10^{10}\;\textrm{m}^{2}

(32)

where the lower bound comes from “early universe cosmology and atomic nuclei”data charmousis2022 ,and the upper bound follows from LAGEOS satellite observations. Regarding the lower bound as negligibly close to zero,the dimensionless version of (32) reads

0<\bar{\alpha}\lesssim 3.2.

(33)

We note that inclusion of preliminary calculations on recent GW data suggest these constraints could potentially tighten to $0<\alpha\lesssim 10^{7}\;\textrm{m}^{2}$ , or alternatively $0<\bar{\alpha}\lesssim 0.0032$ . This would mean that deviations from GR due to 4DEGB would only be detectable in extreme environments such as in the very early universe or near the surface of extremely massive objects.Even tighter bounds were assumed in previous studies of quark stars, where only solutions with $\alpha$ below 6 $\mathrm{km}^{2}$ ( $\bar{\alpha}\leq 0.0019$ ) were considered banerjee_2021_strange ; banerjee_2021_quark ; pretel_2022 .Adopting such a tight bound would make compact stars near the upper end of the mass gap an ideal candidate for investigation the effects of 4DEGB theory.

At this point in time such tighter bounds are not warranted. A proper study of the effects of gravitational radiation in 4DEGB has yet to be carried out. In view of this we shall assume the bound (33), which has strong observational support Clifton:2020xhc ; Fernandes:2022zrq .

III Generalization of the Buchdahl Bound for Charged Compact Stars in 4DEGB Gravity

III.1 Derivation of a Maximal Mass for Charged Compact Objects in 4DEGB Gravity

A generalization of the Buchdahl bound has been derived for charged stars in GR bohmer_2007 and for uncharged objects in 4DEGB theory chakraborty_2020 . In the following we extend these calculations for the case of a charged star in 4DEGB gravity.

By virtue of equation (18) it is straightforward to check that

\frac{d}{dr}\frac{(1-e^{-\Lambda(r)})}{2r^{2}}=\frac{1}{\sqrt{1+4\alpha(\frac{%2m(r)}{r^{3}}-\frac{q(r)^{2}}{r^{4}})}}\frac{d}{dr}\left(\frac{m(r)}{r^{3}}-%\frac{q(r)^{2}}{2r^{4}}\right)

(34)

The next step of the argument hinges on the assumption that the quantity $\frac{(1-e^{-\Lambda(r)})}{r^{2}}$ decreases as radial distance increases chakraborty_2020 . In the uncharged case this is equivalent to assuming decreasing mass density as the star’s surface is approached (though we note that some exotic star models do not have monotonically decreasing density profiles zhang_2021_stellar ).This assumption also seems reasonable in the case of a charged sphere chakraborty_2020 and we employ it in what follows:

Utilizing the conservation equation

(\rho(r)+P(r))=-\frac{2}{\Phi^{\prime}(r)}\left(\frac{\mathrm{d}P}{\mathrm{~{}%d}r}-\frac{q}{4\pi r^{4}}\frac{\mathrm{d}q}{\mathrm{~{}d}r}\right)

(35)

we can write

\frac{dP}{dr}=\frac{q}{4\pi r^{4}}\frac{dq}{dr}-\frac{1}{2}(\rho(r)+P(r))\Phi^%{\prime}(r).

(36)

Recalling the $rr$ component of the field equations, we obtain

		$\displaystyle\frac{1}{r^{2}}\left[r\Phi^{\prime}(r)e^{-\Lambda(r)}-\left(1-e^{%-\Lambda(r)}\right)\right]$		(37)
		$\displaystyle\qquad+\alpha\frac{\left(1-e^{-\Lambda(r)}\right)}{r^{4}}\left[2r%\Phi^{\prime}(r)e^{-\Lambda(r)}+\left(1-e^{-\Lambda(r)}\right)\right]=\kappa%\left(P(r)-\frac{q(r)^{2}}{8\pi r^{4}}\right)$		(37)

and subsequently get

\displaystyle\kappa\left(\frac{dP}{dr}+\frac{q(r)^{2}}{2\pi r^{5}}-\frac{q(r)}%{4\pi r^{4}}\frac{dq}{dr}\right)=\kappa\left(\frac{q(r)^{2}}{2\pi r^{5}}-\frac%{1}{2}(\rho(r)+P(r))\Phi^{\prime}(r)\right)

(38)

by differentiating the left hand side of (37) using(36).

Addition of the $tt$ and $rr$ equations yields

(\rho(r)+P(r))=\frac{1}{\kappa}\left(\frac{1}{r}e^{-\Lambda(r)}\left(\Lambda(r%)^{\prime}+\Phi(r)^{\prime}\right)+\frac{2\alpha}{r^{3}}e^{-\Lambda(r)}\left(1%-e^{-\Lambda(r)}\right)\left(\Lambda(r)^{\prime}+\Phi(r)^{\prime}\right)\right)

(39)

after division by $r^{2}$ . This can be substituted into (38) to obtain

\displaystyle\kappa\left(\frac{dP}{dr}+\frac{q(r)^{2}}{2\pi r^{5}}-\frac{q(r)}%{4\pi r^{4}}\frac{dq}{dr}\right)

\displaystyle=\kappa\frac{q(r)^{2}}{2\pi r^{5}}-\frac{1}{2}\left(\frac{1}{r}e^%{-\Lambda(r)}+\frac{2\alpha}{r^{3}}e^{-\Lambda(r)}\left(1-e^{-\Lambda(r)}%\right)\right)\left(\Lambda(r)^{\prime}+\Phi(r)^{\prime}\right)\Phi^{\prime}(r).

(40)

Making use of equation (37) we find

		$\displaystyle\frac{d}{dr}\left(\frac{1}{r^{2}}\left[r\Phi^{\prime}(r)e^{-%\Lambda(r)}-\left(1-e^{-\Lambda(r)}\right)\right]+\alpha\frac{\left(1-e^{-%\Lambda(r)}\right)}{r^{4}}\left[2r\Phi^{\prime}(r)e^{-\Lambda(r)}+\left(1-e^{-%\Lambda(r)}\right)\right]\right)$		(41)
		$\displaystyle=\kappa\frac{q(r)^{2}}{2\pi r^{5}}-\frac{1}{2}\left(\frac{1}{r}e^%{-\Lambda(r)}+\frac{2\alpha}{r^{3}}e^{-\Lambda(r)}\left(1-e^{-\Lambda(r)}%\right)\right)\left(\Lambda(r)^{\prime}+\Phi(r)^{\prime}\right)\Phi^{\prime}(r).$		(41)

If we define $\beta(r)=\frac{1-e^{-\Lambda(r)}}{r^{2}}$ the above can be rewritten as

\displaystyle\frac{d}{dr}\left(\frac{\Phi^{\prime}(r)e^{-\Lambda(r)}}{r}\right%)(1+2\alpha\beta(r))+Z(r)=\kappa\frac{q(r)^{2}}{2\pi r^{5}}-\frac{\left(\Phi^{%\prime}(r)+\Lambda^{\prime}(r)\right)}{2}\frac{\Phi^{\prime}(r)e^{-\Lambda(r)}%}{r}(1+2\alpha\beta(r))

(42)

where $Z(r)=\left(\frac{2\alpha\Phi^{\prime}(r)e^{-\Lambda(r)}}{r}-1+2\alpha\beta(r)%\right)\beta^{\prime}(r)$ .Equivalently,

	$\displaystyle 2e^{-(\Phi(r)+\Lambda(r))/2}\frac{d}{dr}\left[\frac{e^{-\Lambda(%r)/2}}{r}\frac{de^{\Phi(r)/2}}{dr}\right](1+2\alpha\beta(r))-\kappa\frac{q(r)^%{2}}{2\pi r^{5}}$
	$\displaystyle\qquad\qquad=\left(1-\frac{2\alpha\Phi^{\prime}(r)e^{-\Lambda(r)}%}{r}-2\alpha\beta(r)\right)\beta^{\prime}(r)$		(43)

which can be shown usingthe identity

\left(\frac{d}{dr}\left(\frac{\Phi^{\prime}(r)e^{-\Lambda(r)}}{r}\right)+\frac%{1}{2}(\frac{\Phi^{\prime}(r)e^{-\Lambda(r)}}{r})(\Phi^{\prime}(r)+\Lambda^{%\prime}(r))\right)=2e^{-(\Phi(r)+\Lambda(r))/2}\frac{d}{dr}\left[\frac{e^{-%\Lambda(r)/2}}{r}\frac{de^{\Phi(r)/2}}{dr}\right].

(44)

To proceed further we assume the factors $(1+2\alpha\beta(r))$ and $\left(1-\frac{2\alpha\Phi^{\prime}(r)e^{-\Lambda(r)}}{r}-2\alpha\beta(r)\right)$ are both non-negative, which is valid forsufficiently small, positive $\alpha$ . The right hand side of (43) is then negative due to the monotonically decreasing nature of $\beta(r)$ . We can thus write

2e^{-(\Phi(r)+\Lambda(r))/2}\frac{d}{dr}\left[\frac{e^{-\Lambda(r)/2}}{r}\frac%{de^{\Phi(r)/2}}{dr}\right](1+2\alpha\beta(r))-\kappa\frac{q(r)^{2}}{2\pi r^{5%}}\leq 0

(45)

which upon integration yields

e^{\Phi(r^{\prime\prime})/2}\leq\kappa\int^{r^{\prime\prime}}dr^{\prime}r^{%\prime}e^{\Lambda(r^{\prime})/2}\int^{r^{\prime}}dre^{(\Phi(r)+\Lambda(r))/2}%\frac{q(r)^{2}}{4\pi r^{5}}(1+2\alpha\beta(r))^{-1}\;.

(46)

Making for convenience the following definitions

$\displaystyle\mathcal{Q}$	$\displaystyle\equiv e^{(\Phi(r)+\Lambda(r))/2}\frac{q(r)^{2}}{4\pi r^{5}}(1+2%\alpha\beta(r))^{-1}$	(47)
$\displaystyle\eta(r^{\prime\prime})$	$\displaystyle\equiv\int^{r^{\prime\prime}}e^{\Lambda(r^{\prime})/2}r^{\prime}%\int^{r^{\prime}}\mathcal{Q}(r)drdr^{\prime}$
$\displaystyle\psi$	$\displaystyle\equiv e^{\Phi(r^{\prime\prime})/2}-\eta(r^{\prime\prime})\qquad%\xi\equiv\int^{r^{\prime\prime}}r^{\prime}e^{\Lambda(r^{\prime})/2}dr^{\prime}$

it is straightforward to show that

	$\displaystyle\frac{d\psi}{d\xi}=\frac{1}{r^{\prime\prime}}e^{-\Lambda(r^{%\prime\prime})/2}\frac{d}{dr^{\prime\prime}}e^{\Phi(r^{\prime\prime})/2}-\int^%{r^{\prime\prime}}dr\mathcal{Q}(r^{\prime}),$		(48)
	$\displaystyle\frac{d^{2}\psi}{d\xi^{2}}=\frac{1}{r^{\prime\prime}}e^{-\Lambda(%r^{\prime\prime})/2}\frac{d}{dr}\left[\frac{1}{r^{\prime\prime}}e^{-\Lambda(r^%{\prime\prime})/2}\frac{d}{dr^{\prime\prime}}e^{\Phi(r^{\prime\prime})/2}%\right]-\frac{1}{r^{\prime\prime}}e^{-\Lambda(r^{\prime\prime})/2}\mathcal{Q}(%r^{\prime\prime})\leq 0$		(49)

where integration is from the center ( $r^{\prime\prime}>r^{\prime}>r$ ),or alternatively

\frac{d^{2}\psi}{d\xi^{2}}re^{\Lambda(r)/2}\leq 0

(50)

which follows from (45).Since $e^{\Lambda(r)/2}>0$ and $r>0$ this is equivalent to

\frac{d^{2}\psi}{d\xi^{2}}\leq 0.

(51)

Trivially integrating and applying the mean value theorem bohmer_2007 we see that

\frac{d\psi}{d\xi}\leq\frac{\psi(\xi)-\psi(0)}{\xi-0}\leq\frac{\psi(\xi)}{\xi}

(52)

since at the center of the star $\xi=0$ and $\psi(0)>0$ since $e^{\Phi(r)/2}>0$ inside the star.

Putting everything together we find

		$\displaystyle\left(\frac{1}{r}e^{-\Lambda(r^{\prime\prime})/2}\frac{d}{dr}e^{%\Phi(r^{\prime\prime})/2}-\int^{r^{\prime\prime}}dr^{\prime}\mathcal{Q}(r^{%\prime})\right)\left(\int^{r^{\prime\prime}}dr^{\prime}r^{\prime}e^{\Lambda(r^%{\prime})/2}\right)$
		$\displaystyle\qquad\qquad\qquad\leq e^{\Phi(r^{\prime\prime})/2}-\int^{r^{%\prime\prime}}r^{\prime}e^{\Lambda_{\mathrm{int}}/2}\int^{r^{\prime}}\mathcal{%Q}(r)drdr^{\prime}.$		(53)

Recalling that $\beta(r)=\frac{1-e^{-\Lambda(r)}}{r^{2}}=-\frac{1}{2\alpha}\left[1-\sqrt{1+%\frac{\kappa\alpha}{2\pi}\left(\frac{2m(r)}{r^{3}}-\frac{q(r)^{2}}{r^{4}}%\right)}\right]$ such that

e^{-\Lambda(r)}=1-r^{2}\beta(r),

(54)

and that $\beta(r)\geq\beta(r^{\prime})\geq\beta(r^{\prime\prime})$ (more primes corresponding to further from the center), we obtain

	$\displaystyle\xi\equiv\int^{r^{\prime\prime}}dr^{\prime}r^{\prime}e^{\Lambda(r%^{\prime})/2}=\int^{r^{\prime\prime}}dr^{\prime}r^{\prime}\frac{1}{\sqrt{1-r^{%\prime 2}\beta(r^{\prime})}}$	$\displaystyle\geq\int^{r^{\prime\prime}}dr^{\prime}r^{\prime}\frac{1}{\sqrt{1-%r^{\prime 2}\beta(r^{\prime\prime})}}$		(55)
		$\displaystyle=\frac{1}{\beta(r^{\prime\prime})}\left(1-\sqrt{1-r^{\prime\prime2%}\beta(r^{\prime\prime})}\right)$		(55)

as a bound on $\xi$ .

Similarly we can consider the term

	$\displaystyle\int^{r^{\prime\prime}}\mathcal{Q}(r^{\prime})dr$	$\displaystyle=\kappa\int^{r^{\prime\prime}}e^{(\Phi(r^{\prime})+\Lambda(r^{%\prime}))/2}\frac{q(r^{\prime})^{2}}{4\pi r^{\prime 5}}(1+2\alpha\beta(r^{%\prime}))^{-1}dr$		(56)
		$\displaystyle=\frac{\kappa}{4\pi}\int^{r^{\prime\prime}}\frac{e^{\Phi(r^{%\prime})/2}}{\sqrt{1-r^{\prime 2}\beta(r^{\prime})}\left(1+2\alpha\beta(r^{%\prime})\right)}\frac{q(r^{\prime})^{2}}{r^{\prime 5}}dr^{\prime}.$		(56)

Combining our previous assumption on $\beta$ with the following assumption regarding the behaviour of the charge:

e^{\Phi_{\mathrm{int}}(r)/2}\frac{q(r)^{2}}{r^{5}}\geq e^{\Phi_{\mathrm{int}}(%r^{\prime})/2}\frac{q(r^{\prime})^{2}}{r^{\prime 5}}\geq e^{\Phi_{\mathrm{int}%}(r^{\prime\prime})/2}\frac{q(r^{\prime\prime})^{2}}{r^{\prime\prime 5}}

(57)

the above integral can be bounded as follows:

	$\displaystyle\frac{\kappa}{4\pi}\int^{r^{\prime\prime}}\frac{e^{\Phi(r^{\prime%})/2}}{\sqrt{1-r^{\prime 2}\beta(r^{\prime})}\left(1+2\alpha\beta(r^{\prime})%\right)}\frac{q(r^{\prime})^{2}}{r^{\prime 5}}dr^{\prime}$	$\displaystyle\geq\frac{\kappa}{4\pi}\frac{e^{\Phi_{\mathrm{int}}(r^{\prime%\prime})/2}}{\left(1+2\alpha\beta(0)\right)}\frac{q(r^{\prime\prime})^{2}}{r^{%\prime\prime 5}}\int^{r^{\prime\prime}}\frac{dr^{\prime}}{\sqrt{1-r^{\prime 2}%\beta(r^{\prime\prime})}}$		(58)
		$\displaystyle=\frac{\kappa}{4\pi}\frac{e^{\Phi_{\mathrm{int}}(r^{\prime\prime}%)/2}}{\left(1+2\alpha\beta(0)\right)}\frac{q(r^{\prime\prime})^{2}}{r^{\prime%\prime 5}}\frac{\mathrm{arcsin}\left(\sqrt{\beta(r^{\prime\prime})}r^{\prime%\prime}\right)}{\sqrt{\beta(r^{\prime\prime})}}.$		(58)

Finally we can consider the most complicated integral term:

		$\displaystyle\int^{r^{\prime\prime}}r^{\prime}e^{\Lambda_{\mathrm{int}}/2}\int%^{r^{\prime}}\mathcal{Q}(r)drdr^{\prime}=\frac{\kappa}{4\pi}\int^{r^{\prime%\prime}}r^{\prime}e^{\Lambda_{\mathrm{int}}/2}\int^{r^{\prime}}\frac{e^{\Phi_{%\mathrm{int}}/2}}{\sqrt{1-r^{2}\beta(r)}\left(1+2\alpha\beta(r)\right)}\frac{q%(r)^{2}}{r^{5}}drdr^{\prime}$		(59)
		$\displaystyle\geq\frac{\kappa}{4\pi}\left(1+2\alpha\beta(0)\right)^{-1}\int^{r%^{\prime\prime}}r^{\prime}e^{(\Phi_{\mathrm{int}}(r^{\prime})+\Lambda_{\mathrm%{int}}(r^{\prime}))/2}\frac{q(r^{\prime})^{2}}{r^{\prime 5}}\frac{\mathrm{%arcsin}\left(\sqrt{\beta(r^{\prime})}r^{\prime}\right)}{\sqrt{\beta(r^{\prime}%)}}dr^{\prime}$
		$\displaystyle\geq\frac{\kappa}{4\pi}\left(1+2\alpha\beta(0)\right)^{-1}e^{\Phi%_{\mathrm{int}}(r^{\prime\prime})/2}\frac{q(r^{\prime\prime})^{2}}{r^{\prime%\prime 5}}\int^{r^{\prime\prime}}r^{\prime 2}\frac{1}{\sqrt{1-r^{\prime 2}%\beta(r^{\prime})}}\frac{\mathrm{arcsin}\left(\sqrt{\beta(r^{\prime})}r^{%\prime}\right)}{\sqrt{\beta(r^{\prime})}r^{\prime}}dr^{\prime}.$

Since the function $\mathrm{arcsin}(x)/x$ monotonically increases on the interval $x\in[0,1]$ , the integral is minimal when the smallest $\beta$ is chosen. Thus

		$\displaystyle\int^{r^{\prime\prime}}r^{\prime}e^{\Lambda_{\mathrm{int}}/2}\int%^{r^{\prime}}\mathcal{Q}(r)drdr^{\prime}$		(60)
		$\displaystyle\geq\frac{\kappa}{4\pi}\left(1+2\alpha\beta(0)\right)^{-1}e^{\Phi%_{\mathrm{int}}(r^{\prime\prime})/2}\frac{q(r^{\prime\prime})^{2}}{r^{\prime%\prime 5}}\int^{r^{\prime\prime}}r^{\prime 2}\frac{1}{\sqrt{1-r^{\prime 2}%\beta(r^{\prime\prime})}}\frac{\mathrm{arcsin}\left(\sqrt{\beta(r^{\prime%\prime})}r^{\prime}\right)}{\sqrt{\beta(r^{\prime\prime})}r^{\prime}}dr^{\prime}$
		$\displaystyle=\frac{\kappa}{4\pi}\left(1+2\alpha\beta(0)\right)^{-1}e^{\Phi_{%\mathrm{int}}(r^{\prime\prime})/2}\frac{q(r^{\prime\prime})^{2}}{r^{\prime%\prime 5}}\frac{1}{\beta(r^{\prime\prime})}\left(r^{\prime\prime}-\sqrt{\frac{%1-\beta(r^{\prime\prime})r^{\prime\prime 2}}{\beta(r^{\prime\prime})}}\mathrm{%arcsin}(\sqrt{\beta(r^{\prime\prime})}r^{\prime\prime})\right).$

Inserting these results into (III.1)yields

		$\displaystyle\left(e^{(\Phi(r^{\prime\prime})-\Lambda(r^{\prime\prime}))/2}%\frac{\Phi^{\prime}(r^{\prime\prime})}{2r}-\frac{\kappa}{4\pi}\frac{e^{\Phi(r^%{\prime\prime})/2}}{\left(1+2\alpha\beta(0)\right)}\frac{q(r^{\prime\prime})^{%2}}{r^{\prime\prime 5}}\frac{\mathrm{arcsin}\left(\sqrt{\beta(r^{\prime\prime}%)}r^{\prime\prime}\right)}{\sqrt{\beta(r^{\prime\prime})}}\right)\left(\frac{%\left(1-\sqrt{1-r^{\prime\prime 2}\beta(r^{\prime\prime})}\right)}{\beta(r^{%\prime\prime})}\right)$		(61)
		$\displaystyle\qquad\leq e^{\Phi(r^{\prime\prime})/2}-\frac{\kappa}{4\pi}\left(%1+2\alpha\beta(0)\right)^{-1}e^{\Phi(r^{\prime\prime})/2}\frac{q(r^{\prime%\prime})^{2}}{r^{\prime\prime 5}}\frac{1}{\beta(r^{\prime\prime})}\left(r^{%\prime\prime}-\sqrt{\frac{1-\beta(r^{\prime\prime})r^{\prime\prime 2}}{\beta(r%^{\prime\prime})}}\mathrm{arcsin}(\sqrt{\beta(r^{\prime\prime})}r^{\prime%\prime})\right)$		(61)

and cancelling a factor of $e^{\Phi(r^{\prime\prime})/2}$ leaves

		$\displaystyle\left(\frac{\sqrt{1-r^{\prime\prime 2}\beta(r^{\prime\prime})}}{r%^{\prime\prime}}\frac{\Phi^{\prime}(r^{\prime\prime})}{2}-\frac{\kappa}{4\pi}%\frac{1}{\left(1+2\alpha\beta(0)\right)}\frac{q(r^{\prime\prime})^{2}}{r^{%\prime\prime 5}}\frac{\mathrm{arcsin}\left(\sqrt{\beta(r^{\prime\prime})}r^{%\prime\prime}\right)}{\sqrt{\beta(r^{\prime\prime})}}\right)\left(\frac{\left(%1-\sqrt{1-r^{\prime\prime 2}\beta(r^{\prime\prime})}\right)}{\beta(r^{\prime%\prime})}\right)$		(62)
		$\displaystyle\qquad\qquad\leq 1-\frac{\kappa}{4\pi}\left(1+2\alpha\beta(0)%\right)^{-1}\frac{q(r^{\prime\prime})^{2}}{r^{\prime\prime 5}}\frac{1}{\beta(r%^{\prime\prime})}\left(r^{\prime\prime}-\sqrt{\frac{1-\beta(r^{\prime\prime})r%^{\prime\prime 2}}{\beta(r^{\prime\prime})}}\mathrm{arcsin}(\sqrt{\beta(r^{%\prime\prime})}r^{\prime\prime})\right)$		(62)

which depends explicitly on $\Phi^{\prime}(r^{\prime\prime})$ . This can be determined from the $rr$ field equation (16) to give

\Phi^{\prime}(r)=\frac{\alpha+e^{2\Lambda(r)}\left(\alpha-8\pi r^{4}P(r)+q(r)^%{2}-r^{2}\right)+e^{\Lambda(r)}\left(r^{2}-2\alpha\right)}{2\alpha r-re^{%\Lambda(r)}\left(2\alpha+r^{2}\right)}.

(63)

Making this substitution and simplifying, we find the desired inequality

		$\displaystyle\frac{\left(\sqrt{1-R^{2}\beta(R)}-1\right)\left(q(R)^{2}+R^{4}%\beta(R)(-1+\alpha\beta(R))\right)}{2R^{4}\beta(R)\sqrt{1-R^{2}\beta(R)}(1+2%\alpha\beta(R))}$		(64)
		$\displaystyle\qquad\qquad\qquad\leq 1+\frac{2q(r)^{2}}{R^{4}\beta(R)}\left(%\frac{\mathrm{arcsin}(\sqrt{\beta(R)}R)}{\sqrt{\beta(R)}R}-1\right)$		(64)

which is the generalized Buchdahl bound for charged spheres in 4DEGB gravity.

In the $q\to 0$ limit we recover the following inequality chakraborty_2020

\sqrt{1-\beta(R)R^{2}}\left(1+\alpha\beta(R)\right)>\frac{1}{3}\left(1-\alpha%\beta(R)\right)

(65)

for uncharged objects in 4DEGB.

Alternatively, in the small $\alpha$ limit $\beta_{\mathrm{\alpha\to 0}}(r)=\frac{2m(r)}{r^{3}}-\frac{q(r)^{2}}{r^{4}}$ , and consequently

\frac{\left(\sqrt{1-r^{2}\beta(r)}-1\right)\left(q(r)^{2}-rm(r)\right)}{r^{4}%\sqrt{1-r^{2}\beta(r)}}\leq\beta(r)+\frac{2q(r)^{2}}{r^{4}}\left(\frac{\sin^{-%1}\left(\sqrt{r^{2}\beta(r)}\right)}{\sqrt{r^{2}\beta(r)}}-1\right)

(66)

which matches²²2At the time of publication bohmer_2007 has a factor of 2 missing from their final expression which Dr. Harko has kindly confirmed to be a typo.We note that in bohmer_2007 $\alpha(r)=\beta_{\mathrm{\alpha\to 0}}(r)\frac{r^{3}}{2m(r)}$ in our notation. the final inequality from bohmer_2007 .

III.2 Relation to the Black Hole Horizon

For a given $\alpha$ the radius of a black hole is $R=M+\sqrt{M^{2}-Q^{2}-\alpha}$ . Since the minimum mass $M_{\mathrm{min}}=\sqrt{Q^{2}+\alpha}$ , we obtain $R_{\mathrm{Mmin}}=\sqrt{Q^{2}+\alpha}$ .Substituting this into the Buchdahl bound(64)we find that it is automatically satisfied as an equality; in other words the Buchdahl bound intersects the minimum mass point of the black hole horizon. The $M/R$ curves for nonzero $\alpha$ also smoothly join this point in the limit of large central pressure, as can be seen in figures 1 and 2.

IV Results

Charged Quark Stars and Extreme Compact Objects in Regularized 4D Einstein-Gauss-Bonnet Gravity (1)

Charged Quark Stars and Extreme Compact Objects in Regularized 4D Einstein-Gauss-Bonnet Gravity (2)

Charged Quark Stars and Extreme Compact Objects in Regularized 4D Einstein-Gauss-Bonnet Gravity (7)

Charged Quark Stars and Extreme Compact Objects in Regularized 4D Einstein-Gauss-Bonnet Gravity (8)

The mass/radius ( $M/R$ ) and mass/central density ( $M/\rho_{c}$ ) curves for charged 4DEGB quark stars are presented in figures 1 and 2. As expected from previous workgammon_2024 ; zhang_2021_stellar ,we observe the general trends that for a larger $\alpha$ and/or $Q$ , the mass/radius profiles of the quark stars increase in size. For a fixed nonzero charge the stars have a minimum size,below which the gravitational attraction cannot overcome the self-repulsion of the charge. This is most evident in the right-hand lower panels offigure2.As $\alpha$ gets large the $M/R$ curves converge for all values of fixed $Q$ (and similarly with the $M/\rho_{c}$ curves) with a slight divergence near $\bar{\rho_{c}}=1$ in the $M/\rho_{c}$ plots to account for the structures mentioned above.

As in gammon_2024 , even for small, nonzero $\alpha$ we find curves for the Buchdahl bound that intersect the black hole horizon at the minimum mass point. A new feature, evidently not previously noted, is that this same interaction between the black hole horizon and Buchdahl bound can be observed in GR when $Q\neq 0$ . This is evident in the orange and green dashed curves in the upper left panel of figure1.The difference is that for nonzero $\alpha$ the $M/R$ curves approach this intersection point smoothly, whereas in GR the curves turn away from the horizon/BB and thus never meet.

On the solution curves we have indicated maximum mass points (when this point doesn’t occur at the intersection with the horizon) with a black dot - in general relativity these maximum mass points indicate a transition point from stability against radial perturbations to instability for uncharged stars. In 4DEGB it’s not clear whether this coincidence holds, and offers an interesting avenue for future research.

In the uncharged case a criticality condition was previously noted gammon_2024 wherein for $\bar{\alpha}\geq\frac{3}{4\pi}$ a critical central pressure was present - for central pressures below this critical value, the pressure function diverges with no real roots to define a star’s surface. Above this critical pressure we find numerical solutions which lie at or beyond the black hole horizon, indicating a lack of stable physical stars in this critical regime. Deriving an analogous bound for the charged case is considerably more complicated, however numerically we find that for $Q>0$ (and hence $\gamma>0$ ), such criticality can take place for $\bar{\alpha}<\frac{3}{4\pi}$ . Because of this lack of an analytic bound for the charged case,we do not present numerical solutions in which such behaviour is manifest. With this, we note that we have considered values of $\bar{\alpha}\leq 0.1$ . For larger values of $\bar{\alpha}$ ‘critical central pressure’ behaviour emerges soon after for the charged cases. Examination of these cases requires a more detailed stability analysis, which we leave for future investigation.

Most interestingly we note, similar to the uncharged casegammon_2024 , that the stars described by this theory can be ECCOs (Extreme Compact Charged Objects): charged stellar objects that are smaller (more compact) than the uncharged Buchdahl bound (and sometimes even the Schwarzschild radius) of general relativity. For large enough $\alpha$ the associated $M/R$ curves do not reach their maximum mass until they intersect the black hole horizon - these are particularly interesting candidates for stable ECCOs since in pure uncharged GR we only see stability against radial perturbations when $dM/d\rho>0$ . However, a net charge offsets the stability point from this maximum mass point zhang_2021_stellar - it is likely that the modifications to gravity will have a similar offsetting effect which should be investigated in future studies before declaring these solutions stable against perturbations.

V Stability

We now briefly consider the stability of stars respecting the generalized Buchdahl bound. In general relativity a necessary but insufficient condition for an uncharged compact star to be stable against radial perturbations is $dM/d\rho_{c}>0$ zhang_2021_stellar ; glendenningbook ; arbanil2015 , corresponding to the part of the solution curve before a maximal mass point is reached. In Einstein-Maxwell theory a net charge has been shown to offset the stability point from the maximum mass point in either direction zhang_2021_stellar . In a similar vein, when the coupling to higher curvature gravity is nonzero, it is not obvious whether this coincidence of stability and the maximum mass point will hold for uncharged stars. While we leave a thorough analysis of the fundamental radial oscillation modes for future work (where stability against radial oscillations can be ensured), the speed of sound and effective adiabatic index inside the star can still be discussed.

In the interior of a stable star, the sound speed $c_{s}$ must never exceed the speed of light $c$ . The non-interacting quark equation of state (10) has a constant subluminal sound speed of $c_{s}=\frac{c}{\sqrt{3}}$ , and thus the causality condition is always satisfied. Similarly, the effective adiabatic index

\gamma_{\mathrm{eff}}\equiv\left(1+\frac{\rho}{P}\right)\left(\frac{dP}{d\rho}%\right)_{S}

(67)

is often used as another indicator of stability as it is seen as a bridge between “the relativistic structure of a spherical static object and the equation of state of the interior fluid” moustakidis2017_stability . The subscript $S$ in the above equation indicates that we consider the sound speed at a constant specific entropy. In principle a critical value for $\langle\gamma_{\mathrm{eff}}\rangle$ exists, below which configurations are unstable against radial perturbations. In standard general relativity this critical value can be written chandrasekhar1964 ; moustakidis2017_stability $\gamma_{cr}=\frac{4}{3}+\frac{19}{42}\beta$ , where $\beta=2M/R=R_{S}/R$ is the compactness parameter. If $\beta\to 0$ the well-known classical Newtonian limit is recovered as expected ( $\langle\gamma_{\mathrm{eff}}\rangle\geq\frac{4}{3}$ ).

An equivalent bound has not yet been derived for the 4DEGB theory. Despite this, it is common practice to plot the adiabatic index of the star relative to the Newtonian critical value banerjee_2021_quark ; banerjee_2021_strange ; hansraj_2020_isotropic ; singh_2022_anisotropic . Since we restrict ourselves to a simple, non-interacting equation of state (10) it is straightforward to show that $\gamma_{\mathrm{eff}}=\frac{1}{3}(4+\frac{1}{P(r)})$ and hence the Newtonian bound of 4/3 is always satisfied for non-negative pressure.

VI Summary

In this paper we have investigated the stellar structure of strongly interacting quark stars in the 4D Einstein Gauss-Bonnet theory of gravity for different combinations of the charge $Q$ and 4DEGB coupling constant $\alpha$ . In accord with the lack of a mass gap in the 4DEGB theory charmousis2022 , we find that even for small $\alpha$ the quark star solutions asymptotically approach the 4DEGB black hole horizon radius, and thus have solutions with smaller radii than the GR Buchdahl/Schwarzschild limits. In general, larger $Q$ and $\alpha$ tend to increase the mass-radius profile of quark stars, with large $\alpha$ suppressing the differences between different charges. These findings are generally consistent with what was found in the regime of weak coupling to the 4DEGB theory pretel_2022 .

We have found many additional features in the unexplored regions of parameter space, the most striking of which is that 4DEGB charged quark stars can exist with radiinot only smaller than the general relativistic Buchdahl bound, but also smaller than the $2M$ Schwarzschild radius. Such Extreme Compact Charged Objects represent a possible new state of matter.A full analysis of the stability of such objects would be an interesting topic for future study.

Acknowledgements

This work was supported in part by the Natural Sciences and Engineering Research Council of Canada. We would like to thank Dr. Tiberiu Harko for useful correspondence regarding the existing literature.

References

[1]P.Bueno &P. Cano.Einsteinian cubic gravity.Phys. Rev. D, 94, 2016.
[2]ThomasP. Sotiriou and Valerio Faraoni.f(r) theories of gravity.Reviews of Modern Physics, 82(1):451–497, March 2010.
[3]Shin’ichi Nojiri and SergeiD. Odintsov.Introduction to modified gravity and gravitational alternatives for dark energy.International Journal of Geometric Methods in Modern Physics, 04(01):115–145, 2007.
[4]Timothy Clifton, PedroG. Ferreira, Antonio Padilla, and Constantinos Skordis.Modified gravity and cosmology.Physics Reports, 513(1–3):1–189, March 2012.
[5]J.Ahmed etal.Quintessential Quartic Quasi-topological Quartet, volume 134.2017.
[6]K.S. Stelle.Renormalization of higher-derivative quantum gravity.Phys. Rev. D, 16:953–969, Aug 1977.
[7]RobertH. Brandenberger.A nonsingular universe, 1992.
[8]ROBERTH. BRANDENBERGER.Towards a nonsingular universea.Annals of the New York Academy of Sciences, 688(1):446–453, June 1993.
[9]RobertH. Brandenberger.Implementing markov’s limiting curvature hypothesis, 1995.
[10]D.Lovelock.The einstein tensor and its generalizations.Journal of Mathematical Physics, 12, 1971.
[11]RobieA. Hennigar, David Kubizňák, RobertB. Mann, and Christopher Pollack.On taking the d → 4 limit of gauss-bonnet gravity: theory and solutions.Journal of High Energy Physics, 2020, 07 2020.
[12]Pedro G.S. Fernandes, Pedro Carrilho, Timothy Clifton, and DavidJ. Mulryne.Derivation of Regularized Field Equations for the Einstein-Gauss-Bonnet Theory in Four Dimensions.Phys. Rev. D, 102(2):024025, 2020.
[13]Drazen Glavan and Chunshan Lin.Einstein-gauss-bonnet gravity in four-dimensional spacetime.Physical Review Letters, 124(8), Feb 2020.
[14]Shou-Long Li, Puxun Wu, and Hongwei Yu.Stability of the einstein static universe in $4d$ gauss-bonnet gravity, 2020.
[15]Tsutomu Kobayashi.Effective scalar-tensor description of regularized lovelock gravity in four dimensions.Journal of Cosmology and Astroparticle Physics, 2020(07):013–013, jul 2020.
[16]Arun Kumar, RahulKumar Walia, and SushantG. Ghosh.Bardeen black holes in the regularized 4d einstein gauss bonnet gravity.Universe, 8(4):232, apr 2022.
[17]PedroG.S. Fernandes.Charged black holes in AdS spaces in 4d einstein gauss-bonnet gravity.Physics Letters B, 805:135468, jun 2020.
[18]Rahul Kumar and SushantG. Ghosh.Rotating black holes in 4d einstein-gauss-bonnet gravity and its shadow.Journal of Cosmology and Astroparticle Physics, 2020(07):053–053, jul 2020.
[19]Arun Kumar, Dharmanand Baboolal, and SushantG. Ghosh.Nonsingular black holes in 4d einstein–gauss–bonnet gravity.Universe, 8(4):244, apr 2022.
[20]Daniele Malafarina, Bobir Toshmatov, and Naresh Dadhich.Dust collapse in 4d einstein–gauss–bonnet gravity.Physics of the Dark Universe, 30:100598, dec 2020.
[21]DanielaD. Doneva and StoytchoS. Yazadjiev.Relativistic stars in 4d einstein-gauss-bonnet gravity.Journal of Cosmology and Astroparticle Physics, 2021(05):024, may 2021.
[22]C.Charmousis, A.Lehébel, E.Smyrniotis, and N.Stergioulas.Astrophysical constraints on compact objects in 4d einstein-gauss-bonnet gravity.Journal of Cosmology and Astroparticle Physics, 2022(02):033, Feb 2022.
[23]SushantG. Ghosh and SunilD. Maharaj.Radiating black holes in the novel 4d einstein-gauss-bonnet gravity.2020.
[24]Metin Gürses, TahsinÇağrı Şişman, and Bayram Tekin.Is there a novel einstein–gauss–bonnet theory in four dimensions?The European Physical Journal C, 80(7), jul 2020.
[25]Wen-Yuan Ai.A note on the novel 4d einstein–gauss–bonnet gravity.Communications in Theoretical Physics, 72(9):095402, aug 2020.
[26]Fu-Wen Shu.Vacua in novel 4d einstein-gauss-bonnet gravity: Pathology and instability?Physics Letters B, 811:135907, dec 2020.
[27]RobertB. Mann and S.F. Ross.The $D\rightarrow 2$ limit of general relativity.Class. Quant. Grav., 10:1405–1408, 1993.
[28]H.Lu and YiPang.Horndeski gravity as $D\rightarrow 4$ limit of Gauss-Bonnet.Phys. Lett. B, 809:135717, 2020.
[29]GregoryWalter Horndeski.Second-order scalar-tensor field equations in a four-dimensional space.Int. J. Theor. Phys., 10:363–384, 1974.
[30]Timothy Clifton, Pedro Carrilho, Pedro G.S. Fernandes, and DavidJ. Mulryne.Observational Constraints on the Regularized 4D Einstein-Gauss-Bonnet Theory of Gravity.Phys. Rev. D, 102(8):084005, 2020.
[31]Carola M.A. Zanoletti, BraydenR. Hull, C.Danielle Leonard, and RobertB. Mann.Cosmological constraints on 4-dimensional Einstein-Gauss-Bonnet gravity.JCAP, 01:043, 2024.
[32]Pedro G.S. Fernandes, Pedro Carrilho, Timothy Clifton, and DavidJ. Mulryne.The 4D Einstein–Gauss–Bonnet theory of gravity: a review.Class. Quant. Grav., 39(6):063001, 2022.
[33]Chen Zhang and RobertB. Mann.Unified interacting quark matter and its astrophysical implications.Physical Review D, 103, 03 2021.
[34]Chen Zhang, Michael Gammon, and RobertB. Mann.Stellar structure and stability of charged interacting quark stars and their scaling behaviour.Physical Review D, 104, 12 2021.
[35]Bob Holdom, Jing Ren, and Chen Zhang.Quark matter may not be strange.Physical Review Letters, 120, 05 2018.
[36]Zhiqiang Miao, Jin-Liang Jiang, Ang Li, and Lie-Wen Chen.Bayesian inference of strange star equation of state using the GW170817 and GW190425 data.The Astrophysical Journal Letters, 917(2):L22, aug 2021.
[37]LuizL. Lopes and DeboraP. Menezes.On the nature of the mass-gap object in the GW190814 event.The Astrophysical Journal, 936(1):41, aug 2022.
[38]P.T. Oikonomou and Ch.C. Moustakidis.Colour-flavour locked quark stars in light of the compact object in hess j1731-347 and the gw190814 event, 2023.
[39]Chen Zhang, Yong Gao, Cheng-Jun Xia, and Renxin Xu.Rescaling strangeon stars and its implications on gravitational-wave echoes, 2023.
[40]Michael Gammon, Sarah Rourke, and RobertB. Mann.Quark stars with a unified interacting equation of state in regularized 4d einstein-gauss-bonnet gravity.Phys. Rev. D, 109:024026, Jan 2024.
[41]RobertB. Mann, Sebastian Murk, and DanielR. Terno.Black holes and their horizons in semiclassical and modified theories of gravity.Int. J. Mod. Phys. D, 31(09):2230015, 2022.
[42]RandyS. Conklin and Niayesh Afshordi.Boltzmann Meets Lorentz: A Surrogate Model for Black Hole Echoes.12 2021.
[43]J.E. Horvath, L.S. Rocha, L.M. deSá , P.H. R.S. Moraes, L.G. Barão, M.G.B. deAvellar, A.Bernardo, and R.R.A. Bachega.A light strange star in the remnant HESS j1731-347: Minimal consistency checks.Astronomy & Astrophysics, 672:L11, apr 2023.
[44]Juan M.Z. Pretel, Ayan Banerjee, and Anirudh Pradhan.Electrically charged quark stars in 4d einstein–gauss–bonnet gravity.The European Physical Journal C, 82(2), feb 2022.
[45]Shao-Wen Wei and Yu-Xiao Liu.Testing the nature of gauss-bonnet gravity by four-dimensional rotating black hole shadow, 2020.
[46]Subharthi Ray, AquinoL Espindola, Manuel Malheiro, JoséPS Lemos, and VilsonT Zanchin.Electrically charged compact stars and formation of charged black holes.Physical Review D, 68(8):084004, 2003.
[47]José D. V. Arbañil and M.Malheiro.Equilibrium and stability of charged strange quark stars.Physical Review D, 92(8), oct 2015.
[48]Ayan Banerjee, Takol Tangphati, and Phongpichit Channuie.Strange quark stars in 4d einstein–gauss–bonnet gravity.The Astrophysical Journal, 909(1):14, mar 2021.
[49]Ayan Banerjee, Takol Tangphati, Daris Samart, and Phongpichit Channuie.Quark stars in 4d einstein–gauss–bonnet gravity with an interacting quark equation of state.The Astrophysical Journal, 906:114, 01 2021.
[50]C.G. Böhmer and T.Harko.Minimum mass–radius ratio for charged gravitational objects.General Relativity and Gravitation, 39(6):757–775, mar 2007.
[51]Sumanta Chakraborty and Naresh Dadhich.Limits on stellar structures in lovelock theories of gravity.Physics of the Dark Universe, 30:100658, dec 2020.
[52]NormanK. Glendenning.Compact Stars.1996.
[53]Charalampos Moustakidis.The stability of relativistic stars and the role of the adiabatic index.General Relativity and Gravitation, 49, 04 2017.
[54]S.Chandrasekhar.The Dynamical Instability of Gaseous Masses Approaching the Schwarzschild Limit in General Relativity.Astrophys.J., 140:417, August 1964.
[55]Sudan Hansraj, Ayan Banerjee, Lushen Moodly, and MK Jasim.Isotropic compact stars in 4d einstein–gauss–bonnet gravity.Classical and Quantum Gravity, 38(3):035002, dec 2020.
[56]Ksh.Newton Singh, S.K. Maurya, Piyali Bhar, and Riju Nag.Anisotropic solution for polytropic stars in 4d einstein–gauss–bonnet gravity.The European Physical Journal C, 82(9), sep 2022.

		$\displaystyle\frac{e^{-2\Lambda(r)}\left(r^{2}e^{\Lambda(r)}\left(1-e^{\Lambda%(r)}\right)-r^{3}e^{\Lambda(r)}\Lambda^{\prime}(r)\right)}{r^{4}}$		(15)
		$\displaystyle+\alpha\frac{e^{-2\Lambda(r)}\left(-2re^{\Lambda(r)}\Lambda^{%\prime}(r)+2r\Lambda^{\prime}(r)-e^{\Lambda(r)}\left(1-e^{\Lambda(r)}\right)-e%^{\Lambda(r)}+1\right)}{r^{4}}=-\rho(r)-\frac{E(r)^{2}}{8\pi}$
		$\displaystyle\frac{e^{-2\Lambda(r)}\left(r^{3}e^{\Lambda(r)}\Phi^{\prime}(r)+r%^{2}e^{\Lambda(r)}\left(1-e^{\Lambda(r)}\right)\right)}{r^{4}}$
		$\displaystyle+\alpha\frac{e^{-2\Lambda(r)}\left(2re^{\Lambda(r)}\Phi^{\prime}(%r)-2r\Phi^{\prime}(r)-e^{\Lambda(r)}\left(1-e^{\Lambda(r)}\right)-e^{\Lambda(r%)}+1\right)}{r^{4}}=P(r)-\frac{E(r)^{2}}{8\pi}$		(16)
		$\displaystyle(\rho(r)+P(r))=-\frac{2}{\Phi^{\prime}(r)}\left(\frac{\mathrm{d}P%}{\mathrm{~{}d}r}-\frac{q}{4\pi r^{4}}\frac{\mathrm{d}q}{\mathrm{~{}d}r}\right)$		(17)

$\displaystyle\frac{dq}{dr}$	$\displaystyle=4\pi r^{2}\rho_{e}e^{\frac{\Lambda(r)}{2}}$	(19)
$\displaystyle\frac{dm}{dr}$	$\displaystyle=4\pi r^{2}\rho(r)+\frac{q(r)}{r}\frac{dq}{dr}$	(20)
$\displaystyle\frac{dP}{dr}$	$\displaystyle=(P(r)+\rho(r))\frac{\left(r^{3}(\Gamma+8\pi\alpha P(r)-1)-2%\alpha m(r)\right)}{\Gamma r^{2}\left((\Gamma-1)r^{2}-2\alpha\right)}+\frac{q(%r)}{4\pi r^{4}}\frac{dq}{dr}$	(21)

$\displaystyle\frac{d\bar{q}}{d\bar{r}}$	$\displaystyle=4\pi\bar{r}^{2}\bar{\rho}_{e}e^{\frac{\Lambda(\bar{r})}{2}}$	(25)
$\displaystyle\frac{d\bar{m}}{d\bar{r}}$	$\displaystyle=4\pi\bar{r}^{2}\bar{\rho}(\bar{r})+\frac{\bar{q}(\bar{r})}{\bar{%r}}\frac{d\bar{q}}{d\bar{r}}$	(26)
$\displaystyle\frac{d\bar{P}}{d\bar{r}}$	$\displaystyle=(\bar{P}(\bar{r})+\bar{\rho}(\bar{r}))\frac{\left(\bar{r}^{3}(%\Gamma+8\pi\bar{\alpha}\bar{P}(\bar{r})-1)-2\bar{\alpha}\bar{m}(\bar{r})\right%)}{\Gamma\bar{r}^{2}\left((\Gamma-1)\bar{r}^{2}-2\bar{\alpha}\right)}+\frac{%\bar{q}(\bar{r})}{4\pi\bar{r}^{4}}\frac{d\bar{q}}{d\bar{r}}$	(27)

	$\displaystyle\phi^{\prime}(r)$	$\displaystyle=-\sqrt{\frac{m(0)}{2\alpha r}}-\frac{3m(0)}{4\alpha}+\frac{m(0)%\sqrt{r}\left(\alpha\left(\tilde{q}(0)^{2}-2m^{\prime}(0)\right)-5m(0)^{2}%\right)}{4\sqrt{2}(\alpha m(0))^{3/2}}$		(29)
		$\displaystyle+\frac{r\left(4\alpha\left(-6m^{\prime}(0)+3\tilde{q}(0)^{2}+2%\right)-35m(0)^{2}\right)}{32\alpha^{2}}+\mathcal{O}\left(r^{3/2}\right)$		(29)

	$\displaystyle\lim_{r\to 0}\phi^{\prime}(r)$	$\displaystyle\sim-\sqrt{\frac{\mathcal{M}(0)}{2\alpha}}\sqrt{r}+\frac{r}{4%\alpha}+\mathcal{O}\left(r^{3/2}\right)\approx 0$		(30)
	$\displaystyle\lim_{r\to 0}\phi(r)$	$\displaystyle\sim-\frac{r^{3/2}}{3}\sqrt{\frac{2\mathcal{M}(0)}{\alpha}}+\frac%{r^{2}}{8\alpha}+K\approx K$		(30)