format: text

text:
# Stringy negative-tension branes and the second law of thermodynamics

## 1 Introduction

Branes with negative tension have experienced a recent surge in popularity. In part, this stems from interest in string theoretic objects called orientifolds  as perturbative string calculations tell us that their tension is sometimes less than zero. In addition, negative tension branes have found interesting applications in certain brane-world scenarios (e.g., in  and related work). Indeed, to the authors’ knowledge all such scenarios that address the gauge hierarchy problem contain some form of negative tension object.

It is important to recall that negative tension branes carry negative energy both in the sense that their total energy is less than zero and in the sense that they violate the weak null energy condition . Mention of such objects naturally raises the question of whether they pose a danger to one’s physical theory. On the one hand, the bold physicist may point out the orientifold construction of string theories, relying on string theory’s track record of producing novel pieces of physics which turn out to be surprisingly robust. Past examples of this sort include certain kinds of topology changing processes  and dualities between gravitating theories and non-gravitating gauge theories .

On the other hand, a more conservative physicist may point out excellent reasons to be skeptical of negative energies111While it is understood that negative energy fluxes can arise in quantum field theory , there is significant evidence that such fluxes face fundamental constraints  restricting them to domains where quantum mechanical effects are important and in which their effects are more benign (see e.g. ). In contrast, as will be discussed in detail in section 3, the negative energies mentioned above can be relevant in purely classical regimes.. After all, unbridled negative energies generate disturbing dynamical instabilities. In addition, when gravitational effects are considered, unconstrained negative energies allow the formation of traversable wormholes , “faster than light travel” , naked singularities , and closed timelike curves (aka time machines) . Most importantly for the present work, the foundation of black hole thermodynamics rests on the classical area theorem  which states that, if the energy of matter is non-negative, the area of black hole horizons cannot decrease. In the presence of negative energies, the entropy associated with black holes can decrease, and the conservative physicist would suspect that this will lead to violations of the second law of thermodynamics.

Clearly then, it is important to reconcile the bold and conservative points of view and to discover what (if any) special properties particular negative energy objects may possess that render them immune to such concerns. To date, discussions in the literature have concentrated on the fact that orientifolds have ‘orbifold boundary conditions,’ meaning that the perturbative description of an orientifold involves a quotient of some smooth spacetime by a discrete $`Z_2`$ symmetry. For an orientifold $`p`$-plane, this is a $`Z_2`$ symmetry that leaves a $`(p,1)`$ spacetime surface invariant and reflects the other directions about this surface. Note that this quotient implies that the spacetime fails to be asymptotically flat in the usual sense and that the change of boundary conditions forbids pair production of orientifolds from the vacuum222Here we refer to orientifold points and non-compact or non-contractible orientifolds. Presumably, contractible orientifold planes must have positive energy to avoid this problem since their presence is compatible with the usual boundary conditions. This is a plausible result, though we are not aware of any literature on the subject.. Certain other important aspects of orientifold physics can be captured by working on the two-fold covering space and requiring the physics to be $`Z_2`$ symmetric. In particular, this removes the perturbative negative energy oscillations of the $`p`$-plane. Since the $`Z_2`$ quotient resolves at one stroke two potential instabilities of concern to our conservative physicist, it provides much encouragement to bold physicists wishing to use such objects in brane-world scenarios.

While it may be tempting to regard the creation of closed time-like curves, wormholes, and naked singularities as intriguing possibilities rather than physical inconsistencies, few physicists would be willing to compromise the second law of thermodynamics. Indeed, a famous quote from Albert Einstein  states his belief that thermodynamics was the most likely element of physics to survive without change – in particular, that our understanding of thermodynamics is more reliable than either that of quantum mechanics or relativity. Presumably, it is also more reliable than our understanding of string theory.

It was argued in  that the orbifold boundary condition alone is not sufficient to enforce the second law for negative energy branes. We should first remark that although one can in principle violate the second law using negative energy and ordinary matter by, for example, throwing negative energy into a star and thereby lowering its entropy and temperature, such concerns are in fact removed for ordinary matter if the negative energy is confined to an orbifold plane. The association of the negative energy with the unusual boundary conditions makes it impossible for the orbifold to be destroyed by ordinary matter and thus forbids its negative energy from thermalizing and lowering the temperature of our star.

On the other hand, the issue becomes much more subtle when black holes are involved. Unlike normal matter, black holes absorb the energy of any object that passes inside them. They also have the capability to hide the would-be orbifold plane. The usual tool for considerations of black hole thermodynamics is the Raychaudhuri equation , which shows that a negative energy object crossing the horizon of a black hole in a region of negligible shear necessarily causes null rays to defocus, so that the horizon area (and thus, at first sight, the entropy) tends to decrease. While one does not typically describe objects with orbifold boundary conditions as ‘moving,’ it is straightforward to instead consider a black hole that is thrown toward our negative tension brane. The resulting collision is the same in either case, and corresponds to a symmetric collision of two black holes with a negative tension object in the covering space. The second black hole complicates the analysis somewhat but, depending on the co-dimension $`q`$ of the negative tension brane, a scaling analysis shows that choosing the black holes to be either very large (for $`q=1`$) or very small (for $`q>2`$) compared to the gravitational length scale of the negative tension brane makes the mass of this second black hole negligible in comparison with the amount of negative energy absorbed from the brane. In such regimes, the second black hole should not interfere strongly with the argument above.

In addition to the concern that the negative energy object will cause defocusing of light rays and produce a local decrease in the horizon area, one might have the simpler concern that absorbing a negative energy object would reduce the mass of the black hole so that the final equilibrium configuration would be expected to have smaller horizon area than the initial one. This is essentially a repeat performance of our star argument above, with the notable different that the black hole should in fact be able to thermalize the negative energy of an orientifold. However, one can easily see that this concern is groundless in most cases. A negative energy object is gravitationally repulsive, so that in order to collide, the brane and black hole must be impelled toward each other with a finite kinetic energy333That this kinetic energy is necessarily positive is most easily seen by considering the negative tension brane to be unmoving and thinking of the kinetic energy as associated with the relative motion of our black hole and its image in the two-fold covering space.. Since the energy of any object vanishes when held static at a black hole horizon, this kinetic energy must be as large as the rest mass of the negative tension brane in order for the brane to just reach the horizon of the black hole. In general then, the collision occurs only if the kinetic energy is large enough to counterbalance the negative contribution to the mass of the final object. While there remains a concern about black hole entropy during the collision itself, one therefore expects the final equilibrium state to have larger energy and thus larger entropy than did the initial state.

When the negative tension brane has co-dimension 1 (as considered in  and ), however, this argument breaks down. Such branes are gravitationally attractive. In these cases, one needs to address the question of the entropy of the final state. The O8 plane and AdS soliton provide examples of this type; they are discussed in section 3.2 and 4.

The authors of  explored collisions between a black hole and a negative tension brane in a low dimensional toy model. Interestingly, they found that entropy considerations became moot as such collisions necessarily gave rise to a spacelike (big-crunch like) singularity that engulfed the entire space and prevented any notion of a final equilibrium configuration. Nevertheless, one intuitively suspects that this non-perturbative dynamical instability is linked to the above concerns about the second law. It is interesting to note that the instability of  remains when the system is embedded (as is straightforward to do444We thank Jon Bagger and Raman Sundrum for this observation.) in a supersymmetric system of the sort considered in . An investigation is currently in progress  of whether the instability might be an artifact of the low dimensional setting used in , but preliminary results for a 2+1 dimensional negative tension brane in $`AdS_{3+1}`$ indicate that the instability will survive in higher dimensions .

Such a dynamical instability would be as undesirable for physical models as violations of the second law. We are therefore compelled to seek an alternative resolution of the problems identified above for the stringy negative tension objects. We will see that this study uncovers a new and perhaps surprising property of (most) orientifolds that removes any concern about violations of black hole thermodynamics. This property is related to an unusual set of couplings between the gravitational field and massive modes associated with various classes of non-perturbative stringy effects. Because they involve the curvature, such couplings guarantee that black hole entropy does not reduce strictly to the horizon area  for a black hole in the neighborhood of an orientifold. This allows the entropy of a black hole to increase even though its area decreases as the horizon encounters a negative tension brane.

For each known case of a perturbatively stable stringy negative energy object, we find some novel mechanism that protects the system from violations of the second law. In most cases, this mechanism is the one just mentioned. However, the D8-brane and the AdS soliton  (which has already been used in brane-world constructions ) constitute interesting exceptions. For each of them, we uncover a distinct mechanism for consistency with the second law. One suspects that stable consistent negative energy objects can exist only when they have similarly clever properties and a natural first conjecture is that these three examples provide a complete list. One therefore expects that while more simplified phenomenological descriptions of negative tension objects (e.g., of the sort used in ) can point the way to intriguing and useful new phenomena, they remain at best a rough approximation to the physics involved and when studied in detail are likely to have instabilities of the sort described in .

We begin our analysis in section 2 with the case of the orientifold six-plane (O6-plane). This case has the advantage of a known strong coupling description as a purely gravitational soliton  (the Atiyah-Hitchin manifold ) in M theory. In such a manifestly classical setting it is straightforward to explore issues of black hole entropy and non-perturbative gravitational couplings in detail. As we will describe, at strong coupling the ten-dimensional description of the O6-plane necessarily involves not only the massless fields of type IIA supergravity but also the 11-dimensional Kaluza-Klein modes, which in ten-dimensions are described by an infinite tower of fields at integer multiples of the D0-brane mass. These fields are associated with the well known ‘exponential corrections’ that deform the negative mass Euclidean Taub-NUT manifold  to the smooth Atiyah-Hitchin manifold, and which can be seen in a field theoretic treatment of a D2-brane probe near an O6-plane. We argue that there is clearly no violation of the second law in the M theory description, so curvature couplings to these non-perturbative corrections resolve the puzzle as outlined above.

We then consider other negative energy orientifolds in section 3. This section begins with a scaling analysis to determine whether one would expect these objects to be subject to similar effects. The answer is in the affirmative for most cases, and we find in particular that for the strongly coupled O5-plane they must again come from fields at a non-perturbative mass scale (in this case, the scale set by the D-string tension). We argue that such effects are associated with the onset of a Kosterlitz-Thouless transition  in the field theory of a D-string probe near an O5-plane. At strong coupling, the associated instanton/anti-instanton pairs are unbound and produce non-perturbative corrections to the corresponding moduli space metric.

We turn in section 3.2 to the O8-plane, which forms the single exceptional orientifold case. Instead of relying on non-minimally coupled massive fields, the O8-brane turns out to avoid difficulties through the requirement that it appear in the vicinity of positive tension D8-branes.

A final negative energy object is the ADS soliton, which we address in section 4. This can be considered either as a purely gravitational negative energy soliton or, to a certain extent, as a negative energy brane in one lower dimension. Despite the similarity to the O6-plane and its Atiyah-Hitchin description in M-theory, this negative tension brane avoids thermodynamic difficulties through a completely different effect. We will see that it is simply impossible for such a negative tension brane to cross the event horizon of a black hole! Interestingly, this effect is necessarily associated with a vanishing (Einstein frame) worldvolume metric on the lower-dimensional brane, again leading to a breakdown of a description in terms of pure Einstein-Hilbert gravity coupled to a negative tension brane. See , however, for how to use such solitons to address the hierarchy problem in brane-world settings.

## 2 The orientifold six-plane

For most orientifolds, we wish to argue that the spacetime description receives corrections which break the relationship between entropy and area for nearby black holes. In this section, we will explore in detail the case of an O6-plane555In particular, the O$`6^{}`$-plane in the notation of . in the strongly-coupled IIA theory. The effective description at low energies is then eleven-dimensional supergravity compactified on a circle of radius $`R_{11}gl_s`$ where $`g`$ is the string coupling and $`\mathrm{}_s`$ is the string length. In this case, one expects the dominant corrections to the ten-dimensional geometry to come from the massive Kaluza-Klein modes associated with this compactification (and which become light in the strong coupling limit). The advantage of focusing on this case is that non-perturbative effects associated with $`R_{11}`$ are fully encapsulated in the eleven-dimensional supergravity action. Thus, we gain a much more detailed understanding of the effects of the corrections than in cases where the dominant corrections come from, say, massive string modes. Note that because adding additional orientifolds changes the boundary conditions, one cannot suppress corrections by considering an arbitrary number of coincident O6’s, as one can with D-branes.

Recall that the tension of an O6-plane is

$$\tau _6=\frac{2}{(2\pi )^6g\mathrm{}_{s}^{}{}_{}{}^{7}},$$
(1)
so that the radius of curvature of the ten-dimensional metric is of order

$$RG_{10}\tau _6gl_s,$$
(2)
where $`G_{10}`$ is the ten-dimensional Newton’s constant. Hence, $`R`$ is of order $`R_{11}`$, but $`R`$ is much larger than $`l_s`$, as well as the length scales $`l_p^{(10)}=g^{1/4}l_s`$ and $`l_p^{(11)}=g^{1/3}l_s`$ associated with the ten- and eleven-dimensional Planck lengths. Thus, while Kaluza-Klein corrections are relevant and an eleven-dimensional description is needed to capture the full physics, no further corrections need be considered.

Below, we first review how the eleven-dimensional geometry can be obtained by stringy methods (section 2.1) and then describe the geometry itself (section 2.2) which provides the resolution of our concerns from the eleven-dimensional perspective. We then discuss the entropy (section 2.3) and energy (section 2.4) from the ten-dimensional point of view.

### 2.1 D2-brane probe calculation

We now briefly review how to obtain the complete eleven-dimensional geometry from a probe D-brane calculation , though we refer the reader to the original references for further details.

Consider a probe D2-brane with its worldvolume aligned parallel to two of the directions in the orientifold plane. A D2-brane in the orientifold background will carry a worldvolume field theory with $`SU(2)`$ gauge invariance and $`N=4`$ supersymmetry in $`2+1`$ dimensions. There are seven scalar fields which parametrize the directions orthogonal to the D2-brane in the ten-dimensional geometry which split naturally into four fields corresponding to the directions along the O6-plane (which decouple from the other fields) and three fields $`\varphi _i,i=1,2,3`$ (corresponding to the overall transverse space). These scalars transform in the adjoint of $`SU(2)`$, and their potential energy is minimized when they all commute. At a generic minimum of the potential, the gauge group is broken to $`U(1)`$ and, since we are in 2+1 dimensions, this massless gauge field can be dualized to give an additional scalar $`\sigma `$. Since it arises from a $`U(1)`$ gauge field, $`\sigma `$ is periodically identified, and $`\sigma `$ in fact corresponds to the M-theory circle. We therefore have a four-dimensional Coulomb branch for this field theory for which the metric on moduli space (cross the seven flat directions along the O6-plane) is just the ten-dimensional spatial geometry of the O6-plane666Here we assume that the three-form potential vanishes and that the eleven-dimensional geometry is of the form $`\times _{10}`$ (a direct product of metrics) where $``$ is the time direction. This assumption will be justified a posteriori by the fact that the result satisfies the eleven-dimensional supergravity equations of motion..

At the classical field theory level, the metric on the space of $`\varphi _i,\sigma `$ is just the flat metric on $`(^3\times S^1)/_2`$, where the $`_2`$ identification arises from the action of the Weyl group of $`SU(2)`$. Supersymmetry guarantees that the moduli space cannot be lifted by quantum corrections and that the metric on moduli space remains hyper-Kähler. By studying the action of the symmetries one can show that the one-loop contribution generates a cross term between the $`\varphi _i`$ and $`\sigma `$. In terms of the ten-dimensional string theory, this correction corresponds to the Ramond-Ramond charge carried by the O6-plane. The hyper-Kähler structure implies that there are no new corrections at higher orders in perturbation theory, but there may be non-perturbative corrections. Indeed, one can identify instantons associated with D0-brane exchange between the probe and the orientifold that yield an infinite series of corrections exponentially suppressed by powers of the D0-brane mass. Since there is no Higgs branch for this theory, the metric on the Coulomb branch must be non-singular ; together with the supersymmetry and asymptotic structure, this is enough to determine the metric uniquely. The appropriate metric is the Atiyah-Hitchin metric of . Thus, the flat metric receives a one-loop correction which accounts for the effects of the bulk Ramond-Ramond gauge field, and instanton corrections which remove the singularity at the origin to give the Atiyah-Hitchin metric. The eleven-dimensional metric for the O6-plane is then this metric cross a flat $`^{6,1}`$.

### 2.2 The Atiyah-Hitchin metric

The Atiyah-Hitchin metric is a non-singular $`SO(3)`$ symmetric hyper-Kähler manifold, which also appears in the study of a two-monopole moduli space. The metric takes the form

$$ds^2=\frac{b^2}{r^2}dr^2+a^2\sigma _1^2+b^2\sigma _2^2+c^2\sigma _3^2,$$
(3)
where

$`\sigma _1`$ $`=`$ $`\mathrm{sin}\psi d\theta +\mathrm{cos}\psi \mathrm{sin}\theta d\varphi ,`$
$`\sigma _2`$ $`=`$ $`\mathrm{cos}\psi d\theta +\mathrm{sin}\psi \mathrm{sin}\theta d\varphi ,`$
$`\sigma _3`$ $`=`$ $`d\psi +\mathrm{cos}\theta d\varphi ,`$ (4)
and $`a,b,c`$ are functions only of $`r`$. The ranges of the coordinates are $`0\theta \pi `$, $`0\varphi 2\pi `$, $`0\psi 2\pi `$, and $`\pi R_{11}/2r<\mathrm{}`$. In addition, an identification is imposed under the $`_2`$ symmetry

$$\theta \pi \theta ,\varphi \varphi +\pi ,\psi \psi .$$
(5)
At large $`r`$, the functions $`a,b,c`$ take the form

$`a`$ $`=`$ $`r\left(1{\displaystyle \frac{R_{11}}{r}}\right)^{1/2}{\displaystyle \frac{8r^2}{R_{11}}}\left(1{\displaystyle \frac{R_{11}^2}{8r^2}}\right)e^{\frac{2r}{R_{11}}}+𝒪(e^{\frac{4r}{R_{11}}}),`$
$`b`$ $`=`$ $`r\left(1{\displaystyle \frac{R_{11}}{r}}\right)^{1/2}+{\displaystyle \frac{8r^2}{R_{11}}}\left(1+{\displaystyle \frac{R_{11}}{r}}{\displaystyle \frac{R_{11}^2}{8r^2}}\right)e^{\frac{2r}{R_{11}}}+𝒪(e^{\frac{4r}{R_{11}}}),`$
$`c`$ $`=`$ $`R_{11}\left(1{\displaystyle \frac{R_{11}}{r}}\right)^{1/2}+𝒪(e^{\frac{2r}{R_{11}}}),`$ (6)
where $`2\pi R_{11}`$ is the asymptotic circumference of the $`S^1`$. The above two features together imply that the geometry asymptotically approaches the flat metric on $`(^3\times S^1)/_2`$, with $`\psi `$ the coordinate on the $`S^1`$. Thus, it is the angle $`\psi `$ that corresponds to the M-theory circle and the reduction on this circle to ten dimensions gives a magnetic charge under the Kaluza-Klein gauge field.

The exponentially small corrections at large $`r`$ represent the effect of the instantons described in the preceding section. Because of these corrections, the metric closes off in a regular fashion at $`r=\pi R_{11}/2`$, where $`b=c=\pi `$, and $`a2r\pi R_{11}`$. There is no region with $`r<\pi R_{11}/2`$. For more details, see .

For our purposes, the important point is that the eleven-dimensional geometry is a completely smooth vacuum solution; the above metric is Ricci flat. Despite the vacuum geometry, this spacetime actually has a negative ADM mass (tension). The failure of the positive mass theorems  is associated with the unusual boundary conditions which, for example, rule out the existence of the asymptotically constant spinors used in the Witten proof . The complete description of the collision of a black hole with the O6-plane is then given by a black string wrapped on the compact direction in 11d ‘moving in the Atiyah-Hitchin space’777Such language is strictly appropriate only for a perturbatively small black string, but it is convenient to use this phrase to refer to large black strings as well.. Since this spacetime satisfies the weak energy condition — indeed, there is no local stress-energy of any kind — in the 11d description, the usual 11-dimensional area theorem tells us that the area, and hence entropy, of this black string solution is non-decreasing. Thus, once we take the exponential corrections fully into account, the apparent violation of the second law is completely resolved.

### 2.3 Entropy from the ten-dimensional point of view

The eleven-dimensional picture removes concern about the second law and provides a very satisfying result. However, it is useful to explore the mechanism by which the second law is saved from the ten-dimensional point of view. To flesh out the details, we now re-examine the Atiyah-Hitchin metric from the ten-dimensional point of view and see how the naïve argument in the introduction is invalidated.

First of all, we note that the Atiyah-Hitchin metric (3) depends on $`\psi `$, so that $`_\psi `$ is not a Killing vector field. This implies that the dimensional reduction to ten dimensions will contain excited massive fields from the Kaluza-Klein modes. We wish to argue that the presence of these massive modes implies that the usual relation between the entropy and the area of the black hole horizon is modified. The essential point is that the formula $`S=A/4`$ is valid only in Einstein gravity minimally coupled to matter fields. Wald and Iyer undertook a detailed study of a family of actions with further curvature terms, and showed that the first law of black hole thermodynamics requires the definition of the entropy to be suitably generalized . The full entropy expression is obtained by differentiating the Lagrangian with respect to its dependence on the Riemann tensor, so that each curvature coupling provides a separate contribution to the entropy of any black hole.

In our case, the presence of the massive modes will produce such additional curvature terms. In general, the ten-dimensional action involving these massive modes is a complicated non-linear expression. However, these couplings become tractable when the massive modes can be treated as a small perturbation. This will be enough for our purposes as we seek only a qualitative understanding in ten-dimensional terms.

Writing the 11d metric as $`g_{ab}=g_{ab}^0+h_{ab}`$, the quadratic action for the perturbation $`h_{ab}`$ is

$$S=\frac{1}{32\pi G}d^{11}x\sqrt{g^0}h^{ab}A_{abcd}h^{cd},$$
(7)
where $`g_{ab}^0`$ has been used to raise the indices and

$`A_{abcd}`$ $`=`$ $`{\displaystyle \frac{1}{4}}g_{cd}^0_a_b{\displaystyle \frac{1}{4}}g_{ac}^0_d_b+{\displaystyle \frac{1}{8}}(g_{ac}^0g_{bd}^0+g_{ab}^0g_{cd}^0)_e^e+{\displaystyle \frac{1}{2}}R_{ad}^0g_{bc}^0{\displaystyle \frac{1}{4}}R_{ab}^0g_{cd}`$ (8)
$`+{\displaystyle \frac{1}{16}}R^0g_{ab}^0g_{cd}^0{\displaystyle \frac{1}{8}}R^0g_{ac}^0g_{bd}^0+(ab)+(cd)+(ab,cd).`$
Here $`R_{ab}^0`$ ($`R^0`$) denotes the Ricci tensor (scalar) of the background spacetime $`g_{ab}^0`$. Similarly, $``$ is the covariant derivative associated with $`g_{ab}^0`$. We are mostly interested in the curvature couplings, which can be rewritten as

$$S_{curv}=\frac{1}{32\pi G}d^{11}x\sqrt{g^0}R_{ab}^0\left(2h^{ac}h_c{}_{}{}^{b}h^{ab}h_c{}_{}{}^{c}+\frac{1}{4}h_c{}_{}{}^{c}h_{d}^{}{}_{}{}^{d}\frac{1}{2}h^{ab}h_{ab}\right).$$
(9)
We take the background $`g_{ab}^0`$ to be independent of $`\psi `$ and expand $`h_{ab}`$ in Fourier modes,

$$h_{ab}=\underset{m0}{}h_{ab}^{(m)}e^{im\psi }.$$
(10)
Note that we have taken the perturbation of the zero-mode to vanish, as the perturbation is to represent only the massive fields. It is for this reason that terms in the action that are linear in $`h_{ab}`$ integrate to zero regardless of whether $`g_{ab}^0`$ is taken to be a solution of the equations of motion. Reduction of (9) to ten-dimensions is straightforward once a few conventions are fixed. We choose to work with what we call the “naïve Einstein metric” $`g^{10,E}`$ (and the associated dilaton $`\varphi `$ and 1-form $`A_a`$) in ten-dimensions which is given by the usual algebraic expressions in terms of $`g_{ab}^0`$:

$$g_{ab}^0=e^{\varphi /6}g_{ij}^{E,10}dx^idx^j+e^{4\varphi /3}(d\psi +A_idx^i)^2.$$
(11)
Here we use $`i,j`$ to represent ten-dimensional indices. We use the phrase ‘naïve Einstein metric’ because the associated action is readily seen to contain the curvature couplings

$`S_{curv}`$ $``$ $`{\displaystyle d^{10}xe^{\varphi /6}\sqrt{g^{E,10}}R_{ij}^{(E,10)}x_{,a}^ix_{,b}^j}`$ (12)
$`\times `$ $`{\displaystyle \underset{m}{}}\left(2h^{ac(m)}h_c{}_{}{}^{b(m)}h^{ab(m)}h_c{}_{}{}^{c(m)}+{\displaystyle \frac{1}{4}}h_c{}_{}{}^{c(m)}h_{d}^{}{}_{}{}^{d(m)}{\displaystyle \frac{1}{2}}h^{ab(m)}h_{ab}^{(m)}\right),`$
where $`R_{ij}^{(E,10)}`$ is the Ricci tensor of $`g_{ij}^{E,10}`$ and the matrix of derivatives $`x_{,a}^i`$ implements the pull back of ten-dimensional forms to eleven dimensions.

While this calculation is only valid when the massive modes are a small perturbation, we wish to suggest that similar considerations apply when the massive modes are of order one, as in the region near $`r=\pi R_{11}/2`$ in the Atiyah-Hitchin metric. In particular, we suggest that the terms in (7) and the corresponding higher corrections will continue to source important curvature terms in the ten-dimensional action, which will produce order-one corrections to the entropy formula. The full entropy must after all give just the eleven-dimensional area of the black hole and there is no reason in the presence of massive excitations for this to agree with the area as measured in ten-dimensions.

Let us therefore summarize our picture of the collision of a black $`p`$-brane and an orientifold plane. When the black $`p`$-brane is absorbing the orientifold plane, massive Kaluza-Klein modes will be excited near the horizon. What happens to these modes subsequently? If the black $`p`$-brane completely engulfs the orientifold plane, we would not expect the black $`p`$-brane to be capable of supporting ‘massive hair’. In the linearized analysis, it is easy to see that the massive modes must die off exponentially, so long as the black $`p`$-brane has a radius $`r_0`$ satisfying $`r_0R_{11}`$ (so that it is a truly ten-dimensional solution). This is because in eleven dimensions we are considering just perturbations of the black string of the form $`h_{ab}e^{im\psi }`$. But these are the perturbations considered by Gregory and Laflamme , who showed that they will decay exponentially in time for all $`m>m_{}R_{11}/r_0`$. Thus, the massive modes are absorbed by the black $`p`$-brane, which settles down to a regular solution carrying just the magnetic Kaluza-Klein gauge charge.

From the ten-dimensional point of view, the picture we want to suggest is thus that the orientifold will first fall across the event horizon, lowering the area of the black $`p`$-brane. However, the massive Kaluza-Klein fields will then have order-one excitations in the vicinity of the horizon, so the area will differ substantially from the entropy. These massive modes will then be absorbed by the black $`p`$-brane. Presumably, this influx carries positive energy into the black brane and raises its area so that the final state will be a black brane with energy and horizon area greater than the initial one. Since the massive modes die away, the ten- and eleven-dimensional areas must once again agree in the final state. One can then see that the final area is greater from the eleven-dimensional picture where there is no negative energy. Alternatively, one can use the general argument in the introduction that the kinetic energy must be at least big enough to overcome the negative energy from the orientifold tension.

### 2.4 Energy from the ten-dimensional point of view

We have argued that the description of the horizon dynamics should be rather different from the 10- and 11-dimensional points of view, with the 10-dimensional horizon shrinking during the collision while the 11-dimensional horizon grows. We then suggested that this difference can be traced to the massive modes. In section 2.3 above we showed that these massive modes do induce a series of perturbative corrections to black hole entropy. We now verify the flip-side of our suggestion by showing that the massive modes lead to perturbative negative energy fluxes from the ten-dimensional point of view and so can be associated directly with the defocusing of null rays and the shrinking of black hole horizons. In particular, we argue that at finite coupling the negative energy is not confined to an infinitely thin brane, but instead spreads out over a finite region of space. We again study these excitations in the perturbative regime where there is a clear distinction between the massive fields and the massless IIA fields. For the solution (3), this corresponds to the region $`rR_{11}`$.

We wish to show that Kaluza-Klein reduction of (3) to 10-dimensions yields massive modes which violate the weak null energy condition ($`k^ak^bT_{ab}>0`$) for some null vectors $`k`$. Here, we speak of the stress-energy tensor associated with the naïve Einstein metric $`g^{E,10}`$ introduced above.

If the massive fields were the only sources in ten-dimensions, the sign of the massive mode energy would be related by the Einstein equations to the sign of $`R_{ab}^{E,10}k^ak^b`$, which must be positive to prove the usual black hole area theorem. In the present setting, the massless ten-dimensional matter fields (the dilaton and the Ramond-Ramond 1-form) are also excited and, in fact, dominate over the massive modes at large $`r`$. The stress-energy tensor of these massless fields satisfies the weak energy condition, and therefore so must the full stress-energy tensor in any perturbative treatment. However, our goal here is to analyze the contribution of the massive modes $`T_{ij}^{massive}`$ in particular, which can be obtained by simply subtracting the massless contribution $`T_{ij}^{massless}`$ from the full stress-energy tensor $`T_{ij}^{full}`$. That this difference has a simple representation in terms of eleven-dimensional quantities is readily seen by using the ten-dimensional Einstein equations (which include the effects of the massive modes) to express $`T_{ij}^{massive}`$ in terms of a variation of the massless IIA action:

$`k^ik^jT_{ij}^{massive}`$ $`=`$ $`k^ik^jT_{ij}^{full}k^ik^jT_{ij}^{massless}={\displaystyle \frac{1}{8\pi G_{10}}}k^ik^jR_{ij}^{(E,10)}k^ik^jT_{ij}^{massless}`$ (13)
$`=`$ $`{\displaystyle \frac{2}{\sqrt{g^{E,10}}}}k^ik^j{\displaystyle \frac{\delta }{\delta g^{(10,E)ij}}}S_{masslessIIA}={\displaystyle 𝑑\psi \frac{e^{\varphi /6}\sqrt{g^{11}}}{16\pi ^2G_{10}\sqrt{g^{10,E}}}k^ik^jR_{ij}^0},`$
where we have assumed that $`k`$ is null with respect to both the ten- and eleven-dimensional metrics ($`g_{ab}^{E,10}`$ and $`g_{ab}^0`$) and the symbol $`R_{ab}^0`$ refers to the Ricci tensor of the eleven-dimensional zero-mode metric $`g_{ab}^0`$.

As a result, we need only analyze the sign of $`𝑑\psi k^ak^bR_{ab}^0`$ for an appropriate null vector $`k`$. Note that although the full eleven-dimensional metric $`g_{ab}=g_{ab}^0+h_{ab}`$ is Ricci flat in the Atiyah-Hitchin case, the zero mode piece $`g_{ab}^0`$ need not be and $`k^ak^bR_{ab}^0`$ need not vanish. The setting here is similar to that often used to treat gravitational radiation in which we separate out a gravitational perturbation $`h_{ab}`$ and choose not to regard it not as part of the metric (represented by $`g_{ab}^0`$) but instead as a ‘matter field’ which in turn acts as a gravitational source.

Nonetheless, the vanishing of the full Ricci tensor $`R_{ab}^{(11)}`$ is still quite useful. In particular, we can expand $`R_{ab}^{(11)}`$ as

$$0=R_{ab}^{(11)}=R_{ab}^0+R_{ab}^{0,1}+R_{ab}^{0,2}+\mathrm{},$$
(14)
where $`R_{ab}^{0,1}`$ and $`R_{ab}^{0,2}`$ are respectively linear and quadratic in the perturbation $`h^{ab}`$. Since $`h^{bc}`$ has no zero-mode part, the integral of $`R_{ab}^{0,1}`$ over $`\psi `$ vanishes and to second order in the perturbation we have

$$𝑑\psi k^ak^bR_{ab}^0=𝑑\psi k^ak^bR_{ab}^{0,2}.$$
(15)
We may then consult  to find888The reference  contains the result for perturbations around Minkowski space. However, variations of the Ricci tensor can be written in terms of only the metric perturbation $`h_{ab}`$, its covariant derivatives, and the unperturbed Ricci tensor. Thus, since terms involving $`R_{ab}^0`$ and two $`h`$’s are higher order, $`R_{ab}^{0,2}`$ can be obtained from the flat space result by simply replacing partial derivatives with covariant derivatives and checking the ordering of the second covariant derivatives against (8) above. the relation

$`R_{ab}^{0,2}={\displaystyle \frac{1}{2}}h^{cd}_{(a}_{b)}h_{cd}{\displaystyle \frac{1}{2}}h^{cd}_a_ch_{bd}{\displaystyle \frac{1}{2}}h^{cd}_b_ch_{ad}+{\displaystyle \frac{1}{4}}(_ah_{cd})_bh^{cd}`$ (16)
$`+(^dh_b^c)_{[d}h_{c]a}{\displaystyle \frac{1}{2}}_d(h^{dc}_ch_{ab}){\displaystyle \frac{1}{4}}(^ch)_ch_{ab}(_dh^{cd}{\displaystyle \frac{1}{2}}^ch)_{(a}h_{bc)},`$ (17)
where we have used the fact that terms involving $`R_{ab}^0`$ are higher order in $`h_{ab}.`$ Here indices are raised and lowered with the zero-mode metric $`g_{ab}^0`$ and $`h=g^{0,ab}h_{ab}`$ is the trace of the perturbation in this metric.

It is now straightforward to evaluate $`k^ak^bR_{ab}^{0,2}`$ for various null vectors $`k`$ using the large $`r`$ expansion of (3). We choose $`k=k^t\frac{}{t}\pm k^r\frac{}{r}`$, which is clearly relevant to the horizon of an approaching black hole. The reader may check that such $`k`$ are null simultaneously for both $`g_{ab}^{E,10}`$ and $`g_{ab}^0`$. If we introduce

$$\mathrm{\Delta }^2b^2a^2=\frac{2^7r^3}{R_{11}^3}e^{2r/R_{11}}\left(1+𝒪(\frac{r}{R_{11}})\right),$$
(18)
a somewhat tedious calculation yields

$$\frac{1}{2\pi }𝑑\psi k^ak^bR_{ab}^0=(k^r)^2\frac{3R_{11}^2}{32r^4}\mathrm{\Delta }^4\left(1+𝒪(\frac{r}{R_{11}})\right)<0.$$
(19)
As a result, the stress-energy tensor of the massive modes does indeed violate the weak energy condition in our naïve ten-dimensional Einstein frame.

Of course, this does not mean that an approaching black hole begins to contract while the orientifold is still far away. The energy in the large $`r`$ regime is is dominated by the massless modes and must therefore be positive. However, we expect the massive modes to dominate once $`rR_{11}`$, and there is no reason to expect the sign of the energy carried by such modes to change at that point. Thus, the calculation above indicates that at finite coupling the full stress-energy violates the weak energy condition over a region of size $`R_{11}`$. When the black hole encounters this region, its horizon will begin to contract.

## 3 Other Orientifold-planes

We now want to consider the other Orientifold-planes of type II string theory. These may be classified in much the same manner as the familiar positive tension branes in terms of the bulk gauge fields under which they carry charge. The various categories include the O$`p`$-planes for $`p=0,1,\mathrm{}8`$ (which have the same supersymmetries and the same kinds of charges as the D$`p`$-branes), the ONS5- and OF1-planes which couple to the Neveu-Schwarz B-field, and the OP1 orientifold line which is not charged under a gauge field but does carry momentum.

For the strongly coupled O6-plane, corrections from massive Kaluza-Klein modes were important because the curvature radius $`R`$ of the ten-dimensional metric (the ‘orientifold scale’) was of the same order as $`R_{11}`$ while other corrections (e.g., $`\alpha ^{}`$ corrections) could be ignored (since e.g. $`R\mathrm{}_s`$). To discover which corrections may be important for the other O-planes, we should similarly compare their geometric scales to the fundamental length scales.

In the IIA cases, these fundamental scales are the string length $`\mathrm{}_s`$, the 10-d Planck length $`\mathrm{}_p^{(10)}`$, and the scales of M theory: $`R_{11}`$ and the 11-d Planck length $`\mathrm{}_p^{(11)}`$. While $`R_{11}`$ has no direct analogue in the IIB theory, distance scales associated with the D-brane tensions lead to similar effects. The fundamental scales are related through:

$$\mathrm{}_s=g^{1/4}\mathrm{}_p^{(10)}=g^{1/3}\mathrm{}_p\mathrm{}_p^{(10)}=(\mathrm{}_p^{(10)})^{9/8}R_{11}^{1/8},R_{11}=g\mathrm{}_s=g^{2/3}\mathrm{}_p^{(11)}.$$
(20)
For the (Ramond-Ramond) O$`p`$-planes, the scale is set by the D-brane tension

$$T\frac{m}{L^p}\frac{1}{g\mathrm{}_s^{p+1}}.$$
(21)
As a result, the radius of curvature $`R`$ of the 10-d metric for such cases is

$$R^{7p}G_{10}Tg\mathrm{}_s^{7p}.$$
(22)
Consider first weak coupling. Then the string length is the largest scale, and this is at least as big as the orientifold scale except for $`p=8`$. Thus for $`p8`$, stringy corrections can invalidate the pure supergravity description. Since these orientifolds preserve half the supersymmetry, their interactions may be protected from perturbative corrections. However, we saw explicitly that non-perturbative corrections arise in the O6 case. Dualities then imply that some corrections should be present for the other O$`p`$-planes as well. The O8-plane is an exception to this general picture; it can have no significant corrections at weak coupling and requires some new mechanism to protect the second law. The resolution of this case is discussed in section 3.2 below.

At strong coupling, the situation is reversed and $`R/\mathrm{}_sg^{1/(7p)}`$ is large for $`p<7`$, so that $`O(\alpha ^{})`$ corrections are negligible for these cases. We now need to consider the other fundamental scales, which are larger than the string scale at strong coupling. We can rewrite (22) as

$$R^{7p}g^{(p3)/4}[\mathrm{}_p^{(10)}]^{7p},$$
(23)
so that $`R/\mathrm{}_p^{(10)}`$ is also large at strong coupling and the associated corrections are small for $`p>3`$. In terms of the 11-d Planck length we have

$$R^{7p}g^{(p4)/3}\mathrm{}_p^{7p},$$
(24)
so that 11-d quantum corrections are important for $`p=4`$. Thus, the only cases where quantum corrections are not important are $`p=6`$, considered in the previous section, and $`p=5`$. We will discuss $`p=5`$ in section 3.1 below, and argue that it receives non-perturbative corrections from D-string effects that are similar to those described in section 2 for the orientifold 6-plane.

There are also a few NS-type O-planes to consider. For the ONS5, $`G_{10}T\mathrm{}_s^2`$, and string corrections are always important. For the OF1, we have $`G_{10}Tg^2\mathrm{}_s^6g^{1/2}\mathrm{}_p^{(10)}{}_{}{}^{6}\mathrm{}_{pl}^6`$. Thus, in the IIA theory, 11-d Planck corrections are always important while in IIB string corrections are important at weak coupling. At strong coupling in the IIB theory, the D-string scale $`\mathrm{}_{D1}=g^{1/2}\mathrm{}_s`$ (so that $`T_{D1}\mathrm{}_{D1}^2`$) becomes large enough that we have

$$G_{10}T_{OF1}g^2\mathrm{}_s^6=g^1\mathrm{}_{D1}^6\mathrm{}_{D1}^6,$$
(25)
so that D-string corrections will be important. The T-dual OP1 has a fixed line along a compact direction of radius $`R_{compact}`$, and a tension $`G_{10}Tg^2\mathrm{}_s^8/R_{compact}^2`$, which is smaller than the OF1 tension for $`R_{compact}>\mathrm{}_s.`$ As a result, either the same corrections are important as in the OF1 case above or $`R_{compact}\mathrm{}_s`$ and string scale corrections must again be considered.

Thus, a scaling analysis suggests that for any orientifold but the O8-plane the spacetime description would receive corrections analogous to the discussion in section 2 for the orientifold 6-plane. For the O5-plane, the only relevant corrections should arise at the scale of the D-string tension and we should be able to find them by instanton methods in a probe brane calculation analogous to that of section 2.1. We do so in section 3.1 below and then proceed to analyze the case of the O8-plane in section 3.2.

### 3.1 Nonperturbative corrections and the O5-plane

Let us now turn to the strongly coupled O5-plane. Our discussion below should apply to any of the variants of the O5-plane, as they differ only by discrete Neveu-Schwarz and Ramond-Ramond fluxes. While string and (ten-dimensional) Planck scale corrections are small, the effective description of strongly-coupled IIB is in fact the S-dual weakly-coupled IIB. Recall that there were no uncorrected IIB orientifolds at weak coupling. In fact, the Ramond-Ramond O$`p`$-plane is S-dual to the ONS5-plane, which had the potential for large stringy corrections at weak coupling. From the original strongly-coupled point of view, these should arise somehow from the D-strings which are becoming light. This shows that such corrections are allowed by a scaling analysis and our discussion of the second law predicts that they will be present despite the supersymmetry of the situation.

How do we verify the existence of these corrections? We recall that the non-perturbative corrections for the O6-plane can be seen in the worldvolume field theory of probe branes. A test D2-plane moving in the O6 background is described by an $`SU(2)`$ gauge theory on the D2, and the Atiyah-Hitchin metric is the associated moduli space metric. The massive Kaluza-Klein excitations in the Atiyah-Hitchin metric are related to instanton corrections to the corresponding moduli space.

We now attempt to repeat this analysis for the D1/O5-system. The relevant theory is now $`SU(2)`$ gauge theory in 2 Euclidean dimensions. Since the topology of infinity is $`S^1`$, the space of Euclidean solutions has instanton sectors associated with the winding of $`U(1)`$ subgroups around this $`S^1`$. The lowest action configuration in this sector should be a smooth Euclidean solution. Indeed, a massive version of the theory broken to a $`U(1)`$ subgroup is the abelian Higgs model whose vortex instantons are well-known.

In two Euclidean dimensions, the action of such vortex instantons is in general logarithmically divergent. Because our fields are massless, we expect this to be the case here. Nonetheless, such instantons can still contribute to the partition function. The point is that the action of an instanton/anti-instanton pair will be finite, as the separation $`s`$ of the pair will cut off the logarithmic divergence. Such pairs therefore have an action of size $`S\frac{1}{g}\mathrm{ln}s,`$ where we have explicitly indicated the fact that the action of any D-brane solution is proportional to $`\frac{1}{g}`$. On the other hand, the volume of phase space available for a pair with separation $`s`$ is proportional to the circumference of a circle with radius $`s`$. As a result, such instanton/anti-instanton pairs have an entropy that diverges as $`\mathrm{ln}s`$, without an the accompanying factor of $`\frac{1}{g}`$. Note that at large $`s`$ the pair will provide an approximate solution to the Euclidean equations of motion. It is therefore clear that, for large enough $`g`$, instanton/anti-instanton pairs with large $`s`$ will provide important instanton contributions to the partition function. What we have just described is of course just the usual Kosterlitz-Thouless phase transition  of two-dimensional Euclidean systems.

This analysis verifies our prediction and shows that the moduli space for D-strings moving in the O5-background receives corrections from instantons. This result and the analogy with the O6 system strongly suggests that we think of massive fields at scales set by the D1 tension as being excited near an O5-plane. Presumably, these fields carry negative energy and affect the entropy of black holes in much the same way as in the O6 case studied in section 2.

### 3.2 The O8-plane

Based on the O6-case, we generally expected non-trivial corrections to the (10d) supergravity description of the orientifolds. This is borne out, with the sole exception of the O8-plane at weak coupling. For this case, all of the corrections considered above are negligible. Furthermore, as it has co-dimension one, the O8-plane will attract a positive-mass black hole, so this case is also an exception to the argument in the introduction about kinetic energy being required to impel a black hole toward a negative tension brane. Since all the corrections are suppressed, the resolution of these puzzles must lie within classical (massive) type IIA supergravity.

Let us consider first the question of what happens when an event horizon encounters an O8-brane. We would like to suggest that the resolution of our thermodynamic concerns is related to a well-known difficulty in dealing with O8-planes and the related D8-branes. The solution corresponding to any of these objects in isolation has a singularity at a finite distance of order $`\frac{1}{g}\mathrm{}_s`$ from the brane. The usual interpretation of this singularity is that it is impossible to separate an O8- or D8-brane from an oppositely charged such brane by more than a distance of order $`\frac{1}{g}\mathrm{}_s`$. Thus, when we consider an O8-plane, we should remember that there must necessarily be nearby either an anti O8-plane (in which case the entire universe must be small as these are world-ending branes) or a D8-brane (which has positive tension). As a result, there is simply no setting where one may ask about a black hole (or in fact a black $`p`$-brane) which begins farther away from the O8-plane than the O8 curvature radius $`R\frac{1}{g}\mathrm{}_s`$ and collides with the negative tension orientifold without at the same time encountering a similar amount of positive energy material.

It is instructive to briefly consider black $`p`$-branes that begin close to the orientifold. Such a black hole will act much like its cousin in empty space if the size $`R_{bb}`$ of the black $`p`$-brane is much smaller than the O8-plane length scale $`R_{O8}=\frac{1}{g}\mathrm{}_s`$. If $`R_{bb}/R_{O8}`$ is non-negligible, we expect that distortions of this order will be present in the black brane horizon and consequently that the area $`A`$ of such a black brane will differ from the area $`A_0`$ of its flat space cousin by an amount of order $`A_0\frac{R_{bb}}{R_{O8}}`$. Furthermore, the exact distortion and the corresponding effect on the horizon area will depend on the precise location of the black brane, so that the horizon area of a black brane colliding with the O8-plane will change during the collision by an amount of this order in a manner that is beyond the scope of this work to predict.

Let us also estimate the magnitude of the supposed entropy reduction due to the collision of a black $`p`$-brane with the O8-brane. The amount of negative tension absorbed by the brane is roughly proportional to the tension $`T_{O8}`$ of the orientifold times the cross sectional area $`R_{bb}^{8p}`$ of the black brane, whereas the original tension of the black brane is of size $`\frac{R_{bb}^{7p}}{G_{10}}`$, where $`G_{10}`$ is the ten-dimensional Newton’s constant. Since $`G_{10}T_{O8}1/R_{O8}`$, we have a fractional change in tension of order $`\mathrm{\Delta }T/TR_{bb}/R_{O8}`$ and we would expect a fractional decrease in its area of the same order. But this is of the same size as the dynamical distortion term caused by the proximity of the O8-plane! These finite-distance corrections can therefore resolve the potential decrease in area. We suggest that this is exactly what happens, though it is difficult to study the effect in detail.

In the introduction, we argued that the total mass of the black hole typically could not be reduced by collision with a negative-energy source due to the repulsive gravitational potential. In the present case, this argument breaks down, because the O8-brane is of co-dimension one. However, this problem contains its own solution: since the O8-brane is co-dimension one while black branes always have co-dimension greater than one, a black $`p`$-brane cannot engulf the entire O8-brane. Instead, part of the O8 necessarily remains in the spacetime to source significant distortions of the black $`p`$-brane. Hence, although the mass of the black hole may decrease, one cannot conclude that the black hole’s area decreases when it collides with the O8-brane. Note also that the above scaling analysis suggests that a black hole cannot engulf enough O8 negative tension to bring its total mass to zero or below.

## 4 AdS Soliton

A final negative-energy object in string theory is the AdS soliton . Like the Atiyah-Hitchin metric, the AdS soliton is a pure gravitational solution and is not particularly stringy in and of itself. However, it is present in string theory and its construction was motivated by stringy considerations. As a result, this object is as connected with the consistency of string theory and the second law as the objects that we have previously discussed.

The AdS soliton is an asymptotically AdS solution of the vacuum Einstein equations with a cosmological constant and arises in string theory in the context of sphere compactifications. The metric for the AdS soliton in $`p+2`$ dimensions is

$$ds^2=\frac{r^2}{l^2}\left[\left(1\frac{r_0^{p+1}}{r^{p+1}}\right)d\tau ^2+(dx^i)^2dt^2\right]+\left(1\frac{r_0^{p+1}}{r^{p+1}}\right)^1\frac{l^2}{r^2}dr^2,$$
(26)
where $`i=1\mathrm{}p1`$, $`l`$ is the cosmological length scale, and the coordinate $`\tau `$ is compactified with period

$$\beta =\frac{4\pi l^2}{(p+1)r_0}.$$
(27)
We also take the $`x^i`$ to be compactified on a $`p1`$ dimensional torus of volume $`V_{p1}`$. This is a smooth metric, but it has negative energy,

$$E=\frac{r_0^{p+1}\beta V_{p1}}{16\pi G_{p+2}l^{p+2}},$$
(28)
where a definition of $`E`$ has been used in which the energy of AdS space vanishes.

If one dimensionally reduces along $`\tau `$, this looks like a $`p+1`$ dimensional spacetime with a negative-energy ‘brane’ at $`r=r_0`$. Now, for the metric (26) the proper length of the $`\tau `$ circle diverges at large $`r`$ so that $`p+2`$ dimensional physics is relevant at any value of $`\beta `$. On the other hand, one can make the physics truly $`p+1`$ dimensional by cutting of the spacetime at large $`r`$ as in . A collision with a black hole in the $`p+1`$ dimensional spacetime is then described by considering a black string wrapping the $`\tau `$ direction.

This example differs from the previous cases in that the amount of negative energy is a free parameter (controlled by $`\beta `$), so that it would be difficult to resolve our concerns using fundamental corrections to the supergravity description. In addition, like the O8-brane, the AdS soliton is co-dimension one when considered as a brane and one may check directly that it is gravitationally attractive.

Remarkably, the brane at $`r=r_0`$ cannot cross the black hole horizon in the usual sense. Since $`r=r_0`$ is the origin of the $`r,\tau `$ plane in the $`p+2`$ dimensional geometry, an $`S^1`$ of horizon generators would go to zero size if the black string horizon crossed $`r=r_0`$. Since this geometry is smooth, such a caustic cannot occur on an event horizon . However, the $`p+1`$ dimensional black hole is assuredly attracted to the negative-energy ‘brane’. So what happens? The only possibility is that the brane enters the horizon at a past endpoint, where new generators are entering the horizon and the black hole is still forming. One can see that this is so by again using the rotational symmetry. The point $`r=r_0`$ cannot cross the event horizon along a null generator, as there is no distinguished direction for the corresponding null ray to be travelling. Since the negative-tension brane does not intersect a pre-existing family of horizon generators, the argument in the introduction is simply inapplicable. Once again, in the $`p+2`$-dimensional description, the geometry is everywhere smooth so that the horizon area will not decrease. In particular, the final horizon area must be greater than any reasonable definition of an ‘initial area.’

Note that in contrast to recent work on higher-dimensional solutions , we don’t expect the horizon of our $`p+2`$ “black string” to be toroidal in any invariant sense; the picture here is more like , where one forms a horizon with initially toroidal spacelike slices, but the hole closes up faster than the speed of light. In the 3+1 analogue , there are no causal curves linking the torus999In our higher dimensional setting this question is moot as tori of our dimension link with surfaces and not curves. which extend to infinity (in agreement with topological censorship ). The cosmological constant (i.e., the $`r^2/l^2`$ factor in $`g_{\tau \tau }`$) will accelerate the wrapped black brane toward $`r=r_0`$ even if we place the string far away initially, so the ‘donut hole’ in the middle can close up before observers can pass through it.

Finally, we note an interesting property of this “brane” that is intimately associated with the resolution above. Because in $`p+2`$ dimensions the $`\tau `$ circle shrinks to zero size at $`r=r_0`$, the $`p+1`$ dimensional Einstein frame metric necessarily degenerates at the location of the brane. Such a singularity of the $`p+1`$ metric then forces one to embed the description of this brane in some theory that is more complete than just $`p+1`$ dimensional Einstein-Hilbert gravity in order to describe the brane. Thus we have once again failed to find a stringy negative tension brane which couples only to Einstein-Hilbert gravity.

## 5 Discussion

We have addressed a number of stringy negative tension branes and, in each case, we have identified a plausible mechanism that could make these objects compatible with the generalized second law. We have provided various pieces of evidence, including scaling analyses, perturbative treatments, finding non-perturbative corrections to moduli spaces, and analyses of black hole event horizons.

The two cases of the O8-plane and the AdS soliton were exceptional. For the O8-plane we argued that it is only the required proximity to positive tension D8-branes and the distortions these objects induce in nearby black holes that allows the second law to hold. It would clearly be of interest to compute such distortions, perhaps in the perturbative regime of a small black hole, in order to verify our conjecture that the distortions increase the surface area. On the other hand, we argued for the AdS soliton that it was simply impossible for this brane to meet a ten-dimensional event horizon except at the initial events where the horizon first forms and null generators are added.

In most of the orientifold cases, it was possible to identify a set of corrections that should significantly modify the naïve picture of an infinitely thin negative energy brane interacting with simple (Einstein-Hilbert like) gravitational degrees of freedom. In the O6 case we have seen that these corrections have the interpretation of resulting from massive fields for which the orientifold acts as a source, such that the massive fields have order-one excitations near the orientifold. We believe that this picture holds in general, and the instanton corrections identified for the O5-plane are consistent with this conjecture. These massive fields should have curvature couplings to the naïve Einstein metric, so that in terms of this metric the entropy of a black hole will not be given exactly by $`\frac{1}{4}`$ of the area in Planck units. Such corrections to the entropy should be relevant when an orientifold encounters the horizon of the black hole.

If the black hole completely engulfs the orientifold then these corrections should die out as the system equilibrates. However, for such cases the orientifold has co-dimension greater than 1 and, as described in the introduction, energy considerations suggest that second law violations are of concern only during the collision and not in the final equilibrium state. It is therefore fitting that our resolution should involve corrections of a similar transitory nature.

When we spoke of corrections above, we generally mean ways in which a spacetime description of an orientifold would differ from that obtained by analytic continuation from the positive tension D-branes to appropriate negative tensions. Although we argued that such corrections should exist, we did not map out these corrections in detail for any case other than the O6-plane. Some progress was made for the O5-plane, where we saw that the corrections were related to the Kosterlitz-Thouless unbinding of instanton/anti-instanton pairs.

It would be interesting to find the explicit form of these corrections in other cases as well. Recall, however, that our interest is generally in the region close to the orientifold where we expect the corrections to dominate over the negative tension D-brane geometry. In contrast, natural methods to compute these effects do so in the regime where the corrections are small. While such calculations are therefore unlikely to produce fully satisfactory results, it may be that they reveal interesting properties as we saw in the perturbative calculations of sections 2.3 and 2.4.

The discussion of energy and entropy above was carried out using the naïve Einstein frame. The reader may wonder if one can do better by correcting the metric to obtain a ‘true Einstein frame.’ For definiteness, let us return to the case of the O6-plane and the Atiyah-Hitchin metric in 11-dimensions. In ten dimensions, one may try to find a true Einstein frame by correcting the definition (11) of the ten-dimensional Einstein metric at each order in $`h_{ab}`$. That this is in fact possible may be seen from the fact that the curvature couplings result from the expansion of the 11-dimensional Ricci tensor and are thus linear the ten-dimensional Ricci tensor – the terms do not involve higher powers of the Ricci tensor or other components of the ten-dimensional Riemann tensor. Such terms can always be cancelled by the change in the 10-dimensional Einstein term (essentially the Einstein tensor contracted with the modification of $`g_{ab}^{E,10}`$) induced by an appropriate modification of the metric.

As a result, a true Einstein frame will exist at each perturbative order. Since black hole entropy will be given by the area in this true Einstein frame, this area must match the area as measured by the eleven-dimensional metric. Thus, to the extent that one works in such a true Einstein frame, one expects the description of the physics to match that given in eleven dimensions. One should therefore not find any localized negative energy at all and, to the extent that such a frame makes sense non-perturbatively, the geometry should close off smoothly at the origin with no singular brane to mark any particular points as special. One expects that from this perspective the negative energy is again purely a consequence of the boundary conditions. However, in the weak coupling ($`g0`$) limit all of this structure shrinks to zero size and one is left with an object effectively described as a negative tension brane with orbifold boundary conditions so long as one does not probe this brane too closely.

Note that the ‘true Einstein metric’ described above involves much more than just a conformal rescaling of the metric. In contrast to the more familiar relation between the string and Einstein frames in massless type II supergravity, the present change to the true Einstein frame field will therefore introduce a new notion of null rays, causal structure, and black hole horizons.

Interestingly, in no case did we find a resolution in a regime where the negative tension brane could be described in terms of couplings to pure Einstein-Hilbert gravity. For most of the orientifolds, extra massive fields with complicated curvature couplings seemed to be required for thermodynamic stability. While the O8-brane did not require massive fields in quite the same sense, its resolution resulted from the couplings of massive type IIA supergravity which force the solution to become singular at a finite distance from an isolated brane. When the AdS soliton was viewed as a negative tension brane, the metric necessarily degenerated at the brane and thus passed out of the regime of pure gravity. These are in sharp contrast to the rather benign looking phenomenological negative tension branes used in e.g. , for which the interaction of the metric and brane can be described in terms of the Israel junction conditions. It appears that string theory gives no motivation for stability of simple phenomenological negative tension branes coupled to Einstein-Hilbert like gravity.

While we have by no means ruled out further mechanisms to enforce the second law, it is far from clear what these might be. As a result, one suspects that generic phenomenological negative tension branes face difficulties with the second law and, even when placed at an orbifold, are likely to experience instabilities involving black holes of the sort discussed in . While some more complete treatment will therefore eventually be required, the simplicity of such phenomenological negative tension branes will no doubt maintain their important role in exploring novel and interesting new applications of negative energy objects.

###### Acknowledgments.

We would like to thank Mark Bowick, Clifford Johnson, Dan Kabat, Mark Trodden, and Arkady Tseytlin for useful discussions and to acknowledge the support of the Aspen Center for Physics, where the initial stages of this work were completed. D.M. was supported in part by NSF grant PHY00-98747 and by funds from Syracuse University. S.F.R. is supported by an Advanced Fellowship from EPSRC. S.F.R. also acknowledges the support of the Isaac Newton Institute for Mathematical Sciences during the completion of this work.


text:
# Chandra observations of Abell 2199

## 1 Introduction

Abell 2199 is a nearby rich cluster of galaxies, the central dominant galaxy of which is the cD galaxy, NGC 6166. The intracluster gas is of intermediate temperature (4.7 keV; Edge & Stewart 1991) and gives rise to X-ray emission which is strongly peaked on the central galaxy, leading to the suggestion of a strong cluster cooling flow, e.g., Peres et al. (1998).

NGC 6166 hosts the radio source 3C338 which has been known for nearly 20 years to have a peculiar morphology (Burns, Schwendeman & White, 1983). Ge & Owen (1994) have shown that the radio emission from 3C338 suffers considerable depolarization and Faraday rotation which is consistent with the presence of a dense intracluster medium. Owen & Eilek (1998) have compared the ROSAT HRI X-ray image with a 5GHz radio map and found evidence for a strong interaction between the radio source and the X-ray emitting gas (with the X-ray gas disrupting the radio source). They argue that the intracluster medium is being heated by the radio source.

## 2 Observations

Two observations of Abell 2199 were made using the Advanced CCD Imaging Spectrometer (ACIS-S) on the Chandra satellite. Details are given in Table 1. The standard reprocessed level-two events files have been used for analysis after screening for background flares. Throughout, we adopt a redshift for Abell 2199 of $`z=`$0.0309 and a cosmology with $`H_0=50\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$ and q0=0.5, giving an angular diameter distance of 176 Mpc and a luminosity distance of 187 Mpc. The linear scale is 850pc per arcsec.

In Fig. 1 we show the Chandra X-ray image of the central $``$1 arcmin radius region of Abell 2199 in the 0.3-7.0 keV band. Events from both observations have been combined, taking into account the different aspect solution of the two data sets. The X-ray image was binned to have 1 arcsec pixels. A faint X-ray point source coincident with the active nucleus of NGC 6166 is seen at RA=16:28:38.24, Dec=39:33:04.3 (we use equinox 2000.0 coordinates thoughout). This has been the subject of a previous paper (Di Matteo et al, 2001).

In Fig. 2 we show, as the image, the result of adaptively smoothing (using a three-sigma threshold) the X-ray data and then subtracting the gaussian smoothed (with standard deviation 20 arcsec) data. This emphasises deviations from the large-scale structure. Overlaid on this image is the 1.7 GHz radio map of 3C338 (Giovannini et al., 1998). Depressions in X-ray surface brightness are seen to the east, south and west of the galaxy nucleus. There is a clear correspondence between the position of the eastern radio lobe and a marked depression in the X-ray surface brightness. The south eastern and north western ends of the western radio lobe are also associated with depressions in X-ray surface brightness.

For the case of the eastern radio lobe it may be that the radio-emitting plasma displaces the hot gas as is seen in the Perseus (Fabian et al. 2000) and other clusters (e.g. McNamara et al. 2000). The other depressions to the south may indicate the positions of older radio lobes.

The X-ray data do not account simply for the fact that the radio ridge which joins the two radio lobes is displaced to the south of the nucleus. It is plausible that motion of the nucleus is responsible for this, as discussed by Burns, Schwendeman & White (1983). Probably, it is not the nucleus alone which moves but much of the whole central galaxy and innermost cluster core. If, as seen in projection from Earth, the galaxy has oscillated, or orbited, to the south a distance of about 10 to 20 arcsec (8.5–17 kpc) and has now returned to the north with the radio source off for most of the return trip, then the radio structure may be explained. The bright X-ray emitting region which extends to the south of the nucleus could be a cooling wake, similar to that seen in Abell 1795 (Fabian et al., 2001a).

## 3 Spatially resolved spectroscopy

### 3.1 Analysis of data in projection

The two observations of Abell 2199 were taken with the focal plane at different temperatures, making it inappropriate to combine the event data directly for spectral analysis. We have therefore produced separate spectral files, response data and background files for the two observations, and fitted models to the two observations simultaneously using xspec (Arnaud, 1996). Response matrices and effective area files were made by averaging a $`32\times 32`$ grid of calibrations covering chip 7, using the CIAO tools mkrmf and mkarf, and using weighting factors equal to the number of counts (in the 0.5-7.0 keV band) in the source in the region covered by the calibration. Background spectra were generated from corresponding regions of the standard background fields using software and data from Maxim Markevitch at http://asc.harvard.edu/cal/Links/Acis/acis/Cal\_prods/bkgrnd/current/index.html.

Initially we did not wish to impose a particular geometry on the possible spatial variation of spectral parameters. Therefore, we have used a technique based on the adaptive binning algorithm of Sanders & Fabian (2001) to define tiled square regions of the image where there are at least 1000 counts by progressively increasing the binning of a spatial map until that threshold of counts is reached. Pulse invariant spectra were then extracted from these regions and fitted over the nominal 0.5-8.0 keV band with single temperature mekal (Kaastra & Mewe, 1993; Liedahl et al., 1995) plasma emission models that were allowed to have a freely fitting value for the foreground absorption. (Discrete point sources visible in the Chandra image at positions listed in Table 2 have been excluded from the analysis). We present the results of these fits in Fig. 3.

The temperature map shows significant variations between adjacent pixels. There is a trend in that the temperature drops from values of 4–4.5 keV around the outside edge of our map (R$`75\mathrm{kpc}`$) to values between 1.6–2.2 keV near the nuclear point source. The significance of pixel-to-pixel variations in the metallicity map is much less marked at the edge of our region (where values are between 0.15–0.3 of the solar value), due to the much larger percentage error on this quantity. However, within the inner 45 kpc of the cluster there are significant variations, with individual pixel metalicities varying between 0.25 and 1.2 times the solar value. The absorption map shows that the equivalent neutral hydrogen column density is at least $`2\times 10^{20}\mathrm{cm}^2`$ over most of the mapped region. Only a few pixels are consistent with the value expected from our Galaxy for which $`N_H=8.7\pm 0.6\times 10^{19}\mathrm{cm}^2`$ (unweighted average of surrounding grid points; Dickey & Lockman 1990) in the direction of Abell 2199. Near the centre of the cluster the absorption rises to between $`510\times 10^{20}\mathrm{cm}^2`$. Finally the map of reduced chi-square shows that most of the pixels have values of chi-square between 0.8 and 1.2. Only a few pixels, predominantly near the centre of the cluster have higher values.

We have overlaid the radio contour map on the temperature map, to facilitate comparison. Although no detailed correlation is evident between structure in the temperature or metallicity maps and the radio map, the cooler regions of the cluster are associated with the radio brighter parts of the radio source.

In Fig 4 we have plotted the values of absorption column density, temperature and metallicity as a function of the radial distance of the bin from the cluster centre. Inspection of this plot confirms that most bins at similar radii have similar properties although there are a small number of values inconsistent with the mean for that radius. In particular, we have looked at the region near the two cool pixels centred on (RA=16:28:36.9, Dec=39:32:55.8). We find that a single temperature fit gives kT$`1.3\mathrm{keV}`$, rather cooler then the surrounding regions.

#### 3.1.1 Single temperature radially varying models

In the absence of evidence for gross non-radial variations in the properties of the gas we have made spectra binned azimuthally into eight concentric rings out to a radial distance of 200 kpc from the central point source. One further irregular shaped bin is included which contains all of the data that are outside 200 kpc radius region but which still lie on chip 7. Bins beyond a radial distance of 140 kpc from the nucleus cover slightly different spatial regions in the two data sets, because the cluster is not centred on the chip and the spacecraft roll angle was different for the two observations.

We have fitted a single temperature mekal model absorbed by a foreground screen (to account for absorption from the Galaxy) to each of the concentric spatial regions. Fig. 5 shows the radial profiles of temperature, metallicity, equivalent hydrogen absorbing column density and reduced chi-square. The absorption, parameterised as the equivalent column density of neutral hydrogen, is greater than that expected from 21cm observations of the Galaxy thoughout the region studied. (This may be related to uncertainties in the calibration of the instrumental effective area below 1keV111 http://cxc.harvard.edu/cal/Links/Acis/acis/Cal\_projects/index.html ). Even with this caveat, though, we see an increase in the absorption within the central 20 kpc, relative to that seen further out, by $`2\times 10^{20}\mathrm{cm}^2`$.

The cluster shows a smoothly declining temperature profile that starts at $`4.2`$ keV beyond 200 kpc, reducing to near 2 keV within the central 5 kpc. The broad-beam temperature for this cluster was measured with EXOSAT as $`4.7\pm 0.4`$ keV by Edge & Stewart (1991) indicating that the temperature profile climbs further beyond the spatial extent of our analysis. Our outer measured temperature is in good agreement with the value of $`4.16\pm 0.15`$keV measured from the central 4-5 arcmin radius region by ASCA (Allen et al., 2001a).

Within the central 5kpc we find that the metallicity is close to the ASCA measured value for the central 4-5 arcmin radius region (Allen et al., 2001a) of 0.4 times the solar value. (We use the default solar abundance set given by Anders & Grevesse 1989). By a radial distance of 15kpc, the metallicity has risen to 0.6 times the solar value, before reducing again at larger radii to near 0.25 times the solar value.

#### 3.1.2 Multi-temperature models

We have also examined the annular spectra using multi-temperature models. We first fitted the data with a two-temperature model consisting of freely fitting absorption acting on two single-temperature emission components; the abundances were tied together but the temperatures were allowed to fit freely. The introduction of the second emission component is required (at the 99 per cent confidence level, determined using an F-test) for the three radial bins 10–20kpc, 20–40kpc and 140–200kpc. However, in section 3.2.1 we show that the evidence for multi-temperature components in these annular regions is largely due to projection effects.

Recently, Protassov et al. (2002) have questioned the validity of using the F-test to test for the significance of an additional additive component when fitting models to data because the hypothesis is on the boundary of the parameter space. In our case we have used the F-test to check whether a second mekal component is statistically required by the data. This is on the boundary of parameter space since the null hypothesis is that the second mekal component is not required. In this case the normalization of the second component is zero, and negative normalizations (although possible in practice within the fitting code) are excluded from the fitting procedure as they are not physically meaningful. We have carried out fits to Monte Carlo realizations of our data and find that (in our case) the distribution function of our measured F values matches the formal F distribution very well for values of F above 0.2222 To check the validity of the F-test we have repeatedly run fits with single and two-temperature mekal models to Monte Carlo realizations of our data (generated using the xspec fakeit command) and looked at the distribution of the F-statistic. One thousand realizations of our data were generated from the single temperature plasma model (with freely fit absorption) which best fitted one combined data set, by folding the model through the detector response function and adding random noise appropriate to the total counts generated. First we fitted the same model that was used to generate the data and noted the chi-square value and number of degrees of freedom. We then fitted a two-temperature mekal plasma model in which the metallicity of the second component was fixed to that of the first and the absorption component acted on both mekal components. The temperature and normalization of the second component are then the new fit parameters; we logged the value of the chi-square statistic and number of degrees of freedom in this case. A Kolmogorov-Smirnov test was used to determine whether the distribution of measured F values was significantly different from the true F distribution. For all values of F above 0.2 we find the two distributions to be in excellent agreement (the probability of differences between the two cumulative distribution functions as large as observed (if they are the same distribution) is greater than 20 per cent). Below a measured F-value of 0.2 the two distributions are significantly different (much less than 1 per cent chance of deviations as large as observed if the two distributions are the same). This may be due to either to xspec sometimes not finding the exact minimum value of chi-square, or to problems with the test discussed by Protassov et al. (2002). Since the critical values of F used in our assessment of the significance of extra spectral components in this paper lie near F=4.6, we conclude that the boundary condition problem does not affect our inferences. .

We next explored the effect of replacing the second temperature component with a cooling flow spectrum (in which gas cools to zero degrees from the plasma temperature). We obtain similar results to those found with the two-temperature model and place limits on the size of the cooling flow component within R=100kpc of $`12\pm 3\mathrm{M}_{}\mathrm{yr}^1`$ or $`37_{13}^{+14}\mathrm{M}_{}\mathrm{yr}^1`$ (90 per cent confidence limits), depending whether or not freely fitting intrinsic absorption is included on the cooling flow component.

We have also investigated the region, mostly to the south of the nucleus, defined by the polygon with the following vertices: (16:28:40.03,39:32:54.7), (16:28:40.03,39:32:54.7), (16:28:38.17,39:33:12.9), (16:28:36.25,39:32:54.2), (16:28:40.05,39:32:54.1), but excluding a region of radius 2.5 arcsec centred on the nucleus. This region was targeted for separate analysis because it stands out as having a higher surface brightness than surrounding regions. We find that the gas temperature (fitted with a single-temperature model) is cool (T$`=2.5\pm 0.1`$ keV) and that the data can accommodate a cooling flow of up to $`9\mathrm{M}_{}\mathrm{yr}^1`$. Since this triangular region contains regions of bright radio emission we have checked for the presence of power-law continuum X-ray emission. No significant such emission was found: a two temperature fit is significantly better than a single temperature plus power-law continuum.

### 3.2 Analysis of data accounting for projection effects

#### 3.2.1 Single temperature fits

Assuming that the cluster is spherically symmetric and that the physical conditions change only in the radial coordinate of the cluster, we can account for the projection effects of gas exterior to a particular region. This allows us to determine the physical conditions in the gas as a function of the true 3-dimensional radial coordinate. We do this by fitting spectra from concentric rings simultaneously using multiple mekal models. We allow one mekal model for the outermost annulus, two for the next one in and so on, setting the normalizations of the different projected components according to the relative volumes of the shells. The technique is described in full in Allen et al. (2001b). Another approach, used on the cluster Abell 1795, is described by Ettori et al. (2002).

The temperature and abundance profiles obtained from this method are shown in Fig. 6. The Galactic absorption was allowed to be a free parameter and its value was determined as $`3.05\pm 0.09\times 10^{20}\mathrm{cm}^2`$. The temperature profile shows the central bin to be cooler than that measured in the projected data: $`1.57_{0.16}^{+0.20}`$ keV compared with $`2.08_{0.10}^{+0.12}`$ keV. The profiles both from the simple analysis presented in section 3.1.1 and from this analysis asymptote to similar values, $`4.26\pm 0.07`$keV (deprojected) and $`4.21\pm 0.07`$keV (seen in projection) at radii of 140-200kpc.

The deprojected abundance profile is very similar to that seen in the simple single-temperature fits to the data in projection. The low value seen in the central bin, and the value seen at the abundance maximum near 20 kpc, are more extreme in the deprojected analysis but with larger uncertainties.

#### 3.2.2 Two temperature fits

Finally, we have investigated whether there is a statistical requirement for a second emission component in any annulus, once the projection of the first emission components have been taken into account. To do this, we use the xspec projct model and add in second emission components to each annulus successively (starting at the outside, allowing only the temperature and normalizations to be freely fit), and note how the value of the chi-square statistic reduces. In each case the metal abundance of the second component is tied to the metal abundance of the corresponding first component, and a single absorption component is applied to the whole model. The fiducial fit, in which each region has only one emission component, gave a value $`\chi ^2=4865.3`$ for 3761 degrees of freedom. Table 3 shows the values of the chi-square statistic, number of degrees of freedom and corresponding values of the F-statistic, together with the formal probability that such a large value of F would be obtained if the new emission component was not required. We caution that due to the large value of the reduced chi-square statistic the true probabilities will be larger than listed in the table, so the results are less significant than tabulated.

We find that the outermost region strongly requires a second emission component. This is not surprising since by definition we have not accounted for the projection of emission external on to this region. There is also evidence, but at a lower significance, for a second emission component being required between 20-100kpc.

## 4 Image deprojection and Mass Analysis

The deprojected temperature and surface brightness profiles can be used to constrain the mass profile of the cluster. For this analysis we have used the methods described by Allen et al. (2001b) and Schmidt, Allen & Fabian (2001). In brief, the observed surface brightness profile and a particular parameterized mass model are together used to predict the temperature profile of the X-ray gas, which is then compared with the observed, spectrally-determined results.

Three separate parameterizations for the cluster mass distribution were examined: the Navarro et al. (1997) (NFW) profile, the Moore et al. (1998) profile and a non-singular isothermal sphere (NIS). We find that the NFW profile provides (marginally) the best fit to the data, although all three models give formally unacceptable results. Table 4 lists the best-fit parameters for the three models, together with their $`1\sigma `$ uncertainties, chi-square values, and the number of degrees of freedom for the fit. For the NFW and Moore profiles we note that there are only two fit parameters; the scale radius, $`r_\mathrm{s}`$, and the concentration parameter, $`c`$. However, we also quote values for the effective velocity dispersion $`\sigma =\sqrt{50}r_\mathrm{s}cH`$, which is a function of the other two variables and the Hubble constant.

The temperature profile predicted by the best-fit NFW mass model overplotted on the observed, deprojected temperature profile is shown in the top panel of Fig. 6. The total gravitating mass profile associated with the best-fit NFW model is plotted in Fig. 7, along with its $`1\sigma `$ uncertainties.

The shape of the inner gravitational well around NGC6166 has been probed by the optical spectroscopy of Carter, Bridges & Hau (1999), who find that the velocity dispersion of the galaxy rises outward. Recently Kelson et al. (2002) have confirmed and extended this result. They found that the velocity dispersion continues to rise out to the limit of their data at 60kpc, by which time the stellar velocity dispersion is close to that of the cluster galaxies. Kelson et al. (2002) note that our measurement of the mass within $`0.1h^1`$Mpc (where $`h=H_0/100`$) is in reasonable agreement with their determination, for the same region.

We next used the best-fitting gravitational potential that was determined above to constrain an image deprojection analysis of the Abell 2199 cluster. We follow the method of Fabian et al. (1981) in determining the physical conditions in the cluster as a function of the radius. The Galactic value of the column density is used. Error bars for the derived quantities were determined from one hundred Monte Carlo replications of the original data.

In Fig. 8 we show the major results from the deprojection analysis. The input surface brightness distribution is shown in the top panel together with the implied electron density, cooling time and mass deposition rate in lower panels.

The nominal mass deposition rate integrated out to the point where the cooling time is equal to the age of the universe is $`150\mathrm{M}_{}\mathrm{yr}^1`$. However, the spectral fits (section 3.1.2) place a limit on the mass deposition rate of less than $`50\mathrm{M}_{}\mathrm{yr}^1`$ suggesting that any steady state cooling flow must be younger than 2–3 Gyr (associating the age of the cooling flow with the cooling time of gas where the mass deposition rate from the image deprojection equals that in the spectral fits).

## 5 Heating and cooling rates

We have seen that the radiative cooling time of the gas in Abell 2199 is less than about 7 Gyr within a radius of 100 kpc, less than 2 Gyr within 20 kpc and about 0.1 Gyr within the innermost few kpc. The temperature decreases from above 4 keV beyond 100 kpc to 1.6 keV (deprojected) within 5 kpc. This overall appearance suggests that a cooling flow may be operating (see e.g. Fabian 1994), but the spectral analysis provides only modest support and then only for a small flow of tens $`\mathrm{M}_{}\mathrm{yr}^1`$. XMM-Newton Reflection Grating Spectrometer (RGS) data of other clusters also show little evidence for long-lived, steady-state cooling flows, in which gas cools from the cluster virial temperature to below $`10^6\mathrm{K}`$ by emitting observed X-radiation (Peterson et al. 2001; Tamura et al. 2001; Kaastra et al. 2001). The RGS spectra, which do not spatially resolve any flow, are consistent with the gas cooling from the virial temperature to about one third to one fifth of that value but no lower, possibly in a single-phase cooling flow. The fate of the coolest gas, at 1–2 keV, is unknown. It may be heated, mixed with much cooler gas or expelled (see e.g. Fabian et al. 2001b).

Here we attempt to quantify upper limits to the amount of multiphase gas in the core of Abell 2199, and to see where the steady-state cooling flow assumption breaks down. Using a similar (projct) spectral model, accounting for the projection of outer components on to inner regions, to that discussed in section 3.2.2, limits on the volume filling fraction ($`f_2`$) of any second emission component in the central four spatial bins (out to 40 kpc from the centre of the cluster) were calculated using

$$f_2=\left[1+\frac{e_1}{e_2}\left(\frac{T_1}{T_2}\right)^2\right]^1,$$
$`(1)`$
Sanders & Fabian (2002). $`e_1`$ and $`T_1`$ are the emission measure and temperature of the first temperature component for a particular shell, and $`e_2`$ and $`T_2`$ are the emission measure and temperature of the additional component. Unlike section 3.2.2 we only applied the second emission component to one spatial region at a time. The second component has a freely fitting normalization and a metal abundance that is fixed to the metal abundance of the first component. We stepped the temperature of the second component through a range of values from 0.5keV to 2.5keV, calculating $`f`$ and its uncertainty $`\sigma _f`$ (obtained by propagating the uncertainties in $`e1`$, $`e2`$, and $`T_1`$). Where the value of $`f`$ is less than $`3\sigma _f`$ we plot an upper limit at $`3\sigma _f`$. Our results are shown in Fig. 9. Also shown, in Fig. 10, for comparison, is the continuous differential filling factor distribution expected from a constant pressure cooling flow.

It is clear that there is little evidence for the gas being multiphase at any radius. (Note that the ‘detections’ of a second component at some temperatures in the 10-20kpc region may be in large part due to the projection of the second component present in the 20-40kpc region). For a continuous distribution of temperatures the y-axis is $`df/dT`$ (Fig. 10) and limiting curves lie below the plotted temperature limits which represent the step by step values for one additional temperature component only. (Note also that the test applied in section 3.2.2 is different from the test applied here, as in the former each additional component added to outer rings was left in when evaluating the F-statistic for the interior ring, whereas here we have added the second component in for each region separately and then removed it before adding in the second component for the next region).

Since the temperature distribution appears to be close to single phase, we have then tested whether it resembles a single-phase cooling flow dropping to some minimum temperature $`T_{\mathrm{min}}1\mathrm{keV}`$. Although the temperature profile (Fig. 6) appears consistent with a single-phase cooling flow, the surface brightness profile is not. This is the reason that multiphase flows were introduced in the first place (e.g. Fabian, Nulsen & Canizares 1984). Fig. 11 (solid points) shows the distribution of emission measure against radius. We use the normalizations from xspec, which need to be multiplied by $`3.8\times 10^{68}`$cm2 to convert to $`\mathrm{cm}^3`$.

Also plotted is the distribution expected from a single phase cooling flow of $`100\mathrm{M}_{}\mathrm{yr}^1`$, both with and without the gravitational work done included. To obtain this we use the energy equation of a steady state spherical inflow,

$$\dot{M}\left(\frac{5}{2}\frac{k}{\mu m}\frac{dT}{dr}+\frac{d\varphi }{dr}\right)𝑑r=e\mathrm{\Lambda }+h.$$
$`(2)`$
$`\mu m`$ is the mean mass per particle, $`h`$ the heating rate (assumed zero for Fig. 11), $`e`$ the emission measure and $`\varphi `$ the gravitational potential, obtained from Fig. 7.

We see that the observed emission measure distribution is nowhere as flat as expected from a single-phase flow (the expected values scale up and down with $`\dot{M}`$). There is too little emission at small radii and too much at large radii. The latter can be explained by any flow not yet having reached a steady-state, but the lack of emission at the centre is more difficult to understand. We have already shown that spectral fits to the data exclude any simple distributed mass dropout such as expected from a multiphase flow with a mass deposition rate greater than $`15\mathrm{M}_{}\mathrm{yr}^1`$. However, a flow can be accommodated if the gas, or its emission, can be made to be unobservable. Including absorption is one way to achieve this, and the spectral data allow up to $`50\mathrm{M}_{}\mathrm{yr}^1`$ if the emission is absorbed. The missing soft X-ray luminosity in cooling flows is discussed more extensively by Fabian et al. (2002).

Finally, we have quantified the heating required to maintain the gas in a steady state (Fig. 12), by estimating $`h`$ from equation (2). Here a negative value for the heating means that heat has been lost in a non-radiative manner (for example by mixing with cold gas). Heating must be distributed over a wide range of radii, with a heating rate per unit volume proportional to $`r^{1.5}`$ if there is no mass dropout. The results are only relevant where some steady-state can be assumed. Note that the radio source, at present, affects the gas distribution in the east-west axis out to about 30 kpc.

The peaked abundance distribution determined from the deprojected spectral analysis, Fig. 6, argues against much mixing, or convection of gas across the annulus at 20–30 kpc. The whole of the current gas distribution is stable against convection since the entropy, $`ST/n^{2/3}`$, monotonically increases with radius ($`Sr^{0.8}`$ from a few kpc to about 100 kpc since the deprojected gas temperature is approximately $`1.11r^{0.29}\mathrm{keV}`$ and the electron density is $`0.21r^{0.75}\mathrm{cm}^3,`$ where the radius $`r`$ is in kpc).

## 6 Conclusion

The core of Abell 2199 is complex. Some depressions in the X-ray surface brightness are seen to coincide with the radio lobes. Other depressions to the south may indicate earlier radio activity. The temperature decreases from 4.2 keV to 1.6 keV over radii from 100 kpc to 5 kpc, which is where the radiative cooling time drops from 7 Gyr to 0.1 Gyr. The X-ray surface brightness is not sufficiently peaked however to be consistent with a single phase cooling flow. We find no strong evidence for multiphase gas in the cluster core except in the 20-40kpc region where the radio lobes have disturbed the gas most. We limit the mass deposition rate within the central 100 kpc to $`12\mathrm{M}_{}\mathrm{yr}^1`$. Only if the emission below about 1 keV is somehow undetectable, such as occurs if there is intrinsic absorption, can this rate be significantly increased to around $`40\mathrm{M}_{}\mathrm{yr}^1`$.

We have quantified the heating rate and its distribution in the cluster required to balance cooling if the gas is in a steady state. However, the situation may not be steady. Any cooling flow is only likely to occur and have become steady within the region where the cooling time is a few billion year or less. The radio source is clearly causing disruption to the south, east and west of the nucleus. How much disruption is perhaps indicated by the abundance profile. Our spectral fitting shows that the evidence for two temperatures is almost entirely due to projection effects and Fig. 6 gives the best description of the state of the gas. The abundance shows a maximum at around 30 kpc. The profile (which has a similar shape to that seen in the Centaurus Cluster; Sanders & Fabian 2002) is consistent with some mixing of gas within, but not much beyond, that radius. Abundance gradients can be due to metal injection by SN Ia in the central galaxy (Dupke & White, 2000; Ettori et al., 2002) and are expected to show a monotonic increase inward (unless the metals are poorly mixed and cooling dominates; Morris & Fabian 2002, in preparation). A central decrease can be explained by some exchange of gas between the centre and the region at 40–50 kpc (or further out). Why intermediate radii have not been similarly affected is puzzling. Massive mixing should however have wiped out the gradient. The high abundance of the gas at 10–40 kpc indicates that it has probably resided close to the centre of the galaxy for an extended period.

## Acknowledgments

We thank Stefano Ettori for the use of his software to make exposure maps and to extract surface brightness profiles. SWA and ACF acknowledge support by the Royal Society.