\section{Introduction} \label{sec:intro}

Two--photon interferometry using down converted photons has recently
been extensively used for demonstrating violations of local
\cite{alley,om88,rt89,ci-e94,ps94,IPRA94,Ijosab95,Ipla-4} as well as
non--local \cite{wzm91} hidden--variables theories. However, in spite
of the recent improvement in the efficiencies of single--photon
detectors (close to 80\%) \cite{kw-eb93}, all experiments carried out
to date have had ten or more times fewer coincidence counts than
singles counts and this, in effect, meant a detection efficiency under
10\%. The reason for this lies in a directional uncertainty of signal
with respect to idler photons. On every ten or more detected signal
(idler) photons only one of conjugate idler (signal) photons finds
its way to the other detector. (Detectors in two--photon experiments
must have the same openings.) Therefore, all experiments carried
out so far relied only on coincidence counts and
herewith on additional assumptions---the \it no enhancement
assumption\/ \rm was the most important one---which were considered
very plausible. Then Santos devised \cite{San1,San2,San3} local
hidden--variable models which violate not only the low detection
loophole but also the \it no enhancement assumption\/\rm.
These models, as well as improvements in techniques resulted in an
interest into \it loophole--free\/ \rm experiments. In the past two
years several sophisticated proposals appeared which rely on detailed
elaborations of all experimental details.
\cite{ci-e94,Ijosab95,Ipla-4,Fry,San4,jones} The last three
proposals make use of maximal superpositions and require detection
efficiency of at least 83\%\ \cite{83}, while the first three
references make proposals for nonmaximal superpositions relying on
recent results \cite{eber93,Ipla,Inasa} which require only 67\%\
detection efficiency. In this paper we dispense with all
supplementary assumptions by proposing a feasible selection
method of doing a loophole--free Bell experiment which ideally
requires only 67\%\ detection efficiency and
can reach a realistic detection efficiency as high as 75\%. It is
shown that this means a feasible conclusive experiment with a
realistic visibility of 85\%. The method employs solid
angles of signal and idler photons (in a type--II down conversion
process) which differ five times from each other. This enables a
detection of more than 90\%\ of conjugate photons. We also
consider a method of combining unpolarized independent photons into
spin correlated pairs by means of non--spin observables.  The physics
of such a preparation of spin states by means of non--spin observables
can be paralleled with the polarization correlation between unpolarized
photons discovered by Pavi\v ci\'c \cite{IPRA94} and formulated for
classical light by Paul and Wegmann \cite{paul94}. The main difference
is that in the latter experiments photons cross each other's paths
while in the present proposal they do not. In the end we compare
Hardy's equalities with the Bell inequalities.

\section{The experiment}
\label{sec:exp}

A schematic representation of the experiment is shown in
Fig.~\ref{exp}. An ultra--short laser beam (of frequency $\omega_p$),
split by a beam splitter, simultaneously pumps up two type--II
crystals. We assume that they are beta barium borate (BBO) crystals.
In each type--II crystal the parametric down conversion produces
intersecting cones of mutually perpendicularly polarized signal
(extraordinary linearly polarized within \it e\/\rm--ray plane of the
BBO) and idler (ordinary linearly polarized within \it o\/\rm--ray
plane of the BBO) photons (of frequencies $\omega_e+\omega_o=\omega_p$).
\cite{shih1} Signal and idler photons can be of the same
frequencies ($\omega_e=\omega_o=\omega_p/2$) in which case the cones
are tangent to each other. When we tilt the crystal so as to increase
the angle between the pumping beam and the crystal axis of the BBO
(increasing $\omega_o$ and decreasing $\omega_e$ slightly) the cones
start to intersect each other [see inlet in Fig.~\ref{exp}]. Looking
only at polarization, we see the photons at the intersections of
the cones as entangled, because one cannot know which cone which
photon comes from. \cite{gar} We then entangle one photon from one
pair with one photon from the other pair by an interference of
the fourth order at a beam splitter. Each successful entanglement
(coincidence firing of detectors D1' and D2') selects (opens the
computer gates for their counts) the other two conjugate photons
into a Bell state.

In a real experiment one first has to make photons at the cone
intersections of each BBO indistinguishable, which means one has
to compensate for the finite bandwidths and different group
velocities inside the crystal, i.e., for transversal and
longitudinal (\it e\/\rm--photon
pulls ahead) walkoff effects. Half--wave plates (exchanging
retardation of \it e\/ \rm and \it o\/ \rm photons) and quartz plates
(being positive uniaxial crystal---BBO is negative) do the job.
\cite{ci-e94,shih2,shih3,Kwiat} Then, by rotating the crystal, one can
entangle the photons in a (non)maximal singlet--like or triplet--like
state. \cite{Kwiat} In our proposal we assume both, crystals and
plates, prepared so as to produce maximal singlet--like states. (It
is interesting that starting with two maximal triplet--like states
we arrive at the same final expressions for the probabilities;
Cf.~Ref.~\cite{ps94}.) We also assume that the intensity
of the laser pumping beam is reduced so that the probability
of having more than two down converted singlets in chosen solid
angles within the pumping time is
negligible. \cite{ci-e94} We stress here that we choose a
subpicosecond laser since without such an ultra--short pumping
one would not be able to collect valid coincidence counts of D1'
and D2' simply because there are no detectors which could react
in a time short enough \cite{ps94} to confirm the intensity
interference between two \it independent\/ \rm down converted
photons (from two crystals) whose coherence time lies in a
subpicosecond region. Two successive pumping can take place within
several nanoseconds as determined by the lowest available detector
resolution (recovering) time. For the feasibility of the experiment
it is crucial that the probability of both photon pairs
coming from only one of the two crystals can be made sufficiently
small in comparison with the probability of photon pairs coming
from both crystals by using more and more asymmetric beam splitter
which at the same time lowers the required detection efficiency
more and more towards 67\%. We show this at the end of this section.

As we mentioned in Sec.~\ref{sec:intro}, the main detection
efficiency problem in two photon interferometry is that signal
and idler photons have
to be in equal solid angles and that therefore less than 10\%\
of conjugate photons reach a detector. The present set--up enables
us to use different solid angles for selecting photons (those which
interfere at the beam splitter) and their conjugate photons (whose
counts are passed by the gates). For the purpose, one has to
evaluate the angular width of the conjugate photons once we know
central directions (cone intersections) of both photons.
One can show that the angular width of a conjugate photon depends
on the frequencies of the pump, signal, and idler photons,
on the band widths, on the pump, signal, and idler group refraction
indexes, and on the directions of signal and idler photons with
respect to the pumping beam, but for all combinations of these
terms, the ratio of 1:5 between the solid angle of the
photons we detect by detectors D1' or D2' and the solid angle
centered around central direction of the conjugate photons ($ph$
in  Fig.~\ref{exp}) assures that over 90\%\ of conjugate photons
will be found in the latter solid angle \cite{teich} and that the
probability of detecting ``third party'' photons will be negligible.
Let us now dwell on deriving our probabilities.

We can have three input states, depending on whether the two
pairs come from different crystals or both of them from only
one of the crystals. The probabilities of the pairs being
emitted in any of these three possible ways are equal.
The singlets coming from the crystals are mutually
independent and we therefore formally describe them by
tensor products. The former one, coming from different
crystals is given by
\begin{eqnarray}
|\Psi\rangle={1\over\sqrt2}\bigl(|1_x\rangle_1|1_y\rangle_{1'}\>-
\>\,|1_y\rangle_1|1_x\rangle_{1'}\bigr)
\otimes{1\over\sqrt2}\bigl(|1_x\rangle_2|1_y\rangle_{2'}\>-
\>\,|1_y\rangle_2|1_x\rangle_{2'}\bigr)\,,\label{eq:4-state}
\end{eqnarray}
where $x$ corresponds to the \it e\/\rm--ray planes of the BBO's
and $y$ to the \it o\/\rm--ray planes. The latter ones, coming from
the same crystals are given by
\begin{eqnarray}
|\Psi_{20}\rangle={1\over2}\bigl(|1_x\rangle_1|
1_y\rangle_{1'}-|1_y\rangle_1|1_x\rangle_{1'}\bigr)^2
|0\rangle_{2'}\,,\label{eq:2x0state}
\end{eqnarray}
\begin{eqnarray}
|\Psi_{02}\rangle={1\over2}\bigl(|1_x\rangle_2|1_y\rangle_{2'}-
|1_y\rangle_2|1_x\rangle_{2'}\bigr)^2
|0\rangle_{1'}\,.\label{eq:0x2state}
\end{eqnarray}

To obtain the four photon coincidence probabilities we cannot
superpose these three input states upon one another because
that would violate the principle of indistinguishability. To see
this let us for the moment assume that our detection
efficiency be ideal (100\%) and that the polarizers
P1, P1$^\perp$, P2, and P2$^\perp$ be removed.
Then, with the help of the responses of the detectors D1 and
D2 we could tell $|\Psi\rangle$ (both D1 and D2 would
fire), from $|\Psi_{02}\rangle$ (only D1 would fire) or
from $|\Psi_{20}\rangle$ (only D2 would fire).
If detections in reality had been ideal we would have
used only $|\Psi\rangle$. Since they are not, we have to take
$|\Psi_{02}\rangle$ and $|\Psi_{20}\rangle$ into account as
well, but helpfully it turns out that the Bell inequality
containing their corresponding counts (in addition to
$|\Psi\rangle$ counts) is still violated.
Therefore, we start with $|\Psi\rangle$, i.e., with
two pairs coming out from different crystals, and discuss
$|\Psi_{02}\rangle$ and $|\Psi_{20}\rangle$ later on.
A multi--mode representation of the input state we give
later on.

The outgoing electric--field operators describing
photons---we call them \it selector photons\/\rm---which
pass through beam splitter BS, through polarizers P1' and P2'
(oriented at angles $\theta_{1'}$ and $\theta_{2'}$, respectively),
and are detected by detectors D1' and D2', read (see Fig.~\ref{exp})
\begin{eqnarray}
\hat E_{1'}&=&\left(\hat a_{1'x}t_x\cos\theta_{1'}+\hat
a_{1'y}t_y\sin\theta_{1'}\right)
e^{i\mbox{\scriptsize\bf k}_{1'}\cdot
\mbox{\scriptsize\bf r}_{1'}-
i\omega_1(t-t_{1'}-\tau_{1'})}\nonumber\\
&&+\ i\left(\hat a_{2'x}r_x\cos\theta_{1'}+\hat
a_{2'y}r_y\sin\theta_{1'}\right)
e^{i\tilde{\mbox{\scriptsize\bf k}}_{2'}\cdot
\mbox{\scriptsize\bf r}_{1'}-
i\omega_2(t-t_{2'}-\tau_{1'})}\,,\label{eq:E1'}
\end{eqnarray}
\begin{eqnarray}
\hat E_{2'}&=&\left(\hat a_{2'x}t_x\cos\theta_{2'}+\hat
a_{2'y}t_y\sin\theta_{2'}\right)
e^{i\mbox{\scriptsize\bf k}_{2'}\cdot
\mbox{\scriptsize\bf r}_{2'}-
i\omega_2(t-t_{2'}-\tau_{2'})}\nonumber\\
&&+\ i\left(\hat a_{1'x}r_x\cos\theta_{2'}+\hat
a_{1'y}r_y\sin\theta_{2'}\right)
e^{i\tilde{\mbox{\scriptsize\bf k}}_{1'}\cdot
\mbox{\scriptsize\bf r}_{2'}-
i\omega_1(t-t_{1'}-\tau_{2'})}\,,\label{eq:E2'}
\end{eqnarray}
where $t_x^2$, $t_y^2$ are transmittances, $r_x^2$, $r_y^2$ are
reflectances, $t_j'$ is time delay after which photon $j'$ reaches
BS, $\tau_{j'}$ is time delay between BS and Dj',
$\omega_{j'}$ is the frequency of photon $j'$,
$\mbox{\bf k}_{j'}$ is the wave vector of photon $j'$, and
$\tilde{\mbox{\bf k}}_{j'}$ is the wave vector corresponding
to $\mbox{\bf k}_{j'}$ after reflection at BS.
The annihilation operators act as follows:
${\hat a}_{jx}|1_{x}\rangle_{j'}=|0_{x}\rangle_{j'}$, \
${\hat a}_{jx}|0_{x}\rangle_{j'}=0$, $j'=1,2$.

Operators describing photons---we call them \it Bell
photons\/\rm---which pass through polarizers P1 and P2 (oriented
at angles $\theta_1$ and $\theta_2$, respectively) and are
detected by detectors D1 and D2 read
\begin{eqnarray}
\hat E_1=
(\hat a_{1x}\cos\theta_1+\hat a_{1y}\sin\theta_1)
e^{-i\omega_1t_1}\,,\label{eq:E1}
\end{eqnarray}
\begin{eqnarray}
\hat E_2=
(\hat a_{2x}\cos\theta_2+\hat a_{2y}\sin\theta_2)
e^{-i\omega_2t_2}\,.\label{eq:E2}
\end{eqnarray}

The probability of detecting all four photons by detectors D1,
D2, D1', and D2' is thus
\begin{eqnarray}
P(\theta_{1'},\theta_{2'},\theta_1,\theta_2)
&=&\eta^2\langle\Psi|\hat E_{2'}^\dagger\hat E_{1'}^\dagger\hat E_2^
\dagger\hat E_1^\dagger
\hat E_1^{\phantom\dagger}\hat E_2^{\phantom\dagger}
\hat E_{1'}^{\phantom\dagger}\hat E_{2'}^{\phantom\dagger}|
\Psi\rangle\nonumber\\
&=&{\eta^2\over4}(A^2+B^2-2AB\cos\phi)\,,\label{eq:prob-4}
\end{eqnarray}
where $\eta$ is detection efficiency;
$A=Q(t)_{11'}Q(t)_{22'}$
and $B=Q(r)_{12'}Q(r)_{21'}$; here
$Q(q)_{ij}=q_x\sin\theta_i\cos\theta_j-q_y\cos\theta_i
\sin\theta_j$; $\phi=(\mbox{\boldmath$k$}_1-
\tilde{\mbox{\boldmath$k$}}_2)\cdot\mbox{\boldmath$r$}_1+
(\mbox{\boldmath$k$}_2-\tilde{\mbox{\boldmath$k$}}_1)
\cdot\mbox{\boldmath$r$}_2+(\omega_1-\omega_2)(\tau_{1'}-\tau_{2'})$.

To obtain a realistic estimation of the above result we start
with the multi--mode input states
\begin{eqnarray}
|1\rangle_{1'}|1\rangle_{2'}=\int\int d\omega_1'
d\omega_1'\psi_1(\omega_1)\psi_2(\omega_2)|\omega_1\rangle_{1'}
|\omega_2\rangle_{2'}\,,\label{eq:int}
\end{eqnarray}
which we introduce into Eq.~(\ref{eq:4-state}); $\psi_i(\omega_i)$
($i=1,2$) are both peaked at $\omega={1\over2}\omega_p$:
$\omega_i'=\omega-\omega_i$ ($i=1,2$). We can keep the singlet
state description as given by Eq.~(\ref{eq:4-state}) because
it has recently been proved by Keller and Rubin \cite{rubin}---as
we briefly present below---that a subpicosecond pulse would give
practically the same output as the continuous pumping beams provided
a group velocity condition is matched. In doing so we
rely on the experimental and theoretical results obtained by Kwiat
\it et al.\/ \rm \cite{Kwiat}. We then make a
Fourier decomposition of the electric--field operators
[Eqs.~(\ref{eq:E1'}) and (\ref{eq:E2'})] and obtain the mean
value for $P(\theta_{1'},\theta_{2'},\theta_1,\theta_2)$.
By integrating the latter probability over $\tau_{1'}$, $\tau_{2'}$,
$\omega_1'$, and $\omega_2'$ and using
\begin{eqnarray}
{1\over T}\int_{-T/2}^{T/2}\cos(\omega\tau+a)\,d\tau=
{\sin(\omega T/2)\over\omega T/2}\cos a\,,
\qquad\qquad\qquad
\int_{-\infty}^{\infty}{\sin a\omega\over \omega}\sin b\omega
\, d \omega=0
\,,\label{eq:int1}
\end{eqnarray}
\begin{eqnarray}
\int_{-\infty}^{\infty}{\sin a\omega\over \omega}\cos b\omega
\, d \omega =
\left\{
\begin{array}{ll}
\pi & \mbox{for $b<a$}\\
\pi/2 & \mbox{for $b=a$} \\
0 & \mbox{for $b>a$}\\
\end{array}
\right.\label{eq:int2}
\end{eqnarray}
we obtain
\begin{eqnarray}
P(\theta_{1'},\theta_{2'},\theta_1,\theta_2)
={\eta^2\over4}(A^2+B^2-2ABv_e\cos\Phi)\,,\label{eq:prob-4e}
\end{eqnarray}
where
\begin{eqnarray}
v_e={\int_{-T/2}^{T/2}f_1(\tau-\tau_1)f_2(\tau-\tau_2)\,d\tau\over
\int_{-T/2}^{T/2}f_1^2(\tau-\tau_1)\,d\tau
+\int_{-T/2}^{T/2}f_2^2(\tau-\tau_2)\,d\tau}\,,\label{eq:ve}
\end{eqnarray}
where $f_i(\tau)=\int_{-\infty}^{\infty}\psi_i(\omega)\cos
\omega\tau\,d\omega$, ($i=0,1$), where $T$ is the detection time,
and where $\Phi=2\pi(z_2-z_1)/L$; here $L$ is the spacing of the
interference fringes,  $z_j$ are the coordinates of detectors D$j$'
along $\mbox{\boldmath$k$}_1-\tilde{\mbox{\boldmath$k$}}_2$ and
$\mbox{\boldmath$k$}_2-\tilde{\mbox{\boldmath$k$}}_1$ (see Fig.~1
in reference \cite{Ijosab95}); we dropped
primes from $\tau_{1'}$ and $\tau_{2'}$ for simplicity.
We see that $\Phi$ can be changed by moving the detectors
transversally to the incident beams.

By numerical calculation one can easily show that $v_e$ is not
susceptible to the variation of detection time $T$ provided
$|\tau_1-\tau_2|<|\omega_1-\omega_2|^{-1}$ (for $|\tau_1-
\tau_2|\ll|\omega_1-\omega_2|$ even when $T\gg |\omega_1-
\omega_2|^{-1}$). For $|\tau_1-\tau_2|\ll|\omega_1-\omega_2|^{-1}$ we
have $v_e\rightarrow1$, i.e., the sharpest fringes, and for
$|\tau_1-\tau_2|\gg|\omega_1-\omega_2|^{-1}$ we have $v_e\rightarrow0$
and the fringes disappear. With the experimentally reachable frequency
passband $\Delta\omega$ of the order of magnitude of THz within a
single parametric down conversion with a continuous pumping beam
reaching the condition $|\tau_1-\tau_2|\ll|\Delta\omega|^{-1}$ is not a
problem because the time interval between the idler and signal photons
is of the order of femtoseconds. In our case, when dealing with two
simultaneous down conversions we have to resort to an ultra--short
pumping beam to satisfy the condition. Pulse pumping beam shorter than
1 ps would in general ``make it possible to distinguish pairs of
photons born at sufficiently different depths inside the crystal with
a consequent decrease in two photon interference'' as recently shown
by Keller and Rubin. \cite{rubin} This happens when the center of
momentum of signal and idler photons and the center of pump pulse does
not leave the crystal simultaneously.  When they do, i.e., when, by
choosing appropriate material conditions and pump frequency for a down
conversion within a type--II crystal, we make ``the inverse of the
pump group velocity equals the mean of the idler and signal inverse
group velocity'' \cite{rubin} and therewith we make the photons
indistinguishable again. Singlets appearing from such ``compensated''
crystal therefore keep their description given in Kwiat \it et al.\/
\rm \cite{Kwiat} and that is what we rely on in the afore--given
calculation.

Another realistic detail of the experiment is that the
pinholes of detectors D1' and D2' are not points but have a certain
width $\Delta z$. Therefore, in order to obtain a realistic
probability we integrate Eq.~(\ref{eq:prob-4}) over $z_1$ and $z_2$
over $\Delta z$ to obtain \begin{eqnarray}
P(\theta_{1'},\theta_{2'},\theta_1,\theta_2)
&=&{\eta^2\over4}\int\limits_{z_1-{\Delta z\over 2}}
\limits^{z_1+{\Delta z\over 2}} \int\limits_{z_2-{\Delta
z\over2}}\limits^{z_2+{\Delta z\over 2}}
[A^2+B^2-2ABv_e\cos{2\pi(z_2-z_1)\over L}]dz_1dz_2\nonumber\\
&=&{\eta^2\over4}(A^2+B^2-v2AB\cos\Phi)\,,\label{eq:nu} \end{eqnarray}
where $v=v_e\bigl[\sin(\pi\Delta z/L)/(\pi\Delta z/L)\bigr]^2$ is the
\it visibility\/ \rm of the coincidence counting.

An analysis of Eq.~(\ref{eq:nu}) shows that triggering of
D1' and D2' by selector photons means that their conjugate
\rm Bell photons appear entangled in spite of the fact that they
stem from two independent sources (two crystals) and that
they do not interact in any way (e.g., they do not cross
each other's paths). In general, Bell photons are only
partially entangled as in the case of classical intensity
interferometry. For special cases, however, one can achieve full
quantum nonmaximal entanglement, i.e., one can obtain probability
zero for certain orientations of polarizers P1 and P2. In order
to obtain such an entangled state, which would at the same time
enable a violation of the Bell inequality with only 67\%\
detection  efficiency, it is necessary to use an asymmetrical
beam splitter, to orient polarizers P1' and P2', e.g., along
$\theta_{1'}=90^\circ$ and $\theta_{2'}=0^\circ$, and to put
detectors D1' and D2' in a symmetric position with respect to
BS and with respect to the photons paths from the middle of the
crystals so as to obtain $\Phi=0$. Eq.~(\ref{eq:prob-4}) then
projects out the following nonmaximal singlet--like probability:
\begin{eqnarray}
P(\theta_{1},\theta_{2})&=&
\eta^2s(\cos^2\theta_{1}\sin^2\theta_{2} -
2v\rho\sin{\theta_{1}}\cos\theta_{1}
\sin{\theta_{2}}\cos\theta_{2}\cos\Phi+
\rho^2\sin^2\theta_{1}\cos^2\theta_{2})\,,\nonumber\\
&\equiv&\eta^2p(\theta_{1},\theta_{2})\label{eq:r-prob}
\end{eqnarray}
where we assumed the near normal incidence at BS so as to have
$r_x^2=r_y^2=R$ and $t_x^2=t_y^2=T=1-R$, where we used
$s\equiv T^2/(R^2+T^2)$, $\rho\equiv R/T$, and where we multiplied
Eq.~(\ref{eq:nu}) by 4 for other three possible coincidence detections
[i.e., for ($\theta_{1'}^{\phantom{\perp}}$,$\theta_{2'}^\perp$),
($\theta_{1'}^\perp$,$\theta_{2'}^{\phantom{\perp}}$), and
($\theta_{1'}^\perp$,$\theta_{2'}^\perp$) which we do not
take into account because only
($\theta_{1'}$,$\theta_{2'}$)--triggering opens the gates] and by
$(R^2+T^2)^{-1}$ for photons emerging from the same side of BS
(which also do not open the gates).

The singles--probability of detecting a photon by D1  is
\begin{eqnarray}
P(\theta_1)=
\eta s(\cos^2\theta_1+\rho^2\sin^2\theta_1)\equiv \eta p(\theta_1)\,
.\label{eq:one}
\end{eqnarray}
Analogously, the singles--probability of detecting a photon by D2 is
\begin{eqnarray}
P(\theta_2)=
\eta s(\sin^2\theta_2+
\rho^2\cos^2\theta_2)\equiv \eta p(\theta_2)\,.\hfill\label{eq:one'}
\end{eqnarray}

Introducing the above obtained probabilities into the
Clauser--Horne \cite{cl-h} form of the Bell  inequality
\begin{eqnarray}
B_{CH}\equiv P(\theta_1,\theta_2)-P(\theta_1,\theta_2')+
P(\theta_1',\theta_2')+P(\theta_1',\theta_2)-
P(\theta_1')-P(\theta_2)\leq0\,,\label{eq:loop}
\end{eqnarray}
we obtain the following minimal efficiency for its violation
\begin{eqnarray}
\eta={p(\theta_1')+p(\theta_2)\over
p(\theta_1,\theta_2)-p(\theta_1,\theta_2')
+p(\theta_1',\theta_2')+p(\theta_1',\theta_2)}
\,.\hfill\label{eq:eta}
\end{eqnarray}
We stress here that the probabilities in Eq.~(\ref{eq:loop})
are proper probabilities---not the ratios of coincidence
counts as in the experiments carried out so far. For example,
$P(\theta_2)=\eta p(\theta_2)$ is a total number of counts
detected by detector D2 with the polarizer P2 oriented along
$\theta_2$---it is \it not\/ \rm either $\eta^2p(\infty,a_2)$ or
$\eta^2p(\infty,a_2)/p(\infty,\infty)$.

This efficiency is a function of visibility $v$ and by looking at
Eqs.~(\ref{eq:r-prob}), (\ref{eq:one}), and (\ref{eq:one'})
we see that for each particular $v$ a different set of angles
should minimize it. A computer optimization of angles---presented in
Fig.~\ref{vis}---shows that the lower the reflectivity is, the lower
is the minimal detection efficiency. Also, we see a  rather
unexpected property that a low visibility does not have a significant
impact on the violation of the Bell inequality. For example, with
70\%\ visibility and 0.2 reflectivity of the beam splitter we obtain
a violation of Eqs.~(\ref{eq:loop}) with a lower detection efficiency
than with 100\%\ visibility and 0.5 ($\rho=1$) reflectivity.

In Ref.~\cite{Ijosab95} we have shown that one can select
fully quantum entangled Bell photons even without polarizers
P1' and P2'; i.e., whenever unpolarized selector photons
trigger detectors D1' and D2' they open the gates for maximally
entangled singlet--like state of Bell photons. Now, it is of
interest to find out whether we can use such non--polarization
preparation to prepare full non--maximal polarization--entangled
states. To this aim, we calculate
\begin{eqnarray}
P_{\infty}(\theta_1,\theta_2)=
P(\theta_{1'},\theta_{2'},\theta_1,\theta_2)+
P(\theta_{1'}^\perp,\theta_{2'},\theta_1,\theta_2)+
P(\theta_{1'},\theta_{2'}^\perp,\theta_1,\theta_2)+
P(\theta_{1'}^\perp,\theta_{2'}^\perp,\theta_1,\theta_2)
\label{eq:infty}
\end{eqnarray}
where we obtain the last three probabilities by analogy
with the first one [Eq.~(\ref{eq:prob-4})]; e.g., in order to obtain
$P(\theta_{1'}^\perp,\theta_{2'},\theta_1,\theta_2)$, we
introduce $E_{2'}^\perp$ instead of $E_{2'}$ into Eq.~(\ref{eq:prob-4}),
where we get $E_{2'}^\perp$ from Eq.~(\ref{eq:E2'}) upon substituting
$-\sin\theta_{2'}$ for $\cos\theta_{2'}$ and  $\cos\theta_{2'}$ for
$\sin\theta_{2'}$. Eq.~(\ref{eq:infty}) yields
\begin{eqnarray}
P_{\infty}(\theta_1,\theta_2)=
\eta^2{(1-2r_x^2t_x^2)\sin^2\theta_{1}\sin^2\theta_{2} +
(1-2r_y^2t_y^2)\cos^2\theta_{1}\cos^2\theta_{2}+S-2vW\cos\Phi\over
2(1-2r_x^2t_x^2-2r_y^2t_y^2+t_x^2t_y^2+r_x^2r_y^2)}
\label{eq:infty-1}
\end{eqnarray}
where
\begin{eqnarray}
S=(t_x^2t_y^2+r_x^2r_y^2)(\sin^2\theta_{1}\cos^2\theta_{2}+
\cos^2\theta_{1}\sin^2\theta_{2})\,,
\label{eq:S1}
\end{eqnarray}
\begin{eqnarray}
W=(t_xr_x\sin\theta_{1}\sin\theta_{2}+
t_yr_y\cos\theta_{1}\cos\theta_{2})^2\,.
\label{eq:S2}
\end{eqnarray}
A computer calculation shows that this probability can violate
the Bell inequalities only for a detection efficiency 83\%\ or
higher. It also shows that the probability cannot be used for
obtaining Hardy's equalities \cite{hardy}. On the other hand,
an analysis of Eq.~(\ref{eq:infty-1}) shows that the only
way to obtain a non--partial, i.e., full quantum (non--classical)
entanglement is to use a symmetric beam splitter ($r_x^2=r_y^2=1/2$)
and a symmetric position of detectors D1' and D2' with respect to
BS and with respect to the photons paths from the middle of the
crystals so as to obtain $\Phi=0$. Under these conditions
Eq.~(\ref{eq:infty-1}) yields $P_{\infty}(\theta_1,\theta_2)=
{1\over 2}\sin^2(\theta_{1}-\theta_{2})$,
i.e., a maximal singlet--like state. Thus, by means of non--spin
preparation we can prepare only ``symmetric'' (maximal)
non--classical spin correlated states.

In the end we show that other down conversions which may
occur in the crystals and enter our statistics do not significantly
influence the obtained probabilities. The probability of both
photon pairs coming from only one of the two crystals
and the probability of their coming from both
crystals are of course equal, but for $\rho$ close to 0 the
influence of photon pairs coming from only one of the two crystals
can be made small enough for a conclusive Bell experiment. Let us
see this in detail.

Choosing $\theta_{1'}=90^\circ$, $\theta_{2'}=0^\circ$,
$\Phi=0$,  and rewriting the electric--field
operators [Eqs.~(\ref{eq:E1'}) and (\ref{eq:E2'})] accordingly, we
obtain the following probabilities of detecting the
``intruder'' counts (corresponding to both photons coming from the
same crystal and being both detected by D1 and D2, respectively) while
collecting the singles--probabilities [Eqs.~(\ref{eq:one}) and
(\ref{eq:one'})]
\begin{eqnarray}
P_{20}(90^\circ,0^\circ,\theta_1,-)=P_{20}(\theta_1)
=\eta s \rho (1+v)\sin^2(2\theta_{1})\,,\label{eq:20-prob}
\end{eqnarray}
\begin{eqnarray}
P_{02}(90^\circ,0^\circ,-,\theta_2)=P_{02}(\theta_2)
=\eta s \rho (1+v)\sin^2(2\theta_{2})\,.\label{eq:02-prob}
\end{eqnarray}
We could dispense with these counts only if detectors D1 and D2
could tell one photon from two. It is therefore important to see
whether the Bell inequality Eqs.~(\ref{eq:loop}) is still violated
when we have them in our statistics. In order to include them into
the Bell inequality we should add them to the singles--probabilities
given by Eqs.~(\ref{eq:one}) and (\ref{eq:one'}). By comparing
$P_{20}(\theta_1)$ and $P_{02}(\theta_2)$ with $P(\theta_1)$
and $P(\theta_2)$, respectively, we see that for the angles
close to $\pi\over2$ and $\pi$, for which the asymmetrical states
violate the Bell inequality, the following inequalities hold:
$P_{20}(\theta_1)\ll P(\theta_1)$ and
$P_{02}(\theta_2)\ll P(\theta_1)$.  For example, for
$\rho=0.1$, $\eta=0.75$, $v=0.9$, $\theta_1=104^\circ$, $\theta_1'=
89^\circ$, $\theta_2=181^\circ$, and $\theta_2'=161^\circ$
we obtain the violation $B_{CH}>0$. For the same parameters
we also obtain $B_{CH}-P_{20}-P_{02}>0$. However, this reduces the
value of $B_{CH}$ for which the Bell inequality is violated by 2/3.
On the other hand, we have to use birefringent polarizers P1 and P2
to be able to discard counts which fire both D1 and D1$^\perp$ when
collecting data for the singles--probability $P(\theta_1)$  by
D1 and those which fire both D2 and D2$^\perp$ when
collecting data for the singles--probability $P(\theta_2)$  by
D2. Therefore, in doing a real experiment we should better
split unwanted two--photon wave packets across
additional polarized beam splitters \cite{rt89} or, even better,
by applying \it photon chopping\/ \rm \cite{chop} when
collecting counts for singles--probabilities.
We stress here that by this method we do not affect
the conclusiveness of the Bell experiment but only pick out \it
valid\/ \rm Bell pairs from all those
ones already selected by the D1'--D2' coincidence gates.
That is, we do not discard any counts corresponding to firing of D1
and/or D2 by photons coming from different crystals.

\section{Conclusion}
\label{sec:concl}
Our elaboration shows that the proposed loophole--free Bell
experiment which selects two out of four photons into
nonmaximal singlet--like states can be carried out with the
present technique. The proposal makes use of an asymmetrical
preparation of two \it input\/ \rm photon singlets generated by
two nonlinear type--II crystals. The asymmetry consists in the
fact that we first let one photon from the first singlet interfere
in the fourth order with one photon from the other singlet at a
highly transparent beam splitter. Coincidental detections of the
photons interfering at the beam splitter (we call them \it
selector photons\/\rm) open gates for a
selection of the remaining two conjugate photons, one from each
singlet, into a new correlated state: nonmaximal \it Bell singlet\/\rm.
In other words, since no coincidence detection between signal and idler
photons of the input singlets is needed we can use several times
wider solid angles for the \it Bell photons\/ \rm than for
the \it selector photons\/\rm. With five times wider solid angle
(determined by pinholes \it ph\/ \rm in Fig.~\ref{exp}) we collect
practically all Bell companions of those \it selector
photons\/ \rm which trigger detectors D1' and D2' in coincidence.
In this way we eliminate the main cause of the low detection
in all two--photon experiments so far: loosing of the conjugate
photons (in most cases they ``miss'' the detector opening).
An apparent draw--back to our set--up is that the
probabilities of two pairs coming from both crystals and of
both pairs coming from only one of the crystals are equal.
However, the above calculations show that for reflectivity 0.1,
realizable visibility of 85--90\%, and achievable detector
efficiency of 75\%\ \cite{ci-e94,kw-eb93} the Bell
inequality is violated even when the counts corresponding to
photons emerging from only one of the two crystals are included
into the statistics by which the inequality is fed.

We should mention here that although 67\%\ efficiency result for
Hardy's equalities has been obtained recently as well
\cite{gar2,china} their low marginal violations (of maximal value
0.09 as opposed to 0.41 for the Bell inequalities) make a
loophole--free ``Hardy experiment'' more demanding. However, it
would be worth trying to collect data for it within the proposed
set--up because of its conceptual clarity and because our results
add to the physics of the Hardy experiment. In particular,
an analysis of Eq.~(\ref{eq:infty-1}) shows that Hardy's equalities,
as opposed to the Bell inequalities, cannot be formulated for a
system which is not fully non--classical.  Thus, our set--up
reveals nonlocality of quantum systems as a property of selection
of their subsystems and Hardy's equalities as a test (ideally, some
detectors should always react and some never) of whether the system
is fully quantum or not. It may turn out that quantum nonlocality is
only operationally defined in the same way in which quantum
phase might turn out to be only operationally defined.
\cite{mandel} On the other hand, since Hardy's equalities reach
their maximal violation for $R=0.32$ and not for $R=0.5$, it might
turn out the unwanted effect of both photon pairs coming from the
same crystal on the marginal probabilities can be compensated
sufficiently well to make the experiment feasible.

In the end we mention that the set--up may find its application in
quantum cryptography and quantum computation for its property to
deliver Einstein--Podolsky--Rosen
singlets \cite{cryp} whose ``coherence ... [is] retained over
considerable distances and for long times'' \cite{sten};
actually, our \it Bell singlets\/ \rm stay coherent for ever, i.e.,
until we make use of them and collapse them.

\bigskip
\bigskip
\parindent=0pt
\it Note.\/ \rm\ Some preliminary results to these paper have
been presented within an invited talk at the \sl Adriatico
Research Conference on\/ \it Quantum Interferometry II\/\rm,
held in Trieste, Italy, 4--8 March 1996. \cite{Trst}

\acknowledgments
The author thanks to Harry Paul, Humboldt--Universit\"at zu
Berlin for many valuable discussions. He also acknowledges
supports of the Max--Plack--Geselschaft, Germany.

\parindent=20pt


\section{#1}}
\renewcommand{\theequation}{\arabic{section}.\arabic{equation}}

\def\begin{equation}{\begin{equation}}
\def\end{equation}{\end{equation}}
\def\begin{eqnarray}\samepage{\begin{eqnarray}\samepage}
\def\end{eqnarray}{\end{eqnarray}}
\def\mathrm{\mathrm}
\def\mathcal{\mathcal}
\def\phantom{\phantom}
\def\stackrel{\mr{def}}{=}{\stackrel{\mathrm{def}}{=}}

\font\msbm=msbm10

\topmargin=-0.3in
\headsep=0.0in
\headheight=0.0in
\oddsidemargin=-0.25in
\evensidemargin=-0.25in
\textwidth=6.73in
\textheight=10.0in

\begin{document}

\title{{\bf Dynamic and Thermodynamic Stability and Negative Modes in
Schwarzschild-Anti-de Sitter}}

\author{\\Tim Prestidge \\ \\{\sl T.Prestidge@damtp.cam.ac.uk} \\ \\
\\Department of Applied Mathematics and Theoretical Physics,\\The
University of Cambridge,\\Silver Street,\\Cambridge,\\CB3 9EW. \\ \\
{\bf damtp -- 1999 -- 89}}

\maketitle

\begin{abstract}
The thermodynamic properties of Schwarzschild-anti-de Sitter black
holes confined within finite isothermal cavities are examined. In
contrast to the Schwarzschild case, the infinite cavity limit may be
taken which, if suitably stated, remains double valued. This allows
the correspondence between non-existence of negative modes for
classical solutions and local thermodynamic stability of the
equilibrium configuration of such solutions to be shown in a well
defined manner. This is not possible in the asymptotically flat
case. Furthermore, the non-existence of negative modes for the
larger black hole solution in Schwarzschild-anti-de Sitter provides
strong evidence in favour of the recent positive energy conjecture by
Horowitz and Myers.
\end{abstract}

\renewcommand{\thepage}{ }
\pagebreak
\renewcommand{\thepage}{\arabic{page}}
\setcounter{page}{1}

\sect{Introduction}
It is well known that if the free energy of a system has an imaginary
component, then that system is meta-stable and may decay via some
mechanism in a first-order phase transition. This effect is well
illustrated in, for example, the treatment of hot flat space in four
dimensions given by Gross, Perry, and Yaffe~\cite{GPY}. In their
analysis, the authors demonstrate that flat space at zero temperature
is stable both classically and quantum-mechanically since the flat
space metric $\eta$, which is a true minimum of the action, is the
unique classical solution about which any saddle-point expansion may
be made. In the case of flat space at some non-zero temperature $T$
however, this is no longer true, since there exist other classical
solutions with period $\beta = 1/T$, such as the Euclidean section of
the Schwarzschild solution, which is not a true minimum of the action.

The $n=4$ dimensions Euclidean Schwarzschild solution is seen to
possess a smooth -- i.e. monotonically decreasing in $r$ --
spherically symmetric metric perturbation which decreases the
Euclidean action. The existence of this unique nonconformal `negative
mode' demonstrates that the Schwarzschild solution is a saddle-point
of the action rather than a true minimum, and may be interpreted as a
quantum-mechanical instability of hot flat space since nearby
geometries contribute an imaginary component to the free energy $F = -
1/\beta \: \log Z$ of the system. This permits, at any non-zero
temperature, the spontaneous excitation of hot flat space over an
effective potential barrier and the subsequent nucleation of black
holes with a rate proportional to the imaginary component. This
process may be viewed as a phase transition through topologies with
Euler characteristic $\chi = 0$ to $\chi = 2$, and is purely quantum
in nature with no classical analogue.

The authors conclude, therefore, that asymptotically flat space at
any non-zero temperature $T$ is quantum mechanically unstable due to
the spontaneous nucleation of black holes with mass $M = 1 / 8 \pi
T$. This conclusion is essentially in agreement with the idea that the
canonical ensemble in hot flat space is pathological when gravitational
effects are taken into consideration, although in the form given by
Gross {\sl et al.} the nucleation process is understood to be
thermodynamically inconsistent.

The above results pertain to black hole nucleation in hot flat space
confined to a large ideal isothermal cavity -- essentially in the
limit of infinite volume. However, a far deeper understanding of black
hole thermodynamics and the nucleation process is afforded by instead
considering cavities with finite volume~\cite{SWH1, GP1}. In this
approach, due to York~\cite{YORK1}, the canonical ensemble is
constructed from elements which are ideal isothermal boxes of
invariant surface area $A_B = 4 \pi r_B^2$ and wall temperature
$T_B$. The Euclidean Einstein equations with canonical boundary
conditions then admit either zero black hole solutions if the product
$r_B T_B < \sqrt{27}/8 \pi$, or else two distinct black hole solutions
which are degenerate at equality. The larger of the black hole
solutions is locally thermodynamically stable -- i.e. has positive
heat capacity, and may also have negative free energy, but it does not
have mass $1/8 \pi T_B$. The mass of the smaller black hole does
approach $1/8 \pi T_B$ in the limit of large $r_B$, but it is
thermodynamically unstable and has positive free energy.

These observations allowed York to give a consistent interpretation of
the black hole nucleation process in a finite cavity. In a spontaneous
process such as that proposed the relevant thermodynamic potential --
the free energy in the canonical ensemble -- cannot increase, and so
it is the larger thermodynamically stable black hole which may form by
nucleation. The positive free energy of the smaller hole acts as the
effective potential barrier.

A further advantage of this approach to black hole thermodynamics is
that it demonstrates a clear equivalence between the requirements for
local thermodynamic stability of an equilibrium configuration in the
canonical ensemble and the requirements for local dynamic stability of
the classical solutions. This correspondence may be demonstrated using
the `reduced action' $I_{\star}$ derived by Whiting and
York~\cite{YORK2}. Whiting~\cite{WHITING1} shows explicitly that, at
equilibrium, the second derivative of the reduced action with respect
to the horizon radius is positive --  and so the action is a local
minimum rather than a saddle point --  only in the region where the
heat capacity is positive. This is when the cavity radius satisfies
the constraint $r_+ < r_B < 3 r_+ / 2$ where $r_+$ is the horizon
radius.

Analysis of the Euclidean Schwarzschild nonconformal negative mode is
extended to isothermal cavities of finite volume by
Allen~\cite{ALLEN}, who finds that the mode persists until the cavity
radius $r_B$ falls below some critical radius $r_{\mathrm{crit}} \sim 2.89
\, r_+/2$. The non-existence of a negative mode for a classical
solution is a clear indication of local dynamic stability with respect
to small perturbations of the background, and provides the same
criterion on the action as the second derivative of $I_{\star}$ with
respect to the horizon radius. It should likewise correspond to the
sign of the heat capacity, but there is clearly some small amount of
disagreement here. This is at least partly due to Allen's choice of
boundary conditions: invariant cavity radius $r_B$ rather than
invariant surface area $A_B$, but as York has pointed
out~\cite{YORK1}, this implies that the eigenvalue spectrum of the
negative mode is not calculated at an extremum of the action. In fact,
the boundary conditions for calculating this spectrum depend
explicitly on the magnitude of the metric perturbation in the interior
region. It would appear, therefore, that in a finite isothermal cavity
there can be no well defined way of relating the local dynamics of a
classical solution in terms of the eigenvalue spectrum with the local
thermodynamics of the equilibrium configuration. Of course, the
original calculation of the negative mode eigenvalue in the infinite
volume case is properly defined: there is no choice in the boundary
conditions. However this demonstrates only that a dynamic instability
exists when the heat capacity is negative, not that the instability
vanishes when the solution becomes thermodynamically stable. In the
limit $r_B \rightarrow \infty$ the radius above which this occurs
becomes infinite, and the solution with positive heat capacity does
not exist.

The characteristics of the canonical ensemble in hot flat space
confined to a finite cavity may be contrasted with those in the
presence of a negative cosmological constant. For this purpose, York's
thermodynamic analysis is extended to the case of asymptotically
anti-de Sitter black holes by Brown, Creighton and Mann~\cite{BCM},
who study their properties in both two and three spatial dimensions.
In the case of three spatial dimensions the results are qualitatively
similar to flat space: for a given invariant surface area $A_B$ there
will exist either zero or two black hole solutions to the Euclidean
Einstein equations with canonical boundary conditions. If the constant
boundary temperature $T_B$ is chosen such that two distinct solutions
exist, then the larger black hole will be thermodynamically stable and
the smaller black hole will be unstable.

A key difference between the two cases, however, lies in the behaviour
of the heat capacity. In Schwarzschild-anti-de Sitter this is positive
when the radius of the cavity satisfies the constraint $r_+ < r_B <
f(r_+)$, where $f(r_+) \rightarrow \infty$ as $r_+ \rightarrow b /
\sqrt{3}$ from below, with $f(r_+) = \infty$ for $r_+ > b /
\sqrt{3}$. This means that in the limit of the cavity becoming
infinite, and with an appropriate rescaling of the temperature on the
boundary wall, there will still exist two well-defined black hole
solutions to the Euclidean Einstein equations: one locally
thermodynamically stable and the other unstable. In fact, if suitably
stated, all of the physical characteristics of the solutions are
retained in this limit, and the situation recovered is formally
identical to that which is analysed by Hawking and Page~\cite{PAGE}.

In their approach, these authors identify the temperature fixed at
infinity with the inverse of the periodicity $\beta_{\star}$ in the
imaginary-time direction. The temperature then has a minimum value
$T_0$, below which there are no black hole solutions and above which
there are two. They conjecture that, like Schwarzschild, the
Schwarzschild-anti-de Sitter solution should have a single
nonconformal negative mode, but only when $r_+ < r_0$ where $r_0 = b /
\sqrt{3}$ -- the horizon radius of the unique black hole solution with
temperature $T_0$. They further conjecture that the negative mode
should `pass through zero' at $r_+ = r_0$, and that for solutions with
horizon radius greater than this, no nonconformal negative modes
should exist. This strongly suggests that the action of the larger
solution, should it be negative, is in fact a global minimum of the
action functional for asymptotically anti-de Sitter metrics
periodically identified in imaginary time.

The objective of this work is to confirm these conjectures by
explicitly calculating the eigenvalue spectrum for the nonconformal
negative mode in Schwarzschild-anti-de Sitter. Since the boundary of
the cavity is at infinity this calculation is inherently independent
of the magnitude of the metric perturbation. It will therefore provide
a well defined demonstration of the correspondence between the local
dynamic stability of a classical solution in terms of the existence of
a negative mode and the local thermodynamic stability of the
equilibrium configuration in the canonical ensemble. It would seem
that this is only possible when the cosmological constant is non-zero.

The paper proceeds as follows: the next section is devoted to a review
of the thermodynamic properties of both asymptotically flat and
asymptotically anti-de Sitter black holes confined within a finite
ideal isothermal cavity. For simplicity the analysis is quite
graphical, and attention is restricted to the case of three spatial
dimensions only. Section three briefly discusses the theory behind
one-loop contributions to the partition function in the path-integral
approach to quantum gravity, and defines the physical gauge
independent operator $G$ acting on transverse and trace-free metric
perturbations. The eigenvalue equation for this operator acting on
spherically symmetric modes reduces to a second order O.D.E. in the
radial coordinate. In section four this equation is obtained in a
general form suitable for asymptotically flat and asymptotically
anti-de Sitter solutions in any dimension $n \ge 4$. Section five
discusses the numerical solution of this equation for
Schwarzschild-anti-de Sitter in $n=4$ dimensions and demonstrates the
stated correspondence with the heat capacity. This is then extended to
higher dimensions in section six. The final section discusses the
relevance of these observations to more recent work involving
anti-de Sitter spaces, and in particular to the positive energy
conjecture for periodically identified anti-de Sitter suggested by
Horowitz and Myers~\cite{MYERS}. The non-existence of a negative mode
for the thermodynamically stable Schwarzschild-anti-de Sitter black
hole provides strong evidence in favour of this proposal.

Units in which $G = \hbar = c = k_B = 1$ are used throughout the
paper.

\sect{Thermodynamic Properties of Black Holes in Finite
Isothermal Cavities}
In three spatial dimensions, the Euclidean Schwarzschild solution
takes the well known form
\begin{equation}
ds^2 = V(r) \, d\tau^2 + \frac{1}{V(r)}\, dr^2 + r^2 \, d\Omega^2
\label{metric}
\end{equation}
where the metric function $V(r) = 1 - r_+ / r$, with $r_+ = 2
M$. Clearly this solution is positive-definite only for values of $r$
greater than $r_+$, and a coordinate singularity occurs in the
$(\tau,r)$ plane at $r = r_+$. The $\mathrm{S}^2$ at this point may be
included in the Euclidean section if, in the conventional manner, the
coordinate $\tau$ is identified periodically with a period of
\begin{equation}
\beta_{\star} = \left. \frac{4 \pi}{V'(r)} \right|_{r = r_+} = \: 4
\pi r_+.
\end{equation}
Since the killing vector $\partial / \partial \tau$ is naturally
normalised to 1 in the limit of large $r$, the temperature measured at
infinity may be formally identified with the inverse of this
period. The Tolman law then states that for any static
self-gravitating system in thermal equilibrium, a local observer at
rest will measure a local temperature $T$ which scales as
$g_{00}^{-1/2}$. In the present context, then, the constant of
proportionality is $T_{\infty} = 1 / \beta_{\star}$.

With this in mind, York~\cite{YORK1} defines the elements of the
canonical ensemble for hot flat space to be ideal isothermal cavities
of invariant surface area $A_B = 4 \pi r^2_B$ and wall temperature
$T_B$. One topologically regular solution to the Einstein equations
with these boundary conditions is hot flat space with a uniform
temperature of $T_B$ throughout the cavity. Another solution is the
Schwarzschild metric. If a black hole of horizon radius $r_+ < r_B$
does exist within the cavity then clearly, from the Tolman law, the
wall temperature must satisfy
\begin{equation}
T_B \stackrel{\mr{def}}{=} T(r_B) = (4 \pi r_+)^{-1} \: (1 - r_+ / r_B)^{-1/2}.
\label{sch_surface}
\end{equation}
This equation may be solved for $r_+$ in terms of the constants $r_B$
and $T_B$. If the product $r_B T_B < \sqrt{27}/8 \pi$ then there are
no real positive solutions for $r_+$. In this part of the
$(T_B,r_B)$ plane no black holes exist. If this inequality is not
satisfied then in general two distinct solutions exist, which join
smoothly at equality where $r_+ = 2 r_B / 3$. Figure~\ref{sch} shows
this curve as a constant temperature slice through the $(T_B,r_B,r_+)$
surface.

\begin{figure}
\begin{picture}(0,0)(0,0)
\put(84,36){\small $T_B \! \uparrow$}
\put(121,60){\small $T_B \! \downarrow$}
\put(200,240){\small $r_+ = r_B$}
\put(250,208){\small $r_+ = 8 r_B / 9$}
\put(230,140){\small $r_+ = 2 r_B / 3$}
\put(288,0){\small $r_B$}
\put(-9,265){\small $r_+$}
\put(32,165){\small Exterior Region}
\end{picture}
\centering\epsfig{file=sch.eps,width=10cm}
\caption{A constant temperature slice through the $(T_B,r_B,r_+)$
surface defined by (\ref{sch_surface}) for a Schwarzschild black hole
confined within an ideal isothermal cavity. Also shown are the lines
$r_+ = r_B$ and $2 r_B / 3$, between which the heat capacity $C_A$ is
positive. The broken line is $r_+ = 8 r_B / 9$, above which the free
energy is negative.}
\label{sch}
\end{figure}

The diagram is to be interpreted as follows: for some given wall
temperature $T_B$ the solution curve will appear as shown. If the
scale of $r_B$ along the horizontal axis is kept fixed, then for
higher wall temperatures the solution curve shifts to the left, and
for lower wall temperatures it shifts to the right. The turning point
of the curve, where $r_B T_B = \sqrt{27} / 8 \pi$, remains on the line
$r_+ = 2 r_B / 3$. It is clear therefore that for any cavity radius to
the right of this point there will be two distinct Schwarzschild black
hole solutions -- one above this line and one below it.

For any values of $r_B$ and $T_B$, the entropy of the black hole
solutions to (\ref{sch_surface}) is $S = \pi r_+^2$. The heat capacity
at constant surface area -- the analogue in this context of $C_V$ --
may therefore be calculated for any solution using
\begin{equation}
C_A \stackrel{\mr{def}}{=} T_B \left. \frac{\partial S}{\partial T_B}
\right|_{A_B} = \: - 2 \pi r_+^2 \: (1 - r_+^{\phantom{2}}/ r_B^{\phantom{2}}) \:
(1 - 3 r_+^{\phantom{2}}/ 2 r_B^{\phantom{2}})^{-1}.
\end{equation}
Hence, if $r_+ < r_B < 3 r_+ / 2$ the heat capacity is positive and
the equilibrium configuration is locally thermodynamically stable. As
figure~\ref{sch} shows, this is always the case for the larger black
hole solution.

A far deeper understanding of the thermodynamics of black holes
in finite cavities, and the associated roles of both classical and
non-classical geometries within the canonical ensemble, is afforded by
considering the `reduced action' $I_{\star}$ proposed by Whiting and
York~\cite{YORK2, WHITING1, WHITING2, YORK3}. This is the function
obtained from the Euclidean Einstein action by considering only static
and spherically symmetric geometries which: i) are smooth and have
$\chi=2$ topology; ii) have appropriate boundary data for the
canonical ensemble; and iii) obey the constraints of the Einstein
equations on a family of space-like hyper-surfaces which foliate the
manifold.

With these requirements the reduced action becomes
\begin{equation}
I_{\star} = \beta_B \, E - S
\end{equation}
where $\beta_B$ is the inverse wall temperature, $S$ is the classical
black hole entropy $\pi r_+^2$, and the energy is chosen to be $E =
r_B - r_B \, (1 - r_+ / r_B )^{1/2}$. With this definition, then, both
the energy and the action of hot flat space are zero.

It should be noted that the reduced action has the same form as the
action evaluated for a classical Schwarzschild solution, except that
the variables $\beta_B, r_B,$ and $r_+$ are in this case all
independent rather than being related through (\ref{sch_surface}).
With this function, therefore, it is possible to {\sl derive} both the
thermodynamic properties of the ensemble and the dynamic properties of
the classical solutions within the ensemble by varying with respect to
the single remaining degree of freedom $r_+$. Thus, for example, the
stationary points satisfying $\partial I_{\star} / \partial r_+ = 0$
yield $\beta_B = (4 \pi r_+) \: (1 - r_+ / r_B)^{1/2}$, in agreement
with the necessary classical requirement. Furthermore, at the
stationary points,
\begin{equation}
\frac{\partial^2 I_{\star}}{\partial r_+^2} = \frac{1}{4} \: (1 -
r^{\phantom{0}}_+ / r^{\phantom{0}}_B )^{-1} \: \beta_B^{-2} \: C^{\phantom{0}}_A,
\end{equation}
and so
\begin{equation}
\frac{\partial^2 I_{\star}}{\partial r_+^2} > 0 \Leftrightarrow C_A
> 0.
\end{equation}
The condition for local dynamic stability of the classical solution is
therefore wholly equivalent to the condition for local thermodynamic
stability of the equilibrium configuration, where the latter is shown
to occur only within the range $r_+ < r_B < 3 r_+ / 2$. Should the
partition function in fact be dominated by a classical solution,
i.e. should the action for this solution -- and therefore the free
energy -- be negative, then the equilibrium configuration will be
globally thermodynamically stable. However, this occurs only within
the sub-range $r_+ < r_B < 9 r_+ / 8$.

Following Whiting~\cite{WHITING1}, the preceding analysis may be
written in a gauge invariant form by instead varying the reduced
action with respect to a generalised function $\overline{\beta}(r_+)$.
This function is chosen such that $\partial S / \partial E \: |_{r_B}
= \overline{\beta}$, where the form of $E$ is left undefined.
Stationary points of $I_{\star}$ then occur at $\overline{\beta} =
\beta_B$, and the second variation with respect to $\overline{\beta}$
at the stationary points indicates local dynamic stability as before.

With the introduction of a negative cosmological constant $\Lambda$,
the characteristics of Euclidean solutions to the Einstein equations
are changed in a number of ways. The line element for these
solutions still takes the form (\ref{metric}), but the metric function
becomes $V(r) = 1 - 2M/r + r^2/b^2$ where $b^2 = - 3 / \Lambda$. This
is positive-definite only for values of $r$ greater than the horizon
radius $r_+$ -- the real positive root of $V(r)$, and the mass
parameter may be written in terms of this root as
\begin{equation}
M = \frac{1}{2} \: r^{\phantom{0}}_+ \: ( 1 + r_+^2/b^2_{\phantom{0}} ).
\label{mass_func}
\end{equation}
Again, there is a coordinate singularity in the $(\tau,r)$ plane at $r
= r_+$, and the $\mathrm{S}^2$ at this point may be included in the
Euclidean section if the period of $\tau$ is now fixed to be
\begin{equation}
\beta_{\star} = \frac{4 \pi b^2 r_+}{b_{\phantom{0}}^2 + 3 r_+^2}.
\end{equation}
In this case, however, the temperature measured at infinity may no
longer be identified with the inverse of the periodicity as it can in
asymptotically flat spaces. There is no natural normalisation of the
Killing vector $\partial / \partial \tau$ in the limit of large $r$
due to the $r^2$ term in the metric function, and the locally measured
temperature of any thermal state decreases to zero at infinity.
Nonetheless, it is still consistent to regard $1 / \beta_{\star}$ as
the constant of proportionality in the Tolman equation.

\begin{figure}
\begin{picture}(0,0)(0,0)
\put(87,24){\small $bT_B \! \uparrow$}
\put(119,60){\small $bT_B \! \downarrow$}
\put(200,240){\small $\rho_+ = \rho_B$}
\put(288,0){\small $\rho_B$}
\put(-9,265){\small $\rho_+$}
\put(265,52){\small $\rho_+ = C_0(\rho_B)$}
\put(265,91){\small $\rho_+ = F_0(\rho_B)$}
\put(-6,50){\small 1}
\put(-15,88){\small $\sqrt{3}$}
\put(32,165){\small Exterior Region}
\end{picture}
\centering\epsfig{file=sch_ads.eps,width=10cm}
\caption{A constant temperature slice through the
$(bT_B,\rho_B,\rho_+)$ surface defined by (\ref{sch_ads_surface}) for
a Schwarzschild black hole in anti-de Sitter space confined within an
ideal isothermal cavity. Also shown are the lines $\rho_+ = \rho_B$
and $C_0(\rho_B)$, between which the heat capacity $C_A$ is positive.
Note that $C_0(\rho_B) \rightarrow 1$ in the limit of large $\rho_B$.
The broken line is $\rho_+ = F_0(\rho_B)$, above which the free energy
is negative.}
\label{sch_ads}
\end{figure}

So, the elements of the canonical ensemble for hot anti-de Sitter
space may again be defined as ideal isothermal cavities of invariant
surface area $A_B = 4 \pi r_B^2$ and wall temperature $T_B$, but
should there be a black hole within the cavity, this quantity must
satisfy
\begin{equation}
T_B \stackrel{\mr{def}}{=} T(r_B) = \frac{b_{\phantom{0}}^2+3r_+^2}{4 \pi b^2 r_+} \:
r_B^{1/2} \: \bigl( r_B^{\phantom{0}} - r_+^{\phantom{0}} - r_+^3 / b^2_{\phantom{0}} +
r_B^3 / b^2_{\phantom{0}} \bigr)^{-1/2}.
\end{equation}
Actually, it is convenient when considering anti-de Sitter to
transform to a dimensionless radial coordinate $\rho$, where
\begin{equation}
\rho \stackrel{\mr{def}}{=} \sqrt{3} \, r/b.
\label{rho_def}
\end{equation}
Remaining factors of the constant $b$ may then be absorbed into the
definition of each quantity, although for clarity these are explicitly
written. The wall temperature then becomes
\begin{equation}
b T_B = \frac{3}{4 \pi} \: \frac{1 + \rho_+^2}{\rho_+} \: \rho_B^{1/2}
\: \bigl( \rho_B^3 + 3 \rho_B^{\phantom{0}} - \rho_+^3 - 3 \rho_+^{\phantom{)}}
\bigr)^{-1/2}
\label{sch_ads_surface}
\end{equation}
which may be solved for $\rho_+$ in terms of the constants $\rho_B$
and $bT_B$. Like the asymptotically flat case there will again exist
either zero or two distinct black hole solutions, governed by the
relative values of $\rho_B$ and $bT_B$, although the relation in this
case is not so clear. The two solutions are degenerate along the locus
$\rho_+ = C_0(\rho_B)$, defined implicitly through the real positive
root of
\begin{equation}
\rho_B^3 + 3\rho_B^{\phantom{3}} + \frac{1}{2} \: \left\{ \frac{ \rho_+^5 +
2 \rho_+^3 + 9 \rho_+^{\phantom{0}} }{ \rho_+^2 - 1} \right\} = 0.
\end{equation}
Accordingly, for a given cavity radius, the temperature at and above
which solutions occur may be written as
\begin{equation}
(bT_B)^2 = \frac{9}{8 \pi^2} \: \rho_B \: \frac{1-C_0^2}{C_0^3}.
\end{equation}
This behaviour is illustrated in figure~\ref{sch_ads}, which is to be
interpreted in the same way as the previous diagram for the
asymptotically flat case. It should be noted that the diagrams are
identical in the limit $\rho_B \ll 1$ -- i.e. close to the origin,
since the cosmological constant can have very little effect on the
physics within cavities of radius $r_B \ll b$. For larger cavities,
however, the differences are manifest.

The heat capacity at constant surface area is
\begin{equation}
b^{-2} C_A = \frac{4 \pi}{3} \: \rho_+^2 \: \left\{ -1 + \frac{ 2
\bigl[ \rho_+^3 - 3 \rho_+^{\phantom{0}} \bigr] - (3 \rho_+^2 - 1) \bigl[
\rho_B^3 + 3 \rho_B^{\phantom{0}} \bigr] }{ (\rho_+^2+1) \bigl[ \rho_+^3 + 3
\rho_+^{\phantom{0}} - \rho_B^3 - 3 \rho_B^{\phantom{0}} \bigr] } \right\}^{-1},
\end{equation}
which is negative below the line $\rho_+ = C_0(\rho_B)$ and positive
above it. As before, then, the equilibrium configuration of the larger
black hole solution is always locally thermodynamically stable, while
that for the smaller solution is always unstable. If the limit is
taken in which the cavity radius becomes large while the horizon radii
of the black hole solutions remain fixed~\cite{BCM}, then $b T_B
\rightarrow 0$ and
\begin{equation}
b^{-2} C_A \rightarrow \frac{2 \pi}{3} \: \rho_+^2 \: (\rho_+^2+1)
 / (\rho_+^2-1),
\label{C_A_limit}
\end{equation}
clearly changing sign at $\rho_+ = 1$.

The reduced action approach to black hole thermodynamics may easily be
generalised to the $\Lambda \neq 0$ case. Choosing, for example, a
definition in which both the action and energy of hot anti-de Sitter
space vanish, the energy function becomes
\begin{equation}
b^{-1} E = \frac{1}{3} \: \rho_B^{1/2} \: \left\{ \bigl( \rho_B^3 + 3
\rho_B^{\phantom{0}} \bigr)^{1/2} - \bigl( \rho_B^3 + 3 \rho_B^{\phantom{0}} -
\rho_+^3 - 3 \rho_+^{\phantom{0}} \bigr)^{1/2} \right\}.
\end{equation}
Variation of the corresponding reduced action with respect to $\rho_+$
yields stationary points if the relation (\ref{sch_ads_surface}) is
satisfied, and the second variation at these points again demonstrates
local dynamic stability of the classical solution if and only if $C_A
> 0$. The free energy is negative only between $\rho_+ = \rho_B$ and
$\rho_+ = F_0(\rho_B)$, where $F_0$ is defined implicitly as the
solution curve for $\rho_+$ and $\rho_B$ along which the free energy
vanishes. Within this range the partition function will be dominated
by the classical solution, and so the equilibrium configuration will
be globally thermodynamically stable. Of course, Whiting's gauge
invariant proof of the correspondence between local dynamic and
thermodynamic stability is immediately applicable since it relies only
on the general form of $I_{\star}$, and not on the specific forms of
$E, \beta_B$ and $\overline{\beta}(\rho_+)$ individually.

As the limit (\ref{C_A_limit}) makes clear, an important difference
between this case and the previous one lies in the asymptotic
behaviour of the functions $C_0$ and $F_0$. As the cavity size
increases, these tend to 1 and $\sqrt{3}$ respectively, in contrast to
their analogues in asymptotically flat space which both become
infinite. It is clear, therefore, that for an arbitrarily large cavity
radius  a sufficiently low wall temperature may always be chosen such
that both black hole solutions are of finite size. When the
cosmological constant is zero this is no longer the case since the
turning point of the solution curve, which lies along the line $r_+ =
2 r_B / 3$, occurs at an infinite radius.

In the asymptotically flat case, then, the large black hole solution
is ill-defined as $r_B \rightarrow \infty$, and only the smaller
solution remains. The flat space canonical ensemble simply cannot
survive this limiting process. In the present case, however, with a
suitable rescaling of the wall temperature, a well-defined limit may
be taken in which the elements of the canonical ensemble are both
isothermal and infinite in extent, ostensibly because the functions
$C_0$ and $F_0$ remain finite. Furthermore, all of the physical and
thermodynamic characteristics of the solutions in finite cavities are
retained. This latter observation is essential in the forthcoming
analysis.

Since $\partial / \partial \tau$ possesses no natural normalisation at
infinity, a rescaling of the locally measured temperature merely
amounts to a rescaling of this Killing vector. The temperature is,
after all, only a Lagrange multiplier. With this in mind, if $T_B
\rightarrow \Delta \: T_B$ where the factor $\Delta$ is defined as
\begin{equation}
\Delta \stackrel{\mr{def}}{=} \lim_{r \rightarrow \infty} \: V(r)^{1/2}
\end{equation}
then, in the $\rho_B \rightarrow \infty$ limit, the fixed temperature
at infinity may be identified with $1 / \beta_{\star}$. The solution
recovered through this process is formally identical to that which is
analysed by Hawking and Page~\cite{PAGE}.

\begin{figure}
\begin{picture}(0,0)(0,0)
\put(290,0){\small $\rho_+$}
\put(2,240){\small $b^{-1}\beta_{\star}$}
\put(5,211){\small $\beta_{\mathrm{max}}$}
\put(68,5){\small 1}
\put(95,5){\small $\sqrt{3}$}
\put(-32,135){\small $C_A < 0$}
\put(160,130){\small $C_A > 0$}
\put(55,243){\small $F > 0$}
\put(100,243){\small $F < 0$}
\end{picture}
\centering\epsfig{file=sch_ads2.eps,width=10cm}
\caption{The black hole solution curve in the $(\rho_+,\beta_{\star})$
plane, in the limit of infinite isothermal cavity radius. The turning
point $\beta_{\mathrm{max}} = 1/T_0$ occurs at $\rho_+ = 1 \Leftrightarrow
r_+ = b / \sqrt{3}.$}
\label{sch_ads_beta}
\end{figure}

In the transformed coordinate $\rho$, the inverse temperature takes
the form
\begin{equation}
b^{-1} \: \beta_{\star} = \frac{4 \pi}{\sqrt{3}} \: \rho_+^{\phantom{0}} /
(\rho_+^2 + 1).
\end{equation}
The features retained from the finite cavity case are immediately
apparent: in contrast to the infinite asymptotically flat case, a
minimum temperature still remains below which there are no black hole
solutions. This temperature occurs in the case of the unique solution
with horizon size $\rho_+ = 1$, where $b^{-1} \: \beta_{\mathrm{max}} = 2
\pi / \sqrt{3}$. For temperatures greater than $T_0 = 1 /
\beta_{\mathrm{max}}$ two distinct black hole solutions exist, with
horizon size larger and smaller than 1 respectively. As is manifest
from the limit (\ref{C_A_limit}), the heat capacity for the solutions
with $\rho_+ > 1$ is positive, while for those with $\rho_+ < 1$ it is
negative.

Again this may be related to the local dynamic stability using the
reduced action which, since $E=M$, takes the form
\begin{equation}
I_{\star} = \frac{b}{2 \sqrt{3}} \: \beta^{\phantom{0}}_{\infty} \:
\rho_+^{\phantom{0}} \: (1 + \rho_+^2/3) - \frac{b^2}{3} \pi \rho_+^2.
\end{equation}
Variation with respect to the horizon scale then yields
\begin{equation}
\frac{\partial I_{\star}}{\partial \rho_+} = 0 \: \: \Leftrightarrow
\: \: \beta_{\infty} = \beta_{\star}
\end{equation}
and
\begin{equation}
b^{-2} \: \left. \frac{\partial^2 I_{\star}}{\partial \rho_+^2} \:
\right|_{\beta_{\infty} =  \beta_{\star}} = \mbox{\hspace{5mm}}
\frac{2 \pi}{3} \: (\rho_+^2 - 1)/(\rho_+^2 + 1).
\end{equation}
The sign of the second variation is then clearly the same as the sign
of the heat capacity for the equilibrium configuration. Furthermore,
by forcing the condition $\beta_{\infty} = \beta_{\star}$ in
$I_{\star}$, the action for the classical solutions becomes
\begin{equation}
b^{-2} \: I = - \frac{\pi}{9} \: \rho_+^2 \: (\rho_+^2 - 3) /
(\rho_+^2 + 1),
\end{equation}
from which it may be seen that the free energy indeed changes sign at
$\rho_+ = \sqrt{3}$, in agreement with the limiting value of
$F_0$. Clearly, both the action and the free energy are negative if
the horizon size is greater than $\sqrt{3}$.

These features of the infinite Schwarzschild-anti-de Sitter solution
are summarised in figure~\ref{sch_ads_beta}, and their significance
for the canonical ensemble is exhaustively explained by Hawking and
Page~\cite{PAGE}.

\sect{Nonconformal Negative Modes}
It is evident from the previous section, then, that the correspondence
between local thermodynamic stability of an equilibrium configuration
in the canonical ensemble and local dynamic stability of the classical
solutions is well established. This correspondence may be shown in a
gauge invariant way through variations of the reduced action for the
ensemble with respect to a generalised horizon radius, and is true for
both asymptotically flat and asymptotically anti-de Sitter space
confined within a finite ideal isothermal cavity. In the latter case,
however, both the canonical ensemble and the above correspondence
remain well defined, if suitably stated, even in the limit of infinite
cavity radius.

Nonetheless, it is desirable to have a more direct and intuitive
indication of the dynamic stability of the classical solutions in
terms of a perturbation in the background geometry rather than in a
variation of the reduced action. The measure of stability in this
approach is the eigenvalue spectrum of the nonconformal perturbative
modes for the solution. Should there be a mode with a negative
eigenvalue, then the action for this solution is a saddle-point in its
phase space rather than a true minimum. Consequently, there ought to
be a correspondence between the presence of such a negative mode and
the local thermodynamic stability as governed by the heat capacity. As
will become clear in the forthcoming sections, it is the existence of
a formal limit in the anti-de Sitter case that permits a well defined
prescription for this problem, which seems impossible in the
asymptotically flat case.

Negative modes arise from the analysis of geometric fluctuations about
classical Euclidean solutions of the Einstein field equations.
However, the analysis to confirm their existence must be performed
with care, since the gauge freedom of the Euclidean action will in
general introduce a large number of non-physical negative modes
associated with conformal deformations of the metric. For pure
gravity, the contributions from the conformal and the nonconformal
modes decouple if a suitable gauge is chosen. This procedure is now
well understood.

In the path integral approach~\cite{SWH2}, the partition function $Z$
is generally defined as a functional integral over all metrics with
some fixed asymptotic behaviour on some $n$-dimensional manifold
$\mathcal{M}$,
\begin{equation}
Z = \int_{\mathcal{M}} \! \! D[g] \: \mathrm{e}^{-i I[g]}.
\end{equation}
This integral is formally defined by an analytic continuation to a
Euclidean section of $\mathcal{M}$ to become
\begin{equation}
Z = \int_{\mathcal{M}} \! \! D[g] \: \mathrm{e}^{-\hat{I}[g]},
\label{partition}
\end{equation}
where the integral is performed over all positive definite metrics
$g$. In the case of pure gravity, the Euclidean action $\hat{I}$ is
\begin{equation}
\hat{I} = -k \! \int_{\mathcal{M}} \! \! d^n x \: \sqrt{g} \: \left\{ R - 2
\Lambda \right\} - 2k \! \int_{\partial\mathcal{M}} \! \! d\Sigma \: \:
\mathrm{Tr} K
\end{equation}
where $k$ is a coupling constant and $K$ is the second fundamental
form on the boundary $\partial\mathcal{M}$.

This partition function may be approximated using saddle-point
techniques, by Taylor expanding about the known stationary points
of the Euclidean action -- the solutions to the Einstein field
equations
\begin{equation}
R_{ab} = \frac{2}{n-2} \Lambda \: g_{ab}.
\end{equation}
The expansions are performed by writing the perturbed metric
$\tilde{g}$ as
\begin{equation}
\tilde{g}_{ab} = g_{ab} + \phi_{ab}
\end{equation}
with $\phi$ treated as a quantum field on the classical fixed
background $g$ which vanishes on the boundary $\partial
\mathcal{M}$. Proceeding in this fashion yields
\begin{equation}
\hat{I}[\tilde{g}] = \hat{I}[g] + \hat{I}_2[\phi] + \cdots
\end{equation}
where the linear term $\hat{I}_1$ vanishes precisely because $g$
is a classical solution, $\hat{I}_2$ is quadratic in the field $\phi$,
and `$\cdots$' represents terms of higher than quadratic order.
Ignoring these additional terms, and inserting the remainder into
(\ref{partition}), gives the well-known expansion
\begin{equation}
\log Z = -\hat{I}[g] + \log \int_{\mathcal{M}} \! \! D[\phi] \: \mathrm{e}^{
-\hat{I}_2 [\phi]}.
\end{equation}
The second term on the right is generally called the one-loop
contribution to $\log Z$.

The quadratic contribution to the action is straightforward to
evaluate, and may be written for arbitrary $\phi$ in the
form
\begin{equation}
\hat{I}_2[\phi] = \frac{k}{2} \int \! \! d^n x \: \sqrt{g} \:
\phi^{ab} A_{abcd} \: \phi^{cd}.
\end{equation}
In this context, the operator $A$ takes the rather opaque
form~\cite{SWH2}
\begin{eqnarray}\samepage
\lefteqn{ A_{abcd} = \frac{1}{4} g_{cd} \nabla_a \nabla_b -
\frac{1}{4} g_{ac} \nabla_d \nabla_b + \frac{1}{8} \{ g_{ac} g_{bd} -
g_{ab} g_{cd} \} \nabla_e \nabla^e + \frac{1}{2} R_{ad} g_{bc} } & &
\nonumber \\[2mm]
& & \mbox{\hspace{2cm}} -\frac{1}{4} R_{ab} g_{cd} + \frac{1}{16} R
g_{ab} g_{cd} - \frac{1}{8} R g_{ac} g_{bd} - \frac{1}{8} \Lambda
g_{ab} g_{cd} + \frac{1}{4} \Lambda g_{ac} g_{bd} \nonumber \\[3mm]
& & \mbox{\hspace{4.5cm}} + (a \leftrightarrow b) + (c \leftrightarrow
d) + (a \leftrightarrow b, c \leftrightarrow d).
\end{eqnarray}
Naively, then, the one-loop term may be written as $\frac{1}{2}
\log \mathrm{Det} (\mu^{-2} A)$ where $\mu$ is a regularisation mass, and
the determinant is formally defined as the product of the eigenvalues
of $A$. However, due to the diffeomorphism gauge freedom of the
action, $A$ will in general have a large number of zero eigenvalues,
and so this procedure as stated is ill-defined. The remedy is to add a
gauge fixing term $B$ -- such that the operator $A+B$ has no zero
eigenvalues -- and an associated ghost contribution $C$, to obtain
\begin{equation}
\log Z = -\hat{I}[g] - \frac{1}{2} \log \mathrm{Det} (\mu^{-2} \{A+B\}) +
\log \mathrm{Det} (\mu^{-2} C).
\end{equation}

Such terms may be dealt with by means of generalised zeta functions,
as considered by Gibbons, Hawking, and Perry~\cite{GHP}, and extended
to include a $\Lambda$ term by Hawking~\cite{SWH2}. In order that this
be possible, the terms must be expressed as sums of operators, each
with only a finite number of negative eigenvalues. This may be
achieved by writing $A+B$ as $-F+G$, where $F$ is a scalar operator
acting on the trace of $\phi$. The ghost term $C$ is a spin-1 operator
acting on divergence-free vectors, and $G$ is a physical gauge
invariant spin-2 operator which acts on the transverse and
trace-free part of $\phi$. For any dimension $n \ge 4$, $G$ takes the
simple form
\begin{equation}
G_{abcd} = -g_{ac}g_{bd} \nabla_e \nabla^e - 2 R_{acbd}.
\end{equation}

An observation of primary significance is that, for $\Lambda \le 0$, a
gauge may be chosen in which the operators $F$ and $C$ have no
negative eigenvalues. If the background metric $g$ is flat then, in
addition, $G$ will be positive-definite, but for a non-flat background
this is not the case. In general, $G$ will have some finite number --
generally zero or one -- of negative eigenvalues, which correspond to
the nonconformal negative modes of the solution. The eigenvalues of
$G$ are determined by all solutions to the elliptic equation
\begin{equation}
G^{ab}_{\phantom{ab}cd} \: \phi_{(n)}^{\phantom{n}cd} = \lambda^{\phantom{a}}_{(n)} \:
\phi_{(n)}^{\phantom{n}ab}
\label{G_eq}
\end{equation}
where the eigenfunctions $\phi$ are real, regular, symmetric,
transverse, trace-free, and normalisable tensors. Clearly, should one
of the eigenvalues of $G$ be negative, then the product of all of the
eigenvalues would also be negative. The contribution to $\log Z$ from
fluctuations about the classical solution would then contain an
imaginary component, leading to an instability in the ensemble similar
to the type proposed by Gross, Perry and Yaffe~\cite{GPY}.

\sect{General Spherically Symmetric Solutions}
Following the treatment given in~\cite{GPY}, and subsequently
in~\cite{ALLEN}, for the Schwarzschild case, it is clear that only
spherically symmetric and $\tau$-independent solutions of (\ref{G_eq})
need be considered as candidate nonconformal negative modes. In static
and spherically symmetric backgrounds, modes of higher multipole
moment will necessarily have greater eigenvalues. With this
assumption, it is then straightforward to write down a construction
for such solutions to $G$ valid in {\sl any} $n$-dimensional Euclidean
black hole background of the form (\ref{metric}) for a general metric
function $V(r)$. In the following expressions, a prime denotes $d/dr$.

Since $G$ acts only on symmetric transverse and trace-free tensors,
then clearly the constructed solutions must exhibit all of these
properties. If the mode $\phi_{ab}$ is written in the manifestly
trace-free and symmetric form
\begin{equation}
\phi^a_{\phantom{a}b} = \mathrm{diag} \: \Big( \psi(r), \chi(r),
\underbrace{k(r), \ldots, k(r)}_{(n-2) \: \mathrm{terms}} \Big)
\label{tracefree}
\end{equation}
where the function
\begin{equation}
k(r) = \frac{1}{2-n} \{ \psi(r) + \chi(r) \},
\end{equation}
then the final property, $\nabla^a \phi_{ab}$ = 0, is guaranteed if it
is further assumed that $\psi(r)$ and $\chi(r)$ are related through
the first order equation
\begin{equation}
\psi(r) = \frac{2rV}{rV'-2V} \: \chi'(r) +  \frac{rV'+2(n-1)V}{rV'-2V}
\: \chi(r).
\label{transverse}
\end{equation}

With the ansatz (\ref{tracefree}) and (\ref{transverse}), the
eigenvalue equation (\ref{G_eq}) reduces to a linear second order
ordinary differential equation for the component $\chi(r)$ which, for
general $V(r)$, is
\begin{eqnarray}\samepage
\lefteqn{ -V \chi''(r) +  \left[ \frac{2r^2 [ VV''- V'^2 ] -
r(n-2)VV'+ 2nV^2}{r(rV'-2V)} \right] \chi'(r) } & & \nonumber \\[3mm]
 & & \mbox{\hspace{15mm}} + \left[ \frac{r^2V'V'' + r [ 2(n-1)VV'' -
(n+2)V'^2 ] +4VV'}{r(rV'-2V)} \right] \chi(r) = \lambda \chi(r).
\label{diff_eq}
\end{eqnarray}
The problem of finding the eigenfunctions of $G$ is therefore reduced,
in the cases of interest, to the problem of finding real and
normalisable solutions to this equation that satisfy appropriate
boundary conditions.

For the backgrounds of interest, (\ref{diff_eq}) may have as many as
five singular points, only three of which are relevant. These are at
the horizon radius $r = r_+$, at $r = r_s$ -- the real solution to
$rV' - 2V = 0$, and at $r = \infty$. It is clear from
(\ref{transverse}) that at the first of these, where $V(r)$ vanishes,
the components $\psi(r)$ and $\chi(r)$ are equal and so the metric
perturbation does not alter the periodicity of $\tau$. Furthermore, in
order that $\phi_{ab}$ be  regular everywhere, the product $(rV'-2V)
\: \psi(r)$ must vanish at the second singular point, and hence
\begin{equation}
\chi'(r_s) / \chi(r_s) \: = \left. - \frac{1}{2rV} \: \biggl[rV' +
2(n-1)V \biggr] \: \right|_{r = r_s}.
\label{f_deriv}
\end{equation}

The normalisation of the eigenfunctions may be defined as
\begin{equation}
\mathcal{N}^2 = \int_{\mathcal{M}} \!\! d^nx \: \sqrt{g} \: \phi^{ab} \,
\phi_{ab},
\label{big_norm}
\end{equation}
so that acceptable solutions must be square-integrable functions in
the sense that $\mathcal{N}^2$ be finite, although the extent of $\mathcal{M}$
will depend on the situation under consideration. (Henceforth, any
reference to the solution will imply the normalised function
$\mathcal{N}^{-1} \phi_{ab}$, unless $\mathcal{N}$ is explicitly written.) In
the infinite cavity limit this requirement, and the condition of
regularity at the interior singular points, provide the necessary
boundary conditions to solve for viable solutions. For black holes
confined within finite cavities however, this is no longer the case,
since all regular solutions to (\ref{diff_eq}) will be normalisable.
In such cases the canonical boundary conditions for the ensemble must
be imposed at the cavity wall, and it is this requirement which
presents difficulties.

The perturbed metric $g_{ab} + \epsilon \phi_{ab}$ induces a line
element
\begin{equation}
ds^2  = V(r) \biggl[ 1 + \epsilon \psi(r) \biggr] \: d \tau^2 +
\frac{1}{V(r)} \biggl[ 1 + \epsilon \chi(r) \biggr] \: dr^2 + r^2
\biggl[1 + \epsilon k(r)\biggr] \: d \Omega^2,
\end{equation}
and so the proper length around the $\mathrm{S}^1$ in the $\tau$ direction
at radius $r_B$, which may be identified with $1 / T_B$, becomes
\begin{equation}
1 / T_B = \beta_{\star} \: V(r_B)^{1/2} \: \bigg[1 + \epsilon
\psi(r_B) \biggr]^{1/2}.
\end{equation}
For an isothermal cavity wall at $r = r_B$ the temperature must
remain constant and so the boundary condition, as used by
Allen~\cite{ALLEN}, is
\begin{equation}
\psi(r) \: |_{r = r_B} = 0.
\label{bc1}
\end{equation}
However, this condition is insufficient, since the elements of the
canonical ensemble are defined by both an invariant wall temperature
and an invariant wall area $A_B$. In the perturbed metric the surface
area of a space-like spherical shell at radius $r_B$ is no longer $4
\pi r_B^2$, but becomes
\begin{equation}
\tilde{A}_B = 4 \pi r_B^2 \: \biggl[1 + \epsilon k(r_B) \biggr].
\label{new_area}
\end{equation}
As York points out~\cite{YORK1}, evaluating the eigenvalue spectrum
with the boundary condition (\ref{bc1}) is not physically meaningful
for the classical solutions since the wall area at constant radius is
not invariant. This explains the discrepancy between the cavity radius
$r_B = 3 r_+ / 2$ at which the heat capacity of a Schwarzschild black
hole changes sign, and the critical radius at which the negative mode
vanishes in this case, calculated by Allen to be $r_{\mathrm{crit}} \sim
2.89 \, r_+ / 2$.

So, to apply the correct boundary condition, a modified cavity radius
$r_{B'}$ must be chosen such that $\tilde{A}_{B'} = A_B$, and then the
condition
\begin{equation}
1 + \epsilon \psi(r_{B'}) = V(r_B) / V(r_{B'})
\end{equation}
must be imposed. However, it is clear from (\ref{new_area}) that,
other than in the degenerate case $k(r) = \psi(r) = \chi(r) \equiv 0$,
this procedure is not particularly well defined since the parameter
$\epsilon$ is freely variable. The areas $A_B$ and $\tilde{A}_{B'}$
may therefore be matched at any nearby radius simply by adjusting
$\epsilon$ appropriately. Then for a given choice of $\epsilon$ the
modified cavity radius will be uniquely defined, and the eigenvalue
spectrum of (\ref{diff_eq}) may be evaluated, but this spectrum will
in general be slightly different for each value of $\epsilon$.

This problem disappears when the cavity is infinite. If the solution
$\phi_{ab}$ is normalisable in the sense discussed above, then the
radial functions $\psi(r)$, $\chi(r)$ and $k(r)$ must tend to zero at
some appropriate rate governed by the asymptotic form of
(\ref{diff_eq}). Consequently, the isothermal boundary condition at
infinity is automatically satisfied.

The original calculation of the negative mode eigenvalue for
Schwarzschild in $n = 4$ dimensions, which is manifestly invariant of
the magnitude of $\epsilon$, finds $M^2 \: \lambda_{\mathrm{neg}} \sim
-0.192$. In the case of Schwarzschild-anti-de Sitter, a similar
calculation will achieve two goals. It will confirm the conjecture
made by Hawking and Page regarding the existence of a nonconformal
negative mode and its vanishing point at the critical radius $\rho_+ =
1$. It will also provide a well defined demonstration of the
correspondence between the local dynamic stability of a solution in
terms of its eigenvalue spectrum and the local thermodynamic stability
of the equilibrium configuration. This is only possible since, in
contrast to the flat case, the infinite cavity solution in the
asymptotically anti-de Sitter case retains all of the physical and
thermodynamic characteristics of the finite cavity solutions.

\sect{The Eigenvalue Spectrum in Four Dimensions}
For Schwarzschild-anti-de Sitter in $n=4$ dimensions, the metric
function is
\begin{equation}
V(r) = \frac{1}{r} \: \biggl( r - r_+^{\phantom{0}} - r_+^3/b^2_{\phantom{0}} +
r^3_{\phantom{0}}/b^2_{\phantom{0}} \biggr),
\end{equation}
in which the horizon radius $r_+$ has been substituted for the mass
function $M$ from (\ref{mass_func}). Moving to the dimensionless radial
coordinate $\rho$ defined in (\ref{rho_def}), and defining a new
parameter $\alpha$ such that
\begin{equation}
\alpha \stackrel{\mr{def}}{=} \frac{1}{4} \: \bigl( \rho_+^3 + 3 \rho_+^{\phantom{0}}
\bigr),
\end{equation}
the second order equation for $\chi(\rho)$ may be written in the
homogeneous form
\begin{eqnarray}\samepage
\lefteqn{ \biggl[ \rho \, ( \rho-2\alpha) \, (\rho^3+3\rho-4\alpha)
\biggr] \: \chi''(\rho) + 4 \biggl[ 2 \rho^4 - 5 \alpha \rho^3 + 3
\rho^2 - 11 \alpha \rho + 8 \alpha^2 \biggr] \: \chi'(\rho) } & &
\nonumber \\[2mm]
 & & \mbox{\hspace{5cm}} + \biggl[ (10 + \tilde{\lambda}) \rho^3 - 2
\alpha ( 18 + \tilde{\lambda}) \rho^2 - 16 \alpha \biggr] \:
\chi(\rho) = 0,
\label{4d_diff_eq}
\end{eqnarray}
where a prime now denotes $d/d\rho$, and $\tilde{\lambda} = b^2
\lambda$.

The three real singular points of this equation are at $\rho =
\rho_+$, $2 \alpha$, and $\infty$. The remaining two imaginary
singular points are irrelevant in this analysis. Around $\rho =
\rho_+$, a trial solution of the form $\chi(\rho) = \sum_n \: a_n \:
(\rho - \rho_+)^{k+n}$ yields the indicial roots $k = 0$ and $-1$. The
$k = 0$ case gives the unique regular series since the $k = 1$ series
is not normalisable in the sense that $\mathcal{N}^2$ is infinite. If $a_0$
is chosen to be 1, implying the boundary condition $\chi(\rho) \:
|_{\rho = \rho_+} = 1$, then the next two coefficients become
\begin{equation}
a_1 = - \frac{1}{6} \frac{ \rho_+^2 (18 + \tilde{\lambda}) +
24}{\rho_+^{\phantom{0}} (\rho_+^2 + 1)}
\end{equation}
and
\begin{equation}
a_2 = \frac{1}{108} \frac{ \rho_+^4 ( 648 + 54 \tilde{\lambda} +
\tilde{\lambda}^2_{\phantom{0}}) + \rho_+^2 (1620 + 66 \tilde{\lambda}) +
1008}{ \rho_+^2 ( \rho_+^2 + 1)^2_{\phantom{0}}}.
\end{equation}
A similar expansion with coefficients $b_n$ around $\rho = 2 \alpha$
gives $k = 0$ and 3, both of which produce regular solutions. However,
only the $k = 0$ case contributes terms of order less than
$(\rho-2\alpha)^3$, for which the first three terms are related through
\begin{equation}
b_1 = - \frac{2}{\alpha} \: b_0
\end{equation}
from (\ref{f_deriv}), and
\begin{equation}
b_2 = \frac{1}{2} \frac{ \alpha^2 (18 + \tilde{\lambda}) + 6}{\alpha^2
(4 \alpha^2+1)} \: b_0.
\end{equation}
Using these coefficients a numerical solution to (\ref{4d_diff_eq})
may be generated which is constrained to within $\mathcal{O}(\epsilon^3)$
for $\epsilon \ll 1$ at the interior singular points. Any simple
algorithm is limited to this accuracy by the existence of the $k=3$
series about $\rho = 2 \alpha$.

For a given choice of horizon size $\rho_+$, the aim of such a
numerical method is to find the spectrum of eigenvalues
$\tilde{\lambda}_{(n)}$, and in particular the minimum eigenvalue
$\tilde{\lambda}_{\mathrm{neg}}$, for which regular and normalisable
solutions $\chi(\rho)$ exist. Having fixed $\rho_+$, then, a value for
$\tilde{\lambda}$ is guessed, and an appropriate numerical integrating
routine~\cite{NUM_REC} is used to integrate away from the horizon
using the coefficients $a_n$ as initial data. The solution is
calculated to within a small step $\delta$ of $2 \alpha$, and
then analytically continued through the singular point to $2 \alpha +
\delta$ using the $b_n$. From here, the solution is then calculated
up to some large radius, at which point the criterion of
normalisability is assessed. For most choices of $\tilde{\lambda}$
the solution will fail, but for some choices -- specifically the
eigenvalues -- it will not.

The conditions upon which normalisability is assessed may be obtained
from the specific form in this case of $\mathcal{N}^2$ and the asymptotic
behaviour of (\ref{4d_diff_eq}). With the benefit of relations
(\ref{tracefree}) and (\ref{transverse}), the relation
(\ref{big_norm}) may be translated into an integral condition in one
dimension on $\chi(\rho)$, so that
\begin{equation}
\mathcal{N}^2 = 4 \pi \beta^{\phantom{0}}_{\star} \: \int_{\rho_+}^{\infty} \! \!
d \rho \: \Theta(\rho) \: \chi^2(\rho).
\end{equation}
The function $\Theta(\rho)$ is well defined and finite at the interior
singular points, and has asymptotic behaviour $\sim \rho^6$ as $\rho
\rightarrow \infty$. This behaviour is, of course, the same as that
for the weight function derived from the self-adjoint form of
(\ref{4d_diff_eq}). The solution $\phi_{ab}$ will therefore be
normalisable if $\chi(\rho) \sim \rho^{-(7/2+x)}$ for any
positive $x$.

\begin{figure}
\centering\rotatebox{270}{\epsfig{file=4d.ps,width=11cm}}
\caption{Numerically generated results for
$\tilde{\lambda}_{\mathrm{neg}}$ (vertical) against $\rho_+$ in $n=4$
dimensions Schwarzschild-anti-de Sitter.}
\label{4d_ads}
\end{figure}

The asymptotic form of (\ref{4d_diff_eq}) is readily obtained with the
standard substitution $\rho \rightarrow 1/u$. In the limit $u
\rightarrow 0$ this equation becomes
\begin{equation}
u^2 \: \ddot{\chi}(u) - 6u \: \dot{\chi}(u) + [10 + \tilde{\lambda}]
\: \chi(u) = 0
\label{asymptotic_eq}
\end{equation}
where an over-dot denotes $d/du$ -- the solutions to which describe the
asymptotic behaviour of $\chi(\rho)$. It should be noted that, in
contrast with the asymptotically flat case, the singular point at
$\rho = \infty$ is now regular. In this limit,
\begin{equation}
\chi(\rho) \sim A(\tilde{\lambda}) \: \rho^{-7/2 + \sqrt{9/4 -
\tilde{\lambda}}} \: + \: B(\tilde{\lambda}) \: \rho^{-7/2 - \sqrt{9/4
- \tilde{\lambda}}}.
\label{asymptotic_sol}
\end{equation}

Two points are immediately apparent from the form of this solution.
Firstly, in accordance with the naive expectation that the
eigenfunction corresponding to the lowest eigenvalue is the
`smoothest', the value of $\tilde{\lambda}_{\mathrm{neg}}$ ought to be
strictly less than 9/4. For eigenvalues greater than this the exponent
term acquires an imaginary component, and the solutions oscillate in
$\log \! \rho$. Secondly, the criterion of normalisability is clearly
equivalent to the single requirement $A(\tilde{\lambda}) \equiv 0$.

A simple numerical algorithm for finding $\tilde{\lambda}_{\mathrm{neg}}$
is therefore to proceed as described above and, for some fixed
$\rho_+$ and a guessed $\tilde{\lambda}$, integrate $\chi(\rho)$ out
to the second interior singular point at $2 \alpha$. Then, transform
the solution to a new function $\xi(\rho) = \rho^{7/2} \: \chi(\rho)$
and integrate this out to some large radius. As $\rho$ increases the
new function approaches
\begin{equation}
\xi(\rho) \sim A(\tilde{\lambda}) \: \rho^{\sqrt{9/4 -
\tilde{\lambda}}}
\end{equation}
since the second term must be very small. For some choice of
$\tilde{\lambda}$, then, this remaining term will be negative while
for another it will be positive -- the two choices bracketing the true
value at which $A(\tilde{\lambda})$ vanishes. By subsequently refining
the range over which $\tilde{\lambda}$ is guessed the minimum
eigenvalue can be found to reasonable accuracy. This algorithm is
inelegant but nonetheless quite effective. Results for $\rho_+$ in
the range $0.1 \rightarrow 4.0$ are drawn in figure~\ref{4d_ads}.

The graph clearly shows a single negative mode for $n=4$ dimensions
Schwarzschild-anti-de Sitter vanishing at $\rho_+ = 1$. This confirms
the conjecture of Hawking and Page~\cite{PAGE}, and demonstrates in a
well defined manner the correspondence between this eigenvalue
spectrum and the local thermodynamic stability of the equilibrium
configuration in the canonical ensemble. As an aside, it should be
noted that
\begin{equation}
\lim_{\rho_+ \rightarrow 0} \: M^2 \: \lambda_{\mathrm{neg}} \sim
-0.192
\end{equation}
which demonstrates that the results obtained previously for
Schwarzschild reappear quite naturally in the limit where $r_+ \ll b$.

\sect{Extending the Spectrum to Higher Dimensions}
The various quantities employed in the evaluation of the eigenvalue
spectrum in $n=4$ dimensions Schwarzschild-anti-de Sitter are readily
generalised to any higher dimension. The form of the master equation
(\ref{diff_eq}) then makes it simple to verify that a nonconformal
negative mode will exist for a classical solution in any dimension $n
\ge 4$ only when the heat capacity of the equilibrium configuration is
negative.

In $n$ dimensions, the metric function $V(r)$ for Schwarzschild-anti-de
Sitter becomes
\begin{equation}
V(r) = 1 - 2M / r^{n-3} + r^2/b^2
\end{equation}
where $b^2 = -(n-1)(n-2)/2 \Lambda$ and so, if the horizon is at $r_+$
such that $V(r_+) = 0$, the mass may be expressed as
\begin{equation}
M = \frac{1}{2} \: r_+^{n-3} \: (1 + r_+^2 / b^2_{\phantom{0}}).
\end{equation}
The periodicity then generalises to
\begin{equation}
\beta_{\star} = \frac{4 \pi b^2 r_+}{(n-3)b^2_{\phantom{0}} + (n-1)r_+^2}
\end{equation}
and so the heat capacity becomes
\begin{equation}
C_A = 2 \pi r_+^{n-2} \: \left[ \frac{ (n-1)r_+^2 + (n-3)b^2_{\phantom{0}}
}{ (n-1)r_+^2 - (n-3)b^2_{\phantom{0}} } \right].
\end{equation}
The form of the heat capacity suggests the definition of a
dimensionless radial coordinate $\rho$ such that
\begin{equation}
\rho \stackrel{\mr{def}}{=} (n-1)^{1/2} \: (n-3)^{-1/2} \: r/b
\end{equation}
which clearly reduces to (\ref{rho_def}) when $n=4$. In any dimension
then, there is an infinite discontinuity in $C_A$ at $\rho_+ = 1$,
with $C_A$ negative for $\rho_+ < 1$ and positive for $\rho_+ > 1$.

Inserting into (\ref{diff_eq}) any $n$, and with it the appropriate
form of $V(r)$ obtained from above, yields a second order equation for
$\chi(\rho)$ analogous to (\ref{4d_diff_eq}). Two of the three
relevant singular points of this equation remain fixed -- at $\rho_+$
and $\infty$, while the third is dimension-specific. All three remain
regular, however, regardless of the dimension.

The weight function obtained from the self-adjoint form has
asymptotic behaviour $\sim \rho^{n+2}$, while the general solution to
the asymptotic form, obtained with the substitution $\rho \rightarrow
1/u$, is
\begin{equation}
\chi(\rho) \sim A(\tilde{\lambda}) \: \rho^{-(n+3)/2 + \beta} +
B(\tilde{\lambda}) \: \rho^{-(n+3)/2 - \beta}
\end{equation}
where
\begin{equation}
\beta^2 + \tilde{\lambda} = (n-1)^2 / 4.
\end{equation}
The criterion of normalisability clearly remains the same --
i.e. $A(\tilde{\lambda}) \equiv 0$, and $\tilde{\lambda}_{\mathrm{neg}} <
(n-1)^2 / 4$.

The numerical results for $\tilde{\lambda}_{\mathrm{neg}}$ obtained for
the cases $n=5$ and 7 for $\rho_+$ in the range $0.1 \rightarrow 4.0$
are shown together in figure~\ref{5d_and_7d_ads}. The conjecture of
Hawking and Page, and the correspondence between this eigenvalue
spectrum and local thermodynamic stability, may clearly be extended
to dimensions greater than four.

\begin{figure}
\centering\rotatebox{270}{\epsfig{file=5d_and_7d.ps,width=11cm}}
\caption{Numerically generated results for
$\tilde{\lambda}_{\mathrm{neg}}$ (vertical) against $\rho_+$ in $n=5$ and
$n=7$ dimensions Schwarzschild-anti-de Sitter. The spectrum for $n=7$
gives the `upper' curve (tending to 5), and for $n=5$ the `lower'
curve (tending to 3).}
\label{5d_and_7d_ads}
\end{figure}

\sect{Discussion and Further Applications}
The discovery of a unique nonconformal negative mode for the $n=4$
dimensions Euclidean Schwarzschild solution led Gross, Perry and Yaffe
to propose the instability of flat space at finite temperature due to
the barrier-tunnelling nucleation of black holes. The tunnelling process
is permitted because, in the semi-classical approach, an integration
over quadratic fluctuations of the Schwarzschild solution produces an
imaginary component in the free energy of the system. The existence of
a negative mode is therefore a powerful indication of instability.

The subsequent observation that the Euclidean Schwarzschild solution
becomes double-valued when confined within a box, and that hot flat
space may therefore be stabilised against black hole nucleation,
enabled York to give a thermodynamically consistent interpretation of
the tunnelling process within a finite isothermal cavity. Manifest
within this approach is a correspondence between local dynamic
stability of the classical black hole solutions and local
thermodynamic stability of the equilibrium configurations in the
canonical ensemble. This becomes a global thermodynamic stability if
the action, and so the free energy, is negative. It should also be the
case, therefore, that this correspondence applies to negative modes: a
nonconformal negative mode should exist only in those classical
solutions for which the heat capacity of the equilibrium configuration
is negative. However, there appears to be no well defined way of
demonstrating this explicitly in finite isothermal cavities.

This problem may be overcome by instead considering hot anti-de Sitter
space. Confined again to a finite isothermal cavity, the Euclidean
Schwarzschild anti-de Sitter solution is also double-valued, but in
this case it remains so, if suitably stated, even in the infinite
cavity limit. The result so obtained is formally identical to that
analysed by Hawking and Page, and retains all of the physical and
thermodynamic characteristics of the finite cavity solutions. This
observation not only implies that the canonical ensemble in anti-de
Sitter is well defined, but also presents an ideal opportunity in
which to demonstrate the expected correspondence between negative
modes and heat capacity.

Restricting the master equation (\ref{diff_eq}) derived above to the
case of $n=4$ dimensions Schwarzschild-anti-de Sitter, and solving for
the eigenvalue $\tilde{\lambda}_{\mathrm{neg}}$ of the smoothest
normalisable eigenfunction, results in an eigenvalue spectrum for a
range of classical solutions which indeed has the same sign as the
heat capacity function. This demonstrates the stated correspondence in
a manner invariant of the magnitude of the metric perturbation, and
confirms the conjecture of Hawking and Page that the negative mode
`passes through zero' at the solution with critical horizon radius
$r_+ = b/ \sqrt{3}$. For solutions with horizon radius greater than
this, no nonconformal negative mode exists. Furthermore, the
correspondence may be extended to any spacetime dimension $n >
4$. This latter observation is of particular relevance to more recent
ideas involving anti-de Sitter spaces.

Analogous to that found by Gross {\sl et al.}, a single nonconformal
negative mode exists in the $n=5$ dimensions Schwarzschild solution,
which is described by Witten~\cite{WITT1} as leading to semi-classical
instability of the ground state $\mathrm{M}^4 \times \mathrm{S}^1$ in the
original Kaluza-Klein theory without fermions. In this case the
periodically identified coordinate is not the imaginary time but the
compact fifth spatial direction. In fact the tunnelling instability
phenomenon, and the associated negative mode, seems to be a quite
general characteristic of any metric with periodic identifications --
called a `tachyonic solution'~\cite{GIBB1}. The Lorentzian
continuation of the metric may be interpreted as the result of the
tunnelling instability, or alternatively as providing the Cauchy
development of the time-symmetric initial data set on some plane of
symmetry regarded as $t=0$.

With this in mind, the continuous one-parameter family of initial data
discussed by Horowitz and Myers~\cite{MYERS}, represented by the line
element
\begin{equation}
ds^2 = \biggl[ 1 - M/r^2 + r^2/b^2 \biggr] \: d \tau^2 + \biggl[ 1 - M
/ r^2 + r^2/b^2 \biggr] ^{-1} \: dr^2 + r^2 \: d \Omega^2,
\end{equation}
may be considered as the restriction to the equatorial plane of the
three-spheres in the $n=5$ dimensional Schwarzschild-anti-de Sitter
metric periodically identified in $\tau$. This metric satisfies the
constraint equation $\vphantom{|}^4 R = -12 / b^2$, and the hyper-plane
on which it is defined may be regarded as $t=0$. The radial coordinate
is of course restricted to $r > r_+$ where this is the real positive
root of the metric function as usual.

The Euclidean section of the full $n=5$ dimensions metric is obtained
by setting $t = iz$, and is now known to possess a single nonconformal
negative mode if $r_+ < b / \sqrt{2}$. This mode will `pass through
zero' at equality, and hence solutions with $r_+ > b / \sqrt{2}$ will
have no negative modes. Furthermore, if $r_+$ is great enough, the
action and hence the energy will be negative. By analogy with the
earlier argument for global thermodynamic stability, this strongly
suggests that, should the action be negative, then indeed it will be a
global minimum for metrics with these boundary conditions.

The idea, then, is that the larger Schwarzschild-anti-de Sitter black
hole solution for a given periodicity $\beta_{\star}$, should its
action be negative, is in fact the global minimum of the action for
periodically identified metrics satisfying the scalar curvature
constraint $R = 2 n \Lambda / (n-2)$. This is clearly very similar to
the proposal in~\cite{MYERS} and it is believed that the above work
provides strong evidence in favour of this conjecture, albeit without
recourse to the AdS/CFT correspondence.

It is a pleasure to acknowledge Stephen Hawking for suggesting this
avenue of research, and to thank Malcolm Perry and Gary Gibbons for
their help and advice during the preparation of this work.