text:
# Contents

## 1 Introduction

The properties of higher-dimensional black holes, and more generally spacetimes with event horizons, have come to play an increasingly important role in physics. This is not only for the purely theoretical reason that they may throw light on some hitherto intractable problems of black holes in 3+1 spacetime dimensions, but also because, if current ideas about large extra dimensions are correct, then such black holes may possibly be created by high energy collisions in accelerator experiments, and their behaviour might be accessible to direct observation . As well as having possible applications to laboratory scale physics, higher dimensional black holes, and other types of horizons such as cosmological event horizons , may also have played an important role in the early universe.

Many of the properties of black holes and event horizons in higher dimensions are very similar to their counterparts in 3+1 dimensions. For example, the analogue of the spherically symmetric Schwarzschild black hole exists in all dimensions, with the 2-sphere of the four-dimensional solution replaced by a $`(D2)`$ sphere in $`D`$ dimensions. Likewise, there is an obvious higher-dimensional analogue of the usual four-dimensional de Sitter spacetime. Moreover, subject to the strict condition of asymptotic flatness, the former and their charged versions are unique . However, if one drops the condition of strict asymptotic flatness, by allowing other compact Einstein metrics in place of the usual round sphere in the spatial sections at constant radius, then in higher dimensions there are many more possibilities for black hole solutions, even on manifolds with the same topology as the higher dimensional Schwarzschild solution. This is because of the remarkable fact, discovered by Bohm , that for $`5d9`$, the sphere $`S^d`$ carries infinitely many other inhomogeneous Einstein metrics, in addition to its usual round metric. One might wonder whether the resulting black hole solutions in spacetime dimensions $`7D11`$ could arise during scattering processes. This depends upon their stability. In this paper, we find evidence, and in some cases proofs, that they are in fact unstable. To do so we use methods developed in , which showed that the stability depends on the non-negativity of the spectrum of the operator

$$\mathrm{\Delta }_{\mathrm{stab}}=\mathrm{\Delta }_L+\frac{\mathrm{\Lambda }}{d1}\left(4\frac{(5d)^2}{4}\right),$$
(1)
where $`\mathrm{\Delta }_L`$ is the Lichnerowicz Laplacian on transverse traceless second rank symmetric tensor fields on the $`d`$-dimensional compact Einstein space, which satisfies $`R_{ab}=\mathrm{\Lambda }g_{ab}`$. We obtain numerical results establishing the existence of negative Lichnerowicz modes in Bohm metrics on $`S^5`$, hence demonstrating the instability of seven-dimensional black holes constructed using these metrics. We present general arguments suggesting that all the other Bohm sphere metrics will have negative Lichnerowicz modes too.

Bohm also showed the existence of infinitely many inhomogeneous Einstein metrics on the products of spheres $`S^{N_1}\times S^{N_2}`$, for $`4N_1+N_29`$, with $`N_12`$ and $`N_22`$. The lowest dimensional such examples are on $`S^3\times S^2`$. We have also investigated the associated topologically non-trivial black holes (including the homogeneous product metric), and we prove that these are unstable. It is also known that $`S^3\times S^2`$ admits infinitely many other homogeneous Einstein metrics, the so called $`T^{p,q}Spin(4)/U(1)`$ spaces. Of these, only $`T^{1,1}`$ admits Killing spinors, and in fact only the black hole associated with $`T^{1,1}`$ is stable . There are in addition some other inhomogeneous metrics on $`S^3\times S^2`$ with Killing spinors; these are examples of Einstein-Sasaki metrics. We show here that the associated black holes using these metrics are stable. One might think therefore that an unstable black hole based on the usual product $`S^3\times S^2`$ metric would evolve dynamically into one based on one of the $`S^3\times S^2`$ Einstein-Sasaki spaces. Presumably however, in doing so the area of the event horizon must increase. Now at fixed “mass parameter” we show the areas of the horizons of the Einstein-Sasaki horizons with topology $`S^3\times S^2`$ are less than those of the $`S^3\times S^2`$ Bohm metrics. Thus in the evolution, the mass parameter would have to increase, which seems rather paradoxical. Another way to say this is that at fixed temperature (which is proportional to the inverse of the mass parameter) the entropy of the Einstein-Sasaki metrics is less than that of the Bohm metrics.

In addition to constructing black holes, one may use Bohm’s metrics to obtain static inhomogeneous Lorentzian solutions of the Einstein equations with a positive cosmological constant, which are topologically the same as the static de Sitter metric. Like de Sitter spacetime, they contain a cosmological event horizon . As such they provide counterexamples in $`D9`$ spacetime dimensions to the long standing and hitherto intractable Cosmic Baldness conjecture , which would be a generalisation of Israel’s uniqueness theorem to cover the de Sitter situation. The area of the cosmological event horizon in the Lorentzian Bohm metrics is smaller than that in de Sitter spacetime. We believe therefore that they are unstable, and that under a small perturbation they would evolve, at least within the event horizon of a given observer, to a static de Sitter-like state. If this is the case, then although evading the strict letter of the Cosmic Baldness conjecture they would respect the spirit of the weaker No-Hair conjecture,. The latter asserts that apart from unstable cases of measure zero, the generic solution should settle down to a de Sitter-like state within the horizon of any given observer. This is all that is needed to justify the usual intuition behind inflationary models of the early universe. Another possible application for these generalised de Sitter spacetimes would be as tunnelling metrics.

The plan of this paper is as follows. In section 2 we review the link established in  between black-hole stability and the spectrum of the Lichnerowicz Laplacian on the compact $`d`$-dimensional Einstein space $`M_d`$ that forms the constant-radius spatial sections. We also show that this stability criterion is identical, for all dimensions of $`M_d`$, to that for the stability of Freund-Rubin type AdS$`{}_{n}{}^{}\times M_d`$ solutions of gravity coupled to a $`d`$-form. We then discuss a lower bound on the spectrum of the Lichnerowicz operator, based on considerations of the Weyl curvature, that was considered in .

In section 3 we give a detailed discussion of the Bohm metrics, and we exhibit negative modes of the Lichnerowicz operator in some of these backgrounds. In certain cases, including all the Bohm metrics on $`S^3\times S^2`$ and on $`S^3\times S^3`$, we obtain an analytic proof of the existence of negative modes. Intuition leads one to expect negative modes in all the Bohm metrics, and we back this up with some numerical results in certain examples where an analytic proof is lacking. There are also non-compact examples of Bohm metrics, which are Ricci-flat. We give analytic proofs that the non-compact Bohm metrics on $`^3\times S^2`$ and $`^3\times S^3`$, recently considered by Kol , have negative modes of the Lichnerowicz Laplacian. Again, intuition leads one to expect negative modes for all the non-compact Bohm examples.

In section 4 we discuss Einstein-Sasaki metrics, which may be defined as odd-dimensional Einstein metrics $`ds^2`$ whose cone $`d\widehat{s}^2=dr^2+r^2ds^2`$ is Ricci-flat and Kähler. They admit Killing spinors, and we use this fact to obtain a lower bound on the spectrum of the Lichnerowicz operator. In particular we use this to demonstrate that the associated black holes are always stable. Likewise, this establishes that Freund-Rubin compactifications using Einstein-Sasaki manifolds will always be stable.

Section 5 contains a description of the analytic continuation of the Bohm metrics to give Lorentzian spacetimes that are generalisations of de Sitter spacetime. These provide counterexamples to the Cosmic Baldness conjecture. In section 6 we try to relate the existence of negative modes for the Lichnerowicz Laplacian to the non-uniqueness of the Dirichlet problem for the Einstein equations. Section 7 gives our main conclusions, and points to some other applications of Bohm metrics. For example, they can provide magnetic monopole solutions in Kaluza-Klein theory. An appendix gives further details about our numerical techniques, and includes some graphs illustrating the behaviour of the metric functions in the Bohm solutions.

## 2 Stability and the Lichnerowicz Laplacian

Our first aim will be to study the classical stability of two types of spacetime constructed using a general positive curvature Einstein metric $`M_d`$. The first of these comprises generalisations to higher dimensions of the four-dimensional Schwarzschild black hole, in which the spatial 2-sphere at constant radius is generalised to a higher-dimensional Einstein space $`M_d`$. The second class of examples comprises Freund-Rubin type solutions AdS$`{}_{n}{}^{}\times M_d`$ to a theory of Einstein gravity coupled to a $`d`$-form field strength. As we shall show below, the classical stability criteria for both of these classes of spacetimes are expressible as the same criterion on the spectrum of the Lichnerowicz operator acting on transverse traceless symmetric 2-index tensors on $`M_d`$. To set the stage for this discussion, we begin in section 2.1 with a general discussion of the Lichnerowicz operator, reviewing the manner in which it arises from a consideration of the second variation of the Einstein-Hilbert action. We then review the black-hole stability  and AdS$`{}_{n}{}^{}\times M_d`$ stability  criteria in sections 2.2 and 2.3, and in section 2.4 we review an argument given in  which shows how a lower bound on the Lichnerowicz spectrum can be obtained by considering the eigenvalues of the Weyl tensor.

### 2.1 Stability criteria

We begin by considering the Einstein-Hilbert action in $`d`$ dimensions

$$S=_M\sqrt{g}d^dx(R(d2)\mathrm{\Lambda }),$$
(2)
whose Euler-Lagrange equations give the Einstein equation

$$R_{ab}=\mathrm{\Lambda }g_{ab}.$$
(3)
Under the perturbation

$$g_{ab}g_{ab}+h_{ab},$$
(4)
one finds that up to quadratic order in $`h_{ab}`$, the action $`S`$ is given on-shell by $`S=S_0+S_1+S_2+\mathrm{}`$, with

$`S_0=2\mathrm{\Lambda }{\displaystyle _M}\sqrt{g}d^dx,S_1=0,`$
$`S_2={\displaystyle _M}\sqrt{g}d^dx(\frac{1}{4}h^{ab}\mathrm{\Delta }_2h_{ab}+\frac{1}{4}h\mathrm{\Delta }_0h+\frac{1}{2}(_ah^{ab})^2),`$ (5)
where $`hh_a^a`$, $`\text{    }\text{ }^a_a`$ and

$$\mathrm{\Delta }_0h\text{    }\text{ }h+\frac{1}{2}(d2)\mathrm{\Lambda }h,\mathrm{\Delta }_2\mathrm{\Delta }_L2\mathrm{\Lambda }.$$
(6)
Here $`\mathrm{\Delta }_L`$ is the Lichnerowicz Laplacian operator acting on symmetric rank two tensors,

$$\mathrm{\Delta }_Lh_{ab}\text{    }\text{ }h_{ab}2R_{acbd}h^{cd}+R_{ca}h^c{}_{b}{}^{}+R_{cb}h^c{}_{a}{}^{}.$$
(7)
If we consider a transverse traceless perturbation,

$$^ah_{ab}=0,h^a{}_{a}{}^{}=0,$$
(8)
then the second variation of the action is simply given by

$$S_2=\frac{1}{4}_M\sqrt{g}d^dxh^{ab}\mathrm{\Delta }_2h_{ab}.$$
(9)
In a perturbative classical stability analysis one asks whether there are finite energy solutions to the linearised equations of motion

$$\mathrm{\Delta }_2h_{ab}=0$$
(10)
that grow exponentially in time. This question has been studied in some depth recently for generalised vacuum black hole spacetimes in $`D=d+2`$ dimensions  and also for spacetimes that are a direct product of $`Dd`$ dimensional anti-de Sitter with a $`d`$ dimensional compact Einstein manifold . The anti-de Sitter spacetimes are supported by a gauge field of appropriate rank.

### 2.2 Generalised Schwarzschild-Tangherlini spacetimes

Generalised Schwarzschild-Tangherlini black holes have the form

$$d\widehat{s}^2=\left[1\left(\frac{\mathrm{}}{r}\right)^{d1}\right]dt^2+\frac{dr^2}{\left[1\left(\frac{\mathrm{}}{r}\right)^{d1}\right]}+r^2ds_d^2,$$
(11)
where $`\mathrm{}`$ is a constant and $`ds_d^2`$ is the metric on a $`d`$ dimensional compact Einstein manifold $`B`$ with the curvature normalised to be that of $`S^d`$

$$R_{\alpha \beta }=(d1)g_{\alpha \beta }.$$
(12)
The black hole solution has vanishing cosmological constant.

It was found in  that the dangerous mode for instability is a transverse tracefree eigenfunction of the Lichnerowicz Laplacian on the Einstein manifold $`B`$.

$`\widehat{h}_{0a}=\widehat{h}_{1a}=0,`$
$`\widehat{h}_{\alpha \beta }=h_{\alpha \beta }(x)r^2\varphi (r)e^{\mathrm{i}\omega t},`$ (13)
where $`x`$ are coordinates on $`B`$ and

$$\mathrm{\Delta }_Lh_{\alpha \beta }=\lambda h_{\alpha \beta }.$$
(14)
Here $`\mathrm{\Delta }_L`$ is the Lichnerowicz Laplacian on $`B`$. The stability of the spacetime was found to depend on the spectrum $`\{\lambda \}`$ of the Lichnerowicz Laplacian acting on transverse tracefree modes on $`B`$. Concretely, if the spectrum contains an eigenvalue that is too negative, the spacetime is unstable:

$$\lambda _{min}<\lambda _c4\frac{(5d)^2}{4}\text{instability}.$$
(15)
This result follows from considering the behaviour of the radial dependence of the perturbation (2.2), $`\varphi (r)`$. When the criterion (15) is satisfied, the solution for $`\varphi (r)`$ that decays at infinity also oscillates logarithmically in the interior. This allows it to be matched to a solution that is well behaved at the horizon $`r=\alpha `$. Thus a finite energy mode exists in this case and the spacetime is unstable.

It is perhaps worth pointing out here that a form of Birkhoff theorem holds for the metrics we are considering. In other words if we had assumed a general time dependent metric of the form

$$ds^2=e^{2\varphi (r,t)}dt^2+e^{2\psi (r,t)}dr^2+Y^2(r,t)ds_d^2,$$
(16)
where $`ds_d^2`$ is a $`d`$-dimensional time-independent Einstein metric with scalar curvature $`d(d1)`$, we would have found, on imposing the Einstein equations, that the metrics had to be static. In fact this result also holds if the metric is coupled to a 2-form field strength (in the electric case) or a $`d`$-form field strength the magnetic case. It also holds if one includes a cosmological term. It means that when perturbing the static metric we must consider time dependent perturbations of the transverse or base Einstein metric $`ds_d^2`$. This is another way of seeing why we need information about the the spectrum of the Lichnerowicz operator on this space.

We shall not give a detailed proof of Birkhoff’s theorem here, but merely indicate how to modify an existing treatment of Wiltshire  which assumes $`SO(d+1)`$-invariance, i.e that $`ds_d^2`$ is the unit round metric on $`S^d`$. Wiltshire gives a proof which also covers the case when Gauss-Bonnet terms are present. The argument he presents will not go over to the case of a general Einstein metric, since it makes special use of properties of its Riemann tensor. Thus in what follows we ignore that term, which means we set $`\stackrel{~}{\alpha }=0`$ in his equations. It is an interesting question to ask whether Birkhoff’s theorem remains true when one includes a Gauss-Bonnet term.

The discussion depends upon whether $`_\alpha Y`$ is spacelike, timelike or null. We assume the first case, and make a coordinate choice such that $`Y=r`$. The field equation $`R_t{}_{}{}^{r}=0`$ (equation (6b) in his paper) then yields

$$_t\psi =0.$$

The equation $`R_t{}_{}{}^{t}+R_r{}_{}{}^{r}=0`$ then gives

$$_r\varphi +_r\psi =0.$$

This means that $`\varphi +\psi =f(t)`$, where $`f(t)`$ is an arbitrary differentiable function of $`t`$. By choice of the coordinate $`t`$, $`f(t)`$ may be taken to vanish. It follows that both $`\varphi `$ and $`\psi `$ are independent of time $`t`$, and the metric is therefore static. The remaining field equations show that it takes the Schwarzschild-Tangherlini form.

### 2.3 Anti-de Sitter product spacetimes

These are solutions to a system with a $`d`$-form field strength minimally coupled to gravity

$`d\widehat{s}^2=ds_{AdS_{Dd}}^2+ds_d^2,`$
$`F_d=\left({\displaystyle \frac{2(D2)(d1)}{Dd1}}\right)^{1/2}\text{vol}_B`$ (17)
where we have taken the same normalisation for the curvature of $`B`$ as previously (12).

It was shown in  that the dangerous mode is a transverse tracefree mode on the manifold $`B`$ multiplying a scalar on the AdS spacetime

$`\widehat{h}_{0a}=\widehat{h}_{1a}=0,`$
$`\widehat{h}_{\alpha \beta }=h_{\alpha \beta }(x)\varphi (y),`$ (18)
where $`y`$ are the coordinates on the AdS. The mode $`h_{\alpha \beta }(x)`$ on $`B`$ is an eigentensor as in (14). From the AdS point of view, $`\varphi (y)`$ is seen as a scalar field with mass given by

$$m^2=\lambda 2(d1).$$
(19)
Instability of massive scalars on AdS spacetime is expressed in terms of the Breitenlohner-Freedman bound . In our units, this reads

$$m^2\left(\frac{Dd1}{d1}\right)^2<\frac{(Dd1)^2}{4}\text{instability}.$$
(20)
Using the value of the mass in (19), the criterion (20) is just

$$\lambda _{min}<\lambda _c4\frac{(5d)^2}{4}\text{instability}.$$
(21)
This is immediately seen to be the same as the criterion found for the black hole spacetimes (15). This is an intriguing match.

### 2.4 Lichnerowicz Laplacian and the Weyl tensor

In order to make estimates of the lowest eigenvalue of the Lichnerowicz Laplacian, it is convenient first to rewrite (7) in terms of the Weyl tensor. Since we are assuming that the metric is Einstein, with $`R_{ab}=\mathrm{\Lambda }g_{ab}`$, this is given in $`d`$ dimensions by

$$C_{abcd}=R_{abcd}\frac{\mathrm{\Lambda }}{d1}(g_{ac}g_{bd}g_{ad}g_{bc}).$$
(22)
Thus the Lichnerowicz Laplacian becomes

$$\mathrm{\Delta }_Lh_{ab}=\text{    }\text{ }h_{ab}2C_{acbd}h^{cd}+\frac{2d\mathrm{\Lambda }}{d1}h_{ab}.$$
(23)
A method for obtaining a lower bound on the smallest eigenvalue of $`\mathrm{\Delta }_L`$ was introduced in . One considers the integral of $`3(_{(a}h_{bc)})^2`$, which, after performing an integration by parts and using the transversality and tracelessness of $`h_{ab}`$, gives

$`{\displaystyle _M}h^{ab}\mathrm{\Delta }_Lh_{ab}`$ $`=`$ $`{\displaystyle _M}[3(_{(a}h_{bc)})^24h_{ab}R_{abcd}h^{cd}+4\mathrm{\Lambda }h_{ab}h^{ab}],`$ (24)
$`=`$ $`{\displaystyle _M}[3(_{(a}h_{bc)})^24h_{ab}C_{abcd}h^{cd}+{\displaystyle \frac{4d\mathrm{\Lambda }}{d1}}h_{ab}h^{ab}]`$
$``$ $`{\displaystyle _M}(4h_{ab}C_{abcd}h^{cd}+{\displaystyle \frac{4d\mathrm{\Lambda }}{d1}}h_{ab}h^{ab}).`$
Viewing $`C_{abcd}`$ as a map acting on symmetric traceless tensors $`h_{ab}`$, we can define its eigenvalues $`\kappa `$ by

$$C_{acbd}h^{cd}=\kappa h_{ab}.$$
(25)
For homogeneous spaces, we therefore have a simple inequality

$$\mathrm{\Delta }_L\frac{4d\mathrm{\Lambda }}{d1}4\kappa _{\mathrm{max}},$$
(26)
where $`\kappa _{\mathrm{max}}`$ is the largest eigenvalue of the Weyl tensor. Equality is attained if the corresponding eigentensor $`h_{ab}`$ is a Staeckel tensor, satisfying $`_{(a}h_{bc)}=0`$. The bound (26) was derived for homogeneous Einstein 7-manifolds $`M_7`$ in . It was also shown that in the case of $`U(1)`$ bundles over the Einstein-Kähler product 6-manifolds $`S^2\times ^2`$ and $`S^2\times S^2\times S^2`$, the eigentensor that maximises $`\lambda `$ is in fact Staeckel, and thus the equality in (26) is attained .

For inhomogeneous spaces, such as the Bohm metrics that we shall be studying later in this paper, the bound (24) must be kept in its integrated form, and so we have

$$_Mh^{ab}\mathrm{\Delta }_Lh_{ab}_M\left(\frac{4d\mathrm{\Lambda }}{d1}4\kappa _{\mathrm{max}}(x)\right)h_{ab}h^{ab},$$
(27)
where $`\kappa _{\mathrm{max}}(x)`$ represents the (position dependent) largest eigenvalue of the Weyl tensor. If we define $`\kappa _{\mathrm{max}}^0`$ to be the largest value attained by any of the eigenvalues of the Weyl tensor anywhere in $`M`$, $`\kappa _{\mathrm{max}}^0\text{sup}_x\kappa _{\mathrm{max}}(x)`$, then we obtain the inequality

$$\mathrm{\Delta }_L\frac{4d\mathrm{\Lambda }}{d1}4\kappa _{\mathrm{max}}^0.$$
(28)
It is clear, however, that this bound is not likely to be very sharp, especially if $`\kappa _{\mathrm{max}}`$ depends strongly on position.

We can, nevertheless, extract a general feature of the spectrum of the Lichnerowicz operator from the above considerations, namely that in Einstein spaces of positive Ricci tensor, the minimum eigenvalue tends to be lowered by having a large positive eigenvalue of the Weyl tensor. Indeed, we can see from (28) that if the largest Weyl tensor eigenvalue is not sufficiently positive, then $`\mathrm{\Delta }_L`$ could never be negative or zero. We shall see when we study the Bohm metrics in detail that in these cases the Weyl tensor can in fact have sufficiently large eigenvalues that the Lichnerowicz spectrum includes negative eigenvalues. By using a Rayleigh-Ritz variational method we shall be able to obtain upper bounds on the lowest eigenvalue of the Lichnerowicz Laplacian, allowing us in some cases to give an analytic proof of the existence of negative-eigenvalue modes.

## 3 The Bohm Einstein Metrics on $`S^N`$ and $`S^{N_1}\times S^{N_2}`$

### 3.1 Description of the Bohm construction

Bohm’s construction  gives rise to a countable infinity of Einstein metrics with positive Ricci tensor on the spheres $`S^N`$ for $`5N9`$, and on the product topologies $`S^{N_1}\times S^{N_2}`$ for $`5N_1+N_29`$ with $`N_12`$ and $`N_22`$. We can use these metrics in the black-hole and AdS$`\times M_d`$ spacetimes of the previous section.

The starting point for the construction is the following ansatz for metrics of cohomogeneity one,

$$ds^2=dt^2+a^2d\mathrm{\Omega }_p^2+b^2d\stackrel{~}{\mathrm{\Omega }}_q^2,$$
(29)
where $`a`$ and $`b`$ are functions of the radial variable $`t`$, and $`d\mathrm{\Omega }_p^2`$ and $`d\stackrel{~}{\mathrm{\Omega }}_q^2`$ are the standard metrics on the unit spheres $`S^p`$ and $`S^q`$. An elementary calculation shows that in the orthonormal frame $`e^0=dt`$, $`e^i=a\overline{e}^i`$, $`\widehat{e}^\alpha =b\overline{e}^\alpha `$, where $`\overline{e}^i\overline{e}^i=d\mathrm{\Omega }_p^2`$ and $`\overline{e}^\alpha \overline{e}^\alpha =d\stackrel{~}{\mathrm{\Omega }}_q^2`$, the components of the Riemann tensor are given by

$`R_{0i0j}={\displaystyle \frac{\ddot{a}}{a}}\delta _{ij},R_{0\alpha 0\beta }={\displaystyle \frac{\ddot{b}}{b}}\delta _{\alpha \beta },R_{i\alpha j\beta }={\displaystyle \frac{\dot{a}\dot{b}}{ab}}\delta _{ij}\delta _{\alpha \beta },`$
$`R_{ijk\mathrm{}}={\displaystyle \frac{1\dot{a}^2}{a^2}}(\delta _{ik}\delta _k\mathrm{}\delta _i\mathrm{}\delta _{jk}),R_{\alpha \beta \gamma \delta }={\displaystyle \frac{1\dot{b}^2}{b^2}}(\delta _{\alpha \gamma }\delta _{\beta \delta }\delta _{\alpha \delta }\delta _{\beta \gamma }),`$ (30)
The components of the Ricci tensor are given by

$`R_{00}`$ $`=`$ $`{\displaystyle \frac{p\ddot{a}}{a}}{\displaystyle \frac{q\ddot{b}}{b}},`$
$`R_{ij}`$ $`=`$ $`\left[{\displaystyle \frac{\ddot{a}}{a}}+{\displaystyle \frac{q\dot{a}\dot{b}}{ab}}+{\displaystyle \frac{(p1)(\dot{a}^21)}{a^2}}\right]\delta _{ij},`$
$`R_{\alpha \beta }`$ $`=`$ $`\left[{\displaystyle \frac{\ddot{b}}{b}}+{\displaystyle \frac{p\dot{a}\dot{b}}{ab}}+{\displaystyle \frac{(q1)(\dot{b}^21)}{b^2}}\right]\delta _{\alpha \beta }.`$ (31)
The Einstein equations $`R_{ab}=\mathrm{\Lambda }g_{ab}`$ give rise to two second-order differential equations for $`a`$ and $`b`$,

$`{\displaystyle \frac{\ddot{a}}{a}}+{\displaystyle \frac{q\dot{a}\dot{b}}{ab}}+{\displaystyle \frac{(p1)(\dot{a}^21)}{a^2}}`$ $`=`$ $`\mathrm{\Lambda },`$
$`{\displaystyle \frac{\ddot{b}}{b}}+{\displaystyle \frac{p\dot{a}\dot{b}}{ab}}+{\displaystyle \frac{(q1)(\dot{b}^21)}{b^2}}`$ $`=`$ $`\mathrm{\Lambda },`$ (32)
together with the first-order constraint

$$\frac{p(p1)(\dot{a}^21)}{a^2}+\frac{q(q1)(\dot{b}^21)}{b^2}+\frac{2pq\dot{a}\dot{b}}{ab}+(p+q1)\mathrm{\Lambda }=0.$$
(33)
We shall adopt the conventional normalisation, when considering Einstein metrics with positive $`\mathrm{\Lambda }`$, of taking

$$\mathrm{\Lambda }=p+q,$$
(34)
which is one less than the total dimension of the space.

When $`p>1`$ and $`q>1`$, which we shall be considering here, the general solution of the Einstein equations is not known explicitly. A well-known special solution is

$$a=\mathrm{sin}t,b=\mathrm{cos}t,$$
(35)
in which case the metric (29) becomes just the standard round metric on $`S^{p+q+1}`$,

$$ds^2=dt^2+\mathrm{sin}^2td\mathrm{\Omega }_p^2+\mathrm{cos}^2td\stackrel{~}{\mathrm{\Omega }}_q^2,$$
(36)
written as a foliation by $`S^p\times S^q`$. This can easily be recognised as the metric on the unit $`S^{p+q+1}`$ by introducing coordinates $`x^A`$ on $`^{p+q+2}`$, subject to the unit-radius constraint $`x^Ax^A=1`$, and then introducing orthogonal unit-vectors $`m^A`$ and $`n^A`$ in $`^{p+q+2}`$, such that a general point on the unit $`S^{p+q+1}`$ in $`^{p+q+2}`$ can be written as

$$x^A=m^A\mathrm{sin}t+n^A\mathrm{cos}t.$$
(37)
A second well-known special solution to the Einstein equations is

$$a=\sqrt{\frac{p}{\mathrm{\Lambda }}}\mathrm{sin}\left(\sqrt{\frac{\mathrm{\Lambda }}{p}}t\right),b=\sqrt{\frac{q1}{\mathrm{\Lambda }}},$$
(38)
with, using our conventional choice, $`\mathrm{\Lambda }=p+q`$. This gives the standard homogeneous Einstein metric on $`S^{p+1}\times S^q`$,

$$ds^2=dt^2+\frac{p}{\mathrm{\Lambda }}\mathrm{sin}^2\left(\sqrt{\frac{\mathrm{\Lambda }}{p}}t\right)d\mathrm{\Omega }_p^2+\frac{q1}{\mathrm{\Lambda }}d\stackrel{~}{\mathrm{\Omega }}_q^2.$$
(39)
There is an analogous solution for $`S^p\times S^{q+1}`$ too. Since there is obviously always a discrete transformation under which the roles of the spheres $`S^p`$ and $`S^q`$ are interchanged, we shall not in general bother to mention the symmetry-related possibility.

It is shown in  that the Einstein equations (32) and (33) admit a countably infinite number of solutions giving rise to inequivalent metrics that extend smoothly onto manifold of topology $`S^{p+q+1}`$, and another countable infinity of solutions for which the metrics extend smoothly onto manifolds of topology $`S^{p+1}\times S^q`$. We shall denote these metrics by Bohm$`(p,q)_n`$, where the integer $`n`$ runs over $`n=0,2,4,\mathrm{}`$ for the $`S^{p+q+1}`$ sequence, and $`n=1,3,5\mathrm{}`$ for the $`S^{p+1}\times S^q`$ sequence. The standard unit metric (35) on $`S^{p+q+1}`$ corresponds to Bohm$`(p,q)_0`$, and the standard product Einstein metric (38) on $`S^{p+1}\times S^q`$ corresponds to Bohm$`(p,q)_1`$.

The higher metrics Bohm$`(p,q)_n`$ with $`n2`$ are all inhomogeneous. The radial coordinate $`t`$ runs between endpoints which can be taken to be 0 and $`t_f`$, defined by the vanishing of one or other of the metric functions $`a`$ and $`b`$. The metric extends onto the corresponding degenerate orbit because the associated metric function vanishes like $`t`$ or $`(t_ft)`$, so that one has a regular collapsing of $`p`$-spheres or $`q`$-spheres like in the origin of spherical polar coordinates. For all the $`S^{p+q+1}`$ metrics Bohm$`(p,q)_{2m}`$ one has

$`a(0)=0,\dot{a}(0)=1,b(0)=b_0,\dot{b}(0)=0;`$
$`a(t_f)=a_0,\dot{a}(t_f)=0,b(t_f)=0,\dot{b}(t_f)=1.`$ (40)
On the other hand, for the $`S^{p+1}\times S^q`$ metrics Bohm$`(p,q)_{2m+1}`$ one has

$`a(0)=0,\dot{a}(0)=1,b(0)=b_0,\dot{b}(0)=0;`$
$`a(t_f)=0,\dot{a}(t_f)=1,b(t_f)=\stackrel{~}{b}_0,\dot{b}(t_f)=0.`$ (41)
The functions $`a`$ and $`b`$ are strictly positive for $`0<t<t_f`$, and the quantities $`a_0`$, $`b_0`$ and $`\stackrel{~}{b}_0`$ are certain constants.

Plots of the metric functions $`a`$ and $`b`$ for various Bohm metrics are presented in the Appendix. These have been obtained by performing a numerical integration of the Einstein equations (32). It can be seen that as the index $`n`$ for Bohm$`(p,q)_n`$ increases, the metrics rapidly become approximations to the “double-cone” Einstein metric

$$ds^2=dt^2+\frac{1}{(p+q1)}\mathrm{sin}^2t[(p1)d\mathrm{\Omega }_p^2+(q1)d\stackrel{~}{\mathrm{\Omega }}_q^2],$$
(42)
for most of the range of the radial coordinate. The metric (42) itself is singular at the apexes $`t=0`$ and $`t=\pi `$, since near to each of these points one has a collapse of $`S^p\times S^q`$ surfaces. The actual Bohm$`(p,q)_n`$ metrics with large $`n`$ deviate from (42) just in the vicinity of the apexes, instead approaching the forms given in (40) or (41). It is interesting to note that (42) is in fact the singular limit both of the regular $`S^{p+q+1}`$ sequence Bohm$`(p,q)_{2m}`$ and the regular $`S^{p+1}\times S^q`$ sequence Bohm$`(p,q)_{2m+1}`$.

Note that in the case of the Bohm metrics Bohm$`(p,q)_{2m+1}`$ with the topology $`S^{p+1}\times S^q`$, the fact that the metric function $`b(t)`$ never vanishes means that we can replace the associated round sphere $`S^q`$ with its metric $`d\stackrel{~}{\mathrm{\Omega }}_q^2`$ in (29) by any Einstein space $`Q_q`$ of dimension $`q`$, whose (positive) Ricci curvature is normalised to $`\stackrel{~}{R}_{\alpha \beta }=(q1)\stackrel{~}{g}_{\alpha \beta }`$, and we will again have a complete and non-singular Bohmian metric in $`(p+q+1)`$ dimensions, now with the topology $`S^{p+1}\times Q_q`$. The Einstein space $`Q_q`$ could itself be taken to be a Bohm metric such as Bohm$`(2,2)_n`$, Bohm$`(2,3)_n`$ or Bohm$`(3,2)_n`$.

### 3.2 Estimates and bounds for Lichnerowicz in Bohm metrics

#### 3.2.1 Eigenvalues of the Weyl tensor

We saw earlier, in section 2.4, that positive eigenvalues of the Weyl tensor tend to drive the lowest mode of the Lichnerowicz more negative. Accordingly, we can gain insights into the bounds on the spectrum of the Lichnerowicz operator in the Bohm metrics by studying the Weyl tensor. From (30), and our choice of normalisation where $`R_{ab}=(p+q)g_{ab}`$, we have

$`C_{0i0j}=x_1\delta _{ij},C_{0\alpha 0\beta }=x_2\delta _{\alpha \beta },C_{i\alpha j\beta }=x_3\delta _{ij}\delta _{\alpha \beta },`$
$`C_{ijk\mathrm{}}=x_4(\delta _{ik}\delta _j\mathrm{}\delta _i\mathrm{}\delta _{jk}),C_{\alpha \beta \gamma \delta }=x_5(\delta _{\alpha \gamma }\delta _{\beta \delta }\delta _{\alpha \delta }\delta _{\beta \gamma }),`$ (43)
where

$$x_1=1\frac{\ddot{a}}{a},x_2=1\frac{\ddot{b}}{b},x_3=1\frac{\dot{a}\dot{b}}{ab},x_4=\frac{1\dot{a}^2a^2}{a^2},x_5=\frac{1\dot{b}^2b^2}{b^2}.$$
(44)
Note that these coefficients $`x_i`$ are not all independent, and thus

$$px_1=qx_2=\frac{1}{2}q(q1)x_5\frac{1}{2}p(p1)x_4,x_3=\frac{p(p1)x_4+q(q1)x_5}{2pq}.$$
(45)
It is straightforward to see that with multiplicities $`𝐦`$, the traceless eigenvectors $`h_{AB}`$ and eigenvalues $`\kappa `$ of the Weyl tensor are given by

$`h_{0i}:`$ $`𝐦=p,\kappa =x_1,`$
$`h_{0\alpha }:`$ $`𝐦=q,\kappa =x_2,`$
$`\{h_{ij}|h_{ii}=0\}:`$ $`𝐦=\frac{1}{2}p(p+1)1,\kappa =x_4,`$
$`\{h_{\alpha \beta }|h_{\alpha \alpha }=0\}:`$ $`𝐦=\frac{1}{2}q(q+1)1,\kappa =x_5,`$
$`h_{i\alpha }:`$ $`𝐦=pq,\kappa =x_3,`$ (46)
together with two eigenvectors of the form

$$h_{00}=puqv,h_{ij}=u\delta _{ij},h_{\alpha \beta }=v\delta _{\alpha \beta },$$
(47)
for which the eigenvalues are given by the roots of a quadratic equation,

$`\kappa _\pm `$ $`=`$ $`\frac{1}{2}(p1)x_4+\frac{1}{2}(q1)x_5`$ (48)
$`\pm {\displaystyle \frac{\sqrt{pq[(p1)x_4+(q1)x_5]^2+(p+q+1)[p(p1)x_4q(q1)x_5]^2}}{2\sqrt{pq}}}.`$
The coefficients $`u`$ and $`v`$ are then given by

$$u=2q\kappa _\pm p(p1)x_4+q(q1)x_5,v=2p\kappa _\pm p(p1)x_4+q(q1)x_5.$$
(49)
In total, we have the expected $`\frac{1}{2}(p+q+1)(p+q+2)1`$ symmetric traceless eigenmodes in $`(p+q+1)`$ dimensions.

Using the output of the numerical integration of the Einstein equations for the Bohm metrics, we find that the eigenvalue of the Weyl-tensor that achieves the largest positive value is $`\kappa _+`$ given by (48). It is therefore in this sector that one can expect to find the lowest-lying eigenmodes of the Lichnerowicz operator. We can see from (47) that the associated eigenvector is of a type that may be thought of as a “ballooning mode.” That is to say, it is a mode where one of the spheres $`S^p`$ or $`S^q`$ tends to inflate at the expense of the other. This accords with one’s intuition, which would suggest that the most likely instability for metrics with direct-product orbits would be ballooning modes of this general type.

Some examples of our numerical results for the maximum value of the largest eigenvalue of the Weyl tensor are as follows. For the Bohm$`(2,2)_n`$ metrics on $`S^5`$ and $`S^3\times S^2`$ we find from (48) that $`\kappa _+`$ attains its maximum value at the endpoints of the radial coordinate range, and so

$$\kappa _{\mathrm{max}}^0=\frac{5(1b_0^2)}{3b_0^2}.$$
(50)
In fact, as $`n`$ increases, the function $`\kappa _+`$ peaks more and more strongly around the endpoints. For $`n=0,\mathrm{},6`$, we have approximately

$$\kappa _{\mathrm{max}}^0=\{0,5,24.26,118.45,579.76,3013.72,15106.9\}.$$
(51)
(The results for $`n=0`$ and $`n=1`$ are exact, since these the standard homogeneous metrics on $`S^5`$ and $`S^3\times S^2`$.) Using (28), we obtain the lower bound

$$\mathrm{\Delta }_L\frac{20(14b_0^2)}{3b_0^2}$$
(52)
on the spectrum of the Lichnerowicz operator on Bohm$`(2,2)_n`$. From our numerical results for the values of $`b_0`$ for the first few examples, we find

$`\text{Bohm}(2,2)_0:`$ $`\mathrm{\Delta }_L20,`$
$`\text{Bohm}(2,2)_1:`$ $`\mathrm{\Delta }_L0,`$
$`\text{Bohm}(2,2)_2:`$ $`\mathrm{\Delta }_L77.04,`$
$`\text{Bohm}(2,2)_3:`$ $`\mathrm{\Delta }_L453.8,`$
$`\text{Bohm}(2,2)_4:`$ $`\mathrm{\Delta }_L2342,`$
$`\text{Bohm}(2,2)_5:`$ $`\mathrm{\Delta }_L11972,`$
$`\text{Bohm}(2,2)_6:`$ $`\mathrm{\Delta }_L60407.`$ (53)
The bounds for $`n=0`$ and $`n=1`$ are in fact exactly attained, corresponding to the the cases of the homogeneous $`S^5`$ and $`S^3\times S^2`$ metrics respectively. The zero-mode in the latter case is the ballooning mode on $`S^3\times S^2`$. As we mentioned previously, the lower bounds we obtain for the inhomogeneous Bohm metrics are not expected to be very sharp.

In general for the Bohm$`(p,q)_n`$ metrics we find

$$\kappa _{\mathrm{max}}^0=\frac{(q1)(p+q+1)(1b_0^2)}{(p+1)b_0^2},$$
(54)
and hence we have the lower bound

$$\mathrm{\Delta }_L\frac{4(p+q+1)[q1(p+q)b_0^2]}{(p+1)b_0^2},$$
(55)
where as usual we have normalised the scale so that $`R_{ab}=(p+q)g_{ab}`$.

#### 3.2.2 Transverse tracefree ballooning modes

In order to study transverse tracefree perturbations, we consider a metric of the form

$$ds^2=c^2dt^2+a^2d\mathrm{\Omega }_p^2+b^2d\stackrel{~}{\mathrm{\Omega }}_q^2.$$
(56)
This is similar to (29), except that we have, for convenience, introduce the coordinate gauge function $`c(t)`$ in the metric. Substituting into the Einstein-Hilbert action

$$S=\sqrt{g}d^dx\left[R(d2)(d1)\right],$$
(57)
(where $`d=p+q+1`$) and omitting a constant factor equal to the volume of the product metric on the unit $`S^p\times S^q`$, this gives

$`S`$ $`=`$ $`{\displaystyle }a^pb^qc[2pq{\displaystyle \frac{\dot{a}\dot{b}}{abc^2}}+p(p1){\displaystyle \frac{\dot{a}^2}{a^2c^2}}+q(q1){\displaystyle \frac{\dot{b}^2}{b^2c^2}}+{\displaystyle \frac{p(p1)}{a^2}}+{\displaystyle \frac{q(q1)}{b^2}}`$ (58)
$`(p+q)(p+q1)]dt.`$
Ballooning modes in product metrics, in which one factor contracts and the other expands, are Lichnerowicz zero modes and are typically associated with instabilities . It is reasonable to expect, therefore, as we argued in section 3.2.1, that if instabilities were to arise in the Bohm metrics, they would be associated with modes of a similar type. We are therefore led to seek a generalisation of ballooning modes to the warped product of spheres present in (56). The perturbation

$$aa\sqrt{1+u},bb\sqrt{1+v},c=1\sqrt{1+\gamma },$$
(59)
is tracefree at the linearised level if $`\gamma +pu+qv=0`$, and transverse if

$$\dot{\gamma }+(p\frac{\dot{a}}{a}+q\frac{\dot{b}}{b})\gamma p\frac{\dot{a}}{a}uq\frac{\dot{b}}{b}v=0.$$
(60)
These two conditions can be used in order to solve for $`u`$ and $`v`$ in terms of $`\gamma `$:

$$u=\frac{\dot{\gamma }+\left[(q+1)\dot{b}b^1+p\dot{a}a^1\right]\gamma }{p(\dot{a}a^1\dot{b}b^1)},v=\frac{\dot{\gamma }+\left[(p+1)\dot{a}a^1+q\dot{b}b^1\right]\gamma }{q(\dot{b}b^1\dot{a}a^1)}.$$
(61)
One is free to choose the function $`\gamma `$, which completely determines the perturbation through (61). However, the Bohm$`(p,q)_n`$ metric has $`n`$ interior points at which $`(\dot{a}a^1\dot{b}b^1)`$ vanishes, and hence the expressions for $`u`$ and $`v`$ are singular for generic choices of $`\gamma `$. This problem can be solved by inverting (61) to give $`\gamma `$ in terms of $`u`$,

$$\gamma =\frac{1}{a^pb^{q+1}}a^pb^{q+1}\left[\frac{\dot{a}}{a}\frac{\dot{b}}{b}\right]u𝑑t.$$
(62)
The remaining function, $`v`$, is given by the tracefree condition, $`pu+qv+\gamma =0`$. One is now free to choose a nonsingular function $`u`$ to obtain a perturbation that will at worst be singular at the endpoints $`t=0`$, $`t=t_f`$. These singularities can be avoided as described in the next subsection. One simple choice is

$$u=\frac{pm}{pa^mb^{p+q+1m}},v=\frac{(p+1m)}{qa^mb^{p+q+1m}},\gamma =\frac{1}{a^mb^{p+q+1m}}.$$
(63)
We will use these expressions below in order to exhibit negative Lichnerowicz modes in some of the Bohm metrics.

#### 3.2.3 Rayleigh-Ritz estimates for the lowest Lichnerowicz eigenvalue

In the eigenvalue problem $`\mathrm{\Delta }\varphi =\lambda \varphi `$ for a self-adjoint operator $`\mathrm{\Delta }`$, one can obtain an upper bound on the lowest eigenvalue by the Rayleigh-Ritz method, namely

$$\lambda _{\mathrm{min}}\frac{_M\psi \mathrm{\Delta }\psi }{_M\psi ^2},$$
(64)
with equality being achieved if the trial function $`\psi `$ is actually the eigenfunction corresponding to the lowest eigenvalue. In this section, we shall apply this method to obtain an upper bound on the lowest eigenvalue of the Lichnerowicz operator on TT modes in the Bohm metrics, and in particular, we shall find that there is a negative-eigenvalue mode in some of cases we examine.

The easiest cases to consider are the Bohm metrics Bohm$`(p,q)_{2m+1}`$ on the product topologies $`S^{p+1}\times S^q`$. In these metrics, the function $`b`$ in (29) is nowhere vanishing, and so we can take our trial function to be given by (63) with $`m=0`$:

$$u=\frac{ϵ}{b^{p+q+1}},v=\frac{(p+1)ϵ}{qb^{p+q+1}},\gamma =\frac{ϵ}{b^{p+q+1}}.$$
(65)
Here we have introduced $`ϵ`$ as a small constant order-parameter. As we shall see below, this trial function allows us to prove that certain of the Bohm metrics on products of spheres have negative eigenvalue modes of the Lichnerowicz operator.

### 3.3 Negative Lichnerowicz eigenvalues in Bohm metrics

#### 3.3.1 Analytic results for negative modes for Bohm metrics on $`S^3\times S^2`$

It turns out that the easiest cases to study are the Bohm metrics Bohm$`(2,2)_{2m+1}`$, whose topology is $`S^3\times S^2`$. We are able to obtain completely analytic and explicit results that prove the existence of negative modes of the Lichnerowicz operator for all these examples (for $`m1`$), and so we shall present the details for these metrics here. For these examples, we take $`\gamma =ϵ/b^5`$ as our trial function, as suggested by (65), implying that we have $`u=ϵ/b^5`$ and $`v=3ϵ/(2b^5)`$. From the expansion of $`S`$ given in (58), with the perturbation (59), we can easily extract the terms quadratic in $`ϵ`$, and so by comparing with (6) and (9) we obtain a Rayleigh-Ritz bound for the lowest eigenvalue $`\lambda _{\mathrm{min}}`$ of the Lichnerowicz operator:

$$\lambda _{\mathrm{min}}\frac{_0^{t_f}P𝑑t}{_0^{t_f}Q𝑑t},$$
(66)
where

$`P`$ $`=`$ $`{\displaystyle \frac{3a^2}{b^{10}}}+{\displaystyle \frac{7}{b^8}}{\displaystyle \frac{30a^2}{b^8}}+{\displaystyle \frac{7\dot{a}^2}{b^8}}{\displaystyle \frac{112a\dot{a}\dot{b}}{b^9}}{\displaystyle \frac{91a^2\dot{b}^2}{2b^8}},`$
$`Q`$ $`=`$ $`{\displaystyle \frac{15a^2}{2b^8}},`$ (67)
and $`t_f=2t_c`$ is the upper limit of the range of the radial coordinate $`t`$. Note that in fact it suffices to evaluate the integrals only up to the mid-point $`t=t_c`$ in these cases, since the metric functions are symmetric about $`t=t_c`$ here.

Using the constraint (33), with $`p=q=2`$, we can eliminate the term involving $`\dot{a}^2`$ in $`P`$. We also note that upon use of the second-order equations (32), we can prove the identities

$`{\displaystyle \frac{d}{dt}}\left({\displaystyle \frac{a\dot{a}}{b^8}}\right)`$ $`=`$ $`{\displaystyle \frac{10a\dot{a}\dot{b}}{b^9}}+{\displaystyle \frac{14a^2}{b^8}},`$
$`{\displaystyle \frac{d}{dt}}\left({\displaystyle \frac{a^2\dot{b}}{b^9}}\right)`$ $`=`$ $`{\displaystyle \frac{10a^2\dot{b}^2}{b^{10}}}+{\displaystyle \frac{a^2(14b^2)}{b^{10}}}.`$ (68)
Since $`a\dot{a}/b^8`$ and $`a^2\dot{b}/b^9`$ vanish at both endpoints of the full integration range, we can use these in order to perform integrations by parts in the evaluation of $`P`$. Specifically, we use the former to remove the term in $`P`$ involving $`\dot{a}\dot{b}`$, and then using the latter, we find that

$$_0^{t_f}P𝑑t=\frac{25}{2}_0^{t_f}\frac{a^2\dot{b}^2}{b^{10}}𝑑t,$$
(69)
which is manifestly non-positive. We therefore have the Rayleigh-Ritz bound

$$\lambda _{\mathrm{min}}\frac{5}{3}\frac{_0^{t_f}a^2b^{10}\dot{b}^2𝑑t}{_0^{t_f}a^2b^8𝑑t}$$
(70)
for the lowest eigenvalue of the Lichnerowicz operator. (Recall that we are working in units where $`R_{ab}=4g_{ab}`$.) This proves that the Einstein metrics Bohm$`(2,2)_{2m+1}`$ have a negative eigenvalue for the Lichnerowicz operator on transverse traceless symmetric 2-index tensors, for $`m1`$. (The case $`m=0`$ is the standard product Einstein metric on $`S^3\times S^2`$, with $`b=\frac{1}{2}`$. In this case the numerator gives zero, and in fact we exactly saturate the upper bound, finding the known lowest eigenvalue $`\mathrm{\Delta }_L=0`$ on the product metric.)

Using the numerical results described in the Appendix, we find the following upper bounds on the lowest Lichnerowicz eigenvalue for the Bohm$`(2,2)_3`$ and Bohm$`(2,2)_5`$ Einstein metrics on $`S^3\times S^2`$:

$`\text{Bohm}(2,2)_3:`$ $`\lambda _{\mathrm{min}}7.2766,`$
$`\text{Bohm}(2,2)_5:`$ $`\lambda _{\mathrm{min}}198.008.`$ (71)
(Recall that we are normalising the metrics so that $`R_{ab}=4g_{ab}`$.) As one goes to higher examples Bohm$`(2,2)_{2m+1}`$ with increasing $`m`$, one finds that the upper bound on the lowest Lichnerowicz eigenvalue becomes increasingly negative, tending to $`\mathrm{}`$ in the limit as $`m\mathrm{}`$. Note that our upper bounds (71) are considerably larger than the rather crude lower bounds (53) that we obtained by considering the eigenvalues of the Weyl tensor.

#### 3.3.2 Analytic results for negative modes for Bohm metrics on $`S^3\times S^3`$

For general values of $`p`$ and $`q`$, the analogous expression for the integrands $`Y`$ in (69) and $`Q`$ in (67) that appear in the numerator and denominator of the Rayleigh-Ritz functional (66) turn out to be

$`P`$ $`=`$ $`q^1(p+q+1)^2(p^22p3+pqq)a^pb^{2pq4}\dot{b}^2,`$
$`Q`$ $`=`$ $`q^1(p+q+1)(p+1)a^pb^{2pq2},`$ (72)
if we take the trial function $`\gamma =b^{pq1}`$. In general, this gives us a rather weak positive upper bound on the lowest Lichnerowicz eigenvalue for the Bohm$`(p,q)_{2m+1}`$ metrics on $`S^{p+1}\times S^q`$. In fact only for $`p=q=2`$, which we discussed above, and $`p=2`$, $`q=3`$, does one get a non-positive bound from this choice of trial function. Interestingly, for $`p=2`$, $`q=3`$ the numerator integrand $`P`$ in (72) vanishes identically, and so we obtain the bound

$$\lambda _{\mathrm{min}}0$$
(73)
in this case. It is straightforward to see that for $`m1`$ the trial function $`\gamma =1/b^6`$ does not give an eigenfunction, and hence the inequality in (73) is not saturated. Thus we have an analytic proof that for the Bohm$`(2,3)_{2m+1}`$ Einstein metrics on $`S^3\times S^3`$ with $`m1`$, the lowest eigenvalue of the Lichnerowicz operator on TT symmetric tensors is strictly negative.

For all other cases aside from Bohm$`(2,2)_{2m+1}`$ and Bohm$`(2,3)_{2m+1}`$, the trial function $`\gamma =b^{pq1}`$ does not give a negative upper bound on $`\lambda _{\mathrm{min}}`$. We believe that this is a consequence of a non-optimal choice of trial function, since the qualitative arguments would suggest the existence of negative Lichnerowicz modes for all the Bohm metrics.

#### 3.3.3 Numerical results for negative modes for Bohm$`(2,2)_{2m}`$ metrics

The analytic methods that allowed us to prove the existence of negative modes of the Lichnerowicz operator in the $`S^3\times S^2`$ and $`S^3\times S^3`$ Bohm metrics do not directly extend to any of the Bohm metrics on $`S^{p+q+1}`$. The reason for this is that whilst our trial function $`\gamma =1/b^{p+q+1}`$ is regular everywhere in the metrics on $`S^{p+1}\times S^q`$, it diverges at the right-hand endpoint of the range of the radial coordinate $`t`$ in the metrics on $`S^{p+q+1}`$. A natural modification to the trial function to take account of this is to interpolate smoothly between $`\gamma =1/b^{p+q+1}`$ on the left-hand side and $`\gamma =1/a^{p+q+1}`$ on the right-hand side of the range of $`t`$. We have carried out this procedure numerically in the case of examples of the Bohm$`(2,2)_{2m}`$ metrics on $`S^5`$, and we find that indeed there are negative modes of the Lichnerowicz operator, in accordance with the qualitative arguments. Specifically, we find approximately

$`\text{Bohm}(2,2)_2:`$ $`\lambda _{\mathrm{min}}0.7937,`$
$`\text{Bohm}(2,2)_4:`$ $`\lambda _{\mathrm{min}}38.86,`$
$`\text{Bohm}(2,2)_6:`$ $`\lambda _{\mathrm{min}}1040.6.`$ (74)
These upper bounds are again all considerably larger than the corresponding lower bounds in (53) that we obtained from the eigenvalues of the Weyl tensor.

### 3.4 Non-compact Bohm metrics

A class of complete and non-singular non-compact metrics was also constructed by Bohm . These include examples where the metric ansatz is again taken to be (29), but now the metric is required to be Ricci-flat. These metrics have been considered recently in  in studies of the possibility of topology change. It was shown in  that regular metrics exist in which $`a(t)`$ and $`b(t)`$ satisfy the boundary conditions

$$a(0)=0,\dot{a}(0)=1,b(0)=b_0,\dot{b}(0)=0,$$
(75)
Unlike the previous compact examples, here regularity imposes no constraint on the allowed values for the constant $`b_0`$, and in fact the value of $`b_0`$ now merely sets the overall scale of the metric. Note that $`b`$ is everywhere non-vanishing, and so the there is an $`S^q`$ bolt at $`t=0`$. The metrics are asymptotically conical, approaching cones over the standard product Einstein metric on $`S^p\times S^q`$.

A representative example is presented in the Appendix, for the case of $`p=2`$, $`q=2`$.

The Rayleigh-Ritz method that we described earlier for finding an upper bound on the smallest Lichnerowicz eigenvalue can be applied in these non-compact Bohm metrics too. In fact the trial function $`\gamma =b^{pq1}`$ can be considered here too, since it remains finite everywhere and it falls off rapidly at large $`t`$. We find that the numerator and denominator integrands are then again given by (72), and so again we obtain a negative upper bound on the lowest Lichnerowicz eigenvalue for the cases $`p=q=2`$, and $`p=2`$, $`q=3`$.

Evaluating the integrands numerically for the case $`p=q=2`$, we find that

$$\lambda _{\mathrm{min}}\frac{0.110433}{b_0^2}.$$
(76)
For $`p=2`$, $`q=3`$, we find, as for Bohm$`(2,3)_{2m+1}`$, that the bound is $`\lambda _{\mathrm{min}}0`$, and we can again argue that since the trial function $`\gamma =1/b^6`$ does not give an eigenfunction, we must have $`\lambda _{\mathrm{min}}<0`$. For all other $`p`$ and $`q`$, our choice of trial function does not give a negative upper bound on the lowest eigenvalue of the Lichnerowicz operator. Again, we believe that this is because the trial function is non-optimal in these other cases, since general arguments suggest that the non-compact Bohm metrics should all have negative eigenvalue modes of the Lichnerowicz operator.

## 4 Einstein-Sasaki manifolds

### 4.1 Introduction and definition

In this section we remind the reader that as well as the infinite sequence of cohomogeneity one Bohm metrics that have featured in our discussion, the manifold $`S^3\times S^2`$ admits many other Einstein metrics. For example, it has been known for some time that there are infinitely-many homogeneous but non-supersymmetric $`T^{p,q}`$ spaces, corresponding to $`U(1)`$ bundles over $`S^2\times S^2`$ in which the $`U(1)`$ fibres wind $`p`$ times over one $`S^2`$, and $`q`$ times over the other. These all have the topology $`S^3\times S^2`$ and they all admit an Einstein metric. Only $`T^{1,1}`$ admits Killing spinors.

There are, by contrast, also many inequivalent supersymmetric examples of $`S^3\times S^2`$ Einstein metrics, which do admit Killing spinors. They can thus be used in the AdS/CFT correspondence, replacing $`S^5`$ in the D3 brane metric and its near-horizon limit.

An Einstein-Sasaki metric may be defined as a $`(2m+1)`$-dimensional Einstein metric such that the cone over it is a Calabi-Yau metric

$$ds_{\mathrm{Calabi}\mathrm{Yau}}^2=dR^2+R^2ds_{\mathrm{Einstein}\mathrm{Sasaki}}^2,$$
(77)
or in other words, the cone is a Ricci-flat Kähler metric. The Killing spinors in the Einstein-Sasaki metric come by direct projection from the covariantly-constant spinors of the Calabi-Yau metric. If one uses the complex structure $`J`$ of the Calabi-Yau to act on the Euler vector of the cone, $`R\frac{}{R}`$, one gets a Killing vector on the Einstein-Sasaki manifold with constant magnitude, and thus we may write locally

$$ds_{\mathrm{Einstein}\mathrm{Sasaki}}^2=(d\psi +A)^2+ds_{\mathrm{Einstein}\mathrm{K}\ddot{\mathrm{a}}\mathrm{hler}}^2,$$
(78)
where $`J\frac{}{R}=\frac{}{\psi }`$ and $`ds_{\mathrm{Einstein}\mathrm{K}\ddot{\mathrm{a}}\mathrm{hler}}^2`$ is locally Einstein-Kähler with positive scalar curvature. Globally the $`U(1)`$ action generated by $`\frac{}{\psi }`$ may be free (in which case one speaks of a regular Sasaki structure) and the base is a smooth Einstein-Kähler manifold, or it may have fixed points in which case the Einstein-Kähler base has orbifold singularities. The total space however will still be smooth. We give a more detailed discussion of the relation between the Einstein-Kähler and Einstein-Sasaki spaces below.

Taking $`^2`$ or $`^1\times ^1`$ as the Einstein-Kähler base metric gives the standard homogeneous Sasaski metrics on $`S^5`$ or $`T^{1,1}`$ respectively. If the fibration is regular the only remaining possible base metrics for 5-dimensional Einstein-Sasaki metrics are inhomogeneous metrics on del Pezzo surfaces, i.e. $`^2`$ blown up at $`k`$ points, with $`3k8`$, giving Einstein-Sasaki metrics on the connected sum of $`k`$ copies of $`S^3\times S^2`$. <sup>1</sup><sup>1</sup>1It is worth remarking that any 5-dimensional closed simply-connected spin manifold with no torsion in the second homology group is diffeomorphic to a connected sum of copies of $`S^3\times S^2`$.

Recently Boyer, Galicki et al.  have constructed many inhomogeneous Einstein-Sasaki $`(2n1)`$ metrics on the links $`L_f=C_fS^{2n+1}`$ of weighted homogeneous polynomials $`f`$ on $`^{n+1}`$. The notation is as follows: $`C_f^{n+1}`$ is the zero set $`f(z)=0`$ of the polynomial, and $`S^{2n+1}`$ is the standard sphere. One readily sees that the Hopf fibration descends to $`L_f`$, and this gives the fibration associated to the Sasaki structure. Note that this description is purely topological. The metric is obtained indirectly by means of an existence proof. The present state of the art is that there are at least 14 inequivalent Einstein-Sasaki structures on $`S^3\times S^2`$ . Of these, only $`T^{1,1}`$ is homogeneous, and so if used in the AdS/CFT correspondence the other 13 examples would give supersymmetric vacua with no R-symmetry.

The volume of $`T^{1,1}`$ is well known to be $`\frac{16\pi ^3}{27}`$. According to , the volume of $`L_f`$ is given, in five dimensions, by

$$\frac{\pi ^3}{27w}(|w|d)^3,$$
(79)
where $`d`$ is the degree of $`f`$, $`𝐰=(w_0,w_1,w_2,w_3)`$ are the weights, $`|𝐰|=w_0+w_1+w_2+w_3`$ and $`w=w_0w_1w_2w_3`$. The two inhomogeneous Einstein-Sasaki metrics on $`S^3\times S^2`$ constructed in  are both of degree 256 and have weights $`𝐰=(11,49,69,128)`$ and $`𝐰=(13,35,81,128)`$. They therefore have volumes $`\frac{\pi ^3}{27\times 4760448}`$ and $`\frac{\pi ^3}{27\times 4717440}`$ respectively. These may be compared with the volume of the product metric on $`S^3\times S^2`$, which is $`\frac{\pi ^3}{\sqrt{2}}`$, and of the limiting singular double cone Bohm metric, which is $`\frac{2\pi ^3}{3}`$.

### 4.2 Einstein-Sasaki manifolds as $`U(1)`$ bundles over Einstein-Kähler manifolds

In this subsection, we present what is essentially a review of how Einstein-Sasaki manifolds can be constructed as $`U(1)`$ bundles over Einstein-Kähler manifolds, focusing in particular on the construction of the Killing spinors. The construction can be applied to obtain Einstein-Sasaki manifolds in any odd dimension, and so we shall give the construction for this general case.

Suppose we have an Einstein-Kähler metric $`g_{ab}`$ on a manifold $`M_n`$ of (even) dimension $`n=2m`$. By the standard formulae of Kaluza-Klein reduction, the $`(n+1)`$-dimensional metric

$$d\widehat{s}^2=(d\psi +A)^2+ds^2$$
(80)
has Ricci tensor $`\widehat{R}_{AB}`$ whose frame components are given by

$$\widehat{R}_{ab}=R_{ab}\frac{1}{2}F_a{}_{}{}^{c}F_{bc}^{},\widehat{R}_{00}=\frac{1}{4}F_{ab}F^{ab},\widehat{R}_{0a}=\frac{1}{2}^bF_{ab}$$
(81)
where $`F=dA`$, and $`\widehat{e}^0=d\psi +A`$, $`\widehat{e}^a=e^a`$. Taking $`F_{ab}=\mu J_{ab}`$, where $`J_{ab}`$ is the Kähler form on $`M_n`$, we therefore have

$$\widehat{R}_{ab}=(\mathrm{\Lambda }\frac{1}{2}\mu ^2)\widehat{g}_{ab},\widehat{R}_{00}=\frac{1}{4}n\mu ^2,\widehat{R}_{0a}=0,$$
(82)
where $`R_{ab}=\mathrm{\Lambda }g_{ab}`$ in $`M_n`$, and so the $`(n+1)`$-dimensional metric $`d\widehat{s}^2`$ will be Einstein, $`\widehat{R}_{AB}=\widehat{\mathrm{\Lambda }}\widehat{g}_{AB}`$, and

$$\widehat{\mathrm{\Lambda }}=n,\mathrm{\Lambda }=n+2,$$
(83)
provided that we take $`\mu =2`$.

The covariant exterior derivative on spinors, $`\widehat{D}d+\frac{1}{4}\widehat{\omega }_{AB}\mathrm{\Gamma }^{AB}`$, is easily seen to be given by

$$\widehat{D}=D\frac{1}{4}AJ_{ab}\mathrm{\Gamma }^{ab}+\frac{1}{2}J_{ab}e^b\mathrm{\Gamma }^{0a}\frac{1}{4}J_{ab}d\psi \mathrm{\Gamma }^{ab},$$
(84)
where $`D=d+\frac{1}{4}\omega _{ab}\mathrm{\Gamma }^{ab}`$ is the covariant exterior derivative on spinors in the base space $`M_n`$. Note that since $`n`$ is necessarily even, the spinors in the total space $`\widehat{M}`$ have the same dimension as those in base space $`M_n`$, and so we do not need to make any tensor-product decomposition of the Dirac matrices. The equation for Killing spinors in the $`(n+1)`$-dimensional bundle space $`\widehat{M}`$, in the normalisation $`\widehat{R}_{AB}=n\widehat{g}_{AB}`$ that we established above, is simply $`\widehat{D}_A\eta =\frac{1}{2}\mathrm{i}\sigma \mathrm{\Gamma }_A\eta `$, where $`\sigma =\pm 1`$. From (84), this gives the equations

$`D_a\eta A_a{\displaystyle \frac{\eta }{\psi }}`$ $`=`$ $`\frac{1}{2}J_{ab}\mathrm{\Gamma }^{0b}\eta +\frac{1}{2}\mathrm{i}\sigma \mathrm{\Gamma }_a\eta ,`$
$`{\displaystyle \frac{\eta }{\psi }}`$ $`=`$ $`\frac{1}{4}J_{ab}\mathrm{\Gamma }^{ab}\eta +\frac{1}{2}\mathrm{i}\sigma \mathrm{\Gamma }_0\eta .`$ (85)
As is well known, the Einstein-Kähler space $`M_n`$ admits a gauge-covariantly constant spinor $`\epsilon `$, satisfying

$$𝒟_a\epsilon \mathrm{i}eA_a\epsilon =0,$$
(86)
where as above we have $`dA=F=2J`$, and $`e`$ is the electric charge carried by $`\epsilon `$. This can be determined by examining the integrability condition $`[𝒟_a,𝒟_b]\epsilon =\frac{1}{4}R_{abcd}\mathrm{\Gamma }^{cd}\epsilon 2\mathrm{i}eJ_{ab}\epsilon `$. Multiplying by $`\mathrm{\Gamma }^{ab}`$, this gives $`\mathrm{i}n(n+2)\epsilon =4eJ_{ab}\mathrm{\Gamma }^{ab}\epsilon `$. It is a straightforward exercise to calculate the eigenvalues of the matrix $`J_{ab}\mathrm{\Gamma }^{ab}`$, and to show, in particular, that in general it has only two singlet eigenvalues, which are $`\pm \mathrm{i}n`$. It is these singlets that are associated with the gauge-covariantly-constant spinor $`\epsilon `$ (and its charge conjugate), and so we can deduce that

$$e=\frac{1}{4}(n+2).$$
(87)
One can also then easily show that $`\mathrm{\Gamma }^0\epsilon =\sigma \epsilon `$, where $`\sigma =\pm 1`$. From the second equation in (85) we therefore deduce that if we take $`\eta =f(\psi )\epsilon `$ we shall have

$$f=e^{\frac{1}{4}(n+2)\mathrm{i}\psi },$$
(88)
and then the first equation in (85) confirms that indeed $`\epsilon `$ satisfies

$$D_a\epsilon \frac{1}{4}\mathrm{i}(n+2)A_a\epsilon =0.$$
(89)
In other words, we have proved that if $`\epsilon `$ is the gauge-covariantly constant spinor in the Einstein-Kähler manifold $`M_n`$, then $`\eta =e^{\frac{1}{4}(n+2)\mathrm{i}\psi }\epsilon `$ is a Killing spinor in the $`U(1)`$ bundle over $`M_n`$, which is therefore an Einstein-Sasaki manifold $`\widehat{M}`$. The conjugate spinor satisfies the Killing-spinor equation with the opposite sign on the right-hand side. Lifted up further using (77), one obtains the conjugate pair of covariantly-constant spinors in the Ricci-flat Kähler cone over the Einstein-Sasaki manifold. The situation is summarised in the table below.

$`2m`$-Dimensional $`(2m+1)`$-Dimensional $`(2m+2)`$-Dimensional  Einstein-Kähler Einstein-Sasaki Calabi-Yau Cone    $`\mathrm{\Lambda }>0`$ $`\mathrm{\Lambda }>0`$ $`R_{ab}=0`$  $`(D_a\mathrm{i}eA_a)\epsilon =0`$ $`D_a\eta =\pm \mathrm{i}m\mathrm{\Gamma }_a\eta `$ $`D_a\eta =0`$

The progression from Einstein-Kähler to Einstein-Sasaki to Ricci-flat Kähler cone.

### 4.3 Lichnerowicz bound for Einstein-Sasaki spaces

In any Einstein space $`M`$ that admits Killing spinors, we can prove that the bound (21) that governs the stability of AdS$`\times M`$ solutions, and also the stability of Schwarzschild-Tangherlini black holes, is always satisfied. In other words, we can prove that an Einstein space in $`d`$ dimensions with cosmological constant $`\mathrm{\Lambda }`$ has a Lichnerowicz spectrum such that

$$\mathrm{\Delta }_L\frac{\mathrm{\Lambda }}{d1}\left(4\frac{(5d)^2}{4}\right).$$
(90)
To prove this we shall first, for convenience, make our conventional choice of normalisation $`\mathrm{\Lambda }=d1`$. A Killing spinor therefore satisfies $`D_a\eta =\frac{1}{2}\mathrm{i}\mathrm{\Gamma }_a\eta `$. Suppose that $`h_{ab}`$ is a transverse traceless mode of the Lichnerowicz operator on $`M`$:

$$\mathrm{\Delta }_Lh_{ab}=\lambda h_{ab},^ah_{ab}=0,h^a{}_{a}{}^{}=0.$$
(91)
We now define two vector-spinors:

$$\varphi _ah_{ab}\mathrm{\Gamma }^b\eta ,\chi _a(_bh_{ac})\mathrm{\Gamma }^{bc}\eta .$$
(92)
The assumed properties of $`h_{ab}`$ can easily be seen to imply that

$$D^a\varphi _a=0,\mathrm{\Gamma }^a\varphi _a=0,D^a\chi _a=0,\mathrm{\Gamma }^a\chi _a=0.$$
(93)
We now calculate the action of the Rarita-Schwinger operator on the vector-spinors, finding after some algebra that

$$\mathrm{i}\mathrm{\Gamma }^bD_b\varphi _a=\chi _a\frac{\mathrm{i}}{2}(d2)\varphi _a,\mathrm{i}\mathrm{\Gamma }^bD_b\chi _a=(\lambda d)\varphi _a+\frac{\mathrm{i}}{2}(d4)\chi _a.$$
(94)
Thus by taking an appropriate linear combination of the two vector-spinors, we can form an eigenfunction $`\psi _a=\varphi _a+k\chi _a`$ of the Rarita-Schwinger operator on transverse gamma-traceless spin $`\frac{3}{2}`$ modes,

$$\mathrm{i}\mathrm{\Gamma }^bD_b\psi _a=\mu \psi _a.$$
(95)
It follows immediately from (94) that we shall have an eigenfunction if

$$\mu =\frac{1}{2}(d2)\mathrm{i}k(\lambda d),k\mu =\mathrm{i}\frac{1}{2}k(d4).$$
(96)
These equations determine the constant of proportionality to be $`k=\mathrm{i}/[\mu +\frac{1}{2}(d4)]`$, and hence that the Rarita-Schwinger eigenvalue $`\mu `$ satisfies

$$4\mu ^24\mu d^2+10d8=4\lambda .$$
(97)
Reorganising this we obtain

$$\lambda =\frac{1}{4}(2\mu 1)^2+4\frac{1}{4}(d5)^2.$$
(98)
From the reality of the Rarita-Schwinger eigenvalue $`\mu `$, we therefore deduce that

$$\lambda 4\frac{1}{4}(d5)^2.$$
(99)
Restoring the cosmological constant, we therefore obtain the claimed inequality (90), which must hold for any Einstein space of positive Ricci tensor that admits Killing spinors. In particular, this encompasses the case of all Einstein-Sasaki manifolds, in all odd dimensions.

It is worth remarking that the above proof is a generalisation of an argument that was used in  in the case of seven-dimensional Einstein-Sasaki manifolds. It was argued there that such a manifold $`M_7`$ could be used in order to obtain a supersymmetric solution AdS$`{}_{4}{}^{}\times M_7`$ of eleven-dimensional supergravity. Now it is known that eigenfunctions of the Lichnerowicz operator in the internal space give rise to scalar fields in the AdS spacetime. The supersymmetry of the background implies that these scalars must be members of supermultiplets, including fermions. Since the Kaluza-Klein reduction must necessarily give rise to real masses for the fermions, it follows that the masses of the bosons (making due allowance for the need to define mass carefully in AdS) must be real also. This translates into the statement  that the $`(mass)^2`$ of the scalars must respect the Breitenlohner-Freedman  bound for stability, and hence it follows that the spectrum of the Lichnerowicz operator must be bounded from below by the stability limit, as given in (21), for the case $`d=7`$. The same argument was used recently for AdS$`{}_{5}{}^{}\times M_5`$ compactifications in . Of course our general proof above can be seen to be essentially an extension of the supersymmetry argument of , since in fact the crucial ingredient was not really supersymmetry per se, but rather, the fact that the mass spectrum of scalar fields can be related to the mass spectrum of spin $`\frac{1}{2}`$ fields. Since the Lichnerowicz Laplacian is the mass operator for scalar fields, and the Rarita-Schwinger operator is the mass operator for spin $`\frac{1}{2}`$ fields, our demonstration above that the eigenfunctions of the two operators are related when there are Killing spinors can be seen to reduce to the supersymmetry argument in those special dimensions where supersymmetric AdS vacua can be found. Our argument above is much more general, however, since it dispenses with the excess baggage of supersymmetry, and the need to interpret mass in AdS backgrounds.

In  the Lichnerowicz bounds were investigated for 7-dimensional Einstein metrics on the spaces $`M(m,n)`$ and $`Q(k,\mathrm{},m)`$ which are $`U(1)`$ bundles over $`^2\times S^2`$ and $`S^2\times S^2\times S^2`$ respectively, with the integers specifying the winding numbers of the $`U(1)`$ fibres over the base components. It is was found that for the Einstein-Sasaki examples, namely $`M(3,2)`$ and $`Q(1,1,1)`$, the bound $`\mathrm{\Delta }_L\frac{1}{2}\mathrm{\Lambda }`$ in (90) is strictly exceeded. This has the interesting consequence that for a range of ratios $`m:n`$ or $`k:\mathrm{}:m`$ around the Einstein-Sasaki values, the stability bound is still satisfied despite the absence of supersymmetry . By contrast, it was shown recently in an analogous 5-dimensional calculation for the $`T^{p,q}`$ spaces with Einstein metrics that the bound $`\mathrm{\Delta }_L\mathrm{\Lambda }`$ in (90) is exactly saturated by the Einstein-Sasaki case $`T^{1,1}`$, and that all the non-supersymmetric $`pq`$ spaces have a Lichnerowicz mode lying strictly below the bound .

It is worth remarking that, in view of the equivalence of the criteria for black hole stability and AdS stability described in sections 2.2 and 2.3, we have the immediate consequence that Einstein-Sasaki manifolds will always give stable Schwarzschild-Tangherlini black holes.

A further consequence is that any Einstein metric whose Lichnerowicz spectrum does not respect the lower bound (90) cannot admit Killing spinors, and so it cannot give rise to supersymmetric backgrounds in any supergravity compactification. Examples include not only the case of product metrics, for which it has long been known that there exists a Lichnerowicz zero mode , but also cases such as the Bohm metrics whose negative Lichnerowicz modes we have demonstrated in this paper.

## 5 Lorentzian Bohm Metrics, Real Tunneling Geometries, and Counterexamples to the Cosmic Baldness Conjecture

In this section we discuss metrics obtained by analytic continuation of the Bohm metrics. These metrics, which provide generalisations of de Sitter spacetime as locally static solutions with cosmological horizons, have a number of applications. In particular, they provide counterexamples to a certain form of the Cosmic Baldness conjecture. Furthermore, the Riemannian Bohm metrics have a totally geodesic hypersurface. This allows them to be viewed as real tunnelling geometries for the creation of the Lorentzian Bohm metrics “from nothing.” We first review the geometry by discussing the case of the round $`S^5=\mathrm{Bohm}(2,2)_0`$.

### 5.1 Round $`S^5`$ and $`dS_5`$

The round metric on $`S^5`$ may be written as

$$ds^2=d\rho ^2+\mathrm{sin}^2\rho \left(d\theta ^2+\mathrm{sin}^2\theta d\varphi ^2\right)+\mathrm{cos}^2\rho \left(d\theta _{}^{}{}_{}{}^{2}+\mathrm{sin}^2\theta _{}^{}{}_{}{}^{2}d\varphi _{}^{}{}_{}{}^{2}\right).$$
(100)
This provides an isometric embedding into $`𝔼^6`$, with Cartesian coordinates denoted by the variables $`(X_1,X_2,X_3,X_{}^{}{}_{1}{}^{},X_{}^{}{}_{2}{}^{},X_{}^{}{}_{3}{}^{})`$, via

$$(X_1,X_2,X_3)=\mathrm{sin}\rho (\mathrm{sin}\theta \mathrm{cos}\varphi ,\mathrm{sin}\theta \mathrm{sin}\varphi ,\mathrm{cos}\theta ),$$
(101)
and

$$(X_{}^{}{}_{1}{}^{},X_{}^{}{}_{2}{}^{},X_{}^{}{}_{3}{}^{})=\mathrm{cos}\rho (\mathrm{sin}\theta ^{}\mathrm{cos}\varphi ^{},\mathrm{sin}\theta ^{}\mathrm{sin}\varphi ^{},\mathrm{cos}\theta ^{}).$$
(102)
Thus, as required, one has

$$X_1^2+X_2^2+X_3^2+X_{}^{}{}_{1}{}^{2}+X_{}^{}{}_{2}{}^{2}+X_{}^{}{}_{3}{}^{2}=1.$$
(103)
Note also, for later use, that

$$X_1^2+X_2^2+X_3^2=\mathrm{sin}^2\rho ,X_{}^{}{}_{1}{}^{2}+X_{}^{}{}_{2}{}^{2}+X_{}^{}{}_{3}{}^{2}=\mathrm{cos}^2\rho .$$
(104)
The range of $`\rho `$ is seen to be $`\rho [0,\frac{\pi }{2}]`$.

To get the locally static Lorentzian de Sitter solution $`dS_5`$, one can set the angle $`\varphi =\mathrm{i}t`$ with $`t`$ real. This means that $`X_2=\mathrm{i}T`$ with $`T`$ real, and the embedding is into (but not onto) the quadric

$$X_1^2T^2+X_3^2+X_{}^{}{}_{1}{}^{2}+X_{}^{}{}_{2}{}^{2}+X_{}^{}{}_{3}{}^{2}=1.$$
(105)
Equations (104) now become

$$X_1^2T^2+X_3^2=\mathrm{sin}^2\rho ,X_{}^{}{}_{1}{}^{2}+X_{}^{}{}_{2}{}^{2}+X_{}^{}{}_{3}{}^{2}=\mathrm{cos}^2\rho .$$
(106)
It is clear that there can be points on $`dS_5`$ for which $`X_1^2T^2+X_3^2`$ is negative, and therefore for which $`X_{}^{}{}_{1}{}^{2}+X_{}^{}{}_{2}{}^{2}+X_{}^{}{}_{3}{}^{2}`$ exceeds unity. It follows that we need to use a different parameterisation. We set $`\mathrm{sin}^2\rho =1b^2`$ and get the metric

$$ds^2=\frac{db^2}{1b^2}+(1b^2)\left(d\theta ^2\mathrm{sin}^2\theta dt^2\right)+b^2\left(d\theta _{}^{}{}_{}{}^{2}+\mathrm{sin}^2\theta ^{}d\varphi _{}^{}{}_{}{}^{2}\right).$$
(107)
In this metric $`b[0,\mathrm{})`$. However, $`b=1`$ is a coordinate singularity, and for $`b>1`$ the orbits of $`_t`$ are spacelike. We use (107) in the region $`0b<1`$, where $`_t`$ is timelike. There are Killing horizons of $`_t`$ at $`\theta =0`$ and $`\theta =\pi `$. The metric on the horizon is

$$ds^2=\frac{db^2}{1b^2}+b^2\left(d\theta _{}^{}{}_{}{}^{2}+\mathrm{sin}^2\theta ^{}d\varphi _{}^{}{}_{}{}^{2}\right),$$
(108)
which is the standard metric on $`S^3`$. It is important to note that in order to get all of the $`S^3`$ we need both of the copies of this metric with $`0b<1`$ that from each of the two values $`\theta =0`$ and $`\theta =\pi `$, which each cover half of the horizon. This is best seen from the embedding of (101) and (102). The cosmological event horizon is not unique, since one may act with the $`SO(5,1)`$ isometry group. This fact is connected with the observer dependence of the associated Hawking thermal radiation . The coordinate singularity at $`b=1`$ arises because the 2-dimensional orbit of $`SO(2,1)`$ on the $`dS_2`$ factor changes from being timelike to spacelike as it crosses the surface $`b=1`$.

It is interesting to observe that with the metric on $`S^5`$ written as in (100), the map $`\varphi \varphi `$ is an isometry, which fixes pointwise a separating totally-geodesic hypersurface $`\mathrm{\Sigma }`$, given by $`\varphi =0`$ and $`\varphi =\pi `$ (because $`\pi \pi `$ here). In terms of the embedding described earlier, $`\mathrm{\Sigma }`$ is the hypersurface $`X_2=0`$. The $`_2`$ isometry implies the vanishing of the second fundamental form on $`\mathrm{\Sigma }`$, i.e. $`K\frac{1}{2}_ng=0`$, and hence the totally-geodesic property. Here $`n`$ is the normal to $`\mathrm{\Sigma }`$. The metric on $`\mathrm{\Sigma }`$ is

$$ds^2=\frac{db^2}{1b^2}+(1b^2)d\theta ^2+b^2\left(d\theta _{}^{}{}_{}{}^{2}+\mathrm{sin}^2\theta ^{}d\varphi _{}^{}{}_{}{}^{2}\right),$$
(109)
which is in fact just the round metric on $`S^4`$. We get all of $`S^4`$ because we have two copies of this metric, each with $`\theta [0,\pi ]`$, corresponding to $`\varphi =0`$ and $`\varphi =\pi `$. Again, this is seen most immediately in terms of the embedding (101)and (102), where $`X_2=0`$ manifestly defines an $`S^4`$. Thus we have a real tunnelling geometry in the sense of . That is, we have a compact gravitational instanton with totally-geodesic boundary, such as one might use to approximate a proposed wave function for the universe. The Riemannian metric may be grafted onto the Lorentzian $`dS_5`$ metric (107) at $`t=0`$, where it is clear that (107) also has the same totally-geodesic hypersurface with metric(109) defined as the fixed point set of $`tt`$.

### 5.2 Lorentzian Bohm metrics

The setup of the previous subsection 5.1 generalises straightforwardly to Bohm metrics. An isometric embedding is no longer possible, because spheres are the only positive curvature Einstein metrics that may be embedded isometrically into Euclidean space of one dimension higher. However, the topological statements go through. The analytic continuation of the five dimensional Bohm metrics gives

$$ds^2=d\rho ^2+a^2(\rho )\left(d\theta ^2\mathrm{sin}^2\theta dt^2\right)+b^2(\rho )\left(d\theta _{}^{}{}_{}{}^{2}+\mathrm{sin}^2\theta ^{}d\varphi _{}^{}{}_{}{}^{2}\right).$$
(110)
Note that we could have analytically continued the second sphere instead. In some cases this gives two inequivalent Lorentzian metrics; we shall discuss making this interchange $`ab`$ below. Again, we have Killing horizons at $`\theta =0`$ and $`\theta =\pi `$. The range of $`\rho `$ depends on the specifics of the Bohm metric.

Consider first the cases where the corresponding Riemannian Bohm metric has topology $`S^5`$. The topology of the horizon is $`S^3`$, as it was for the round $`S^5`$ of the previous subsection 5.1. This follows from the fact that the metric functions $`a(\rho )`$ and $`b(\rho )`$ behave near the two endpoints $`\rho =0`$ and $`\rho =\rho _f`$ in the same way as $`\mathrm{sin}\rho `$ and $`\mathrm{cos}\rho `$ respectively behave near the endpoints $`\rho =0`$ and $`\rho =\frac{\pi }{2}`$ of the round $`S^5`$ metric. As we shall see shortly, the area $`A`$ of this cosmological event horizon is always less than in the round case. Interchanging $`a`$ and $`b`$ in the $`S^5`$ Bohm examples gives the same Lorentzian metric.

When the Riemannian Bohm metric has topology $`S^3\times S^2`$, an exchange of functions $`ab`$ in (110) will change the topology of the Lorentzian manifold, and in particular the topology of the event horizon. This is because in these cases $`a(\rho )`$ goes to zero at both endpoints, whilst $`b(\rho )`$ never goes to zero. For topological purposes, we may think of $`a`$ as behaving like $`\mathrm{sin}\rho `$ with endpoints $`\rho =0,\rho =\pi `$ and $`b`$ behaving like a constant function (just as in the “trivial” Bohm metric Bohm$`(2,2)_1`$, which is simply the product Einstein metric on $`S^3\times S^2`$). Thus we have two possibilities. The metric (110) has a horizon with topology $`S^1\times S^2`$. If we exchange $`a`$ and $`b`$, the horizon will have topology $`S^3`$. These topologies are seen in the same way as in the previous subsection 5.1, and as always we should take care to include the two values $`\theta =0`$ and $`\theta =\pi `$. The area of the horizons are

$$A_1=8\pi a^2𝑑t$$
(111)
and

$$A_2=8\pi b^2𝑑t,$$
(112)
respectively. Note that there is an extra factor of 2 because there are contributions from both $`\theta =0`$ and $`\theta =\pi `$. The arguments of the following sections suggest that these two quantities should be equal, and less than the horizon area for the de Sitter spacetime $`dS_5`$. These are nontrivial conditions on the functions $`a`$ and $`b`$.

The non-uniqueness of these cosmological horizons is reduced compared to the de Sitter case, because the relevant isometry group is now only $`SO(2,1)`$.

The $`_2`$ isometry of the previous subsection 5.1 is also present in the Bohm metrics, and therefore we recover a totally-geodesic submanifold. Thus one might consider using Bohm metrics in tunnelling calculations for the creation of a Lorentzian Bohm universe. In that application the number of negative modes of $`\mathrm{\Delta }_2`$ should be an odd number, so as to get an imaginary part for the free energy when one evaluates the functional integral. Usually one expects just one negative mode, and the contribution of instantons with more than one is often ignored.

Another interesting question is what is the volume of $`\mathrm{\Sigma }`$, the totally-geodesic boundary. For tunnelling geometries constructed from hyperbolic tunnelling manifolds  <sup>2</sup><sup>2</sup>2Note that $`\mathrm{\Delta }_2`$ has no negative modes in this case, and so one might worry about the tunnelling interpretation., $`\mathrm{Vol}(\mathrm{\Sigma })`$ is a measure of the complexity of $`\mathrm{\Sigma }`$ and it is possible to bound the volume of the tunnelling geometry in terms of the volume of the boundary such that larger complexity, i.e. larger $`\mathrm{Vol}(\mathrm{\Sigma })`$, means larger volume . In the case of positive scalar curvature, and with boundaries $`\mathrm{\Sigma }`$ of simple topology ( $`S^4`$ or $`S^2\times S^2`$ in our case), the notion of complexity is not relevant. However, it is still interesting to know how the volume of the manifold is related to the volume of the totally-geodesic boundary. Specifically, the volume of the five-dimensional manifold is

$$V=16\pi ^2a^2b^2𝑑t,$$
(113)
while

$`\mathrm{Vol}(\mathrm{\Sigma }_2)=8\pi ^2{\displaystyle a^2b𝑑t},`$
$`\mathrm{Vol}(\mathrm{\Sigma }_1)=8\pi ^2{\displaystyle ab^2𝑑t}.`$ (114)
The two different values for $`\mathrm{\Sigma }`$ correspond to interchanging $`a`$ and $`b`$ in the metric. In the $`S^5`$ cases these will be the same, but for the $`S^3\times S^2`$ cases they will be different. Some examples are illustrated in the following table. The table also collects information about the corresponding values of $`b_0`$ and $`t_{fin}`$.

In this table the results for the double cone are found analytically, and the Bohm$`(2,2)_0`$ metric (the round 5-sphere) may also be calculated analytically as a check on the numerics. The volumes decrease from the round sphere to the double cone, as expected from Bishop’s theorem. The volumes of the $`\mathrm{\Sigma }`$’s decrease with the volume in the $`S^5`$ cases, but not in the $`S^3\times S^2`$ cases.

Similar calculations may of course be done with higher-dimensional Bohm metrics, as illustrated in the following table. The generalisation of $`\mathrm{\Sigma }`$ to higher dimensions is slightly more involved, and will not be discussed here.

### 5.3 Cosmological event horizons

In the Riemannian metrics, the circle action on $`S^5`$ generated by $`_\varphi `$ rotates the $`X_1X_2`$ plane. The action has an $`S^3`$’s worth of fixed points for which $`X_1=0`$ and $`X_2=0`$, corresponding to $`\theta =0`$ and $`\theta =\pi `$. Because the reversal of $`\varphi `$ is also an isometry, we have in fact an $`O(2)`$ action, which allows the analytic continuation to a locally static, i.e. time-reversal invariant, metric with a hypersurface-orthogonal locally timelike Killing vector field with a Killing horizon. These are the cosmological horizons of the previous subsections. In such cases the Lorentzian metric may be written locally as

$$ds^2=U^2dt^2+g_{ij}dx^idx^j,$$
(115)
where $`g_{ij}`$ is the metric on the orbit space Q of $`_t`$, and the gravitational field equations imply in particular that

$$^2U=\mathrm{\Lambda }U.$$
(116)
The quantity $`U`$ vanishes on the horizon, and its normal derivative $`\frac{U}{n}`$ on the horizon $`U=0`$ is a constant, which is called the surface gravity $`\kappa `$. The period in imaginary time, i.e. the real period of $`\tau \mathrm{i}t`$, is $`2\pi /\kappa `$. If $`V`$ is the volume of the corresponding Riemannian manifold, and $`A`$ is the area of the event horizon in the Lorentzian manifold, one has

$$V=\frac{2\pi }{\kappa }_QU\sqrt{g}d^{n1}x.$$
(117)
The boundary $`Q`$ is the event horizon, where the orbits degenerate. For example, when the Bohm metric has topology $`S^5`$ the boundary $`Q`$ is a 4-ball, $`B^4`$, with boundary the even horizon $`S^3`$. The metric on $`Q`$ is in fact given by (109). However, the crucial difference is that we take only $`\varphi =0`$ to intersect all the orbits of $`_\varphi `$ once, whilst for $`\mathrm{\Sigma }`$ we needed to take both $`\varphi =0`$ and $`\varphi =\pi `$. The $`S^4`$ that we had before corresponded to gluing two copies of $`B^4`$s across their boundary $`S^3`$.

Integration of $`^2U`$ gives

$$\mathrm{\Lambda }_QU\sqrt{g}d^{n1}x=\kappa A,$$
(118)
whence

$$V\mathrm{\Lambda }=2\pi A.$$
(119)
This argument shows that when there are two possible inequivalent analytic continuations, such as for the Bohm metrics on $`S^3\times S^2`$, the horizon areas should be the same, $`A_1=A_2`$. This is illustrated in the following table, which also illustrates the relationship (119), showing that it works for the various topologies. The values of $`V`$ are repeated from the previous table.

Note that $`\kappa `$ cancels in (119), as it must since it depends on the normalisation of the length of the Killing field, which is arbitrary. This relation between area and volume is quite universal and holds for any Einstein metric admitting an $`O(2)`$ action. It allows us to relate the on-shell action to the area of the horizon, and hence to show that formally at least, the entropy $`S`$ is given by

$$S=\frac{1}{4}A,$$
(120)
just as in four dimensions. This general argument was first given in four dimensions in . Now, a theorem of Bishop  tells us that for fixed $`\mathrm{\Lambda }`$, the volume $`V`$ never exceeds its value for the round sphere, with equality only in the case of roundness. It follows that this area or entropy is always less than the area or entropy of the corresponding horizon in $`dS_5`$.

Since dynamically one expects the area to increase, and thermodynamically one expects the entropy to increase, there seem to be some physical grounds for believing that the static Lorentzian Bohm metrics are dynamically unstable. Indeed, one might conjecture that if they are perturbed slightly at some initial time, then they will evolve to an asymptotically de Sitter state, and that this evolution will be such that the area of the cosmological horizon increases monotonically from a value near, but smaller than, its value for the initial Bohm metric, to its value for the de Sitter spacetime. It should not be impossible to investigate this conjecture numerically.

### 5.4 Consequences for Cosmic Baldness

It has been conjectured for some time  (in four dimensions) that there should exist only one regular static solution of the Einstein equations with a cosmological constant that has only a single cosmological horizon (so that $`Q=S^{n2}`$). In fact one usually thinks of $`Q`$ as being topologically an $`(n1)`$-ball, as described above. This is a much stronger statement than that locally, within the event horizon of every observer or most observers, the metric will settle down to the static de Sitter form. In fact for generic initial data one cannot hope that the metric will settle down globally to the de Sitter state, as was originally made clear in  (see also  for a detailed discussion using the exact Lorentzian Taub-NUT metrics).

It is now clear that the Bohm metrics provide infinitely many counterexamples to the Cosmic Baldness conjecture in dimensions $`5n9`$. The situation in four dimensions remains unclear. It is still possible to believe an even stronger conjecture, the truth of which would imply the Cosmic Baldness conjecture, namely that there is only one Einstein metric on $`S^4`$. At present all that is known is that if there is another Einstein metric on $`S^4`$, then its volume must be less than that of the round metric by a factor of 3 , and that the magnitude of the Weyl tensor must exceed a certain threshold . This is interesting in the light of the fact that it is the magnitude of the Weyl tensor which appears to play a role in controlling the spectrum of the Lichnerowicz operator. It may perhaps suggest that any counterexample will have a negative mode of $`\mathrm{\Delta }_2`$.

Curiously, there are some proofs of a form of the Cosmic No-Hair conjecture in the literature , but these proofs require a smooth structure at future spacelike infinity $`I^+`$. It seems likely that in our examples, the future timelike infinity will not be of the sort envisaged in those proofs. It would be interesting to investigate this point further, but this would seem to require analytic formulae for $`a(\rho )`$ and $`b(\rho )`$. As mentioned in the previous sub-section it seems likely that these static metrics will be dynamically unstable, and they may well evolve into an asymptotically de Sitter-like state. If this is true then the main physical spirit of the No-Hair conjecture will hold, even though the letter of the Baldness conjecture is broken.

## 6 Negative modes and non-uniqueness of the Dirichlet problem

The existence of negative Lichnerowicz modes has been a major theme in this work. This section contains speculative comments on a generic connection between negative Lichnerowicz modes and the non-uniqueness of solutions to the Dirichlet problem. In particular, we argue that there will be infinitely many negative Lichnerowicz modes on the non-compact Ricci-flat Bohm metrics. The existence of $`L^2`$ negative modes for the Lichnerowicz operator for non-compact Ricci flat manifolds such as the Riemannian Schwarzschild solution first came to light when considering the negative specific heat of black holes. Since then a considerable literature has grown up, analysing that and related cases.

In the non-compact Ricci flat case, it seems that the general picture is as follows. One has a class of metrics on a manifold $`M_n`$ depending on some parameters, in the simplest case just one overall scaling parameter $`\mu `$, such as the mass in the Schwarzschild case. One asks whether this metric can fill in a given boundary $`\mathrm{\Sigma }_{n1}`$ that has a given metric $`h_{ij}`$. That is, one tries to solve the Dirichlet problem for the Einstein equations. In the four-dimensional Riemannian Schwarzschild case for example, the boundary is taken to be $`S^2\times S^1`$ with the product metric. This is specified by the radius $`R`$ of the 2-sphere and the period $`\beta `$ of the circle. Physically, we are putting a black hole in a spherical box of radius $`R`$, and fixing the temperature on the boundary of the box to equal $`T=\beta ^1`$. If the metric and the boundary data are to agree then we must have

$$8\pi \mu \sqrt{1\frac{2\mu }{R}}=\beta .$$
(121)
The number of solutions of this equation for $`\mu `$ depends on the ratio $`\beta /R`$ that specifies, up to a scale, the boundary metric $`h_{ij}`$. One finds that if the ratio is small there are two solutions for $`\mu `$. The Einstein action $`I`$ of the two solutions differs. The action for the smaller value of $`\mu `$ is the smaller. We shall refer to these two solutions as branches. If the ratio is large there are no solutions for $`\mu `$. At the critical value the two solutions for $`\mu `$ coincide, and give $`\mu =\frac{1}{3}R`$.

Now consider the operator $`\mathrm{\Delta }_2`$ (which equals $`\mathrm{\Delta }_L`$ in this Ricci-flat situation), subject to Dirichlet boundary conditions, which gives the Hessian of the action $`I`$. For a large box (in relation to the scale set by the temperature $`T=\beta ^1`$), $`\mathrm{\Delta }_2`$ has a single negative mode for the branch with the smaller value of $`\mu `$, and it a positive but no negative mode for the branch with the larger value. As one reaches the critical value, the two branches coincide and so do the two eigenvalues. At the critical point there is an eigenmode of $`\mathrm{\Delta }_L`$ with zero eigenvalue. In other words, at the critical point there is a marginally stable mode. Because the specific heat is given essentially by the Hessian of the action (i.e. the free energy) considered as a function of the boundary data, it changes sign at this value.

Next, consider what happens as the radius of the box increases to infinity, so that the ratio $`\beta /R`$ goes to zero while keeping $`\beta `$ a constant, and take the branch on which $`\mu `$ is smaller. In the limit, one finds that $`\mu \beta /(8\pi )`$. On the other branch, one finds that $`\mu \mathrm{}`$. Thus, in the limit $`R\mathrm{}`$ one gets a non-compact Ricci-flat Einstein manifold, namely the standard Schwarzschild solution, and by following the mode which first appears as a zero-mode at the critical value, one gets an $`L^2`$ negative mode for $`\mathrm{\Delta }_L`$ on that manifold. This process is illustrated at fixed temperature in the following figure.

The picture described above has been vindicated by detailed numerical calculations in this and related cases . For example, Hawking and Page  studied black holes in anti-de Sitter spacetime. The classical solution is the Kottler or Schwarzschild anti-de Sitter solution, an Einstein metric with negative scalar curvature $`4\mathrm{\Lambda }`$. The role of the radius $`R`$ is now played by the cosmological constant $`\mathrm{\Lambda }`$, and the manifolds considered are always non-compact, but the general picture is similar.

The arguments given above are heuristic rather than being completely rigorous, but they suggest the following generalisation. One considers a one-parameter family of Dirichlet problems for the Einstein equations. As the parameter varies one finds a discrete non-uniqueness, with more and more branches appearing, generically in pairs, and as each new branch appears a zero mode of $`\mathrm{\Delta }_2`$ occurs, which then splits into a pair of modes, one with positive eigenvalue and one with negative eigenvalue. In the limit that one gets a non-compact manifold, one should have found, on the correct branch, as many $`L^2`$ negative modes as the number of critical values one has passed.

An obvious example on which to try this argument is the noncompact metric of Bohm on $`^3\times S^2`$ , recently considered by Kol  and discussed above. In fact it exhibits a feature not seen previously, which is that even within the restricted framework of cohomogeneity one metrics, the Dirichlet problem may have infinitely many solutions.

These metrics are determined by a single (scale) parameter, which may be taken to be the radius $`b_0`$ of the 2-sphere bolt. Now the geometry of a boundary at some fixed value of the radius $`R`$ is given by the ratio $`a(R)/b(R)`$. Thus the possible filling solutions are determined by the intersections of the curve of $`a/b`$ with a straight line through the origin in the $`(a,b)`$ plane, as shown in figure 11. Clearly, if we let $`a/b`$ tend to one there are more and more intersections, which appear at critical points in pairs when the straight line touches the slightly wiggly but almost straight $`a/b`$ curve. It seems reasonable to suppose that an additional zero-mode of $`\mathrm{\Delta }_L`$ appears at this point, and then as the slope of the straight line gets closer to unity, a pair of eigenvalues, one positive and one negative, branches off. If this intuition is correct, and if the branches are connected at $`a/b=1`$ say, then one expects that the noncompact Bohm metrics on $`^3\times S^3`$ should have infinitely many $`L^2`$ Lichnerowicz negative modes. This process is illustrated in the following figure.

Some evidence for this viewpoint comes from examining $`SO(p)\times SO(q)`$-invariant transverse tracefree perturbations of the singular Ricci-flat cone on $`S^p\times S^q`$ obtained in the limit $`b_00`$. That is to say, the cone is

$$ds^2=dr^2+\frac{r^2}{p+q1}\left((p1)d\mathrm{\Omega }_p^2+(q1)d\stackrel{~}{\mathrm{\Omega }}_q^2\right),$$
(122)
and the perturbation is

$$h_{\alpha \beta }=r^2\varphi (r)\left(\begin{array}{cc}\frac{1}{p}g_p& \\ & \frac{1}{q}\stackrel{~}{g}_q\end{array}\right)_{\alpha \beta },$$
(123)
where $`g_p`$ and $`\stackrel{~}{g}_q`$ are the round metrics on $`S^p`$ and $`S^q`$, respectively. This is a zero mode on $`S^p\times S^q`$. If we want a mode on the cone with Lichnerowicz eigenvalue $`\lambda `$, the equation for $`\varphi `$ is

$$\frac{d^2\varphi }{dr^2}+\frac{d}{r}\frac{d\varphi }{dr}+\frac{2d2}{r^2}\varphi =\lambda \varphi .$$
(124)
The solutions to this equation for $`\lambda 0`$ are, writing $`d=p+q`$,

$$\varphi (r)=r^{(1d)/2}\left[AI_\mu (\lambda ^{1/2}r)+BK_\mu (\lambda ^{1/2}r)\right],\mu =\frac{1}{2}\sqrt{(d1)(d9)},$$
(125)
where $`A`$ and $`B`$ are constants, and $`I_\mu `$ and $`K_\mu `$ are the modified Bessel functions. When $`\lambda =0`$ we get

$$\varphi (r)=r^{(1d)/2\pm i\sqrt{(d1)(9d)}/2}.$$
(126)
This expression was written down in  (his $`d`$ is shifted by one). It is also the behaviour of the Bessel function solutions in (125) as $`r0`$. As $`r\mathrm{}`$, the Bessel function solutions go as $`r^{d/2}e^{\pm r}`$. The $`K_\mu `$ Bessel function is the better behaved.

We are interested in the Bohm cases $`4<d<9`$ that coincide with oscillatory behaviour in (125) and (126). For these dimensions, the zero-mode solutions are not normalisable as $`r\mathrm{}`$ or as $`r0`$, although there are bounded as $`r\mathrm{}`$. The negative mode solution with $`K_\mu `$ is normalisable at infinity.

Thus the $`K_\mu `$ solutions decay at infinity, and as we move in towards the origin, start oscillating at $`r\lambda ^{1/2}`$. In the singular cone limit, the modes are not normalisable at the origin. However, suppose we are in a rounded-off cone. For $`rb_0`$ the metric is essentially that of the singular cone, and we may use our solutions (125). If further we have $`\lambda ^{1/2}b_0`$, logarithmic oscillations will set in, within this asymptotic regime. We should then expect to be able to match this solution to a solution in the inner regions that is well-behaved at the origin (cf. ), for a certain discrete set of values for $`\lambda `$. This will give us a spectrum of negative Lichnerowicz modes. It would seem that there will be an infinity of such modes, accumulating at zero.

## 7 Conclusions and discussion

The principal focus of this paper has been to study applications of the countable infinities of inhomogeneous Einstein metrics on certain spheres and products of spheres, which were discovered recently by Bohm . These occur for the topologies $`S^{p+q+1}`$ and $`S^{p+1}\times S^q`$, for $`5p+q9`$ and $`p2`$, $`q2`$. They may be used in place of the usual round-sphere Einstein metrics in a variety of constructions including black holes and Freund-Rubin solutions, and after a Wick rotation to a Lorentzian section, they may be interpreted as spacetime metrics in their own right.

The stability of generalised Schwarzschild-Tangherlini black holes, where the $`d`$-dimensional constant-radius spatial sections $`M_d`$ in the $`(d+2)`$-dimensional spacetime are taken to be positive Ricci curvature Einstein spaces, was studied recently in . It was shown that a solution will be classically stable if the spectrum of eigenvalues of the Lichnerowicz operator on transverse traceless symmetric two-index tensors in $`M_d`$ is bounded below by a value corresponding to $`\mathrm{\Delta }\mathrm{stab}0`$ in (1). One of our results in this paper has been to show that this stability criterion is identical to one obtained in  for the stability of Freund-Rubin solutions AdS$`{}_{n}{}^{}\times M_d`$ of gravity coupled to a $`d`$-form (or $`n`$-form) field strength. Thus it becomes of considerable interest to try to obtain bounds on the spectrum of the Lichnerowicz operator on Einstein spaces $`M_d`$.

The Bohm metrics are sequences of cohomogeneity one Einstein metrics of the form (29), which more and more nearly approach a “double-cone” form (42) as one progresses along the sequence. One therefore intuitively expects that the “ballooning” instabilities associated with direct products of the sphere metrics forming the principal orbits will give rise to negative-eigenvalue modes of the Lichnerowicz operator. Such modes would a fortiori violate the stability criterion described above, implying that Schwarzschild-Tangherlini black holes or AdS$`{}_{n}{}^{}\times M_d`$ solutions constructed using the Bohm metrics would be unstable. In certain cases, including all the Bohm metrics on $`S^3\times S^2`$ and $`S^3\times S^3`$, we constructed analytic proofs that indeed show the existence of negative modes of the Lichnerowicz operator. Numerical calculations for other examples, namely Bohm metrics on $`S^5`$, have confirmed that these too have negative-eigenvalue modes of the Lichnerowicz operator. We believe that in fact all the Bohm metrics have negative Lichnerowicz modes.

One can perform analytic coordinate continuations in the Bohm metrics in order to obtain spacetimes with positive cosmological constant that generalise de Sitter spacetime. If one does this for the Bohm metrics that are themselves topologically spheres, then the resulting spacetimes have the same topology and global structure as de Sitter spacetime itself. These metrics provide infinitely many counterexamples, in dimensions $`5n9`$, to the Cosmic Baldness conjecture, which asserted the uniqueness of regular static solutions of the Einstein equations with a single cosmological horizon. However, although the Bohmian analogues of de Sitter spacetime are regular, we have argued that they are unstable and that they are likely to decay into a de Sitter-like state. This would mean that the No-Hair conjecture would remain inviolate.

In order to explore possible endpoints for the decay of spacetimes constructed using Bohm metrics, we were also led to consider other geometries for Einstein spaces $`M_d`$ that would satisfy the criteria for stability. In particular, we considered compact Einstein spaces of positive Ricci curvature that admit Killing spinors. We showed that in all such spaces, in any dimension, one can derive a lower bound on the spectrum of the Lichnerowicz operator which implies that the stability criterion $`\mathrm{\Delta }_{\mathrm{stab}}0`$ is satisfied. These examples include all the Einstein-Sasaki spaces, which may be defined as odd-dimensional Einstein spaces whose cones give Ricci-flat Kähler spaces in one higher dimension. It is straightforward to see that the covariantly-constant spinors on the Ricci-flat Kähler cone project down as Killing spinors on the Einstein-Sasaki base. A by-product of our results is that it demonstrates that the Bohm metrics, for which $`\mathrm{\Delta }_L`$ (and hence a fortiori $`\mathrm{\Delta }_{\mathrm{stab}}`$) can be negative, cannot admit Killing spinors.

There also exist Ricci-flat Bohm metrics, with non-compact topology. The structure of these metrics at short distance looks very similar to that near one of the two endpoints of the compact metrics. However, lacking the cosmological term that causes the metric functions to turn over and recollapse in the compact examples, the non-compact spaces asymptotically approach the cone over the an $`S^p\times S^q`$ direct-product base. The non-compact spaces have the topology $`^{p+1}\times S^q`$. The only parameter in the non-compact metrics is the overall scale. We constructed an analytic proof that the non-compact Bohm metrics for $`p=q=2`$ and $`p=2`$, $`q=3`$ have negative Lichnerowicz modes, and we presented a general argument that indicates that all the non-compact Bohm metrics will have infinitely many $`L^2`$ normalisable negative-eigenvalue Lichnerowicz modes.

Another possible use of the Bohm metrics is to construct four-dimensional gravitating monopoles and black holes by dimensional reduction, as studied in . This is possible because many of the Bohm metrics have $`S^3`$ factors. Recalling that $`S^3`$ is isomorphic to $`SU(2)`$, one can quotient by a $`U(1)`$ action on $`S^3`$ to end up with $`SU(2)/\left[U(1)\times _2\right]S^2`$. Thus the resulting lower-dimensional space will have an $`S^2`$ factor, i.e. it will look like a monopole or black hole, and it will come with a $`U(1)`$ gauge field. Explicitly, one can write the metric on $`S^3`$ using Euler angles as

$$ds^2=d\theta ^2+\mathrm{sin}^2\theta d\psi ^2+(d\varphi +\mathrm{cos}\theta d\psi )^2,$$
(127)
where $`0\theta \pi `$, $`0\varphi <2\pi `$ and $`0\psi 4\pi `$. The quotient by the $`_\varphi `$ isometry, and the $`_2`$ quotient $`\psi \psi +2\pi `$, leaves us with the standard metric on $`S^2`$, namely $`d\theta ^2+\mathrm{sin}^2\theta d\psi ^2`$, and a charge-two Dirac monopole $`A=\mathrm{cos}\theta d\psi `$. One can also quotient by the whole $`SU(2)`$, and in this case because there is no fibration over the $`S^3`$ the $`SU(2)`$ gauge fields obtained will be trivial.

Thus, for example, take the seven-dimensional noncompact Bohm metric over two copies of $`S^3`$, supplement the metric by an eighth timelike direction $`dt^2`$, and dimensionally reduce on $`SU(2)\times U(1)`$. One will obtain a gravitating $`U(1)`$ monopole with four scalar fields. Another possibility would be to take the generalised black hole in eight dimensions over a compact Bohm metric in six dimensions with topology $`S^3\times S^3`$. Again quotient by $`SU(2)\times U(1)`$, where the $`SU(2)`$ is acting on the round $`S^3`$ in the Bohm metric. We will obtain a $`U(1)`$ magnetically charged black hole in four dimensions. Because the $`S^3`$ that the $`U(1)`$ was acting on is not round, the black hole will not be spherically symmetric.

## Acknowledgements

We should like to thank Christoph Bohm for useful discussions. C.N.P. is grateful to the Isaac Newton Institute for Mathematical Sciences and the Centre for Mathematical Studies in Cambridge, and CERN, Geneva, for hospitality and financial support during the course of this work. C.N.P. is supported in part by DOE DE-FG03-95ER40917. G.W.G. acknowledges partial support from PPARC through SPG#613. S.A.H. is funded by the Sims scholarship.

APPENDICES

## Appendix A Numerical Solutions for Bohm Metrics

### A.1 Numerical techniques

The Einstein equations (32) and (33) cannot be solved explicitly when $`p>1`$ and $`q>1`$. It was shown in  that countable infinities of smooth solutions satisfying the $`S^{p+q+1}`$ or $`S^{p+1}\times S^q`$ boundary conditions (40) or (41) exist. It is quite a straightforward matter to obtain these solutions by numerical methods, since it turns out that the second-order equations (32) are quite stable.

We have constructed numerical solutions by first obtaining Taylor-series expansions for the metric functions $`a`$ and $`b`$ near to $`t=0`$, imposing the $`t=0`$ boundary conditions given in (40) (or, equivalently, in (41)). To the first couple of orders, these give

$`a`$ $`=`$ $`t{\displaystyle \frac{q(q1)+b_0^2(p+q)(pq+1)}{6b_0^2p(p+1)}}t^3+O(t^5),`$
$`b`$ $`=`$ $`b_0+{\displaystyle \frac{q1b_0^2(p+q)}{2b_0(p+1)}}t^2+O(t^4).`$ (128)
Note that $`a`$ is an odd function of $`t`$, whilst $`b`$ is an even function.

Using a Taylor expansion of the form (128) (which we actually evaluated up to order $`t^9`$), we then set initial data just outside the singular point, for a very small positive value of $`t`$. These data are then evolved forward in $`t`$ numerically, using the second-order equations (32). The exercise then becomes a “shooting problem,” in which one seeks to adjust the one free initial parameter $`b_0`$ so as to achieve a smooth termination of the evolved data at a point $`t=t_f`$ where $`a`$ and $`b`$ satisfy one or other of the $`t=t_f`$ boundary conditions given in (40) or (41).

In cases where $`p=q`$, the numerical analysis is simpler, since the regular solutions are all symmetric under reflection about the midpoint $`t=t_c=t_f/2`$. Thus one can avoid the need to handle the integrations in the region near $`t=t_f`$ where one or other metric function is tending to zero. Instead, the shooting problem reduces to finding a $`b_0`$ for which either $`\dot{a}=\dot{b}=0`$ at some point $`t=t_c`$ (for the Bohm$`(p,p)_{2m}=S^{2p+1}`$ metrics), or else for which $`a=b`$ and $`\dot{a}=\dot{b}`$ (for the Bohm$`(p,p)_{2m+1}=S^{p+1}\times S^p`$ metrics).

It is known from the results in  that there is a countable infinity of values $`b_0`$ for which a regular termination at some $`t_f`$ occurs. The largest $`b_0`$ yielding a regular solution is $`b_0=1`$, leading to the standard unit $`S^{p+q+1}`$ metric (35), which is called Bohm$`(p,q)_0`$. The next value is $`b_0=\sqrt{(q1)/(p+q)}`$, giving the direct-product Einstein metric (38) on $`S^{p+1,q}`$ that we call Bohm$`(p,q)_1`$. There is then a monotonically-decreasing sequence of $`b_0`$ values, giving the Bohm$`(p,q)_n`$ sequence of Einstein metrics, alternating between terminating $`t_0`$ boundary conditions given by (40) and (41). The limit point of the sequence is $`b_0=0`$, giving the double-cone singular metric (42).

Plots for the five-dimensional Einstein metrics Bohm$`(2,2)_n`$ on $`S^5`$ and $`S^3\times S^2`$ are given below, for $`0n6`$, with $`b_0=(1,\frac{1}{2},0.253554255,0.117794,0.053054,0.023571,0.010503)`$. We also give a plot for the case of Bohm$`(2,3)_2`$, which is topologically $`S^6`$, to illustrate an example where there is no symmetry between $`a`$ an $`b`$. This has $`b_00.297647`$, and the endpoint is at $`t_f2.68296`$ (there is no natural significance to the midpoint of the radial coordinate range in the $`pq`$ examples).

A few isolated examples for other values of $`p`$ and $`q`$ are as follows. The Bohm$`(3,3)_2`$ metric on $`S^7`$ has $`b_00.3055210896`$, with the midpoint occurring at $`t=t_c1.34689859293`$. The Bohm$`(3,3)_3`$ metric on $`S^4\times S^3`$ has $`b_00.14291337`$ and $`t_c1.4691901856`$. The Bohm$`(4,4)_2`$ metric on $`S^9`$ has $`b_00.2851829`$ and $`t_c1.376730624`$, whilst the Bohm$`(4,4)_3`$ metric on $`S^5\times S^4`$ has $`b_00.09135`$ and $`t_c1.5099148`$.

We can also treat the analysis of the non-compact Bohm metrics described in section 3.4 in a similar fashion. Since these are Ricci-flat solutions of the Einstein equations, the terms involving the cosmological constant will be absent in (32) and (33), but otherwise all the formulae are analogous. The short-distance Taylor expansions (128) now become

$`a`$ $`=`$ $`t{\displaystyle \frac{q(q1)}{6b_0^2p(p+1)}}t^3+O(t^5),`$
$`b`$ $`=`$ $`b_0+{\displaystyle \frac{q1}{2b_0(p+1)}}t^2+O(t^4).`$ (129)
Using this, taken to order $`t^9`$, to set initial data just outside the $`S^q`$ bolt at $`t=0`$, we again performed numerical integrations. The plots for the functions $`a`$ and $`b`$ in the representative example $`p=q=2`$ are given at the end of this section.