Flat Horizon D-Dimensional Black Holes in the Cubic Form of f(Q) Theory

1. Introduction

One notable and worrisome observation from the last 20 years is the universe’s acceleration caused by Dark Energy (DE). This cosmological event is confirmed by recent developments in observational cosmology: Cosmic Microwave Background Radiation[3], type Ia supernovae [4,5,6], the Lyman-forest power spectrum from the Sloan Digital Sky Survey[7], large-scale structure observations [8,9,10], and the investigation of high-energy DE models with weak lensing data [11]. However, GR is unable to account for observations of massive pulsars [12,13] and white dwarfs [14,15,16] with masses higher than the Chandrasekhar mass limit, or $1.44M_{⊚}$, which is the customary maximum limit. Moreover, GR is contradicted by the strong gravitational field and fresh findings [17,18,19]. Consequently, scientists search for appropriate modifications to the GR, such as gravitation $f\left(R\right)$[20,21,22,23,24,25,26], $f(R,G)$, where the Gauss-Bonnet and Ricci scalar expressions are denoted by R and G, respectively [27], $f\left(T\right)$ gravity, with the torsion scalar T[28,29,30], $f\left(G\right)$ gravity [31,32,33], Brans-Dicke (BD) gravity [34,35], and so on. The concept of higher-order curvature, or more precisely $f\left(R\right)$ gravity, is the most successful adaptation of GR, which uses data to refute the theory of gravity and explain the existence of dark matter[36]. Recently, $f\left(Q\right)$ gravity, a well-motivated theory of gravity, was put forth by Jim’enez et al. [37]. Lagrangian density, on which it is based, produces a general function of the non-metricity scalar Q. Non-metricity drives the gravitational interaction in space-time in this theory. The modified theory of $f\left(Q\right)$ gravity leads to intriguing cosmic phenomenology at the background level [38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57].

Moreover, it has effectively been tested against diverse observational data related to background and perturbations, Type Ia Supernovae (SNIa), including the Cosmic Microwave Background (CMB), Redshift Space Distortion (RSD), growth data, Baryonic Acoustic Oscillations (BAO), and similar datasets [58,59,60,61,62,63,64]. Ultimately, the constraints of Big Bang Nucleosynthesis (BBN) are simply transcended by $f\left(Q\right)$ gravity[65]. Different yet comparable theories of gravity are produced depending on T (torsion) or Q (nonmetricity). These are known as the teleparallel equivalent of GR, or $\left(TEGR\right)$[66,67] and symmetric teleparallel GR (STGR) [68,69,70]. Instead of curvature and torsion, nonmetricity underpins the concept of gravity in STGR. Inspired by the interesting qualities of $f\left(Q\right)$ gravity, we are going to derive $D-dimensions$ flat horizons black hole using the cubic form of $\mathrm{f}\left(Q\right)$, i.e., $\mathrm{f}\left(Q\right)=Q+1/2\gamma Q^2+1/3{\gamma }_1{Q}^3-2\Lambda$, where $\gamma$ and ${\gamma }_1$ are two dimensional constants of $\mathrm{l}engt{h}^2$$\mathrm{l}engt{h}^3$ and $\Lambda$. $\Lambda$ in this study represents the cosmological constant.

The structure of this investigation is outlined as follows: Section 2, deals with the examination of the field equations and a brief summary of the non-mitricity formalism. Subsequently, we present the equation of motion for gravity within the framework of $f\left(Q\right)$. The ansatz of the metric with a flat horizon in D-dimensions is utilized to the equations of motion of $f\left(Q\right)$ gravity in Section 3. Applying this approach leads to deriving a new solution in D-dimensions. The asymptotic behavior of the solution corresponds to Anti-de-Sitter (AdS) space. The relevant physical properties of these solutions are discussed in Section 4. We explore the black holes thermodynamics in Section 5. Finally, Section 6 has closing comments.

2. The theory of $f\left(Q\right)$

This section covers some of the generic characteristics of $f\left(Q\right)$-gravity. We will restrict the scope of this explanation to components (the reader can turn to Ref. [71,72,73] for a more rigorous derivation in terms of forms).

For a parallelizable and differentiable manifold, the affine connection can be written in the form:

$${\Gamma }_{\mu \nu }^{\sigma }={\tilde{\Gamma }}_{\mu \nu }^{\sigma }+K_{\mu \nu }^{\sigma }+L_{\mu \nu }^{\sigma },$$

(1)

where the Levi-Civita connection is represented by ${\tilde{\Gamma }}_{\mu \nu }^{\sigma }$ which has the following definition:

$${\tilde{\Gamma }}_{\mu \nu }^{\sigma }=\frac{1}{2}{g}^{\sigma \rho }\left({\partial }_{\mu }{g}_{\rho \nu }+{\partial }_{\nu }{g}_{\rho \mu }-{\partial }_{\rho }{g}_{\mu \nu }\right).$$

(2)

Also, the contortion ${K^{\sigma }}_{\mu \nu }$ is defined as:

$$K_{\mu \nu }^{\sigma }=\frac{1}{2}{T}_{\mu \nu }^{\sigma }+T_{\left(\mu \nu \right)}^{\sigma },$$

(3)

with the torsion tensor $T_{\mu \nu }^{\sigma }=2{\Gamma }_{\left[\mu \nu \right]}^{\sigma }$. Lastly, the deformation is $L_{\mu \nu }^{\sigma }$, which reads,

$$L_{\mu \nu }^{\sigma }=\frac{1}{2}{Q}_{\mu \nu }^{\sigma }-Q_{\left(\mu \nu \right)}^{\sigma },$$

(4)

where the non-metricity tensor, $Q_{\mu \nu }^{\sigma }$, is provided by

$$Q_{\sigma \mu \nu }={\nabla }_{\sigma }{g}_{\mu \nu }={\partial }_{\sigma }{g}_{\mu \nu }-{\Gamma }_{\sigma \mu }^{\rho }{g}_{\nu \rho }-{\Gamma }_{\sigma \nu }^{\rho }{g}_{\mu \rho }.$$

(5)

Consequently, the scalar of non-metricity is

$$\begin{array}{cc}\hfill Q=& g^{\mu \nu }{(L_{\beta \nu }^{\alpha }{L}_{\mu \alpha }^{\beta }-L_{\alpha \beta }^{\beta }L)}_{\mu \nu }^{\alpha }\equiv Q_{\sigma \mu \nu }{P}^{\sigma \mu \nu },\hfill \end{array}$$

(6)

where $P^{\sigma \mu \nu }$, the conjugate of non-metricity, is given by

$$P_{\mu \nu }^{\sigma }=\frac{1}{4}\left(-Q_{\mu \nu }^{\sigma }+2Q_{\left(\mu \nu \right)}^{\sigma }+Q^{\sigma }{g}_{\mu \nu }-{\tilde{Q}}^{\sigma }{g}_{\mu \nu }-{\delta }_{(\mu }^{\sigma }{Q}_{\nu )}\right),$$

(7)

where $Q_{\sigma }=Q_{\sigma \mu }^{\mu }$, and ${\tilde{Q}}_{\sigma }=Q_{\sigma \mu }^{\mu }$.

In the absence of torsion and non-metricity, the connection takes on the same form as the metrically compatible Levi-Civita connection. Curvature and torsion are both zero in STG, non-metricity is contingent on the interaction between the metric and the connection.

In Ref. [37] authors introduced modified symmetric teleparallel gravity, with the action reading,

$$I=-\frac{1}{2{\kappa }^2}{\int }_{\mathcal{M}}f\left(Q\right)\sqrt{-g}{d}^{4}x+{\int }_{\mathcal{M}}{\mathcal{L}}_m\sqrt{-g}{d}^{4}x,$$

(8)

where g is the determinant of the metric tensor $g_{\mu \nu }$, $\mathcal{M}$ is the space-time manifold, ${\mathcal{L}}_m$ is the Lagrangian density of matter contents, and $f\left(Q\right)$ is a generic function of the non-metricity scalar Q.

One applies to Eq. (8) independent variations with respect to both the metric and the connection in order to obtain the field equations of the theory, having so

$$\xi \equiv \frac{2}{\sqrt{-g}}{\nabla }_{\alpha }\left(\sqrt{-g}{f}_Q{P}_{\mu \nu }^{\alpha }\right)+\frac{1}{2}{g}_{\mu \nu }f+f_Q\left(P_{\mu \alpha \beta }{Q}_{\nu }^{\alpha \beta }-2P_{\alpha \beta \mu }{Q}_{\nu }^{\alpha \beta }\right)={\kappa }^2{\mathcal{T}}_{\mu \nu },$$

(9)

$${\nabla }_{\mu }{\nabla }_{\nu }\left(\sqrt{-g}{f}_Q{P}_{\mu \nu }^{\alpha }\right)=0,$$

(10)

where ${\mathcal{T}}_{\mu \nu }$, as is traditional, represents the energy-momentum tensor of matter

$${\mathcal{T}}_{\mu \nu }=-\frac{2}{\sqrt{-g}}\frac{\delta \left(\sqrt{-g}{\mathcal{L}}_m\right)}{\delta g^{\mu \nu }}.$$

(11)

The above expression has two parts: $f_Q=\frac{df\left(Q\right)}{dQ}$ and $f\equiv f\left(Q\right)$. We observe that there is no hyper-momentum because the Lagrangian density of matter is calculated without consideration of the connection. Furthermore, it is well known that by presenting $f\left(Q\right)=Q$, the Lagrangian density $\mathcal{L}=-\frac{Q}{2{\kappa }^2}+{\mathcal{L}}_m$ may be produced, yielding the results of GR (in the STEGR framework).

3. Static anti-de-Sitter black hole solution

We investigate the cylindrical D-dimensional spacetime using the field equations of $f\left(Q\right)$ gravity, given by Eq. (9). The line element that emerges from this analysis is shown in cylindrical coordinates (t, r, ${\zeta }_1$, ${\zeta }_2$,⋯, ${\zeta }_{D-2}$), as elaborated in [74]:

$$\begin{array}{c}\hfill d{s}^2=\mu \left(r\right)d{t}^2-\frac{d{r}^2}{\nu \left(r\right)}-r^2\sum _{i=1}^{D-2}d{\zeta }_i^2.\end{array}$$

(12)

In this context, $\mu \left(r\right)$ and ${\nu }_1\left(r\right)$ denote two variables that depend on the radial coordinate. Additionally, the form of non-mitricity, Q, given by Eq. (12), yields the following form in D-dimension:

$$\begin{array}{c}\hfill Q=-\frac{(D-2)\nu [(D-3)\mu +{\mu }^{\prime }]}{r^2\mu }.\end{array}$$

(13)

By applying Eq. (12) to the equations of motion (9), we derive the following non-zero components:

$$\begin{array}{cc}\hfill & {{\xi }_t}^t\equiv \frac{1}{2r^2{\mu }^2}\left[2{(D-2)}^2{\nu }^2{r}^2{f}_{QQ}[{\mu }^{\prime 2}-\mu {\mu }^{\prime \prime }]+(D-2)r\mu \nu {\mu }^{\prime }[2(D-2)f_{QQ}\{\nu -r{\mu }_1^{\prime }\}+r^2{f}_Q]+\mu \left\{(D-2)r\mu {\mu }_1^{\prime }\left[r^2{f}_Q\right.\right.\right.\hfill \\ \hfill & \left.\left.\left.-2(D-3)(D-2)f_{QQ}\nu \right]+\mu (r^{4}f\left(Q\right)+2(D-2)(D-3)r^2\nu f_Q+4(D-3){(N-2)}^2{\nu }^2{f}_{QQ})\right\}\right]=0,\hfill \\ \hfill & {{\xi }_r}^r\equiv \frac{f\left(Q\right)r^2\mu +2(D-2)r{f}_Q\nu {\mu }^{\prime }+2(D-2)(D-3)f_Q\mu \nu }{2r^2\mu }=0,\hfill \\ \hfill & {{\xi }_{{\zeta }_1}}^{{\zeta }_1}={{\xi }_{{\zeta }_2}}^{{\zeta }_2}=\cdots \cdots ={{\xi }_{{\zeta }_{D-2}}}^{{\zeta }_{D-2}}\equiv \frac{1}{4r^4{\mu }^3}\left\{2r^2\nu \mu {\mu }^{\prime \prime }[\mu (r^2{f}_Q-2(D-2)(D-3)f_{QQ}\nu )-(D-2)r\nu {\mu }^{\prime }{f}_{QQ}]\right.\hfill \\ \hfill & +2(D-2)f_{QQ}{r}^3{\nu }^2{\mu }^{\prime 3}-r^2\mu \nu {\mu }^{\prime 2}[r^2{f}_Q+2f_{QQ}(D-2)(r{\mu }_1^{\prime }-(2D-5)\nu )]+r{\mu }^2{\mu }^{\prime }\left[r{\mu }_1^{\prime }(r^2{f}_Q-6(D-2)(D-3)f_{QQ}\nu )\right.\hfill \\ \hfill & \left.+2(2D-5)r^2{f}_Q\nu +8(D-2)(D-3)f_{QQ}{\nu }^2\right]+2{\mu }^2\left[(D-3)r\mu {\mu }_1^{\prime }(r^2{f}_Q-2(D-2)(D-3)\nu f_{QQ})+\mu \left(r^{4}f\left(Q\right)\right.\right.\hfill \\ \hfill & \left.\left.\left.+2{(D-3)}^2{r}^2\nu f_Q+4(D-2){(D-3)}^2{f}_{QQ}{\nu }^2\right)\right]\right\}=0.\hfill \end{array}$$

(14)

Following that, we will find a complete solution to Eqs. (14) by using a specific expression for $f\left(Q\right)$, namely:

$$\begin{array}{c}\hfill \mathrm{f}\left(Q\right)=Q+\frac{1}{2}\gamma Q^2+\frac{1}{3}{\gamma }_1{Q}^3-2\Lambda ,\end{array}$$

(15)

where $\gamma$ and ${\gamma }_1$ are dimensional constants that have the unites of $\mathit{l}engt{h}^2$, $\mathit{l}engt{h}^3$, $\Lambda$ is the cosmological constant. Given this specific $f\left(Q\right)$ configuration, the following results are obtained from Eqs. (13):

$$\begin{array}{cc}\hfill & {{\xi }_t}^t\equiv \frac{1}{2r^2{\mu }^2}\left\{2{(N-2)}^2{\nu }^2{r}^2[\gamma +2{\gamma }_{1}Q][{\mu }^{\prime 2}-\mu {\mu }^{\prime \prime }]+(N-2)r\mu \nu {\mu }^{\prime }[2(N-2)[\gamma +2{\gamma }_{1}Q]\{\nu -r{\mu }_1^{\prime }\}+r^2(1+\gamma Q+{\gamma }_1{Q}^2)]\right.\hfill \\ \hfill & +\mu \left[(N-2)r\mu {\mu }_1^{\prime }\left[r^2[1+\gamma Q+{\gamma }_1{Q}^2]-2(N-3)(N-2)[\gamma +2{\gamma }_{1}Q]\nu \right]+\mu \left\{r^4[Q+\frac{1}{2}\gamma Q^2+\frac{1}{3}{\gamma }_1{Q}^3-2\Lambda ]\right.\right.\hfill \\ \hfill & \left.\left.\left.+2(N-2)(N-3)r^2\nu [1+\gamma Q+{\gamma }_1{Q}^2]+4(N-3){(N-2)}^2{\nu }^2[\gamma +2{\gamma }_{1}Q]\right\}\right]\right\}=0,\hfill \\ \hfill & {{\xi }_r}^r\equiv \frac{[Q+\frac{1}{2}\gamma Q^2+\frac{1}{3}{\gamma }_1{Q}^3-2\Lambda ]r^2\mu +2(N-2)r[1+\gamma Q+{\gamma }_1{Q}^2]\nu {\mu }^{\prime }+2(N-2)(N-3)[1+\gamma Q+{\gamma }_1{Q}^2]\mu \nu }{2r^2\mu }=0,\hfill \\ \hfill & {{\xi }_{{\zeta }_1}}^{{\zeta }_1}={{\xi }_{{\zeta }_2}}^{{\zeta }_2}=\cdots \cdots ={{\xi }_{{\zeta }_{D-2}}}^{{\zeta }_{D-2}}\equiv \frac{1}{4r^4{\mu }^3}\left\{2r^2\nu \mu {\mu }^{\prime \prime }\left[\mu \left(r^2[1+\gamma Q+{\gamma }_1{Q}^2]-2(N-2)(N-3)[\gamma +2{\gamma }_{1}Q]\nu \right)\right.\right.\hfill \\ \hfill & \left.-(D-2)r\nu {\mu }^{\prime }(\gamma +{\gamma }_{1}Q)\right]+2(D-2)[\gamma +2{\gamma }_{1}Q]r^3{\nu }^2{\mu }^{\prime 3}-r^2\mu \nu {\mu }^{\prime 2}\left[r^2[1+\gamma Q+{\gamma }_1{Q}^2]+2[\gamma +2{\gamma }_{1}Q][r{\mu }_1^{\prime }-(2D-5)\nu ]\right.\hfill \\ \hfill & \left.\times (N-2)\right]+r{\mu }^2{\mu }^{\prime }\left[r{\mu }_1^{\prime }(r^2[1+\gamma Q+{\gamma }_1{Q}^2]-6(N-2)(N-3)[\gamma +{\gamma }_{1}Q]\nu )+2(2D-5)r^2[1+\gamma Q+{\gamma }_1{Q}^2]\nu \right.\hfill \\ \hfill & \left.+8(D-2)(D-3)[\gamma +{\gamma }_{1}Q]{\nu }^2\right]+2{\mu }^2\left[(N-3)r\mu {\mu }_1^{\prime }(r^2[1+\gamma Q+{\gamma }_1{Q}^2]-2(D-2)(D-3)\nu [\gamma +2{\gamma }_{1}Q])\right.\hfill \\ \hfill & \left.\left.+\mu \left(r^4[Q+\frac{1}{2}\gamma Q^2+\frac{1}{3}{\gamma }_1{Q}^3-2\Lambda ]+2{(N-3)}^2{r}^2\nu [1+\gamma Q+{\gamma }_1{Q}^2]+4(D-2){(D-3)}^2[\gamma +2{\gamma }_{1}Q]{\nu }^2\right)\right]\right\}=0.\hfill \end{array}$$

(16)

For Eq. (16), a general D-dimensional case solution is:

$$\begin{array}{cc}\hfill & \mu \left(r\right)=\nu \left(r\right)=\frac{1}{3}\left(\frac{1}{20}\frac{\sqrt[3]{600{{\gamma }_1}^2\Lambda -90{\gamma }_1\alpha +27{\gamma }^3+10\sqrt{80{\gamma }_1-27{\gamma }^2-1080{\alpha }_1\gamma \Lambda +3600{{\gamma }_1}^2{\Lambda }^2+324\Lambda {\gamma }^3}{\gamma }_1}}{{\gamma }_1}\right.\hfill \\ \hfill & \left.+\frac{r^2}{20}\frac{9{\gamma }^2-20{\gamma }_1}{{\gamma }_1\sqrt[3]{600{{\gamma }_1}^2\Lambda -90{\gamma }_1\gamma +27{\gamma }^3+10\sqrt{80{\gamma }_1-27{\gamma }^2-1080{\gamma }_1\gamma \Lambda +3600{{\gamma }_1}^2{\Lambda }^2+324\Lambda {\gamma }^3}{\gamma }_1}}+\frac{3}{20}\frac{\gamma }{{\gamma }_1}\right)+\frac{c_1}{r}.\hfill \end{array}$$

(17)

The general behavior of the above solution shows that $\mu$ and $\nu$ behave generally as anti-de-Sitter(AdS) or de-Sitter (dS) spacetime. Here, $c_1$ stands for a dimensional integration constant.In an effort to streamline the computations, we shall assume that

$$\begin{array}{c}\hfill \Lambda =\frac{1}{12\gamma },{\alpha }_1=\frac{2{\gamma }^2}{5}.\end{array}$$

(18)

This presumption results in a special solution that has the following form:

$$\begin{array}{ccc}\hfill & & \mu \left(r\right)=\nu \left(r\right)=\frac{r^2}{(D-1)(D-2)\alpha }+\frac{c_1}{r^{N-3}}.\hfill \end{array}$$

(19)

It is evident from Eq. (19) that in the case of cubic form, the higher order of $f\left(Q\right)$ serves as a cosmological constant.

4. The fundamental features of the black hole solutions (19)

Let’s investigate now certain relevant facets of the solution discussed in the previous section. The formulation of the solution’s line element (19) is as follows:

$$\begin{array}{cc}\hfill & d{s}^2=\left[r^2{\Lambda }_1-\frac{2M}{r^{D-3}}\right]d{t}^2-\frac{d{r}^2}{r^2{\Lambda }_1-\frac{2M}{r^{D-3}}}-r^2\sum _{i=1}^{D-2}d{\zeta }_i^2,\hfill \\ \hfill \end{array}$$

(20)

with ${\Lambda }_1$ being the cosmological constant that related to the theory of $f\left(Q\right)$ and is defined as ${\Lambda }_1=\frac{1}{(D-1)(D-2)\gamma }$ and $c_1=-2M$. Equation (19) clearly signifies that the line element of the solution approaches AdS geometry. It’s important to highlight that there is no counterpart for the non-metricity solution in its linear form as we approach the limit of $\gamma \to 0$.

Singularity:

Physical singularities in this framework are identified by assessing all possible invariants within the domain of $f\left(Q\right)$ theory. The ansatz $\mu \left(r\right)$ might exhibit roots, represented as $r_h$. Consequently, one must investigate the invariant behavior around these roots. Following the evaluation from the different invariants, we get:

$$\begin{array}{cc}\hfill & Q=-\frac{1}{\gamma },\hfill \\ \hfill & R^{\mu \nu \lambda \rho }{R}_{\mu \nu \lambda \rho }=\frac{2N}{(N-1){(N-2)}^2{\gamma }^2}+\frac{(D-1){(D-2)}^2(D-3)M^2}{r^{2(D-1)}},\hfill \\ \hfill & R^{\mu \nu }{R}_{\mu \nu }=\frac{D}{{(D-2)}^2{\gamma }^2},\mathbb{R}=\frac{D}{(D-2)\gamma },\hfill \\ \hfill & Q^{\mu \nu \rho }{Q}_{\mu \nu \rho }\approx -\frac{{(D-3)}^2}{(D-1)\alpha }+\frac{10(D-3)M}{r^{D-1}}+\mathcal{O}\left(r^{D-1}\right),\hfill \\ \hfill & Q_{\mu }{Q}^{\mu }=\frac{4(D-2)}{(D-1)\gamma }+\frac{4{(D-2)}^2}{r^{(D-1)}},\hfill \end{array}$$

(21)

where $R^{\mu \nu \lambda \rho }R\mu \nu \lambda \rho$, $R^{\mu \nu }Q\mu \nu$, Q, $Q^{\mu \nu \lambda }{Q}_{\mu \nu \lambda }$, $Q^{\mu }{Q}_{\mu }$, ${\tilde{Q}}^{\mu }{\tilde{Q}}_{\mu }$, and Q represent all the conceivable invariants that can be formulated within this theory1 which demonstrates the singularity of the invariants at $r=0$, which is described as a singularity in curvature.

5. The black holes thermodynamic properties as expressed by Eq. (19)

Using the recently found solution given in Eq. (20), we explored the thermodynamic properties by introducing the concept of the Hawking temperature [75,76] as:

$$T_h=\frac{{\mu }^{\prime }\left(r_h\right)}{4\pi }.$$

(22)

The notation ′ in this scenario signifies a derivative in relation to the event horizon, $r_h$, which represents the most significant positive root of $\nu \left(r_h\right)=0$, while ensuring that ${\nu }^{\prime }\left(r_h\right)$ is not equal to zero. The $f\left(Q\right)$ theory’s Bekenstein-Hawking entropy is expressed as [77,78]2:

$$S\left(r_h\right)=\frac{1}{4}A{f}_Q\left(r_h\right).$$

(23)

The event horizon’s surface area in this frame is denoted by A. As per the heat capacity indicator $C_h$, the black hole will be thermodynamically stable; if $C_h>0$, it will be stable, and if $C_h<0$, it will not be stable. In the next study, we determine if these black hole solutions are thermally stable by looking at how each of their distinct heat capacity behaves.[79,80]

$$H_h=\frac{d{E}_1}{d{T}_h}=\frac{\partial M}{\partial r_h}{\left(\frac{\partial T}{\partial r_h}\right)}^{-1}.$$

(24)

In this frame, $E_1$ describes the quasilocal energy. Within the framework of four dimensions and in relation to the solution given in Eq. (20), the horizons are derived as follows:

$$\begin{array}{cc}\hfill & r_h=\sqrt[3]{12M\gamma }.\hfill \end{array}$$

(25)

Furthermore, we can obtain the following mass equation from Eq. (20):

$$\begin{array}{cc}\hfill & M=\frac{{r_h}^3}{12\gamma }.\hfill \end{array}$$

(26)

The black hole’s total mass is influenced by the horizon, as demonstrated by equation (26). Figure 1(a) illustrates the relationship between $\nu \left(r\right)$ and r, illustrating the potential horizons.

The entropy of solution (20) takes the following manner:

$$\begin{array}{c}\hfill S_h=\frac{2\pi {r_1}^2}{5}.\end{array}$$

(27)

The patterns of entropy are depicted in Figure 1 (b), revealing a consistent behavior of the entropy.

The following formula is used to get the Hawking temperature of Eq. (20):

$$\begin{array}{cc}\hfill & T_h=\frac{{r_h}^3+6M\gamma }{2\pi r_h(12M\gamma -{r_h}^3)}.\hfill \end{array}$$

(28)

where, $T_h$ represents the Hawking temperature. The Fig. (c), shows the temperature revealing that it is always positive.

Equations (26) and (28) are substituted into (24) to yield:

$$\begin{array}{c}\hfill H_h=\frac{\pi {r_h}^4{({r_h}^3-12M\gamma )}^2}{2\gamma ({r_h}^6+48{r_h}^{3}M\gamma -72M^2{\gamma }^2)}.\end{array}$$

(29)

Figure 1(d) shows the patterns of the heat capacity for solution (20) for several values of the model parameters. The heat capacity is consistently positive as long as $r>h_h$, suggesting higher global stability, as seen in Figure 1(d).

6. Conclusions and discussion

In this study, we delved into the intricate realm of the cubic form of the $f\left(Q\right)$ gravity theory, aiming to unravel its implications and characteristics. The cubic form, encapsulated by the function $f\left(Q\right)$ where Q represents the non-metricity, introduces a compelling dimension to our understanding of gravitational dynamics. Throughout our investigation, we scrutinized various aspects, including the solutions to the field equations, the behavior of invariants, and the implications for spacetime geometry.

One of the pivotal findings of our study pertains to the reality that the dimensional quantities relate to the quadratic and cubic higher-order theories, i.e., $\gamma$ and ${\gamma }_1$ are finally unify to show behavior of cosmological constant. This sheds light on the nuanced interplay between the cubic form of $f\left(Q\right)$ gravity and the fundamental aspects of gravitational physics which indeed ensure that $f\left(Q\right)$ will generally do not different from the first-order approximation of $f\left(Q\right)$ in the case of static geometry.

It’s crucial to acknowledge the limitations of our study, such as the charged geometry. This aspect warrant further investigation and refinement in future research endeavors.

In conclusion, our study on the cubic form of $f\left(Q\right)$ gravity theory represents a step forward in comprehending the complexities of alternative gravitational theories. The intriguing patterns and phenomena uncovered in this exploration pave the way for continued research and offer valuable contributions to the broader landscape of gravitational physics.