Boundary Anomalies and Correlation Functions

It was shown recently that boundary terms of conformal anomalies recover the universal contribution to the entanglement entropy and also play an important role in the boundary monotonicity theorem of odd-dimensional quantum field theories. Motivated by these results, we investigate relationships between boundary anomalies and the stress tensor correlation functions in conformal field theories. In particular, we focus on how the conformal Ward identity and the renormalization group equation are modified by boundary central charges. Renormalized stress tensors induced by boundary Weyl invariants are also discussed, with examples in spherical and cylindrical geometries.


Introduction
As a non-trivial extension of Poincaré symmetry, Weyl invariance imposes significant constraints on the structure of correlation functions. For a generally non-conformal quantum field theory (QFT), the correlation functions approach those of a conformal field theory (CFT) sitting at one end point of the renormalization group (RG) flow at short distance, and those of a CFT sitting at another end point of the RG flow at long distance. It is therefore of great importance to search for a general principle constraining the flows between the two end points. In d = 2, the well-known c-theorem was proved in [1] thirty years ago. Conjectured first by [2], a proof of the a-theorem in d = 4 was given only recently by [3] [4], using the so-called dilaton anomaly effective action. The proof has not been found for d = 6 QFTs. See [5] [6] [7] [8] for the recent progress. In proving these monotonicity theorems, the central charges, defined as the coefficients of the trace anomaly when embedding the theory in a curved background, play the central role. In particular, it is the central charge of the topological Euler density that satisfies the irreversibility of the RG flow.
There is no Weyl anomaly in odd-dimensional (compact) manifolds, and hence a definition for central charge becomes elusive. An alternative candidate that satisfies the monotonic behaviour along the RG flow in d = 3 was suggested by [9] as the Euclidean path integral of the CFT conformally mapped to S 3 . See [10] [11] for further discussions on such an F-theorem; [12] points out an interpolation between the a-theorem and the F-theorem at the fixed points. The universal part of the partition function on S 3 can be further identified as the constant piece of the vacuum entanglement entropy (EE) across a disk. Moreover, with the help of strong subadditivity [13], the irreversibility of the RG flow in d = 3 can be proved [14]. It remains an important question whether one can link the fundamental properties of EE to the monotonicity theorem for spacetime dimensions higher than three. (See [15] [16] for examining the general RG flow using holography.) However, in many physical systems, in particular in condensed matter physics, edge effects are important. It has also been shown recently that the universal part of EE can be understood purely as a boundary effect [17] (so in some sense the "area/boundary law" of EE is extended to also include the UV cut-off-independent log-term), which solves the main puzzle left from an earlier attempt [18]. The study of EE in d = 3 with a spacetime boundary is discussed in a more recent work [19]. Most important for this paper is the fact that there are additional boundary invariants in the presence of a boundary, for any-dimensional CFTs. To make a connection with the monotonicity theorem, we notice that the boundary central charge orders the boundary RG flow in d = 3 [20]; this b-theorem is then a generalization of the boundary g-theorem [21]. See [22] for a related discussion.
Motivated by these results, we here consider how the boundary central charges affect the conformal Ward identities and the RG equation of stress tensor correlation functions. We will largely focus on d = 4, but will also discuss d = 3 where only boundary anomalies exist. (In d = 2, there is no new boundary anomaly and only a boundary term of the Euler density is needed, which we discuss briefly in the appendix.) To set the stage for our discussion, in the next section we first review the conformal invariance and the correlation function. The explicit expression for the stress tensor three-point function is rather bulky and we will refer the reader to [23] [24]. Our main focus here instead is the additional contribution to correlation functions from the local counter-terms (conformal anomaly), in particular when a manifold has a boundary. In Sec. 3, we revisit several main identities of correlation functions. We are interested in correlation functions in flat bulk spacetime but we allow a generally curved, codimension-1 boundary. The boundary will be assumed to be compact and smooth (no corners). When revisiting these identities, we will not set the Weyl anomaly to be zero in the flat limit since the boundary terms survive even in the flat space. The boundary local counter-terms also contribute non-trivially to the stress tensor in the flat limit, so we will not drop the stress tensor either. These identities then generalize the ones given in the literature (for example, [23] [24]).
The d = 4 CFTs will be considered in Sec. 4. We first discuss the RG equation of the threepoint function and the conformal Ward identity for a compact manifold. Then, we generalize these results by introducing a boundary. It is found that the RG equation is modified by a special boundary central charge defined by a Weyl invariant constructed solely from the extrinsic curvature. The Ward identity is also modified due to the boundary counter-terms. Moreover, we obtain stress tensors in the vicinity of a boundary whose values are determined by boundary central charges. Examples in ball and cylindrical geometries are given. A similar analysis for d = 3 CFTs is in Sec. 5, where the story becomes simpler because there is no bulk conformal anomaly and there are only two simple Weyl anomalies living on the boundary.
Among our new results, the most important ones are: a formal expression of the trace conformal Ward identity for general d (17); an RG equation for the three-point function (65), the Ward identity (66), a b 1 -type stress tensor for a 4-cylinder (82), a-type stress tensors for a 4-ball (88) and for a 4-cylinder (91) for d = 4; the Ward identity (97), an RG equation for the two-point function (98), a c-type stress tensor for a 3-cylinder (104), a-type stress tensors for a 3-ball (108) and for a 3-cylinder (111) for d = 3.
In the conclusion we point out related questions. The appendix contains useful formulae for metric variation, and also a brief discussion on d = 2 CFTs with a boundary.

Conformal Invariance and Anomalous Terms
For a conformal transformation g, x µ → x ′ µ = (gx) µ , we can define a local orthogonal matrix on where Ω g (x) is the scale factor showing up in the transformed line element ds ′2 = Ω g (x) −2 ds 2 . In the Euclidean signature we replace η µν by δ µν on R d and R µν (x) becomes a local rotation matrix belonging to the group O(d). One can generate a conformal transformation by a combination of translations, rotations and inversions, although the conformal group SO(d + 1, 1) can only be formed by even numbers of inversions since the inversion itself is not connected to the identity. If we consider an inversion through the origin, The matrix I µν (x − y) transforms like a vector and it can be regarded as a parallel transport for the conformal transformations. For three points, x, y, z, one can also define a covariant vector that transforms homogeneously. At the point z such a vector is given by Similarly, at points x, y one can define covariant vectors X µ and Y ν via cyclic permutation. Denote a conformal primary operator as O i (x) with i representing components in some representation of the rotation group. To construct the correlation functions, it is useful to adopt the induced representations [25] to write the conformal transformation as where ∆ is the conformal dimension; D i j (R(x)) is the matrix in the associated representation. The conjugate representation of (4) for the conjugate field. The two-point function for conformal primary operators can be written by in the irreducible representations of O(d). The overall constant C O might be set to one by operator redefinition.
Our main subject of interest is the correlation functions of the stress tensor T µν , which is a symmetric tensor satisfying ∂ µ T µν = 0 and T µ µ = 0. The latter property implies conformal invariance at the classical level. T µν has the scale dimension d required also by the conformal invariance. Now applying (5) to the stress tensor gives where C T is a constant coefficient and Here L λρ αβ is a projection operator onto the space of a symmetric, traceless tensor. The expression of the three-point function becomes rather involved [23] [24]. To keep the expression simple we will not list its full expression but simply write it as The structure of Γ µναβσρ (x, y, z) is determined by the conformal symmetry and the stress tensor conservation. Its explicit form will not be too relevant in the present discussion. However, it is important to know that there are three independent forms (in d = 3 there are two and only one in d = 2) for the stress tensor three-point function [23].
Here we focus on the additional contribution to (8) coming from possible local counter-terms in the effective action. These counter-terms give precisely the Weyl anomaly. We will adopt dimensional regularization working in d → d + ǫ dimensions and refer to these counter-terms as the anomaly effective action denoted by W ≡ µ ǫ ǫW . By including possible counter-terms we have In d = 4, there are two types of local counter-terms: one type-A anomaly and one type-B anomaly. (We will remove the scheme-dependent type-D total derivative anomaly.) The contributions from the bulk anomalies were first calculated by [23] and we will review them. Notice that the scale µdependent part in the second term of (9) determines the RG equation of the three-point function.
Our main input is the new correction from the boundary terms of the effective action, i.e. the second term in The effective action must reproduce the boundary anomalies via the conformal transformation. As we will discuss in more detail, in addition to the Euler characteristic boundary term needed to preserve topological invariance, there exists new boundary Weyl invariants. To our knowledge, the corrections to the stress tensor correlation functions from the boundary anomalies have not been discussed in the literature.

Ward Identities and Anomalies Revisited
Considering a field theory coupled to a non-dynamical, curved background, we define the stress tensor in the Euclidean signature by where the functional W ≡ W (g(x)) is the effective action (the generating functional for connected Green's functions). The classical stress tensor trace vanishes for CFTs but at the quantum level, the regularization and renormalization introduce conformal symmetry breaking counter-terms that result in T µ µ (x) = 0; W denoted earlier is a part of W . We assume the theory is regulated in a diffeomorphism-invariant way and we will focus on the vacuum. The trace anomaly is a function of intrinsic curvatures (for a compact manifold) that must vanish in the flat limit. However, the flat space correlation functions depend on central charges: a result reminding us to adopt the correct order of limits, that the flat limit is imposed only after performing the metric variation.
We have the following definitions of the functional differentiation: where W | g=0 is the standard procedure getting rid of the source term in the action after performing the variation. Notice that we have restricted the results to flat space at the end of the computation. From these definitions, we can derive the following identity: where vanishes for a manifold without a boundary because the expectation value of the stress tensor vanishes in flat space. However, we will find a non-trivial contribution to T µν (x) | δµν when we include boundary Weyl anomalies. The conservation of the stress tensor implies An important relation that we can also derive from these definitions is the following trace conformal Ward identity relating the three-point function and the two-point function: where For CFTs in a compact manifold the last two lines of (17) do not contribute, and therefore these terms are not included in the literature. (Take d = 4 CFTs for example. After performing a metric variation on the stress tensor trace, the result vanishes in the flat limit since the trace anomaly is a function of curvature squared. The same argument applies for higher dimensions.) We see that in general the existence of the boundary Weyl anomaly modifies the trace Ward identity.

Compact Manifold
In this section we reproduce several main results for d = 4 CFTs without boundary terms. We follow [23] [24] closely but the way we adopt the dimensional regularization will be slightly different in the details. The conformal anomaly reads The type-A anomaly is defined by the Euler density, where δ αβγδ µνσρ is the fully anti-symmetrized product of four Kronecker delta functions. The Euler number χ is given by M d 4 x √ gE 4 = 4πχ. The type-B anomaly defines the central charge c; the Weyl tensor in d-dimensions (for d > 3) is given by which has basic properties: The last term in (19), referred to as the type-D anomaly, is scheme-dependent: adding an R 2 local counter-term shifts its anomaly coefficient. Here we use this freedom to remove this anomaly. See [26] and [27] for related discussions on the type-D anomaly. 1 Following expression (9), we write the stress tensor three-point function as where the µ-dependent part is determined by the counter-terms. We take In the dimensional regularization scheme, we consider a different treatment on the totally-antisymmetic tensor compared to [23] [24]. We write the Euler density as Defining σ(x) as the Weyl transformation parameter, g µν (x) → e 2σ g µν (x), we have µνλρ ) which confirms that the counter-terms produce the Weyl anomaly. Let us first consider a metric variation on (25). Performing integration by parts and recalling the metricity and the Bianchi identity, in varying the integrated Euler density we in fact only need to vary the vielbeins. From the relation 2δ/δg ν µ = e a (ν δ/δe a µ) , we obtain After performing two more metric variations and imposing the flat space limit, the result is However, the a-charge does not contribute to the RG equation of the three-point function because lim d→4 D E µνσραβ (x, y, z) = 0 due to the five totally anti-symmetized indices. It is the c-charge that controls the RG equation While the explicit expression of D W µνσραβ (x, y, z) is not essential for our purpose, we remark that, from the following re-writing one can also say the RG flow is dictated by the Q-curvature 2 : The above third-order metric variation can be performed straightforwardly. However, we do not find a compact expression so we do not list the result here. Finally, given the conformal anomaly (19), the Ward identity (17) can also be computed and the result is given by [23] T The anomalous terms are determined by the second-order expansion of the anomalies: where a projector is introduced that shares the same symmetries as the Weyl tensor: 2 In d = 4 the notion of Q-curvature was introduced by Branson and Ørsted [29] as a global conformal invariant.
However, note that the curvature Q4 in (30) is different from the Q-curvature defined in [29] by a total derivative, R. For an analysis of Q-curvature in d = 6, see [30].
(To the first order in the metric expansion, one has W µσρν ∼ 2P µσρν,αγδβ ∂ γ ∂ δ δg αβ .) From the general statement [23] that there are three independent coefficients in the three-point function, the relation (32) implies that there exists a linear relation between the three coefficients and the central charges a and c. (Assume that the type-D anomaly is removed.)

Boundary Terms of Anomalies and Effective Action
The complete classification of possible boundary terms based on the Wess-Zumino consistency for d = 4 CFTs was carried out recently by [17]. Denoting the induced boundary metric as h µν = g µν − n µ n ν where n µ is the unit-length, outward-pointing normal vector to ∂M, the full Weyl anomaly with boundary terms is given by where δ(x ⊥ ) is a Dirac delta function with support on the boundary. The Chern-Simons-like boundary term of the Euler Characteristic, is used to supplement the bulk density E 4 to preserve the topological invariance. (In the literature, it is also referred to as the boundary term for the Lovelock gravity. See [31] or [17] for a review.) We have denotedE αβ as the boundary Einstein tensor. There are two additional boundary Weyl invariants in the anomaly (36) and we refer to b 1 and b 2 as "boundary" central charges. The Weyl curvature in the last term of (36) and the Riemann curvature in (37) are pulled-back tensors. The traceless part of the extrinsic curvature is given bŷ with d = 4 here. (We still adopt the Greek indices for these boundary tensors, but note that their normal component is empty.) It is an important property thatK αβ transforms covariantly under the Weyl scaling. The boundary b 1 -type anomaly to our knowledge first appeared in [32]; [33] first pointed out the boundary b 2 -type anomaly. Spelling out the expressions, we have Notice that, because of these boundary terms, the Weyl anomaly (the log-divergent term of the partition function) of CFTs does not vanish even in flat space. In particular, the boundary b 1 -type anomaly constructed solely from the extrinsic curvature exists in any dimensions, and it might play a role to order the (boundary) RG flows in any-dimensional QFTs.
The conformal anomaly has its origin tracing back to the local counter-terms. Adopting the dimensional regularization, we should include the corresponding boundary terms in the effective action. We have the following identities: αβγδ .
The first identity is the consequence of the topological invariance when the boundary term of the Euler Characteristic is included; the last two identities reflect that they are covariant Weyl tensors. We find a compact expression for the boundary type-b 1 anomaly given by using three Kronecker delta functions. In the dimensional regularization, we write similar to the way we express the bulk Euler density in d dimensions using vielbeins in (25). We will work in Gaussian normal coordinates, x µ = {x ⊥ , x i }, such that x ⊥ = 0 is the local function for the boundary. The metric is given by Focusing on the response from varying the bulk metric of M, we keep the boundary metric h ij (x ⊥ = 0, x i ) fixed while we perform the variation. However, the normal derivative of the metric variation on the boundary can be non-zero in general. Under the coordinate (46), on the boundary we require The metric variation of extrinsic curvatures on the boundary then can be written by 3 We will not discuss the correlation functions of the boundary stress tensor obtained by varying the boundary metric, h. In this case the contribution instead comes from the tangential fluctuation of the metric along the boundary, and it will be independent of the bulk stress tensor correlation functions considered here.
It is instructive to verify the identity (42) using expression (45). We should show In the metric variation we receive contributions from varying the extrinsic curvatures and also from varying the vielbeins. 4 The latter procedure results in terms ∼K ω ρK σ ηK λ δ δ ρηδµ ωσλν and identity (42) can be recovered by contracting µ with ν. But one then needs to show that the contribution from varying the extrinsic curvatures with respect to metric is traceless. We can see this is indeed the case by writing the explicit metric expansion aŝ where a tensor is introduced, with the desired traceless property: H λαρ α = H βλ αλ = 0. The higher-order expansions of (52) vanish due to the boundary condition. Therefore, all we need to include on the boundary is the first-order part.

Renormalization Group Flow with a Boundary
Taking into account the boundary terms of the conformal anomaly (36), we would like to see how the RG equation (29) might be affected by boundary central charges. We shall focus on potential µ-dependent divergences coming from the boundary local counter-terms.
First we discuss the boundary term of the a-type anomaly. With the boundary term, the fact that the Lovelock gravity has a well-defined variational principle in the presence of a boundary implies that there will be no boundary contribution to the metric variation left over. In dimensions d = 4 + ǫ, the additional variation comes from varying the vielbeins, which only appear through the wedge products. Varying these vielbeins gives additional indices in the totally antisymmetric tensor that vanishes in the limit d → 4 5 . In short, the topological nature of Euler characteristic is still preserved in a non-compact manifold and the RG flow remains to be independent of the a-charge. The argument applies to any even dimensions.
The story is more subtle for the boundary Weyl invariants. Let us first consider the boundary b 2 -type. In a recent paper [34], it is shown that the metric variation of the following action: 4 We impose the boundary condition on the metric only in the physical dimensions. In the dimensional regularization scheme, the extra dimensions introduce additional vielbeins, whose variation on the boundary should be kept non-vanishing. It is the variation of these vielbeins that allows us to correctly reproduce the anomaly by varying an effective action. 5 For certain geometries one could obtain finite stress tensors. We will discuss them in the next section. However, for the RG equation only the µ-dependent poles are relevant.
does not have boundary terms containing normal derivatives of the metric variation left over. The tensor P αβµν has the same symmetries as the W αβµν and is defined by where P (0) αβµν = n α n ν K βµ − n β n ν K αµ − n α n µ K βν + n β n µ K αν , Note that P α µαν = 0. Under the conformal transformation, P αβµν → e 3σ P αβµν . To relate this formulation with the effective action of the type-B anomaly, we take n = 2 in (54). Using tr(P W ) = tr(P (0) W ) = 4K να W µναβ n µ n β we writẽ The effective action then might be defined by multiplying the above result by the µ ǫ ǫ factor. It is interesting to notice that the coefficient "8" of the boundary term suggests b 2 = 8c [34].
There are, however, some suspicions about this derivation that leads to the relation b 2 = 8c, which would mean that b 2 is not an independent central charge. First of all, from the classical viewpoint, in order to integrate the equations of motion of Weyl gravity, one needs to specify more boundary data. The theory of Weyl gravity might be still consistent (at the classical level) even if on the boundary there are some normal derivatives of the metric variation left over. (The Lovelock gravity is special in that the action is quadratic in time derivatives.) On the other hand, at the quantum level, it becomes not clear that the anomaly effective action must have a nice variational principle. In other words, one might allow some (traceless) delta function distributions on the boundary as a part of the contribution to the type-B anomaly induced stress tensor.
Because of these potential issues, one might adopt the following decomposition of the type-b anomalies related effective action: From the conformal invariance requirement, we are allowed to add more hKW on the boundary with the coefficient measuring the difference b ′ 2 = b 2 − 8c. Indeed, if the argument using the variational principle is invalid, one needs another derivation to verify b ′ 2 . However, we notice that the spin 0, 1 2 , 1 case's free field calculations given recently in [35] suggest the following universal result (independent of boundary conditions) This is the scenario we adopt in what follows. It is still interesting, though, to find a general proof of the result (61) without referring to the variational method.
Let us go back to the correlation function. As the consequence of the decomposition (60) and the result (61), the b 2 -type boundary term basically plays the role of a Gibbons-Hawking-like term for the bulk type-B anomaly effective action. Notice that (59) applies directly in d = 4 + ǫ dimensions. Therefore, this boundary b 2 -type does not generate a µ-dependent pole in the threepoint function, and hence the RG flow is not touched by the b 2 -charge.
Finally we consider the boundary b 1 -type. We will see this type of boundary Weyl invariant contributes to the RG equation in the vicinity of the boundary. This central charge depends on boundary conditions; see [35] for related heat kernel computation. Focusing on the µ-dependent, singular contribution in the d → 4 limit, we denote the correction as where Using expression (45), we obtain The normal component does not contribute and the correction only exists near the boundary. In summary, as the generalization of (29), the RG equation of the three-point function for d = 4 QFTs in a flat manifold with a boundary is given by µνσραβ (x, y, z) .
D W µνσραβ is determined by the Q-curvature via (31) and the boundary correction is given in (64).

Conformal Ward Identity with a Boundary and Stress Tensor
We next revisit the conformal Ward identity by including boundary anomalies. As the generalization of (32), we find where is the flat limit of the conformal anomaly. The results A σλρη (x, y, z) and B σλρη (x, y, z) are given by (33) and (34), respectively; the tensor S σραβ is defined in (18). We also have We have used the fact that when the Weyl tensor or the Riemann tensor is pulled-back on the boundary, there are no normal derivatives acting on the metric variation left over.
Notice that there is a µ-dependent pole from the two-point function proportional to the ccharge. As noticed by [23], this singular behaviour can be reproduced (or derived from definition (12)) by requiring compatibility between the three point function and the Ward identity (16). The result is given by where S µν = (∂ µ ∂ ν − δ µν ). (For the contact term, see (6).) Similarly, the boundary charge b 1 also contributes a pole to the two-point function. We have The correction does not have the normal component and only exists near the boundary. Let us also discuss the expectation value of the stress tensor in the flat limit. Normally the stress tensor vanishes in flat space, but here a curved boundary generating boundary Weyl anomalies can lead to non-vanishing stress tensors.
The boundary b 2 -type anomaly, as we have discussed earlier, does not touch the resulting metric variation because of the Gibbons-Hawking mechanism (we take b 2 = 8c). Therefore, there is no stress tensor correction by the b 2 -charge.
Next we consider the contribution from the boundary b 1 -type. We find where the second term comes from varying the vielbeins and is given by Notice that the metric h µ ν showing up inK µ ν contracts with the generalized delta function and it then generates the factor ǫ = d − 4, cancelling the pole. We have One might argue that the first term simply vanishes in (d − 1) → 3 (the physical dimensionality of the boundary) limit, due to the additional indices in the totally-antisymmetric tensor. However, if the boundary geometry has a special form of the extrinsic curvature, for example, where c is a constant, the contraction of indices can still lead to a finite result. In other words, the expression of the finite part of the stress tensor depends on the structure of the boundary 6 . For a ball whereK α β vanishes in d-dimensions, the stress tensor (73) vanishes directly. This type of boundary anomaly does not exist for the ball-like geometry.
The cylinder however provides a non-trivial example. The corresponding metric and extrinsic curvatures are given by It is useful to denote a delta functionδ µ ν = diag 0, 0, 1, 1, 1, 1, ... where r and z components are empty. The extrinsic curvature then can be written as K ν µ = 1 rδ ν µ . We consider a tube that has a radius r 0 so the boundary is set by δ(x ⊥ ) = δ(r − r 0 ). The first term of (75) gives which generates an overall d − 4 factor cancelling the pole. The second term of (75) evaluated in d = 4 contains the following contribution: where in the last equality we have used the cylindrical geometry. However, we see that this nonvanishing result does not generate a d − 4 factor, and therefore it leads to an infinite contribution to the stress tensor. To have better behaviour, we add the following regulator: with the coefficient c ′ = − b 1 16π 2 being adjusted to cancel the divergence when taking ǫ → 0. 7 This regulator does not touch the anomaly coefficient since it is manifestly Weyl invariant. We obtain the final (renormalized) stress tensor The stress tensor contributes near the boundary and there is no r-component contribution, as expected. Taking the trace on this stress tensor gives which reproduces the Weyl anomaly evaluated for a d = 4 cylinder where We have mentioned that the Euler characteristic boundary term provides the Gibbons-Hawking mechanism so that the RG flow is untouched by the boundary counter-term a ∂M E (bry) . However, there can be finite contribution to the stress tensor from varying the vielbeins. (Varying the vielbeins on the bulk Euler characteristic counter-term gives the a-type stress tensor in curved space [26] and the result vanishes in the flat limit.) Including the boundary term, the Euler characteristic gives the following contribution to the stress tensor in the flat limit Let us first consider a ball with a radius r 0 . Starting with bulk dimensionality d, we have ds 2 = dr 2 + r 2 dθ 2 + r 2 sin 2 θdφ 2 + r 2 sin 2 θ sin 2 φdψ 2 + ... , We obtain a finite result where we have sent d → 4 in the last equality. The trace of this stress tensor is which recovers the Weyl anomaly for a 4-ball where Let us next consider a cylinder. Notice that E (bry) 4 vanishes for a 4-cylinder and the anomaly comes solely from the boundary b 1 -charge. One might think that, because of the vanishing E , there should not have any a-type stress tensor correction near the boundary. However, it is important to recall that the geometry can be fixed only after working out the variation. We find for a 4-cylinder. The result is traceless, as expected.
It is of great interest to reproduce these new stress tensors in the vicinity of a boundary whose values are determined by boundary central charges using a different approach.

Boundary Weyl Anomalies
Folklore has it that there is no Weyl anomaly in an odd dimension. This statement is based on the fact that it is impossible to construct a scalar with the correct dimension using curvatures. However, for an odd-dimensional manifold with a boundary, there can be Weyl anomalies localizing on the boundary.
For a d = 3 manifold, the anomaly is given by We denoteR as the intrinsic Ricci scalar on the boundary. The trace-free extrinsic curvature in this case is given byK In this notation, we have (a 3 , c 3 ) = (−1, 1) for a conformal scalar with the Dirichlet boundary condition and (a 3 , c 3 ) = (1, 1) for the conformal Robin condition [36]. (See [34] for a recent discussion on d = 5 CFTs.) In this section we discuss how these boundary charges contribute to stress tensor correlation functions.

Correlation Function and Stress Tensor
We again focus on the correlation functions obtained by the bulk metric variation. We will not perturb the boundary metric, h, and the boundary conditions are (47), (48) and (49). The following identity relates the boundary counter-term to the a 3 -type anomaly: The bulk metric variation does not touch the boundary Riemann curvatureR a 1 a 2 µ 1 µ 2 which is intrinsic on ∂M. The contribution comes only from varying the vielbeins. The a 3 -charge does not generate a µ-dependent pole, and hence the RG flows of (bulk) stress tensor correlation functions are independent from this boundary central charge.
For the boundary Weyl invariant, trK 2 , we first re-write it as Working now in d = 3 + ǫ dimensions, we have the following identity: that relates this conformal anomaly to the boundary counter-term. Using the general expression (17), we obtain the Ward identity where the flat limit of the anomaly is the same as (107). We have defined S αβσρ in (18). We can also compute the RG equations of the correlation functions. The third-order metric expansion on the c 3 -type effective action (which contains only twoK µν here) does not generate a µ-dependent pole, and therefore there is no correction to the RG equation of the three-point function. For the two-point function, a relevant contribution comes from varying the extrinsic curvatures. (Varying the vielbeins will not give a µ-dependent pole.) The RG flow is given by This result applies for d = 3 QFTs in a flat manifold with a boundary. The correction appears only near the boundary.
The stress tensor trace, recovers the Weyl anomaly evaluated for a 3-ball where For a 3-cylinder with the radius r 0 , we have This result is traceless, as expected, since the a 3 -type anomaly vanishes for a 3-cylinder.

Conclusion
The general motivation of this paper comes from the constraints of the RG flows in both evenand odd-dimensional QFTs in flat manifolds with a boundary, and also comes from the universal contribution to the entanglement entropy. It is certainly of great interest to find more physical quantities characterized by boundary terms of conformal anomalies. Here we investigate the stress tensor correlation functions in flat spacetime with a generally curved boundary, focusing on the contribution from the boundary counter-terms. In particular, in d = 4, we find that the conformal Ward identity is modified by boundary central charges and the charge b 1 gives an additional correction to the RG equation near the boundary. Moreover, the boundary counterterms induce new stress tensors near the boundary. We have considered examples using a ball and a cylinder. We also discussed the similar story for d = 3 CFTs, where the Weyl anomaly exists only on the boundary. It will be interesting to generalize these results to five and six dimensions. Let us conclude by listing some questions that relate to boundary Weyl anomalies: (1) How do the boundary central charges modify the n-point functions of the stress tensor at non-zero separation?
(2) For four-dimensional QFTs with a boundary, (i) is it generally true that the edge central charge satisfies b 1(UV) > b 1(IR) ? (ii) Does the bulk a-charge still satisfy monotonicity under the RG flow in the presence of a boundary? (The boundary terms of the dilaton effective action calculated in [17] might be useful.) It will be also interesting to interpret these boundary Weyl anomalies in the AdS/CFT correspondence.

Metric Variation
For the convenience of the reader, here we list formulae for the metric variation of curvatures.

Two-dimensional CFTs with a Boundary
For a compact d = 2 manifold, the Weyl anomaly is given by The central charge notation c is more common in the literature. One has where χ is the Euler number. Using the expansion on the Ricci scalar around flat space, one obtains the Ward identity In the presence of a boundary one has The Ward identity is modified by the boundary term and we obtain T µ µ (x)T σλ (y) = − c 12π (∂ σ ∂ λ − δ σλ )δ 2 (x − y) − δ(x ⊥ )h σλ ∂ n δ 2 (x − y) +δ σλ δ 2 (x − y) T µ µ (x) | δµν + 2δ 2 (x − y) T σλ (x) | δµν .
The anomaly in the flat limit is given by Kδ(x ⊥ ) .
We define the vacuum stress tensor of a plane to vanish. For a disk with a radius r 0 , we have (Now we are working in d = 2 + ǫ dimensions.) Only the angular component is non-vanishing. The trace of this stress tensor reproduces the anomaly in the flat limit. Notice that in d = 2 curved space, the bulk stress tensor has to be expressed in terms of the Weyl factor σ(x) defined via g µν = e 2σ δ µν , because the Einstein tensor vanishes in d = 2. See [17] for more discussions on d = 2 CFTs.