Boundary Holographic Witten Diagrams

In this paper we discuss geodesic Witten diagrams in generic holographic conformal field theories with boundary or defect. Boundary CFTs allow two different decompositions of two-point functions into conformal blocks: boundary channel and ambient channel. Building on earlier work, we derive a holographic dual of the boundary channel decomposition in terms of bulk-to-bulk propagators on lower dimensional AdS slices. In the situation in which we can treat the boundary or defect as a perturbation around pure AdS spacetime, we obtain the leading corrections to the two-point function both in boundary and ambient channel in terms of geodesic Witten diagrams which exactly reproduce the decomposition into corresponding conformal blocks on the field theory side.


Introduction
Using the operator product expansion (OPE), correlation functions in conformal field theories naturally organize themselves into contributions from "conformal blocks". These blocks sum up the contributions of all the descendants associated with a given primary operator arising in the OPE of a given pair of operators. In field theories with holographic dual, correlation functions in the bulk can be calculated from Witten diagrams [1], that is position space Feynman diagrams in asymptotically anti-de Sitter (AdS) space. While the full correlation function must of course respect its decomposition into conformal blocks, the individual Witten diagrams do not nicely separate the contribution from any given block. It has been argued in [2] that one can isolate the contribution of a single conformal block to a conformal 4-pt function by calculating so called "geodesic Witten diagrams". Geodesic Witten diagrams differ from their standard cousins in that the bulk interaction vertices are only integrated along a geodesic connecting two boundary operator insertions instead of over all of AdS space. One can show that these geodesic Witten diagrams represent the contribution of a single block by explicit calculation. But a more elegant method is to demonstrate that they obey the defining Casimir differential equations that blocks must obey together with the correct boundary conditions.
In conformal field theories with boundaries (bCFTs) the notion of conformal blocks becomes more interesting [3]. The presence of a conformally invariant boundary reduces the conformal group in d spacetime dimensions from SO(d, 2) to SO(d − 1, 2). The reduced symmetry allows for the appearance of a non-trivial function depending on a single conformally invariant cross-ratio already at the level of the 2-pt function. This correlator can be decomposed into conformal blocks in two distinct ways: in the "ambient space channel" one uses the standard ambient space 1 OPE to re-express the two-point function as a sum over one-point functions (which need not vanish in a bCFT). The contribution of a given primary and its descendants gets summed up into an "ambient block". In the "boundary channel" one uses a novel operators expansion, the BOPE or boundary operator expansion, to expand any ambient operators in terms of boundary localized operators [3]. This way the full ambient space 2-pt functions gets reduced to a sum over 2-pt functions of boundary localized operators. Once again, the contribution of a given primary and its descendants can be summed up into a "boundary block". Demanding equality of the decompositions into ambient blocks and boundary blocks gives interesting constraints on the bCFT data, as encapsulated in the boundary bootstrap program [5].
It is natural to expect that the conformal block decompositions of the 2-pt function in a bCFT with holographic dual can once again be captured by geodesic Witten diagrams. First steps in this direction have been taken in [6]. bCFTs are dual to a spacetime with d + 1 non-compact directions that allow a slicing in terms of AdS d [7]. The simplest models have a d + 1 dimensional bulk given by a metric ds 2 = e 2A(r) ds 2 AdS d + dr 2 (1.1) potentially times some internal space. If e A = cosh(r/L) (hereafter we take the AdS radius L = 1) this metric is simply AdS d+1 . For the holographic dual of genuine bCFTs the standard holographic dictionary requires the warpfactor A to approach this asymptotic form for large r. Examples in this class include AdS sliced Randall-Sundrum models [7,8], the very closely related AdS/bCFT proposal by Takayanagi [9], as well as the d = 4 Janus solution of type IIB supergravity [10] together with its cousins in other d. The former two are toy models, based on Einstein gravity coupled to branes with tension. They have no known embedding in string theory and no explicitly known dual field theory. The latter is an explicit top-down solution; its dual field theory consists of N = 4 super Yang-Mills (SYM) theory with a step function defect across which the coupling constant jumps 2 . The non-trivial metric (1.1) is supported 1 As in [4] we use the term ambient space for the d-dimensional space-time (labelled by indices µ, ν) in which the d − 1 dimensional defect (labelled by i, j) is embedded. The direction transverse to the defect is called w. We reserve the term "bulk" for the d + 1 dimensional spacetime of the holographic dual given by (1.1). 2 Solutions like the Janus solution in which ambient space extends on both sides of the defect with different properties are often referred to as (holographic duals of) interface conformal field theories, or iCFTs. If the by matter fields whose profile is independent of the coordinates on the slice. We'll collectively denote these background fields as X(r). For example, in the Janus solution there is a single scalar turned on (the dilaton) with an r-dependent profile. In more general holographic bCFT constructions, such as the dual of N = 4 SYM on half-space with supersymmetry preserving boundary conditions [11][12][13], the warpfactor A depends non-trivially also on the compact internal space. In this case we can still use the metric (1.1) with the understanding that in this D dimensional metric r stands for the set of all internal variables, dr 2 is the D − d dimensional metric on the internal space, and A(r) really is, in general, a function of all these D − d internal coordinates.
In [6] only the simplest case of a holographic dCFT was addressed. In the case where the bCFT is really a dCFT with a small number of matter fields localized on a defect in a large N gauge theory, one can neglect the backreaction of the matter fields on the ambient space field theory. A simple top-down example of such a "probe brane dCFT" is the D3/D5 system of [8,14], representing N = 4 SYM coupled to a 2 + 1 dimensional hypermultiplet in the fundamental representation of the gauge group. In this case, the dual geometry remains AdS 5 throughout (that is e A = cosh(r) for all values of r). The defect is dual to a probe D5 brane living on one of the AdS 4 slices 3 . In this case the prescription for geodesic Witten diagrams is fairly straightforward as both geodesic and propagator retain their standard AdS form and all one has to account for are the extra brane localized matter and interactions. The resulting proposals of [6] can once again be confirmed by explicit calculation as well as by the Casimir method. But the prescriptions as phrased in this work heavily rely on the special probe brane scenario and it remained far from clear of how to implement the idea of geodesic Witten diagrams in generic boundary conformal field theories. It is our aim in this work to fill this gap.
In fact, the bulk manifestation of boundary conformal blocks has been understood in a slightly different context before. In [4] it was shown that the BOPE manifests itself as a mode decomposition in the bulk. This construction was used in [15] to show that particular integrals of bulk scalar fields along geodesics, the so called weighted X-ray transforms, are the correct bulk duals to the boundary conformal blocks. They naturally live in "boundary kinematic space". We will review both these constructions in the next section, as they will play a crucial role in deriving the correct geodesic Witten diagram prescriptions for generic holographic bCFTs. In section 3 we will give this derivation of the diagrams associated to the contribution of a single conformal block in the boundary and ambient channel respectively. After presenting some simple examples in section 4 we will conclude in section 5.
field theory on both sides of the interface is the same but extra degrees of freedom are localized on it, the system is often referred to as a defect CFT or dCFT. Both iCFTs and dCFTs can be seen as special cases of bCFT by employing the folding trick: the interface/defect can be viewed as a boundary in a theory whose ambient space contains two decoupled copies of the original CFT on half-space with interactions localized on the boundary. 3 For the simplest case of a D5 brane intersecting N D3 branes the probe is located at r = 0, but it can move to a different r * if we let some of the D3 branes end on the D5 brane [8].

Holographic Boundary Operator Expansion and Boundary Blocks
Underlying the decomposition of the 2-pt function into boundary conformal blocks is the notion of a boundary operator expansion (or BOPE). As was demonstrated in [4], the BOPE in the bulk can naturally be understood as a mode-decomposition of fields living on an AdS d sliced geometry with metric (1.1). We want to solve the equation of motion for a scalar field φ d+1 (r, y), dual to an ambient space operator O of dimension ∆. Here y and r stand collectively for the coordinates along the AdS d slice and the transverse directions respectively. Let us for concreteness focus on the case with only on internal variable r; the generalization to many r is straight forward as we will see as we go along. Using the d + 1 dimensional geometry (1.1) with background fields X(r) turned on the equation of motion reads g is the Laplacian in the full Janus background geometry (1.1) and ∂ 2 d stands for the AdS d Laplacian on the slice. The radial operator D 2 r is defined as The potential term M 2 (r) includes the bulk scalar mass M 2 0 , but also all the interactions with the background fields X(r) 4 . For example a quartic X 2 φ 2 coupling in the Lagrangian will give rise to an extra X 2 (r) term in M 2 . The only important property we need from M is that it does not depend on the y coordinates as guaranteed by the SO(d − 1, 2) defect conformal symmetry. We can make a separation of variables ansatz φ d+1 (r, y) = n ψ n (r)φ d,n (y) (2.3) so modes φ d,n obey a standard scalar wave equation on the slice The eigenvalues m 2 n are then determined by the internal mode equation: This 2nd order differential equation can easily be brought into the form of a 1d Schrödinger equation by a simple change of variables 5 ; correspondingly the modefunctions can be chosen 4 The dilaton is special because its coupling appears in front of a kinetic term. We will treat this case in section 4 as an example of our prescription. 5 A change of variables from r to a conformal coordinate z with dr = e A dz removes the e −2A factor in front of the eigenvalues m 2 n and a further rescaling ψn = e −(d−1)A/2 Ψn eliminates the first derivatives acting on the mode-function so that we are left with a standard Schrödinger equation for Ψn(z) together with its usual norm and an effective potential of [16] V (z) = 1 2 6) to be complete and orthonormal with respect to the Schrödinger norm, which in the original variables implies n ψ n (r)ψ n (r ) = e −(d−2)A(r) δ(r − r ) , dr e (d−2)A(r) ψ m ψ n = δ mn , (2.7) and the eigenvalues m 2 n are real. For cases with more than one internal coordinate or couplings of X(r) to terms involving derivatives of φ d+1 , it is straightforward to write down the corresponding eigenvalue problem. This is the only ingredient that needs to be changed in these cases.
Since each mode φ d,n on the slice obeys a scalar wave equation on AdS d , it is dual to an SO(d − 1, 2) primary operator o n localized on the defect. It was shown in [4] that these d − 1 dimensional operators are exactly the ones appearing in the BOPE of O, that is they appear in the expansion of O in terms of boundary localized operators We broke out the contribution c O 1 of the identity operator for clarity. In the BOPE as written in (2.8) the operators on the right hand side, labeled by k, are both primaries and descendants. However the contribution of the descendants are completely determined by that of the primaries, labeled by n, and they can be summed into non-local block operators B n ( x, w) so that the BOPE reads: The functional form of the blocks is uniquely fixed by symmetry [3]. According to [4] the primaries appearing in the BOPE are in 1-to-1 correspondence with the modes φ d,n , their dimensions are given by the eigenvalues m 2 n by the usual AdS d relation and the OPE coefficients c O n are encoded in the asymptotic fall-off of the modefunctions ψ n [4].
The decomposition of an ambient operator O into conformal blocks can be inserted into any correlation function which reduces the ambient space correlator into correlators of the nonlocal blocks. In order to implement this procedure for the special case of the 2-pt function 6 , we would like to get a bulk representation of the conformal block itself. This was provided 6 In the literature the name "conformal block" is often used both for the non-local operator B itself, but also for the non-trivial function of cross-ratios it contributes to a particular correlation function (usually the 4-pt function in CFTs without boundary and the 2-pt function in bCFTs). We will try to be careful in the following to reserve the name for the operator itself and will refer to the non-trivial function of the cross-ratio appearing in the 2-pt function as the contribution from a particular block.
in [15] using the construction of boundary kinematic space, which generalized the previously introduced kinematic space of [17][18][19][20] to bCFTs. It was shown that from the bulk field φ d+1 one can construct kinematic space operators R n φ by a weighted X-ray transform R n φ(y) = γ dr e (d−2)A ψ n (r)φ d+1 (r, y) . (2.11) Here γ is a symmetry enhanced geodesic, as defined in [15], emanating from the boundary point y at r → ∞. Of course there are many geodesics anchored at this point and their detailed properties depend on the warp-factor A. But it was shown in [15] that the line is a geodesic for any choice of warpfactor A(r) and is in fact singled out to be the only one compatible with the expected symmetries of the block. This is the geodesic γ that appears in the weighted X-ray transform. The weights ψ n are exactly the mode-functions appearing in (2.5). Using the orthogonality of the modefunctions as well as the mode decomposition (2.3) we see that the weighted X-ray transform exactly pulls out the AdS d modes φ d,n : Writing the AdS d metric on the slice parametrized by the d coordinates y as 7 it was further shown in [15] that R n φ( x, w) was in fact equal to the conformal block B n ( x, w), that is the radial direction along the slice takes the role of the direction orthogonal to the defect. In the following section we will use this insight to derive the contribution of a given boundary conformal block to the 2-pt function of two ambient space operators in terms of geodesic Witten diagrams.
3 Block decomposition of the 2-pt function In a conformal field theory without boundary, the conformal block decomposition is usually applied to 4-pt functions since they are the simplest correlators that allow any non-trivial functional dependence on the position of the insertion points. With 4 insertion points one can form 2 conformally invariant cross-ratios and the correlator can be a non-trivial function of both of them which can be decomposed (in different ways) into contributions from blocks appearing in the OPE applied to two operators at a time. For the case of a bCFT, a nontrivial cross-ratio already appears in the 2-pt function is conformally invariant and correspondingly the general form of the 2-pt function is Note that, unlike in the case without boundary, this 2pt-function need not vanish in the case when the two operators have different dimensions ∆ 1 and ∆ 2 . Both ambient and boundary block expansions can in principle be applied to the general case of two different operators. In the following we will however restrict ourselves to the case where both insertions are the same operator O with dimension ∆. Applying the BOPE to O( x 1 , w 1 ) and O( x 2 , w 2 ) separately, the 2-pt function can be written as a sum of contributions from 2-pt functions of the blocks. The residual boundary conformal group insures that the only non-vanishing 2-pt functions arise from one and the same primary o n appearing in both BOPEs. Furthermore, for 2-pt functions of scalar operators angular momentum conservation implies that only blocks build from scalar operators can contribute. This allows one to give a boundary-channel expansion for f (η) of the form c O 1 denotes the contribution of the identity operator. The contribution of the n-th boundary block, f ∂ (∆ n , η), is fixed by conformal invariance 8 . The explicit form of f ∂ can either be obtained by summing up the contributions of the descendants [3] or, more elegantly, by the Casimir method [5]: since all the descendants of a given primary sit in the same representation of the conformal algebra, they all have to correspond to the same eigenvalue of the Casimir operator L 2 ∂ of the conformal group. In terms of the generators D (dilatation), K i (special conformal), P i (translations), and M ij (boosts and rotations) in their standard representation as differential operators acting on a scalar at ( x, w) one has The differential equation together with the boundary condition that as the ambient operators approach the defect (η → ∞) the contribution of the blocks is dominated by the primary or in other words gives the contribution of the n-th block as [5] f A different expansion of f (η) can be obtained by using the standard ambient space OPE on the product O( . This reduces the calculation of the 2-pt function to a sum over 1-pt functions of all the operators appearing in the OPE of O with itself. In a standard CFT the only non-vanishing 1-pt function comes from the identity operator. So the 2-pt function is completely determined by the identity block. In a bCFT any scalar primary operator can have a non-trivial 1-pt function and so the 2-pt function of ambient operators can be reduced to a sum over 1-pt functions. This give us an expansion of f of the form where the sum over N is a sum over ambient space primaries (that is primaries under the full SO(d, 2)), λ N are the OPE coefficients and a N the constants determining their 1-pt functions according to (3.9). The contribution of the N -th ambient block can once more be obtained by explicit summation or the Casimir method [5]: Like f ∂ , f ambient can be found as an eigenfunction of a conformal Casimir operator, but this time it is the full SO(d, 2) Casimir (acting on x 1 , w 1 , x 2 , w 2 ) and the eigenvalue is Equating the two expansion gives rise to the boundary bootstrap equation [5], which will however not play a major role in this work.

Simple Examples: Dirichlet, Neumann and "No-brane"
The simplest examples to illustrate the various block decompositions are free field theories. The three cases one wants to distinguish are a free scalar bCFT with Dirichlet boundary conditions (the "Dirichlet theory"), a free scalar bCFT with Neumann boundary conditions (the "Neumann theory") or the "free No-braner" theory: a free scalar ϕ without boundary or interface in which one randomly picks the w = 0 surface to be treated as an interface. As emphasized in [4], in a free theory the BOPE is essentially a Taylor expansion. Correspondingly, the defect operators o n are built from O( x, w = 0) and ∂ w O( x, w = 0) and so on 9 . Due to the equations of motion, in a free theory ∂ 2 w O for O = ϕ is related to the on-slice Laplacian and hence is already a descendent of the defect conformal algebra. Correspondingly, the only operators that can appear in the BOPE are O (dimension ∆ ϕ = d/2 − 1) and ∂ w O (with dimension d/2). The Dirichlet theory only has the former, the Neumann theory the latter and the free no-braner has both. Correspondingly the boundary channel expansion yields [5] f which means that the full correlator takes the expected form one would get from a method of images construction Instead of interpreting the mirror charge terms as the contribution of boundary blocks with dimension ∆ ϕ + 1 and ∆ ϕ in the Dirichlet and Neumann case respectively, we can also give them an ambient channel representation: they are the ambient block of dimension 2∆ ϕ associated with the operator ϕ 2 appearing in the operator product of ϕ with itself. The difference in sign comes from the difference in vacuum expectation value of ϕ 2 . Note that the "no-braner" construction can be applied to any CFT, be it interacting or not. The fact that the BOPE truncates to only two primaries is special to the case of a free no-braner. In particularly, in the "no-braner" of a theory with holographic duals, one finds the entire tower of fields with dimension ∆ n = ∆ + n associates to the primaries built from

Holographic calculation
In principle it is easy to calculate a dCFT 2-pt function from Witten diagrams: we simply need to obtain the bulk-to-boundary propagator in the full holographic dCFT geometry of w w r w (x 1 ,w 1 ) (x 2 ,w 2 ) Figure 1. Full 2-pt function in an iCFT. The red triangle indicates the region in which the metric strongly deviates from AdS d+1 . In a bCFT space would truncate smoothly inside the triangle. We also indicated how the r and w coordinates used in the main text parametrize the space depicted in this figure. All following figures also depict the same r-w plane.
(1.1). This amounts to calculating a single Witten diagram as depicted in figure 1. So to some extend there is much less urgency in this case to organize the contributions according to conformal blocks. Nevertheless, it may sometimes be convenient to do so. We can use the exact representation in terms of the diagram of figure 1 in order to derive the geodesic Witten diagrams associated with the block decomposition.
In principle the Witten diagram in figure 1 is very easy to calculate. In terms of the bulk-to-boundary propagator K ∆,d+1 (r, (3.17) The constant C encodes the prefactor of the action 11 . The bulk-to-boundary propagator is, as usual, defined as a solution to the scalar equations of motion approaching the appropriate delta function at the boundary As in (2.5), M 2 (r) is a non-trivial function of the radial direction not just encoding the mass of the scalar field but also its interactions with all the matter fields that have a non-trivial 10 For the propagators and their relation to the correlation functions we are following the conventions of [21]. 11 For a scalar field, we take the action to be where M, N label d + 1 coordinates of the bulk spacetime.
profile in the background geometry. Asymptotically the metric approaches that of AdS d+1 and all matter fields go to constants so that M 2 (r) approaches M 2 0 which is related to the dimension ∆ by the usual The bulk-to-boundary propagator can also be obtained as a limit of the full bulk-to-bulk propagator G ∆,d+1 which obeys We can recover K ∆,d+1 via 12 (3.23)

Boundary Channel
Despite the simplicity of the Witten diagrams for the full 2-pt function leading to the expression (3.23) it can be helpful to decompose this full answer into a sum over blocks. Here we will use the full form of the 2-pt function to derive an expression in terms of boundary channel blocks. The contribution of a single block will be shown to be given by the diagrams in figure 2(b). In the special case that the holographic bulk dual is simply given by a probe brane our prescription reduces to the Witten diagram of 2(a) as it appeared previously in [6].
To derive the block decomposition, let us start with a representation of the bulk-to-bulk propagator in the full geometry in terms of mode functions. Instead of directly solving (3.21) for the propagator, we first find a set of appropriate eigenfunctions of the scalar wave operator. Usually one would look for modes of the form and then write the Green's function as (3.25) 12 As one approaches the boundary of asymptotically AdS d+1 space, the metric diverges as f −2 where the "defining function" f has a single zero at the boundary. To extract the asymptotic behavior of the various fields on this space, one multiplies with the appropriate powers of f . In the metric (1.1) one naively may have expected to use cosh(r) ∼ e r /2 as the defining function and hence multiply K ∆,d+1 simply with e ∆r /2 ∆ . However in this case one would obtain answers relevant for a field theory on AdS d . If we are interested in extracting correlators for a flat space bCFT, we need to use f = e r /(2w) as the defining function.  Here X stands for all the coordinates and K labels the modes. As long as the modes are a complete set, the right hand side automatically gives a delta function when acted upon with the wave operator. In our coordinate system it is easier to follow a slight variation of this strategy. We are looking for modes obeying In the Poincaré patch slicing of (2.14) k is a continuous label. As standard eigenfunctions of the wave equation, φ d,n,k form a complete set on the AdS d slice: where g 0 is the determinant of the AdS d metric on the slice. In terms of these we can write Using the defining equation (3.26) as well as completeness (2.7) and (3.28) of the modes it is straightforward to check that this representation indeed obeys the equation defining the bulk-to-bulk propagator: An important observation is to note that this representation can be used in order to reexpress the full d + 1 dimensional bulk-to-bulk propagator in terms of d dimensional bulkto-bulk propagators G 0 ∆n,d of the modes on the slice (that is, these are propagators for a d dimensional scalar on AdS d with mass squared given by m 2 n ). The 0 superscript here is reminding us that, unlike G ∆,d+1 , these d dimensional propagators on the slice are calculated using an AdS d geometry and so are completely known. Given a mode representation of G 0 and the separation of variables ansatz (2.3) we can easily see that the full d + 1 dimensional propagator can also be represented as 13 G ∆,d+1 (r 1 , x 1 , w 1 , r 2 , x 2 , w 2 ) = n ψ n (r 1 )ψ n (r 2 )G 0 ∆n,d ( x 1 , w 1 , x 2 , w 2 ) .

(3.32)
Our claim is that this is exactly the decomposition into boundary blocks. This immediately follows from the prescription of [15] for identifying the blocks. The bulk-to-bulk propagator 13 As a final check of our construction, note that it is easy to confirm that the propagator in this representation indeed obeys the defining equation (3.21): where we used the completeness of the radial modes (2.7) as well as is the correlation function of two bulk scalar insertions. We can perform a weighted X-ray transform (2.13) on both insertion points to extract the contribution of a given block: dr 2 e (d−2)(A(r 1 )+A(r 2 )) ψ n (r 1 )ψ n (r 2 )G ∆,d+1 (r 1 , x 1 , w 1 , r 2 , x 2 , w 2 ) .
(3.33) Furthermore the orthogonality of the mode-functions (2.7) yields That is, the conformal boundary block is supposed to be equal to the AdS d bulk-to-bulk propagator up to some normalization constant which we will determine shortly. The explicit form of the AdS d bulk-to-bulk propagator in the conventions of [21] is and ξ the chordal distance, which in the coordinates of (2.14) becomes or in other words Using a quadratic hypergeometric identity 14 we can also write this as Comparing with the expression for f ∂ (3.8) we see that they indeed have exactly the same functional form. This can also easily be verified by the Casimir method [6]; the form of the block is entirely fixed by conformal invariance. To confirm that the G n we obtained via the weighted X-ray transform indeed contribute to the CFT correlation functions exactly like a boundary block we need to plug in our representation of the propagator (3.32) into the formula for the full 2-pt function (3.23). To do so we need the asymptotic behavior of the mode functions ψ n (r). According to the analysis in [4] at large r one finds that is the fall-off is universally determined by the mass M 2 0 of the bulk scalar irrespective of n. The numerical factors C n are related to the BOPE coefficients as we will also make explicit below. Correspondingly we find for the 2-pt function Collecting all the numerical factors and comparing with the desired boundary channel decomposition (3.3) one finds (3.42) Recall that C was the prefactor of the action, C n the coefficient governing and fall-off of the mode functions, C ∆n,d the standard normalization factors (3.36) appearing in the bulk-tobulk propagator and c O n the BOPE coefficient. The factor of (2∆ − d) 2 is inherited from the relation between correlation function and bulk-to-bulk propagator (3.23).
The prescription of (3.41) with (3.42) provides a full non-perturbative decomposition of the 2-pt function into conformal blocks valid in any bCFT with holographic dual. It constitutes the main result of our work. If we want to describe the bCFT in terms of Witten diagrams in which the effect of the defect itself is treated order by order in a diagrammatic expansion, we need to look for a situation where the effects of the defect/interface can be taken into account perturbatively so that a calculation of the wavefunctions order by order in a diagrammatic expansion is meaningful to begin with. That is, we are interested in a situation in which the metric and background fields obey with ε a small parameter 15 . The mode functions encode the bulk geometry via the analog Schrödinger equation with potential (2.6) as emphasized in [15]. As long as the background has the form of (3.43), the potential fixing the Hamiltonian of the analog Schrödinger system takes the form where V 0 is the potential associated with an AdS d+1 geometry. Correspondingly we can expand the energy and the eigenfunctions as m 2 n = (m 0 n ) 2 + ε δm 2 n and ψ n = ψ 0 n + ε δψ n . (3.45) We will study the eigenfunctions of V 0 in all detail in the next section, for now it suffices to say that, using our result from (3.32) the propagator, is just a non-standard representation of the usual AdS d+1 bulk-to-bulk propagator. Now let us turn to the leading correction to the 2-pt function at small ε. Since ψ n = ψ 0 n + ε δψ n we also have C n = C 0 n + ε δC n . According to standard quantum mechanical perturbation theory we can write with (3.49) The leading correction to the 2-pt function hence becomes Eq. (3.50) can be read as our result for the geodesic Witten diagram prescription for calculating the contribution of a given boundary block to the scalar 2-pt function. Symmetry enhanced geodesics are used to identify points on the AdS d slices associated to a given boundary point. The first term corresponds to a diagram that represents the bulk-to-bulk propagator on the slice weighted by m =n γ mn C 0 m C 0 n . The second term is proportional to the correction of the propagator on the slice due to the shifted mass. Diagrammatically this could be represented by a mass insertion on the slice connected with two propagators, even though it is not clear this representation would be very illuminating. Note that the integrand of δm 2 n and γ mn is only non-vanishing in the region where the geometry differs from AdS d+1 as depicted in figure 2(b).
The case of probe brane considered in [6] and depicted in figure 2(a) looks similar but different in detail. In [6] the on-slice propagator is weighted by two bulk-to-boundary propagators as well as two interaction vertices. We will now show that this result in fact also naturally arises from our more general prescription. What is new in the probe brane case is that we have in the holographic theory matter fields which are completely localized on the brane. They are dual to defect localized operators in the CFT which do not arise from restriction of an ambient space operator to the defect but instead arise from matter that only lives on the defect. In terms of our bulk prescription such brane localized matter fields correspond to special modes ψ M of dimension ∆ M which, before accounting for the interactions, have no support near the boundary, meaning C 0 M = 0. As such, they make no contribution to the zeroth order 2-pt function. Since C 0 M vanishes, we also have δC n = γ nM C 0 M = 0. The leading contribution of these extra field comes from the correction to ψ M (and hence C M ): ψ 0 M has no support near r → ∞ and so does not contributed to correlation functions. The leading contribution of this mode to the 2-pt function hence is This result appears already much closer to the answer quoted in [6] in that it is order ε 2 , that is it involves two coupling insertions, and is proportional to the propagator of a brane localized field which did not contribute to the zeroth order answer. In [6] this propagator was decorated with two bulk-to-boundary propagators of the bulk field dual to the operator of dimension ∆. These bulk-to-boundary propagators connected the boundary operator insertion point with coordinates ( x, w) to the bulk point (r = 0, x, w) to which it is tied via a symmetry enhanced geodesic. For this special case the bulk-to-boundary propagators simply amount to inserting 16 the factors of (2w 1 ) ∆ and (2w 2 ) ∆ which arise automatically in our answer (3.53). So our answer is indeed completely equivalent to the one found in [6].

An explicit example: The Holographic No-braner
The simplest theory to test our decomposition (3.41) on is the holographic no-braner. That is, as described in subsection 3.1.2, we take our field theory to be a field theory without boundary with the plane at w = 0 treated as a defect. In this case the BOPE simply becomes a Taylor series. It was already found in [4] that this seemingly trivial example appears actually quite non-trivial from the point of view of the BOPE. Here we will see that also from the point of view of the conformal block expansion we require some seemingly miraculous cancellations. Let us start from the field theory point of view. Consider an operator O of dimension ∆ in a d dimensional CFT without any brane, boundary or defect. In this case, the full SO(d, 2) conformal group demands that (3.55) 16 The bulk to boundary propagator on AdS d+1 in the standard Poincaré coordinates with ds 2 = We can change variables to r, w, x by setting cosh(r)w −1 = W −1 , w tanh(r) = X d−1 , x = X . With this the boundary point (w, x) lives at (X d−1 = w, X = x) where as the bulk point (r = 0, w, x) at (W = w, X d−1 = 0, x). The bulk-to-boundary propagator connecting the two indeed just gives a factor of (2w) −∆ .
Here N is an overall normalization constant which could be chosen to be 1. Comparing this with the form (3.2) of the 2-pt function in a bCFT we see that this corresponds to This point is obvious from the ambient channel decomposition of the 2-pt function in (3.10).
In the "no-braner" theory only the identity operator has a non-trivial 1-pt function and so this is the only block that contributes. From the boundary channel point of view, we however have an infinite tower of boundary operators with dimension ∆ + n for n = 0, 1, 2, . . . contributing to the 2-pt function. The fact that O doesn't have a vev means there is no contribution from the identity operator in the boundary channel. The coefficients c O n in the decomposition (3.3) have to conspire in such a way that the contribution of all the blocks sums up to the simple expression in (3.56). In the special case of a free no-brane theory we saw that only two boundary operators contributed, but for a general bCFT (in particular one with a holographic dual) we genuinely need an infinite tower to sum up into the simple power law we are looking for.
This simple example can then be seen as one explicit solution to the conformal boundary bootstrap where an infinite number of boundary blocks reproduces the ambient channel decomposition in which only the identity contributes: with f ∂,n given by (3.8) with dimension ∆ n = ∆ + n. This requires a rather involved identity of hypergeometric functions and provides a non-trivial check for the holographic calculation of the coefficients in (3.42). Below we will prove the identity for the special case that d = ∆ = 4, but we expect (3.57) to hold for general d and ∆ 17 . The dynamics of a holographic bCFT is encoded in the bulk geometry via the warpfactor e A as well as the non-trivial background fields X(r). As emphasized above, the only place this data enters into the BOPE coefficients in (3.42) is via the C n , that is via the asymptotic fall-off of the mode functions. For the no-brane theory, the geometry is AdS d+1 , which means e A = cosh(r), and all other background fields are turned off. In this case the mode equation (2.5) can be solved analytically [7,24]. For d = ∆ = 4 one finds The overall prefactor can be determined by requiring the orthogonality condition (2.7). The 17 The identity (3.57) is a special case of eq. (A.7b) in [22]. If one inserts h = d/2, 1 = 2 = −∆ and ρ = η/(η + 4) into eq. (A.7b) in [22], one will obtain (3.57). Closely related identities have also been used in [23]. We would like to thank Christopher Herzog for pointing out [22,23] to us.

(3.64)
Now the sum over n is straightforward to do and so is the subsequent integral over t. We find, as hoped for, in perfect agreement with (3.56). 18 Ref. [4] chose to work with un-normalized wavefunctions and an n-independent fall-off Cn. In this case the non-trivial coefficient was obtained from the norm C −2 n ∼ dr e 2A ψ 2 n . 19 The general expression is . 20 We can simplify the product of Gamma functions by using the Legendre Duplication formula Γ(2z) = 2 2z−1 π 1/2 Γ(z)Γ(z + 1/2) on Γ(2n + 6).

Ambient Channel
The boundary channel expansion relied heavily on a new feature of bCFTs: the BOPE. In contrast, the ambient channel expansion uses the standard OPE for operators. The defect only makes its presence known via allowing non-vanishing 1-pt functions for scalar operators. Correspondingly the structure of the ambient channel expansion closely follows the pattern in the theory without defects [2]. First note that we can always think of the non-trivial metric (1.1) in the framework of being AdS d+1 plus deformations, as formally displayed in (3.43) with ε not necessarily small. The simple Witten diagram for the full 2-pt function of figure 1 can be thought of as an infinite sum of Witten diagrams in AdS d+1 with δg and δX insertions. A generic contribution is displayed in figure 3. For any expansion in terms of Witten diagrams to make sense, we need to focus on the case of small ε. In this case we consider diagrams in pure AdS d+1 and treat the deviations of the warpfactor from the e A = cosh(r) form as well as all matter fields with non-trivial profile X(r) as extra sources. We can derive the ambient channel expansion from the identity Here X and Z stand for bulk points (r 1 , x 1 , w 1 ) and (r 2 , x 2 , w 2 ) respectively. Y is a bulk point and integration region for Y is the whole of AdS d+1 spacetime. g denotes the Laplacian at Y . Expanding the identity around the background and taking the r 1,2 → ∞ limit, we obtain Here, as before, the 0 superscript indicates that these quantities take their un-deformed AdS d+1 values. The two bulk-to-boundary prop-  agators can be decomposed as [2] Here (x) n = Γ(x+n)/Γ(x) represents the Pochhammer symbol. The geodesicγ 0 parametrized by λ is the usual AdS d+1 geodesic connecting the boundary points ( x 1 , w 1 ) and ( x 2 , w 2 ). Plugging this decomposition into (3.67), we obtain conformal block expansions of the ambient channel,

(3.69)
This shows that there is basically no difference between a probe brane setup (in which case the sources are localized at a single locus r = r * as depicted in figure 4(a)) and the generic holographic bCFT geometry of (1.1) (in which case the sources have support over some region in r as depicted in figure 4(b)). Note that this integral only involves the uncorrected AdS d+1 propagators. The proof in [6] that this is indeed an eigenfunction of the conformal Casimir with the right boundary conditions and hence corresponds to a contribution of a single ambient channel block only relied on the properties of the geodesicγ 0 and the bulk-to-bulk propagator G 0 ∆,d+1 . It applies immediately to our case as well.
The one interesting upshot of this analysis is that to leading order the only aspect of the bulk geometry that affects the 2-pt function of two scalar operator O are the bulk scalars X(r) with non-trivial profile. This is due to the fact that only scalar blocks contribute in the 2-pt function of ambient space scalars. This follows immediately from angular momentum conservation.

Equivalence of two different decompositions
In subsections 3.2 and 3.3 we decomposed 2-pt functions in two different ways; the boundary channel (3.50) and the ambient channel (3.69). These decompositions should be the same. In this subsection we give an explicit proof of this equivalence. At leading order we confirmed this in subsection 3.2.1 when discussing the no-brane case.
Our mode decomposition of the AdS d+1 propagator (3.46) implies a similar representation for the bulk-to-boundary propagator via the limiting procedure of (3.22). Using this representation for the bulk-to-boundary propagators in our ambient channel result (3.67), we get The first term is perfect agreement with that of the leading correction of boundary channel (3.50). To show the equivalence of the second term, some further computations are required. Using (3.31) and the completeness relation of φ d,n,k , we obtain 21 This computation shows that the second terms are also same and we proved that the boundary channel (3.50) and the ambient channel (3.69) are the same exactly. In this proof we do not assume any specific form for δV (r). We can apply this proof to any case.

Example: The Janus iCFT
In this section, we consider the Janus iCFT as an example. The coupling constant of the Janus iCFT jumps across the interface. From the bulk point of view, this is because the dilaton field is not constant and has a non-trivial profile. Regarding the dilaton as a source, we consider the 2-pt function of operators dual to axions. The bulk dual of the Janus iCFT is a solution of type IIB supergravity. The dilaton field φ of Janus has a non-trivial profile and depends only on the radial coordinate. The dilaton satisfies the equation of motion, For a dilaton with only dependence on the radial coordinate r this implies where ε an integration constant and will be assumed to be small when we consider perturbation theory. ε is proportional to the jump in coupling constant in the dual iCFT. The corresponding correction to the metric is of order ε 2 and so, as in (3.43), can be neglected to obtain the leading order correction to the 2-pt function. 21 In terms of diagrams what we are saying is that the change of the AdS d propagator of a field with shifted mass (m 0 n ) 2 + ε δm 2 n can be obtained from a Witten diagram with an interaction vertex δm 2 n , integrated over all of AdS d , connecting two propagators associated to mass (m 0 n ) 2 .
The action of type IIB supergravity in the Einstein frame contains a coupling term between the axion field a and the dilaton field, Since the dilaton field is not constant, the 2-pt function of the dimension 4 operator TrF ∧ F dual to the axion field will be modified by the dilaton field. Plugging a mode expansion of the axion field as in (2.3) a = n ψ n (r)φ d,n (y) into the equation of motion of the axion it reduces to e −2φ D 2 r + e −2A−2φ m 2 n + ∂ r e −2φ · ∂ r a = 0 . Naively, if we expand the dilaton term in the above equation in ε and regard the leading term of order ε as δV , we may obtain leading correction terms to modefunctions and energy. But the problem is not so simple because this naive potential δV contains a first derivative term about r and the dilaton term appears in front of the energy term. Thus, a more careful treatment is required to discuss perturbation theory in this case.
As noted in the previous footnote 5, to use standard quantum mechanical perturbation results we should first change the variable from r to z with dr = e A dz. Furthermore, we rescale the field as ψ n = e −(d−1)A/2+φ Ψ n to remove the first derivative term. Then the problem reduces to standard quantum mechanics with an energy E n = m 2 n /2, a kinetic term −(1/2)d 2 /dz 2 , a potential (where we introduce a mass term though this term does not appear in the equation of motion of the axion) and a standard normalization 22 Let's return to the discussion of perturbation due to the dilaton profile. Expanding the potential to leading order in ε, we obtain When the dilaton term is in front of the kinetic term, the original orthogonality relation is dr e (d−2)A−2φ ψmψn = δmn .

Conclusion and Discussion
In this paper we discussed the conformal block expansion of 2-pt functions in general holographic bCFTs. In section 3.2, we provided the decomposition of the 2-pt function in the boundary channel. This was accomplished by decomposing the bulk-to-bulk propagator on the full d + 1 dimensional geometry into the radial direction and AdS d space. It was shown that conformal blocks in the boundary channel are given by bulk-to-bulk propagators on the AdS d slice. We also obtained the leading correction of 2-pt functions by the perturbation around pure AdS and background fields.
We also confirmed that our conformal block expansion works in the case without boundary. The summation of conformal blocks can be written as a 1-pt function of just the identity operator. This is an expected result because all 1-pt functions except that of the identity operator vanish without boundary, but reproducing this in the boundary channel expansion proved to be surprisingly tedious.
In section 3.3, we discussed the ambient channel. We provided the leading correction due to conformal blocks in the ambient channel from first principles. The contribution of a given conformal block contains products of two bulk-to-boundary propagators and one bulk-tobulk propagator. They intersect at points on the geodesic as in the probe brane case [6]. The remaining point connected to its bulk-to-bulk propagator couples to a source term. When the source term is a delta function in the radial direction, our decomposition reduces to that of [6]. We also proved the equivalence between the two decompositions, boundary channel and ambient channel, in section 3.4.
In section 4, we considered d dimensional Janus solutions as an example. Since the d = 4 Janus solutions is constructed from type IIB supergravity, the dual CFT is known explicitly. So, the Janus geometry is a good example. In Janus, the source is an odd function with respect to the radial coordinate r. Hence, the conformal dimension is not affected by the source and only γ mn is non-trivial. We explicitly computed the potential and obtained γ mn as integrals over Legendre polynomials . Our prescription can be easily generalized to other cases.