
==== Front
Eur Phys J C Part Fields
Eur Phys J C Part Fields
The European Physical Journal. C, Particles and Fields
1434-6044
1434-6052
Springer Berlin Heidelberg Berlin/Heidelberg

39175581
13203
10.1140/epjc/s10052-024-13203-9
Regular Article – Theoretical Physics
Exact dynamical black hole solutions in five or higher dimensions
Fahim Bardia H.
http://orcid.org/0000-0002-7882-9172
Ghezelbash A. M. masoud.ghezelbash@usask.ca

https://ror.org/010x8gc63 grid.25152.31 0000 0001 2154 235X Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SK S7N 5E2 Canada
20 8 2024
20 8 2024
2024
84 8 8375 6 2024
3 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3.
We construct new classes of the dynamical black hole solutions in five or higher dimensional Einstein–Maxwell theory, coupled to a dilaton field, in the presence of an arbitrary cosmological constant. The dilaton field interacts non-trivially with the Maxwell field, as well as the cosmological constant, with two arbitrary coupling constants. The solutions are non-stationary, and almost conformally regular everywhere. To construct the solutions, we use the four-dimensional Bianchi type IX geometry, as the base space. We find three different classes of solutions, based on the values of the coupling constants. We notice that our solutions could be asymptotically de-Sitter, anti-de-Sitter or flat. We find the relevant quantities of the solutions, and discuss the properties of the solutions.

http://dx.doi.org/10.13039/501100000038 Natural Sciences and Engineering Research Council of Canada issue-copyright-statement© EDP Sciences, Societa Italiana di Fisica (SIF) and Springer-Verlag GmbH, DE, part of Springer Nature 2024
==== Body
pmcIntroduction

Finding the exact solutions to the Einstein gravity, especially in the presence of matter fields in different dimensions, is truly the main aim of gravitational physics. In this regard, the dimensional compactification of higher dimensional gravity has been studied broadly in different references [1, 2]. In other area of research on holography between two different models of physics, especially AdS/CFT and Kerr/CFT correspondences, constructing the exact solutions to the asymptotically de-Sitter and Anti-de-Sitter Einstein gravity is of utmost importance [3]. Moreover, the exact solutions to the Einstein gravity coupled to the different matter fields, such as Maxwell field, dilaton field and NUT charges are constructed in [4–6]. The Einstein–Maxwell-dilaton theory with different types of interaction between the fields, can describe physical phenomenons, such as slowly rotation black holes [7], topological charged hairy black holes [8], cosmic censorship [9], gravitational radiation [10], hyperscaling violation [11] and compactification of M-theory in generalized Freund-Rubin theory [12]. The Einstein–Maxwell-dilaton theory and its extensions have also been used in other areas of research on black holes [13–26].

In a recent work [27], the authors considered the Einstein–Maxwell-dilaton theory with two extra vector fields. The first field supports the non-trivial topology, and the second field supports states with the finite charge density. They found a new class of charged black holes with hyper-scaling violating asymptotics. The black holes have non-trivial horizon topology, for arbitrary Lifshitz exponent and a hyper-scaling violation parameter [27].

Moreover, in [28–31], the authors constructed time-dependent charged black hole solutions in different dimensions in Einstein–Maxwell theory. Some of the time- dependent solutions can describe the coalescence of the extremal charged black holes in different dimensions.

After direct detection of coalescing black holes [32], there is a huge interest in finding the exact analytical dynamical black hole solutions. In general, finding the exact analytical dynamical black hole solutions is a difficult task in general relativity (and its modified versions which include the matter fields).

Inspired by above considerations, in this article we generalize the Kastor–Traschen-like black holes [33, 34] in any dimensions D≥5, based on an embedded four-dimensional Bianchi type IX space. The Kastor–Traschen black holes on Gibbons-Hawking space were constructed in [35, 36], which describe a system of coalescing black holes.

The organization of the paper is as follows. In Sect. 2, we briefly review the Bianchi type IX space and discuss its properties. Then, we consider the Einstein–Maxwell-dilaton theory in N+1-dimensions with two non-equal coupling constants. We employ special ansatzes for the metric, the Maxwell field and the dilaton, and explicitly solve all the field equations. We find analytical exact solutions for all the metric functions, the Maxwell field and the dilaton field. Moreover, we find a constraint on the two non-equal coupling constants and also another constraint on the cosmological constant. We discuss and plot the quantities related to the geometry. In Sect. 3, we consider the Einstein–Maxwell-dilaton theory in N+1-dimensions with two non-zero equal coupling constants. We employ another different set of ansatzes for the metric, the Maxwell field and the dilaton, and explicitly solve all the field equations. We find analytical exact solutions for all the metric functions, the Maxwell field and the dilaton field. Moreover, we find a constraint on the cosmological constant. We discuss and plot the quantities related to the geometry. In Sect. 4, we consider the Einstein–Maxwell-dilaton theory in N+1-dimensions with two zero coupling constants. The theory reduces to the Einstein–Maxwell theory. We employ special ansatzes for the metric and the Maxwell field, and explicitly solve all the field equations. We find analytical exact solutions for all the metric functions and the Maxwell field. We discuss and plot the quantities related to the geometry. In Sect. 5, we show that the solutions to the Einstein–Maxwell-dilaton theory in Sect. 2 can be embedded to a higher-dimensional gravity theory coupled to a form field. We find that the dimension of the the internal space is related to the coupling constant b in the Einstein–Maxwell-dilaton theory. We wrap up the article by an appendix and the concluding remarks and comments on the future works.

Bianchi type IX geometry and Einstein–Maxwell-dilaton theory with two non-equal coupling constants

In Bianchi’s classification of the homogeneous spaces, Bianchi type IX is a self-dual and asymptotically Euclidean space and includes two important sub-spaces such as Eguchi–Hanson I and II. The Bianchi type IX geometry with an SU(2) isometry group is given by1 ds2=e2f(η)σ12+e2h(η)σ22+e2g(η)σ32+e2(f(η)+h(η)+g(η))dη2,

where the Maurer-Cartan one-forms σi satisfy dσi=12ϵijkσjσk, and are given by2 σ1=dψ+cosθdϕ,

3 σ2=cosψsinθdϕ-sinψdθ,

4 σ3=sinψsinθdϕ+cosψdθ.

The periodicity of the Euler angles θ, ϕ and ψ are π, 2π and 4π, respectively. From the self-duality property of Bianchi type IX, we find5 2dfdη=e2h+e2g-e2f-2λ1eh+g,

6 2dhdη=e2g+e2f-e2h-2λ2ef+g,

7 2dgdη=e2f+e2h-e2g-2λ3ef+h,

where the constants λi satisfy λiλj=ϵijkλk, with i=1,2,3. By choosing (λ1,λ2,λ3)=(0,0,0) and solving the vacuum Einstein equations, we find the functions f(η), h(η) and g(η) in terms of the standard Jacobi elliptic functions sn, cn and dn as [37, 38]8 f(η)=12lnc2cn(c2η,k2)dn(c2η,k2)sn(-c2η,k2),

9 h(η)=12lnc2cn(c2η,k2)dn(c2η,k2)sn(-c2η,k2),

10 g(η)=12lnc2dn(c2η,k2)cn(c2η,k2)sn(-c2η,k2).

Changing the coordinate η to r=2csn(c2η,k2), we find the triaxial Bianchi type IX geometry11 dsB.IX2=dr2g(r)+r24g(r)(dψ+cosθdϕ)21-a14r4+(-sinψdθ+cosψsinθdϕ)21-a24r4+(cosψdθ+sinψsinθdϕ)21-a34r4,

where12 g(r)=1-a14r41-a24r41-a34r4,

and we can choose the parameters a1=0, a2=2kc and a3=2c, with constants c>0 and 0≤k≤1. In order to preserve the positive definiteness of the metric (11), we need to impose r≥a3. Bianchi type IX is an important geometry and has been used in different theories such as supergravity, loop quantum cosmology and string theory [39–42].

Having constructed a background space for the theory, we focus on the action of the Einstein–Maxwell-dilaton theory, in N+1 dimensions, where the dilaton field is coupled to the electromagnetic field and the cosmological constant, with two different coupling constants a and b [43]13 S=∫dN+1x-g{R-4N-1(∇ϕ)2-e-4/(N-1)aϕF2-e4/(N-1)bϕΛ},

where R is the Ricci scalar, Fμν=∇μAν-∇νAμ is the electromagnetic field tensor, ϕ is the dilaton field and Λ the cosmological constant. We consider the following ansatz for the N+1-dimensional metric as a background geometry14 dsN+12=-1H(r,θ)2dt2+H(r,θ)2(N-2)R(t)2[dsB.IX2+Σi=1N-4dxi2],

where R(t) and H(r,θ) are two metric functions, Σi=1N-4dxi2 is the extended Euclidean spaces, and dsB.IX2 is the four-dimensional Bianchi type IX metric given by (11).

Variation of the action (13) with respect to the metric tensor gμν, electromagnetic gauge field Aμ and the dilaton field ϕ leads to the Einstein field equations, electromagnetic and dilaton field equations in N+1 dimensions [43],15 Gμν≡Rμν-12gμνR-4N-1∇μϕ∇νϕ-12gμν(∇ϕ)2-e-4aϕN-12FμρFνρ-12gμν(F)2+12e4bϕN-1gμνΛ=0,

16 Mμ≡∇ν(e-4/(N-1)aϕFμν)=0,

17 D≡∇2ϕ-b2e4/(N-1)bϕΛ+a2e-4/(N-1)aϕF2=0.

respectively.

We consider the following ansatz for the dilaton field18 ϕ(t,r,θ)=-(N-1)4aln(HU(r,θ)RV(t)),

where U and V are two constants that will be determined throughout solving the field equations of motion. We also consider the following ansatz for the electromagnetic gauge field19 At(t,r,θ)=αRX(t)HY(r,θ),

where we assumed At(t,r,θ) to be the only non-zero component of the electromagnetic gauge field, which generates an electric field in r and θ directions. In (19), α, X and Y are constants.

The Mr component of the electromagnetic field equation (16) gives the following differential equation20 Mr=-r8-16c4(k4+1)r4+256c8k4r8(∂rH)(∂tR)Yα(X+V+N-2)×RX+V-3HN-4N-2+Y+U=0,

where N is the number of spatial dimension. This Eq. (20) leads to the following relation between the constants X and V21 X+V=2-N.

Moreover, from the Gtr component of the Einstein field equation22 Gtr=-(N-1)∂tR∂rH(UV+4a2)8a2HR,

we find23 UV=-4a2.

Substituting these constraints (21) and (23) in Grθ component of the Einstein field equation, we can determine the constants that appear in dilaton field (18)24 U=2a2N-2,V=-2(N-2),

and the constants in the gauge ansatz (19)25 X=N-2,Y=-1-a2N-2,α2=N-12(a2+N-2).

Substituting these results into the Mt component of the Maxwell field equation (which due to its length, we show it explicitly in the appendix) gives a differential equation for the metric function H(r,θ), that can be solved as26 H(r,θ)=(g+r2cosθ+g-)N-2a2+N-2,

where g± are arbitrary constants. This suggest that for a fixed dimension, increasing the coupling constant |a|, decreases the value of H(r,θ). We also note that based on the line element ansatz in Eq. (14), the function H(r,θ) needs to be a real-positive function. In Fig. 1, we represent the behaviour of the metric function H(r,θ) with respect to the coordinates r and θ, for three different dimensions, where we set the coupling constant a=1.Fig. 1 The metric function H(r,θ) in terms of the coordinates r and θ for three different spatial dimensions N=4, N=5 and N=6, which correspond to the lower, middle and upper surface, respectively. We assumed the constants g+=0.5, g-=15 and a=1

Fig. 2 The behaviour of the cosmological constant in terms of a the spatial dimension N for two different values of the coupling constant a, and b in terms of N for two different values of a, where we set η=1

By solving Grr and Gtt components of Einstein field equations, we find the following solutions for the metric function R(t) and the cosmological constant Λ27 R(t)=(ηt+ν)a2(N-2)2,

28 Λ=N-1(N-2)2a2η2Na2(N-2)2-1,

where η and ν are arbitrary constants, N is the number of spatial dimensions and a is the coupling constant. We can avoid R(t)=0 by imposing η≥0 and ν>0. The solution (27) for the metric function R(t), indicates that for a fixed number of dimension, the value of R(t) increases by increasing the norm of the coupling constant |a|.

From Eq. (28), we realize that depending on the value of the coupling constant a, and the number of spatial dimension N, the cosmological constant Λ can become positive, negative or zero. For example, in (4+1)-dimensions, Λ>0 (dS space) when the coupling constant a>1, Λ=0 for a=±1, and Λ<0 (AdS space) for -1<a<1. As an example, in Fig. 2, we illustrate the changes in the cosmological constant Λ with respect to the spatial dimensions N for different values of the coupling constant a, and its changes with respect to the coupling constant a for different values of N.

We find the following relation between the coupling constants a and b by solving the other non-zero components of the field equations29 ab=-(N-2),

which indicates that the coupling constants a and b cannot be equal to each other. Moreover, we can rewrite the action of Einstein–Maxwell-dilaton theory (13), as30 S=∫dN+1x-gR-4N-1(∇ϕ)2-e-4/(N-1)aϕF2-e-4(N-2)a(N-1)ϕΛ,

from which we realize that increasing the coupling constant a, leads to an increase in the strength of the interaction between the dilaton field and the cosmological constant, and the decrease in the strength of the interaction between the dilaton field and the electromagnetic field. It is worth noting that our N+1-dimensional results are all independent of the constants k and c that appear in the Bianchi type IX geometry (11). Therefore, our results satisfy all the field equations for any sub-classes of the Bianchi type IX geometry. Substituting the result for the metric functions H(r,θ) and R(t) into the dilaton field Φ(t,r,θ) (18), we get31 Φ(t,r,θ)=N-1-N+2lng+r2cosθ+g-+lnηt+νa2+N-2a2N-2a2+N-2.

We show the behaviour of the dilaton field Φ(t,r,θ) for three different spatial dimensions in Fig. 3. We notice that for a specific time slice, the dilaton field decreases by increasing the dimension.Fig. 3 The behaviour of the dilaton field Φ(t,r,θ) for three different spatial dimensions N=4, N=5 and N=6, which correspond to the upper, middle and lower surface, respectively. We set the constants g+=0.5, g-=15, η=1, ν=2 and a=1, and assume the time slice t=1

In order to study the singularities of this spacetime, we calculate the Ricci scalar and the Kretschmann invariant of the solutions in N+1-dimensions32 RN=f(N)(t,r,θ,ψ)r10R(t)2H(r,θ)sin2θ,

33 KN=g(N)(t,r,θ,ψ)r14R(t)2H(r,θ)sin2θ,

where we show the numerator of these expressions by f(N)(t,r,θ,ψ) and g(N)(t,r,θ,ψ), respectively, which are functions of the coordinates. These functions change based on the number of dimensions. We realize that in any dimension N≥4, the Ricci scalar and the Kretschmann invariant diverge at r=0, sinθ=0 and on the hypersurface H(r,θ)=0. We can avoid the singularities at R(t)=0, by restricting the constants η≥0 and ν>0.

There is a correspondence in asymptotically AdS/dS spacetimes, between the UV/IR physics in dual conformal field theory and the near boundary of the spacetime. The c-theorem states that for an expanding dS spacetime, the renormalization group flows to the ultraviolet, and for a contracting dS spacetime, to the infrared [44–46]. The c-function in N+1-dimensions is given by34 c∼1(Gtt)N-12,

where Gtt is the effective Einstein field tensor. In Fig. 4, we show the behaviour of the c-function in N=4.Fig. 4 The behaviour of the c-function with respect to the time coordinate in N=4, for a set of values for the constants

The electromagnetic gauge field Aμ given in Eq. (19) is explicitly given by35 At(t,r,θ)=α(ηt+ν)a2N-2(g+r2cosθ+g-)-1,

which yields the electric field in r and θ directions. Furnished with (35), we find the components of the electric field are given by36 Er=-rcosθg+(ηt+ν)a2N-22N-2(g+r2cosθ+g-)2a2+N-2,

37 Eθ=r2sinθg+(ηt+ν)a2N-22N-22(g+r2cosθ+g-)2a2+N-2.

We show the behaviour of these electric fields in Fig. 5 for assumed values for the constants.Fig. 5 The electric fields a Er and b Eθ in (4+1)-dimensions at t=5, where we set a=3, η=1, ν=2, g+=1 and g-=15

We show in Sect. 5 that our solutions in this section can be uplifted to a higher dimensional Einstein-form theory, however they can’t be uplifted into a higher dimensional Einstein–Maxwell theory with cosmological constant, or just higher dimensional Einstein gravity with a cosmological constant. In fact, in the first uplifting situation, the uplifting works if a=b, and the number of extra spatial directions is given by [5, 47]38 d=3a21-a2.

However as we noticed before, the consistency of the solutions implies ab=-(N-2) and so, there is no value for the coupling constant a, such that a=b. The latter uplifting, i.e. the case of uplifting the N+1-dimensional solutions to Einstein–Maxwell-dilaton theory with two coupling constants to the solutions of N+2-dimensional gravity with the cosmological constant, is only possible if the coupling constants are equal to a=±2 and b=±12, [48, 49]. However, the product of these coupling constants are equal to +1, and so they are not in agreement with the relation between the coupling constants (29). We actually explicitly check out that the N+2-dimensional metric39 dsN+22=e∓4312ϕ(t,r,θ)dsN+12+e±412ϕ(t,r,θ)(dz+2At(t,r,θ)dt)2,

where dsN+12, ϕ(t,r,θ) and At(t,r,θ) are given by (14), (31), and (35), respectively, does not satisfy the field equations of the N+2-dimensional gravity with a cosmological constant. We note that in (39), z is the uplifted coordinate. In fact, a careful analysis of the uplifting process in [49], shows that the N+1-dimensional metric for the Einstein–Maxwell-dilaton theory, always is diagonal; though our metric ansatz (14) has off-diagonal elements due to the Bianchi type IX geometry. It looks like that we may need a new anstaz for the uplifting of the N+1-dimensional Einstein–Maxwell-dilaton theory into a N+2-dimensional gravity with the cosmological constant.

In the next section, we propose a different set of ansatzes to find a new exact solution to the theory, which includes the case where the coupling constants are equal to each other a=b.

Einstein–Maxwell-dilaton theory with two non-zero equal coupling constants

One of the main results of the previous section is the restriction on the coupling constants ab=-(N-2), which indicates that with the considered ansatzes (14), (18), and (19), the coupling constants can’t be equal to each other. In order to explore a new class of exact solution to the Einstein–Maxwell-dilaton theory with equal coupling constants a=b, we assume a new ansatz for the N+1-dimensional spacetime as40 dsN+12=-1H(t,r,θ)2dt2+H(t,r,θ)2(N-2)R(t)2[dsB.IX2+Σi=1N-4dxi2],

where the metric function H(t,r,θ) is now a function of the spatial coordinates r and θ, as well as the time coordinate t. We also consider a new ansatz for the electromagnetic gauge41 At(t,r,θ)=αRX(t)HY(t,r,θ),

and the dilaton field42 ϕ(t,r,θ)=-(N-1)4aln(HU(t,r,θ)RV(t)),

where X, Y, U and V are constants that will be determined by the field equations.Fig. 6 The metric function H(t,r,θ) in terms of the coordinates r and θ for three different spatial dimensions N=4, N=5 and N=6 (a), which correspond to the upper, middle and lower surface, respectively. The metric function H(t,r,θ) in terms of the coordinates r and θ for N=5 (b), which is the zoom in of middle surface in a. We assumed the constants η=1, ν=2, g+=0.5, g-=5, a=1.5 and the time slice t=2

Combining the electromagnetic field equations Mr and Mt, leads to a differential equation for the metric function H(t,r,θ), which can be solved as43 H(t,r,θ)=R-(N-2)(t)(R(t)p+g+r2cosθ+g-)N-2N-2+a2,

where g± and p are arbitrary constants. Substituting the result for H(t,r,θ) in the other non-zero components of the electromagnetic and Einstein field equations, we determine the constants X, Y and α in the Maxwell gauge field (41)44 X=-a2,Y=-1-a2(N-2),α2=(N-1)2(a2+(N-2)),

and the constants U and V in the dilaton field (42)45 U=2a2(N-2),V=2a2.

In Fig. 6, we represent the behaviour of the metric function H(t,r,θ) with respect to the coordinates r and θ, for three different dimensions, where we set the coupling constant a=1.5. As we notice from Fig. 6a, the fluctuations in the metric function is not visible especially for N=5 and N=6, because of the vertical axis values. In Fig. 6b, we zoom in on the metric function for N=5, and notice non-trivial fluctuations versus the coordinates r and θ. We notice that for a specific time slice, the metric function decreases by increasing the dimension of spacetime.

Substituting these results in the Einstein equation Gtt and Grr, we find the metric function R(t) and the cosmological constant as46 R(t)=(ηt+ν)Δ,

47 Λ=-(N-1)η2(a2-N)(a2+N-2)2a2+2a2+4-N2,

where η and ν are arbitrary constants, and Δ is a constant that has the following relation with p that appears as the power of R(t) in Eq. (43)48 p=(N-2)Δ+1Δ.

In order to reproduce the same results of the reference [50] for the five-dimensional Einstein–Maxwell-dilaton theory, based on the Bianchi type IX geometry, we consider p=a2+2 with no loss of generality. This choice makes Δ=1a2+4-N. We note that the divergence in the Δ can be removed by choosing different values for p.

Substituting the results for the metric functions H(r,θ) and R(t), into the dilaton field Φ(t,r,θ) (42), we get49 Φ(t,r,θ)=-2a(N-1)N-2+a2ln{(ηt+ν)a2+2+g+r2cosθ+g-}.

We show the behaviour of the dilaton field Φ(t,r,θ) for three different spatial dimensions in Fig. 7. As we notice from Fig. 7a, the fluctuations in the dilaton field is not visible for different values of N=4, N=5 and N=6, because of the vertical axis values. In Fig. 7b, we zoom in on the dilaton field for N=5, and notice non-trivial fluctuations versus the coordinates r and θ. We notice that for a specific time slice, the dilaton field increases by increasing the dimension of spacetime.Fig. 7 The dilaton field Φ(t,r,θ) in terms of the coordinates r and θ for three different spatial dimensions N=4, N=5 and N=6 (a), which correspond to the lower, middle and upper surface, respectively. The dilaton field Φ(t,r,θ) in terms of the coordinates r and θ for N=5 (b), which is the zoom in of middle surface in a. We assumed the constants η=1, ν=2, g+=0.5, g-=5, a=1.5 and the time slice t=2

Fig. 8 The behaviour of the cosmological constant Λ in terms of a the number of spatial dimension N for two different values of the coupling constant a, and b in terms of a in two different dimensions, where we assumed η=1

We realize that the cosmological constant Λ in N+1-dimensions can be positive, negative or zero, based on the coupling constant a and the number of the spatial dimension N. In Fig. 8, we show the behaviour of the cosmological constant Λ in terms of the number of dimension N in Fig. 8a, and in terms of the coupling constant a, for different values of N, in Fig. 8b.

It’s worth mentioning that for the case where the coupling constant a=1, in N=4 dimensions, the action of Einstein–Maxwell-dilaton theory reduces to the low-energy effective action for heterotic string theory [51].

In Fig. 9, we represent the c-function of the solutions, where N=4.Fig. 9 The behaviour of the c-function for N=4 for a set of values for the constants

The electromagnetic gauge field Aμ given in Eq. (41), is explicitly given by50 At(t,r,θ)=α(ηt+ν)-a2a2+4-N×{(ηt+ν)N-2+a2a2+4-N((ηt+ν)2+a2a2+4-N+g+r2cosθ+g-)-1},

which yields the electric field in r and θ direction. Furnished with (50), we find the components of the electric field are given by51 Er=-rcosθg+(ηt+ν)2a4-N+2-a2+N-42N-2(1+(g+r2cosθ+g-)(ηt+ν)a4-a2+N-4)2a2+N-2,

52 Eθ=r2sinθg+(ηt+ν)2a4-N+2-a2+N-42N-2(1+(g+r2cosθ+g-)(ηt+ν)a4-a2+N-4)22a2+N-2.

In Fig. 10, we show the electric fields Er and Eθ in terms of the coordinates r and θ, for a set values for the constants.Fig. 10 The electric fields a Er and b Eθ in (4+1)-dimensions for a=b=1, where we set t=10, η=1, ν=2, g+=1 and g-=15

We also note that the N+1-dimensional metric (40) can be uplifted to the solution to the higher dimensional Einstein–Maxwell theory with a cosmological constant, if the coupling constant a is greater than or equal to 12 and less than 1 [49]. The Einstein–Maxwell theory with a cosmological constant is a N+1+D-dimensional theory with the cosmological constant equal to Λ=9μ22(D+3)(D+4)D2, where D=3a21-a2. The N+1+D-dimensional metric is given by53 dsN+1+D2=e43DD+3ϕ(t,r,θ)dsN+12+e-4DDD+3ϕ(t,r,θ)dy→·dy→,

where y→=(y1,…,yD) parametrizes an Euclidean ED space, and dsN+12 and ϕ(t,r,θ) are given by (40) and (49), respectively. We explicitly check that (53) satisfies all the Einstein–Maxwell field equations for a=12,25 and 12 (where D=1,2 and 3) with the cosmological constants Λ=90μ2,1354μ2 and 21μ2, respectively.

We should note that the solution for the dilaton field (42) is not well defined, if the coupling constant a, is equal to zero. In other words, if we are interested in the exact solutions, in the limit where a→0, we shall impose the limit in the action (13), and solve the corresponding field equations. In the next section, we propose a different set of ansatzes to find a new exact solution to the theory, which includes the case where the coupling constants are equal to zero a=b=0.

Einstein–Maxwell-dilaton theory with two zero coupling constants

In this section, we consider the N+1-dimensional Einstein–Maxwell theory in the presence of the cosmological constant. Considering a=b=0 in the last section makes the dilaton field (42) diverges. Moreover, this assumption leads to a divergent cosmological constant (47) when N=4. Therefore, we assume the following ansatz for the line element, and the electromagnetic gauge54 dsN+12=-1H(t,r,θ)2dt2+H(t,r,θ)2(N-2)R(t)2[dsB.IX2+Σi=1N-4dxi2],

55 At(t,r,θ)=αH(t,r,θ),

respectively. Solving the field equations, we find the metric function H(t,r,θ) as56 H(t,r,θ)=1+(g+r2cosθ+g-)RY(t),

where g± are two arbitrary constants, and Y is a constant that we will determine by the field equations. Moreover, solving the Gtt component of Einstein equation, and substituting the result and the form of H(t,r,θ) (56) into Grr, we find a differential equation for R(t), which gives us57 R(t)=(ηexpϵΛX1/2t)p,

where η is an arbitrary constant, ϵ=±1 and X and p will be determined throughout the equations. From Grθ component of Einstein equation58 Grθ=-4r3g+2sinθcosθ(α2-(N-1)2(N-2))(g+r2cosθ+g-+exp{(N-2)ϵ(ΛX)1/2t}),

we find the constant α2 that appears in electromagnetic gauge field (55) as59 α2=(N-1)2(N-2).

Substituting these results into the rest of the equations of motion, we determine the constants Y and X as60 Y=-(N-2)p,X=N(N-1)ϵ2.

Assuming p=1, which makes Y=-(N-2), recovers the same results proposed in [50] for the five-dimensional Einstein–Maxwell-dilaton theory, based on the Bianchi type IX geometry.

In Fig. 11, we show the behaviour of the metric function H(t,r,θ) in 5-dimensional spacetime (N=4) for three different time slices.

The c-function in N+1-dimensions is given by61 c∼1(Gtt)N-12,

where Gtt is the effective Einstein field tensor. For N=4, we show the behaviour of the c-function in terms of the time coordinate, for both ϵ=+1 and ϵ=-1 in Fig. 12.Fig. 11 The metric function H(t,r,θ) in (4+1)-dimensions, where the upper surface, middle surface and the lower surface correspond to specific time slices t=1, t=2 and t=3, respectively. We consider the constants as a=1, g+=0.5, g-=15 and Λ=1

Fig. 12 c-function in terms of time for a ϵ=+1 and b ϵ=-1, for a set of values for the constants

Furnished with our results, we find the electric fields produced by the electromagnetic gauge field (55) as62 Er=-rg+cosθ2N-2exp(N+2)N2-NϵΛ1/2tN-2((g+r2cosθ+g-)exp2N2-NϵΛ1/2t+expNN-1ϵΛ1/2t)2,

63 Eθ=r2g+sinθ2N-2exp(N+2)N2-NϵΛ1/2t2N-2((g+r2cosθ+g-)exp2N2-NϵΛ1/2t+expNN-1ϵΛ1/2t)2.

Fig. 13 The electric fields a Er and b Eθ in (4+1)-dimensions, where we set t=2, Λ=1, ϵ=1, g+=0.5 and g-=15

In Fig. 13, we show the behaviour of the electric fields Er and Eθ in (4+1)-dimensions, where we set t=2 and Λ=1.

Embedding of the solutions in higher-dimensional theory

In this section, we consider the D-dimensional gravity coupled to a form field, in presence of a cosmological constant ΛD [52],64 SD=∫dDx-gR-12(q+2)!F[q+2]2+2ΛD,

where D=p+q+1. The B[q+1] denotes a q+1-potential. We note that in Eq. (64), R is the D-dimensional Ricci scalar, and F[q+2] is given as,65 F[q+2]=dB[q+1].

Note that in (65), F[q+2] is the q+2-field strength form, and dB[q+1] is the exterior derivation of the potential B[q+1].

We now consider the dimensional reduction of the D-dimensional theory (64) to p+1-dimensions. We consider an internal curved q-dimensional space, with the line element dKq2 [52]. We consider the following ansatz for the D-dimensional metric,66 dsD2=e-δϕ′dsp+12+eϕ′2δ(p-1)-δdKq2,

and also the q+1-potential B[q+1],67 B[q+1]=A[1]∧dKq.

It turns out that the theory in p+1 dimensions, is given by [52],68 Sp+1=∫dp+1xR′-12(∇ϕ′)2-14eγϕ′F[2]2+2ΛDe-δϕ′+2Λ′e-2δ(p-1)ϕ′.

We note that in Eq. (68), R′ is the p+1-dimensional Ricci scalar, and Λ′=R′′/2, where R′′ is the q-dimensional Ricci scalar of the internal space. We also note that in the action (68), δ and γ are the dilaton coupling constants, which are given by [47],69 δ=2q(p-1)(p+q-1)1/2,

70 γ=δ(2-p).

We compare Eq. (13) with (68), and find that the dilaton fields are related by71 ϕ′=8N-1ϕ.

We also find that72 ΛD=0,

and73 2Λ′=-Λ,

and74 A[1]=2Atdt.

Moreover, the coupling constants in (68) are given by75 δ=-N-12(p-1)b8N-1,

76 γ=-4aN-1N-18,

in terms of coupling constants a and b in (13). Equations (29), (70), (75) and (76) yield77 p=N.

Moreover from Eq. (69), we find the dimension of the compact space is given by78 q=N-1b2-1.

So, we can conclude that our N+1-dimensional solutions in Sect. 2 with the metric (14), the dilaton field (18), and the electromagnetic field (19) can be uplifted to the D=N+1+N-1b2-1-dimensional theory (64) with no cosmological constant ΛD=0. Of course, Eq. (78) imposes another constraint on the coupling constant b, to have an integer value for q.

Conclusions

In this article, we present new classes of exact solutions to the D-dimensional Einstein–Maxwell-dilaton theory which describes the dynamical black holes in D≥5 dimensions. The solutions are based partially on the Bianchi type IX geometry, where there are two couplings between the dilaton field and the electromagnetic field, as well as the dilaton with the cosmological constant. We consider three different cases where the coupling constants are not equal, are non-zero and equal, and finally both are zero. In each case, we use different ansatzes for the D-dimensional metric, the dilaton and the Maxwell’s field. We analytically solve all the field equations and find unique solutions for the metric functions, the dilaton and electromagnetic fields. By solving the field equations, we find that for the case of not equal coupling constants, there is a constraint on the coupling constants. Moreover, we find that the cosmological constant can take any positive, zero or negative values. We also show that for special values for the coupling constant b, the solutions with two non-equal coupling constants, can be uplifted to a higher dimensional Einstein-form theory with no cosmological constant. Moreover, We show that for special values for the coupling constant a, the solutions with two equal coupling constants, can be uplifted to a higher dimensional Einstein theory with a cosmological constant. We also should mention that the different classes of exact solutions are completely unique and analytical, and we do not use any approximations to find them, at all. We conclude with the observation that the well-known holography between the rotating black holes and the conformal field theories (CFTs) enjoys the independence of the central charges of the CFT on the non-gravitational matter fields [53, 54]. In this article, we found some exact solutions to the five and higher dimensional gravity coupled to non-gravitational fields. It would be an interesting project to find the rotating versions of the exact solutions, presented in the article. The rotating solutions provide a treasure trove of solutions, including black holes, where we can find and study their holographic dual CFTs. Moreover, we can test the independence of the central charges of the CFT on the non-gravitational fields, for a broader class of gravitational theories. One other interesting line of research is to seek the possible hidden symmetry in the solutions space of a probe field, in the background of rotating versions of the exact solutions, presented in the article. These symmetries, in general, lead to finding the possible dual hidden CFT to the black holes [55]. Moreover extending the dual hidden CFT by introducing a deformation parameter in the radial equation of the probe field, as well as finding the different pictures for the dual hidden CFT [56], are some of other applications of the rotating versions of the exact solutions, presented in the article. We leave studying the rotating versions of the exact solutions and their above-mentioned applications in holography for a future article.

Appendix

For the case where the coupling constants are not equal a≠b, we show the Mt component of the electromagnetic field equation in N+1-dimensions, which gives a differential equation for the metric function H(r,θ)79 Mt=-64(H-N+a2N-2α(a2+N-2)(N-2)2r5R4(t)sinθ16k4c4-r416c4-r4(((N-2)sinθ(c4k4-c4)cos2ψ-k4c4+r4/16)Hr3∂2∂θ2H+(N-2)(r9/32-c4(k4+1)r5/2+8c8k4r)sinθH×∂2∂r2H+sinθa2(c4(k4-1)cos2ψ-k4c4+r4/16)r3∂∂θH2-(N-2)cosθH×(c4(k4-1)cos2ψ+c4-r4/16)r3∂∂θH+4sinθ∂∂rH((a2r9/256-a2c4(k4-1)r5/16+a2c8k4r)∂∂rH-(N-2)(-3r8/256+c4(k4-1)r4/16+c8k4)H)).

Solving this equation, we find a solution for H(r,θ)80 H(r,θ)=(g+r2cosθ+g-)N-2a2+N-2,

where g± are arbitrary constants.

Acknowledgements

This work was supported by the Natural Sciences and Engineering Research Council of Canada.

Data Availability Statement

This manuscript has no associated data. [Authors’ comment: Data sharing not applicable to this article as no datasets were generated or analysed during the current study.]

Code Availability Statement

Code/software will be made available on reasonable request. [Authors’ comment: The code/software generated during and/or analysed during the current study is available from the corresponding author on reasonable request.]
==== Refs
References

1. Freund PGO Rubin MA Phys. Lett. B 1980 97 233
P.G.O. Freund, M.A. Rubin, Phys. Lett. B 97, 233 (1980)
2. Th. Kaluza, arXiv preprint arXiv:1803.08616
3. Guica M Hartman T Song W Strominger A Phys. Rev. D 2009 80 124008
M. Guica, T. Hartman, W. Song, A. Strominger, Phys. Rev. D 80, 124008 (2009)
4. Ghezelbash AM Phys. Rev. D 2015 91 084003
A.M. Ghezelbash, Phys. Rev. D 91, 084003 (2015)
5. Goutéraux B Smolic J Smolic M Skenderis K Taylor M JHEP 2012 1 89
B. Goutéraux, J. Smolic, M. Smolic, K. Skenderis, M. Taylor, JHEP 1, 89 (2012)
6. Mahapatra S Roy P JHEP 2018 11 138
S. Mahapatra, P. Roy, JHEP 11, 138 (2018)
7. Stetsko MM Eur. Phys. J. 2019 C79 244
M.M. Stetsko, Eur. Phys. J. C79, 244 (2019)
8. Mahapatra S Priyadahshinee S Reddy G Shukla B Phys. Rev. D 2020 102 024042
S. Mahapatra, S. Priyadahshinee, G. Reddy, B. Shukla, Phys. Rev. D 102, 024042 (2020)
9. P. Goulart, arXiv preprint arXiv:1809.06533
10. Julié F J. Cosmol. Astropart. Phys. 2018 10 033
F. Julié, J. Cosmol. Astropart. Phys. 10, 033 (2018)
11. Li L Phys. Lett. B 2017 767 278
L. Li, Phys. Lett. B 767, 278 (2017)
12. Torii T Shiromizu T Phys. Lett. B 2003 551 161
T. Torii, T. Shiromizu, Phys. Lett. B 551, 161 (2003)
13. Rougemont R Grefa J Hippert M Noronha J Noronha-Hostler J Portillo I Ratti C Prog. Part. Nucl. Phys. 2024 135C 104093
R. Rougemont, J. Grefa, M. Hippert, J. Noronha, J. Noronha-Hostler, I. Portillo, C. Ratti, Prog. Part. Nucl. Phys. 135C, 104093 (2024)
14. Badía J Eiroa EF Phys. Rev. D 2023 107 124028
J. Badía, E.F. Eiroa, Phys. Rev. D 107, 124028 (2023)
15. Heydari-Fard M Heydari-Fard M Sepangi HR Phys. Rev. D 2022 105 124009
M. Heydari-Fard, M. Heydari-Fard, H.R. Sepangi, Phys. Rev. D 105, 124009 (2022)
16. Tripathi A Zhou B Abdikamalov AB Ayzenberg D Bambi C JCAP 2021 2107 002
A. Tripathi, B. Zhou, A.B. Abdikamalov, D. Ayzenberg, C. Bambi, JCAP 2107, 002 (2021)
17. Tan W Phys. Rev. D 2020 102 044054
W. Tan, Phys. Rev. D 102, 044054 (2020)
18. Lü H Mao P Wu J JHEP 2019 11 005
H. Lü, P. Mao, J. Wu, JHEP 11, 005 (2019)
19. Rocha JV Tomašević M Phys. Rev. D 2018 98 104063
J.V. Rocha, M. Tomašević, Phys. Rev. D 98, 104063 (2018)
20. Li S Wei H Phys. Rev. D 2019 99 064002
S. Li, H. Wei, Phys. Rev. D 99, 064002 (2019)
21. Azreg-Aïnou M Ahmed AK Jamil M Class. Quantum Gravity 2018 35 235001
M. Azreg-Aïnou, A.K. Ahmed, M. Jamil, Class. Quantum Gravity 35, 235001 (2018)
22. Khalil M Sennett N Steinhoff J Vines J Buonanno A Phys. Rev. D 2018 98 104010
M. Khalil, N. Sennett, J. Steinhoff, J. Vines, A. Buonanno, Phys. Rev. D 98, 104010 (2018)
23. Brito R Pacilio C Phys. Rev. D 2018 98 104042
R. Brito, C. Pacilio, Phys. Rev. D 98, 104042 (2018)
24. Azreg-Aïnou M Haroon S Jamil M Rizwan M Int. J. Mod. Phys. D 2019 28 1950063
M. Azreg-Aïnou, S. Haroon, M. Jamil, M. Rizwan, Int. J. Mod. Phys. D 28, 1950063 (2019)
25. Pacilio C Phys. Rev. D 2018 98 064055
C. Pacilio, Phys. Rev. D 98, 064055 (2018)
26. Butler M Ghezelbash AM Int. J. Mod. Phys. 2019 A34 1950061
M. Butler, A.M. Ghezelbash, Int. J. Mod. Phys. A34, 1950061 (2019)
27. Pedraza JF Sybesma W Visser MR Class. Quantum Gravity 2019 36 054002
J.F. Pedraza, W. Sybesma, M.R. Visser, Class. Quantum Gravity 36, 054002 (2019)
28. Ishihara H Kimura M Matsuno K Phys. Rev. D 2016 93 024037
H. Ishihara, M. Kimura, K. Matsuno, Phys. Rev. D 93, 024037 (2016)
29. Matsuno K Ishihara H Kimura M Class. Quantum Gravity 2015 32 215008
K. Matsuno, H. Ishihara, M. Kimura, Class. Quantum Gravity 32, 215008 (2015)
30. Kanou Y Ishihara H Kimura M Matsun K Tatsuoka T Phys. Rev. D 2014 90 084004
Y. Kanou, H. Ishihara, M. Kimura, K. Matsun, T. Tatsuoka, Phys. Rev. D 90, 084004 (2014)
31. Ida D Ishihara H Kimura M Matsuno K Morisawa Y Tomizawa S Class. Quantum Gravity 2007 24 3141
D. Ida, H. Ishihara, M. Kimura, K. Matsuno, Y. Morisawa, S. Tomizawa, Class. Quantum Gravity 24, 3141 (2007)
32. Abbott BP Phys. Rev. Lett. 2016 116 061102 26918975
B.P. Abbott et al., Phys. Rev. Lett. 116, 061102 (2016)26918975
33. Kastor D Traschen JH Phys. Rev. D 1993 47 5370
D. Kastor, J.H. Traschen, Phys. Rev. D 47, 5370 (1993)
34. London LAJ Nucl. Phys. B 1995 434 709
L.A.J. London, Nucl. Phys. B 434, 709 (1995)
35. Ishihara H Kimura M Tomizawa S Class. Quantum Gravity 2006 23 L89
H. Ishihara, M. Kimura, S. Tomizawa, Class. Quantum Gravity 23, L89 (2006)
36. Yoo CM Ishihara H Kimura M Matsuno K Tomizawa S Class. Quantum Gravity 2008 25 095017
C.M. Yoo, H. Ishihara, M. Kimura, K. Matsuno, S. Tomizawa, Class. Quantum Gravity 25, 095017 (2008)
37. Ghezelbash AM Phys. Rev. D 2008 78 126002
A.M. Ghezelbash, Phys. Rev. D 78, 126002 (2008)
38. Meyer KR Am. Math. Mon. 2001 108 729
K.R. Meyer, Am. Math. Mon. 108, 729 (2001)
39. Bali R Dave S Pramana 2001 56 513
R. Bali, S. Dave, Pramana 56, 513 (2001)
40. Ghezelbash AM Phys. Rev. D 2006 74 126004
A.M. Ghezelbash, Phys. Rev. D 74, 126004 (2006)
41. Ashtekar A Wilson-Ewing E Phys. Rev. D 2009 80 123532
A. Ashtekar, E. Wilson-Ewing, Phys. Rev. D 80, 123532 (2009)
42. Starobinsky AA Sushkov SV Volkov MS Phys. Rev. D 2020 101 064039
A.A. Starobinsky, S.V. Sushkov, M.S. Volkov, Phys. Rev. D 101, 064039 (2020)
43. Maki T Shiraishi K Class. Quantum Gravity 1993 10 2171
T. Maki, K. Shiraishi, Class. Quantum Gravity 10, 2171 (1993)
44. Ghezelbash AM Phys. Rev. D 2010 81 044027
A.M. Ghezelbash, Phys. Rev. D 81, 044027 (2010)
45. Strominger A JHEP 2001 10 034
A. Strominger, JHEP 10, 034 (2001)
46. Balasubramanian V De Boer J Minic D Phys. Rev. D 2002 65 123508
V. Balasubramanian, J. De Boer, D. Minic, Phys. Rev. D 65, 123508 (2002)
47. Goutéraux B Kiritsis E JHEP 2011 12 036
B. Goutéraux, E. Kiritsis, JHEP 12, 036 (2011)
48. Charmousis C Gouteraux B Soda J Phys. Rev. D 2009 80 024028
C. Charmousis, B. Gouteraux, J. Soda, Phys. Rev. D 80, 024028 (2009)
49. Charmousis C Langlois D Steer D Zegers R JHEP 2007 0702 064
C. Charmousis, D. Langlois, D. Steer, R. Zegers, JHEP 0702, 064 (2007)
50. Fahim BH Ghezelbash AM Eur. Phys. J. C 2021 81 1
B.H. Fahim, A.M. Ghezelbash, Eur. Phys. J. C 81, 1 (2021)
51. Rocha JV Tomašević M Phys. Rev. D 2018 98 104063
J.V. Rocha, M. Tomašević, Phys. Rev. D 98, 104063 (2018)
52. Ghezelbash AM Phys. Rev. D 2017 95 064030
A.M. Ghezelbash, Phys. Rev. D 95, 064030 (2017)
53. Compere G Murata K Nishioka T JHEP 2009 09 077
G. Compere, K. Murata, T. Nishioka, JHEP 09, 077 (2009)
54. Ghezelbash AM JHEP 2009 09 045
A.M. Ghezelbash, JHEP 09, 045 (2009)
55. Castro A Maloney A Strominger A Phys. Rev. D 2010 82 024008
A. Castro, A. Maloney, A. Strominger, Phys. Rev. D 82, 024008 (2010)
56. Ghezelbash AM Siahaan H Gen. Relat. Gravit. 2014 46 1783
A.M. Ghezelbash, H. Siahaan, Gen. Relat. Gravit. 46, 1783 (2014)
