Nuclear Physics B 744 (2006) 1–13
Symmetry and entropy of black hole horizons
Olaf Dreyer, Fotini Markopoulou, Lee Smolin ∗
Perimeter Institute for Theoretical Physics, 35 King Street North, Waterloo, ON N2J 2W9, Canada Department of Physics, University of Waterloo, Waterloo, ON N2L 3G1, Canada Received 31 August 2005; accepted 3 February 2006 Available online 31 March 2006
Abstract We argue, using methods taken from the theory of noiseless subsystems in quantum information theory, that the quantum states associated with a Schwarzschild black hole live in the restricted subspace of the Hilbert space of horizon boundary states in which all punctures are equal. Consequently, one value of the Immirzi parameter matches both the Hawking value for the entropy and the quasi normal mode spectrum of the Schwarzschild black hole. The method of noiseless subsystems allows us to understand, in this example and more generally, how symmetries, which take physical states to physical states, can emerge from a diffeomorphism invariant formulation of quantum gravity. © 2006 Published by Elsevier B.V.
1. Introduction
This paper is concerned with two distinct problems in background independent formulations of quantum gravity. First, how one can recover the semiclassical Bekenstein–Hawking results [1,2] in the context of results about the entropy of horizons and surfaces in loop quantum gravity [3–6]. Second, how the description of spacetime in terms of classical general relativity emerges from the quantum theory. Up till now, these problems have been treated in isolation. Here we propose a perspective that relates them, and also addresses a technical issue concerning the black hole entropy. This technical issue concerns the value of a certain free parameter, γ , in the theory called the Immirzi parameter [7]. This parameter, not present in the classical theory, arises from an
* Corresponding author. E-mail addresses: [email protected] (O. Dreyer), [email protected] (F. Markopoulou), [email protected] (L. Smolin).
0550-3213/$ – see front matter © 2006 Published by Elsevier B.V. doi:10.1016/j.nuclphysb.2006.02.045 2 O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 ambiguity in the quantization procedure because there is a choice in the connection variable that is used in the construction of the quantum observable algebra. As it does not appear in the classical theory, its value should be fixed by some physical requirement. In previous work [3–5], it was shown that γ is fixed by the requirement that the quantum gravity computation reproduce the Hawking value of the entropy. Since the story is very simple and physical, we review it for non-experts. The Immirzi para- meter comes into the formula for the area of a surface. When a surface is punctured by a set of spin network edges with spin labels {jα}, then the area of the surface is given by A {jα} = 8πγ jα(jα + 1), (1) α 2 which we take here to be measured in units of the Planck area lPl. In loop quantum gravity, there is an exact description of the quantum geometry of black hole horizons [5], as well as a more general class of boundaries [3]. The entropy of the horizon is defined in terms of the Hilbert space of the boundary theory. For reasons that we spell out shortly, it is taken to be the logarithm of the dimension of the Hilbert space of boundary states. The Immirzi parameter does not come into the entropy. Hence it appears in the ratio of the entropy S[A] of a black hole to its area: S[A] c R = = , (2) A[{jα}] γ where c is a constant to be determined by calculating the entropy. At the same time, we believe from Hawking’s semiclassical calculation of black hole entropy that [2] 1 R = . (3) 4 For the two equations to agree, γ = 4c. (4) Recently, one of us, following a clue uncovered by Hod [8], described a semiclassical argu- ment which fixes γ by an argument that appears to be independent of considerations of entropy or radiation but rather makes use of the quasinormal mode spectrum [6]. As it is relevant, we repeat here the basic argument. The quasinormal mode spectrum of the Schwarzschild black hole turns out to have an asymp- totic form ln(3) i 1 Mω = + n + , (5) 8π 4 2 for arbitrarily large integer n.Asn goes to infinity, the decay times of the modes goes to zero in units of the light crossing time of the horizon. This means that the excitations involve more and more local regions of the horizon. At the same time, the real part of the frequency goes to the asymptotic value [9] ln(3) ω = . (6) qnm 8πM This frequency is then associated with an arbitrarily short lived, and hence local, excitation of the horizon. By the correspondence principle, this must correspond to an energy ln(3)h¯ M = hω¯ = . (7) qnm qnm 8πM O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 3
In the classical theory, the region of the horizon excited may be arbitrarily small, but in the quantum theory there is a smallest region that can be excited, which is a minimal puncture with minimal area A = 8πγ jmin(jmin + 1). (8) However, for a Schwarzschild black hole,
A = 16πM2, (9) from which it follows that, for a quasinormal mode, = = 2 Aqnm 32πMM 4ln(3)lPl. (10) The result is a prediction for the Immirzi parameter, which is ln(3) γ qnm = √ . (11) 2π jmin(jmin + 1) Does this match the value needed to get R = 1/4? It turns out to depend on what we take for the Hilbert space of the black hole horizon. By arguments given in [5], we know that the horizon Hilbert space, for fixed area A in Planck units, is related to the Hilbert space of the 2 { } V 1 invariant states of U(1) Chern–Simons theory on an S with punctures jα , denoted by {jα}. The level k is related to the mass of the black hole, and is assumed here to be large. Given a set of punctures {j}, V = H {jα} Inv jα (12) α where Hj is the (2j + 1)-dimensional spin space for spin j and Inv means the Hilbert states contains only states with total spin zero. We know that the Schwarzschild black hole has definite area A[{j}] and definite energy:
A({j}) M {j} = . (13) 16π As the area is quantized, we see that the energy is quantized as well. V We then expect that at least some of the Hilbert spaces {jα} contain states which correspond to Schwarzschild black holes. The state is expected to be thermal, because it is entangled with the states of radiation that has left the black hole. We are, however, ignorant of the exact state so, instead, the black hole entropy is usually taken to be Schwarzscild = HSchwarzschild S ln dim A , (14) HSchwarzschild Schwarzscild where A is a Hilbert space on which the density matrix ρA describing the ensemble of Schwarzschild black holes may be expected to be non-degenerate. One proposal is to take HSchwarzschild = HAll = V A A {j}, (15) A({j})=A
1 There is an equivalent description in terms of SU(2) Chern–Simons theory [3]. 4 O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 where we use the superscript All to remind us that in this case we sum over all sets of punctures that give an area between A and A + A, for some small A. However, this proposal can be criticized on the basis that the only information about the Schwarzschild black hole that is used is its area. We may expect that the actual state of the Schwarzschild black hole is determined by additional physical considerations and hence is non- degenerate only in a subspace HSchwarzschild HAll A A (16) that is picked out by additional physical input. In this article, we will argue that the physical Hilbert space associated to the horizon of a Schwarzschild black hole is dominated by the V{j} for which all punctures have equal minimal spin jmin. The difference is important, because in the first case we have ln dim HAll RAll = A , (17) A({j}) from which recent calculations [10] have deduced the value2
γ All = 0.23753295796592 ... . (18) On the other hand, in the second case
ln dim V{j ,...,j } Rmin = min min . (19) A({j}) This is easy to compute [3,6] and leads to a value ln(3) γ min = √ . (20) 2π jmin(jmin + 1) We note that
γ qnm = γ min = γ All. (21) Hence, we have the following situation. If indeed the right choice for the Hilbert space of a Schwarzschild black hole consists only of states with equal and minimal punctures, then there is a remarkable agreement between the two computations. This supports the case that there is something physically correct about the description of black holes in loop quantum gravity. HAll On the other hand, if A is the right Hilbert space for the horizon of a Schwarzschild black hole, loop quantum gravity is in deep trouble, because there is no choice of γ that will agree with both the Hawking entropy and the quasi-normal mode spectrum. It is clear that to resolve this problem we need additional physical input. The only property of the Schwarzschild black hole used in previous work to constrain the corresponding quantum state is its area. But a more complete treatment should take into account other characteristics of the Schwarzschild black hole such as its symmetry and stability. To do so, we need to understand how those classical properties can emerge from a quantum state, described so far in the background independent language of loop quantum gravity.
2 But see [17] for a criticism of that calculation. O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 5
Recently, one of us proposed a new perspective on the general problem of the emergence of particles and other semiclassical states from background independent approaches to quantum gravity [13]. This makes use of the concept of noiseless subsystems, developed in the context of quantum information theory to describe how particle-like states may emerge there [12]. We shall see in the following that the black hole problem provides a nice example of the general strategy proposed in [13], while it resolves the present problem.3
2. Symmetry and noiseless subsystems
We begin by asking two questions:
HSchwarzschild (1) How do we find the subspace A corresponding to a Schwarzschild black hole, as opposed to a general surface of a given area that satisfies the appropriate boundary condi- tions? (2) How do we recognize in that subspace the excitations that, in the classical limit, correspond to the quasinormal modes?
We should here make an important comment, which clarifies the sense in which these ques- tions are asked. We note that the boundary conditions we require for our analysis are more general and apply to all horizons, not just Schwarzschild black holes. Thus, our framework differs from the isolated horizon framework [5] in one crucial aspect. We posit that the states of all quantum black holes live in a single Hilbert space, so that angular momentum and other multiple moments are to be measured by quantum operators. In the isolated horizon picture the angular momentum and multiple moments are treated classically, and the Hilbert space is defined for each value of them. Only in the our framework does it make sense to ask how the states corresponding to a non-rotating black hole emerge from a more general Hilbert space of horizon states. Our two questions are then analogous to problems concerned with the emergence and stability of persistent quantum states in condensed matter physics and quantum information theory. For the present purposes, we will use the following idea from the theory of noiseless subsystems in quantum information theory. We have a complete quantum experiment, which we want to divide into the system S and environment E. We want to understand what quantum properties of the system may survive stably in spite of continual and uncontrollable interactions with the environment. The joint Hilbert space decomposes into the product of system and environment, S E Htotal = H ⊗ H , (22) while the Hamiltonian decomposes into the sum S E H = H + H + H int (23) where H S acts only on the system, H E acts only on the environment and all the interactions int between them are contained in H . The reduced dynamics of S is given by a completely positive operator φ : HS → HS , S = S = S ⊗ E † ρf φ ρi trE U ρi ρi U , (24)
3 While this paper was in draft, a paper was posted that reaches a similar conclusion by a different argument [11].More recently, additional arguments for it have been given by [18]. 6 O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 where the joint state of S and E evolves unitarily and then the environment is traced out. While a generic “noise” φ will affect the entire state space HS , there may be a noiseless subsystem of S, and possibly even a subspace of HS that is left unchanged by φ, i.e., evolves unitarily. What follows is a necessary and sufficient condition for the existence of a noiseless subsystem. Eq. (24) can be rewritten as S = S † φ ρ Akρ Ak, (25) k where
ψi|Ak|ψf =ψi|⊗ek|U|ψf ⊗|ek , (26) | | ∈HS {| } HE for any ψi , ψf and ek an orthonormal basis on . One can check that † = k AkAk 1. If B(HS ) is the algebra of all operators acting on HS ,theinteraction algebra Aint ⊆ B(HS ) int is the subalgebra generated by the Ak (assuming that A is closed under † and φ is unital). int Up to a unitary transformation, A can be written as a direct sum of dj × dj complex matrix algebras, each of which appear with multiplicity μj :
Aint ⊗ B dj 1μj C , (27) j
μj where 1μj is the identity operator on C . The commutant, Aint of Aint is the set of all operators in B(HS ) that commute with every element of Aint, i.e., Aint B μj ⊗ C 1dj . (28) j This decomposition induces a natural decomposition of HS :
S H = Cμj ⊗ Cdj . (29) j
int Note now that any state ρ in A is a fixed point of φ since it commutes with all the Ak: S = S † = † S = S φ ρ Akρ Ak AkAkρ ρ . (30) k k It can be shown [12] that the reverse also holds, i.e., S φ ρ = ρ ⇔ ρ ∈ Aint . (31) Hence, the noiseless subsystem can be identified with the Cμj in Eq. (29). These are relevant for the physical description of the system, because the interactions with the environment will not disturb them and hence are useful for describing the long term behavior of the system, because they are conserved. The relevance of this argument for quantum gravity was proposed in [13]. If we divide the quantum state of the gravitational field arbitrarily into subsystems, those properties which are conserved under interactions between the subsystems are going to characterize the low energy limit of the spacetime geometry. If classical spacetime physics emerges in the low energy limit, in which these regions are arbitrarily large on Planck O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 7 scale, the commutant of the interaction algebra should include the symmetries that characterize classical spacetime in the ground state (presumably the Poincaré or the de Sitter group, or some deformation of them). In the case of vanishing cosmological constant, we then expect that the low energy limit is characterized by representations of the Poincaré group. But these, of course, are the elementary particles, as described in quantum field theory on Minkowski spacetime. The noiseless subsystem idea is then a way to understand how particle states, and the low energy vacuum could emerge from a background independent quantum theory of gravity. In the following section, the system state space is that of the states of loop quantum gravity that correspond to a black hole of a given area. This is coupled to an environment of bulk spin network states and, while we do not know the quantum dynamics, we do know that a Schwarzschild black hole has a classical SO(3) symmetry. In the next section, we make the reasonable assumption that there is a microscopic analogue of SO(3) that acts on the system, a symmetry that should be contained in the commutant of the interaction algebra for black holes.
3. Application to black holes
HSchwarzschild In this section we shall look for the subspace A of states that correspond to a HSchwarzschild classical Schwarzschild black hole. We want to show that A is a proper subspace of HAll ± the space A of states of horizons with area A A, i.e., HSchwarzschild HAll A A . (32) The geometry of a horizon is only partially fixed by its area A. A Schwarzschild horizon is a very special example of such a geometry that is distinguished by its SO(3) symmetry. Since all the HAll other horizon geometries have quantum counterparts in A as well, we expect that the Hilbert HSchwarzschild HAll space A corresponding to a Schwarzschild black hole is a proper subspace of A . In contrast to the smooth horizon geometry of a classical Schwarzschild black hole, the geom- etry of the corresponding quantum horizon is concentrated at the punctures and is thus discrete. HSchwarzschild Consequently, we cannot expect an action of the smooth Lie group SO(3) on A . What we can expect is an action of some discrete symmetry group Gq that only on the classical level coincides with the action of SO(3) (see Fig. 1). What can reasonably be assumed about the group Gq ? Since Gq is a group of symmetries we assume it commutes with the action of the Hamiltonian H int on the horizon Hilbert space HS , i.e., int Gq ⊂ A . (33) We do not know the precise action of the Hamiltonian H int but it is possible to restrict it suffi- ciently to compute the commutant Gq . We first review the basics of the description of black hole horizons in loop quantum gravity [5]. The Hilbert space is at the kinematical level of the form (22), where the environmental Hilbert space has a basis consisting of all spin networks that end on the horizon, in any number of punctures. Given that the energy of a black hole is a function of its area, it makes sense to decompose the Hilbert space in eigenspaces of the area operator. The eigenvalues are given by sets of spins at the punctures {ji}, where i labels the punctures. Gauss’s law acting at the boundary enforces that the labels on the punctures on the horizon match the spins of the edges of spin network in the bulk, incident on the horizon. We then have S E Htotal = H ⊗ H . (34) {ji } {ji } {ji } 8 O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13
Fig. 1. The geometry of a quantum horizon is concentrated at a discrete set of punctures. The classical symmetry group SO(3) will thus not act directly on the quantum states representing a Schwarzschild black hole. A discrete group Gq ,on the other hand, is expected to act on these states. In the classical limit, the action of this discrete group coincides with that of SO(3). For a large number of punctures, the discrete group action shown above would approximate a rotation of the horizon.
The environment (sometimes called “bulk”) Hilbert spaces HE have bases consisting of dif- {ji } feomorphism classes of spin networks in the bulk, with edges ending on the boundary at the punctures with the specified spin labels. We note that the diffeomorphisms on the surface are fixed. (There may also be an exterior boundary at which the diffeomorphisms are also restricted.) The system, or surface, Hilbert spaces HS are direct sums of the one-dimensional Hilbert {ji } spaces of U(1) Chern–Simons theory on the sphere, with punctures labeled by charges mi at 2 points σi ∈ S , subject to the conditions | | mi ji, (35) mi = 0. (36) i
The level k of each Chern–Simons theory depends on the ji ’s and is given by [5] a [{j }] k = 0 i (37) 4πγ where a0[{ji}] is the nearest number to A[{ji}] such that the level k is an integer. We note that the level k is very large for black holes large in Planck units, so that for some estimates the limit k →∞can be taken. This is normally done in the computation of the entropy, for which the difference between the classical and quantum dimensions may be neglected. The U(1) connection on the boundary satisfies 4π F (σ ) = m δ2(σ, σ ). (38) ij k i i i We then can use a basis of the boundary theory, which is labeled (j1,m1), (j2,m2), . . . , (jN ,mN ) . (39) We can now turn to the specification of the interaction algebra. We begin by specifying a general interaction algebra, which consists of all possible interactions between an environment and a quantum horizon. Then we will specialize to a subalgebra which excludes the case of an O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 9 asymmetric environment or external geometry. The general interaction algebra, Agint,mustbe generated by operators that
(1) act simultaneously on the environment and system Hilbert space; (2) act locally.
A sufficient set of generators for Agint consists of:
• addition or removal of a puncture; • braidings of the punctures.
The Gauss’ law constraint that ties the labels on punctures to the labels on edges of spin networks that meet them implies that each of these involve changes to both the surface and environment state. We note that the braidings generate a group, called the braid group and together with adding and removing punctures, these generate the tangle algebra discussed by Baez in [15]. Let us consider, for example, a generator in Agint that corresponds to braiding two punctures, which we will label i = 1, 2. It has the effect of rotating, with a positive orientation, the two punctures around each other, returning them to their original positions. This braids the two edges of the bulk spin network, which changes the diffeomorphism class of the bulk state. The action of H int on HS due to a braiding of the punctures can be shown to modify each basis state (39) by a phase. As the punctures are unchanged, it is represented by an operator Bˆ in HS . A simple computation in the quantum Chern–Simons theory [14] shows that this is 12 {ji } realized by ˆ B12 ◦ (j1,m1), (j2,m2),...,(jN ,mN ) + 2πı(m1 m2 = e k (j1,m1), (j2,m2),...,(jN ,mN ) . (40) The general interaction algebra does not take into account any symmetries, hence it describes the general case, in which the black hole horizon may have non-vanishing multiple moments.4 If we want to specialize to the case of a black hole with symmetries, we must impose condi- tions on the interaction algebra, for a symmetric subsystem cannot exist stably in an asymmetric environment. In principle we could code various kinds of symmetries in the choice of interac- tion algebra. But it suffices to consider the simplest choice, which is the subalgebra Aint ⊂ Agint which consists of adding and removing punctures plus a symmetric sum of braiding operations, ˆT ˆ B = Bij , (41) i P1. The elements of Gq do not change the area of the horizon. 4 In the isolated horizon approach [5] angular momentum and multiple moments are coded as classical parameters. Here we seek instead to define different subspaces of the horizon Hilbert space associated with black holes with different macroscopic properties. 10 O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 5 We take this to mean that the elements of Gq commute with the area operator. This is trivially true for the action of SO(3) on the classical spacetime and we assume that it is also true for the action of G . Because of the property P1, the elements of G will map the spaces HS into q q {ji } themselves. A necessary condition for Eq. (33) to hold, i.e., for Gq to commute with the action of H int, is that G commutes with the action of H int on each of the spaces HS . q {ji } On each HS the braiding generators in H int act as shown in Eq. (40). The commutant of {ji } int A is then easy to find. It is generated by all those operators that just permute the mi values for given spins ji : ...,(j,m1), . . . (j, mK ), . . . −→ ...,(j,mπ(1)),...(j,mπ(K)),..., (42) for some π ∈ PK , the permutation group of K objects. Assuming Gq to be unitary we then have ⊂ Gq PKj , (43) j where the product is over all the different j’s in {ji}. We now assume that Gq also shares the following property with SO(3): P2. Gq acts transitively on the punctures, i.e., for every two punctures there is an element in Gq that connects the two. Because of Eq. (43), we know that Gq is a subgroup of = P PKj . (44) j For Gq to act transitively, it is necessary that the bigger group P acts transitively on the punctures. Since P is the product of permutation groups, this is only possible if P coincides with just one permutation group. It follows that there is one j such that Kj = N. (45) HSchwarzschild V All the punctures have the same spin j and the Hilbert space A consists of those {j} for which all the j’s coincide. It is then straightforward to show that the dimension of the V{j} varies sharply with the spin j and is dominated by the lowest spin jmin. 4. Conclusions In this paper, we used a microscopic analogue of the classical SO(3) symmetry of a Schwarz- schild black hole to restrict the loop quantum gravity state space HA of black holes of area A to a smaller subspace, left invariant by the symmetry. The method we used is the noiseless subsys- tems of quantum information theory, in which the symmetry of the dynamics implies a non-trivial commutant of the interaction algebra of the system. The construction of the state space HA in loop quantum gravity is a hybrid construction that requires the imposition of black hole horizon conditions at the classical level. One eventually 5 In principle we could consider a weaker requirement which is that the elements of Gq do not change the expectation value of the area. This would be an interesting extension of the usual problem to study. O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 11 wants to be able to identify black hole states directly in the quantum theory and derive the clas- sical geometry in the appropriate classical limit. While the present work has the same hybrid character (we start from HA), it indicates that properties of the quantum states can be inferred from the dynamics algebraically, without need for a classical geometry. This approach also clarifies the important but subtle question of how symmetries can arise from a diffeomorphism invariant state. All the states discussed here are invariant under spatial diffeomorphisms. The point is that the commutant, Gq takes physical, diffeomorphism invariant states, to other distinct physical, diffeomorphism invariant, states. The transitive permutations that we discuss translate the black hole horizon with respect to the spin networks that define the external geometry. This is physically meaningful because these transformations are not sub- groups of the diffeomorphisms, instead they act on the space of diffeomorphism invariant states. Thus, we have a realization of the proposal in [13] that in a diffeomorphism invariant theory a symmetry can only arise as a motion of one subsystem with respect to another, and its generators must live in the commutant of the interaction algebra, defined by the splitting of the universe into subsystems. This proposal has ramifications that might lead to a better understanding of how the classical properties of black holes arise from the quantum geometry of the horizon. We also see in a sim- ple example, how properties of classical spacetime geometries can emerge from exact quantum geometries, by making use of the insights gained from the study of similar questions in quantum information theory. We close by mentioning several queries that may be made regarding the construction used here. • What about the connection to quasi-normal modes? Can we understand it as other than a coincidence? It may be possible to understand how the quasi-normal modes arise from the states which transform non-trivially under those generators of the commutant that become rotations in the A →∞limit. • What about rotating black holes? It is of interest to use these methods to characterize the states of a rotating black hole, or a horizon with multiple moments, in terms of the same language. We expect that this will involve other choices of the interaction algebra than that made here. We may note that attempts to extend the connection between quasi-normal modes and entropy to rotating black holes have not succeeded; this is something that needs explanation. Note also that while we employ here the same boundary conditions as the isolated horizon for- malism, the intension is different. In particular, in the isolated horizon picture, rotation and other multiple moments are fixed in the classical description. The isolated horizon boundary conditions do not change, but the functions on the phase space that correspond to different observables will depend on the multiple moments. In our picture, we seek to recover the generators of rotation, and hence the conjugate conserved quantities, as operators in a single Hilbert space. Hence here, unlike the pure isolated horizon case, we expect that a single Hilbert space contains the states corresponding to black holes with all values of angular momentum. • HAll HSchwarzchild What about the states in A that are not in A ? These are states whose pictures give the right area, but are not in the noiseless subsystem we identify. But they arose from a quantization of the isolated horizon boundary conditions. We cannot give a definitive answer without more dynamical input, particularly a Hamiltonian on the horizon states. But we can argue that they must correspond to either non-static, classical metrics or they correspond to no classical geometry. 12 O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 Given that the boundary conditions do not distinguish rotating from non-rotating black holes, some of theses states must be associated with rotating black holes. Apart from these, it must be that most of these states are associated with semiclassical states of non-stationary horizons, corresponding to distorted black holes. These will be excited, transient states of the black hole, which classically decay to stationary states, emitting gravitational radiation. How do we know HAll HAll that such states are included in A ? All that distinguishes states in A is that there is a horizon of a particular area. This will include all such configurations, whether belonging to stationary states, semiclassical states corresponding to non-stationary horizons as well as quantum states that have no semiclassical approximation. In particular, these must include states corresponding to quasi normal modes. Beyond those categories, our argument implies that the remaining states couple strongly to the random noise coming from the environment of the black hole. They are then states that cannot play a role in the semiclassical limit. The only definitive way to distinguish the different kinds of states is dynamically. If we had an effective Hamiltonian for the horizon and near horizon geometry, then the non-stationary states would have to have higher energy than the stationary states, corresponding to their potential for radiating energy to infinity in gravitational radiation. Hence, we can deduce from the fact that in the classical theory non-stationary states radiate, that those states would be suppressed in a HAll Boltzman weighted, equilibrium ensemble. In the ensemble A they appear unsuppressed, but that is only because no dynamics are invoked, so states of different energies are counted as having the same weight. Given that there are many more non-stationary classical configurations than stationary ones, it is reasonable that, for a given horizon area, those states that correspond to static black holes must HAll HSchwarzchild appear to be a minority, as is the case if we compare A with its subspace A .But, HAll HSchwarzchild it is incorrect to argue just on the basis of the fact that A is larger than A , that a typical state corresponding to an equilibrium, static black hole will be in the latter rather than the former, because dynamics has yet to be invoked. The point of the construction of this paper is then that, even in the absence of an explicit Hamil- tonian, we can use symmetry properties of the Hamiltonian to apply an argument from quantum information to deduce which subspace of states will emerge in the low energy limit as corre- sponding to stationary configurations. The fact that, in the absence of dynamics, these are a minority of the ensemble of states with a given area must be expected and is not an objection against the construction presented here. • Finally, we note that the same method could be applied within quantum cosmology, to study how states that correspond to symmetric universes emerge out of a more general Hilbert space, corresponding to the full theory [16]. Acknowledgements We are very grateful to Paolo Zanardi for pointing out an oversight in an earlier version of this paper. We also thank Abhay Ashtekar and Jerzy Lewandowski for very useful discussions clarifying the treatment of symmetries in the isolated horizon formalism. References [1] J.D. Bekenstein, Black holes and entropy, Phys. Rev. D 7 (1973) 2333. [2] S.W. Hawking, Particle creation by black holes, Commun. Math. Phys. 43 (1975) 199. O. Dreyer et al. / Nuclear Physics B 744 (2006) 1–13 13 [3] L. Smolin, Linking topological quantum field theory and nonperturbative quantum gravity, J. Math. Phys. 36 (1995) 6417, gr-qc/9505028. [4] K. Krasnov, On quantum statistical mechanics of a Schwarzschild black hole, Gen. Relativ. Gravit. 30 (1998) 53–68, gr-qc/9605047; C. Rovelli, Black hole entropy from loop quantum gravity, gr-qc/9603063. [5] A. Ashtekar, J. Baez, K. Krasnov, Quantum geometry of isolated horizons and black hole entropy, gr-qc/0005126; A. Ashtekar, J. Baez, A. Corichi, K. Krasnov, Quantum geometry and black hole entropy, Phys. Rev. Lett. 80 (1998) 904–907, gr-qc/9710007. [6] O. Dreyer, Quasinormal modes, the area spectrum, and black hole entropy, Phys. Rev. Lett. 90 (2003) 081301, gr-qc/0211076. [7] G. Immirzi, Nucl. Phys. B (Proc. Suppl.) 57 (1997) 65, gr-qc/9701052. [8] S. Hod, Phys. Rev. Lett. 81 (1998) 4293, gr-qc/9812002. [9] L. Motl, An analytical computation of asymptotic Schwarzschild quasinormal frequencies, gr-qc/0212096; L. Motl, A. Neitzke, Asymptotic black hole quasinormal frequencies, hep-th/0301173. [10] K.A. Meissner, Black hole entropy in loop quantum gravity, gr-qc/0407052; M. Domagala, J. Lewandowski, Black hole entropy from quantum geometry, gr-qc/0407051. [11] S. Alexandrov, On the counting of black hole states in loop quantum gravity, gr-qc/0408033. [12] P. Zanardi, M. Rasetti, Phys. Rev. Lett. 79 (1997) 3306; E. Knill, R. Laflamme, L. Viola, Phys. Rev. Lett. 84 (2000) 2525; J. Kempe, D. Bacon, D.A. Lidar, K.B. Whaley, Phys. Rev. A 63 (2001) 042307. [13] F. Markopoulou, D. Poulin, in preparation. [14] D.P. Arovas, R. Schrieffer, F. Wilczek, A. Zee, Statistical mechanics of anyons, Nucl. Phys. B 251 (1985) 117; S. Rao, Anyon primer, hep-th/9209066. [15] J. Baez, Quantum gravity and the algebra of tangles, Class. Quantum Grav. 10 (1993) 673–694, hep-th/9205007. [16] E. Poisson, personal communication. [17] T. Tamaki, H. Nomura, Ambiguity of black hole entropy in loop quantum gravity, hep-th/0508142. [18] E.R. Livine, D.R. Terno, Quantum black holes: Entropy and entanglement on the horizon, gr-qc/0508085. Nuclear Physics B 744 (2006) 14–33 Finite volume effects for the pion mass at two loops Gilberto Colangelo ∗, Christoph Haefeli Institut für Theoretische Physik, Universität Bern, Sidlerstr. 5, 3012 Bern, Switzerland Received 16 February 2006; accepted 9 March 2006 Available online 30 March 2006 Abstract We evaluate the pion mass in finite volume to two loops within chiral perturbation theory. The results are compared with a recently proposed extension of the asymptotic formula of Lüscher. We find that contribu- tions, which were neglected in the latter, are numerically very small at the two-loop level and conclude that for Mπ L 2, L 2 fm the finite volume effects in the meson sector are analytically well under control. © 2006 Elsevier B.V. All rights reserved. 1. Introduction Numerical simulations in lattice QCD are bound to rather small volumes. When one deter- mines the hadron spectrum and other low energy observables one has to understand and properly account for the volume dependence in order to correctly interpret the numerical data. Analytical methods which allow one to predict the size of the finite volume effects are particularly useful in this respect. In the case of hadron masses, there are two different methods to analytically evaluate the finite volume effects: the asymptotic formula derived by Lüscher [1], and chiral perturbation theory (ChPT) in finite volume [2–4]. Lüscher’s formula relates the volume dependence of the mass of a hadron to an integral over the π-hadron forward scattering amplitude. Knowledge of the latter scattering amplitude immediately translates into an estimate of finite volume effects for the hadron mass. Since the integral is dominated by the low-energy region, one can rely on the chiral representation of the scattering amplitude in order to numerically evaluate the integral. In this manner one makes use of ChPT only in infinite volume and obtains an estimate only of the leading exponential term in the finite volume dependence, the term of the order exp(−Mπ L), where L is the box size. Alternatively one can perform the calculation of the hadron mass in * Corresponding author. E-mail address: [email protected] (G. Colangelo). 0550-3213/$ – see front matter © 2006 Elsevier B.V. All rights reserved. doi:10.1016/j.nuclphysb.2006.03.010 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 15 ChPT in finite volume and obtains an estimate also of the terms which are exponentially sub- leading. A detailed study and comparison of the two approaches has been performed in [5,6] for the case of the pion mass. This case is particularly interesting both because it is a simple hadron to be studied on the lattice and because the ππ scattering amplitude, which is needed in the Lüscher formula, is now known to next-to-next-to-leading order [7–10] in the chiral expansion. It was therefore possible to study numerically how well the leading exponential term dominates the series, and how fast the chiral expansion (of the leading exponential term) converges. The somewhat surprising result was that the leading term in both series receives large corrections both from subleading exponentials as well as from the next-to-leading order in the chiral expansion. The next-to-next-to-leading order chiral correction, on the other hand, was found to be rather small, which indicates that—at this level—the chiral series has started converging well. The fact that the Lüscher formula appeared to give the numerically dominating term only for volumes so large that the finite volume correction itself has become negligible appeared to be more worrisome. At first sight this might have lead one to conclude that in cases of practical interest one could not rely on this extremely convenient formula in order to evaluate reliably the finite volume effects. We believe, however, that such a negative conclusion is unjustified and have proposed a resummation of the Lüscher formula [11,12] which retains all the convenience of the original one but does not suffer from the same large corrections. This resummation can be understood in simple terms: Lüscher has shown that the leading exponential correction comes from a radiative correction in which the emitted virtual light particle (the pion, in QCD) goes around the world once (thanks to the periodic boundary conditions) before being reabsorbed. The resummed formula takes into account all other possible ways which the pion has to go around the world (go along the diagonal, or go around more than once, etc.). As was shown in [12] this resummed formula exactly reproduces the one-loop ChPT calcula- tion of the pion mass in finite volume if one inserts in the integral the leading chiral representation of the ππ scattering amplitude. Inserting the next-to-leading chiral representation of the scatter- ing amplitude one reproduces the full two-loop calculation√ of the pion√ mass in ChPT in finite ¯ ¯ volume up to corrections of order exp(−ML) with M ( 3 + 1)/ 2Mπ . Despite this im- provement of the algebraic accuracy of the resummed formula with respect to the Lüscher’s one, it is essential to show that the resummation is numerically effective and that the corrections to it are small. In order to make this check we have now performed a full two loop calculation of the pion mass in finite volume in ChPT. As we will demonstrate in this article, the corrections which are not captured by the resummed asymptotic formula are negligibly small for Mπ L 2. The results of the present two-loop calculation were anticipated in [13,14]. Here, we give further details about the calculation and the results. One of the reasons for performing this further investigation of finite volume effects for hadron masses is that one may view Lüscher’s formula (or its resummed version) as a way to determine on the lattice a scattering amplitude, cf. [12]. While it is true that the scattering amplitude is seen here in an exponentially suppressed effect, it is also true that the direct calculation of a scattering amplitude on the lattice is much more difficult than that of a mass. This indirect method to extract a scattering amplitude may in some cases turn out to be more practical, and we find it worthwhile to discuss its theoretical feasibility. Indeed if the part of the finite volume corrections which is related to the scattering amplitude is not strongly dominating with respect to the rest, this is not a viable method to determine the scattering amplitude. The conclusion reached in this paper is therefore encouraging in this respect and can be used to estimate in which region of the (Mπ ,L) 16 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 plane the finite volume correction to a hadron mass is given to a good approximation by the integral containing the scattering amplitude. The two-loop calculation is interesting in its own right, also from a technical viewpoint. To date, a number of finite volume calculations have been performed at one-loop order [15–21],but as far as we know the only finite volume two-loop calculation which has been performed until now is for the quark condensate [22], which is a much simpler calculation. In the related context of finite temperature field theory there are a few examples of calculations beyond one loop [23, 24]. We wish to briefly mention related work. Most notably, the asymptotic formula may also be applied to the nucleon mass [25], see Koma and Koma [26] as well as [14] for recent work. Braun, Pirner and Klein have evaluated the volume dependence of the pion mass based on a quark–meson model [27]. In the framework of a lattice regularized ChPT finite volume effects have been addressed by Borasoy et al. in Ref. [28]. The outline of the article is as follows. In Section 2 we set the notation and remind of the basic assumptions for an application of ChPT in finite volume. Section 3 is devoted to outline the calculation and state the main results. In Section 4 we show the explicit expressions which have been used for the numerical analysis in Section 5. We conclude with a summary. 2. Preliminaries In this section we shall set the notation and the basic definitions for the two-loop calculation. 2.1. ChPT in finite and in infinite volume Chiral perturbation theory (ChPT) is the effective theory for QCD at low energies. If we first restrict ourselves to the infinite volume case, the effective Lagrangian of QCD for two light flavours at low energies consists of an infinite number of terms [7], Leff = L2 + L4 + L6 +···. (1) As we wish to calculate the pion mass, an on-shell quantity, external fields can be dropped in Leff. We work in the isospin symmetry limit mu = md in Euclidean space–time, and for the choice of the pion fields we use the nonlinear sigma model parameterization. We have F 2 L = u u − χ+, (2) 2 4 μ μ with φ φ2 U = σ + i ,σ2 + = , 2 1 F √ F π 0 2π + φ = √ = φiτ i, 2π − −π 0 = † † =− † = † = † † + † uμ iu ∂μUu iu∂μU u uμ,χ+ u χu uχ u, 1 χ = 2Bmˆ 1, mˆ = (m + m ), (3) 2 u d G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 17 2 with u = U. The symbol A denotes the trace of the two-by-two matrix A. The Lagrangian L4 can be written as [7] 4 L4 = iPi +···, (4) i=1 where 1 1 P =− u u 2,P=− u u 2, 1 4 μ μ 2 4 μ ν 1 2 i P =− χ+ ,P= u χ− , (5) 3 16 4 4 μ μ with † † † χ−μ = u ∂μχu − u∂μχ u. (6) The ellipsis in Eq. (4) denotes terms that do not contribute to the pion mass. The low energy constants i are divergent and remove the ultraviolet divergences generated by one-loop graphs from L2. The complete effective Lagrangian L6 with its divergence structure at d = 4 has been con- structed in [29,30]. As will be discussed in Section 3, terms from L6 merely renormalize the pion mass in infinite volume and do not contribute to finite volume corrections which we are interested in. Thus, we refrain from showing it here. Given the effective Lagrangian and the parameteriza- tion for the pion fields, it is straightforward to calculate the pion mass to two-loops. We refer to [9], where one also finds a detailed discussion of the renormalization procedure (two-loop diagrams in ChPT are discussed in [31]). The effective framework is still appropriate, when the system is enclosed by a large box of size V = L3. We refer to the literature for the foundations [2–4] and a recent review [11]. Here, we only recall a few fundamental results which guided the present calculation: the volume has to be large enough, such that ChPT can give reliable results, 2Fπ L 1. The value of the parameter Mπ L determines the power counting for the perturbative calculation: if Mπ L 1 one is in the “p-regime” in which 1/L counts as a small quantity of order Mπ .IfMπ L 1 one is in the “ - 2 regime” and 1/L is a quantity of order Mπ . In both cases the effective Lagrangian is the same as in the infinite volume. In this article we only consider the “p-regime”, where the system is distorted mildly and the only change brought about by the finite volume is a modification of the pion propagator due to the periodic boundary conditions of the pion fields1 0 0 G x , x = G0 x , x + nL , (7) n∈Z3 with G0(x) the propagator in infinite volume. 2.2. Basic definitions In Euclidean space–time the propagator is defined through the connected correlation function ab = a b G(x)δ φ (x)φ (0) L, (8) 1 Throughout we denote by the volume the three-dimensional volume V = L3, whereas the time direction is not compactified. 18 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 where a, b stand for isospin indices of the pion fields and the subscript L in the correlation function denotes that it is evaluated in finite volume. We also have 0 − dp φ1(x)φ1(0) = L 3 eipx G p0, p , L 2π p 0 −1 2 2 0 G p , p = M + p − ΣL p , p , (9) where the momenta p can only have discrete values 2π p = n, n ∈ Z3, (10) L and M2 = 2Bmˆ is the tree-level pion mass in infinite volume. The pion mass in finite volume MπL is now defined by the pole equation −1 G(pˆL) = 0, for pˆL = (iMπL, 0). (11) 3. Outline of calculation and statement of results For a large volume, the finite size effects are expected to be small, such that the pole equation can be solved perturbatively. We outline the calculation and proceed with the main results. More details about the calculation are relegated to later sections or Appendices A and B. 3.1. One-loop result Since the effective Lagrangian remains unchanged when going to the finite volume, we can immediately write down the Feynman diagrams which contribute to the self-energy at two-loops, see Fig. 1, the only difference with respect to an infinite volume calculation is that the propagators Fig. 1. Self-energy graphs to two-loops in ChPT. A spline corresponds to a periodified propagator, whereas those without correspond to an infinite volume propagator, cf. Eq. (7). Normal vertices come from L2, squared vertices from L4 and the circle-crossed from L6. G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 19 need to be periodified, cf. Eq. (7). At leading order, the graphs 1(a) and 1(e) need to be evaluated and the self-energy admits the form 1 3 M4 1 Σ p0, = G( ) − M2 − p2 − π + O , L p 2 0 π 2 3 2 4 (12) Fπ 2 Fπ Fπ with G(0) the value of the finite volume propagator at the origin. It contains a logarithmic diver- gence due to the contribution from the term n = 0 in Eq. (7). In dimensional regularization 1 1 (1 − d/2) − G (0) = dd p = Md 2. (13) 0 (2π)d p2 + M2 (4π)d/2 The remaining terms with n = 0 are finite and may be expressed in terms of a kinematical func- 2 2 tion g1(M , 0,L), 2 G(0) = G0(0) + g1 M , 0,L , ∞ 2 2 dτ − 2 n L g M2, ,L = e τM − . 1 0 d/2 exp (14) (4πτ) = 4τ 0 n 0 For a derivation of Eq. (14), we refer to [32].InAppendix A we provide a different derivation, based on a contour integration analysis. The pole of the gamma function in Eq. (13) in four dimensions is absorbed in a renormalization of the low-energy constant 3. One readily verifies that inserting Eq. (14) into Eq. (12) yields for the leading finite volume shift [2] 1 1 M = M + g M2 , ,L + O . πL π 1 2 1 π 0 4 (15) 4Fπ Fπ The separation of the cut-off and the volume-dependence is as expected: finite volume corrections do not generate new UV-divergences. At leading order the finite volume corrections could be isolated immediately. This will not be the case at the two-loop level. The graph (d) in Fig. 1 does not factorize in pure one-loop integrals and further steps need to be performed. We refer to Appendix A for further details. 3.2. Minimal set of periodified propagators In the evaluation of the Feynman diagrams, it is mandatory to make use of the following result. In order to evaluate a graph consisting of L loops, one needs to consider only a certain set of L propagators as the periodified, finite-volume ones—the others can be taken as infinite-volume propagators. In the following we briefly discuss this statement (see also p. 18ff in Ref. [1]). Consider an arbitrary self-energy graph with L loops, I internal lines and V vertices. Since the number of loops is the number of independent integrations over momenta, we have L = I − V + 1. (16) For every line the propagator is an infinite sum of terms characterized by an integer vector n(). The loop graph itself is an infinite sum of terms which can be identified by a set of integer vectors {n()} assigned to all the internal lines. As remarked by Lüscher this structure can be 2 The second argument of g1 is the temperature, which we keep zero. Notation as in Ref. [2]. 20 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 viewed as a gauge field on the self-energy graph. Gauge transformations are defined as3 n () = n() + f() − i() , (17) where (v) is some field of integer vectors and i() (f ()) is the initial (final) vertex of the line . One verifies immediately that two contributions of a Feynman graph which differ only by a gauge transformation yield the same mathematical expression. In the calculation of the loop graph what matters is the sum over representatives of gauge equivalence classes. A convenient representa- tive can be found, e.g., by adjusting (i()) and (f ()) iteratively, such that n() = 0 for as many internal lines as possible. This can be achieved for V − 1 lines. Therefore, there remain I − V + 1 internal lines where the periodified propagator has to be inserted—which coincides with the number of loops of the graph. This shall be our minimal set of periodified propagators. In Fig. 1, we have attached a spline to a periodified propagator, whereas lines without a spline correspond to an infinite volume propagator. 3.3. Two-loop result It is convenient to split the sum over the equivalence classes into three parts4: Σ(0): n() = 0 ∀ (pure gauge), Σ(1): n() = 0 ∀ except for one line ¯ (simple gauge), (2) ¯ ¯ Σ : n() = 0 ∀ except for two lines 1, 2. (18) The pion mass in finite volume to two-loops then admits the form 2 = 2 − (1) − (2) MπL Mπ Σ Σ , 2 = 2 − (0) Mπ M Σ , (19) where, using λπ = Mπ L, we get (1) 3 Σ = Ip + Ic + O ξ , (20) ∞ 2 ∞ √ M m(n) − + 2 I = π √ dy F (iy)e n(1 y )λπ , p 2 π (21) 16π λπ n n=1 −∞ ∞ ∞ √ ∞ − ˜+ 2 iM2 m(n) e n(s y )λπ I =− π √ dy ds˜ disc F (s,˜ 1 + iy) , (22) c 3 ˜ + π 32π λπ = n s 2iy n 1 −∞ 4 2 with m(n) ≡ number of integer vectors z with z = n. Ip denotes the pole and Ic the cut contri- bution, whose meaning and precise definitions will be explained below. The expression for Σ(2) is more cumbersome: (2) = 2 2 9 ˜ 2 − 1 ˜ ∂ ˜ + + O 3 Σ Mπ ξ g1(λπ ) λπ g1(λπ ) g1(λπ ) Δ ξ , (23) 8 8 ∂λπ 3 We borrow the notation from Ref. [1]. 4 Note the slight difference in the definition of simple fields with respect to [1]. We do not require |n()|=1. G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 21 ˜2 and we merely note that it can be split into products of one-loop contributions terms with g1 and a pure two-loop part, denoted by Δ. We have introduced the abbreviation M2 ξ = π . 2 2 (24) 16π Fπ A detailed derivation of these results will be given in the subsequent sections and we confine ourselves to a few comments at this stage: In Ip one recovers the asymptotic formula of Lüscher, if one restricts the sum to the first addendum. Its extension to the present form has already been suggested in [11] and was applied in [12]. The function Fπ (ν) denotes the isospin zero ππ scattering amplitude in the forward kinematics. It contains cuts due to the two-pion intermediate state. In the derivation of the asymptotic formula, Lüscher consistently dropped contributions arising from the cuts, since in a large volume expansion they are beyond the order of accuracy he aimed at. Here, we take them into account up to two-loops, relying on the chiral representation of ππ scattering amplitude. The outcome of the analysis is summarized in Ic in which the same symbol Fπ , this time with two arguments appears: with the latter we mean the scattering ampli- tude of two on-shell pions into two off-shell pions in the forward kinematics configuration. The contributions from the cuts can still be written as integrals over the ππ scattering amplitude with an exponential weight, as in the asymptotic formula. Notice however, that this formula explicitly relies on the chiral representation of the scattering amplitude and is only valid up to two-loops, whereas the Lüscher formula holds to all orders. Contributions from two pion propagators in finite volume are ultimately captured in Σ(2), being expressed in terms of a dimensionless function g˜1(λπ ) and a numerically small correction Δ arising from graph 1(d). Both are explicitly given in Appendix A and Eq. (44), respectively. The pure gauge contributions are not volume dependent and merely renormalize the pion mass, cf. Eq. (19). A detailed discussion of this calculation can be found in [33,9], with which we agree—this was a useful, nontrivial check for our calculation. Fig. 2. Skeleton diagrams representing (a) Eq. (25) and (b) Eq. (B.3). The blob in (a) stands for the 4-point function of ππ scattering in infinite volume and in (b) for a subtracted 6-point function of πππ scattering in infinite volume. The 2 double-line with the spline is a finite volume propagator with the physical pion mass Mπ . 22 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 3.4. Self-energy to first order: Σ(1) The simple fields can be summed up in closed form and may be represented by a skeleton diagram, see Fig. 2(a) ∞ 1 d4q √ (1) iq1 nL 2 Σ = m(n)e G q Γππ(p,q,ˆ −ˆp,−q), (25) 2 (2π)4 0 n=1 with pˆ = (iMπ , 0) and Γππ(p,q,ˆ −ˆp,−q) the 4-point function of ππ scattering in the forward scattering kinematics. Note that the 4-point function in Eq. (25) is an off–shell amplitude and that (2) we evaluate it at pˆ and not at pˆL. The difference is taken into account in Σ . The representation of Σ(1) in Eq. (25) has already been used by Lüscher and eventually led him to the asymptotic 2 formula [1]. He then discussed the contribution of the pole of the propagator G0(q ) that one meets at = 2 + 2 + 2 = q1 i Mπ q⊥ q0 ,q⊥ (q2,q3). (26) Above this pole the singularities come from the cuts of the propagator and the 4-point function Γππ(p,q,ˆ −ˆp,−q). These start from =− ˆ + 2 2 =− ˆ − 2 2 − 2 2 s (p q) 4Mπ ,u (p q) 4Mπ , q 9Mπ . (27) Lüscher showed that the pole contribution is dominating with respect to those coming from the cuts and neglected the latter. In fact, his discussion involved no further assumptions about the 4-point function and remains true at every order of the perturbative expansion. Since our goal is to test the asymptotic formula beyond the leading exponentials, we wish to work out the impact of the contributions that were dropped by Lüscher. Doing this at the two-loop level is rather straightforward, as we will now show. Up to O(p4), the four-point function can be decomposed into a combination of functions which have either a singularity in s or in u. Since Σ(1) is symmetric in s and u, we may write ¯ Γππ(p,q,ˆ −ˆp,−q) in terms of a function Γππ(s, pq)ˆ , which depends only on the variables s ˆ 2 and pq, and which has cuts at s 4Mπ , 1 Γ (p,q,ˆ −ˆp,−q) = Γ¯ (s, pq)ˆ + O . ππ ππ 6 (28) Fπ Fig. 3. Integration contour in the complex q1 plane with the pole from the pion propagator and the branch cut from the ππ scattering amplitude. G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 23 = + ∈ R 2 We write q1 x iy and for the domain s , s 4Mπ we have 2 − 2 − 2 − 2 2 =− y x q0 q⊥ 3Mπ ,xyq0Mπ . (29) In particular we observe that in the complex q1 plane the cut starts above the pole of Eq. (26), as illustrated in Fig. 3. Performing now the contour integration in the upper half plane we obtain two terms: the first one, Ip, comes from the pole and is simply its residuum. The other comes from the integral along the new integration path. The latter contribution vanishes as we push the integration lines to infinity, except for the path which goes around the cut, to be denoted by Ic in the following, (1) Σ = Ip + Ic. As was shown by Lüscher [1], the former can be simplified considerably and can be brought in the usual form of the asymptotic formula ∞ 2 ∞ √ M m(n) − + 2 I = π √ dy F (iy)e n(1 y )λπ , p 2 π 16π λπ n n=1 −∞ with Fπ (iy) the ππ forward scattering amplitude. Restricting the sum to the first term we recover Lüscher’s formula [1].ForIc we find ∞ √ 2 iq1 nL 1 dq d q⊥ e − I = m(n) 0 dq disc Γ¯ (s, pq)ˆ + O F 6 , (30) c 4 1 2 + 2 ππ π 2 = (2π) Mπ q n 1 [ 2 ] s 4Mπ ¯ ¯ where disc[Γππ(s, pq)ˆ ] denotes the discontinuity of Γππ along the cut. It is now convenient to shift the integration path in q0 from Im(q0) = 0toIm(q0) =−iMπ . Along this path we have =¯ − =−¯2 − 2 − 2 ¯ ∈ R q0 q0 iMπ ,sq0 q1 q⊥, q0 , (31) and the integration over q1 = x + iy falls onto the imaginary axis, = = 2 +¯2 + 2 x 0,yy0 4Mπ q0 q⊥, ∞ ∞ ¯ 2 √ [ ¯ ˆ ] i dq0 d q⊥ −y nL disc Γππ(s, pq) Ic = m(n) dy e . (32) 2 (2π)4 M2 + q2 n=1 π y0 Next, we change the integration variable from q1 to s, ∞ ∞ − +¯2+ 2 n(s q q⊥)L ¯ i dq¯ d2q⊥ e 0 disc[Γ (s, pq)ˆ ] I =− m(n) 0 ds ππ , (33) c 4 +¯2 + 2 1/2 + ¯ 2 = (2π) 2(s q q⊥) s 2iMπ q0 n 1 2 0 4Mπ and make use of 2 √ d q⊥ 1 − 2+ 2 1 − e n(μ q⊥)L = √ e μ nL, (34) 2 2 2 1/2 (2π) 2(μ + q⊥) 4π nL to carry through the integration over q⊥ and to end up with Eq. (22). 24 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 3.5. Self-energy to second order: Σ(2) The term Σ(2) is more complicated to manipulate, although it is suggestive to think that it may be related to an integral over the 3π → 3π amplitude. A similar question has been asked in the framework of finite temperature QCD [23]. We discuss the situation in finite volume in Appendix B. Since these considerations have not lead us to a nice and compact representation for Σ(2) we have written it as in (23) where we only have split the factorizable two-loop contributions from the rest. The latter, which we have denoted by Δ we have evaluated numerically. We conclude this section with a remark concerning the large L behaviour of the two-loop results. In the large volume limit, the contributions from Eq. (21) behave according to 1 −λπ lim Ip ∼ e . (35) L→∞ 3/2 λπ Similarly, we may evaluate the large volume behaviour of the terms occurring in Eqs. (23), (22).It turns out that√ these, besides√ being exponentially suppressed and behaving at least as exp(−αλπ ), with α = ( 3 + 1)/ 2, are also suppressed√ by a power of λπ . Unfortunately this suppression is not particularly strong, being of 1/ λπ with respect to Eq. (35). While we can conclude that√ the resummed Lüscher formula is dominating in comparison to other two-loop diagrams in a 1/ λπ expansion, i.e., ¯ −αλπ MπL = MπL 1 + O e / λπ , ¯ 1 MπL = Mπ − Ip, (36) 2Mπ we are convinced that the most important test of the resummed formula is the numerical one, which we will discuss in the following. 4. Summary of analytical results In this section we shall give the analytical results in explicit form, i.e., insert the chiral rep- resentation of the amplitudes which appear in the formulae given in the previous section and express our results in terms of a few basic integrals. In Eqs. (19)–(24) we have split the finite size (1) (2) effects of the pion mass into Σ = Ip + Ic and Σ . For the former two we find ∞ = 2 m(n)√ 1 (2) + (4) Ip Mπ ξ IM ξIM , n λπ π π n=1 ∞ = 2 m(n)√ 1 2 ˜ Ic Mπ ξ Ic, (37) n λπ n=1 where ξ is the chiral expansion parameter defined in Eq. (24) and the expressions I (2) , I (4) have Mπ Mπ already been given in Refs. [6,12], and we reproduce them here for completeness: I (2) =−B0, Mπ (4) 0 55 ¯ 8 ¯ 5 ¯ ¯ 2 112 8 ¯ 32 ¯ I = B − + 41 + 2 − 3 − 24 + B − 1 − 2 Mπ 18 3 2 9 3 3 13 16 40 + R0 − R1 − R2, (38) 3 0 3 0 3 0 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 25 where √ 0,2 0,2 0 2 B ≡ B ( nλπ ), B (x) = 2K1(x), B (x) = 2K2(x)/x, (39) k and the integrals R0 are defined as ∞ √ Re − +˜2 k even, Rk = dy˜ y˜k e n(1 y )λπ g(2 + 2iy)˜ for (40) 0 Im k odd, −∞ ¯ 2 where g(x) is related to the standard one-loop function J(xMπ ) through = 2 ¯ 2 g(x) 16π J xMπ , (41) ¯ 2 and J(xMπ ) given in Eq. (A.10). For instance for x<0, we have σ − 1 g(x) = σ log + 2, with σ = 1 − 4/x, (42) σ + 1 and elsewhere defined through analytic continuation. In Ref. [12] also the next-to-next-to-leading order term in I , namely ξ 2I (6) has been given and we will use it in our numerical analysis. The p Mπ ˜ coefficient Ic can be expressed as a combination of basic integrals 1 I˜ = 112C0,2 + 37C1,0 − 40C1,2 − 4C2,0 , c 3 ∞ ∞ √ − ˜+ 2 e n(s y )λπ 4 1/2 Cj,k = dy ds˜ 1 − s˜j yk. (43) s˜2 + 4y2 s˜ −∞ 4 Notice that the latter expression is obtained from the chiral representation of the off-shell ππ scattering amplitude for forward kinematics. This is parameterization dependent, and the result given here is obtained for the parameterization discussed in Section 2.1. The dependence on the parameterization must cancel in the full result. The self-energy to second order has already been introduced in Eq. (23) and reproduced here for convenience (2) = 2 2 9 ˜ 2 − 1 ˜ ∂ ˜ + + O 3 Σ Mπ ξ g1(λπ ) λπ g1(λπ ) g1(λπ ) Δ ξ , 8 8 ∂λπ with ˆ 2 2 7 p Δ = 16π 4p˜μp˜νHμν + 4p˜μHμ + H , p˜ = , (44) 6 Mπ where H , Hμ and Hμν are related to the sunset-type integrals of Fig. 1(d). We have not been able to find a compact representation for these integrals and only elaborate on their numerical analysis in Appendix A in some detail. 26 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 5. Numerics 5.1. Setup The numerical analysis is performed in line with the setup of Ref. [12]. The quantity of interest is − ≡ MπL Mπ RMπ , (45) Mπ whose quark mass dependence shall be evaluated numerically for different sizes L. The parame- ters of RMπ are(seeEqs.(19)–(24) and Eqs. (37)–(44)) the pion mass Mπ and the pion decay constant F in infinite volume as well as (implicitly in I (4/6))theSU(2) low energy constants. π Mπ The quark mass dependence of the pion decay constant may be taken into account by expressing Fπ as a function of the pion mass Mπ (see, e.g., [6]). Regarding the low energy constants, we use the ones determined in [9,10] which are the same as in our previous finite size studies [6,12]. 5.2. Results = We plot our results for RMπ in Figs. 4 and 5 both for L 2, 3, 4 fm as a function of Mπ and for Mπ = 100, 300, 500 MeV as a function of L. We show the one-loop result (LO) as well as the two-loop result (NLO). These shall be compared with the resummed asymptotic formula with LO/NLO/NNLO input for the ππ scattering amplitude. Note that the one-loop result and the resummed asymptotic formula to LO coincide. The best estimate for RMπ is finally obtained by adding to the asymptotic pure three-loop contribution the two-loop result (NNLO asympt. + full NLO). At NLO, the finite size effects contain low energy constants, see, e.g., diagrams (f) and (g) in Fig. 1, leading to a nonnegligible error band. We take up a point which was already observed in [5], namely the large contributions when going from LO to NLO in the asymptotic formula (dotted to thin-dash-dotted). Compared to this gap, the additional contributions from the full two-loop result (thick-dash-dotted) are very small. The two-loop and the NLO result from the asymptotic formula only drift away, when we go be- yond the region where the p-regime can be safely applied. In Table 1 we wish to underline these statements with a numerical example: we show the relative numerical impact of the LO, the pure NLO and the pure NNLO contribution for RMπ . In the second column we give the source of the effect. The fifth line, e.g., contains only the contribution of the dispersive terms Ic to RMπ , once the kinematical prefactors are properly accounted for. We observe that Ic is strongly suppressed, irrespective of the value Mπ L. Consider the column with Mπ L 1.4. Although the bulk of the subleading effects is still due to the asymptotic contributions, the terms of Σ(2) play a significant role and cannot be neglected. This behaviour was expected, since with Mπ L 1.4 we might have already crossed the border of the p-regime. As we increase this parameter to Mπ L 2, the asymptotic regime begins to set in. The contributions from the resummed asymptotic formula (2) are now dominating with respect to those from Σ . The numerical results for Mπ L 2.5 con- firm this trend. The fact that even the additional NNLO asymptotic terms (a partial three-loop result) are larger than the NLO nonasymptotic contributions is in nice agreement with the analyt- ical expectation found in Eq. (36). However, we find that the suppression seen in the numbers is stronger than what one could have expected on the basis of the latter analytical argument. Finally, the discussion of the subleading effects allows us to give a reliable estimate for the lower bound G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 27 = − = Fig. 4. RMπ MπL/Mπ 1vs.Mπ for L 2, 3, 4 fm. The result of the resummed asymptotic Lüscher formula (21) with LO/NLO/NNLO chiral input (with attribute “asympt.” for NLO and NNLO in legenda) is compared to the one-loop (LO) and two-loop (NLO) result. The best estimate for RMπ is obtained by adding the pure three loop contribution from the asymptotic formula to the two-loop result (NNLO asympt. + full NLO). The error band comes from the uncertainties in the low energy constants and is only shown for the best estimate. In the region above the Mπ L = 2 line, one is not safely in the p-regime and our results should not be trusted. of Mπ L for the p-regime, Mπ L 2: lower bound for p-regime. (46) 6. Summary (i) We have evaluated the finite volume corrections for the pion mass to two-loops within the framework of chiral perturbation theory (ChPT) in the p-regime (Mπ L 1, L>2fm, Mπ < 500 MeV). (ii) We have compared the two-loop result with the resummed version of the asymptotic for- mula of Lüscher. We have found that whenever the effects are calculated for masses and volumes such that Mπ L 1, such that one is safely within the p-regime of ChPT, the contributions which are not included in the resummed asymptotic formula are very small. The result gives us confi- dence in the claim that the resummed asymptotic formula is a convenient and efficient way to reliably calculate finite volume effects for hadronic masses [12]. (iii) The derivation of the asymptotic formula for decay constants [34] follows closely the original one for the masses [1]—the only new, subtle point, concerns the amplitude which appears in the integrand of the asymptotic formula, which has a pole in the integration region which first needs to be subtracted (see [34] for details). This is, however, a technical point—the physics 28 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 = − = Fig. 5. RMπ MπL/Mπ 1vs.L for Mπ 100, 300, 500 MeV. The rest of the legenda is as in Fig. 4. Table 1 = − = RMπ MπL/Mπ 1 for a 2 fm volume and pion masses Mπ 140, 200, 250 MeV. We show the relative numerical impact of the LO, the pure NLO and the pure NNLO contribution for RMπ . In the second column we give the source of the effect. The fifth line, e.g., contains the contribution to RMπ of the dispersive terms Ic alone, once kinematical prefactors are added accordingly. The numerical results for the Ip contributions, to LO, NLO and NNLO have already been given in [6,12] = = = RMπ Mπ 140 MeV Mπ 200 MeV Mπ 250 MeV L = 2fm Mπ L 1.4 Mπ L 2.0 Mπ L 2.5 Ip LO 0.0463 0.0199 0.0105 Ip ΔNLO 0.0291 0.0149 0.0088 Ic ΔNLO 0.0005 0.0001 0.0000 Σ(2) ΔNLO 0.0206 0.0038 0.0011 Ip ΔNNLO 0.0076 0.0053 0.0035 Total 0.1041(55) 0.0440(50) 0.0239(42) of the finite volume corrections for masses and decay constants appears to be rather similar. In view of this we believe that the present results speak also in favour of the resummed asymptotic formula for decay constants [12]. A check of this claim could be obtained by evaluating the pion decay constant to two loops. (iv) The two-loop calculation performed here allows us to better estimate the region of valid- ity of ChPT in the p-regime. Our conclusion is that in the case of the pion mass it is necessary to have Mπ L 2. Again, this serves as a guideline also for decay constants and masses of other hadrons. G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 29 Acknowledgements We thank Stephan Dürr for useful discussions and a careful reading of the manuscript. This work has been supported by the Schweizerischer Nationalfonds and partly by the EU “Euridice” program under code HPRN-CT2002-00311. Appendix A. Finite volume integrals In this appendix we give further details on the finite volume integrals which occurred in the two-loop calculation. Throughout we applied dimensional regularization, as is common in ChPT. Since the finite box breaks Lorentz invariance, tensor simplifications have to be performed with care. Consider, e.g., the integral dd q eiqnL pˆ pˆ A = (pq)ˆ 2, (A.1) μ ν μν (2π)d M2 + q2 n∈Z3 π which arises when evaluating diagram (f) in Fig. 1. Remember pˆ = (iMπ , 0). The prime in the sum denotes that the term with n = 0 is excluded. The ansatz Aμν = δμνA (A.2) with A a scalar integral which is valid in infinite volume, leads to an incorrect result in finite volume. A direct calculation of the integral with the help of Eq. (A.15) yields ∞ (M2 )1+d/2 m (n) √ pˆ pˆ A =− π d K λ n . μ ν μν d/2 d/4 d/2 π (A.3) (2πλπ ) n n=1 In four dimensions the result agrees with the derivation outlined in the next section, where one proceeds along the same lines as in Section 3.4. A.1. Tadpole In Fig. 1, the finite volume corrections of the diagrams (a)–(c), (f) and (g) factorize into one- loop integrals which are of the generic form 4 iqnL d q e P pqˆ ; q2 , (A.4) (2π)4 (M2 + q2)k n∈Z3 π with k an integer positive number and P(pq,qˆ 2) a polynomial of pqˆ and q2. It suffices to discuss the case for k = 1, since k = 2, 3,... are obtained through appropriate derivatives with respect 2 to Mπ . The only singularity of the integrand is the pole of the propagator. Therefore, the contour integration analysis applied in Section 3.4 yields the result ∞ ∞ 2 √ m(n) M − + 2 √ π dy e n(1 y )λπ P iM2 y;−M2 . 2 π π (A.5) n 8π λπ n=1 −∞ In particular, in the text we have used the dimensionless function g˜1(λπ ), 16π 2 g˜ (λ ) = g M2 , ,L , 1 π 2 1 π 0 Mπ 30 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 dd q eiqnL g M2 , 0,L = , (A.6) k π (2π)d (M2 + q2)k n∈Z3 π = 2 evaluated at d 4, and where g1(Mπ , 0,L) was introduced a long time ago by Gasser and Leutwyler [2,3], see also Eq. (14). A.2. Sunset The only real two-loop diagram is the sunset graph Fig. 1(d), and we will comment on it in some detail. It is convenient to split the finite volume integrals as in Eq. (18) only after the tensor simplifications. As alluded in the beginning of the appendix, tensor simplifications in finite volume have to be performed with care. Even though one cannot rely on Lorentz invariance, the sunset tensor integrals may still be reduced to the structures d ipx 2 {H, Hμ, Hμν}= d xe G(x) {1,i∂μ, −∂μ∂ν}G(x) , (A.7) with G(x) from Eq. (7). A rather direct way to perform these steps is to work in coordinate space and to make use of partial integrations, i.e., 1 dd xeipx ∂ G(x)∂ G(x)G(x) =− dd xeipx G(x)2 p i∂ G(x) + ∂ ∂ G(x) . μ ν 2 μ ν μ ν (A.8) Notice that the same identities in momentum space are derived with the help of translational invariance which is still respected in finite volume due to the periodic boundary conditions. In the following, we will only elaborate on the scalar integral H in more detail. We expand the integral in terms of number of finite volume propagators as motivated in Section 3.3, H = H(0) + 3H(1) + H(2), (A.9) where the first (second) addend corresponds to the pure (simple) gauge fields contribution. For H(0) we refer to [31]. Further, √ ∞ d4q eiq1 nL H(1) = m(n) J¯ (pˆ − q)2 + g M2 , 0,L J(0), (A.10) (2π)4 M2 + q2 1 π n=1 π with dd 1 1 J k2 = = J¯ k2 + J(0). (A.11) d [ 2 + 2] [ 2 + − 2] (2π) Mπ Mπ (k ) The first term of Eq. (A.10) is finite, the second carries an UV-divergence which is absorbed by a counterterm. This shows that although finite volume effects do not generate new UV-divergences, they still appear at intermediate steps of the calculation. It is a thorough check on our calculation that these nonanalytic divergences cancel. Finally, d4q d4k eiqnL eikrL 1 H(2) = . (A.12) (2π)4 (2π)4 [M2 + q2] [M2 + k2] [M2 + (pˆ − q − k)2] n,r∈Z3\0 π π π n=r G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 31 In the text, we used its dimensionless version 4 4 ˜ iqn˜ λ ikr˜ λ ˜ ˜ ˜ d q˜ d k e π e π {1; kμ; kμkν} {H ; Hμ; Hμν}= , (A.13) (2π)4 (2π)4 [1 +˜q2] [1 + k˜2] [1 + (p˜ −˜q − k)˜ 2] n,r∈Z3\0 n=r which are UV-finite and do not need to be renormalized. We only have to find a convenient representation for the numerical analysis. We shall restrict ourselves to the scalar integral in the following. Introducing a Feynman parameter by combining the second and third denominator, we find 1 ˜ − + ˜ d4q˜ d4k˜ 1 eiq(r nx)λπ inkλπ H = dx . (A.14) (2π)4 (2π)4 [1 +˜q2] [1 + (p˜ −˜q)2x(1 − x) + k˜2]2 0 n,r∈Z3\0 n=r We use Schwinger’s trick for both remaining denominators ∞ 1 − + 2 = dαe α(1 x ). (A.15) 1 + x2 0 The integrals over k˜ and q˜ are then of the Gaussian type and can be performed analytically. We are then left with three integrations over a rather lengthy expression which shall not be written down here. Despite their unhandy form, the integrations may still safely be performed numerically. The accuracy of the determination of Eq. (A.14) is not limited by the integration routine, but by the rather slow convergence of the sum in n and r for moderate λπ . Consequently, the evaluation of the sunset integrals going into Δ is restricted to three significant digits. (The last digit given for Σ(2) in Table 1 is not significant.) Note that the uncertainty of the H -type integrals is not a serious matter. Firstly, it could be lowered by brute force and secondly it only plays a (minor) role, in case when the p-regime cannot be safely applied anymore. Appendix B. Self-energy to second order: Σ(2) In this section we ask ourselves whether the self-energy to second order can be represented in a similar compact form as in the case of the self-energy to first order in Eq. (25). As will be discussed, it is indeed possible to relate the self-energy to second order to a 3π scattering amplitude in the forward scattering kinematics, as illustrated in Fig. 2(b). This has already been observed by Schenk in a related context [23]. He investigated the dynamics of pions in a cold heat bath of temperature T and examined the effective mass of the pion Mπ (T ) within the framework of ChPT. The close relation between finite temperature field theory and finite volume effects becomes apparent in the imaginary time formalism, where one treats the inverse temperature as a finite extension in the imaginary time direction. The expansion in terms of the number of finite volume propagators is then in one-to-one correspondence with the expansion in terms of the number of finite temperature propagators. In the following, Schenk’s approach shall be adapted to the finite volume scenario. We first establish the relation between Σ(2) and the three-to-three particle scattering amplitude and proceed with various remarks. 32 G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 Consider the 6-point function in d dimensions in the nonlinear sigma model parameterization Eq. (3) in the forward kinematics 3 − − − − − − dx ···dx e ip(x1 x4) ik(x2 x5) iq(x3 x6)0|Tϕ1 ϕa ϕb ϕ1 ϕa ϕb |0 1 5 x1 x2 x3 x4 x5 x6 a,b=1 Z3 = T (p,k,q), (B.1) 2 + 2 2 2 + 2 2 2 + 2 2 πππ (Mπ p ) (Mπ k ) (Mπ q ) with Z the wave function renormalization constant and ϕx ≡ ϕ(x). The amplitude Tπππ(p,k,q) ˆ2 =− 2 contains a pole at p Mπ which needs to be subtracted R(p,k,q) T (p,k,q)= Tˆ (p,k,q)+ . (B.2) πππ πππ 2 + 2 Mπ p The self-energy to second order can then be written in terms of the subtracted amplitude ˆ Tπππ(p,k,q)ˆ 4 4 ˆ 1 d q d k + T (p,k,q)ˆ 1 Σ(2) = eiqnL ikrL πππ + O , (B.3) 8 (2π)4 (2π)4 (M2 + q2)(M2 + k2) F 6 n,r π π π where the prime restricts the sum to integer vectors n and r obeying n = 0 = r, and for dia- gram (c) in Fig. 6 in addition n = r. The latter restriction avoids double counting as this term is already accounted for in a simple gauge field of Σ(1). (i) The reason for the subtraction is easily accounted for. In Fig. 6 we decompose the 3π–3π scattering amplitude in terms of 1-particle irreducible parts. The subtraction removes the con- tribution from the diagram in Fig. 6(e) which does not correspond to a 1-particle irreducible self-energy diagram, once the external legs are appropriately closed. A thorough discussion on the physical interpretation of the pole term can be found in [23]. (ii) Notice that Eq. (B.2) defines a subtracted off-shell amplitude which depends on the regu- larization scheme as well as on the parameterization of the pion fields. At the order we are work- ˆ ing, the regularization dependence is not an issue, since the subtracted amplitude Tπππ(p,k,q)ˆ is only needed at tree level. However, the dependence on the parameterization of the pion fields is of concern. While the momentum integrations in Eq. (B.3) for the diagrams Fig. 6(b) and (d) put the momenta k and q on-shell and the ambiguity due to the parameterization therefore drops out, this does not happen in the case of diagram 6(c). Note that the same parameterization ambiguity already occurred in Ic in the dispersive analysis (1) 2 =− 2 of Σ (cf. Eq. (31) where q Mπ ). In order to understand the close relation between these Fig. 6. Decomposition of the 3π–3π amplitude in terms of 1 particle irreducible parts. The characters 1, a and b on the external legs denote isospin indices. G. Colangelo, C. Haefeli / Nuclear Physics B 744 (2006) 14–33 33 two terms, we first note that only diagram Fig. 1(d) contributed to Ic. Further, the simple gauge field of Fig. 1(d) with n = 0 for one propagator can immediately be written as a contribution with two finite volume propagators: one periodifies a second propagator, however with the same n = 0 as already for the first one. Since MπL does not depend on the parameterization of the pion fields, the ambiguities of the two terms have to cancel each other. To explicitly show this is however nontrivial. (2) In summary, while our representations for Ic in Eq. (22), respectively Eq. (37) and for Σ in Eq. (B.3), respectively Eq. (23) do depend on the off-shell dependence of the scattering ampli- tudes, the sum in Eq. (19) does not. (iii) Even though the self-energy to second order can be expressed in a compact form, fur- ther simplifications (similar to those performed in Section 3.4) seem not to be straightforward. Instead, for the (numerical) integrations we had to discuss the various terms contributing to Σ(2) one by one. A general discussion of finite volume integrals has been provided in Appendix A. After treating the finite volume integrals, we end up with the basic functions given in Eq. (23). References [1] M. Luscher, Commun. Math. Phys. 104 (1986) 177. [2] J. Gasser, H. Leutwyler, Phys. Lett. B 184 (1987) 83. [3] J. Gasser, H. Leutwyler, Phys. Lett. B 188 (1987) 477. [4] J. Gasser, H. Leutwyler, Nucl. Phys. B 307 (1988) 763. [5] G. Colangelo, S. Durr, R. Sommer, Nucl. Phys. B (Proc. Suppl.) 119 (2003) 254, hep-lat/0209110. [6] G. Colangelo, S. Durr, Eur. Phys. J. C 33 (2004) 543, hep-lat/0311023. [7] J. Gasser, H. Leutwyler, Ann. Phys. 158 (1984) 142. [8] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser, M.E. Sainio, Phys. Lett. B 374 (1996) 210, hep-ph/9511397. [9] J. Bijnens, G. Colangelo, G. Ecker, J. Gasser, M.E. Sainio, Nucl. Phys. B 508 (1997) 263, hep-ph/9707291; J. Bijnens, G. Colangelo, G. Ecker, J. Gasser, M.E. Sainio, Nucl. Phys. B 517 (1998) 639, Erratum. [10] G. Colangelo, J. Gasser, H. Leutwyler, Nucl. Phys. B 603 (2001) 125, hep-ph/0103088. [11] G. Colangelo, Nucl. Phys. B (Proc. Suppl.) 140 (2005) 120, hep-lat/0409111. [12] G. Colangelo, S. Durr, C. Haefeli, Nucl. Phys. B 721 (2005) 136, hep-lat/0503014. [13] C. Haefeli, hep-lat/0509078. [14] G. Colangelo, A. Fuhrer, C. Haefeli, hep-lat/0512002. [15] S.R. Sharpe, Phys. Rev. D 46 (1992) 3146, hep-lat/9205020. [16] D. Becirevic, G. Villadoro, Phys. Rev. D 69 (2004) 054010, hep-lat/0311028. [17] QCDSF-UKQCD Collaboration, A. Ali Khan, et al., Nucl. Phys. B 689 (2004) 175, hep-lat/0312030. [18] D. Arndt, C.J.D. Lin, Phys. Rev. D 70 (2004) 014503, hep-lat/0403012. [19] S.R. Beane, Phys. Rev. D 70 (2004) 034507, hep-lat/0403015. [20] S.R. Beane, M.J. Savage, Phys. Rev. D 70 (2004) 074029, hep-ph/0404131. [21] P.F. Bedaque, I. Sato, A. Walker-Loud, hep-lat/0601033. [22] J. Bijnens, N. Danielsson, K. Ghorbani, T. Lahde, hep-lat/0509042. [23] A. Schenk, Phys. Rev. D 47 (1993) 5138. [24] D. Toublan, Phys. Rev. D 56 (1997) 5629, hep-ph/9706273. [25] M. Luscher, DESY 83/116, Lecture given at Cargese Summer Institute, Cargese, France, 1–15 September 1983. [26] Y. Koma, M. Koma, Nucl. Phys. B 713 (2005) 575, hep-lat/0406034; Y. Koma, M. Koma, hep-lat/0504009. [27] J. Braun, B. Klein, H.J. Pirner, Phys. Rev. D 71 (2005) 014032, hep-ph/0408116. [28] B. Borasoy, G.M. von Hippel, H. Krebs, R. Lewis, hep-lat/0509007. [29] J. Bijnens, G. Colangelo, G. Ecker, JHEP 9902 (1999) 020, hep-ph/9902437. [30] J. Bijnens, G. Colangelo, G. Ecker, Ann. Phys. 280 (2000) 100, hep-ph/9907333. [31] J. Gasser, M.E. Sainio, Eur. Phys. J. C 6 (1999) 297, hep-ph/9803251. [32] P. Hasenfratz, H. Leutwyler, Nucl. Phys. B 343 (1990) 241. [33] U. Burgi, Nucl. Phys. B 479 (1996) 392, hep-ph/9602429. [34] G. Colangelo, C. Haefeli, Phys. Lett. B 590 (2004) 258, hep-lat/0403025. Nuclear Physics B 744 (2006) 34–58 Pseudoscalar meson decay constants and couplings, the Witten–Veneziano formula beyond large Nc, and the topological susceptibility ✩ G.M. Shore Department of Physics, University of Wales, Swansea, Swansea SA2 8PP, UK Received 13 January 2006; received in revised form 9 March 2006; accepted 10 March 2006 Available online 31 March 2006 Abstract The QCD formulae for the radiative decays η,η → γγ, and the corresponding Dashen–Gell-Mann– Oakes–Renner relations, differ from conventional PCAC results due to the gluonic U(1)A axial anomaly. This introduces a critical dependence on the gluon topological susceptibility. In this paper, we revisit our earlier theoretical analysis of radiative pseudoscalar decays and the DGMOR relations and extract explicit experimental values for the decay constants. This is our main result. The flavour singlet DGMOR relation is the generalisation of the Witten–Veneziano formula beyond large Nc, so we are able to give a quantitative assessment of the realisation of the 1/Nc expansion in the U(1)A sector of QCD. Applications to other aspects of η physics, including the relation with the first moment sum rule for the polarised photon structure γ function g , are highlighted. The U(1)A Goldberger–Treiman relation is extended to accommodate SU(3) 1 flavour breaking and the implications of a more precise measurement of the η and η -nucleon couplings are discussed. A comparison with the existing literature on pseudoscalar meson decay constants using large-Nc chiral Lagrangians is also made. © 2006 Elsevier B.V. All rights reserved. 1. Introduction The phenomenology of the pseudoscalar mesons opens a window on many interesting aspects of the non-perturbative dynamics of QCD, including spontaneous chiral symmetry breaking, the electromagnetic and gluonic axial anomalies, the OZI rule, the gluon topological susceptibility, ✩ This research is supported in part by PPARC grants PP/G/O/2002/00470 and PP/D507407/1. E-mail address: [email protected] (G.M. Shore). 0550-3213/$ – see front matter © 2006 Elsevier B.V. All rights reserved. doi:10.1016/j.nuclphysb.2006.03.011 G.M. Shore / Nuclear Physics B 744 (2006) 34–58 35 and so on. Nevertheless, until comparatively recently, phenomenological analyses did not take fully into account the role of the gluonic U(1)A anomaly in this sector, with the result that most of the existing determinations of the pseudoscalar meson decay constants and couplings are based on over-simplified theoretical formulae which miss the most interesting anomaly-sensitive physics. In a previous paper [1] (see also [2]) we have analysed in detail the radiative decays of the neutral pseudoscalar mesons, π 0,η,η → γγ, and provided a set of formulae describing these processes together with modified Dashen–Gell-Mann–Oakes–Renner (DGMOR) [3,4] relations which fully include the effect of the anomaly and the gluon topological susceptibility in the flavour singlet sector as well as explicit flavour SU(3) breaking. In this paper, we confront these formulae with experimental data to extract a set of results for the four pseudoscalar meson decay constants, f 0η ,f0η,f8η and f 8η. Our approach, which is based on a straightforward generalisation of conventional PCAC to include the anomaly in a renormalisation group consistent way, is unique in that it is theoretically consistent yet does rely on using the 1/Nc expansion. The decay formulae and DGMOR relations we derive are valid for all Nc. At some point, however, the fact that the η is not a Nambu– Goldstone (NG) boson means that to make predictions from our formulae we need to augment the standard dynamical approximations of PCAC or chiral perturbation theory with a further dynamical input. Here, we make a judicious use of 1/Nc ideas to justify the use of the lattice calculation of the topological susceptibility in pure Yang–Mills theory as an input into our flavour singlet DGMOR relation. The self-consistency of this approach allows us to test the validity of 1/Nc methods in the U(1)A sector against real experimental data. In particular, the flavour singlet DGMOR relation is the finite-Nc generalisation of the well-known Witten–Veneziano formula [5,6], and we can use our results to study how well the large-Nc limit is realised in real QCD. This is especially important in view of the extensive use of large Nc in modern duality-based approaches to non-perturbative gauge theories such as the AdS–CFT correspondence. The relation of our analysis of the low-energy decays η (η) → γγ to high-energy physics in γ the form of the first moment sum rule for the polarised photon structure function g1 [7,8] is also discussed. The dependence of this sum rule on the photon momentum encodes many aspects of non-perturbative QCD physics and its full measurement is now within the reach of planned high-luminosity e+e− colliders [9]. Our methods can also be readily extended to a variety of other reactions involving the pseudoscalar mesons, including η (η) → Vγ, where V is one of the flavour singlet vector mesons ρ,ω,φ, η (η) → π +π −γ , a variety of electro and photoproduction reactions such as γp → pη (η) and γp → pφ, and η and η production in pp collisions, pp → ppη (η). The latter are relevant for determining the couplings gηNN and gη NN which occur in gen- eralisations of the Goldberger–Treiman relation [10]. In this paper, we present a new version of the U(1)A Goldberger–Treiman relation [11–13] incorporating flavour octet–singlet mixing in the η–η sector. We are then able to use our results for the decay constants to show how the physical interpretation of this relation, which is closely related to the ‘proton spin’ problem in p high-energy QCD (i.e., the first moment sum rule for the polarised proton structure function g1 ), depends critically on the experimental determination of gη NN. Finally, in Appendix A, we give a brief comparison of the similarities and differences be- tween our approach and an alternative theoretically consistent framework for describing η and η physics, viz. the large-Nc chiral Lagrangian formalism of Refs. [14–16] (see also [17–21]). A comprehensive review of the phenomenology of the pseudoscalar mesons, incorporating the chiral Lagrangian results, is given in Ref. [22] 36 G.M. Shore / Nuclear Physics B 744 (2006) 34–58 2. Radiative decay and DGMOR formulae In this section, we present the radiative decay and DGMOR formulae and discuss how to extract phenomenological quantities from them, emphasising their 1/Nc dependence. A brief review of their derivation in the language of ‘U(1)A PCAC’ is given in Section 3, but we refer to the original papers [1,2,23] for full details, including a more precise formulation using functional methods. Throughout, we consider the case of three light flavours, nf = 3. The assumptions made in deriving these formulae are the standard ones of PCAC. The for- mulae are based on the zero-momentum chiral Ward identities. It is then assumed that the decay aα 2 2 ‘constants’ (f (k ) in the notation below) and couplings (gηαγγ(k )) are approximately con- stant functions of momentum in the range from zero, where the Ward identities are applied, to the appropriate physical particle mass. It is important that this approximation is applied only to pole- free quantities which have only an implicit dependence on the quark masses. This is equivalent to assuming pole dominance of the propagators for the relevant operators by the (pseudo-)NG bosons, and naturally becomes exact in the chiral limit. Notice that these dynamical approxima- tions are equally necessary in the chiral Lagrangian approach, where they are built in to the basic structure of the model—the effective fields are chosen to be in one-to-one correspondence with the NG bosons and pole dominance is implemented through the low-momentum expansion with momentum-independent parameters corresponding to the decay constants. The relations for the π0 decouple from those for the η and η in the realistic approximation of SU(2) flavour invariance of the quark condensates, i.e., ¯uu=dd¯ . The decay formula for π0 → γγ is the standard one, α f g = a3 em (2.1) π πγγ em π 3 = 1 with aem 3 Nc, while the DGMOR formula is simply 2 2 =− ¯ + ¯ fπ mπ mu uu md dd . (2.2) The coupling gπγγ is defined as usual from the decay amplitude: 0 =− α β λ ρ γγ π igπγγ λραβ p1 p2 (p1) (p2) (2.3) in obvious notation [1]. The r.h.s. of Eq. (2.1) arises from the electromagnetic U(1)A anom- aly and is of course sensitive to the number of colours, Nc. The decay constant has a sim- ple theoretical interpretation as the coupling of the pion to the flavour triplet axial current, | 3 | = 0 Jμ5 π ikμfπ , which is conserved in the chiral limit. In the η,η sector, however, explicit SU(3) flavour breaking means there is mixing. The decay constants therefore form a 2 × 2 matrix: 0η 0η aα f f f = , (2.4) f 8η f 8η where the index a is an SU(3) × U(1) flavour index (including the singlet a = 0) and α denotes the physical particle states π 0,η,η . These four decay constants are independent [14,15,24,25], in contradiction to the original phenomenological parametrisations which erroneously expressed f aα as a diagonal decay constant matrix times an orthogonal η–η mixing matrix, giving a total of only three independent parameters. Sometimes the four decay constants are expressed in terms of two constants and two mixing angles but, while we also express our results in this form in Section 4, there is no particular reason to parametrise in this way. G.M. Shore / Nuclear Physics B 744 (2006) 34–58 37