Cosmological Solutions with Charged Black Holes

We consider the problem of constructing cosmological solutions of the Einstein-Maxwell equations that contain multiple charged black holes. By considering the field equations as a set of constraint and evolution equations, we construct exact initial data for N charged black holes on a hypersphere. This corresponds to the maximum of expansion of a cosmological solution, and provides sufficient information for a unique evolution. We then consider the specific example of a universe that contains eight charged black holes, and show that the existence of non-zero electric charge reduces the scale of the cosmological region of the space. These solutions generalize the Majumdar-Papapetrou solutions away from the extremal limit of charged black holes, and provide what we believe to be some of the first relativistic calculations of the effects of electric charge on cosmological backreaction.


Introduction
Inhomogeneous cosmological models are an important tool for investigating the foundational assumptions that go into the standard model of cosmology. One such assumption is that all space-times that are statistically homogeneous and isotropic on large scales should evolve in precisely the same way as the perfectly homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker (FLRW) solutions of Einstein's equations. However, it has never been proven that this assumption should be true under general circumstances, and there are in fact good reasons to question its veracity. This is because averaging and evolution operations do not commute in Einstein's theory, which opens the possibility that small-scale inhomogeneities could affect the large-scale dynamics of the Universe [1][2][3]. Inhomogeneous cosmologies that do not rely on averaging are ideal for investigating this backreaction question, and in this paper we will extend a class of such models to include massive bodies with electric charge.
There are several motivations for including electric charge in these models. The first is that it offers a generalisation of existing models, analogous to the generalisation from Schwarzschild to Reissner-Nordström in the study of static black holes. This should be expected to give new mathematical and physical insights into the nature of these inhomogeneous cosmologies. Moreover, in an FLRW universe the homogeneity of space forbids regions of charged spacetime by definition. There is also relatively little known about the cosmological consequences of electric charge in inhomogeneous models. It therefore remains a somewhat open question as to what the effects of allowing electric charge in cosmology should be in general. The lattice models studied in refs.  allow us the opportunity to study this problem in a well-defined and precise way. Moreover, this work offers the opportunity to generalise the well-known Majumdar-Papapetrou and Kastor-Traschen solutions for multiple extremal black holes with |q| = m [35][36][37] to the case where the black holes have |q| = m. Finally, the inclusion of electric charge offers a way to break the discrete rota-tional symmetries that might otherwise exist within these lattice cosmologies, as we find that the black holes must necessarily contain non-identical charges.
This rest of this paper is organised as follows: in Section 2 we review the Einstein-Maxwell constraint equations on a time-symmetric hypersurface. Then in Section 3 we derive initial data for a cosmological model that contains electrically charged black holes. The example of a universe that contains eight black holes is considered in Section 4, which we then compare with a Friedmann cosmology before calculating the apparent horizons of the black holes. Finally we conclude in Section 5. We find that as the magnitude of the electric charge is increased, the discrepancy between the uncharged lattice models and the corresponding FLRW solutions reduces, and that the horizon size of each black hole decreases. Throughout this paper we adopt geometrised units where 8πG = c = 1.

Geometrostatics with an electric field
The Einstein-Maxwell equations that govern electromagnetic fields in the presence of gravitation are given by where R µν is the Ricci tensor, and R is the Ricci scalar. The Faraday tensor is denoted by F µν , and obeys the following differential relations: where the indices µ, ν, ... run from 0 to 3. It is these equations that we wish to solve, to obtain the geometry of a universe that contains N electrically-charged black holes. To find the relevant initial data for this system, we start by choosing a space-like hypersurface with metric γ ij and normal u µ , where i, j, ... run from 1 to 3. The four constraints from equation (1) can then be found to be where R (3) is the Ricci curvature scalar of the intial 3-dimensional space, K ij is the extrinsic curvature of this 3-space in the 4-dimensional space-time, and where ρ ≡ T µν u µ u ν and q i ≡ −δ j i u µ T µj are the energy and momentum densities, respectively. All indices in these expressions should be understood to be raised and lowered with the metric of the 3-space γ ij , and the covariant derivative symbolized with a colon is with respect to γ ij .
If our initial hypersurface and matter fields are time-symmetric then we immediately have K ij = 0, while F iµ = −E i u µ and F ij = 0. This means that ρ = E i E i and q i = 0, so that the energy constraint equation (3) becomes and the momentum constraint equation (4) is identically satisfied. Simultaneously, the differential relations from equation (2) imply If equations (5) and (6) are both satisfied, then we have a solution to the constraint equations. This provides the initial data for a unique evolution under the six remaining Einstein-Maxwell equations in (1) and the three evolution equations in (2) [38].
Let us now make the following ansätze: firstly, let us write the metric as where ψ and χ are functions of all spatial coordinates, and where h ij is the metric of a conformal 3-dimensional space of curvature R = constant. Secondly, let us choose, as in [31], The time-symmetric constraint equations (5) and (6) then becomē where∇ 2 indicates a covariant Laplacian with respect to the metric h ij . These are two copies of Helmholtz equation, for the two functions χ and ψ. The solutions to this equation are well known, and it can be seen that E i = 0 as long as χ and ψ are not directly proportional. The reader may note that for periodic lattices of black holes the equations in (9) can have positive energy solutions only if R is positive (see [23] for proof).

A universe full of charged black holes
To find cosmologically interesting configurations of black holes, it is convenient to choose h ij to be the metric of a 3-sphere: so that R = 6. By noting that equations (9) are both linear, and by recognising that f ∝ 1/ sin(r/2) is a solution, it can be seen that we can write solutions to these equations as follows: where i runs from 1 to N , and where c i and d i are two sets of N constants (yet to be determined). These solutions contain N poles on the conformal 3-sphere, located at arbitrarily chosen locations, and the symbol r i is intended to denote the value of the r coordinate after rotating the sphere so that the ith pole is located at r i = 0. Each of these poles will correspond to an electrically-charged mass point, at the centre of a black hole in a cosmological model.

A slice through Reissner-Nordström
In principle, equations (11) and (12) represent exact initial data for N electrically charged black holes in a hyperspherical cosmology. In practise, we must relate the constants c i and d i to the charge and mass of each of the black holes, in order for the geometry to be fully specified. To make this connection, it is useful to recall that the Reissner-Nordström solution with mass m and charge q can be written as By making the transformation we can therefore write a time-symmetric slice of the metric from equation (13) as A further transformation, to cast this metric in the form of equations (7), (11) and (12), is given by ρ = k tan(r/2), where k is a constant. This transforms the conformal 3-space from a plane to a sphere, and results in where ds 2 is the line-element on a 3-sphere, as given in equation (10). Finally, if we choose 2k = √ m − q √ m + q, then we get This result is of exactly the form that we would have obtained from considering equations (7), (11) and (12) with one mass positioned at r = 0 and a second mass positioned at r = π. If we further choose c 1 = d 2 = √ m + q and c 2 = d 1 = √ m − q then the two geometries are formally identical. This means that the time-symmetric initial data for the Reissner-Nordström solution can be considered as a special case of our more general intial data, with one mass of charge +q at r = 0 and a second mass of charge −q at r = π.

N arbitrarily positioned, charged black holes
Let us now consider the general case, where N masses are positioned at arbitrarily chosen locations. In the vicinity of a mass point, and by using equation (16), we can relate c i and d i to the mass m i of each of the black holes in any arbitrary distribution of points. The first step in this is rotating coordinates so that the ith mass appears at r = 0. In the limit when r → 0 we can then expand the terms in equation (7) to find where r ij is the coordinate distance between the positions of the ith and jth masses. Comparing this equation with the expanded version of equation (16), again around r = 0, allows us to read off the following from the coefficients of the leading-order terms: and The former of these equations simply relates the coordinate systems of the Reissner-Nordström solution and our multi-black hole solution in the vicinity of one of our mass points. The latter, however, allows the mass of each of our black holes to be identified. That is, if a multi-black hole solution is specified (with a full set of c i s and d i s), then the above analysis shows that the geometry of space in the vicinity of any one of the black holes will be similar to a Reissner-Nordström black hole with mass m, where m is given by equation (20).
It now remains to identify the electric charge on each of our black holes. A general definition of electric charge within a region Ω can be given by [31] where E i is the electric field and n i is the unit inward-pointing normal. For our black holes, it is convenient to take the boundary ∂Ω to correspond to asymptotic infinity, on the far side of the Einstein-Rosen bridge. This gives and dS = ψ 2 χ 2 r 2 sin θdθdφ, where in both of these expressions χ and ψ should be taken to be evaluated in the limit r → 0. Evaluating equation (21) in this limit, and using equations (11), (12) and (22), gives This equation gives the charge on the mass located at r i = 0, and can equally be used to evaluate the charge on every other black hole in our solution. It is interesting to note that equation (23) immediately gives i q i = 0, meaning that the total charge in the universe must be zero, independently of how the black holes are distributed and their masses. This makes sense physically, as lines of flux can only end on masses in a closed cosmology. Adding and subtracting the equations (20) and (23) gives a somewhat simpler pair of equations: While relatively simple to write down in this form, this is still a system of non-linear equations, which makes it difficult to solve for the {c i } and {d i } directly, after the desired {m i } and {q i } have been specified.

A periodic universe with eight charged black holes
Probably the simplest solution to equations (24) and (25) is the case where all but one of the black holes is extremal, with m i = q i . If this is the case, and the one exceptional black hole is labelled i = 1, then vanishing total charge means that we must have q 1 = − ı =1 q i = − ı =1 m i . The only solutions to equations (24) and (25) then have m 1 = −q 1 , and we have the initial data for a Majumdar-Papapetrou [35,36] solution with spatial infinity transformed into the black hole with i = 1. In the next section we will present an example with non-extremal black holes, in which an exact solution for {c i } and {d i } can be found once the mass and charge of each of the black holes has been specified.

Determining the sets of constants {c i } and {d i }
To find a regular arrangement of black holes in a closed cosmological model, we can choose to tile the conformal hypersphere from our 3-dimensional initial data with eight cubic lattice cells. This is one of the six possible tilings of a 3-sphere with regular convex polyhedra, and corresponds to the structure one would obtain by placing a hypercube within the hypersphere in a 4-dimensional Euclidean embedding space, and by ensuring that the vertices of the hypercube are all touching the hypersphere simultaneously (i.e. by circumscribing the cube with the sphere). The lines that connect the points where these two structures touch then form the edges of eight equal-sized cubes, which can be used as the primitive lattice cells of our tiling. If we place one mass point at the centre of each of these cells, then we have a completely regular lattice in which each mass is exactly equidistant from each of its nearest neighbours. One of the nice features of the regular eight-mass universe just described is that each of the black holes will be antipodal to another black hole, just as in the time-symmetric slice through the Reissner-Nordström solution discussed in the previous section. This can be verified by looking at the final column of Table 1, which gives the coordinates of the location of each mass in both a set of Cartesian coordinates in the 4-dimensional Euclidean embedding space, as well as in a set of hyperspherical polar coodinates intrinsic to the 3-space itself. Taking a hint from the existence of the time-symmetric Reissner-Nordström geometry, we can now choose the {c i } and {d i } of these eight masses such that the c i from each mass is equal to the d i of the antipodal mass, i.e. so that where e and f are constants (yet to be determined). This choice then very conveniently allows us to determine from equations (24) and (25) that so that every black hole in the universe has an identical mass to all of the others, while having an equal and opposite electric charge to its antipodal partner. This seems to be the simplest way to satisfy the requirement that the total electric charge of all black holes must vanish, whilst maintaining at least some degree of regularity.
mass number, i (w, x, y, z) (r, θ, φ) Finally, it is now straightforward to solve equations (28) and (29) for e and f to get and f = −q/2(1 + 3 √ 2)e, where m = m i is the mass of each black hole and |q| = |q i | is the magnitude of the charge on each black hole (with sign chosen so that q = q 1 ). We now have an explicit solution where all black holes have equal mass, and half of them have positive charge while the other half have negative charge.

Comparison with Friedmann cosmology
With a cosmological solution in hand, it is interesting to check whether the large-scale properties of our space bear any resemblance to the predictions of the commonly used homogeneous and isotropic FLRW models. The first thing that one could compare, between these two cases, is the scale of the cosmological region at the maximum of expansion. The initial data constructed above can immediately be seen to be at such a maximum, as it was chosen to be time-symmetric. The scale factor at the corresponding moment of time in a spatially closed and dust dominated FLRW model, which we take to be the most similar to our configuration of black holes, is given by where M is the total mass in dust in the entire space-time. By comparing this number to the scale of the cosmological region in our black hole space-time we can then deduce the effect of condensing the mass into a finite number of points, with discretely distributed electric charges. This will extend previous results concerning cosmologies filled with uncharged black holes, and will allow us to consider the cosmological effects of electrical charge. The most natural way to determine the scale of the cosmological region of the black hole space-time is to calculate the proper length of one of the edges of one of our lattice cells. If we rotate our solution until a cell edge lies along a curve with θ = constant and φ = constant, then the proper length of the edge is given by where r 1 and r 2 are the coordinates of the end-points of the edge being studied (i.e. the location of the vertices of the lattice cell). The proper length of a curve that subtends the same angle in a closed FLRW model is given by where a(t max ) is the scale factor from equation (31). In general, the value of d will depend on both the mass of the black holes and their charge (m and q), while the value of d FLRW will depend only on the total mass in the space M (once the locations of the cell vertices have been chosen). If we take M = 8m, so that each cosmology has the same total proper mass, then the relevant comparison of scales in the two cosmologies will be given by where d and d FLRW are given by equations (32) and (33). This ratio is only a function of the charge to mass ratio of the black holes q/m, and reduces to the measure of cosmological back-reaction studied previously in ref. [16] in the limit q → 0. The reader may note that this comparison between FLRW and our black hole lattice neglects the contribution of interaction energies between black holes, which may be considerable. A study of the effect of these quantities on the global cosmology will appear elsewhere [39]. The value of the ratio of scale factors as a function of q/m is straightforward to calculate, and is displayed graphically in figure 1. When q/m → 0 we recover the results of the uncharged case, and as q/m → ±1 the black holes become extremal. It can be seen that the difference from the predicted scale factor from FLRW cosmology is greatest when the black holes are uncharged. Increasing the charges on the black holes decreases the discrepancy with FLRW, but the scale of the cosmological region in our black hole cosmology is always greater than that of the corresponding FLRW model, even when the black holes approach extremality. The reader may note that there are two curves in figure 1. These correspond to the situations where the three cells that meet along the cell edge under investigation all contain black holes that have the same charge (orange line), or when two of the cells have the same charge and the other has the opposite charge (blue line). The results do not depend on whether the same charges in these two situations are positive or negative; the ratio of scale factors is the same in either case. These two curves therefore describe every cell edge that exists within our eight-black hole model. As well as the scale of the cosmological region in our black hole universe, it is also of interest to determine the location of the horizons around each of the black holes. This is partly to ensure that there is no overlap in the horizons corresponding to different black holes, so that the solution can reasonably be referred to as a "cosmological model". For our purposes, the location of the event horizon is impossible to locate, as we do not have knowledge of the global space-time. The position of the apparent horizon, however, can be determined from knowledge of the initial data alone. We will therefore use this surface to approximate the location of the horizon of each of our black holes.
The apparent horizon is the outermost marginally outer trapped surface, by definition. As our initial data is time-symmetric, this surface must be an extremal surface in this 3-space [40]. We can approximate its location by ro-  Fig. 2 The position of the apparent horizon around one of our black holes, r min , displayed as a fraction of the coordinate distance to the halfway point between black holes, r half . We find r min → 0 as q/m → ±1.
tating our 3-sphere until one of the masses appears at r = 0. The area of any sphere centered on this mass must then be given by where r = constant. Integrating this quantity numerically, and determining the value of r that minimizes it, then gives us a good approximation to the location of the black hole horizon. This method of assuming a sphere of constant coordinate radius to approximate the apparent horizon necessarily produces a slight overestimate of the horizon area, as the true area of the apparent horizon is that of a minimal sphere in the initial geometry. However, as long as the horizons are well separated, this error is expected to be very small [17]. The results of calculating the position of the apparent horizon in this way are shown in figure 2, where the r coordinate of the apparent horizon as a fraction of the coordinate distance between neighbouring black holes is displayed. When q/m → 0 it can be seen that the horizon extends about 27% of the way to the halfway point between black holes, as expected from previous work. As the charge on the black hole is increased, to either positive or negative values, the horizon then shrinks back towards the centre of the black hole, until it reaches zero in the limit in which q/m → ±1. This is similar behaviour to what should be expected from a maximal slice through a Reissner-Nordström black hole, and verifies that the black hole horizons in our eight-mass model never touch.

Conclusion
We have presented an analysis of an inhomogeneous cosmological model that consists of a lattice of regularly arranged, charged black holes at maximum of expansion. This generalises existing work in the literature and provides insight into the effect that electric charge has on cosmological physics. We have found that the only solutions that exist are those that have zero net charge. We then applied this analysis to a specific model containing eight black holes, and found a relatively simply way of obtaining an uncharged universe overall -by demanding that black holes exist in antipodal pairs, and that each black hole has equal and opposite charge to its antipodal partner.
We then compared this model to an FLRW model with positive spatial curvature containing the same total mass in dust. We did this by comparing the scale between the two types of universe at maximum of expansion as a function of the charge to mass ratio of the black holes. We found that as the magnitude of the charge on the black holes was increased, the discrepancy between lattice and FLRW was reduced, with the largest discrepancy corresponding to the case where there is no charge. Lastly we investigated the size of the apparent horizon in this initial data, as a function of this charge to mass ratio. We found that as the magnitude of the charge of the black holes is increased the size of the apparent horizon decreases, and that in the limit where the black hole charges become extremal the apparent horizons recede to zero radius.