Area-scaling of quantum fluctuations Amos Yarom, Ram Brustein Department of Physics, Ben-Gurion University, Beer-Sheva 84105, Israel E-mail: [email protected], [email protected], We show that fluctuations of bulk operators that are restricted to some region of space scale as the surface area of the region, independently of its geometry. Specif- 4 0 ically, we consider two point functions of operators that are integrals over local 0 2 operator densities whose two point functions falls off rapidly at large distances, and n a does not diverge too strongly at short distances. We show that the two point func- J 3 tion of such bulk operators is proportional to the area of the common boundary of 1 the two spatial regions. Consequences of this, relevant to the holographic principle 1 v and to area-scaling of Unruh radiation are briefly discussed. 1 8 0 1 I. INTRODUCTION 0 4 0 / The discovery of the non extensive nature of black hole entropy [1, 2] has lead to entropy h t - bounds on matter, and to the proposal of the holographic principle (see [3] for a review) – a p e conjecture regarding the reduction of the number of degrees of freedom needed to describe h : v a theory of quantum gravity. Later, it was discovered that other thermodynamic quantities i X of fields in a black-hole background also scale as the surface area of its horizon, and con- r a sequently, it was hypothesized that it is the entangled nature of the quantum state of the system inside and outside the horizon, which leads to area-scaling. This view has received some support from the area-scaling properties of entanglement entropy in flat space [4, 5, 6]. In this paper we study area-scaling of two point functions of a certain class of quantum field theory bulk operators in Minkowski space, to be defined shortly. We shall show, under weak assumptions, that quantum expectation values of operators restricted to a region of space will scale as its surface area, regardless of the region’s geometry. This area-scaling property may then be used to establish area-scaling of thermodynamic quantities [7], and in some sense, a bulk-boundary correspondence [8]. The two point functions we wish to consider are the expectation values of a product of two operators which are restricted to some regions V and V of Minkowski space. In 1 2 2 order to restrict operators to a region V , we use operator densities (~x), and define OVi = i Oj j (~x)ddx. ViOj R We show that if the connected two point function of the operator densities F( ~x ~y ), | − | (~x) (~y) = F( ~x ~y ) 2g satisfies the following conditions: i j C hO O i | − | ≡ ∇ (I) g(ξ) is short range: at large distances it behaves like g(ξ) 1/ξa, with a d 1, ∼ ≥ − (II) g(ξ) is not too singular at short distances. Explicitly, we require that for small ξ, g(ξ) 1/ξa with a < d 2, ∼ − and if V and V are finite, then the connected two point function OV1OV2 is proportional 1 2 h i j iC to the area of the common boundary of the two regions V and V : OV1OV2 S(B(V ) 1 2 h i j iC ∝ 1 ∩ B(V )). Here B(V) is the boundary of V, and S(B(V)) is its area. 2 This implies that the fluctuations (or variance: var(OV) = (OV2 OV 2 ) of the op- i h i − h i i i erator OV scale as the surface area. In particular, the energy fluctuations inside V will be i proportional to S(B(V)). This is discussed at length in [7]. In section II we give a detailed proof of the area-scaling property of two point functions of bulk operators satisfying conditions (I) and (II). Section III contains an explicit calculation of energy fluctuations, and fluctuations of the boost operator. We discuss these results in section IV. II. AREA-SCALING OF TWO POINT FUNCTIONS We shall first give a general explanation as to why conditions (I) and (II) are required for area-scaling, and outline the proof for area-scaling of two point functions. A more detailed discussion will then follow. In order to evaluate OV1OV2 , we may express it as follows, h i j iC OV1OV2 = F ( ~x ~y )ddxddy h i j iC i,j | − | ZV ZV = D(ξ)F (ξ)dξ. (1) i,j Z The integral has been factored into a product of a purely geometric term D(ξ), and an operator dependent term F (ξ). Using 2g (ξ) F (ξ), we may integrate eq.(1) by i,j i,j i,j ∇ ≡ parts. The surface term then vanishes due to conditions (I) and (II), and we get d 1 d OV1OV2 = D(ξ) ξd−1 g(ξ)dξ. (2) h i j iC − dξ ξd−1 dξ Z (cid:18) (cid:19) 3 FIG. 1: A general case when the regions have some overlap with each other, and neither is fully contained in the other. The regions have common ‘boundaries from within’, and ‘boundaries from the outside’. To proceed we note that the geometric term is of the form: D (ξ) = δ(ξ ~r ~r )ddr ddr . (3) V1,V2 −| 1 − 2| 1 2 ZV1ZV2 We show in the next subsection that D = G Vξd−1+G Sξd+ (ξd+1), with V a volume V1,V2 V S O term, S a surface area term, and G numerical coefficients. There are some geometries for which this scaling is more obvious: for example, if V and V are disjoint volumes with a 1 2 common boundary which may be approximated as flat, then one may carry out the integral in eq.(3) by switching to ‘center of mass coordinates’ R~ = ~r1+~r2, and~r =~r ~r . Restricting 2 2− 1 the value of ~r to be equal to ξ, results in confining R~ to be, at most, a distance of ξ/2 | | ~ from the boundary. The integral over the R coordinate will give a term proportional to Sξ, S being the area of the boundary. The remaining integral over the ~r coordinate will give a term proportional to ξd−1. Combining these results, we get the claimed form of D (the V1,V2 volume term vanishes in this case). Another possibility is that V = V . In this case an integral over the center of mass 1 2 coordinate R~ will yield, to leading order, a volume dependent term. The remaining integral over the ~r coordinate will give the appropriate powers of ξ. Using these two simple cases, one can now construct D for general geometries by V1,V2 dividing one of the volumes, say V , into sub-volumes which are either contained in V , just 1 2 touching V , or disconnected from V . In order for the appropriate form of D to follow, 2 2 V1,V2 4 one would need that the coefficient G be exactly the same for both cases described above. S This requires a more detailed calculation and is shown in the following subsection. Going back to eq.(2), we observe that the integral comes from the region ξ 0, and that → the contribution of the volume dependent term of D will vanish ∂ D(ξ) S + (ξ). V1,V2 ∂ξ ξd−1 ∼ O The leading contribution to the integral will come from the surface term(cid:16), wh(cid:17)ich gives area- scaling behavior of two point functions, as claimed. The vanishing of the contribution of the volume term is a result of properties (I) and (II) of g, which make the surface term vanish when integrating eq.(1) by parts, and of the special polynomial dependence of D (ξ) on V1,V2 ξ. A. The geometric term. As we have shown schematically, the area-scaling properties of correlations depend on the properties of the geometric term defined in (eq.(3)). In this subsection we shall study it in some detail. First we note that there exists a ξ and ξ , such that D(ξ) = 0 for ξ ξ or min max max ≥ ξ ξ : define ξ = inf ~x ~y ~x V ,~y V . For ξ < ξ there are no values of ~x min min 1 2 min ≤ {| − | ∈ ∈ } and ~y which will have a non zero c(cid:12)ontribution to the integral, and therefore D(ξ) = 0 for (cid:12) this region. Similarly ξ = sup ~x ~y ~x V ,~y V . max 1 2 {| − | ∈ ∈ } For the cases where ξ = 0 (that is, (cid:12)the closed sets V and V are not disjoint), we shall min 1 2 (cid:12) show that D (ξ) satisfies V1,V2 D (ξ) = G Vξd−1 +G (S(B ) S(B ))ξd + (ξd+1). V1,V2 V S in − out O G and G are constants which depend only on the dimensionality of space, explicitly V S G = d πd/2 . B is the common boundary of V and V . B is the part of the S (d−1)Γ(d/2+1) in/out 1 2 out boundary when V and V are on opposite sides of the boundary, whereas B is the part 1 2 in of the boundary when V and V are on the same side of the boundary. The first term is a 1 2 volume dependent term. To solve the integral (3) we switch to a ‘center of mass’ coordinate system (shown in Fig. 2): R~ = ~r1+~r2, and ~r = ~r ~r , so that ~r points from ~r to ~r , and R~ points to the 2 2 − 1 1 2 middle of ~r (since R~ =~r +1/2~r). In this new coordinate system, 1 D (ξ) = δ(ξ r)ddRddr. (4) V1,V2 − Z Z 5 V1 V2 r r1 R r 2 FIG. 2: A pictorial description of the (~r ,~r ) coordinates, and the (R~,~r) coordinates. 1 2 In order to carry out the integration, we wish to use a generalized ‘radial’ coordinate ρ such that ρ = ρ will define the boundary B(V ) of V , and generalized ‘angular’ coordinates 0 1 1 α which will define solid angles on the boundary. To define such a coordinate system we i foliate space into surfaces which, when very close to B(V ), are parallel to it. ρ is chosen to 1 be the coordinate which points to different leaves of the foliation. ρ = ρ , as stated, defines 0 the surface B(V ). We also choose ρ such that for a point R~ for which ρ ρ is small 1 0 | − | enough, then ρ ρ will specify the distance of R~ from the boundary. α are generalized 0 i | − | angles on each hyper-surface. The unit volume in this coordinate system is ddR = J(ρ,α )dρ dα . i i i Y For the vector ~r we choose a polar coordinate system: ddr = rd−1dΩ = rd−1sind−2θdθdΩ . ⊥ ~ For a given point R, the integration over the angular coordinates of ~r will give us the solid angle subtended by all allowed values of ~r ~r , when it is centered at R~. In preparation 2 1 − for performing the integral in eq.(4) for a general geometry (as shown in Fig. 1) we first consider several particular cases. Case 1: Regions which have a common boundary, but no common interior. Here we expect to have D (ξ) = G S(B(V ) B(V ))ξd + (ξd+1). V1,V2 S 1 ∩ 2 O 6 ξ V1 V 2 FIG. 3: Region of integration of the R~ coordinate: The boundaries of regions 1 and 2 are given by the thick lines. The allowed region of R~ is given by the thin lines. Fixing ξ to have some small value, we look at all values ~r is allowed to have for a fixed value of R~. Defining ξ = r~ r~ , then since ~r and ~r are on different sides of the boundary 1 2 1 2 | − | B B(V ) B(V ), R~ is restricted to a distance of ξ/2 from B (See Fig.3). Therefore 1 2 ≡ ∩ ρ0+ξ/2 ddR = J(ρ,α )dρ dα . (5) i i Z Zρ0−ξ/2 ZA(ρ) i Y In this case A(ρ) is the region of angular integration for each leaf of the foliation. We denote the range of angles which define the common boundary as A , so that B B = (ρ ,α ) α A . 0 i i B { | ∈ } We note that with this definition we may also write J(ρ ,α ) dα = S(B)+ (ξ). 0 i i O ZA(ρ0) i Y For a certain point R~ = (ρ,α ) close to the boundary: ρ ρ < ξ/2 (shown in Fig. 4), i 0 | − | ¯ the angular integration over ~r will be restricted to a region Ω(ρ,α ;ξ). i Therefore we have: ρ0+ξ/2 D (ξ) = J(ρ,α ) ξd−1dΩ dα dρ. (6) V1,V2 i i Zρ0−ξ/2 ZA(ρ) ZΩ¯(ρ,αi;ξ) i Y 7 ξ V1 V 2 r R FIG. 4: The allowed region for ~r in the case where the volumes are ‘just touching’. The dark end of the circle shows the region where ~r is allowed to point to. Changing variables to the dimensionless variable ζ = (ρ ρ )/(ξ/2), and integrating over 0 − the angular coordinate Ω we obtain, 1 ξ D (ξ) = ξd−1 J(ρ +ξ/2ζ,α )Ω¯(ρ +ξ/2ζ,α ;ξ) dα dζ. (7) V1,V2 0 i 0 i i2 Z−1ZA(ρ0+ξ/2ζ) i Y This is an exact expression. Expanding the integrand on the r.h.s. of eq.(7) in powers of ξ, we obtain J(ρ +ξ/2ζ,α )Ω¯(ρ +ξ/2ζ,α ;ξ) dα 0 i 0 i i ZA(ρ0+ξ/2ζ) i Y = J(ρ +ξ/2ζ,α )Ω¯(ρ +ξ/2ζ,α ;ξ) dα + (ξ) 0 i 0 i i ZA(ρ0+ξ/2ζ) Yi (cid:12)ξ=0 O (cid:12) = J(ρ ,α )Ω(ζ) dα + (ξ), (cid:12) (8) 0 i i O ZA(ρ0) i Y where ¯ Ω(ζ) = limΩ(ρ +ξ/2ζ,α ;ξ). 0 i ξ→0 We will show that Ω is a function of ζ only, and obtain an explicit expression for it. For very small ξ, we look at a point R~, a distance ρ ρ < ξ/2 from the boundary, and 0 | − | calculate the solid angle subtended by a vector ~r centered at R~ whose one end is in V , and 1 8 θ0 ρ θ0 ρ0 ρ0 ρ FIG. 5: The allowed region for ~r (dark region) in the case where ξ is very small, so that the boundary may be approximated as flat. ρ specifies the location of the common boundary, and ρ 0 the leaf of the foliation to which R~ is pointing to. The left diagram is for the case where R~ V , 1 ∈ and the right diagram for R~ V . 2 ∈ other end is in V . Since we are taking the limit where ξ 0, the shape of the boundary 2 → close to the point R~ may be considered flat (shown in Fig. 5). Therefore Ω is the solid angle generated by ~r when very close to a flat surface. Defining the z-axis as the axis perpendicular to the surface, we get that the angle θ between the z-axis and ~r can range over values from 0 to θ , with cos(θ ) = (ρ ρ)/(ξ/2) 0 0 0 − for ρ ρ and cos(θ ) = (ρ ρ )/(ξ/2) for ρ ρ . Note that in any case 0 θ π/2, 0 0 0 0 0 ≤ − ≥ ≤ ≤ implying that cosθ 0. 0 ≥ Therefore, for d > 2, θ0 Ω = sind−2θdθ dΩ . ⊥ Z0 Z We note that θ0 1 1 d 1 1 d 1 sind−2θdθ = B , − B , − (9) ζ2 2 2 2 − 2 2 Z0 (cid:18) (cid:18) (cid:19) (cid:18) (cid:19)(cid:19) where B (a,b) is the partial beta function. The full solid angle of a d-dimensional sphere is x given by π C = sind−2θdθ dΩ d ⊥ Z0 Z 1 d 1 = B , − dΩ . (10) ⊥ 2 2 (cid:18) (cid:19)Z Putting eq.(10) into eq.(8), we get J(ρ ,α )Cd 1 Bζ2(12, d−21) dα + (ξ). 0 i 2 − B(1, d−1) i O ZA(ρ0) 2 2 ! i Y 9 Integrating over the angular coordinates we are left with Cd 1 Bζ2(21, d−21) S(B)+ (ξ), 2 − B(1, d−1) ! O 2 2 and plugging this into eq.(7), will give us: D (ξ) = ξξd−1 Cd 1 1 Bζ2(12, d−21) S(B)dζ + (ξ) V1,V2 2 2 − B(1, d−1) O " Z−1 2 2 ! # = ξdC S(B) 1 1 Bζ2(21, d−21) dζ + (ξd+1). (11) 2 d − B(1, d−1) O Z0 2 2 ! So now we have to evaluate the integral 1Bζ2 1, d−1 dζ = 1 ζ2t−21(1 t)d−23dtdζ, 2 2 − Z0 (cid:18) (cid:19) Z0 Z0 which can be done by changing the order of integration. The result is that C B(1, d−1) D (ξ) = ξd dS(B) 2 + (ξd+1). V1,V2 2 B(1, d−1) O 2 2 Expressing the Beta function as a product of Gamma functions, and using the explicit expression C = dπd/2 , we get the final result d Γ(d/2+1) dπd/2 D (ξ) = S(B)ξd + (ξ). (12) V1,V2 (d 1)Γ(d/2+1) O − Case 2: Next we consider geometries where V = V V, for which we expect to have: 1 2 ≡ D (ξ) = G Vξd−1 G S(B(V))ξd + (ξd+1). V1,V2 V − S O We consider the same foliation of space as before. In this case, even at small ξ, we integrate over points inside V. The angular integration over ~r will be constrained only when the region where R~ is a distance of ξ/2 from the boundary. Therefore: ρ0 D (ξ) = ξd−1J(ρ,α )dΩ dα dρ (13) V1,V2 i i Z0 ZA(ρ)ZΩ¯ i Y ρ0−ξ/2 = ξd−1J(ρ,α )C dα dρ i d i Z0 ZA(ρ) i Y ρ0 + ξd−1J(ρ,α )dΩ dα dρ. i i Zρ0−ξ/2ZA(ρ)ZΩ¯ i Y 10 The first expression is the unconstrained part, and the second expression is the constrained part. Considering first the integral for which ~r is not constrained, we have ρ0−ξ/2 ξd−1J(ρ,α )C dα dρ i d i Z0 ZA(ρ) i Y ρ0 π = VC ξd−1 ξd−1J(ρ,α ) sind−2θdθdΩ dα dρ. d i ⊥ i − Zρ0−ξ/2ZA(ρ) Z0 i Y Proceeding as before, we change variables of integration to ζ = (ρ ρ )/(ξ/2), and expand 0 − in small ξ. ξd 0 π VC ξd−1 S(B) sind−2θdθdΩ dζ + (ξd+1). (14) d ⊥ − 2 O Z−1Z0 Under the same approximation, the constrained part reduces to ξd 0 S(B) Ω(ζ)dζ + (ξd+1). (15) 2 O Z−1 In this case, θ, the angle between ~r and the z-axis, is restricted to θ , and π θ where 0 0 − cosθ = ζ. Hence 0 − π−θ0 Ω(ζ) = sind−2θdθ. Zθ0 Combining the contributions of the surface terms of (14) and (15) to those of the surface term of (13), we have ξd 0 θ π S(B) + sind−2θdθ Ω dζ. ⊥ − 2 Z−1(cid:18)Z0 Zπ−θ(cid:19) Z Since 0 θ π/2, both angular integrals are equal, and the above equation simplifies to ≤ ≤ ξd 0 θ S(B) 2 sind−2θdθ Ω dζ. ⊥ − 2 Z−1 Z0 Z Using eq.(9) we may carry out the integral over the θ coordinate: ξd 0 1 1 d 1 1 d 1 S(B) 2 B , − B , − dζ ζ2 − 2 2 2 2 − 2 2 Z−1 (cid:18) (cid:18) (cid:19) (cid:18) (cid:19)(cid:19) ξd 1 1 1 d 1 1 d 1 = S(B) B , − B , − dζ. ζ2 − 2 2 2 2 − 2 2 Z−1 (cid:18) (cid:18) (cid:19) (cid:18) (cid:19)(cid:19) Comparing this with eq.(11), we find that the surface term contribution has exactly the same magnitude as the leading order contribution to D (ξ) in Case 1, its sign, however, V1,V2