The extended minimal geometric deformation of SU($N$) dark glueball condensates

The extended minimal geometric deformation (EMGD) procedure, in the holographic membrane paradigm, is employed to model stellar distributions that arise upon self-interacting scalar glueball dark matter condensation. Such scalar glueballs are SU($N$) Yang-Mills hidden sectors beyond the Standard Model. Then, corrections to the gravitational wave radiation, emitted by SU($N$) EMGD dark glueball stars mergers, are derived, and their respective spectra are studied in the EMGD framework, due to a phenomenological brane tension with finite value. The bulk Weyl fluid that drives the EMGD is then proposed to be experimentally detected by enhanced windows at the eLISA and LIGO.


I. INTRODUCTION
The membrane paradigm of brane-world scenarios can be deployed into the fluid/gravity correspondence as the low energy regime of AdS/CFT conjecture [1], which, briefly, states that a theory of gravity in an anti-de Sitter (AdS) space of d + 1 dimensions is dual to a strongly coupled conformal field theory (CFT) -whose low energy underlying hydrodynamics corresponds to the Navier-Stokes equations -at the d-dimensional border of the AdS space [2][3][4]. In the membrane paradigm of the fluid/gravity correspondence [5], black holes in the brane, as well as black branes and black strings, studied from the viewpoint of fluid membranes [6,7] in the long wavelength limit [8,9]. Other successful paradigms are also described in the seminal Refs. [10,11].
In the membrane paradigm and beyond the general relativity (GR), the so called method of geometrical deformation (MGD) places itself as an important procedure to generate new solutions of the effective Einstein's field equations on the brane [8,[12][13][14], describing black holes and compact stellar distributions as well, in a 5D bulk Weyl fluid bath [15,16]. The MGD and its extensions take into account the brane Einstein's field equations [17,18], where the effective stress-energy tensor has additional terms, in particular regarding the Gauss-Codazzi equations from the bulk stress-energy tensor into the brane [10]. Important terms are the bulk dark radiation, the bulk dark pressure, the electric part of the Weyl tensor and quadratic terms on the brane stress-energy tensor. This last one is derived for regimes of energy that are further than the (finite) brane tension in the theory. The MGD and its extensions [6] have been recently equipped with experimental, phenomenological, and ob-servational very precise bounds on the running parameters. MGD gravitational lensing effects were explored in Ref. [19] and the classical tests of GR imposed the most strict bound on the brane tension in Ref. [16]. Besides, the configurational entropy was used to provide account for the critical stellar densities, in the MGD, in full compliance with the Chandrasekhar's critical stellar densities that are extremal points of the system associated informational entropy [20]. Besides, MGD black holes analogues were explored in Ref. [21] where sound waves through de Laval nozzles in a laboratory were shown to be suitable to derive experimental data regarding the 5D bulk Weyl fluid, being acoustic perturbations indeed the analogue MGD quasinormal modes. MGD black branes was also studied in Ref. [22].
On the other hand, gauge forces, that are not encompassed by the Yang-Mills fields of the Standard Model, might play a prominent role in describing features of the observed universe. Indeed, the Standard Model may present a coupling to hidden sectors ruled by a pure Yang-Mills setup, in the low energy regime. Non-Abelian dark forces can be then implemented by a (pure) Yang-Mills setup, with simple gauge group, confining at the energy scale of the theory with dark gluons. This system is then bound into some spectrum of dark glueballs, that are colour-neutral, that have mass or the order of the energy scale [23][24][25]. The scalar glueball dark matter is then a prime candidate emulating SU(N ) dark gauge hidden sectors [23,[26][27][28], whose cross section is very large. When bosons and fermions in the Standard Model are precluded to interact with the SU(N ) scalar glueball whatsoever, a Bose-Einstein condensation of glueballs can occur, originating thus compact stellar distributions. In particular, in the holographic AdS/CFT interpretation, Bose-Einstein condensates of scalar glueballs occupy a relevant place [29]. Ref. [28] discusses relevant elastic scatterings among glueballs, manifesting their feasibility as a self-interacting dark matter candidate. Dark SU(N ) glueball stars were studied on fluid Eötvös branes in Ref. [7], with proposed experimental signatures provided by gravitational waves their mergers. Hereon, the extended MGD (EMGD) procedure [6] is proposed to explore physical signatures of SU(N ) dark glueball condensates. This paper is organised as follows: Sect. II is devoted to introduce the EMGD method, reviewing the MGD one, through the deformation of the radial component of the metric, presenting solutions for the brane Einstein's field equations. The EMGD is implemented by taking into account the deformation of the metric temporal component, given by a k parameter, phenomenologically driving the extension of the MGD. The Schwarzschild solution is obtained for k = 0 and, for k > 0, terms up to order O(r −(k+1) ) are considered. In Sect. III, the self-interacting system of glueball dark matter is studied, with a self-interacting scalar glueball potential that induces condensation into SU(N ) dark EMGD stellar distributions. Then, the corrections to the gravitational wave radiation that is emitted by SU(N ) dark EMGD glueball mergers are derived, being their spectra obtained for two important cases in the EMGD setup, due to a phenomenological brane tension with finite value, opening a wider range probed by the eLISA and the LIGO. In Sect. IV the concluding remarks are outlined with the perspectives.

II. EMGD STELLAR DISTRIBUTIONS
The Minimal Geometric Deformation (MGD) procedure can be realized as a mechanism employed to derive high energy corrections to the general relativity (GR), in such a way that the Einstein's field equations nonlinearity does not produce inconsistencies in the solutions. In fact, the AdS/CFT correspondence can bind warped extra-dimensional models to 4D theories that are strongly-coupled. According to the membrane paradigm, which has been used to realize the deformation method, our 4D universe is a brane with a finite positive tension (or energy density) σ [10]. Systems with energy E σ neither feel the self-gravity effects nor the bulk effects, which allows then the recovery of GR in such regime, corresponding to an infinitely rigid brane, wherein σ → ∞. The most strict brane tension bound, σ 3.18 × 10 6 MeV 4 , was obtained in the MGD context [20].
The study the EMGD, the 5D Einstein Field Equations (EFE) in the bulk (5) must be computed, where (5) G AB denotes the 5D Einstein tensor and (5) T AB is the 5D stress-energy tensor. Natural units 8πG = c = 1 = are used hereon, and M Pl shall denote the 4D Planck mass. Greek indexes µ, ν represent the 4D brane indexes running from 0 to 3, whilst A, B = 0, . . . , 4 denote bulk indexes. The Gauss-Codazzi equations can be used to implement the projection of Eq. (1) onto the brane [17], re-sulting in the effective EFE on the brane. The projected EFE presents high energy corrections evinced by the bulk permeated by a 5D Weyl fluid, whose projection onto the brane yields the electric part of the Weyl tensor E µν . The interaction between the brane and the bulk is strictly gravitational, with no exchange of further fields. Therefore, the projection of the 5D Weyl fluid on the brane encrypts the effects due to 5D gravitons for the linearized case, namely, the Kaluza-Klein modes. Effects provided by E µν are therefore considered non-local, since it does not depend on any data on the brane, having the form [10], where h µν = g µν +u µ u ν is the orthogonal projector to the 4-velocity u µ onto the brane, being g µν the brane metric components. Besides, U = − 1 6 σE µν u µ u ν is the effective energy density; P µν = − 1 6 σ h ρ (µ h σ ν) − 1 3 h ρσ h µν E µν is the effective non-local anisotropic stress-tensor; and the effective non-local energy flux on the brane, Q µ = − 1 6 σh ρ µ E ρν u ν , is originated from the bulk free gravitational field. Local corrections are encoded in the tensor: where T µν is the brane matter stress-tensor, in the MGD setup represented by a perfect fluid: with ρ being the fluid density and p represents the fluid pressure field. Notice that S µν is a correction arising from the bulk, quadratic in T µν . The 4D effective EFE then read: Since E µν ∼ σ −1 , it is straightforward to notice that in the infinitely rigid brane limit, σ → ∞, GR is recovered and G µν = T µν . The brane tension has such important role in the MGD, which provides a way to check that the deformation chosen in the method recovers the standard GR equations. Compact stars in 4D, which must be solutions of Eq. (5), can be described by a static, spherically symmetric metric, written in Schwarzschild-like coordinates as: By solving Eq. (1), using the metric of Eq. (6) and the stress-energy tensor (4), one obtains the effective pressure components of the Weyl fluid: where either a = 0, counting the temporal component, or a = 1, if the radial component is taken into account, and where P = P µ µ . Effects of the 5D Weyl fluid makes the fluid projected on the brane to present a certain degree of anisotropy [12][13][14], encoded into the equation: Eq. (8) asserts that the difference between the components of the fluid pressure is a function of σ −1 , being the fluid isotropic when the GR limit σ −1 → 0 is recovered. Considering the inner region of a stellar distribution and the outer vacuum, the radial coordinate r ranges from the center of the star (r = 0) to the star surface (r = R), and then beyond into the outer vacuum, where ρ = p = 0 in Eq. (4). We shall encode such a condition into the following equation, which is the general standard GR expression for the radial component of the metric: where m(r) = 8πG 4 r 0ř 2 ρdř stands for the GR mass function. In fact, in the GR limit 1/σ → 0, the function M = M (σ) in Eq. (9) is led to the standard GR mass value, M 0 . This limit yields M | σ→∞ = M 0 = m(R).
For the metric of Eq. (6) to be a solution of Eq. (5), one assumes the ansatz for its radial component [13]: where the geometric deformation given by f (r) is what an observer on the brane experiences due to the projected 5D bulk gravity effects onto the brane. In fact, f (r) is a distorting function added to the standard GR solution, given by ξ(r), therefore yielding a deformation, reading [13]: where one denotes ∂ r = d/dr, H = H(p, ρ, ν), χ = χ(σ) is an integration constant, and with r 0 chosen according to the region of interest. The most important contribution in the geometric deformation comes from the function H, In fact, the field H equals zero for the standard GR solution, yielding an infimum for the MGD, since H > 0 [16]. Thus, in this limit, the geometrical deformation shall be solely driven by the pressure and the energy density of the source. When it sets in, the MGD,f (r), is explicitly given by: Therefore, starting from the choice ν = ν GR , one can derive the deformed radial component of the metric by evaluating Eq. (10) with the minimal deformationf (r) [8]. Afterwards, the deformed temporal metric component is derived with the remaining expressions from the effective EFE relating λ and ν. The field χ in Eq. (11) must be zero in the standard GR limit. Besides, the integration constant χ must equal to zero in the star inner region (r < R), for the metric to be smooth at r = 0. Notwithstanding, at r > R, χ can assume values that are different of zero. Hence, a geometrical deformation can be associated with the standard GR Schwarzschild ("S") solution: e νS = e −λ S = 1 − 2M r , in which case the MGD field becomesf (r)| p=ρ=0 = χe −I(r) , implying [12]:f . In Ref. [16] such a parameter was phenomenologically weaker bound as |b(M )| 5 × 10 −11 . Therefore, the deformed exterior temporal and radial metric components are respectively given by: MGD black holes, described by this metric, has two event horizons, r 1 = 2M and r 2 = 3M 2 − ψ. Thus, up to the second order on σ −2 , the deformation term in Eq. (10) yields f (r) = − ψ r . For a stellar distribution, 5D bulk effects are maximal, near the MGD stellar distribution surface r = R. The parameter |ψ| stands for the star compactness, being higher for higher values of ψ. Observational data imply the strongest bound |ψ| 3.18 × 10 −11 [16], hence matching good results up to the order O(σ −2 ). The MGD procedure was used in Ref. [7] to model dark SU(N ) stars and to study their experimental signatures. Furthermore, another application was accomplished in the context of de Laval nozzles, which can be associated to MGD black hole analogues, so that 5D effects might be observed in the laboratory [21]. Now, an extension of the MGD can be regarded [6], generalizing the MGD for the outer region. It can be accomplished by considering a deformation not only on the radial but also on the temporal metric component. This determines what is called the extended minimal geometric deformation (EMGD) method, which will be discussed in this section, as introduced in [6].
Let us take a geometric deformation of the temporal metric component ν(r) = ν S + h(r) in eq. (6) [12]. The ν S term defines, as previously, the Schwarzschild temporal metric component, and h(r) is the temporal deformation, which is proportional to σ −1 , thus providing the GR limit. Using the effective vacuum EFE, the radial geometric deformation f (r) can be written in terms of h(r) as [6]: where I is given by Eq. (12), and now: Therefore, the exterior deformed radial metric component becomes with the extended geometric deformation f (r) redefined according to Eq. (18). Notice that a constant h implies G = 0, which produces an exterior MGD as one had before. On the other hand, it is also possible to achieve a MGD with a h(r), which is given by setting Eq. (19) equal to zero. The solution of the resulting differential equation reads [6]: With the assumption of asymptotic flatness, lim r→∞ e ν = 1, implying h → 0, the integration constants a and d, both dependent of the brane tension, are constrained by a = 1 − d 2M . Hence, the minimallydeformed temporal metric component becomes [6]: The minimally-deformed radial metric component, on the other hand, reads: A more general solution for the exterior radial metric component of Eq. (20) can be derived [6], under a geometric deformation such that G(h) = 0. The choice h(r) = k log 1 − 2M r yields: where k is a constant known as the deformation parameter. Naturally, k = 0 results no temporal geometric deformation, being directly associated with the Schwarzschild metric. For k = 1, one has which then allows the calculation of the radial metric component, through Eq. (20), which yields [6]: . Now, in order to the radial metric component asymptotically approach the Schwarzschild behavior with ADM mass M = 2M , e −λ(r) ∼ 1 − 2M r + O(r −2 ), one must necessarily have κ 1 = −2M . In this case, the temporal and spatial components of the metric shall be inversely equal to each other (as it is the case of the Schwarzschild solution), containing a tidal charge Q = 4M 2 reproducing a solution that is tidally charged by the Weyl fluid [30]: It is worth to emphasize that the metric of Eq. (27) has a degenerate event horizon at r h = 2M = M. Since the degenerate horizon lies behind the Schwarzschild event horizon, r h = M < r s = 2M, 5D bulk effects are then responsible for decreasing the gravitational field strength on the brane. Now the exterior solution for k = 2 can be constructed, making Eq. (24) to yield: for where Q = 12M 2 and M = 3M . The radial component, on the other hand, reads: where The asymptotic Schwarzschild behavior is then assured when κ 2 = −M/32. In this case, the degenerate event horizon is at r e ≈ 1.12M [6], which makes clear that here as well the 5D bulk effects induce a weaker gravity.
The classical tests of GR applied to the EMGD metric provide the following constraints on the value of the deformation parameter: k 4.5 for the perihelion precession of Mercury, k 4.3, for the deflection of light by the Sun, and k 4.2 for the gravitational redshift of light. Hence, observational and experimental data from the classical tests of GR make us hereon to opt to study the cases k = 1 and k = 2, as above detailed.

III. SU(N ) EMGD GLUEBALL DARK STELLAR DISTRIBUTIONS AND THEIR GRAVITATIONAL WAVE RADIATION SPECTRA
Hidden SU(N ) gauge sectors may be described by the (scalar) glueball dark matter parading [23,28,31]. In fact, glueballs can interact by exchanging gravitons, constituting a self-interacting system. The number, N , of colors driving the SU(N ) gauge sectors and the scalar glueball mass, m, are the parameters that model the glueballs self-interaction. These two parameters drive the 2 → 2 elastic scattering cross sections of the lightest glueball state in the hidden dark sector, given by σ ∼ m −2 N −4 [28]. Dark SU(N ) glueball stars can be formed when the self-gravity unbalances the energy density of the system, inducing a Bose-Einstein condensation of the glueball system into stellar distributions. This phenomenon can occur in the parameters ranges 10 3 N 10 6 and 10 eV m 10 KeV [24,28].
The action for the glueball system is given by The resulting Klein-Gordon-Einstein system can be then scrutinized in what follows, for a φ 4 self-interacting glueball potential, condensating into an EMGD stellar distribution due to a finite brane tension, for λ > 0, making the self-interacting force to counterbalance the gravitational force. Hence, a stable configuration for a SU(N ) dark stellar distribution can be naturally achieved [7,28,32] and therefore the range of frequencies for gravitational waves generated from SU(N ) EMGD dark stars mergers, will be derived and studied. The Klein-Gordon-Einstein system, derived from Eqs. (6,31,32), with the assumption that spherically symmetric glueball fields are also periodic in time, φ(r, t) = Φ(r) exp(iωt) [28], read where [28,32] M(x) m represents the star mass within a radius x m , yielding dimensionless M and x parameters. Besides, one also uses the notation λ = 12λ m 2 and Ω = ω m . The glueball (dark) matter is then modelled by Eq. (33a), yielding Figs. 1 -3 depict the numerical computations, where below we consecutively present the results for the EMGD, for k = 1 and k = 2, comparing with the standard MGD (k = 0) in Eq. (24), and its subsequent GR limit. In what follows, the brane tension value σ ∼ 10 6 MeV 4 shall denote the exact current bound σ 3.18 × 10 −6 GeV 4 derived in Ref. [20]. σ∼∞ σ∼10 6 MeV 4 (k=0) , with respect to e ν(0) Ω 2 . The black line regards the GR limit σ → ∞; the other lines depict the phenomenological bound for the brane tension σ 3.18 × 10 −6 GeV 4 : the standard MGD procedure is a particular case of the EMGD for k = 0 (gray line); the EMGD case is plot for k = 1 (black dashed line) and for k = 2 (gray dashed line).   , with respect to e ν(0) Ω 2 . The dashed black line regards the GR limit σ → ∞; the other lines depict the phenomenological bound for the brane tension σ 3.18 × 10 −6 GeV 4 : the standard MGD procedure is a particular case of the EMGD for k = 0 (black line); the EMGD case is plot for k = 1 (gray line) and for k = 2 (gray dashed line). , with respect to e ν(0) Ω 2 . The dashed black line regards the GR limit σ → ∞; the other lines depict the phenomenological bound for the brane tension σ 3.18 × 10 −6 GeV 4 . The standard MGD procedure is a particular case of the EMGD for k = 0 (black line); the EMGD case is plot for k = 1 (gray line) and for k = 2 (gray dashed line). Fig. 1 shows that a SU(N ) EMGD dark stellar distribution can accrete the dark matter that surrounds it [32], having the stellar mass increased up to a maximal mass, shown in Table I. As general SU(N ) dark glueball stellar distributions have typical mass and radius respectively given by [32], the glueball dark SU(N ) EMGD dark star have, then, the following values for its mass and radius, where M denotes, as usual, the Solar mass. The range of frequencies emitted from SU(N ) (Schwarzschild) dark mergers read f max = 1 2π GM R 3 1/2 [28,33]. SU(N ) EMGD dark stars mergers represent candidates to enlarge such range of frequencies, as we shall show in what follows. Eqs. (39, 40), and Table I, thus imply that the peak gravitational wave frequency reads where the function γ(k, σ) = 4m 2 √ λ c 10 4 /GeV 2 is represented by a factor for σ → ∞ (GR limit), 1.262, for k = 0, 1.521, for k = 1, 1.960, for k = 2.
(42) that measures the corrections further the unit, which regards the σ → ∞ GR limit, for different values of the brane tension. Such a parameter c hence provides corrections to peak wave frequencies from SU(N ) EMGD dark star mergers, in a finite tension membrane paradigm. The range of gravitational wave frequencies is clearly a function that depends on the m glueball dark matter mass and on the number N of colors driving the hidden gauge sectors, and may be detected by the eLISA and the LIGO experiments [33]. The relevant stars m and N parameters ranges imply a maximum EMGD star mass that lies between 10 6 M and 10 9 M , where M denotes the solar mass. On the other hand, the lowest EMGD star radius ranges between 10 2 R and 10 5 R , where R denotes the solar radius. SU(N ) EMGD dark stars hence present a completely distinct experimental signature that are quite different of Schwarzschild black hole mergers, due to the subsequent analysis of Table I, as well as Eqs. (42).
The peak frequency f max , emitted from dark SU(N ) MGD stellar mergers, may be apportioned into a range between 30 µHz and 100 mHz, at the eLISA experiments [34], whereas the LIGO experiment can currently detect frequencies between 50 Hz and 1 KHz. The figures below represent the N -m parameter space. First, Fig. 4 represents the σ → ∞ GR limit, and Figs. 5-7 take into account the MGD (namely, the EMGD for k = 0), the EMGD for k = 1, and the EMGD for k = 2, with finite brane tension lying in the bound σ 3.2 × 10 −6 GeV 4 . LIGO range Self-interacting glueball dark matter eLISA range   [20], for the EMGD procedure, for k = 2. The black [light gray] band shows the gravitational wave peak frequency to be detected by the LIGO [eLISA] experiment, whereas the gray band shows the gravitational wave peak frequency emitted from SU(N ) EMGD dark star mergers.
The bulk Weyl fluid in the membrane paradigm can be, thus, detected, in a scenario where hidden sectors to the Standard Model condensate into SU(N ) glueball stars, described by EMGD stellar distributions. Their limiting cases can be further attained, as k = 0 in Eq. (24) corresponds to the standard MGD and, subsequently, its GR limit σ → ∞. In fact, Fig. 7 represents the detection window for gravitational wave radiation emitted by SU(N ) EMGD dark star mergers, for k = 2 in Eq. (24), that is wider than in Fig. 6, for k = 1. Comparing also with Figs. 4 and 5, the bigger the k, that drives the extension of the MGD, the wider the window for experimental detection is. In the next section we more precisely analyze our results.

IV. CONCLUSIONS
SU(N ) hidden gauge sectors were employed to study the scalar glueball dark matter model, when condensation into SU(N ) EMGD dark stars takes part, in the holographic membrane paradigm. Such a process then generates stable compact stellar distributions, modelled by the extended MGD. After briefly reviewing the MGD and the EMGD procedures, we focus on the deformation of the metric temporal component by the parameter k, in Eq. (24). The phenomenological upper bound k 4.5 makes us to study the cases k = 0 corresponding to the MGD (Fig. 5), whose limit σ → ∞, that regards the Schwarzschild solution was also investigated (Fig. 4). For k = 1 and k = 2, we illustrate our results in Figs. 6 and 7, respectively. Although the case for k = 1 resembles a generalized form for the Reissner-Nordström metric, with tidal charge from a 5D bulk Weyl fluid in Eq. (27) [6], for k = 2 the obtained metric (30) does not take any known form yet. Hence, it may be interesting to study it to estimate the behavior of EMGD solutions in higher orders of k, up to the phenomenological bound. Therefore, SU(N ) EMGD dark glueball stars present both radius and effective mass that are corrected by a bulk 5D Weyl fluid. Figs. 6 and 7 show that observational signatures evinced from EMGD mergers are more probable to be experimentally detectable at LIGO [33] and eLISA project [34]. In fact, the window for detection in Figs. 6 and 7 are wider than for the SU(N ) MGD dark glueball stars (Fig. 5), corresponding to the EMGD for k = 0, and even wider than for SU(N ) Schwarzschild dark glueball stars (Fig. 4).
It shall be still worth to study the effects of the dynamical corrections and the anomalous dimension effects on the scalar glueballs here studied [31], as well as glueballs in finite-temperature AdS/QCD models [35], to further refine the spectra of emitted gravitational wave radiation and its consequences in the membrane paradigm.
Here we proposed that a 5D Weyl fluid, in the membrane paradigm, can induce experimental signatures, at LIGO and eLISA, that are amplified by the EMGD procedure, with respect to the MGD and its σ → ∞ GR limit. Hence, our method is alternative to collider data that depends on confirmation at the LHC or future col-liders [36].