# The functional bootstrap for boundary CFT

Apratim Kaviraj$^{a,b}$ and Miguel F. Paulos$^a$

$^a$Laboratoire de Physique Théorique de l'École Normale Supérieure, PSL University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, 24 rue Lhomond, 75231 Paris Cedex 05, France

$^b$Institut de Physique Théorique Philippe Meyer, de l'École Normale Supérieure, PSL University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, 24 rue Lhomond, 75231 Paris Cedex 05, France

E-mail: apratim.kaviraj@lpt.ens.fr, miguel.paulos@lpt.ens.fr

**ABSTRACT:** We introduce a new approach to the study of the crossing equation for CFTs in the presence of a boundary. We argue that there is a basis for this equation related to the generalized free field solution. The dual basis is a set of linear functionals which act on the crossing equation to give a set of sum rules on the boundary CFT data: the functional bootstrap equations. We show these equations are essentially equivalent to a Polyakov-type approach to the bootstrap of BCFTs, and show how to fix the so-called contact term ambiguity in that context. Finally, the functional bootstrap equations diagonalize perturbation theory around generalized free fields, which we use to recover the Wilson-Fisher BCFT data in the $\epsilon$-expansion to order $\epsilon^2$.

**KEYWORDS:** Conformal Field Theory, Nonperturbative Effects, AdS-CFT Correspondence

ARXIV ePRINT: 1812.04034
# Contents

1 Introduction 1

2 Review of the boundary bootstrap 2

3 Functionals and a simple basis 5

3.1 Functional ansatz and constraints 5

3.2 Warm-up: a simple basis 8

3.2.1 Special values of $\Delta_{\phi}$ and orthonormal functionals 8

3.2.2 Shifted functionals and general $\Delta_{\phi}$ 10

3.3 Functional bootstrap equations 11

4 Bootstrap equations and BCFT basis 14

4.1 Witten diagrams and functional equations 14

4.2 Constraining the kernels 17

4.3 Boundary functionals 18

4.3.1 General solution 18

4.3.2 Neumann contact term ambiguity 20

4.3.3 Shifted functionals 21

4.4 Bulk functionals 22

4.5 Functional actions 23

5 Witten diagrams for BCFT 25

5.1 Setup 25

5.2 Boundary exchange 26

5.3 Bulk exchange 28

5.4 Contact diagrams 30

6 Polyakov and functional bootstrap 31

6.1 Polyakov blocks 31

6.2 Dirichlet case 33

6.3 Neumann case 33

7 Application — Wilson-Fisher at $O(\epsilon^2)$ 35

8 Discussion and outlook 39

A Functional actions 41

A.1 $\alpha$ and $\beta$ functionals 41

A.1.1 Action on boundary blocks 41

A.1.2 Acting on bulk blocks 42

A.2 Bulk functionals 43
B Computation of bulk exchange Witten diagram 44

C Witten diagram block decompositions 46

C.1 Boundary exchange diagram 46

C.2 Bulk exchange diagram 48

C.3 Contact diagram 50

# 1 Introduction

The bootstrap philosophy [1, 2] was first applied to boundary conformal field theories (BCFTs) in [3] building on general results of [4, 5] and the methods introduced in [6]. In that work, unitarity and crossing symmetry was used to derive numerical bounds for general BCFTs, as well as exact results for the Wilson-Fisher fixed points in the presence of a boundary. Since then, many other works have studied and generalized this to conformal defects in a number of situations [7–14].

The usual bootstrap starting point are the two point functions, which in the presence of a defect are non-trivial functions of one or more cross-ratios. A notable feature is that unlike in the usual CFT bootstrap, the crossing equation for such correlators involves coefficients which lack a definite sign, even assuming unitarity. This makes ordinary numerical methods harder or impossible to justify. Approaches based on sparseness of extremal solutions [15, 16] have so far obtained limited, but interesting, results [12, 15], which should be extended and explored.

Very recently, the works [17–19] have introduced a new analytic approach to studying 1d CFT correlators, or more generally CFTs restricted to a line.<fnref n="1" /> In this approach one introduces a complete set of linear functionals associated to the generalized free field solutions to crossing. They are the extremal functionals [16, 21] that correspond to these solutions, and are associated to optimal bounds on the CFT data. By acting on the crossing equation with these functionals one obtains a set of sum rules that this data must satisfy. These sum rules are a complete reformulation of the original crossing equation, and have many nice properties. Among them they provide upper and lower bounds on the Operator Product Expansion (OPE) data, diagonalize perturbation theory around generalized free solutions, and also allow a rigorous derivation of the (Mellin-)Polyakov bootstrap [2, 22–25]. In particular, in this 1d context it fixes the polynomial (or contact-term) ambiguity, that limits the applicability of the Polyakov bootstrap in higher dimensions [26, 27].

The success of this approach suggests we should attempt to generalize it to other settings. In this work we will initiate the study of the functional approach to the boundary bootstrap. This is particularly nice as our results will be immediately relevant for higher-dimensional CFTs in the presence of a boundary, providing a useful stepping stone towards a functional bootstrap analysis of the higher-d CFT crossing equation.

<footnote n="1">¹See also [20] for related work.</footnote>
Concretely, in this paper we will explain how to construct suitable bases of linear functionals that act on the crossing equation for a scalar CFT correlator $\langle\phi\phi\rangle$ on a half-space. These functionals are associated to the generalized free boson CFT in the presence of a boundary with Neumann or Dirichlet boundary conditions. Technically the construction shares sufficient features of the 1d CFT case that will allow us to make progress in a relatively straightforward manner.

We act with our set of functionals on the crossing equation to obtain a set of functional bootstrap equations contraining the OPE and Boundary Operator Expansion (BOE) data of the BCFT. As we explain, these equations can be equivalently thought of as arising from a Polyakov-type approach to bootstrapping the BCFT crossing equation. In particular, the functional actions on conformal blocks compute the conformal block decompositions of Witten diagrams in an $AdS_{d+1}$ geometry in the presence of a radially extended $d$-dimensional brane. The functional bootstrap equations then guarantee that two-point functions can be written in a manifestly crossing symmetric manner in terms of Polyakov blocks for BCFTs.

By perturbing the functional bootstrap equations around the generalized free solution we can continuously deform it to more interesting interacting solutions. In particular, if we assume the states in the spectrum to be essentially deformations of those of the generalized free case, apart from eventually some fixed number of states which we control, the deformation is uniquely fixed. This is the idea of extremality: sparseness of a solution to crossing can be used to uniquely determine it. Here we will use this approach to rederive Wilson-Fisher BCFT data up to $O(\epsilon^2)$ [28].

The outline of this paper is as follows. In the next section we briefly review the key points of the boundary bootstrap which are relevant to this work. Section 3 introduces the general ansatz for the functionals we will consider, and constructs a simpler, pedagogical basis of such functionals which can be used to derive a toy version of the functional bootstrap equations and the Polyakov bootstrap. The construction of the main functional basis takes place in section 4, which is motivated by considering Witten diagrams in $AdS$ space in the presence of a brane. This picture is confirmed in 5 by explicitly computing the relevant Witten diagrams, where we show the functional actions appear as coefficients of bulk double trace operators and boundary derivative operators. In section 6 we explain in detail how the functionals we construct are related to a Polyakov-like approach to the bootstrap of BCFT correlators. As an application of the functional bootstrap equations, we study them in perturbation theory around free theory in $4-\epsilon$ dimensions to recover the Wilson-Fisher BCFT OPE and BOE data to $O(\epsilon^2)$. We conclude with a discussion and outlook. The paper is complemented by technical appendices.

**Note added.** While this work was being completed we became aware of the work [29] which overlaps with ours.

## 2 Review of the boundary bootstrap

In this section we will briefly review the bootstrap approach to BCFTs. We will be brief, referring the reader to the literature for further details, e.g. [3–5].
We consider d-dimensional CFTs on the Euclidean half-space with coordinates:

$$ (\vec{x}, x^{\perp}) = (x^1, \dots, x^{d-1}, x^{\perp}), \quad \text{with } x^{\perp} \ge 0. \qquad (2.1) $$

A boundary CFT arises when it is possible to choose boundary conditions at $x^{\perp} = 0$ in such a way that a SO(1, d) subgroup of the conformal group is unbroken, namely the group of conformal transformations which preserve the hyperplane $x^{\perp} = 0$, which we will call the “defect”. Under these circumstances, the theory contains both the original bulk operators $\mathcal{O}$ as well as operators $\hat{\mathcal{O}}$ which live on the defect. The latter can be thought of as arising by taking bulk operators sufficiently close to the defect. Correlators containing only defect operators satisfy all the axioms of an ordinary CFT with the exception of the existence of a stress-tensor.<fnref n="2" />

Since the original conformal group is partially broken, bulk scalar operators can now develop expectation values while preserving the SO(1, d) subgroup:

$$ \langle \mathcal{O}(\vec{x}, x^{\perp}) \rangle = \frac{a_{\mathcal{O}}}{(x^{\perp})^{\Delta_{\mathcal{O}}}}. \qquad (2.2) $$

In turn this implies that the bulk two point function can become non-trivial. Indeed, with two points we can make up a cross-ratio which is invariant under the unbroken conformal transformations:<fnref n="3" />

$$ z = \frac{4x_1^{\perp}x_2^{\perp}}{|\vec{x}_{12}|^2 + (x_1^{\perp} + x_2^{\perp})^2}, \quad \vec{x}_{ij} \equiv \vec{x}_i - \vec{x}_j. \qquad (2.3) $$

In terms of this cross-ratio we can write a scalar two-point function as

$$ \langle \phi_1(\vec{x}_1, x_1^{\perp}) \phi_2(\vec{x}_2, x_2^{\perp}) \rangle = z^{\frac{\Delta_{\phi_1} + \Delta_{\phi_2}}{2}} \frac{\mathcal{G}(z)}{(x_1^{\perp})^{\Delta_{\phi_1}} (x_2^{\perp})^{\Delta_{\phi_2}}}. \qquad (2.4) $$

We have pulled out a power of $z$ for convenience. In the bulk we have an OPE which takes the schematic form:

$$ \phi_1 \times \phi_2 = \sum_{\mathcal{O} \in \phi_1 \times \phi_2} \lambda_{\phi_1 \phi_2 \mathcal{O}} (\mathcal{O} + \text{descendants}) \qquad (2.5) $$

We can use this to derive an expansion for the correlator around $z=1$:

$$ \mathcal{G}(z) = \sum_{\mathcal{O} \in \phi_1 \times \phi_2} a_{\mathcal{O}} \lambda_{\phi_1 \phi_2 \mathcal{O}} G_{\Delta_{\mathcal{O}}}(z | \Delta_{\phi_1}, \Delta_{\phi_2}), \qquad (2.6) $$

where we have introduced the bulk conformal block,

$$ G_{\Delta}(z | \Delta_{\phi_1}, \Delta_{\phi_2}) = \frac{(1-z)^{\frac{\Delta - \Delta_{\phi_1} - \Delta_{\phi_2}}{2}}}{z^{\frac{\Delta_{12}}{2}}} {}_2F_1 \left( \frac{\Delta}{2} + \frac{\Delta_{12}}{2}, \frac{\Delta}{2} - \varepsilon + \frac{\Delta_{12}}{2}, \Delta - \varepsilon, 1-z \right), \qquad (2.7) $$

with

$$ \varepsilon \equiv \frac{d-2}{2}, \quad \Delta_{12} \equiv \Delta_{\phi_1} - \Delta_{\phi_2} \qquad (2.8) $$

<footnote n="2">²Since the half-space is conformally invariant to AdS, this can be thought of as a special case of the construction in [30], where the bulk QFT happens to be a CFT.</footnote>

<footnote n="3">³Note that this is related to the more usual $\xi = \frac{|\vec{x}_{12}|^2 + (x_1^{\perp} - x_2^{\perp})^2}{4x_1^{\perp} x_2^{\perp}}$ by $\xi = \frac{1-z}{z}$.</footnote>
which captures contributions to the two-point function from a bulk primary and all its descendants.

It is possible to do a different decomposition of the same two point function [31]. This arises from the fact that any bulk operator can be expanded in a basis of defect operators. This is a direct consequence of radially quantizing the theory around a point contained in the defect. This boundary operator expansion (BOE) takes the form

$$ \phi(\vec{x}, x^{\perp}) = \sum_{\hat{\phi}} \mu_{\hat{\phi}}(x^{\perp})^{\Delta_{\hat{\phi}} - \Delta_{\phi}} C[(x^{\perp})^2 \partial^2] \hat{O}(\vec{x}) \quad (2.9) $$

where $C[(x^{\perp})^2 \partial^2]$ captures contributions of boundary descendants [4]. Using the BOE we now have

$$ \mathcal{G}(z) = \sum_{\hat{\phi}} \mu_{\hat{\phi}}^2 \hat{G}_{\Delta_{\hat{\phi}}}(z|\Delta_{\phi_1}, \Delta_{\phi_2}) \quad (2.10) $$

where we've defined the boundary blocks:

$$ \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi_1}, \Delta_{\phi_2}) = z^{\hat{\Delta} - \frac{\Delta_{\phi_1} + \Delta_{\phi_2}}{2}} {}_2F_1(\hat{\Delta}, \hat{\Delta} - \varepsilon, 2\hat{\Delta} - 2\varepsilon, z). \quad (2.11) $$

From now on we will be interested in the case where $\Delta_{\phi_1} = \Delta_{\phi_2} = \Delta_{\phi}$. Equating the two different expansions for the correlator we obtain the boundary bootstrap crossing equation:

$$ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi}) = \sum_{\Delta} (a\lambda)_{\Delta} G_{\Delta}(z|\Delta_{\phi}). \quad (2.12) $$

From now on we'll want to think of the above as a mathematical equation for the OPE and BOE data, whose solutions may or not correspond to a physical theory. Accordingly the sums now run over pure labels $\hat{\Delta}$, $\Delta$ and we have made several self-evident abreviations.

The BCFT crossing equation above is very similar to the $D=1$ CFT bootstrap crossing equation equating s- and t-channel expansions of a four-point function of identical operators. Indeed, for $\varepsilon = 0$ and setting $(a\lambda)_{\Delta} = \mu_{\Delta/2}^2$ the two become equivalent. More generally however, not only are the two channels now independent, but crucially the bulk channel does not have positivity, since the $(a\lambda)_{\Delta}$ can have either sign.

Although the original equation was defined only for $z \in (0, 1)$ in the Euclidean region, we actually have convergence in the much wider region $\mathcal{R} \equiv \mathbb{C} \setminus (-\infty, 0) \cup (1, \infty)$. This region is obtained by complexifying $z$ and avoiding the branch cuts of the bulk and boundary blocks.<fnref n="4" />

## The generalized free solution

An important example of a CFT in the presence of a boundary is the generalized free field (GFF) with Dirichlet or Neumann boundary conditions. This is a theory with an elementary field $\phi$ whose correlators are given by Wick contractions. The two point function

<footnote n="4">⁴It is possible to prove this in the same way as the analogous statement for the $D=1$ CFT, namely by introducing a radial coordinate $\rho(z)$ which maps the cut complex plane to the disk [32].</footnote>
in turn is given as

$$ \begin{aligned} \langle \phi(\vec{x}_1, x_1^{\perp}) \phi(\vec{x}_2, x_2^{\perp}) \rangle &= \frac{1}{(|\vec{x}_{12}|^2 + (x_1^{\perp} - x_2^{\perp})^2)^{\Delta_{\phi}}} + \nu \frac{1}{(|\vec{x}_{12}|^2 + (x_1^{\perp} + x_2^{\perp})^2)^{\Delta_{\phi}}} \\ \Leftrightarrow \mathcal{G}(z) &= \frac{1}{(1-z)^{\Delta_{\phi}}} + \nu, \end{aligned} \quad (2.13) $$

where $\nu = 1, -1$ for Neumann, Dirichlet boundary conditions respectively. In the bulk channel the correlator contains the identity and double trace operators $\phi \Box^n \phi$ with dimension $\Delta_n = 2\Delta_\phi + 2n$. In the boundary channel one finds operators with dimension $\hat{\Delta}_n^+ = \Delta_\phi + 2n$ for Neumann boundary conditions, and $\hat{\Delta}_n^- = \Delta_\phi + 1 + 2n$ for Dirichlet, corresponding to (transverse) derivative operators $\partial_\perp^{2n} \phi(\vec{x}, 0)$ or $\partial_\perp^{1+2n} \phi(\vec{x}, 0)$ respectively. The OPE and BOE data in this case are [3]:

$$ (a\lambda)_{\Delta_n} = \nu(-1)^n \frac{(\Delta_\phi)_n (2\Delta_\phi - 1 - \varepsilon + 2n)_{-n}}{(\Delta_\phi - \varepsilon + n + 1)_{-n} n!} \quad (2.14a) $$

$$ \mu_{\hat{\Delta}_n^+}^2 = \frac{(\Delta_\phi)_{2n} (\Delta_\phi - \varepsilon - \frac{1}{2} + 2n)_{-n}}{2^{2n-1} (2n)! (\Delta_\phi - \varepsilon + n)_{-n}}, \quad \text{Neumann} \quad (2.14b) $$

$$ \mu_{\hat{\Delta}_n^-}^2 = \frac{(\Delta_\phi)_{2n+1} (\Delta_\phi - \varepsilon + \frac{1}{2} + 2n)_{-n}}{2^{2n} (2n+1)! (\Delta_\phi - \varepsilon + n)_{-n}}, \quad \text{Dirichlet} \quad (2.14c) $$

The GFF theory can be thought of as a free scalar field with mass $m^2 = \Delta_\phi(\Delta_\phi - d)$ living in $AdS_{d+1}$ [33]. In the present context, we have placed a brane in AdS perpendicular to its boundary, partially breaking the AdS isometries. The induced metric on this brane is $AdS_d$, which makes manifest that a subgroup of conformal transformations is preserved. On this brane the field is chosen to have Neumann or Dirichlet boundary conditions.<fnref n="5" /> The resulting theory trivially reproduces the above correlators once AdS fields are pushed out to the boundary. It is natural however to think of deformations of the solution above, introduced by adding interactions in $AdS_{d+1}$. This will play an important role in our story later on.

# 3 Functionals and a simple basis

## 3.1 Functional ansatz and constraints

We propose that the optimal way to extract information from the crossing equation is to act on it with a suitable basis of linear functionals. Note that we will use the word ‘basis’ somewhat loosely in this work. In particular a basis of functionals should be understood as a set of functionals which establishes both necessary and sufficient conditions for crossing symmetry to hold. However in this paper we focus on constructing a set of functionals without proving their sufficiency.

The functionals should have the property that they commute with the infinite sums over states in the crossing equation — the swapping condition [34] — and lead to finite

<footnote n="5">⁵This is not to be confused with the condition at the radial boundary of $AdS_{d+1}$.</footnote>
results. Let us denote one such functional by $\omega$. Acting on the crossing equation we obtain a sum rule

$$ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \omega(\hat{\Delta}|\Delta_{\phi}) = \sum_{\Delta} (a\lambda)_{\Delta} \omega(\Delta|\Delta_{\phi}), \qquad (3.1) $$

which should be thought of as a constraint on the BCFT data. In the above we have defined the useful abbreviations:

$$ \omega[\hat{G}_{\hat{\Delta}}(z|\Delta_{\phi})] \equiv \omega(\hat{\Delta}|\Delta_{\phi}), \qquad \omega[G_{\Delta}(z|\Delta_{\phi})] \equiv \omega(\Delta|\Delta_{\phi}). \qquad (3.2) $$

We ask for the indulgence of the reader for the slight abuse of notation, which will pay off in a lot less hats. We will also often drop the explicit $\Delta_{\phi}$ dependence and use $\omega(\Delta)$, $\omega(\hat{\Delta})$.

The two basic properties we demand of a linear functional are:

### Finiteness

$$ \begin{align} \omega(\hat{\Delta}) < \infty \quad &\text{for } \hat{\Delta} \ge 0 \nonumber \\ \omega(\Delta) < \infty \quad &\text{for } \Delta \ge 0 \tag{3.3} \\ \omega[G] < \infty \quad &\text{for } \Delta \ge 0 \nonumber \end{align} $$

### Swapping

$$ \omega[G] = \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \omega(\hat{\Delta}|\Delta_{\phi}) = \sum_{\Delta} (a\lambda)_{\Delta} \omega(\Delta|\Delta_{\phi}). \qquad (3.4) $$

These conditions are broad enough that they allow for many different sets of functionals. Clearly the original crossing constraints are recovered by considering the (trivial) set of functionals which evaluate a function at a specific point. A more interesting basis is the set of derivatives at $z = \frac{1}{2}$, which is the approach used in the numerical bootstrap [6]. In this paper our main focus is on an interesting set of functionals which in a sense are dual to the generalized free field solutions to crossing discussed at the end of the previous section. However, to illustrate the main ideas we will first consider a toy basis of functionals in the next subsection.

For now, let us introduce our general ansatz for the functionals and examine in detail how it is constrained by the above requirements. A natural starting point is:

$$ \omega[\mathcal{F}] := \int_{-\infty}^{0} dz h_{-}(z) \mathrm{Im}[\mathcal{F}(z)] + \int_{1}^{\infty} dz h_{+}(z) \mathrm{Im}[\mathcal{F}(z)] \qquad (3.5) $$

with $\mathrm{Im}[\mathcal{F}(z)] := \lim_{\epsilon \to 0^+} \frac{1}{2i} [\mathcal{F}(z + i\epsilon) - \mathcal{F}^*(z - i\epsilon)]$. It is very natural to consider such an ansatz, since we will be acting on functions which are holomorphic in the crossing region $\mathcal{R}$, and such functions are essentially determined by their values near the cuts $(-\infty, 0) \cup (1, \infty)$. We can think of the above as integrating the kernels $h_+$, $h_-$ against test functions on a contour that wraps these cuts. Assuming for now sufficiently soft behaviour at $z = \infty$ as well as holomorphicity of $h_{\pm}$ away from the real axis, we can deform the contour of integration to obtain instead:

$$ \omega[\mathcal{F}] = \int_{\frac{1}{2}+i\mathbb{R}} \frac{dz}{2\pi i} [h_{+}(z) - h_{-}(z)]\mathcal{F}(z) \\ + \int_{\frac{1}{2}}^{1} \frac{dz}{\pi} \mathrm{Im}[h_{+}(z)]\mathcal{F}(z) - \int_{0}^{\frac{1}{2}} \frac{dz}{\pi} \mathrm{Im}[h_{-}(z)]\mathcal{F}(z). \qquad (3.6) $$
<figure><img src="image_1.png" /><figcaption>Figure 1. Contour deformation leading from (3.5) to (3.6).</figcaption></figure>

This deformation can be visualized in figure 1. One should really think of this form of the functional action as the more fundamental one. In particular, when deforming the contour to go back to (3.5) we might not be allowed to wrap the contours around the branch points at $z = 0, 1$.

From this expression we find that finiteness of the action on the conformal blocks (2.7), (2.11) demands:

• Finiteness:

$$h_+(z) - h_-(z) \stackrel{|z| \to \infty}{\sim} O(|z|^{\Delta_{\phi} - \varepsilon - 1 - \eta}), \qquad (3.7a)$$

$$\mathrm{Im}[h_+(z)] \stackrel{z \to 1^-}{\sim} O[(1-z)^{\Delta_{\phi}-1+\eta}], \qquad (3.7b)$$

$$\mathrm{Im}[h_-(z)] \stackrel{z \to 0^+}{\sim} O[(-z)^{\Delta_{\phi}-1+\eta}]. \qquad (3.7c)$$

for some $\eta > 0$. Now consider swapping. Since the crossing sums are uniformly convergent in the crossing region, by deforming the contour of integration we see that we need only worry about the behaviour of the kernels in the regions near $\infty, 0$ and $1$. However, it is clear that for the last two no additional constraints come from the swapping condition. For instance, it is true that the infinite sum over boundary blocks behaves worse near $z = 1$ than any single term in the sum. However, we can use crossing to bound this behaviour by the bulk channel expansion, for which we are already assuming the functional action is finite. The same argument goes through mutatis mutandis for $z = 0$. So, we must worry only about the behaviour at infinity. Here, there is indeed a potential problem, since the tail of the infinite sum over blocks falls off more slowly near $z = \infty$ than any individual block.

Let us see this explicitly. By crossing, it is sufficient to consider one particular channel so let us pick for definiteness the boundary one. We will use the fact that the boundary block can be written in terms of a radial coordinate [32]:

$$\hat{G}_{\hat{\Delta}}(z) = \hat{G}_{\hat{\Delta}}^{\rho}(\rho_z) = (4\rho_z)^{\hat{\Delta}} {}_2F_1\left(\frac{1}{2}, \frac{1}{2}; \hat{\Delta}, \frac{1}{2} + \hat{\Delta}; \rho_z^2\right), \qquad \rho_z := \frac{z}{(1 + \sqrt{1 - z})^2}. \qquad (3.8)$$

We will abuse notation and denote the expression on the right by $\hat{G}_{\hat{\Delta}}^{\rho}(\rho_z|\Delta_{\phi})$. The important point is that everywhere on the crossing region, and in particular along the functional
kernel integration contour, we have $|\rho_z| < 1$. We then have:

$$ \left| \sum_{\hat{\Delta} > \hat{\Delta}^*} \mu_{\hat{\Delta}}^2 \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi}) \right| \le \sum_{\hat{\Delta} > \hat{\Delta}^*} \mu_{\hat{\Delta}}^2 \frac{\hat{G}_{\hat{\Delta}}^{\rho}(|\rho_z|)}{|z|^{\Delta_{\phi}}} = \sum_{\hat{\Delta} > \hat{\Delta}^*} \mu_{\hat{\Delta}}^2 \frac{\hat{G}_{\hat{\Delta}}(z_e)}{|z|^{\Delta_{\phi}}} \quad (3.9) $$

where we have introduced an effective z variable:

$$ z_e = \frac{4|\rho_z|}{(1 + |\rho_z|)^2}, \quad z_e \underset{|z| \to \infty}{\sim} 1 - \frac{1}{|z|} \quad (3.10) $$

Hence we see that in the above the large z limit is dominated by the bulk channel identity. Therefore we can bound:

$$ \sum_{\hat{\Delta} > \hat{\Delta}^*} \mu_{\hat{\Delta}}^2 \frac{\hat{G}_{\hat{\Delta}}(z_e)}{|z|^{\Delta_{\phi}}} \le \frac{C|z_e|^{\Delta_{\phi}}}{|z|^{\Delta_{\phi}}(1-z_e)^{\Delta_{\phi}}} \underset{|z| \to \infty}{=} O(1) \quad (3.11) $$

This should be compared with the behaviour of an individual block, which goes as $\sim |z|^{-\Delta_{\phi}}$. The contribution from the tail will be negligible as long as we make the contribution from the integration region near $z = \infty$ vanishingly small. Given the bound above, we see that this actually requires a stronger falloff condition on $h_+(z) - h_-(z)$ than (3.7).

The conclusion is that overall, finiteness and swapping conditions constrain our ansatz in the following way:

### Finiteness+Swapping

$$ h_+(z) - h_-(z) \underset{|z| \to \infty}{\sim} O(|z|^{-1-\eta}), \quad (3.12a) $$

$$ \mathrm{Im}[h_+(z)] \underset{z \to 1^-}{\sim} O[(1-z)^{\Delta_{\phi}-1+\eta}], \quad (3.12b) $$

$$ \mathrm{Im}[h_-(z)] \underset{z \to 0^+}{\sim} O[(-z)^{\Delta_{\phi}-1+\eta}]. \quad (3.12c) $$

The first condition above is the most non-trivial. Relaxing or strengthening it allows for functionals that act on objects with softer or harder behaviour near $z = \infty$. In particular, perturbative corrections to correlators may not share the same large z-behaviour as fully-fledged correlators.

## 3.2 Warm-up: a simple basis

### 3.2.1 Special values of $\Delta_{\phi}$ and orthonormal functionals

We will now consider a simple basis of functionals which nevertheless has several nice properties. In this basis, the functionals will be dual to blocks in each channel with specific integer-spaced scaling dimensions that appear in the generalized free solutions. We will take a pedagogical route to explain general features which will be useful later, but the final result will be extremely simple.

The construction starts from the simple observation that the discontinuities of both bulk and boundary blocks are made up of two pieces. One of the pieces contains a sine
function, which vanishes automatically for those values of $\Delta, \hat{\Delta}$ corresponding to the generalized free field solutions. The other piece is, for both cases, a Jacobi function of the first kind, which reduces to a polynomial when $\hat{\Delta} = 1 + \varepsilon + n$ and $\Delta = 2 + 2\varepsilon + 2n$ with $n$ integer. This suggests using the orthogonality properties of Jacobi polynomials to construct the following functionals:

$$ \tau_n^{\text{pre}}: \quad h_+(z) = 0, \quad h_-(z) \propto (1-z)^{\Delta_\phi - 2-\varepsilon} P_n^{(0,\varepsilon)}\left(\frac{1+z}{1-z}\right) \quad (3.13) $$

$$ \hat{\tau}_n^{\text{pre}}: \quad h_-(z) = 0, \quad h_+(z) \propto z^{\Delta_\phi - 2 - \varepsilon} P_n^{(-\varepsilon, \varepsilon)}\left(\frac{2-z}{z}\right). \quad (3.14) $$

Here the superscript stands for 'pre-functionals'. We call these objects pre-functionals, since we have not yet imposed on them the finiteness and swapping properties. With a suitable choice of normalization factor we have:

$$ \tau_n^{\text{pre}}(\Delta_m) = \delta_{nm}, \qquad \tau_n^{\text{pre}}(\hat{\Delta}_m) = 0 \quad (3.15a) $$

$$ \hat{\tau}_n^{\text{pre}}(\Delta_m) = 0, \qquad \hat{\tau}_n^{\text{pre}}(\hat{\Delta}_m) = \delta_{nm} \quad (3.15b) $$

with $\hat{\Delta}_n = 1 + \varepsilon + n$ and $\Delta_m = 2 + 2\varepsilon + 2n$ as before. In particular, we see that as long as $\Delta_\phi = \varepsilon + 1 + k$ with $k \ge 0$ some integer, these functionals are dual to those blocks appearing in the crossing equation for the generalized free field solutions.

In general we have more prefunctionals than we need, in the sense that some of these are not dual to any blocks appearing in the decomposition of the generalized free two point function. For $\Delta_\phi = \varepsilon + 1 + k$ with $k \ge 0$, we see that the bulk channel prefunctionals with $n = 0, 1, \dots, k-1$ are not required. As for the boundary channel, we do not need the functionals with $n = 0, 1, \dots, k-1$ for Neumann boundary condition and in the Dirichlet case those with $n = 0, 1, \dots, k$. These unwanted functionals can be put to good use to construct bona fide functionals (satisfying the swapping property) from the prefunctionals above.

The construction is clearer to understand with an example. Let us set $\Delta_\phi = \varepsilon + 1$. The prefunctionals $\hat{\tau}_n$ have kernels $h_+(z)$ that do not falloff sufficiently fast at infinity to satisfy conditions (3.12). Hence we must proceed with caution. There is a simple thing to do which is merely to consider linear combinations of prefunctionals that do falloff faster. One particular choice is to set

$$ \hat{\tau}_n = \hat{\tau}_n^{\text{pre}} - c_n \hat{\tau}_0^{\text{pre}} \quad (3.16) $$

for some suitable coefficient $c_n$, so as to improve the falloff at infinity from the corresponding kernels from $1/z$ to $1/z^2$. It is clear many other choices of subtractions will do, but this is canonical in a sense, since recall that we can definitely afford to lose $\hat{\tau}_0^{\text{pre}}$ for the Dirichlet boundary condition. For Neumann it looks like we are apparently missing a functional, but there is a good reason for this, as will be explained further below.

After this subtraction procedure we have perfectly good boundary functionals. We can proceed in the same way for the bulk channel functionals, defining

$$ \tau_n = \tau_n^{\text{pre}} - d_n \hat{\tau}_0^{\text{pre}} \quad (3.17) $$
with $d_n$ chosen such that $h_+(z) - h_-(z)$ decays like $1/z^2$ at infinity. In this way we construct a set of functionals satisfying finiteness, swapping and the duality conditions:

$$ \begin{aligned} \tau_n(\Delta_m) &= \delta_{nm}, & \tau_n(\hat{\Delta}_m) &= -d_n\delta_{m,0}, & n, m \ge 0 \\ \hat{\tau}_n(\Delta_m) &= 0, & \hat{\tau}_n(\hat{\Delta}_m) &= \delta_{n,m} - c_n\delta_{m,0}, & n \ge 1, m \ge 0. \end{aligned} \qquad (3.18) $$

where now $\Delta_m = 2\Delta_\phi + 2m$, $\hat{\Delta}_n = \Delta_\phi + n$. This procedure can be repeated for other values of $\Delta_\phi = 1 + \epsilon + k$ in a similar manner. Everytime we increase $\Delta_\phi$ by one unit the asymptotic behaviour of the prefunctionals is worse by an extra factor of $z$, but this is compensated by there being more prefunctionals to use. Not all of these are independent, and in the end the counting of degrees of freedom always works out that we end up with just enough functionals to obtain the duality relations written above.

### 3.2.2 Shifted functionals and general $\Delta_\phi$

Let us go back to our simple example at $\Delta_\phi = 1 + \epsilon$. In this case it is easy to check that for $n=1$ one gets:

$$ \hat{\tau}_1: \qquad h_+(z) \propto \frac{1}{z^2}. \qquad (3.19) $$

In fact, one can check that for any $\Delta_\phi = 1 + \epsilon + k$ one always finds this very same result. The reason for this is easy to understand. The Jacobi polynomials can be written as polynomials in $z^{-1}$, and so by taking linear combinations it must be that a particularly simple set of functionals satisfying finiteness and swapping can be defined as:

$$ s\hat{\tau}_n: \qquad h_-(z) = 0, \qquad h_+(z) \propto \frac{1}{z^{1+n}}, \qquad n \ge 1. \qquad (3.20) $$

The letter ‘s’ stands for “simple” or “shifted”, since we see all kernels are simply determined from the $n=1$ case by shifting the exponent or multiplication by a power of $z^{-1}$.

The shifted functionals satisfy the duality conditions:

$$ s\hat{\tau}_n(\Delta_m) = 0 \qquad \text{for } m \ge 0, \quad n \ge 1, \qquad (3.21a) $$

$$ s\hat{\tau}_n(\hat{\Delta}_m) = \delta_{nm}, \qquad \text{for } m \ge n, \quad n \ge 1, \qquad (3.21b) $$

where $\hat{\Delta}_m = \Delta_\phi + m = \epsilon + 2 + m$ and $\Delta_m = 2\Delta_\phi + 2m$.

The discussion goes through mutatis mutandis for the bulk channel, for which we can choose:

$$ s\tau_n: \qquad h_+(z) = 0, \qquad h_-(z) \propto \frac{1}{(1-z)^{1+n}}, \qquad n \ge 1 \qquad (3.22) $$

$$ s\tau_0: \qquad h_+(z) \propto -\frac{1}{z}, \qquad h_-(z) \propto \frac{1}{(1-z)}. \qquad (3.23) $$

where for $s\tau_0$ we must choose proportionality factors that $h_+(z) - h_-(z)$ falls off as $1/z^2$ at infinity. These functionals now satisfy

$$ s\tau_n(\Delta_m) = \delta_{nm}, \qquad \text{for } m \ge n, \quad n \ge 0, \qquad (3.24a) $$

$$ s\tau_n(\hat{\Delta}_m) = -\delta_{m,0}e_m, \qquad \text{for } m, n \ge 0, \qquad (3.24b) $$

for some constants $e_m$.
This particular functional basis has extremely simple kernels $h_{\pm}(z)$, and for this we paid the relatively small price that the orthonormality conditions are more involved. It is clear however that starting from this basis, for any $n$ a finite Gram-Schmidt decomposition can be made to go back to the orthonormal basis.

Now comes the crucial point. Although we have arrived at this basis via a somewhat circuitous route, the final functionals are extremely simple, and in particular all details related to the fact that we started with Jacobi polynomials have disappeared. This is not an accident, since the functionals constructed above will actually have the same properties for any $\Delta_{\phi}$.

To see this, consider for definiteness the action of the boundary functionals. The bulk case is analogous. When acting on bulk blocks we get

$$ s\hat{\tau}_n(\Delta) \propto \int_1^\infty \frac{dz}{z^{1+n}} \mathrm{Im}[G_\Delta(z|\Delta_\phi)] \\ = \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_\phi)\right] \int_1^\infty \frac{dz}{z^{1+n}} (z-1)^{\frac{\Delta}{2}-\Delta_\phi} \frac{2F_1\left(\frac{\Delta}{2}, \frac{\Delta}{2}, \Delta - \varepsilon, \frac{z-1}{z}\right)}{z^{\frac{\Delta}{2}}} \quad (3.25) $$

Thanks to the oscillating factor we see the functional is automatically zero for $\Delta = 2\Delta_{\phi}+m$ independently of $\Delta_{\phi}$, with $m \ge 0$.<fnref n="6" /> Now consider the action on the boundary blocks:

$$ s\hat{\tau}_n(\hat{\Delta}) \propto \int_1^\infty \frac{dz}{z^{1+n}} \mathrm{Im}[\hat{G}_{\hat{\Delta}}(z|\Delta_\phi)]. \quad (3.26) $$

Starting from this expression we deform the contour of integration to wrap around the discontinuity of $\hat{G}_{\hat{\Delta}}$ for negative $z$. We will see momentarily what are the conditions for this manipulation. If we do this we now get

$$ s\hat{\tau}_n(\hat{\Delta}) \propto \sin\left[\pi(\hat{\Delta} - \Delta_\phi)\right] \int_{-\infty}^{0} \frac{dz}{z^{1+n}} \frac{\hat{G}_{\hat{\Delta}}\left(\frac{z}{z-1}|\Delta_\phi\right)}{(1-z)^{\Delta_\phi}} \quad (3.27) $$

We see that automatically the functional vanishes on boundary blocks with dimensions $\hat{\Delta} = \Delta_{\phi} + m$, with $m$ integer. But this is true only down to $m = n$, since at this point the integral develops a divergence from the region near $z = 0^{-}$. The divergence cancels the zero of the sine function, leading to a finite answer. For $\hat{\Delta}$ smaller than $\Delta_{\phi} + n$ we can no longer trust the contour deformation and must go back to the original definition of the functional action. We conclude that $s\hat{\tau}_n$ satisfies indeed the duality relations (3.21). A similar argument goes through for the bulk functional basis. The messy details involving $\Delta_{\phi}$ and $\varepsilon$ only become important if we insist in having a fully orthonormal basis.

### 3.3 Functional bootstrap equations

To summarize, we have shown that an orthonormal basis of functionals $\tau_n, \hat{\tau}_n$ exist with the properties:

$$ \begin{aligned} \tau_n(\Delta_m) &= \delta_{nm}, & \tau_n(\hat{\Delta}_m) &= -d_n\delta_{m,0}, & n \ge 1, m \ge 0 \\ \hat{\tau}_n(\Delta_m) &= 0, & \hat{\tau}_n(\hat{\Delta}_m) &= \delta_{nm} - c_n\delta_{m,0}, & n, m \ge 0 \end{aligned} \quad (3.28) $$

<footnote n="6">⁶For $\Delta \le 2\Delta - 2$ the integral diverges, and we must use instead expression (3.6) for the functional action.</footnote>
with $\hat{\Delta}_m = \Delta_\phi + m$ and $\Delta_m = 2\Delta_\phi + 2m$. This basis can be obtained for any $\Delta_\phi$ by starting with the shifted functionals $s\tau_n$ and $s\hat{\tau}_n$ and performing a finite Gram-Schmidt orthonormalization procedure. Applying this functional basis to a generic crossing equation this gives a set of necessary conditions on the BCFT data:

$$ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \omega_n(\hat{\Delta}) = \sum_{\hat{\Delta}} (a\lambda)_{\Delta} \omega_n(\Delta), \quad \omega = \tau, \hat{\tau}, \qquad (3.29) $$

which are the functional bootstrap equations associated to this particular functional basis.

There is a different way to understand the origin of these equations. The orthonormality conditions suggest one may write

$$ G_{\Delta}(z) = \sum_{n=0}^{\infty} \left[ \tau_n(\Delta) G_{2\Delta_\phi+2n}(z) + \hat{\tau}_n(\Delta) \hat{G}_{1+\Delta_\phi+n}(z) \right] \qquad (3.30a) $$

$$ \hat{G}_{\hat{\Delta}}(z) = \sum_{n=0}^{\infty} \left[ \tau_n(\Delta) G_{2\Delta_\phi+2n}(z) + \hat{\tau}_n(\Delta) \hat{G}_{1+\Delta_\phi+n}(z) \right] \qquad (3.30b) $$

where we have omitted the dependence of the blocks (and functionals) on $\Delta_\phi$ for simplicity. These expressions can be thought of as decompositions of generic bulk and boundary blocks into a basis, the basis being constituted by bulk and boundary blocks with specific, generalized free, scaling dimensions. Again, we use the word basis loosely here: roughly speaking we have in mind the vector space made up of linear combinations of bulk blocks and boundary blocks.<fnref n="7" />

If these decompositions are actually true (we will not attempt to prove this here), then by plugging them into the original crossing equation and assuming we can commute sums we get

$$ \sum_{n=0}^{\infty} \left\{ \left[ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \tau_n(\hat{\Delta}) - \sum_{\hat{\Delta}} (a\lambda)_{\Delta} \tau_n(\Delta) \right] G_{2\Delta_\phi+2n}(z) + \left[ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \hat{\tau}_n(\hat{\Delta}) - \sum_{\hat{\Delta}} (a\lambda)_{\Delta} \hat{\tau}_n(\Delta) \right] \hat{G}_{1+\Delta_\phi+n}(z) \right\} = 0, \qquad (3.31) $$

and the functional equations now follow demanding that the coefficient of each “basis element” is vanishing.

A closely related perspective on the same equations allow us to define a toy version of the Polyakov bootstrap. If we define the (toy) Polyakov blocks

$$ \begin{aligned} P_{\Delta}(z) &:= G_{\Delta}(z) - \sum_{n=0}^{\infty} \tau_n(\Delta) G_{2\Delta_\phi+2n}(z), \\ \hat{P}_{\hat{\Delta}}(z) &:= \hat{G}_{\hat{\Delta}}(z) - \sum_{n=0}^{\infty} \hat{\tau}_n(\hat{\Delta}) \hat{G}_{1+\Delta_\phi+n}(z), \end{aligned} \qquad (3.32) $$

<footnote n="7">⁷A more precise statement would be that we would have a basis if the dual functionals $\tau_n$, $\hat{\tau}_n$ would lead to a set of constraints which fully capture those of the crossing equation.</footnote>
then it is not hard to show using the basis decompositions above that the functional bootstrap equations arise from demanding that correlators may be expressed as

$$ \mathcal{G}(z) = \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \hat{P}_{\hat{\Delta}}(z) + \sum_{\Delta} (a\lambda)_{\Delta} P_{\Delta}(z). \qquad (3.33) $$

Hence, this simple basis have allowed us to derive a toy version of the Polyakov bootstrap for BCFTs. Even though we have not attempted to, it is likely all the above statements can be made rigorous. In the next sections we will construct a different basis of functionals, leading to different functional bootstrap equations and different Polyakov blocks, but the basic ideas are as above.

The functionals equations can be used for instance to reconstruct the generalized free solutions to crossing. Assuming the correct spectrum for each boundary condition, one can read off the BCFT data:

$$ \begin{aligned} \mu_{\Delta_{\phi}+1+2p}^2 &= \hat{\tau}_{1+2p}(0), & (a\lambda)_{2\Delta_{\phi}+2p} &= -\tau_p(0), & \text{Dirichlet} \\ \mu_{\Delta_{\phi}+2p}^2 &= \hat{\tau}_{2p}(0) + c_{2p}\mu_{\Delta_{\phi}}^2, & (a\lambda)_{2\Delta_{\phi}+2p} &= -\tau_p(0) - d_p\mu_{\Delta_{\phi}}^2, & \text{Neumann} \end{aligned} \qquad (3.34) $$

Note we also have the constraints:

$$ \begin{aligned} \hat{\tau}_{2p}(0) &= 0, & \text{Dirichlet} \\ \hat{\tau}_{1+2p}(0) &= -c_{1+2p}\mu_{\Delta_{\phi}}^2, & \text{Neumann} \end{aligned} \qquad (3.35) $$

which follow from applying these functionals to the Dirichlet and Neumann solutions respectively. From these relations we can reconstruct the whole OPE and BOE data for both solutions given $c_n, d_n$ which were fixed by demanding the correct asymptotics of the functional kernels.

It may seem puzzling that we are apparently missing a functional, namely $\hat{\tau}_0$, which we had to remove to fix the behaviour of the kernels at $z = \infty$. However, this is precisely how it must be. Indeed, notice that by combining the above two solutions defining

$$ \mathcal{G}(z) = 2x - 1 + \frac{1}{(1-z)^{\Delta_{\phi}}} = x\mathcal{G}^{N}(z) + (1-x)\mathcal{G}^{D}(z), \qquad (3.36) $$

with $0 \le x \le 1$, leads to a perfectly reasonable family of solutions to crossing which generically contain all operators in both GFF solutions. Had we been able to find a new functional such that the duality conditions would hold even for $\hat{\Delta} = \Delta_{\phi}$, i.e. with $c_n, d_n = 0$. Then we could read off the OPE and BOE data simply by acting with the functionals to obtain

$$ \mu_{\hat{\Delta}_n}^2 = \hat{\tau}_n(0), \quad (a\lambda)_{\Delta_n} = \tau_n(0) \qquad (3.37) $$

Since on the r.h.s. of these equations we just have some fixed numbers, this would be in contradiction with the existence of the family of solutions described above. The fact that we have this degree of freedom in the choice of crossing solution which keeps the set $\hat{\Delta}_m = \Delta_{\phi} + m$ and $\Delta_m = 2\Delta_{\phi} + 2m$ fixed in the block decomposition, translates into the absence of a functional in the dual space. Indeed, $x$ controls the $O(1)$ behaviour of the
<figure><img src="image_2.png" /><figcaption>**Figure 2.** The bulk exchange, boundary exchange and contact Witten diagrams respectively. The semicircular area represents $AdS_{d+1}$ bulk in the presence of a brane, which is represented by the dashed line.</figcaption></figure>

correlator at $\infty$, and our functionals were carefully constructed such that they do not care about this overall constant.

We could now go on to investigate the consequences that these equations imply for generic solutions to crossing. However, it seems likely this will not be very useful. Consider for instance starting from the Dirichlet GFF solution and deforming it in a very slight manner. Unless this deformation is trivial like in the $x$- dependent family of correlators above, one should expect that at least the boundary spectrum should change non-trivially. While to zeroth order the equations are nicely diagonal in the spectrum thanks to orthonormality, allowing us to just read off the OPE and BOE data, the minute we expand in the spectrum the equations will become completely coupled. This is because $\partial_{\Delta}\omega_n(\Delta_m) \neq \delta_{nm}$ for these functionals. Clearly it would be nice to have a basis that has these properties, since it would allow us to do perturbation theory in an efficient manner. Such a basis indeed exists, and we will construct it in the next section

Before we do so, we would just like to comment that in [28], the authors were able to bootstrap the Wilson-Fisher fixed point to $O(\epsilon^2)$ by essentially using the basis of functionals constructed in this section, although they did not formulate their approach in this language.<fnref n="8" /> The reason this was possible is that the full set of boundary and bulk operators appear only at $O(\epsilon^2)$ and only a finite set of them appear with dimensions away from their free integer values. It is clear however that this approach is doomed to fail at order $O(\epsilon^3)$ where an infinite number of operators pick up anomalous dimensions. With the basis that we will construct next, there are no such limitations and it is possible in principle to go to arbitrarily high orders in perturbation theory.

# 4 Bootstrap equations and BCFT basis

## 4.1 Witten diagrams and functional equations

As we have discussed in section 2, the generalized free solutions to crossing can be understood as arising from the theory of a free massive scalar field in $AdS_{d+1}$ in the presence of a radially extended brane, on which Neumann or Dirichlet boundary conditions are imposed. In such a theory it is very natural to consider perturbations. We will focus on those terms

<footnote n="8">⁸To be precise they used our initial Jacobi polynomial basis and only for the bulk channel.</footnote>
in perturbation theory which modify how the field interacts with the brane, while keeping the bulk operator dimensions fixed. At tree level there are three sets of perturbations that we can consider. We represent them diagramatically in figure 2. They correspond to bulk and boundary exchange Witten diagrams, as well as contact interactions. Such diagrams have both bulk and boundary conformal block expansions. Crossing symmetry for each such diagram leads to the following equations:

$$G_{\Delta}(z) = \sum_{n=0}^{\infty} \left[ a_n(\Delta) \hat{G}_{\hat{\Delta}_n}(z) + b_n(\Delta) \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n}(z) + t_n(\Delta) G_{\Delta_n}(z) \right] \quad (4.1a)$$

$$\hat{G}_{\hat{\Delta}}(z) = \sum_{n=0}^{\infty} \left[ \hat{a}_n(\hat{\Delta}) \hat{G}_{\hat{\Delta}_n}(z) + \hat{b}_n(\hat{\Delta}) \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n}(z) + \hat{t}_n(\hat{\Delta}) G_{\Delta_n}(z) \right] \quad (4.1b)$$

$$0 = \sum_{n=0}^{\infty} \left[ d_n^{(k)} \hat{G}_{\hat{\Delta}_n}(z) + e_n^{(k)} \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n}(z) + f_n^{(k)} G_{\Delta_n}(z) \right] \quad (4.1c)$$

with $\Delta_n = 2\Delta_\phi + 2n$ running over bulk double trace operators, $\hat{\Delta}_n = \frac{1-\nu}{2} + \Delta_\phi + 2n$ over boundary derivative operators and $k$ labels the contact interaction. We have omitted the dependence of the blocks on $\Delta_\phi$ for brevity.

By analogy with the results of the last section, these expressions suggest that it may be possible to use the sets $G_{\hat{\Delta}_n}$, $\partial_{\hat{\Delta}} G_{\hat{\Delta}_n}$ and $G_{\Delta_n}$ (modulo some elements to account for the contact interactions above) as a basis. Notice that unlike the toy basis of the previous section, here the boundary operator dimensions have spacing of two units, but this is compensated by the extra terms involving derivatives of these blocks. Our job in this section is then to construct functionals which are dual to this basis, i.e. $\alpha_n, \beta_n, \theta_n$ such that:

$$\alpha_n(\hat{\Delta}_m) = \delta_{nm}, \qquad \partial_{\hat{\Delta}} \alpha_n(\hat{\Delta}_m) = 0, \qquad \alpha_n(\Delta_m) = 0, \quad (4.2a)$$

$$\beta_n(\hat{\Delta}_m) = 0, \qquad \partial_{\hat{\Delta}} \beta_n(\hat{\Delta}_m) = \delta_{nm}, \qquad \beta_n(\Delta_m) = 0, \quad (4.2b)$$

$$\theta_n(\hat{\Delta}_m) = 0, \qquad \partial_{\hat{\Delta}} \theta_n(\hat{\Delta}_m) = 0, \qquad \theta_n(\Delta_m) = \delta_{nm} \quad (4.2c)$$

To be precise, there are really two different bases one can consider, depending on whether we choose functionals dual to the Neumann or Dirichlet solution. When necessary we will make this explicit by adding a superscript $+$, $-$ to the functionals for Neumann, Dirichlet respectively. We will see that in the latter case we will be able to indeed satisfy the above conditions, while for Neumann there will be a slight modification, as for the toy basis of the previous section, cf. equation (3.28).<fnref n="9" />

By acting with the functional bases on bulk and boundary blocks, one expects that decompositions of the form (4.1) should exist, with:

$$a_n(\Delta) \rightarrow \alpha_n(\Delta), \qquad b_n(\Delta) \rightarrow \beta_n(\Delta), \qquad t_n(\Delta) \rightarrow \theta_n(\Delta) \quad (4.3)$$

$$\hat{a}_n(\hat{\Delta}) \rightarrow \alpha_n(\hat{\Delta}), \qquad \hat{b}_n(\hat{\Delta}) \rightarrow \beta_n(\hat{\Delta}), \qquad \hat{t}_n(\hat{\Delta}) \rightarrow \theta_n(\hat{\Delta}) \quad (4.4)$$

<footnote n="9">⁹This is also very similar to what happens in the 1D bootstrap case, with our Dirichlet and Neumann cases being analogous to the fermionic and bosonic basis [19].</footnote>
In other words, these functionals should allow us to bootstrap Witten diagrams in $AdS_{d+1}$. We will compute these Witten diagrams explicitly in the next section, and study the precise relation to the functional bases in section 6.

Just like with the toy basis, we can act with these functionals on the crossing equation (2.12), to obtain a set of functional bootstrap equations:

$$ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \alpha_n(\hat{\Delta}) = \sum_{\Delta} (a\lambda)_{\Delta} \alpha_n(\Delta) \qquad (4.5a) $$

$$ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \beta_n(\hat{\Delta}) = \sum_{\Delta} (a\lambda)_{\Delta} \beta_n(\Delta) \qquad (4.5b) $$

$$ \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \theta_n(\hat{\Delta}) = \sum_{\Delta} (a\lambda)_{\Delta} \theta_n(\Delta) \qquad (4.5c) $$

for all $n \in \mathbb{N}$.<fnref n="10" /> The functional bootstrap equations are a set of fully non-perturbative constraints on the OPE and BOE data. Unlike the toy basis, the $\alpha, \beta, \theta$ basis now provide an ideal starting point for a perturbative expansion around the generalized free solution.

To see this explicitly, consider starting with the free correlator and deforming it by adding a new bulk operator with some small coefficient $g$. Expanding all relevant quantities,

$$ \mathcal{G}(z) = \mathcal{G}^{(0)}(z) + g \mathcal{G}^{(1)}(z) + O(g^2), \qquad (4.6a) $$

$$ \Delta_n = \Delta_n^{(0)} + g\Delta_n^{(1)} + O(g^2), \qquad \hat{\Delta}_n = \hat{\Delta}_n^{(0)} + 2n + g\hat{\Delta}_n^{(1)} + O(g^2) \qquad (4.6b) $$

$$ (a\lambda)_{\Delta_n} = t^{(0)} + g t^{(1)} + O(g^2), \qquad \mu_{\hat{\Delta}_n}^2 = a_n^{(0)} + g a_n^{(1)} + O(g^2) \qquad (4.6c) $$

where $t^{(0)}, a_n^{(0)}$ are the generalized free BCFT data given in (2.14), $\Delta_n^{(0)} = 2\Delta_\phi + 2n$, $\hat{\Delta}_n^{(0)} = \frac{1-\nu}{2} + 2\Delta_\phi + 2n$ and

$$ \mathcal{G}^{(0)}(z) = G_0(z) + \sum_{n=0}^{\infty} t_n^{(0)} G_{\Delta_n^{(0)}}(z) = \sum_{n=0}^{\infty} a_n^{(0)} \hat{G}_{\hat{\Delta}_n^{(0)}}(z) \qquad (4.7) $$

$$ \mathcal{G}^{(1)}(z) = G_{\Delta}(z) + \sum_{n=0}^{\infty} \left[ t_n^{(1)} G_{\Delta_n^{(0)}}(z) + \Delta_n^{(1)} \partial_{\Delta} G_{\Delta_n^{(0)}}(z) \right] \\ = \sum_{n=0}^{+\infty} \left[ a_n^{(1)} \hat{G}_{\hat{\Delta}_n^{(0)}}(z) + a_n^{(0)} \hat{\Delta}_n^{(1)} \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n^{(0)}}(z) \right] \qquad (4.8) $$

Acting with the functionals we read off:

$$ a_n^{(0)} = \alpha_n(0), \qquad a_n^{(1)} = \alpha_n(\Delta) - \sum_{n=0}^{+\infty} \Delta_n^{(1)} \partial_{\Delta} \alpha_n(\Delta_n^{(0)}) \qquad (4.9a) $$

$$ 0 = \beta_n(0), \qquad a_n^{(0)} \hat{\Delta}_n^{(1)} = \beta_n(\Delta) - \sum_{n=0}^{+\infty} \Delta_n^{(1)} \partial_{\Delta} \beta_n(\Delta_n^{(0)}) \qquad (4.9b) $$

$$ t_n^{(0)} = -\theta_n(0), \qquad t_n^{(1)} = -\theta_n(\Delta) - \sum_{n=0}^{+\infty} \Delta_n^{(1)} \partial_{\Delta} \theta_n(\Delta_n^{(0)}) \qquad (4.9c) $$

<footnote n="10">¹⁰Our definition of $\mathbb{N}$ is such that it includes zero.</footnote>
We see that all the BCFT data, including boundary operator anomalous dimensions may be determined from knowledge of the functional actions, as long as the quantities $\Delta_n^{(1)}$ are given as inputs: those anomalous dimensions belong purely to the bulk CFT and must be determined independently. A particular case corresponds to setting them to zero, whereupon we would recover the bulk exchange Witten diagram decomposition.<fnref n="11" /> It is not hard to see that this behaviour persists to arbitrary orders in perturbation theory, in other words, the $\alpha, \beta, \theta$ functional basis diagonalises perturbation theory around the generalized free solution as long as the bulk spectrum is kept fixed (or known). We will discuss a particular application in section 7, where we will bootstrap the Wilson-Fisher fixed point with a boundary with Neumann and Dirichlet boundary conditions.

After these motivating remarks, we will now turn to the actual construction of the functional basis. We begin by discussing general constraints on the functional kernels, before solving these constraints appropriately, first for the boundary channel functionals $\alpha, \beta$, and then the bulk functionals $\theta$.

## 4.2 Constraining the kernels

We start with the ansatz introduced in the previous section:

$$ \omega[\mathcal{F}] := \int_{-\infty}^{0} dz h_{-}(z) \mathrm{Im}[\mathcal{F}(z)] + \int_{1}^{\infty} dz h_{+}(z) \mathrm{Im}[\mathcal{F}(z)] \quad (4.10) $$

The orthonormality properties (4.2) imply that we want to find functionals which have double zeros in the boundary channel. So, unlike for our warm-up basis, it is much less obvious how to choose appropriate kernels. However, one of the lessons of that construction is that it is useful to rotate contours so as to go to the region where discontinuities of the blocks are simple.

Consider the action of a generic functional $\omega$ on the bulk blocks. Rotating contours we can write

$$ \omega(\Delta) = \int_{1}^{\infty} dz h_{+}(z) \mathrm{Im}[G_{\Delta}(z|\Delta_{\phi})] - \int_{0}^{1} dz \mathrm{Im}[h_{-}(z)] G_{\Delta}(z|\Delta_{\phi}) \\ - \int_{1}^{\infty} dz \{\mathrm{Im}[h_{-}(z)] \mathrm{Re}[G_{\Delta}(z|\Delta_{\phi})] + \mathrm{Re}[h_{-}(z)] \mathrm{Im}[G_{\Delta}(z|\Delta_{\phi})]\} \quad (4.11) $$

We want functionals to annihilate the generalized free boson bulk spectrum. In the region $z > 1$ we have that the real and imaginary parts of the bulk block are proportional to $\cos[\frac{\pi}{2}(\Delta - \Delta_{\phi})]$ and $\sin[\frac{\pi}{2}(\Delta - \Delta_{\phi})]$ respectively. Hence it is natural to keep only the latter terms in the expression above, which can be achieved by demanding

$$ \mathrm{Im}[h_{-}(z)] = 0, \quad \forall z \in \mathbb{R}. \quad (4.12) $$

This is then one of our conditions on the kernels. For future reference, we note that in this case the functional action becomes

$$ \omega(\Delta) = \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_{\phi})\right] \int_{1}^{\infty} dz [h_{+}(z) - h_{-}(z)] \frac{2F_{1}\left(\frac{\Delta}{2}, \frac{\Delta}{2}, \Delta - \varepsilon, \frac{z-1}{z}\right)}{(z-1)^{\Delta_{\phi}-\frac{\Delta}{2}} z^{\frac{\Delta}{2}}}. \quad (4.13) $$

<footnote n="11">¹¹It would be interesting to follow the deformation to second order, as we would expect then that bulk operators should be forced to pick up anomalous dimensions.</footnote>
This is valid only for $\Delta$ sufficiently large, depending on the singularity structure near $z = 1$ of $h_-, h_+$.

Next we look at the functional action on the boundary channel. In this case we want to do the contour rotation in a different way, swapping the roles of $h_+, h_-$ to obtain

$$ \omega(\hat{\Delta}) = \int_{-\infty}^{0} dz h_{-}(z) \mathrm{Im}[\hat{G}_{\hat{\Delta}}(z|\Delta_{\phi})] - \int_{0}^{1} dz \mathrm{Im}[h_{+}(z)] \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi}) \\ - \int_{-\infty}^{0} dz \left\{ \mathrm{Im}[h_{+}(z)] \mathrm{Re}[\hat{G}_{\hat{\Delta}}(z|\Delta_{\phi})] + \mathrm{Re}[h_{+}(z)] \mathrm{Im}[\hat{G}_{\hat{\Delta}}(z|\Delta_{\phi})] \right\} \quad (4.14) $$

We would like for the functional action to have double zeros. Since the imaginary and real part of the blocks oscillate as $\sin[\pi(\hat{\Delta} - \Delta_{\phi})]$ and $\cos[\pi(\hat{\Delta} - \Delta_{\phi})]$, we can achieve this by combining the non-oscillating and cosine pieces and setting the sine pieces to zero. That is to say, we demand:

$$ h_{-}(z) = \mathrm{Re}[h_{+}(z)], \qquad 0 > z \qquad (4.15a) $$

$$ \mathrm{Im}[h_{+}(z)] = -\nu \mathrm{Im} \left[ h_{+} \left( \frac{z}{z-1} \right) \right] (1-z)^{\Delta_{\phi}-2}, \qquad 0 < z < 1. \qquad (4.15b) $$

which constitutes a second set of constraints on the kernels. Under these conditions we get

$$ \omega(\hat{\Delta}) = -(1 - \nu \cos[\pi(\hat{\Delta} - \Delta_{\phi})]) \int_{0}^{1} dz \mathrm{Im}[h_{+}(z)] \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi}) \quad (4.16) $$

which will have generically double zeros for $\hat{\Delta} = \frac{1-\nu}{2} + \Delta_{\phi} + 2n$, again at least for $\hat{\Delta}$ sufficiently large that allows us to do the contour deformations.

In the following we will find solutions to equations (4.12), (4.15) subject to boundary conditions appropriate for $\alpha_n, \beta_n$ and $\theta_n$ functionals satisfying the orthonormality conditions (4.2), as well as the swapping and finiteness conditions (3.12).

Before we do this, there are important simplifications to keep in mind. Firstly, since $h_-(z)$ is real for real $z$, these equations as well as the finiteness and boundary conditions (3.12) only care about the difference $h_+(z) - h_-(z)$. Hence by shifting $h_+(z) \to h_+(z) + h_-(z)$ we can get rid of $h_-(z)$ altogether and work only with $h_+(z)$.

Secondly, and more importantly, these equations are invariant under the shifts:

$$ \begin{aligned} h_{+}(z) &\to z^{-p}(1-z)^{-q} h_{+}(z), \\ \Delta_{\phi} &\to \Delta_{\phi} - 2q - p, \\ \nu &\to (-1)^p \nu, \end{aligned} \qquad (4.17) $$

with $p, q \in \mathbb{Z}$. This means that if we know a solution with specified singular behaviour near $z = 0, z = 1$ it is easy to obtain all remaining solutions.

## 4.3 Boundary functionals

### 4.3.1 General solution

Boundary functionals are comprised of the $\alpha_n$ and $\beta_n$. The orthonormality conditions (4.2) tell us the corresponding functional actions have to be zero on all double trace blocks in the
$k \le (m+1)/2$ contain combinations of $\beta_1^+, \dots, \beta_{(m+1)/2}^+$, and for $k > (m+1)/2$ we get combinations of $\alpha_0^+, \dots, \alpha_{(m-1)/2}^+$ with the same $\beta$ functionals. In particular we find

$$ \beta_0^-: \quad h_+(z) = h_{0,0}(z) = \frac{\Gamma(\Delta_\phi + 1)}{\pi^{3/2}} z^{-3/2} {}_2\tilde{F}_1 \left(-\frac{1}{2}, \frac{3}{2}; \Delta_\phi + \frac{1}{2}; \frac{1}{z}\right) \quad (4.24a) $$

$$ \beta_1^+: \quad h_+(z) = h_{1,0}(z) = \frac{\Gamma(\Delta_\phi + 2)}{4\pi^{3/2}} z^{-3/2} {}_2\tilde{F}_1 \left(-\frac{3}{2}, \frac{5}{2}; \Delta_\phi + \frac{1}{2}; \frac{1}{z}\right) \quad (4.24b) $$

$$ \alpha_0^-: \quad h_+(z) = \frac{8}{(1 + \Delta_\phi)(2 + \Delta_\phi)} (h_{2,2}(z) - h_{2,0}(z) + n_{\Delta_\phi}h_{0,0}(z)) \quad (4.24c) $$

$$ \alpha_0^+: \quad h_+(z) = \frac{8}{\Delta_\phi(1 + \Delta_\phi)} (h_{1,0}(z) + h_{1,1}(z)) \quad (4.24d) $$

with

$$ n_{\Delta_{\phi}} = \frac{(\Delta_{\phi} + 1)(\Delta_{\phi} + 2)}{32} (1 + 8 \log(2) - 4 H_{\Delta_{\phi}}) \quad (4.25) $$

and $H_n$ the harmonic number function.

### 4.3.2 Neumann contact term ambiguity

There is a slight problem in the above. Although it is not clear from the notation, as we have defined it $\partial_\Delta(\alpha_0^+)(\Delta_\phi) \ne 0$ as required by (4.2), and it would only be by adding a suitable multiple of $\beta_0^+$ that we could make this vanish. But $\beta_0^+$ has not made an appearance in the set of kernels above. It should correspond to something like a kernel with $m = k = -1$ but this then does not have the correct asymptotic behaviour at infinity.

In fact, such a functional does not exist, and the reason is very similar to what happened for the easy basis in section 3.2. Here the reason why such a functional cannot exist, is that it would preclude the existence of a relevant contact interaction in $AdS_{d+1}$ between the bulk field and the brane, namely $\int_{AdS_d} d^dx \sqrt{g} \Phi^2(\vec{x}, x^\perp, y=0)$. This decays at large $z$ and is only possible for Neumann boundary conditions. Such a contact interaction deforms the Neumann GFF solution while introducing anomalous dimensions to boundary operators that decay sufficiently fast for large $\hat{\Delta}$.<fnref n="12" /> It is easy to check that had we been able to construct a full orthonormal basis of Neumann functionals, this contact term would've been ruled out, since it is a modification of a solution to crossing that does not introduce new states.

Hence, a fully orthonormal basis cannot exist for the Neumann case. Instead, the best we can do is to find a set of $\alpha_n, \beta_n$ with $\beta_0 \equiv 0$, such that

$$ \alpha_n^+(\hat{\Delta}_m) = \delta_{nm}, \quad \partial_\Delta \alpha_n^+(\hat{\Delta}_m) = -\delta_{m0}d_n, \quad \alpha_n^+(\Delta_m) = 0, \quad (4.26a) $$

$$ \beta_n^+(\hat{\Delta}_m) = 0, \quad \partial_{\hat{\Delta}} \beta_n^+(\hat{\Delta}_m) = \delta_{nm} - \delta_{m0}e_n, \quad \beta_n^+(\Delta_m) = 0, \quad (4.26b) $$

with $e_0 \equiv 1$. The parameters $d_n, e_n, f_n$ are functions $\Delta_\phi$ and spacetime dimension and follow from the construction of the functional kernels. They can be determined by explicit computation of the functional actions, as shown in appendix A. The constants $f_n$ will make their appearance in the next subsection when we construct bulk functionals.

<footnote n="12">¹²This is entirely analogous to what happens in the 1d functional construction when one attempts to construct a basis associated to the generalized free boson [19].</footnote>
Thanks to these modified orthonormality conditions, the functionals now allow for a homogeneous solution to crossing:

$$0 = \sum_{n=0}^{\infty} \left[ d_n \hat{G}_{\hat{\Delta}_n}(z|\Delta_\phi) + e_n \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n}(z|\Delta_\phi) + f_n G_{\Delta_n}(z|\Delta_\phi) \right]. \quad (4.27)$$

This is of course nothing but the contact interaction described above. This will be checked by direct computation of the Witten diagram in section 5.

To conclude this subsection, let us explain why we refer to this as an “ambiguity” in the choice of basis. The point is that, unlike for the Dirichlet case, for Neumann it is simply impossible to have a fully orthonormal set of basis functionals with the correct large z behaviour, and once we have given up on full orthonormality, there is no good reason to hold relations (4.26) to be sacrosanct. Equivalently there is now no canonical choice of basis, unlike for the Dirichlet functionals. In particular in the above we have decided to set $\beta_0 \equiv 0$, but it is clear we could have chosen differently, e.g. by shifting every functional by $\beta_k$ we could have set $\beta_k = 0$ instead, or we could have even set some finite linear combination of functionals to vanish. It is in this sense that the choice of basis is ambiguous. Nevertheless, the choice above is particularly simple and we shall stick to it.

### 4.3.3 Shifted functionals

We will now show a different way of constructing the functional basis, by using the solution-generating shifts introduced in (4.17).

Let us set:

$$\beta_0^{-} : \quad h_+(z) = h_{\beta_0^{-}}(z|\Delta_{\phi}), \quad (4.28a)$$

$$\alpha_0^{-} : \quad h_+(z) = h_{\alpha_0^{-}}(z|\Delta_{\phi}), \quad (4.28b)$$

$$\beta_1^{+} : \quad h_+(z) = h_{\beta_1^{+}}(z|\Delta_{\phi}), \quad (4.28c)$$

$$\alpha_0^{+} : \quad h_+(z) = h_{\alpha_0^{+}}(z|\Delta_{\phi}), \quad (4.28d)$$

i.e. the kernels on the righthand side can be read off from (4.24). Then we can define

$$s\beta_n^{-} : \quad h_+(z) = z^{-2n}h_{\beta_0^{-}}(z|\Delta_{\phi} + 2n) \quad (4.29)$$

$$s\alpha_n^{-} : \quad h_+(z) = z^{-2n}h_{\alpha_0^{-}}(z|\Delta_{\phi} + 2n) \quad (4.30)$$

These functionals satisfy all the desired constraints and have the correct singular behaviour near $z = 0$. By taking finite linear combinations of $\omega_n, \omega_{n-1}, \dots, \omega_0$ with $\omega = s\alpha^{-}, s\beta^{-}$ these can be chosen to satisfy the orthonormality conditions (4.2). This completes the construction of the Dirichlet boundary functionals. As for the Neumann case, first we can check that indeed:

$$h_{\beta_1^{+}}(z|\Delta_{\phi}) = z^{-1}h_{\beta_0^{-}}(z|\Delta_{\phi} + 1) \quad (4.31)$$

Then we can define

$$s\beta_n^{+} : \quad h_+(z) = z^{1-2n}h_{\beta_0^{-}}(z|\Delta_{\phi} - 1 + 2n) \quad (4.32)$$

$$s\alpha_n^{+} : \quad h_+(z) = z^{1-2n}h_{\alpha_0^{-}}(z|\Delta_{\phi} - 1 + 2n) \quad (4.33)$$
for $n \ge 1$. For $n = 0$ this method fails, since we would break the required asymptotic condition near $z = \infty$. For $\alpha_0^+$ however there is a way out: we can form a suitable linear combination of $h_{\beta_0^-}, h_{\alpha_0^-}$ which falls off not as $z^{-3/2}$ as each individual kernel, but as $z^{-5/2}$. We can then set

$$ \alpha_0^+ : \quad h^+(z) = z \left[ h_{\alpha_0^-}(z|\Delta_\phi - 1) + A(\Delta_\phi)h_{\beta_0^-}(z|\Delta_\phi - 1) \right]. \quad (4.34) $$

The choice $A(\Delta_\phi) = H_{\Delta_\phi} - \log(4)$ does the trick, and we can check that the resulting kernel matches the one we obtained in (4.24). This procedure does not work for obtaining $\beta_0^+$, since we are not allowed to use $\alpha_0^-$, so we conclude again that such a functional cannot be constructed.

## 4.4 Bulk functionals

We now turn our attention to the bulk functionals $\theta_n$. The orthonormality conditions (4.2) tell us the corresponding functional actions must have double zeros on all boundary derivative operators. In the bulk channel they must be zero for all $\Delta_m = 2\Delta_\phi + 2m$ larger than $\Delta_n$. Given the functional actions (4.16), (4.13) we conclude that we must demand the behaviours:

$$ \theta_n : \quad \mathrm{Im}[h_+(z)] \stackrel{z \to 0^+}{\sim} O[z^{-1-\frac{1-\nu}{2}+\eta}], \quad \eta > 0 \quad (4.35a) $$

$$ \theta_n : \quad h_+(z) \stackrel{z \to 1^+}{\sim} (z-1)^{-1-n}, \quad (4.35b) $$

where the first condition guarantees that representation (4.16) is valid all the way down to $\hat{\Delta} = \frac{1-\nu}{2} + \Delta_\phi$, and in particular that the functional action is zero for all boundary operators. As may be expected from the discussion of the Neumann case in the previous subsection, we should expect that something should go wrong when $\nu = 1$, and that's indeed what we will find.

After a somewhat circuitous route<fnref n="13" /> we were led to the following simple result for the lowest Dirichlet bulk functional:

$$ \theta_0^- : \quad h_+(z) = h_{\theta_0^-}(z|\Delta_\phi) := N_{\Delta_\phi} \sqrt{z} \left[ \frac{(\Delta_\phi - \frac{1}{2})}{2(z-1)} {}_2\tilde{F}_1 \left( -\frac{1}{2}, \frac{3}{2}, \frac{1}{2} + \Delta_\phi; \frac{1}{z} \right) \right. \\ \left. + {}_2\tilde{F}_1 \left( -\frac{1}{2}, \frac{1}{2}, \Delta_\phi - \frac{3}{2}; \frac{1}{z} \right) - {}_2\tilde{F}_1 \left( -\frac{1}{2}, -\frac{1}{2}, \Delta_\phi - \frac{3}{2}; \frac{1}{z} \right) \right], \quad (4.36) $$

with the normalization factor chosen such that $\theta_0^-(2\Delta_\phi) = 1$,

$$ N_{\Delta_\phi} = \frac{4\Gamma(\Delta_\phi - 1)\Gamma(\Delta_\phi + 1)}{\pi\Gamma(\frac{1}{2} + \Delta_\phi)}. \quad (4.37) $$

The functional kernel is a combination of two pieces which separately satisfy (4.15), namely the top and bottom lines in the above. In particular the first piece is nothing but $\frac{z^2}{z-1}$ times

<footnote n="13">¹³This involved constructing particular solutions for various integer $\Delta_\phi$, obtaining series expansions for large $z$, and finally guessing the general $\Delta_\phi$ result.</footnote>
the boundary functional with kernel $h_{0,0}(z)$ cf. (4.19) (which thanks to (4.17) is again a solution). The specific combination above has the properties

$$h_+(z) \stackrel{z \to \infty}{\sim} O(z^{-\frac{3}{2}}) \qquad (4.38a)$$

$$\mathrm{Im}[h_+(z)] \stackrel{z \to 1^-}{\sim} O[(1-z)^{\Delta_{\phi}-\frac{1}{2}}] \qquad (4.38b)$$

$$h_+(z) \stackrel{z \to 1^+}{\sim} \frac{2}{\pi z - 1} \qquad (4.38c)$$

$$\mathrm{Im}[h_+(z)] \stackrel{z \to 0^+}{\sim} O[\log(z)] \qquad (4.38d)$$

The first two conditions guarantee finiteness and swapping, and the bottom two characterize it as $\theta_0^-$. From this expression it is now easy to obtain all other Dirichlet functionals and (almost) all Neumann ones, by using the shift property (4.17):

$$s\theta_n^-: \quad h_+(z) = \frac{1}{(z-1)^n} h_{\theta_0^-}(z|\Delta_\phi + 2n), \quad n \ge 0 \qquad (4.39)$$

$$s\theta_n^+: \quad h_+(z) = \frac{z}{(z-1)^n} h_{\theta_0^-}(z|\Delta_\phi - 1 + 2n), \quad n \ge 1 \qquad (4.40)$$

We cannot allow $n = 0$ in the last equation since this would destroy the behaviour near $z = \infty$. As discussed in the previous subsection this was expected: for Neumann boundary conditions we cannot build fully orthonormal functionals since this would rule out the contact term interaction. This also give us the clue on how to obtain the $\theta_0^+$: we must combine it with a boundary functional. This will mean that $\theta_0^+$ will no longer have double zeros on all boundary double trace blocks, and upon the Gram-Schmidt decomposition to orthonormalise the $\tilde{\theta}_n^+$ this will be the case for all functionals. A suitable linear combination is

$$\theta_0^+: \quad h_+(z) = z \left[ h_{\theta_0^-}(z|\Delta_\phi - 1) - \frac{\sqrt{\pi \Gamma(\Delta_\phi)}}{\Gamma(\Delta_\phi + \frac{1}{2})} h_{\theta_0^-}(z|\Delta_\phi - 1) \right] \qquad (4.41)$$

where $h_{\beta_0^-}(z|\Delta_\phi)$ is a boundary functional kernel, cf. (4.24a), where we've added an extra label to show the dependence on $\Delta_\phi$. The shift and multiplication by $z$ was necessary to obtain a Neumann type boundary kernel. The linear combination above falls off at infinity as $z^{-3/2}$ unlike each individual term. One can check all other requirements on the functional are satisfied, except that now near $z = 0$ we have $\mathrm{Im}[h_+(z)] \sim O(1/z)$. Hence this functional will have a first order zero at $\hat{\Delta} = \Delta_\phi$.

The conclusion is that, after orthonormalisation, the best we can do for the Neumann basis is:

$$\theta_n(\hat{\Delta}_m) = 0, \qquad \partial_{\hat{\Delta}} \theta_n(\hat{\Delta}_m) = -\delta_{m0} f_n, \qquad \theta_n(\Delta_m) = \delta_{nm} \qquad (4.42)$$

for $m, n \ge 0$. The constants $f_n$ are related to a contact diagram and given in appendix C.3.

## 4.5 Functional actions

Given the functional kernels $h_+(z)$, we are interested in determining the functional actions. We will now summarize how this is done, leaving specific results to appendix A.
For the purposes of computing the $\alpha$ and $\beta$ functional actions, we will use

$$ \omega_{m,k}: \quad h_{m,k}(z) = -\frac{\Gamma(1+\Delta_{\phi}+m)}{\Gamma(1+m-k)} \frac{2\tilde{F}_1\left(k-m-\frac{1}{2}, \frac{3}{2}; m-k, \frac{1}{2}+\Delta_{\phi}+k, \frac{1}{z}\right)}{4^{m-k}\pi^{\frac{3}{2}}z^{\frac{3}{2}+k}} \quad (4.43) $$

For the bulk functionals we will compute the shifted functionals from the kernels

$$ \begin{aligned} s\theta_n^+ &: \quad h_+(z) = \frac{z}{(1-z)^n} h_{\theta_0^-}(z|\Delta_\phi - 1 + 2n) \\ s\theta_n^- &: \quad h_+(z) = \frac{1}{(1-z)^n} h_{\theta_0^-}(z|\Delta_\phi + 2n) \end{aligned} \quad (4.44) $$

with

$$ h_{\theta_0^{-}}(z|\Delta_{\phi}) := N_{\Delta_{\phi}} \sqrt{z} \left[ \frac{(\Delta_{\phi} - \frac{1}{2})}{2(1-z)} {}_2\tilde{F}_1 \left( -\frac{1}{2}, \frac{3}{2}, \frac{1}{2} + \Delta_{\phi}, \frac{1}{z} \right) \right. \\ \left. - {}_2\tilde{F}_1 \left( -\frac{1}{2}, \frac{1}{2}, \Delta_{\phi} - \frac{3}{2}, \frac{1}{z} \right) + {}_2\tilde{F}_1 \left( -\frac{1}{2}, -\frac{1}{2}, \Delta_{\phi} - \frac{3}{2}, \frac{1}{z} \right) \right], \quad (4.45) $$

and the normalization factor:

$$ N_{\Delta_{\phi}} = \frac{4\Gamma(\Delta_{\phi} - 1)\Gamma(\Delta_{\phi} + 1)}{\pi\Gamma(\frac{1}{2} + \Delta_{\phi})}. \quad (4.46) $$

With our choice of kernels we have in particular

$$ s\theta_n^{\pm}(2\Delta_{\phi} + 2n) = 1. \quad (4.47) $$

To evaluate the functional actions we have proceeded as follows. For the $\omega_{m,k}$ functionals we used the representations

$$ \begin{aligned} \omega_{m,k}(\Delta) &= \int_1^\infty dz h_{m,k}(z) \operatorname{Im}[G_\Delta(z|\Delta_\phi)] \\ \omega_{m,k}(\hat{\Delta}) &= \int_1^\infty dz h_{m,k}(z) \operatorname{Im}[\hat{G}_{\hat{\Delta}}(z|\Delta_\phi)] \end{aligned} \quad (4.48) $$

The discontinuities of the blocks are given in terms of Gaussian hypergeometric functions, and so is the kernel $h_{m,k}$. We have evaluated these expressions by replacing the hypergeometric functions by their series representations and commuting the integration with the double series. We check a posteriori equality of the resulting expressions with direct numerical evaluation.

For the shifted bulk functionals we found it necessary to use different representations depending on which action was considered. For the action on bulk blocks we did as for the $\omega_{n,k}$. However, on boundary blocks it was necessary to use instead

$$ s\theta_n(\hat{\Delta}) = -\left(1 - \nu \cos[\pi(\hat{\Delta} - \Delta_\phi)]\right) \int_0^1 dz \operatorname{Im}[h_{s\theta_n}(z)] G_{\hat{\Delta}}(z|\Delta_\phi). \quad (4.49) $$

In both cases the computation of the functional action proceeds as for $\omega_{m,k}$. However, for the bulk functionals the final answer leaves one of the series unperformed.

To conclude note that given the functional actions above, one can go on to obtain the orthonormal functionals by performing a finite Gram-Schmidt orthonormalization, and in some cases we found it possible then to go on to guess the general orthonormal answer.
# 5 Witten diagrams for BCFT

## 5.1 Setup

We will now explicitly compute the Witten diagrams which motivated the functional basis constructed in the previous section. The Witten diagrams for external scalars with interface boundary conditions, were constructed in [33] (see also [35] for the defect case). For Neumann and Dirichlet boundary conditions, we will use similar techniques in our computations. We consider a free scalar field propagating in AdS$_{d+1}$ in the presence of a brane extended along the radial direction. In the Poincaré patch the action is given by

$$S = \int_{x^{\perp}, y \ge 0} \frac{dy}{y^{d+1}} d^{d-1}x \, dx^{\perp} \left[ \frac{1}{2} \nabla \Phi \cdot \nabla \Phi + \frac{1}{2} m^2 \Phi^2 \right] \quad (5.1)$$

with $m^2 = \Delta_{\phi}(\Delta_{\phi} - d)$ and boundary conditions:

$$ \begin{aligned} \Phi(\vec{x}, x^{\perp} = 0, y) &= 0, & \nu &= -1 & (\text{Dirichlet}) \\ \partial_{\perp} \Phi(\vec{x}, x^{\perp} = 0, y) &= 0, & \nu &= 1 & (\text{Neumann}) \end{aligned} \quad (5.2) $$

The bulk field $\Phi$ is dual to the boundary operator $\phi$, which will satisfy the same boundary condition at $x^{\perp} = 0$:

$$ \begin{aligned} \Phi(\vec{x}, x^{\perp}, y) &\stackrel{y \to 0}{\sim} y^{\Delta_{\phi}} \phi(\vec{x}, x^{\perp}) \\ \phi(\vec{x}, x^{\perp} = 0) &= 0, & \nu &= -1 & (\text{Dirichlet}) \\ \partial_{\perp} \phi(\vec{x}, x^{\perp} = 0) &= 0, & \nu &= 1 & (\text{Neumann}) \end{aligned} \quad (5.3) $$

We emphasize that throughout this section, we are not interested in this AdS theory for its own sake, rather it merely serves as a way of generating individual Witten exchange diagrams of interest.

In the absence of the brane, the bulk-to-boundary and bulk-to-bulk propagators would be<fnref n="14" />

$$ \begin{aligned} G_{B\partial}(P, X|\Delta_{\phi}) &= \frac{1}{(-2P \cdot X)^{\Delta_{\phi}}} & (5.4) \\ G_{BB}(X_1, X_2|\Delta_{\phi}) &= \left(-\frac{1}{\zeta}\right)^{\Delta_{\phi}} {}_2F_1\left(\Delta_{\phi}, \Delta_{\phi} - \frac{d-1}{2}, 2\Delta_{\phi} - d + 1, -\frac{4}{\zeta}\right) & (5.5) \end{aligned} $$

with $\zeta = \frac{(X_1 - X_2)^2}{4}$. For economy of space we've introduced the embedding formalism — for a review and more details see [36–38]. In practice, we just need to know that $P^M, X^M$ are $d+2$-dimensional vectors satisfying $X^2 = -1, P^2 = 0$, with

$$ X^M = \{X^+, X^-, X^i, X^{\perp}\} = \frac{1}{y} \{1, |\vec{x}|^2 + (x^{\perp})^2 + y^2, x^i, x^{\perp}\} \quad (5.6) $$

$$ P^M = \{P^+, P^-, P^i, P^{\perp}\} = \{1, |\vec{x}|^2 + (x^{\perp})^2, x^i, x^{\perp}\} \quad (5.7) $$

<footnote n="14">¹⁴Up to normalization factors which will be irrelevant for us.</footnote>
with $i = 1, \dots, d-1$. This gives for instance

$$-2P_1 \cdot P_2 = |\vec{x}_{12}|^2 + (x_1^{\perp} - x_2^{\perp})^2 \qquad (5.8)$$

$$-2P_1 \cdot X_2 = \frac{1}{y_2} \left[ |\vec{x}_{12}|^2 + (x_1^{\perp} - x_2^{\perp})^2 + y_2^2 \right] \qquad (5.9)$$

$$\zeta = \frac{|\vec{x}_{12}|^2 + (x_1^{\perp} - x_2^{\perp})^2 + (y_1 - y_2)^2}{4y_1 y_2}. \qquad (5.10)$$

In the presence of the brane, the propagators get modified. Since this is a free theory, it is straightforward to work out the correct propagators via the method of images. The result is that:

$$G_{BB}^{\nu}(X_1, X_2 | \Delta_{\phi}) = G_{BB}(X_1, X_2 | \Delta_{\phi}) + \nu G_{BB}(X_{1r}, X_2 | \Delta_{\phi}), \qquad (5.11a)$$

$$G_{B\partial}^{\nu}(P, X | \Delta_{\phi}) = G_{B\partial}(P, X | \Delta_{\phi}) + \nu G_{B\partial}(P_r, X | \Delta_{\phi}) \\ = G_{B\partial}(P, X | \Delta_{\phi}) + \nu G_{B\partial}(P, X_r | \Delta_{\phi}), \qquad (5.11b)$$

where a subscript 'r' indicates a reflection $x^{\perp} \to -x^{\perp}$.

Below we will consider three kinds of Witten diagram. These correspond to exchanges of fields propagating in the AdS$_{d+1}$ bulk (bulk exchange), fields which propagate on the AdS$_d$ brane (boundary exchange) and finally we will also consider the simplest possible contact interaction with the brane. These diagrams will compute corrections to a BCFT$_d$ two point function $\langle \phi_1(\vec{x}_1, x_1^{\perp}) \phi_2(\vec{x}_2, x_2^{\perp}) \rangle$ of fields with dimensions $\Delta_1, \Delta_2$, although we will eventually be interested in setting $\phi_1 = \phi_2 = \phi$ and $\Delta_1 = \Delta_2 = \Delta_{\phi}$. We will write for the two point functions:

$$W(P_1, P_2) = \frac{W(z)}{(x_1^{\perp})^{\Delta_1} (x_2^{\perp})^{\Delta_2}}, \qquad (5.12)$$

where $z$ is the cross-ratio we've been using so far, and which can be written in terms of $P_1$ and $P_2$ explicitly as:

$$z = - \frac{2P_1^{\perp} P_2^{\perp}}{P_1 \cdot P_2}. \qquad (5.13)$$

We will stick to the general rule that hatted quantities relate to defect and unhatted to bulk. Hence Witten diagrams corresponding to boundary (bulk) exchanges will be represented by hatted (unhatted) quantities.

## 5.2 Boundary exchange

The boundary exchange Witten diagram is represented in figure 3. The diagram is then computed for the two different boundary conditions, as follows:

Let us put Neumann boundary conditions on the AdS$_{d+1}$ scalars $\Phi_1$ and $\Phi_2$, which are the external legs of the diagram. We assume the vertices $\Phi_i \hat{\Psi}_{\Delta}$, with $i = 1, 2$, on the AdS$_d$ defect. Here $\hat{\Psi}_{\Delta}$ is a freely propagating scalar field on AdS$_d$, of dimension $\hat{\Delta}$. The boundary-bulk propagators end on the AdS$_d$ brane, where they become the same as ordinary AdS$_{d+1}$ bulk to boundary propagators (up to a factor of 2). The diagram above then reduces to the same boundary exchange Witten diagram already evaluated in [33],

$$\widehat{W}^{+}(P_1, P_2) = \int_{AdS_d} dW_1 dW_2 G_{B\partial}^{+}(P_1, W_1 | \Delta_1) G_{B\partial}^{+}(P_2, W_2 | \Delta_2) G_{BB}(W_1, W_2 | \hat{\Delta}). \qquad (5.14)$$
<figure><img src="image_3.png" /><figcaption>Figure 3. The boundary exchange Witten diagram.</figcaption></figure>

The integration runs over the AdS$_d$ brane, and the measure is simply the ordinary one on AdS$_d$. The $G_{BB}(W_1, W_2|\hat{\Delta})$ is an ordinary bulk-to-bulk propagator in AdS$_d$.

For Dirichlet boundary conditions, the vertex is instead $\partial_{\perp} \Phi_i \hat{\Psi}_{\hat{\Delta}}$ on the AdS$_d$ brane. The diagram has the following form,

$$ \widehat{W}^{-}(P_1, P_2) = \int_{AdS_d} dW_1 dW_2 \partial_{\perp}^{W_1} G_{B\partial}^{-}(P_1, W_1|\Delta_1) \partial_{\perp}^{W_2} G_{B\partial}^{-}(P_2, W_2|\Delta_2) G_{BB}(W_1, W_2|\hat{\Delta}). \qquad (5.15) $$

Here we denote $\partial_{\perp}^{W} G(P, W) = \lim_{X \to W} \partial_{x^{\perp}} G(P, X)$. Near the brane we have:

$$ G_{B\partial}^{-}(P, X) \xrightarrow{X \to W + O(x^{\perp})} \frac{x^{\perp}}{(-2P \cdot W)} G_{B\partial}^{+}(P, W). \qquad (5.16) $$

Thanks to the derivatives in the integral (5.15), the integral then becomes similar to the Neumann case, with boundary-bulk propagators whose external scalar dimensions have effectively increased by one.

For unequal external scalar fields $\Phi_1, \Phi_2$, and Neumann boundary conditions the result is given as follows:

$$ \widehat{W}_{\hat{\Delta}}^{+}(z) = \mathcal{N}_{\text{bdy}}^{N}(\hat{\Delta}) \int_{-i\infty}^{i\infty} d\tau \left[ \frac{z}{4(1-z)} \right]^{\tau} \frac{\Gamma(\tau)\Gamma(h' - \tau)}{\Gamma(2h' - 2\tau)} \int_{-i\infty}^{i\infty} dc \frac{\hat{f}(c, \tau)\hat{f}(-c, \tau)}{(\hat{\Delta} - h')^2 - c^2}, \qquad (5.17) $$

where $h' = (d-1)/2$ and,

$$ \hat{f}(c, \tau) = \frac{\Gamma\left(\frac{c+\Delta_1-h'}{2}\right) \Gamma\left(\frac{c+\Delta_2-h'}{2}\right) \Gamma(c+h' - \tau)}{2\Gamma(c)}. \qquad (5.18) $$

The normalization factor is chosen such that when decomposing into boundary blocks we find:

$$ z^{-\frac{\Delta_1+\Delta_2}{2}} \widehat{W}_{\hat{\Delta}}^{+}(z) = \widehat{G}_{\hat{\Delta}}(z|\Delta_1, \Delta_2) - \sum_{n=0}^{\infty} \left[ \hat{a}_n^{+}(\Delta_1, \Delta_2, \hat{\Delta}) \widehat{G}_{\Delta_1+2n}(z|\Delta_1, \Delta_2) + \hat{a}_n^{+}(\Delta_2, \Delta_1, \hat{\Delta}) \widehat{G}_{\Delta_2+2n}(z|\Delta_1, \Delta_2) \right]. \qquad (5.19) $$
Setting $\Delta_1 = \Delta_\phi + \frac{\Delta_{12}}{2}$, $\Delta_2 = \Delta_\phi - \frac{\Delta_{12}}{2}$ and sending $\Delta_{12} \to 0$ we find

$$ \hat{a}_n^+(\Delta_1, \Delta_2, \hat{\Delta}) \stackrel{\Delta_{12} \to 0}{\sim} \frac{\hat{b}_n^+(\hat{\Delta})}{\Delta_{12}} + \frac{1}{2} \hat{a}_n^+(\hat{\Delta}) + O(\Delta_{12}). \quad (5.20) $$

We have omitted the dependence of $\hat{b}_n^+, \hat{a}_n^+$ on $\Delta_1 = \Delta_2 = \Delta_\phi$ for simplicity. The boundary block expansion becomes:

$$ z^{-\Delta_{\phi}} \hat{W}_{\hat{\Delta}}^{+}(z) = \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi}) - \sum_{n=0}^{\infty} \left[ \hat{a}_{n}^{+}(\hat{\Delta}) \hat{G}_{\Delta_{\phi}+2n}(z|\Delta_{\phi}) + \hat{b}_{n}^{+}(\hat{\Delta}) \partial_{\hat{\Delta}} \hat{G}_{\Delta_{\phi}+2n}(z|\Delta_{\phi}) \right]. \quad (5.21) $$

As for the bulk channel we get instead

$$ \begin{aligned} \hat{W}_{\hat{\Delta}}^{+}(z) &= z^{\frac{\Delta_1 + \Delta_2}{2}} \sum_{n} \hat{t}_{n}^{N}(\Delta_1, \Delta_2, \hat{\Delta}) G_{\Delta_1 + \Delta_2 + 2n}(z|\Delta_1, \Delta_2). \\ \stackrel{\Delta_{12}=0}{=} z^{\Delta_{\phi}} \sum_{n} \hat{t}_{n}^{+}(\hat{\Delta}) G_{2\Delta_{\phi} + 2n}(z|\Delta_{\phi}). \end{aligned} \quad (5.22) $$

The Dirichlet Witten diagram $\hat{W}_{\hat{\Delta}}^{-}$ is obtained from $\hat{W}_{\hat{\Delta}}^{+}$, by shifting the dimensions $\Delta_1, \Delta_2 \to \Delta_1 + 1, \Delta_2 + 1$. The corresponding boundary block coefficients can be directly obtained this way. The bulk block coefficients can be obtained by inspecting the expansion around $z \to 1$, where the expansion coefficients can be obtained from the above shifts. Explicit expressions for all coefficients appearing here are relegated to appendix C.1.

## 5.3 Bulk exchange

The bulk exchange Witten diagram is shown in figure 4. For this diagram, we assume a vertex $\Phi_1\Phi_2\Psi_\Delta$ in the bulk, and a linear vertex $\Psi_\Delta$ on the brane. The $\Psi_\Delta$ is a freely propagating scalar operator in $AdS_{d+1}$ that gets exchanged. We put Neumann/Dirichlet boundary conditions on $\Phi_1$ and $\Phi_2$, but on $\Psi_\Delta$ we always put Neumann boundary conditions.<fnref n="15" /> The diagram is then given by:

$$ W_{\Delta}^{\nu}(P_1, P_2) = \int_{\frac{AdS_{d+1}}{2}} dX \int_{AdS_d} dW \, G_{B\partial}^{\nu}(P_1, X|\Delta_1) G_{B\partial}^{\nu}(P_2, X|\Delta_2) G_{BB}(X, W|\Delta). \quad (5.23) $$

The integral label $\frac{AdS_{d+1}}{2}$ indicates that the integration is restricted to $x^{\perp} \ge 0$. Note that the bulk-to-bulk propagator is the ordinary one. This is easy to see since irrespective of the boundary conditions on $\Phi_1$ and $\Phi_2$, the exchange operator always has Neumann boundary condition. So the propagator $G_{BB}^{\nu}(X_1, X_2|\Delta)$ from (5.11a) reduces to (twice of) $G_{BB}(X_1, X_2|\Delta)$ when one of the points $X_1$ and $X_2$ lies on the $AdS_d$ brane.

Using the propagators (5.11) and rearranging we can get:

$$ W_{\Delta}^{\nu}(P_1, P_2) = \int_{AdS_{d+1}} dX \int_{AdS_d} dW \left[ G_{B\partial}(P_1, X|\Delta_1) G_{B\partial}(P_2, X|\Delta_2) G_{BB}(X, W|\Delta) \right. \\ \left. \qquad + \nu G_{B\partial}(P_{1r}, X|\Delta_1) G_{B\partial}(P_2, X|\Delta_2) G_{BB}(X, W|\Delta) \right]. \quad (5.24) $$

<footnote n="15">¹⁵We could have chosen a different boundary condition at the cost of modifying the vertex on the brane.</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 4.** The bulk exchange Witten diagrams can be written as a sum over two diagrams on full AdS$_{d+1}$. In the figure we show the Neumann case, Dirichlet would correspond to a minus sign on the r.h.s.</figcaption></figure>

These integrals can be computed using standard tools summarized nicely in [33]. The first integral is straightforward, while the second is more involved. We give the details of this computation in appendix B. The final expression is given by:

$$ W_{\Delta}^{\nu}(z) = \mathcal{N}_{\text{bulk}}^{\nu}(\Delta) \int_{-i\infty}^{+i\infty} \frac{d\tau}{2\pi i} \left[ \frac{4(1-z)}{z} \right]^{\tau - \frac{\Delta_1 + \Delta_2}{2}} \frac{\Gamma(\tau)\Gamma(\tau - \frac{\Delta_{12}}{2})\Gamma(\tau - \frac{\Delta_{21}}{2})}{\Gamma(2\tau)\Gamma(\frac{1}{2} - \tau)} \\ \times \int_{-i\infty}^{+i\infty} \frac{dc}{2\pi i} \frac{f(c, \tau)f(-c, \tau)}{((\Delta - h)^2 - c^2)} \left[ 1 + \nu \frac{2F_1(\frac{1-h+c}{2}, \frac{1-h-c}{2}, \frac{1}{2} - \tau, z)}{(1-z)^{\frac{d-\Delta_1-\Delta_2-1}{2}}} \right], \quad (5.25) $$

where $h = d/2$ and,

$$ f(c, \tau) = \frac{\Gamma(\frac{\Delta_1+\Delta_2-h+c}{2})\Gamma(\frac{1+c-h}{2})\Gamma(\frac{h+c}{2}-\tau)}{2\Gamma(c)}. \quad (5.26) $$

The above expression<fnref n="16" /> has the following boundary channel decomposition (as $z \to 0$),

$$ W_{\Delta}^{\nu}(z) = z^{\frac{\Delta_1 + \Delta_2}{2}} \sum_{n=0}^{\infty} \quad (5.27) \\ \times \left[ a_{n}^{\nu}(\Delta_1, \Delta_2, \Delta) \hat{G}_{\Delta_1+2n+\frac{1-\nu}{2}}(z|\Delta_1, \Delta_2) + a_{n}^{\nu}(\Delta_2, \Delta_1, \Delta) \hat{G}_{\Delta_2+2n+\frac{1-\nu}{2}}(z|\Delta_1, \Delta_2) \right]. $$

The fact that we only get e.g. $\Delta_1$ plus an even integer depends crucially on the combination of terms in (5.25). As before, when we approach $\Delta_1 \to \Delta_\phi + \frac{\Delta_{12}}{2}$ and $\Delta_2 \to \Delta_\phi - \frac{\Delta_{12}}{2}$ we find

$$ a_n^\nu(\Delta_1, \Delta_2, \Delta) \stackrel{\Delta_{12} \to 0}{\rightleftharpoons} \frac{b_n^\nu(\Delta)}{\Delta_{12}} + \frac{1}{2} a_n^\nu(\Delta) + O(\Delta_{12}), \quad (5.28) $$

and the boundary block expansion becomes:

$$ W_{\Delta}^{\nu}(z) = z^{\Delta_{\phi}} \sum_{n=0}^{\infty} \left[ a_{n}^{\nu}(\Delta) \hat{G}_{\Delta_{\phi}+2n+\frac{1-\nu}{2}}(z|\Delta_{\phi}) + b_{n}^{\nu}(\Delta) \partial_{\Delta} \hat{G}_{\Delta_{\phi}+2n+\frac{1-\nu}{2}}(z|\Delta_{\phi}) \right]. \quad (5.29) $$

<footnote n="16">¹⁶Note the presence of the extra hypergeometric $_2F_1$ function, which was not there for the interface case considered in [33]. It is needed to implement the boundary condition.</footnote>
<figure><img src="image_5.png" /><figcaption>Figure 5. The contact Witten diagram.</figcaption></figure>

Finally, in the bulk channel we find

$$ W_{\Delta}^{\nu}(z) = z^{\frac{\Delta_1+\Delta_2}{2}} \left[ G_{\Delta}^{\Delta_1 2}(z|\Delta_1, \Delta_2) - \sum_n t_n^{\nu}(\Delta_1, \Delta_2, \Delta) G_{\Delta_1+\Delta_2+2n}^{\Delta_1 2}(z|\Delta_1, \Delta_2) \right] \\ \stackrel{\Delta_1 2=0}{=} z^{\Delta_{\phi}} \left[ G_{\Delta}(z|\Delta_{\phi}) - \sum_n t_n^{\nu}(\Delta) G_{2\Delta_{\phi}+2n}(z|\Delta_{\phi}) \right]. \qquad (5.30) $$

All expansion coefficients are given later in appendix C.2.

## 5.4 Contact diagrams

The final tree level Witten diagrams we will consider are contact diagrams. They can be represented by two $AdS_{d+1}$ boundary-bulk propagators meeting at the $AdS_d$ boundary where they interact via a quadratic vertex. There are several kinds of such diagrams depending on the choice of this vertex. Here we will consider the simplest possible cases.

With Neumann boundary condition, the simplest such vertex is $\Phi_1\Phi_2$ vertex on the $AdS_d$ brane. We have:

$$ W_C^+ = \int_{AdS_d} dW G_{B\partial}^+(P_1, W|\Delta_1) G_{B\partial}^+(P_2, W|\Delta_2). \qquad (5.31) $$

For Dirichlet, we consider the vertex $\partial_{\perp}\Phi_1\partial_{\perp}\Phi_2$. There is a simplification here similar to the Dirichlet boundary exchange diagram. Following the notations of subsection 5.2, the contact diagram can be written as follows,

$$ W_C^- = \int_{AdS_d} dW \partial_{\perp}^W G_{B\partial}^-(P_1, W|\Delta_1) \partial_{\perp}^W G_{B\partial}^-(P_2, W|\Delta_2). \qquad (5.32) $$

Due to (5.16) and the derivatives on the propagators, the above integral becomes equivalent to the Neumann case (5.31) with effective external dimensions shifted as $\Delta_1 \to \Delta_1 + 1$ and $\Delta_2 \to \Delta_2 + 1$.

For Neumann boundary conditions, the boundary-bulk propagators are just the ordinary $AdS_{d+1}$ propagators (times 2). So the above integral simplifies to the interface case diagram computed in [33]. For both boundary conditions, the contact diagram can be
written in a compact way as follows,

$$W_C^\nu(z) = \int_{-i\infty}^{i\infty} d\tau \left(\frac{1-z}{z}\right)^{-\tau} \frac{\Gamma(\Delta_1 - \tau + \frac{1-\nu}{2})\Gamma(\Delta_2 - \tau + \frac{1-\nu}{2})\Gamma(\tau)}{\Gamma(\frac{\Delta_1+\Delta_2-\nu}{2} - \tau + 1)}. \quad (5.33)$$

As $z \to 0$ this has the decomposition,

$$W_C^\nu(z) = z^{\frac{\Delta_1+\Delta_2}{2}} \sum_{n=0}^{\infty} [d_n^\nu(\Delta_1, \Delta_2) \hat{G}_{\Delta_1+2n+\frac{1-\nu}{2}}(z|\Delta_1, \Delta_2) + d_n^\nu(\Delta_2, \Delta_1) \hat{G}_{\Delta_2+2n+\frac{1-\nu}{2}}(z|\Delta_1, \Delta_2)] \\ \Delta_1 = \Delta_2 = \Delta_\phi \qquad z^{\Delta_\phi} \sum_{n=0}^{\infty} \left[ d_n^\nu \hat{G}_{\Delta_\phi+2n+\frac{1-\nu}{2}}(z|\Delta_\phi) + e_n^\nu \partial_{\hat{\Delta}} \hat{G}_{\Delta_\phi+2n+\frac{1-\nu}{2}}(z|\Delta_\phi) \right], \quad (5.34)$$

where as before, when we approach $\Delta_1 \to \Delta_\phi + \frac{\Delta_{12}}{2}$ and $\Delta_2 \to \Delta_\phi - \frac{\Delta_{12}}{2}$, we have

$$d_n^\nu(\Delta_1, \Delta_2) \stackrel{\Delta_{12} \to 0}{\rightleftharpoons} \frac{e_n^\nu}{\Delta_{12}} + \frac{1}{2} d_n^\nu + O(\Delta_{12}). \quad (5.35)$$

As $z \to 1$ we have the bulk channel decomposition,

$$W_C^\nu(z) = -z^{\frac{\Delta_1+\Delta_2}{2}} \sum_n f_n^\nu(\Delta_1, \Delta_2) G_{\Delta_1+\Delta_2+2n}^{\Delta_{12}}(z|\Delta_1, \Delta_2) \\ \stackrel{\Delta_{12}=0}{\rightleftharpoons} -z^{\Delta_\phi} \sum_n f_n^\nu G_{2\Delta_\phi+2n}^{\Delta_\phi}(z|\Delta_\phi). \quad (5.36)$$

The expansion coefficients for the contact diagram are given in appendix C.3.

# 6 Polyakov and functional bootstrap

## 6.1 Polyakov blocks

In section 4 we constructed two interesting bases of linear functionals that act on the BCFT crossing equation. This construction was partly motivated by considering the general structure of Witten diagrams in a bulk dual to the generalized free field BCFT. In the previous section we have explicitly computed these diagrams. We now close the circle and show in detail how these are related.

The main claim of this section is that there exist Polyakov bulk and boundary blocks $P_\Delta(z|\Delta_\phi)$ and $\hat{P}_{\hat{\Delta}}(z|\Delta_\phi)$, for both Neumann and Dirichlet boundary conditions, which satisfy:

$$\begin{aligned} P_\Delta(z|\Delta_\phi) &= G_\Delta(z|\Delta_\phi) - \sum_{n=0}^{\infty} \theta_n(\Delta|\Delta_\phi) G_{\Delta_n}(z|\Delta_\phi) \\ &= \sum_{n=0}^{\infty} \left[ \alpha_n(\Delta|\Delta_\phi) \hat{G}_{\hat{\Delta}_n}(z|\Delta_\phi) + \beta_n(\Delta|\Delta_\phi) \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n}(z|\Delta_\phi) \right] \end{aligned} \quad (6.1a)$$

$$\begin{aligned} \hat{P}_{\hat{\Delta}}(z|\Delta_{\phi}) &= \hat{G}_{\hat{\Delta}}(z|\Delta_{\phi}) - \sum_{n=0}^{\infty} \left[ \alpha_n(\hat{\Delta}|\Delta_{\phi}) \hat{G}_{\hat{\Delta}_n}(z|\Delta_{\phi}) + \beta_n(\hat{\Delta}|\Delta_{\phi}) \partial_{\hat{\Delta}} \hat{G}_{\hat{\Delta}_n}(z|\Delta_{\phi}) \right] \\ &= \sum_{n=0}^{\infty} \theta_n(\hat{\Delta}) G_{\Delta_n}(z). \end{aligned} \quad (6.1b)$$
Concretely, we will show that the Polyakov blocks may be obtained from the Witten dia-
grams computed in the previous section. Indeed, note that the crossing relations satisfied
by the Polyakov blocks look remarkably similar to those arising from Witten exchange
diagrams, cf. (4.1a), (4.1b), but with the notable difference that here the coefficients in
the expansion are functional actions. We will show that it is possible to combine exchange
diagrams with contact terms in a way that in those equations we effectively set:

$$
a_n(\Delta) \rightarrow \alpha_n(\Delta), \quad \hat{a}_n(\hat{\Delta}) \rightarrow \alpha_n(\hat{\Delta}) \qquad (6.2a)
$$

$$
b_n(\Delta) \rightarrow \beta_n(\Delta), \quad \hat{b}_n(\hat{\Delta}) \rightarrow \beta_n(\hat{\Delta}) \qquad (6.2b)
$$

$$
c_n(\Delta) \rightarrow \theta_n(\Delta), \quad \hat{c}_n(\hat{\Delta}) \rightarrow \theta_n(\hat{\Delta}) \qquad (6.2c)
$$

In other words, the functional bases allow us to bootstrap Witten exchange diagrams.

Before we do this, let us briefly comment on why these Polyakov blocks are interesting
objects. The role in life of the Polyakov blocks is that they provide us with a way of
writing a generic two-point function in an explicitly crossing symmetric way. Concretely,
for a given two point function $\mathcal{G}(z)$, we have:

$$
\begin{align}
\mathcal{G}(z) &= \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \hat{G}_{\hat{\Delta}}(z) |\Delta_{\phi}) = \sum_{\Delta} (a\lambda)_{\Delta} G_{\Delta}(z|\Delta_{\phi}) \nonumber \\
&= \sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \hat{P}_{\hat{\Delta}}(z|\Delta_{\phi}) + \sum_{\Delta} (a\lambda)_{\Delta} P_{\Delta}(z) \tag{6.3}
\end{align}
$$

This remarkable statement is true only if the contributions from the unphysical bulk double
trace operators and boundary derivative operators appearing in the conformal block decom-
positions of the Polyakov blocks drop out from the total sum. If we plug in the expression
for the Polyakov blocks into the above equation and *commute the two series*, one finds that
the unphysical states drop out precisely if the functional bootstrap equations (4.5a) hold,
which we repeat here for convenience:

$$
\sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \alpha_n(\hat{\Delta}) = \sum_{\Delta} (a\lambda)_{\Delta} \alpha_n(\Delta) \qquad (6.4a)
$$

$$
\sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \beta_n(\hat{\Delta}) = \sum_{\Delta} (a\lambda)_{\Delta} \beta_n(\Delta) \qquad (6.4b)
$$

$$
\sum_{\hat{\Delta}} \mu_{\hat{\Delta}}^2 \theta_n(\hat{\Delta}) = \sum_{\Delta} (a\lambda)_{\Delta} \theta_n(\Delta) \qquad (6.4c)
$$

This establishes the link between the so-called Polyakov bootstrap and our functional
bootstrap approach.

Of course, we have not shown that it is actually possible to commute the series. This
would amount to showing that the functional bases we have constructed is in fact complete
in some sense, i.e. that the functional bootstrap equations are completely equivalent to the
original crossing equation. We will not attempt to prove this here. However, given the fact
that these equations allow us to bootstrap arbitrary deformations of the generalized free
solution, it seems likely that this is correct.<fnref n="17" />

<footnote n="17">¹⁷A similar problem arose in [19]. There the proof of completeness crucially relied on establishing upper bounds on the OPE density, so we expect the same must be true here. However, establishing such bounds seems more difficult here due to lack of positivity of the bulk channel expansion.</footnote>
## 6.2 Dirichlet case

Let us then show how Polyakov blocks and Witten diagrams are related, beginning with the Dirichlet case. In this case, we expect the situation to be simpler than for Neumann, since there were no ambiguities in the definition of our basis of functionals. As we’ve discussed before, there is a simple reason for this.

Ambiguities in the basis correspond to the existence of solutions to crossing which preserve the form of the conformal block decomposition and which decay sufficiently fast at large $z$. This is crucial, since we have constructed our functional basis by assuming the functionals should satisfy the swapping property, but this is dependent on the large $z$ behaviour of a correlator. It can be checked by direct computation that the contact interaction with smallest number of derivatives and Dirichlet boundary conditions behaves like $\sim \sqrt{z}$ for large $z$. Since this is worse behaviour than for a generic two point function, which is only constant in the same limit, this explains the absence of ambiguities in this case, since such contact terms cannot be bootstrapped with our basis.<fnref n="18" />

On the other hand, in our computation of the Witten diagrams we made specific choices of vertices on the brane. For Dirichlet boundary conditions it is straightforward to check that the boundary exchange Witten diagram falls off as $1/\sqrt{z}$ (as can be determined e.g. from its conformal block decomposition). The bulk exchange diagram also has the same fall-off behaviour. This can be checked by taking special values of the exchange dimension e.g. $\Delta = 2\Delta_\phi - 2$.<fnref n="19" /> Then we can infer that 1) we cannot shift the Witten exchange correlators by a contact interaction without changing the large $z$ behaviour, and 2) that there should be no issue in acting directly with our functional basis on these correlators. By using the orthonormality properties of the functionals, it is natural to conclude:

$$z^{-\Delta_{\phi}} \mathcal{W}_{\Delta}^{-}(z) = \mathcal{P}_{\Delta}^{-}(z) \qquad (6.5)$$

$$z^{-\Delta_{\phi}} \widehat{\mathcal{W}}_{\Delta}^{-}(z) = \widehat{\mathcal{P}}_{\Delta}^{-}(z). \qquad (6.6)$$

That is, the Dirichlet functionals directly compute the conformal block decompositions of Dirichlet Witten exchange diagrams. To check this we should test relations (6.2). We have evaluated the functional actions for several $n$ using (4.16), (4.13) and/or the original definition (3.6), as well as the kernels constructed in section 4. To the extent we were able to determine these, the results are presented in appendix A. We then compared them with the conformal block decompositions of the Witten diagrams given in appendix C. In some cases it is possible to do this analytically (in particular for the $\alpha, \beta$ funcitonals), in others it is more convenient to do numerical checks. In all cases we found perfect agreement.

## 6.3 Neumann case

The Neumann case is more interesting. Here the simplest contact interaction falls off as $z^{-1/2}$ near infinity, while the boundary exchange Witten diagram computed in the section 5 falls off like $z^{-3/2}$.<fnref n="20" /> This is interesting for two reasons. Firstly, this means that our basis

<footnote n="18">¹⁸They can however be bootstrapped by demanding stronger falloff conditions on the functional basis.</footnote>

<footnote n="19">¹⁹In general for $\Delta = 2\Delta_\phi - 2m$ with integer $m$ the expression for $W_\Delta^\nu$ simplifies (see [33]).</footnote>

<footnote n="20">²⁰Recall the boundary channel decomposition in the Neumann case is related to the Dirichlet one by shifting everything by $\Delta_\phi \to \Delta_\phi - 1$ and dividing by $z$.</footnote>
of functionals can act on any linear combination of the exchange diagram and the contact interaction, since it is suitable for correlators which grow as $z^{1/2-\eta}$ for any $\eta > 0$. Secondly, recall that there is a missing Neumann functional precisely because we demanded that the falloff at infinity was sufficiently fast. This prefunctional was such that the decay at infinity was only $1/\sqrt{z}$. But we see that this is sufficient to act on the boundary exchange Witten diagram, since overall we get a $1/z^2$ decay. The missing functional is simply

$$ \tilde{\beta}_0^+ : \qquad h_+(z) = z h_{\beta_0^-}(z|\Delta_\phi - 1) \qquad (6.7) $$

We put a tilde to emphasize these functionals have poorer asymptotic behaviour. We can combine it with the remaining basis functionals constructed in section 4 to get a complete set of fully orthonormal prefunctionals, all of which decay at infinity as $z^{-1/2}$. From the orthonormality relations (4.26) and (4.42) the correct combination is

$$ \begin{aligned} \tilde{\alpha}_n^+ &= \alpha_n^+ + d_n \tilde{\beta}_0^+ \\ \tilde{\beta}_n^+ &= \beta_n^+ + e_n \tilde{\beta}_0^+, \\ \tilde{\theta}_n^+ &= \theta_n^+ + f_n \tilde{\beta}_0^+. \end{aligned} \qquad (6.8a) $$

The claim then is that these prefunctionals directly compute the conformal block decomposition of the Witten diagram, for example

$$ z^{-\Delta_{\phi}} \tilde{W}_{\Delta}^{+}(z) = \hat{G}_{\Delta}(z) - \sum_{n=0}^{\infty} \left[ \tilde{\alpha}_{n}^{+}(\hat{\Delta}) \hat{G}_{\Delta_{\phi}+2n}(z) + \tilde{\beta}_{n}^{+}(\hat{\Delta}) \partial_{\Delta} \hat{G}_{\Delta_{\phi}+2n}(z) \right]. \qquad (6.9) $$

The $\tilde{\alpha}^+, \tilde{\beta}^+$ functional actions are determined in appendix A, and they can be indeed matched with the conformal block decompositions of the Neumann boundary exchange diagram. It works analogously for the bulk channel decomposition.

The bulk exchange diagram (5.25) however falls off as $1/\sqrt{z}$, so the above prefunctionals cannot be used directly. One can however modify the diagram by adding a contact diagram to (5.25), since the resulting expression will still have the correct block decomposition as in (5.28) and (5.29). Since the contact diagram also has the $1/\sqrt{z}$ fall-off, adding it with the right factor can cancel the leading large $z$ behavior, so that the modified diagram falls off as $z^{-3/2}$. We can write the modified bulk exchange diagram as,

$$ \tilde{W}_{\Delta}^{+}(z) = W_{\Delta}^{+}(z) - \frac{w_1}{w_2} W_{C}^{+}(z), \qquad (6.10) $$

where the functions $w_1$ and $w_2$ are defined by the respective fall-off behaviours $z^{-\Delta_\phi} W_\Delta^{+} \xrightarrow{z \gg 1} w_1 z^{-1/2}$ and $z^{-\Delta_\phi} W_C^{+} \xrightarrow{z \gg 1} w_2 z^{-1/2}$. The above object now allows the action of the prefunctionals (6.8), and hence has a block decomposition,

$$ \begin{aligned} z^{-\Delta_{\phi}} \tilde{W}_{\Delta}^{+}(z) &= \sum_{n=0}^{\infty} \left[ \tilde{\alpha}_{n}^{+}(\Delta) \hat{G}_{\Delta_{\phi}+2n}(z) + \tilde{\beta}_{n}^{+}(\Delta) \partial_{\Delta} \hat{G}_{\Delta_{\phi}+2n}(z) \right] \\ &= G_{\Delta}(z) - \sum_{n=0}^{\infty} \tilde{\theta}_{n}^{+}(\Delta) G_{2\Delta_{\phi}+2n}(z). \end{aligned} \qquad (6.11) $$
We have not obtained the Polyakov block yet, since we want the coefficients to be the good functional basis with stronger falloff near $z = \infty$. Hence we expect that we need to take a suitable linear combination of the Witten exchange and the contact diagram to recover it. This combination is easily found. Since our basis does not contain $\beta_0^+$, we should have

$$z^{\Delta_{\phi}} \hat{P}_{\hat{\Delta}}^{+}(z) = \hat{W}_{\hat{\Delta}}^{+}(z) - \tilde{\beta}_{0}(\hat{\Delta}) W_{C}^{+}(z) \qquad (6.12)$$

$$z^{\Delta_{\phi}} \mathcal{P}_{\Delta}^{+}(z) = \tilde{W}_{\Delta}^{+}(z) - \tilde{\beta}_{0}(\Delta) W_{C}^{+}(z). \qquad (6.13)$$

It can be checked that for both the above objects the coefficient of $\partial_{\hat{\Delta}} \hat{G}_{\Delta_{\phi}}(z)$ cancels in the boundary conformal block expansions. This is compatible with (6.8) if and only if $d_n, e_n, f_n$ compute the contact diagram, which we have already shown in the previous section. Note that we can alternatively define $\mathcal{P}_{\Delta}^{+}$ from $W_{\Delta}^{+}$ by adding $W_{C}^{+}$ and requiring the coefficient of $\partial_{\hat{\Delta}} \hat{G}_{\Delta_{\phi}}(z)$ to vanish from the boundary channel expansion (which is natural to expect given the decomposition of $\hat{P}_{\Delta}^{+}$). This way one arrives at the same expression as above. It is straightforward to check that the condition (6.2) holds.

To conclude, as pointed out in section 4 our choice of basis of functionals is ambiguous in the Neumann case, since we cannot have fully orthonormal functionals. This ambiguity is reflected in the fact that unlike for Dirichlet boundary conditions, there is no longer a canonical Polyakov block because of the existence of the contact interaction. In particular, if $\omega$ is a good functional satisfying the correct fall-off condition, we are always free to shift:

$$\hat{P}_{\hat{\Delta}}^{+} \rightarrow \hat{P}_{\hat{\Delta}}^{+} + \omega(\hat{\Delta}) z^{-\Delta_{\phi}} W_{C}(z). \qquad (6.14)$$

$$\mathcal{P}_{\Delta}^{+} \rightarrow \mathcal{P}_{\Delta}^{+} + \omega(\Delta) z^{-\Delta_{\phi}} W_{C}(z). \qquad (6.15)$$

# 7 Application — Wilson-Fisher at $O(\epsilon^2)$

The functional bootstrap equations (4.5a) reformulate the crossing equation in a form which is ideally suited for perturbative expansions. As a particular application, we will reconstruct the BCFT data of the Wilson-Fisher fixed point, i.e. $\lambda\phi^4$ theory in $d = 4 - \epsilon$ dimensions tuned to the critical point) on a half-space with Neumann or Dirichlet boundary conditions, by bootstrapping the two point function of the fundamental field $\phi$. These boundary conditions preserve conformality on the defect and are usually called in the literature the *special* and *ordinary* transitions respectively.<fnref n="21" />

The above critical points have been studied in literature with a traditional RG approach, as formulated in [39] and extended in [40, 41] for the ordinary transition, and in [42] for the special transition. The main results in perturbation theory include the critical exponents related to the surface-surface and bulk-surface correlation of the field $\phi$, and the scaling of $\phi^2$ (order parameter or energy density) near the boundary up to an order $O(\epsilon^2)$. We direct the reader to [43] for a concise review.

<footnote n="21">²¹We will briefly comment in the next section on a third possible boundary condition, corresponding to the extraordinary transition.</footnote>
To obtain BCFT data one often considers the dimensions in bulk spectrum as input, and we will use the same principle in the BCFT bootstrap approach. The bulk dimensions are known from the extensive study of usual CFT without boundary [23, 44–46]. We will treat the Wilson-Fisher fixed point as a continuous deformation of the (generalized) free solution, i.e. one which preserves the structure of the spectrum of states in the two point function. This is certainly incorrect non-perturbatively since this spectrum is very different in both cases: for instance we expect that the number of operators appearing in the bulk channel increases exponentially fast for sufficiently large scaling dimension. In perturbation theory however, almost all such operators have coinciding scaling dimensions, and so we expect to get correct results for the $\langle\phi\phi\rangle$ correlator, up to the order where degeneracies among operators should get lifted. For the Wilson-Fisher BCFTs this turns out to first happen at $O(\epsilon^3)$. More precisely at this order one will require extra information to disentangle the anomalous dimensions of degenerate operators but one should still recover the full correlator.<fnref n="22" /> Here we will only consider perturbation theory up to $O(\epsilon^2)$.

The functional bootstrap equations take the form:

$$ \sum_{n=0}^{\infty} \mu_{\hat{\Delta}_n}^2 \omega_m(\hat{\Delta}_n | \Delta_{\phi}, \epsilon) = \omega_m(0 | \Delta_{\phi}, \epsilon) + \sum_{n=0}^{\infty} (a\lambda)_{\Delta_n} \omega_m(\Delta_n | \Delta_{\phi}, \epsilon) \quad (7.1) $$

where $\omega = \alpha, \beta, \theta$ and we've explicitly the dependence of the functional actions on both $\Delta_{\phi}$ and the spacetime dimension via $\epsilon$. We will expand parameters as follows:

$$ \begin{aligned} \hat{\Delta}_n &= \hat{\Delta}_n^{(0)} + \epsilon \hat{\Delta}_n^{(1)} + \epsilon^2 \hat{\Delta}_n^{(2)}, & \Delta_n &= \Delta_n^{(0)} + \epsilon \Delta_n^{(1)} + \epsilon^2 \Delta_n^{(2)}, \\ \mu_n^2 &= a_n^{(0)} + \epsilon a_n^{(1)} + \epsilon^2 a_n^{(2)}, & (a\lambda)_n &= t_n^{(0)} + \epsilon t_n^{(1)} + \epsilon^2 t_n^{(2)}, \\ \Delta_{\phi} &= \Delta_{\phi}^{(0)} + \epsilon \Delta_{\phi}^{(1)} + \epsilon^2 \Delta_{\phi}^{(2)}. \end{aligned} \quad (7.2) $$

with

$$ \hat{\Delta}_n^{(0)} = \frac{1-\nu}{2} + \Delta_{\phi}^{(0)} + 2n, \quad \Delta_n^{(0)} = 2\Delta_{\phi}^{(0)} + 2n, \quad \Delta_{\phi}^{(0)} = 1. \quad (7.3) $$

It is also convenient to define

$$ \hat{\gamma}_n^{(i)} = \hat{\Delta}_n^{(i)} - \Delta_{\phi}^{(i)}, \quad \gamma_n^{(i)} = \Delta_n^{(i)} - 2\Delta_{\phi}^{(i)}, \quad i = 1, 2, \dots \quad (7.4) $$

Note these are anomalous dimensions with respect to the generalized free field solution.

Upon inserting these expressions into the functional equations and expanding in $\epsilon$ we will obtain the functional actions and its derivatives with respect to $\Delta, \hat{\Delta}, \Delta_{\phi}$ and $\epsilon$ evaluated at $\Delta_{\phi} = 1$ and $\epsilon = 0$. We will omit this dependence below.

### $O(\epsilon^0)$

At this order we merely recover the free solution. In particular we get

$$ a_n^{(0)} = \alpha_n^{-}(0) = \delta_{n,0}, \quad t_n^{(0)} = -\theta_n^{-}(0) = -\delta_{n0}, \quad \text{(Dirichlet)} \quad (7.5) $$

$$ a_n^{(0)} = \alpha_n^{+}(0) = 2\delta_{n,0}, \quad t_n^{(0)} = -\theta_n^{+}(0) = \delta_{n0}, \quad \text{(Neumann)} \quad (7.6) $$

<footnote n="22">²²Beyond this order, it is not that the functional bootstrap equations do not fail, but rather they give rise to some other solution to crossing, namely one which has the sparsest possible spectrum of operators. That said, the functional bootstrap equations give solutions to crossing, not theories, so it is quite likely such a solution is actually unphysical.</footnote>
These equations tell us that the free solution to crossing contains one block in both boundary and bulk channels, on top of the bulk identity. The fact that $a_n^{(0)}$, $\theta_n^{(0)}$ vanish for $n \ge 1$ will lead to simplifications below. Notice that we have not imposed the free equation of motion at any point: rather the infinite decouplings of double trace and boundary derivative operators follows automatically from the functionals.

### $O(\epsilon)$

At this order we get:

$$ \sum_{n=0}^{\infty} \left\{ a_n^{(1)} \omega_m(\hat{\Delta}_n^{(0)}) + a_n^{(0)} \hat{\gamma}_n^{(1)} \partial_{\hat{\Delta}} \omega_m(\hat{\Delta}_n^{(0)}) \right\} = D_{\epsilon} \omega_m(0) + \sum_{n=0}^{\infty} \left\{ t_n^{(1)} \omega_m(\Delta_n^{(0)}) + t_n^{(0)} \gamma_n^{(1)} \partial_{\Delta} \omega_m(\Delta_n^{(0)}) \right\}. \quad (7.7) $$

In order to derive this expression we have used the orthonormality relations satisfied by $\alpha_m, \beta_m, \theta_m$ to simplify results, for instance, $\partial_{\Delta_{\phi}} \omega_m(\Delta_n^{(0)}) = -2\partial_{\Delta} \omega_m(\Delta_n^{(0)})$. Using the results from the previous order we can get:

$$ a_m^{(1)} - D_{\epsilon}\alpha_m(0) = t_0^{(0)}\gamma_0^{(1)}\partial_{\Delta}\alpha_m(\hat{\Delta}_0^{(0)}) - a_0^{(0)}\hat{\gamma}_0^{(1)}\partial_{\hat{\Delta}}\alpha_m(\hat{\Delta}_0^{(0)}) \quad (7.8a) $$

$$ t_m^{(1)} + D_{\epsilon}\theta_m(0) = -t_0^{(0)}\gamma_0^{(1)}\partial_{\Delta}\theta_m(\hat{\Delta}_0^{(0)}) + a_0^{(0)}\hat{\gamma}_0^{(1)}\partial_{\hat{\Delta}}\theta_m(\hat{\Delta}_0^{(0)}) \quad (7.8b) $$

$$ 0 = t_0^{(0)}\gamma_0^{(1)}\partial_{\Delta}\beta_m(\hat{\Delta}_0^{(0)}) - a_0^{(0)}\hat{\gamma}_0^{(1)}\partial_{\hat{\Delta}}\beta_m(\hat{\Delta}_0^{(0)}) \quad (7.8c) $$

where we have defined:

$$ D_{\epsilon} = \Delta_{\phi}^{(1)} \partial_{\Delta_{\phi}} + \partial_{\epsilon}. \quad (7.9) $$

On the other hand by explicit computation we find the relations:

$$ \partial_{\Delta} \alpha_m(\Delta_0^{(0)}) = -\partial_{\hat{\Delta}} \alpha_m(\hat{\Delta}_0^{(0)}), \quad m \ge 1 \quad (7.10) $$

$$ \partial_{\Delta} \theta_m(\Delta_0^{(0)}) = -\partial_{\hat{\Delta}} \theta_m(\hat{\Delta}_0^{(0)}), \quad m \ge 2 \quad (7.11) $$

$$ \partial_{\Delta} \beta_m(\Delta_0^{(0)}) = -\partial_{\hat{\Delta}} \beta_m(\hat{\Delta}_0^{(0)}), \quad m \ge 1. \quad (7.12) $$

These expressions are valid for both Dirichlet and Neumann functionals. In the former case, the r.h.s. of these equations is automatically zero. Recall also that all derivatives are evaluated for $\Delta_{\phi} = 1, \epsilon = 0$. We find therefore

$$ a_m^{(1)} = 0, \quad \text{for } m \ge 1, \quad (7.13) $$

$$ t_m^{(1)} = 0, \quad \text{for } m \ge 2, \quad (7.14) $$

for either boundary condition, where we have used that for the Wilson-Fisher fixed point we have $\Delta_{\phi}^{(1)} = -\frac{1}{2}$. As for the remaining quantities we get

$$ a_0^{(1)} = \frac{1}{2}\gamma_0^{(1)} - \frac{1}{2}, \quad t_0^{(1)} = \frac{1}{2}\gamma_0^{(1)}, \quad t_1^{(1)} = \frac{1}{4}\gamma_0^{(1)}, \quad \hat{\gamma}_0^{(1)} = -\frac{1}{2}\gamma_0^{(1)}, \quad (\text{Dirichlet}) \quad (7.14) $$

$$ a_0^{(1)} = 0, \quad t_0^{(1)} = \frac{1}{2}\gamma_0^{(1)}, \quad t_1^{(1)} = \frac{1}{4}\gamma_0^{(1)}, \quad \hat{\gamma}_0^{(1)} = -\frac{1}{2}\gamma_0^{(1)}, \quad (\text{Neumann}) \quad (7.14) $$

This is in perfect agreement with the results of [3, 4].
• O(ε²)

At this order the equations become more involved. We obtain:

$$
\begin{equation}
\begin{aligned}
& \sum_{n=0}^{\infty} \left\{ a_n^{(2)} \omega_m(\hat{\Delta}_n^{(0)}) + (a_n^{(0)} \hat{\gamma}_n^{(2)} + a_n^{(1)} \hat{\gamma}_n^{(1)}) \partial_{\hat{\Delta}} \omega_m(\hat{\Delta}_n^{(0)}) + \right. \\
& \qquad \left. + \frac{1}{2} a_n^{(0)} \left[ (\hat{\Delta}_n^{(1)})^2 - (\Delta_{\phi}^{(1)})^2 \right] \partial_{\hat{\Delta}}^2 \omega_m(\hat{\Delta}_n^{(0)}) + a_n^{(0)} \hat{\gamma}_n^{(1)} D_{\epsilon} \partial_{\hat{\Delta}} \omega_m(\hat{\Delta}_n^{(0)}) \right\} \\
& = \left( \frac{1}{2} D_{\epsilon}^2 + \Delta_{\phi}^{(2)} \partial_{\Delta_{\phi}} \right) \omega_m(0) + \sum_{n=0}^{\infty} \left\{ t_n^{(2)} \omega_m(\Delta_n^{(0)}) + (t_n^{(0)} \gamma_n^{(2)} + t_n^{(1)} \gamma_n^{(1)}) \partial_{\Delta} \omega_m(\Delta_n^{(0)}) \right. \\
& \qquad \left. + \frac{1}{2} t_n^{(0)} \left[ (\Delta_n^{(1)})^2 - 4(\Delta_{\phi}^{(1)})^2 \right] \partial_{\Delta}^2 \omega_m(\Delta_n^{(0)}) + t_n^{(0)} \gamma_n^{(1)} D_{\epsilon} \partial_{\Delta} \omega_m(\Delta_n^{(0)}) \right\}
\end{aligned}
\tag{7.15}
\end{equation}
$$

where again we used the orthonormality properties of $\omega_m$ to simplify things. Specializing
to the $\alpha, \beta, \theta$ functionals and doing some simplifications we get

$$
\begin{align}
& a_m^{(2)} + \frac{1}{2} a_0^{(0)} (\hat{\gamma}_0^{(1)})^2 \partial_{\hat{\Delta}}^2 \alpha_m(\hat{\Delta}_0^{(0)}) - \nonumber \\
& \quad - (a_0^{(0)} \hat{\gamma}_0^{(2)} + a_0^{(1)} \hat{\gamma}_0^{(1)}) d_m - a_0^{(0)} \hat{\gamma}_0^{(1)} D_{\epsilon} d_m \tag{7.16} \\
&= \Delta_{\phi}^{(2)} \partial_{\Delta_{\phi}} \alpha_m(0) + (t_0^{(0)} \gamma_0^{(2)} + t_0^{(1)} \gamma_0^{(1)}) \partial_{\Delta} \alpha_m(\Delta_0^{(0)}) + t_1^{(1)} \gamma_1^{(1)} \partial_{\Delta} \alpha_m(\Delta_1^{(0)}) \nonumber \\
& \quad + \frac{1}{2} t_0^{(0)} \left[ (\Delta_0^{(1)})^2 - 4(\Delta_{\phi}^{(1)})^2 \right] \partial_{\Delta}^2 \alpha_m(\Delta_0^{(0)}) + t_0^{(0)} \gamma_0^{(1)} D_{\epsilon} \partial_{\Delta} \alpha_m(\Delta_0^{(0)}) \nonumber
\end{align}
$$

$$
\begin{align}
& a_m^{(0)} \hat{\gamma}_m^{(2)} + a_m^{(1)} \hat{\gamma}_m^{(1)} + \frac{1}{2} a_0^{(0)} (\hat{\gamma}_0^{(1)})^2 \partial_{\hat{\Delta}}^2 \beta_m (\hat{\Delta}_0^{(0)}) - \nonumber \\
& \quad - (a_0^{(0)} \hat{\gamma}_0^{(2)} + a_0^{(1)} \hat{\gamma}_0^{(1)}) e_m - a_0^{(0)} \hat{\gamma}_0^{(1)} D_\epsilon e_m \tag{7.17} \\
&= (t_0^{(0)} \gamma_0^{(2)} + t_0^{(1)} \gamma_0^{(1)}) \partial_{\Delta} \beta_m (\Delta_0^{(0)}) + t_1^{(1)} \gamma_1^{(1)} \partial_{\Delta} \beta_m (\Delta_1^{(0)}) \nonumber \\
& \quad + \frac{1}{2} t_0^{(0)} \left[ (\Delta_0^{(1)})^2 - 4(\Delta_\phi^{(1)})^2 \right] \partial_{\Delta}^2 \beta_m (\Delta_0^{(0)}) + t_0^{(0)} \gamma_0^{(1)} D_\epsilon \partial_{\Delta} \beta_m (\Delta_0^{(0)}) \nonumber
\end{align}
$$

$$
\begin{align}
& \frac{1}{2} a_0^{(0)} (\hat{\gamma}_0^{(1)})^2 \partial_{\hat{\Delta}}^2 \theta_m (\hat{\Delta}_0^{(0)}) - \nonumber \\
& \quad - (a_0^{(0)} \hat{\gamma}_0^{(2)} + a_0^{(1)} \hat{\gamma}_0^{(1)}) f_m - a_0^{(0)} \hat{\gamma}_0^{(1)} D_{\epsilon} f_m \tag{7.18} \\
& = \Delta_{\phi}^{(2)} \partial_{\Delta_{\phi}} \theta_m (0) + t_m^{(2)} + (t_0^{(0)} \gamma_0^{(2)} + t_0^{(1)} \gamma_0^{(1)}) \partial_{\Delta} \alpha_m (\Delta_0^{(0)}) + t_1^{(1)} \gamma_1^{(1)} \partial_{\Delta} \theta_m (\Delta_1^{(0)}) \nonumber \\
& \quad + \frac{1}{2} t_0^{(0)} \left[ (\Delta_0^{(1)})^2 - 4(\Delta_{\phi}^{(1)})^2 \right] \partial_{\Delta}^2 \theta_m (\Delta_0^{(0)}) + t_0^{(0)} \gamma_0^{(1)} D_{\epsilon} \partial_{\Delta} \theta_m (\Delta_0^{(0)}) \nonumber
\end{align}
$$

In each of these three equations, note that the second line is only non-zero for the Neumann
boundary condition. Recall that the coefficients $d_m, e_m, f_m$ provide the conformal block
decomposition of the simplest contact interaction with Neumann boundary conditions, cf.
appendix C.3. In this case, the second equation is also only valid for $m \ge 1$, as $\beta_0^+ = 0$.

Determining the second order coefficients is now a matter of evaluating derivatives of
functional actions at specific values. In some cases this can be done analytically, but for
the bulk functional actions especially this is difficult. For instance from the second set of equations one finds for Dirichlet boundary conditions:

$$ \hat{\gamma}_0^{(2)} = -\gamma_0^{(2)} + \frac{1}{4}\gamma_0^{(1)}(\gamma_0^{(1)} - 1) \quad (\text{Dirichlet}) \quad (7.19) $$

For the purposes of this work, we were content to check analytically where possible, and more generally numerically, that the results following from these equations are in perfect agreement with those of [28].

### * O(ε³) and beyond

It is easy to see that there are no obstacles to going to higher orders in ε expansion. In particular at p-th order the equations look schematically as

$$ \sum_{n=0}^{\infty} \left\{ a_n^{(p)} \omega_m(\hat{\Delta}_n^{(0)}) + a_n^{(0)} \hat{\gamma}_n^{(p)} \partial_{\hat{\Delta}} \omega_m(\hat{\Delta}_n^{(0)}) + (\dots) \right\} \qquad (7.20) \\ = \Delta_{\phi}^{(p)} \partial_{\Delta_{\phi}} \omega_m(0) + \sum_{n=0}^{\infty} \left\{ t_n^{(p)} \omega_m(\Delta_n^{(0)}) + t_n^{(0)} \gamma_n^{(p)} \partial_{\Delta} \omega_m(\Delta_n^{(0)}) + (\dots) \right\} $$

where on both sides of the equation the dots involve only results that are lower order in ε. In particular, hidden in the dots we have contributions from infinite sets of operators in both bulk and boundary channels starting from p = 3. Since the functionals satisfy orthornormality properties, it follows that by choosing ω = α, β, θ we can isolate specific a_m(p), γ_m(p) or t_m(p) in the infinite sums above, e.g.:

$$ a_n^{(p)} = \Delta_{\phi}^{(p)} \partial_{\Delta_{\phi}} \omega_m(0) + \sum_{n=0}^{+\infty} t_n^{(0)} \gamma_n^{(p)} \partial_{\Delta} \alpha_m(\Delta_n^{(0)}) + (\dots) \quad (7.21) $$

The order εp data will be expressed only in terms of results from previous orders in ε, as well as the infinite set of bulk scaling dimensions γn(p) and Δφ(p) which must be given as inputs.

## 8 Discussion and outlook

In this work we have constructed two functional bases for the BCFT crossing equation, which are associated to generalized free field solutions with Neumann and Dirichlet boundary conditions. The functionals lead to a set of sum rules, the functional bootstrap equations, which constrain the BCFT data. We have explained the relation of these equations to a Polyakov-like approach to the BCFT bootstrap, and used them to rederive perturbation theory results for the Wilson-Fisher fixed point to O(ε²). There are a number of open problems and possible applications and extensions of these results.

The functional sum rules derived in this paper form a set of necessary conditions on the BCFT data. It would be nice to see if they can be used to derive bounds on this data as was done for 1D CFTs in [19]. We expect this to be possible at least in the boundary channel where there is manifest positivity. Such bounds are crucial if we want to show that
the equations are also sufficient for crossing. In particular, until we do so, our proof of the validity of the Polyakov bootstrap will remain incomplete.

The sum rules lend themselves nicely not only to perturbative but also numerical studies of the crossing equation. This can be done by truncating the sum over the spectrum and considering only a finite set of functional equations. One can think of this as applying Gliozzi's [15] (or extremal flows [16]) method to the BCFT crossing equation in a particularly nice basis. Since the same basis reproduces the correct $\epsilon$-expansion coefficients, at least to some order, it is reasonable to expect that repeating the analysis numerically should also give good results. Right now the chief technical obstacle both for perturbative and numerical analysis is better control over the functional actions.

We have discussed two kinds of bases for the BCFT crossing equation and their associated functionals, in sections 3 and 4. For both kinds, generic dimension bulk and boundary blocks were expanded in terms of double trace bulk blocks and boundary derivative blocks. In the first, “easy” basis, the boundary channel was comprised of blocks $\hat{G}_{\hat{\Delta}_n}$ of dimension $\Delta_n = \Delta_{\hat{\phi}} + n$ for $n = 0, 1, \dots$, whereas in the main bases of section 4 we had $\Delta_n = \Delta_{\hat{\phi}} + 2n$. This halving was compensated by the appearance of derivatives of boundary blocks $\partial_{\hat{\Delta}}\hat{G}_{\hat{\Delta}_n}$. It was thanks to the functionals dual to these latter objects (i.e. the $\beta_n$) that we were able to control deformations of the generalized free solution where the boundary operator dimensions were allowed to change.

In both kinds of bases however, the bulk channel is comprised of blocks with dimension $G_{\Delta_n}$ with $\Delta_n = 2\Delta_\phi + 2n$, but not their derivatives. This implies that in this basis we do not have good control over deformations where the bulk operator dimensions change, which means in practice these must be inputs. It seems then that there should be a third kind of basis, where also in the bulk channel we halve the set of operators and introduce compensating $\partial_{\Delta}G_{\Delta_n}$. In general however, we expect that such a basis should be very non-trivial, since it would not be associated to generalized free fields. The reason why we should expect this basis to exist at all is because the associated functionals can be constructed numerically: indeed, the original numerical bootstrap for BCFTs implemented in [3] can be thought of as constructing exactly such functionals. It would be very interesting to see what can be said of such basis, as they might then allow us to access scaling dimensions of higher-d bulk CFTs.

In this paper we have constructed explicit expressions for Witten diagrams for BCFTs, subject to Neumann and Dirichlet boundary conditions. It would be worthwhile to explore their structure more, to see if there are further simplifications and interesting features in their block decomposition [27, 47].

We restricted the application of the functionals to the perturbative analysis of a Wilson-Fisher CFT with a boundary with Neumann and Dirichlet boundary conditions. On similar lines one may investigate other boundary CFT examples that allow a perturbative approach. One straightforward but interesting generalization to our analysis is the case of extraordinary transition for Wilson-Fisher theory in $d = 4 - \epsilon$ dimensions. As shown in appendix B.4 of [3], it is defined by the existence of a one-point function for the bulk field $\phi$ i.e. the existence of a boundary identity operator. We leave study of this interesting case for the future.
Let us conclude by noting there are several possible generalizations of the results in this paper which should be obtainable with moderate efforts. Among these stand out bootstrapping two point functions of spinning fields [48] and including various amounts of supersymmetry [49]. The latter is particularly interesting in light of applications to line defects in superconformal CFTs where many exact results are available. We look forward to exploring these and other developments in the future.

## Acknowledgments

We would like to thank Adel Bilal, Parijat Dey, Kausik Ghosh, Rajesh Gopakumar, Dalimil Mazáč, Leonardo Rastelli, Slava Rychkov, Aninda Sinha, Emilio Trevisani and Xinan Zhou for useful discussions and/or comments on the draft. AK and MFP acknowledge the role of the Simons bootstrap collaboration workshops in Azores and Caltech where parts of this work were completed. AK would like to thank YITP, Stony Brook for hospitality during the course of this work. AK is supported by the Simons Foundation grant 488655 (Simons Collaboration on the Nonperturbative Bootstrap).

## A Functional actions

In this section we compute the functional actions starting from the expressions for the kernels determined in section 4

### A.1 α and β functionals

The α and β functionals can be obtained from the $\omega_{m,k}$ by Gram-Schmidt decomposition. Concretely, we have:

$$ \beta_n^- = \frac{\omega_{2n,0} - \sum_{i=0}^{n-1} \partial_{\hat{\Delta}} \omega_{2n,0}(\hat{\Delta}_i^-) \beta_i^-}{\partial_{\hat{\Delta}} \omega_{2n,0}(\hat{\Delta}_n^-)} \quad (A.1) $$

$$ \beta_n^+ = \frac{\omega_{2n-1,0} - \sum_{i=1}^{n-1} \partial_{\hat{\Delta}} \omega_{2n-1,0}(\hat{\Delta}_i^+) \beta_i^+}{\partial_{\hat{\Delta}} \omega_{2n-1,0}(\hat{\Delta}_n^+)} \quad (A.2) $$

$$ \alpha_n^- = \frac{\omega_{2n+2,2n+2} - \sum_{i=0}^{n+1} \partial_{\hat{\Delta}} \omega_{2n+2,2n+2}(\hat{\Delta}_i^-) \beta_i^- - \sum_{i=1}^{n-1} \omega_{2n+2,2n+2}(\hat{\Delta}_i^-) \alpha_i^-}{\omega_{2n+2,2n+2}(\hat{\Delta}_n^-)} \quad (A.3) $$

$$ \alpha_n^+ = \frac{\omega_{2n+1,2n+1} - \sum_{i=1}^{n+1} \partial_{\hat{\Delta}} \omega_{2n+1,2n+1}(\hat{\Delta}_i^+) \beta_i^+ - \sum_{i=1}^{n-1} \omega_{2n+1,2n+1}(\hat{\Delta}_i^+) \alpha_i^-}{\omega_{2n+1,2n+1}(\hat{\Delta}_n^+)} \quad (A.4) $$

with $\hat{\Delta}_n^- = 1 + \Delta_\phi + 2n$ and $\hat{\Delta}_n^+ = \Delta_\phi + 2n$.

#### A.1.1 Action on boundary blocks

One finds:

$$ \omega_{m,k}(\hat{\Delta}) = - \frac{4^{-\Delta_{\phi}+\Delta-m-\frac{1}{2}} \Gamma(\Delta+\frac{3-d}{2}) \Gamma(k+\Delta_{\phi}+\frac{3-d}{2}) \Gamma(m+\Delta_{\phi}+1)}{(m-k)! \Gamma(\Delta) \Gamma(\frac{m-\Delta+\Delta_{\phi}+3}{2}) \Gamma(\frac{m+\Delta+\Delta_{\phi}+4-d}{2}) \Gamma(\frac{2k-m-\Delta+\Delta_{\phi}+1}{2}) \Gamma(\frac{2k-m+\Delta+\Delta_{\phi}+2-d}{2})} \quad (A.5) $$
Using this expression we find

$$ \beta_n^+(\hat{\Delta}) = \tilde{\beta}_n^+(\hat{\Delta}) - e_n \tilde{\beta}_0^+(\hat{\Delta}) \quad (A.6) $$

where

$$ \tilde{\beta}_n^+(\hat{\Delta}) = \frac{4^{\hat{\Delta}-\hat{\Delta}_n+\frac{3}{2}} \Gamma(-\frac{d}{2}+\hat{\Delta}+\frac{3}{2}) \Gamma(\hat{\Delta}_n) \Gamma(-\frac{d+\hat{\Delta}_n+\hat{\Delta}_0+1}{2})^2}{\Gamma(\hat{\Delta})(\hat{\Delta}-\hat{\Delta}_n)(-d+\hat{\Delta}+\hat{\Delta}_n+1) \Gamma(\frac{\hat{\Delta}_0-\hat{\Delta}}{2})^2 \Gamma(-\frac{d+\hat{\Delta}+\hat{\Delta}_0+1}{2})^2 \Gamma(-\frac{d}{2}+\hat{\Delta}_n+\frac{1}{2}) \Gamma(\frac{\hat{\Delta}_n-\hat{\Delta}_0+2}{2})^2} \quad (A.7) $$

In this expression $\hat{\Delta}_n = \Delta_\phi + 2n$. Recall that the $\tilde{\beta}$, $\tilde{\alpha}$ are the fully orthonormal Neumann prefunctionals, i.e. functionals whose kernels fall off as $z^{-\frac{1}{2}}$ at infinity, and not as $z^{-\frac{3}{2}}$ as they should. The $e_n$ are found to be

$$ e_n = \frac{4^{\Delta_\phi - \hat{\Delta}_n} \Gamma(\hat{\Delta}_n) \Gamma\left(\frac{-d + \hat{\Delta}_n + \Delta_\phi + 1}{2}\right)^2}{\Gamma(\Delta_\phi) \Gamma\left(-\frac{d}{2} + \Delta_\phi + \frac{1}{2}\right) \Gamma\left(-\frac{d}{2} + \hat{\Delta}_n + \frac{1}{2}\right) \Gamma\left(\frac{\hat{\Delta}_n - \Delta_\phi + 2}{2}\right)^2} \quad (A.8) $$

The Dirichlet functionals are then obtained simply by

$$ \beta_n^-(\hat{\Delta}|\Delta_\phi) = \tilde{\beta}_n^+(\hat{\Delta}|\Delta_\phi + 1). \quad (A.9) $$

In particular this leads to exactly the same expression as the r.h.s. of (A.7) but now with $\hat{\Delta}_n = 1 + \Delta_\phi + 2n$.

As for the $\alpha_n$ functionals, after computation we observe that:

$$ \alpha_n = \frac{1}{2} \partial_n \beta_n. \quad (A.10) $$

This holds for tilded and untilded Neumann functional actions as well as Dirichlet. In particular we also get

$$ d_n = \frac{1}{2} \partial_n e_n. \quad (A.11) $$

### A.1.2 Acting on bulk blocks

In this case we find

$$ \omega_{m,k}(\Delta) = N_{m,k}(\Delta) \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_{\phi})\right] 4 \tilde{F}_3 \left[ \frac{\frac{\Delta - 2\Delta_{\phi} + 2}{2}}{4k - 2m + \Delta}, m + \frac{\Delta}{2} + 2, \frac{2-d}{2} + \Delta \middle| 1 \right] \quad (A.12) $$

with

$$ N_{m,k}(\Delta) = \frac{4^{k-m} \Gamma\left(\frac{2-d}{2} + \Delta\right) \Gamma\left(\frac{\Delta}{2} - \Delta_{\phi} + 1\right) \Gamma\left(k + \frac{\Delta+1}{2}\right) \Gamma\left(m + \Delta_{\phi} + 1\right)}{\pi^{3/2}(m-k)!} \quad (A.13) $$

After Gram-Schmidt decomposition we were able to find relatively simple expressions for the $\beta_n(\Delta)$:

$$ \beta_n^-(\Delta) = \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_\phi)\right] \sum_{k=0}^{n} b_{n,k}^-(\Delta) 4 \tilde{F}_3 \left[ \frac{\frac{\Delta - 2\Delta_\phi + 2}{2}}{-4k + \Delta}, \frac{4+4k+\Delta}{2}, \frac{2-d}{2} + \Delta \middle| 1 \right] \quad (A.14) $$

$$ \beta_n^+(\Delta) = \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_\phi)\right] \sum_{k=1}^{n} b_{n,k}^+(\Delta) 4 \tilde{F}_3 \left[ \frac{\frac{\Delta - 2\Delta_\phi + 2}{2}}{2-4k+\Delta}, \frac{\Delta}{2}, \frac{1+\Delta}{2}, \frac{2-d}{2} + \Delta \middle| 1 \right] \quad (A.15) $$
with

$$ b_{n,k}^{-}(\Delta) = \frac{(1+2k)}{4^{2n}\pi^{\frac{3}{2}}} \frac{(\Delta_{\phi} + \frac{1-d}{2})_{n-k}}{(\Delta_{\phi} + n + k + \frac{3-d}{2})_{n-k}} \frac{\Gamma(\hat{\Delta}_{n}^{-}) \Gamma(1-\Delta_{\phi} + \frac{\Delta}{2}) \Gamma(\frac{1+\Delta_{\phi}}{2}) \Gamma(\Delta + \frac{2-d}{2})}{(n-k)!(n+k)!} \quad (A.16) $$

$$ b_{n,k}^{+}(\Delta) = \frac{2k}{4^{2n-1}\pi^{\frac{3}{2}}} \frac{(\Delta_{\phi} + \frac{1-d}{2})_{n-k}}{(\Delta_{\phi} + n + k + \frac{1-d}{2})_{n-k}} \frac{\Gamma(\hat{\Delta}_{n}^{+}) \Gamma(1-\Delta_{\phi} + \frac{\Delta}{2}) \Gamma(\frac{1+\Delta_{\phi}}{2}) \Gamma(\Delta + \frac{2-d}{2})}{(n-k)!(n+k)!} \quad (A.17) $$

To obtain the tilded functionals one just needs to know that

$$ \tilde{\beta}_0 = \omega_{-1,-1}, \quad (A.18) $$

since $\tilde{\beta}_n = \beta_n + e_n \tilde{\beta}_0$.

As for the $\alpha_n(\Delta)$ functional actions, these two are given in terms of finite sums of $_4F_3(1)$ as above, but we were now not able to derive a closed form expression for their coefficients.

## A.2 Bulk functionals

The orthonormal bulk functionals may be obtained from the shifted bulk functionals by Gram-Schmidt orthonormalization, which in this case takes the form

$$ \theta_n^{-} = s\theta_n - \sum_{i=0}^{n-1} s\theta_n^{-}(\Delta_i)\theta_i^{-} \quad (A.19) $$

$$ \tilde{\theta}_n^+ = s\theta_n^+ - \sum_{i=0}^{n-1} s\theta_n^+(\Delta_i)\tilde{\theta}_i^+ \quad (A.20) $$

To recover the untilded functionals, i.e. those with the correct falloff in the Neumann case there is a slight modification:

$$ \theta_n^+ = s\theta_n^+ - \sum_{i=0}^{n-1} s\theta_n^+(\Delta_i)\theta_i^+, \quad n \ge 1 \quad (A.21) $$

$$ \theta_0^+ = \tilde{\theta}_0^+ - f_0 \tilde{\beta}_0, \quad \tilde{\theta}_0^+ = s\theta_0^+ \quad (A.22) $$

where $f_0$ was previously determined to be

$$ f_0 = \frac{\sqrt{\pi}\Gamma(\Delta_{\phi})}{\Gamma(\frac{1}{2} + \Delta_{\phi})} \quad (A.23) $$

Note we must have also

$$ \theta_n^+ = \tilde{\theta}_n^+ - f_n \tilde{\beta}_0 \quad (A.24) $$

where

$$ f_n = \frac{(-1)^n \sqrt{\pi} (2\Delta_{\phi} + 1 - d) \Gamma(n + \Delta_{\phi})^2 (\Delta_{\phi} - h + \frac{3}{2})_{n-1}}{2(n!) \Gamma(\Delta_{\phi}) \Gamma(n + \Delta_{\phi} + \frac{1}{2}) (n + 2\Delta_{\phi} - h)_n}, \quad (A.25) $$
is determined in appendix C.3. Consistency requires

$$ \sum_{i=0}^{n} s\theta_{n}^{+}(\Delta_{i}) f_{i} = 0 \qquad (A.26) $$

and we have checked this is the case numerically to great accuracy.

As for the shifted functional actions we find

$$ s\theta_{n}^{-}(\Delta) = \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_{\phi})\right] \mathcal{N}_{n}^{-} \sum_{p=0}^{+\infty} \sum_{i=0}^{1} t_{n,p,i}^{-}(\Delta) {}_3\tilde{F}_2 \left[ \begin{array}{l} p - \frac{1}{2}, \frac{\Delta+2-d}{2}, i + \frac{2\Delta_{\phi}+2n+2p-1-d}{2} \\ \frac{\Delta+1-d+p}{2}, i + \frac{2\Delta_{\phi}+\Delta+2n+2p-1-d}{2} \end{array} ; 1 \right] \qquad (A.27) $$

$$ s\theta_{n}^{+}(\Delta) = \sin\left[\frac{\pi}{2}(\Delta - 2\Delta_{\phi})\right] \mathcal{N}_{n}^{+} \sum_{p=0}^{+\infty} \sum_{i=0}^{1} t_{n,p,i}^{+}(\Delta) {}_3\tilde{F}_2 \left[ \begin{array}{l} p - \frac{3}{2}, \frac{\Delta+2-d}{2}, i + \frac{2\Delta_{\phi}+2n+2p-3-d}{2} \\ \frac{\Delta-1-d+p}{2}, i + \frac{2\Delta_{\phi}+\Delta+2n+2p-3-d}{2} \end{array} ; 1 \right] \qquad (A.28) $$

and

$$ s\theta_{n}^{-}(\hat{\Delta}) = \sin^2\left[\frac{\pi}{2}(\Delta - \hat{\Delta}_{0}^{-})\right] \hat{\mathcal{N}}_{n}^{-} \sum_{p=0}^{+\infty} \sum_{i=0}^{1} \hat{t}_{n,p,i}^{-}(\hat{\Delta}) {}_3\tilde{F}_2 \left[ \begin{array}{l} \Delta, \Delta_{\phi} + \Delta + 1 - d, \frac{2\Delta + 2 - d}{2} \\ 2\Delta + 2 - d, i + \Delta_{\phi} + \Delta + n + p + \frac{1-d}{2} \end{array} ; 1 \right] \qquad (A.29) $$

$$ s\theta_{n}^{+}(\hat{\Delta}) = \sin^2\left[\frac{\pi}{2}(\Delta - \hat{\Delta}_{0}^{+})\right] \hat{\mathcal{N}}_{n}^{+} \sum_{p=0}^{+\infty} \sum_{i=0}^{1} \hat{t}_{n,p,i}^{+}(\hat{\Delta}) {}_3\tilde{F}_2 \left[ \begin{array}{l} \Delta, \Delta_{\phi} + \Delta - d, \frac{2\Delta + 2 - d}{2} \\ 2\Delta + 2 - d, i + \Delta_{\phi} + \Delta + n + p - \frac{1-d}{2} \end{array} ; 1 \right] \qquad (A.30) $$

The coefficients appearing in these and the above expressions can be obtained from the authors upon request.

## B Computation of bulk exchange Witten diagram

In this appendix we give the details of the computation of the bulk Witten diagram (5.24), which we show here for convenience,

$$ W_{\Delta}^{\nu}(P_1, P_2) = \int_{AdS_{d+1}} dX \int_{AdS_d} dW \left[ G_{B\partial}(P_1, X|\Delta_1) G_{B\partial}(P_2, X|\Delta_2) G_{BB}(X, W|\Delta) \right. \\ \left. + \nu G_{B\partial}(P_{1r}, X|\Delta_1) G_{B\partial}(P_2, X|\Delta_2) G_{BB}(X, W|\Delta) \right]. \qquad (B.1) $$

The first part of (5.24) is the same bulk exchange integral evaluated in [33], for the interface boundary condition. The second integral is slightly different, so let us work it out in detail. The bulk-boundary propagator $G_{B\partial}(P, X)$ is given by (5.4). The bulk-bulk propagator in its spectral representation is given by,

$$ G_{BB}(X, W|\Delta) = \int_{-i\infty}^{i\infty} \frac{dc}{(c^2 - (\Delta - h)^2)} \frac{\Gamma(h+c)\Gamma(h-c)}{2\pi^{2h}\Gamma(c)\Gamma(-c)} \int_{\partial AdS_{d+1}} dP (-2P.X)^{h+c} (-2P.W)^{h-c}. \qquad (B.2) $$
Carrying out the integral over $X$, which is over bulk $AdS_{d+1}$ gives us the following three-point function,

$$ \langle O_1(P_1)O_2(P_{2r})O_{h+c}(P) \rangle = \frac{\pi^h \Gamma\left(\frac{\Delta_1+\Delta_2-h-c}{2}\right) \Gamma\left(\frac{\Delta_1+\Delta_2-h+c}{2}\right) \Gamma\left(\frac{\Delta_1-\Delta_2+h+c}{2}\right) \Gamma\left(\frac{\Delta_2-\Delta_1+h+c}{2}\right)}{2\Gamma(\Delta_1)\Gamma(\Delta_2)\Gamma(h+c)(-2P_1.P_{2r})^{\frac{\Delta_1+\Delta_2-h-c}{2}}(-2P_1.P)^{\frac{h+c+\Delta_{12}}{2}}(-2P_{2r}.P)^{\frac{h+c+\Delta_{21}}{2}}}. \quad (B.3) $$

Here we have

$$ \begin{aligned} -2P_1.P_{2r} &= (\vec{x}_1 - \vec{x}_2)^2 + (x_1^{\perp} + x_2^{\perp})^2 \\ -2P_1.P &= (\vec{x}_1 - \vec{x})^2 + (x_1^{\perp} - x^{\perp})^2 \\ -2P_{2r}.P &= (\vec{x}_2 - \vec{x})^2 + (x_2^{\perp} + x^{\perp})^2. \end{aligned} \quad (B.4) $$

Carrying out the $W$ integral gives the one point function,

$$ \langle O_{h-c}(P) \rangle = \frac{\pi^{h-\frac{1}{2}} \Gamma\left(\frac{h-c}{2}\right) \Gamma\left(\frac{-h-c+1}{2}\right)}{2\Gamma(h-c)(x^{\perp})^{\frac{h-c}{2}}}. \quad (B.5) $$

The integral over $P$ can be done with Schwinger parameters, that yields,

$$ \int_{-\infty}^{\infty} \frac{dx^{\perp} d^{d-1}x}{(-2P_{2r}.P)^{\frac{h+c+\Delta_1-\Delta_2}{2}} (-2P_1.P)^{\frac{h+c+\Delta_2-\Delta_1}{2}} (x^{\perp})^{\frac{h-c}{2}}} = \frac{\pi^h}{\Gamma\left(\frac{h+c+\Delta_{12}}{2}\right) \Gamma\left(\frac{h+c+\Delta_{21}}{2}\right) \Gamma\left(\frac{h-c}{2})} \\ \times \int_0^{\infty} \frac{ds}{s} \frac{dt}{t} \frac{du}{u} \frac{s^{\frac{h-c}{2}} t^{\frac{h+c+\Delta_{12}}{2}} u^{\frac{h+c+\Delta_{21}}{2}}}{(t+u)^{h-\frac{1}{2}} (s+t+u)^{\frac{1}{2}}} \exp \left[ -\frac{tu(-2P_1.P_{2r})}{t+u} - \frac{s(tx_1^{\perp} - ux_2^{\perp})}{(t+u)(s+t+u)} \right]. \quad (B.6) $$

Here $x_1^{\perp}, x_2^{\perp} > 0$. To evaluate the integral, the usual steps are to insert the following,

$$ 1 = \int d\lambda \, \delta(\lambda - (s + t + u)) \int d\rho \, \delta(\rho - (t + u)), \quad (B.7) $$

and do the $s$ integral. Then we rescale $\lambda = \rho\lambda$, $t = \rho t$ and $u = \rho u$ and carry out the $t$ integral. This gives,

$$ \int_{1}^{\infty} d\lambda \int_{0}^{1} du \int_{0}^{\infty} d\rho \, (\lambda - 1)^{\frac{h-c}{2}-1} (1-u)^{\frac{h+c+\Delta_{12}}{2}-1} u^{\frac{h+c+\Delta_{21}}{2}-1} \rho^{\frac{h-c}{2}-1} \\ \times \exp \left[ -\rho u(1-u)(-2P_1.P_{2r}) - \frac{\rho(\lambda-1)}{\lambda}((1-u)x_1^{\perp} - ux_2^{\perp})^2 \right]. \quad (B.8) $$

Note that the limits of the integrals follow from the above rescalings. Here we have to do an inverse Mellin transform, to rewrite one part inside the exponential as a power law. For this we use the following,

$$ e^{-x} = \int_{-i\infty}^{i\infty} \Gamma(\tau)x^{-\tau}d\tau. \quad (B.9) $$
We choose the arguments of exponentials and power laws carefully such that the integrals over $\lambda$ and $u$ can be carried out. Rearranging the exponential slightly we can write,

$$ 
\begin{aligned}
& \exp \left[ -\rho u(1-u)(-2P_1.P_{2r}) - \frac{\rho(\lambda-1)}{\lambda}((1-u)x_1^{\perp} + u x_2^{\perp})^2 + 4\frac{\rho(\lambda-1)}{\lambda}(u(1-u)x_1^{\perp}x_2^{\perp}) \right] \\
& = \exp \left[ -\rho u(1-u)(-2P_1.P_{2r}) + 4\frac{\rho(\lambda-1)}{\lambda}(u(1-u)x_1^{\perp}x_2^{\perp}) \right] \\
& \quad \times \int_{-i\infty}^{i\infty} d\tau \Gamma(\tau) \left( \frac{\rho(\lambda-1)}{\lambda}((1-u)x_1^{\perp} + u x_2^{\perp})^2 \right)^{-\tau}.
\end{aligned} 
\quad (B.10) $$

One can easily verify that the quantities inside exponential and in power of $\tau$ are of definite signs, which allows us to do the $\rho, u$ and $\lambda$ integrals. This gives a hypergeometric $_2F_1$. The final answer, including the first term of (5.24), is given by (5.25).

## C Witten diagram block decompositions

In this appendix we are going to give the details of the coefficients of expansion of the various Witten diagrams, and how to compute them. Since the Witten diagram expressions in section 5 are written as functions of $(1-z)/z$, it is convenient to use the variable $\xi$,

$$ \xi = \frac{1-z}{z}. \quad (C.1) $$

In terms of this, the bulk and boundary blocks are respectively given by,

$$ 
\begin{aligned}
G_{\Delta}(z|\Delta_1, \Delta_2) &= \xi^{\frac{\Delta-\Delta_1-\Delta_2}{2}} (1+\xi)^{\frac{\Delta_1+\Delta_2}{2}} {}_2F_1\left(\frac{\Delta+\Delta_1}{2}, \frac{\Delta+\Delta_2}{2}; \Delta-h+1, -\xi\right) \\
\hat{G}_{\hat{\Delta}}(z|\Delta_1, \Delta_2) &= \xi^{-\hat{\Delta}} (1+\xi)^{\frac{\Delta_1+\hat{\Delta}}{2}} {}_2F_1\left(\hat{\Delta}, \hat{\Delta}-h+1, 2\hat{\Delta}+2-d, -\frac{1}{\xi}\right).
\end{aligned} 
\quad (C.2) $$

It is clear that the limit $z \to 0$, which corresponds to boundary channel, will be given by $\xi \to \infty$. In the bulk channel, where $z \to 1$, we have $\xi \to 0$.

### C.1 Boundary exchange diagram

In the diagram (5.17) with Neumann boundary conditions, we have the following set of poles in $\tau$,

1. $\tau = h' \pm c + m_1$, for $m_1 \in \mathbb{Z}$ and $m_1 \ge 0$

2. $\tau = -m_2$, for $m_2 \in \mathbb{Z}$ and $m_2 \ge 0$

and the following poles in $c$,

1. $c = \pm(\hat{\Delta} - h')$

2. $c = \pm(\Delta_1 - h' + 2m_3)$ for $m_3 \in \mathbb{Z}$ and $m_3 \ge 0$

3. $c = \pm(\Delta_2 - h' + 2m_4)$ for $m_4 \in \mathbb{Z}$ and $m_4 \ge 0$

4. $c = \pm(h' - \tau + m_5)$, for $m_5 \in \mathbb{Z}$ and $m_5 \ge 0$
In the boundary channel decomposition (5.19), we have $\xi \to \infty$. Here the various $\xi$ dependencies arise from the collisions of $\tau$ and $c$ poles. The $\tau$-contour is closed on the right. To obtain the “physical term”, we take the residues at the poles $\tau = h' \pm c$. For the ‘+’ (‘-’) sign we close the $c$-contour on the right (left) and then choose $c = \hat{\Delta} - h'$ ($c = h' - \hat{\Delta}$). The total gives us the first term in the large $\xi$ expansion of the boundary block $\hat{G}_{\hat{\Delta}}((1+\xi)^{-1}|\Delta_1, \Delta_2)$. Since we normalize it to unit coefficient, this gives us $\mathcal{N}_{\text{bdy}}^{+}(\hat{\Delta})$ that reads,

$$ \mathcal{N}_{\text{bdy}}^{+}(\hat{\Delta}) = \frac{4^{\Delta+1}\Gamma(\hat{\Delta} + \frac{3}{2} - h)}{\Gamma(\hat{\Delta})\Gamma(\frac{\Delta_1-\hat{\Delta}}{2})\Gamma(\frac{\Delta_2-\hat{\Delta}}{2})\Gamma(\frac{\hat{\Delta}+\Delta_1+1}{2} - h)\Gamma(\frac{\hat{\Delta}+\Delta_2+1}{2} - h)}. \quad (C.3) $$

Let us now compute the coefficient of the boundary block $\hat{G}_{\Delta_1+2n}$, which is related to the coefficients of $\xi^{-\Delta_1-2n}$. It is obtained by taking residues at $\tau = h' \pm c+m_1$. Once again the choice of the $c$-contour depends on the sign. Then we take residues at $c = \Delta_1-h+2[n-\frac{m_1}{2}]$ with $m_1 = 1, 2, \dots, n$ (or $c = -\Delta_1+h'-2[n-\frac{m_1}{2}]$ as appropriate from the choice of contour). This gives all the contribution to $\xi^{-\Delta_1-2n}$, from which we subtract the contributions from the blocks $\hat{G}_{\Delta_1+2k}$ with $k < n$. This gives us the coefficient $a_n^+(\Delta_1, \Delta_2, \hat{\Delta})$,

$$ \hat{a}_n^+(\Delta_1, \Delta_2, \hat{\Delta}) = \mathcal{N}_{\text{bdy}}^+(\hat{\Delta}) \frac{(-1)^n \Gamma(2n + \Delta_1) \Gamma(\frac{\Delta_2-2n-\Delta_1}{2}) \Gamma(n + \Delta_1 + \frac{1}{2} - h) \Gamma(\frac{\Delta_1+\Delta_2+1}{2} - h + n)}{4^{\Delta_1+2n} n! (\hat{\Delta} - \Delta_1 - 2n) (\hat{\Delta} + \Delta_1 + 2n + 1 - d) \Gamma(2n + \Delta_1 + \frac{1}{2} - h)}. \quad (C.4) $$

The third set of poles in $c$ simply gives us $\hat{a}_n^+(\Delta_2, \Delta_1, \hat{\Delta})$ similarly as above. The fourth set has already been used up from the $\tau$ poles. For equal scalars the coefficients $\hat{b}_n^+$ and $\hat{a}_n^+$ can be obtained from (C.4) as shown in the main text.

In the bulk channel expansion of the diagram (5.17), we have $\xi \to 0$. So the $\tau$-contour is closed on left. The coefficient $\hat{t}_n^+$ of the bulk block $G_{\Delta_1+\Delta_2+2n}$ can be obtained by taking the residues at the poles $\tau = -m_2$ of $\Gamma(\tau)$ and rearranging the powers of $\xi$ as expected in the bulk blocks. The coefficient of $\xi^n$, coming from the residue at $\tau = -n$ is given by,

$$ \tilde{t}_n^+(\Delta_1, \Delta_2, \hat{\Delta}) = \mathcal{N}_{\text{bdy}}^+(\hat{\Delta}) \frac{(-4)^n \Gamma(h' + n)}{n! \Gamma(2h' + 2n)} \int_{-i\infty}^{i\infty} \frac{dc}{2\pi i} \frac{\hat{f}(c, -n) \hat{f}(-c, -n)}{(\hat{\Delta} - h')^2 - c^2}, \quad (C.5) $$

where we remind the reader,

$$ \hat{f}(c, \tau) = \frac{\Gamma\left(\frac{c+\Delta_1-h'}{2}\right) \Gamma\left(\frac{c+\Delta_2-h'}{2}\right) \Gamma(c+h'-\tau)}{2\Gamma(c)}. \quad (C.6) $$

We get $\hat{t}_n^+$, by subtracting the contributions to $\tilde{t}_n^+$ from lower blocks $G_{\Delta_1+\Delta_2+2k}$ with $k < n$. The general expression is too tedious to present, so we just indicate how to obtain them recursively,

$$ \hat{t}_n^+ = \tilde{t}_n^+ - \sum_{k=0}^{n-1} \left[ \frac{(-1)^{n-k} (\Delta_1 + k)_{n-k} (\Delta_2 + k)_{n-k}}{(n-k)!(\Delta_1 + \Delta_2 + 1 - h + 2k)_{n-k}} \right] \hat{t}_k^N, \quad (C.7) $$

and for $n=0$ it is simply $\hat{t}_0^+ = \tilde{t}_0^+$.
For the Dirichlet case, the analogous relations are simply obtained from the Neumann case as shown below,

$$ \mathcal{N}_{\text{bdy}}^{-}(\hat{\Delta}) = \left[ \mathcal{N}_{\text{bdy}}^{+}(\hat{\Delta}) \right]_{\substack{\Delta_1 \to \Delta_1+1 \\ \Delta_2 \to \Delta_2+1}} \qquad (C.8) $$

$$ \hat{a}_{n}^{-}(\Delta_{1}, \Delta_{2}, \hat{\Delta}) = \hat{a}_{n}^{+}(\Delta_{1} + 1, \Delta_{2} + 1, \hat{\Delta}), \qquad (C.9) $$

$$ \tilde{t}_{n}^{-}(\Delta_{1}, \Delta_{2}, \hat{\Delta}) = \tilde{t}_{n}^{+}(\Delta_{1} + 1, \Delta_{2} + 1, \hat{\Delta}). \qquad (C.10) $$

## C.2 Bulk exchange diagram

Here we give the details of the block expansion coefficients from the bulk exchange operators. The Witten diagram is given by (5.25). The following are the poles in $\tau$:

$$ 1. \tau = \frac{h \pm c}{2} + m_1, \text{ for } m_1 \in \mathbb{Z} \text{ and } m_1 \ge 0 $$

$$ 2. \tau = \pm \frac{\Delta_{12}}{2} - m_2, \text{ for } m_2 \in \mathbb{Z} \text{ and } m_2 \ge 0 $$

and the poles in $c$ are,

$$ 1. c = \pm(\Delta - h) $$

$$ 2. c = \pm(\Delta_1 + \Delta_2 - h + 2m_3) \text{ for } m_3 \in \mathbb{Z} \text{ and } m_3 \ge 0 $$

$$ 3. c = \pm(1 - h + 2m_4) \text{ for } m_4 \in \mathbb{Z} \text{ and } m_4 \ge 0 $$

$$ 4. c = \pm(h - 2\tau + 2m_5) \text{ for } m_5 \in \mathbb{Z} \text{ and } m_5 \ge 0. $$

In the bulk channel, we have $\xi \to 0$. So the $\tau$ contour is closed on the right. First we will compute the normalization $\mathcal{N}_{\text{bulk}}^{\nu}(\Delta)$. To get this, first we take the residue at the poles $\tau = \frac{h \pm c}{2} + m$,

$$ \text{Res}_{\tau = -\frac{h \pm c}{2} + m} [W_{\text{bulk}}^{\nu}] = \mathcal{N}_{\text{bulk}}^{\nu} \int \frac{dc \Gamma(\frac{1+c-h}{2}) \Gamma(\frac{\Delta_1+\Delta_2+c-h}{2}) \Gamma(\frac{2m+c+h}{2}) \Gamma(\frac{h+c+2m-\Delta_{12}}{2}) \Gamma(\frac{h+c+2m-\Delta_{21}}{2})}{(-1)^m 2 \times m! \Gamma(\pm c) (c^2 - (h-\Delta)^2) \Gamma(\frac{1 \mp c-h-2m}{2}) \Gamma(h \pm c+2m)} \\ \times f(-c, \frac{c+h+2m}{2}) (4\xi)^{\frac{h \pm c - \Delta_1 - \Delta_2 + 2m}{2}} \left[ 1 + \nu \frac{2F_1[\frac{1-c-h}{2}, \frac{1+c-h}{2}, \frac{1 \mp c-h-2m}{2}, z]}{(1-z)^{h-\frac{\Delta_1+\Delta_2+1}{2}}} \right]. \qquad (C.11) $$

Let us remind the reader that,

$$ f(c, \tau) = \frac{\Gamma(\frac{\Delta_1+\Delta_2-h+c}{2}) \Gamma(\frac{1+c-h}{2}) \Gamma(\frac{h+c}{2} - \tau)}{2\Gamma(c)}. \qquad (C.12) $$

The normalization $\mathcal{N}_{\text{bulk}}^{\nu}(\Delta)$ comes from requiring coefficient of $G_{\Delta}$ in the bulk channel decomposition to be 1. In the $\xi \to 0$ expansion of $G_{\Delta}((1+\xi)^{-1}|\Delta_1, \Delta_2)$ the first term is $\xi^{\frac{\Delta_1-\Delta_2}{2}}$. This is obtained from (C.11), by putting $m=0$ and taking the residue at $c=\Delta-h$ ($c=h-\Delta$) for the upper (lower) sign with the c-contour closed on right (left). This gives,

$$ \mathcal{N}_{\text{bulk}}^{\nu}(\Delta) = \frac{2^{\Delta_1+\Delta_2-\Delta+2} \times \Gamma(\Delta)\Gamma(\Delta-h+1)}{\Gamma(\frac{\Delta}{2})\Gamma(\frac{\Delta+\Delta_1-\Delta_2}{2})\Gamma(\frac{\Delta-\Delta_1+\Delta_2}{2})\Gamma(\frac{\Delta_1+\Delta_2-\Delta}{2})\Gamma(\frac{\Delta+1-d}{2})\Gamma(\frac{\Delta+\Delta_1+\Delta_2-d}{2})}. \qquad (C.13) $$
We will now look at the coefficients of the blocks $G_{\Delta_1+\Delta_2+2n}$. For simplicity of the expressions we will take $\Delta_1 = \Delta_2 = \Delta_\phi$, which is the case we are interested in anyway. First we have to obtain the power of $\xi^n$. This is done again from (C.11) with the following steps

**Step 1.** For the first term in parentheses, for the upper (lower) sign, we take the residue at $c = \Delta_1 + \Delta_2 - h + 2(n-m)$ ($c = -\Delta_1 - \Delta_2 + h - 2(n-m)$) and do the sums over $m = 0, 1, \dots, n$. Finally we expand in $\xi$ and get the coefficient of $\xi^n$.

**Step 2.** For the second term (i.e. the hypergeometric) we expand in $\xi$, and integrate the coefficient of $\xi^n$ over $c$ using residues. Then we sum over all $m \ge 0$.

Let us denote the total coefficient of $\xi^n$ obtained in this way by $\tilde{t}_n^\nu$. It can be expressed as an infinite sum of $_3F_2$ hypergeometric functions, but their expression is too tedious to display here. The coefficient of the block $t_n^\nu$ is then obtained by using a recursion formula like (C.7):

$$ -t_n^\nu = \tilde{t}_n^\nu + \sum_{k=0}^{n-1} \left[ \frac{(-1)^{n-k} ((\Delta_\phi + k)_{n-k})^2}{(n-k)!(2\Delta_\phi + 1 - h + 2k)_{n-k}} \right] t_k^\nu. \quad (C.14) $$

Let us now calculate the boundary channel decomposition coefficient $a_n(\Delta_1, \Delta_2, \Delta)$ of the boundary block $G_{\Delta_1+2n+\frac{\nu+1}{2}}$ from (5.27). The leading term of the block is $\xi^{-\Delta_1-2n-\frac{1-\nu}{2}}$ (times the prefactor $(1+\xi)^{\frac{\Delta_1+\Delta_2}{2}}$ coming from (C.2)). Since the limit is $\xi \to \infty$, the $\tau$ contour is closed on the left. A certain power $\xi^{-\Delta_1-m}$ is obtained by taking residues of $\tau = \frac{\Delta_2}{2} - k$, with $k = 0, 1, \dots, m$, then expanding the terms at large $\xi$ for each $k$ and summing. Let us call this coefficient $\tilde{a}_m^\nu$, where

$$ \tilde{a}_m^\nu = \sum_{k=0}^{m} \left[ \frac{\mathcal{N}_{\text{bulk}}^\nu}{2^{2\Delta_1+2k}} \int_{-i\infty}^{+i\infty} \frac{dc}{2\pi i} \left[ \left( \delta_{km} + \nu \sum_{p=0}^{m-k} \binom{h-\frac{\Delta_1+\Delta_2+1}{2}}{m-k-p} \left( \sum_{q=1}^{p} \frac{(-q) \binom{1-c-h}{2-q} \binom{1+c-h}{2-q}}{q! (k+\frac{\Delta_1+1}{2})_q} + \delta_{p,0} \right) \right) \right. \\ \left. \frac{f(c, \frac{\Delta_2}{2}-k) f(-c, \frac{\Delta_2}{2}-k)}{((\Delta-h)^2-c^2)} \right] \times \frac{(-1)^k \Gamma(\frac{\Delta_2}{2}-k) \Gamma(\Delta_2-k)}{k! \Gamma(\frac{1-\Delta_2}{2}+k) \Gamma(\Delta_2-2k)} \right]. \quad (C.15) $$

Integrals of this type can be evaluated using the following identity [27]

$$ \int_{-i\infty}^{i\infty} \frac{dc}{2\pi i} \frac{\prod_{i=1}^{3} \Gamma(a_i - \frac{c}{2}) \Gamma(a_i + \frac{c}{2})}{(4a_4^2 - c^2) \Gamma(c) \Gamma(-c)} = \frac{\Gamma(a_1 + a_2) \Gamma(a_1 + a_3) \Gamma(a_2 + a_4) \Gamma(a_3 + a_4)}{(a_1 + a_4) \Gamma(2a_4 + 1)} \\ \times {}_3F_2 \left( \begin{matrix} a_1 + a_4, a_4 - a_2 + 1, a_4 - a_3 + 1 \\ a_1 + a_4 + 1, 2a_4 + 1 \end{matrix} ; 1 \right). \quad (C.16) $$

This gives the following,

$$ \tilde{a}_m^\nu(\Delta_1, \Delta_2, \Delta) = \sum_{k=0}^{m} \frac{\mathcal{N}_{\text{bulk}}^\nu \Gamma(\frac{\Delta_2}{2}-k) \Gamma(\frac{\Delta+\Delta_1+\Delta_2}{2}-h) \Gamma(k+\frac{\Delta-\Delta_2}{2})}{2^{2\Delta_1+2k+1} (-1)^k k! (\Delta_2-k)_{-k} \Gamma(\Delta+1-h)} \left[ \left( \delta_{k,m} + \nu \binom{h-\frac{\Delta_1+\Delta_2+1}{2}}{m-k} \right) \right. \\ \times \frac{\Gamma(\frac{\Delta_1+\Delta_2+1-d}{2})}{(\Delta-d+1)} {}_3F_2 \left( \begin{matrix} \frac{\Delta+1-d}{2}, \frac{\Delta-\Delta_1-\Delta_2}{2}+1, \frac{\Delta+\Delta_2-1-d}{2}+1-k \\ \frac{\Delta+3-d}{2}, \Delta-h+1 \end{matrix} ; 1 \right) + \nu \sum_{p=1}^{m-k} \sum_{q=1}^{p} \binom{h-\frac{\Delta_1+\Delta_2+1}{2}}{m-k-p} \\ \left. \times \left[ \frac{(-q) (k-\frac{\Delta_2}{2}+\frac{1}{2})_q \Gamma(\frac{\Delta_1+\Delta_2+1}{2}-h+q)}{q! (k+\frac{\Delta_1+1}{2})_q (\Delta-2h+2q+1)} \right] {}_3F_2 \left( \begin{matrix} \frac{\Delta+1-d}{2}+q, \frac{\Delta-\Delta_1-\Delta_2}{2}+1, \frac{\Delta+\Delta_2-1-d}{2}+1-k \\ \frac{\Delta+3-d}{2}+q, \Delta-h+1 \end{matrix} ; 1 \right) \right] \quad (C.17) $$
The $a_n^\nu$-s are obtained by subtracting from $a_{2n+\frac{1-\nu}{2}}^\nu$ the coefficients of the lower boundary blocks $a_{m<n}^\nu$. We give $a_n^\nu$ in the recursion formula,

$$a_n^\nu(\Delta_1, \Delta_2, \Delta) = \tilde{a}_{2n+\frac{1-\nu}{2}}^\nu - \sum_{k=0}^{n-1} \left[ \frac{(\Delta_1 + 2k')_{2n-2k}(\Delta_1 - h + 2k' + 1)_{2n-2k}}{(2n-2k)!(2\Delta_1 + 4k' + 2-d)_{2n-2k}} \right] a_k^\nu(\Delta, \Delta_1, \Delta_2), \quad (C.18)$$

where $k' = k + \frac{1-\nu}{4}$.

### C.3 Contact diagram

The contact Witten diagram (5.33) is rather easy to decompose into boundary and bulk blocks. There is only the $\tau$-integral, whose contour is closed on the right for boundary channel, and on the left for bulk channel. As we expand around $\xi \to \infty$ and $\xi \to 0$ the coefficients of $\xi^{-\Delta_1-n}$ and $\xi^n$ are respectively given by,

$$\begin{gathered} \tilde{d}_n^\nu(\Delta_1, \Delta_2) = \frac{(-1)^n \Gamma(\Delta_{21} - n) \Gamma(n + \frac{1-\nu}{2} + \Delta_1)}{n! \Gamma(\frac{\Delta_{21}+1-2n}{2})}, \\ \text{and} \\ \tilde{f}_n^\nu(\Delta_1, \Delta_2) = \frac{\Gamma(n + \frac{1-\nu}{2} + \Delta_1) \Gamma(n + \frac{1-\nu}{2} + \Delta_2)}{n! \Gamma(\frac{2n+\Delta_1+\Delta_2+1}{2})}. \end{gathered} \quad (C.19)$$

From the above, obtaining the block decomposition coefficients is straightforward. For the Neumann case, we have the boundary channel expansion coefficients,

$$d_n^+(\Delta_1, \Delta_2) = \frac{(\Delta_1 + \Delta_2 + 1 - d) \Gamma(2n + \Delta_1) \Gamma(\Delta_{21} - 2n) \Gamma(n + \frac{1}{2} + \Delta_1 - h) \left( \frac{\Delta_1 + \Delta_2 + 3 - d}{2} \right)_{n-1}}{2^{2n+1} (-1)^n n! \Gamma(\frac{1+\Delta_{21}}{2} - n) \Gamma(2n + \Delta_1 + \frac{1}{2} - h)} \quad (C.20)$$

The corresponding coefficients for Dirichlet are simply given by

$$d_n^-(\Delta_1, \Delta_2) = d_n^+(\Delta_1 + 1, \Delta_2 + 1). \quad (C.21)$$

The bulk decomposition coefficients, for the Neumann case, are given by,

$$f_n^+(\Delta_1, \Delta_2) = \frac{(-1)^n (d - \Delta_1 - \Delta_2 - 1) \Gamma(n + \Delta_1) \Gamma(n + \Delta_2) \left(\frac{\Delta_1 + \Delta_2 + 3}{2} - h\right)_{n-1}}{2(n!) \Gamma\left(n + \frac{\Delta_1 + \Delta_2}{2} + \frac{1}{2}\right) \left(n + \Delta_1 + \Delta_2 - h\right)_n} \quad (C.22)$$

For Dirichlet the corresponding coefficients $f_n^-$ can be obtained from $\tilde{f}_n^-$, which are given by,

$$\tilde{f}_n^- = \tilde{f}_n^+(\Delta_1 + 1, \Delta_2 + 1). \quad (C.23)$$

One can approach the equal scalar limit, $\Delta_1 \to \Delta_\phi + \frac{\Delta_{12}}{2}$ and $\Delta_2 \to \Delta_\phi - \frac{\Delta_{12}}{2}$, to get

$$d_n^\nu(\Delta_1, \Delta_2) \stackrel{\Delta_{12} \to 0}{\rightleftharpoons} \frac{e_n^\nu}{\Delta_{12}} + \frac{1}{2} d_n^\nu + O(\Delta_{12}). \quad (C.24)$$

For later use, let us show these quantities for the Neumann case. Choosing the normalization $e_n = e_n^+/e_0^+$, $d_n = d_n^+/e_0^+$ and $f_n = f_n^+/e_0^+$, we get the following,

$$e_n = \frac{16^{-n} (\Delta_\phi)_{2n} \Gamma\left(n + \Delta_\phi - h + \frac{1}{2}\right)^2}{\Gamma(n+1)^2 \Gamma\left(\Delta_\phi - h + \frac{1}{2}\right) \Gamma\left(2n + \Delta_\phi - h + \frac{1}{2}\right)}, \quad \text{and} \quad d_n = \frac{1}{2} \partial_n e_n. \quad (C.25)$$
The bulk expansion coefficients are given by,

$$f_n = \frac{(-1)^n \sqrt{\pi}(2\Delta_\phi + 1 - d)\Gamma(n + \Delta_\phi)^2 (\Delta_\phi - h + \frac{3}{2})_{n-1}}{2(n!)\Gamma(\Delta_\phi)\Gamma(n + \Delta_\phi + \frac{1}{2}) (n + 2\Delta_\phi - h)_n} \quad (C.26)$$

**Open Access.** This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.

## References

[1] S. Ferrara, A.F. Grillo and R. Gatto, *Tensor representations of conformal algebra and conformally covariant operator product expansion*, *Annals Phys.* **76** (1973) 161 [INSPIRE].

[2] A.M. Polyakov, *Nonhamiltonian approach to conformal quantum field theory*, *Zh. Eksp. Teor. Fiz.* **66** (1974) 23 [Sov. Phys. JETP **39** (1974) 9] [INSPIRE].

[3] P. Liendo, L. Rastelli and B.C. van Rees, *The Bootstrap Program for Boundary CFT*, *JHEP* **07** (2013) 113 [arXiv:1210.4258] [INSPIRE].

[4] D.M. McAvity and H. Osborn, *Conformal field theories near a boundary in general dimensions*, *Nucl. Phys. B* **455** (1995) 522 [cond-mat/9505127] [INSPIRE].

[5] J.L. Cardy, *Conformal Invariance and Surface Critical Behavior*, *Nucl. Phys. B* **240** (1984) 514 [INSPIRE].

[6] R. Rattazzi, V.S. Rychkov, E. Tonni and A. Vichi, *Bounding scalar operator dimensions in 4D CFT*, *JHEP* **12** (2008) 031 [arXiv:0807.0004] [INSPIRE].

[7] F. Gliozzi, P. Liendo, M. Meineri and A. Rago, *Boundary and Interface CFTs from the Conformal Bootstrap*, *JHEP* **05** (2015) 036 [arXiv:1502.07217] [INSPIRE].

[8] M. Billò, V. Gonçalves, E. Lauria and M. Meineri, *Defects in conformal field theory*, *JHEP* **04** (2016) 091 [arXiv:1601.02883] [INSPIRE].

[9] M. Lemos, P. Liendo, M. Meineri and S. Sarkar, *Universality at large transverse spin in defect CFT*, *JHEP* **09** (2018) 091 [arXiv:1712.08185] [INSPIRE].

[10] P. Liendo, C. Meneghelli and V. Mitev, *Bootstrapping the half-BPS line defect*, *JHEP* **10** (2018) 077 [arXiv:1806.01862] [INSPIRE].

[11] M. Hogervorst, *Crossing Kernels for Boundary and Crosscap CFTs*, arXiv:1703.08159 [INSPIRE].

[12] F. Gliozzi, *Truncatable bootstrap equations in algebraic form and critical surface exponents*, *JHEP* **10** (2016) 037 [arXiv:1605.04175] [INSPIRE].

[13] A. Gadde, *Conformal constraints on defects*, *JHEP* **01** (2020) 038 [arXiv:1602.06354] [INSPIRE].

[14] E. Lauria, M. Meineri and E. Trevisani, *Radial coordinates for defect CFTs*, *JHEP* **11** (2018) 148 [arXiv:1712.07668] [INSPIRE].

[15] F. Gliozzi, *More constraining conformal bootstrap*, *Phys. Rev. Lett.* **111** (2013) 161602 [arXiv:1307.3111] [INSPIRE].

[16] S. El-Showk and M.F. Paulos, *Extremal bootstrapping: go with the flow*, *JHEP* **03** (2018) 148 [arXiv:1605.08087] [INSPIRE].
[17] D. Mazáč and M.F. Paulos, The analytic functional bootstrap. Part I: 1D CFTs and 2D S-matrices, JHEP 02 (2019) 162 [arXiv:1803.10233] [INSPIRE].

[18] D. Mazáč, Analytic bounds and emergence of AdS$_2$ physics from the conformal bootstrap, JHEP 04 (2017) 146 [arXiv:1611.10060] [INSPIRE].

[19] D. Mazáč and M.F. Paulos, The analytic functional bootstrap. Part II. Natural bases for the crossing equation, JHEP 02 (2019) 163 [arXiv:1811.10646] [INSPIRE].

[20] D. Mazáč, A Crossing-Symmetric OPE Inversion Formula, JHEP 06 (2019) 082 [arXiv:1812.02254] [INSPIRE].

[21] S. El-Showk and M.F. Paulos, Bootstrapping Conformal Field Theories with the Extremal Functional Method, Phys. Rev. Lett. 111 (2013) 241601 [arXiv:1211.2810] [INSPIRE].

[22] K. Sen and A. Sinha, On critical exponents without Feynman diagrams, J. Phys. A 49 (2016) 445401 [arXiv:1510.07770] [INSPIRE].

[23] R. Gopakumar, A. Kaviraj, K. Sen and A. Sinha, Conformal Bootstrap in Mellin Space, Phys. Rev. Lett. 118 (2017) 081601 [arXiv:1609.00572] [INSPIRE].

[24] R. Gopakumar, A. Kaviraj, K. Sen and A. Sinha, A Mellin space approach to the conformal bootstrap, JHEP 05 (2017) 027 [arXiv:1611.08407] [INSPIRE].

[25] K. Ghosh, Polyakov-Mellin Bootstrap for AdS loops, JHEP 02 (2020) 006 [arXiv:1811.00504] [INSPIRE].

[26] P. Dey, K. Ghosh and A. Sinha, Simplifying large spin bootstrap in Mellin space, JHEP 01 (2018) 152 [arXiv:1709.06110] [INSPIRE].

[27] R. Gopakumar and A. Sinha, On the Polyakov-Mellin bootstrap, JHEP 12 (2018) 040 [arXiv:1809.10975] [INSPIRE].

[28] A. Bissi, T. Hansen and A. Söderberg, Analytic Bootstrap for Boundary CFT, JHEP 01 (2019) 010 [arXiv:1808.08155] [INSPIRE].

[29] D. Mazáč, L. Rastelli and X. Zhou, An analytic approach to $BCFT_d$, JHEP 12 (2019) 004.

[30] M.F. Paulos, J. Penedones, J. Toledo, B.C. van Rees and P. Vieira, The S-matrix bootstrap. Part I: QFT in AdS, JHEP 11 (2017) 133 [arXiv:1607.06109] [INSPIRE].

[31] J.L. Cardy and D.C. Lewellen, Bulk and boundary operators in conformal field theory, Phys. Lett. B 259 (1991) 274 [INSPIRE].

[32] M. Hogervorst and S. Rychkov, Radial Coordinates for Conformal Blocks, Phys. Rev. D 87 (2013) 106004 [arXiv:1303.1111] [INSPIRE].

[33] L. Rastelli and X. Zhou, The Mellin Formalism for Boundary CFT$_d$, JHEP 10 (2017) 146 [arXiv:1705.05362] [INSPIRE].

[34] J. Qiao and S. Rychkov, Cut-touching linear functionals in the conformal bootstrap, JHEP 06 (2017) 076 [arXiv:1705.01357] [INSPIRE].

[35] V. Goncalves and G. Itsios, A note on defect Mellin amplitudes, arXiv:1803.06721 [INSPIRE].

[36] S. Weinberg, Six-dimensional Methods for Four-dimensional Conformal Field Theories II: Irreducible Fields, Phys. Rev. D 86 (2012) 085013 [arXiv:1209.4659] [INSPIRE].

[37] J. Penedones, High Energy Scattering in the AdS/CFT Correspondence, Ph.D. thesis, Porto University, (2007), [arXiv:0712.0802] [INSPIRE].
[38] M.F. Paulos, *Towards Feynman rules for Mellin amplitudes*, JHEP 10 (2011) 074
[arXiv:1107.1504] [INSPIRE].

[39] T.C. Lubensky and M.H. Rubin, *Critical phenomena in semi-infinite systems. I. $\epsilon$ expansion for positive extrapolation length*, Phys. Rev. B 11 (1975) 4533.

[40] H.W. Diehl and S. Dietrich, *Field-theoretical approach to static critical phenomena in semi-infinite systems*, Z. Phys. B 42 (1981) 65 [INSPIRE].

[41] J.S. Reeve and A.J. Guttman, *Critical behaviour of the $n$-vector model with a free surface*, Phys. Rev. Lett. 45 (1980) 1581.

[42] J.S. Reeve, *Renormalisation group calculation of the critical exponents of the special transition in semi-infinite systems*, Phys. Lett. A 81 (1981) 237.

[43] H.W. Diehl, *The theory of boundary critical phenomena*, Int. J. Mod. Phys. B 11 (1997) 3503 [cond-mat/9610143] [INSPIRE].

[44] H. Kleinert and V. Schulte-Frohlinde, *Exact five loop renormalization group functions of $\theta^4$ theory with $O(N)$ symmetric and cubic interactions: Critical exponents up to $\epsilon^5$*, Phys. Lett. B 342 (1995) 284 [cond-mat/9503038] [INSPIRE].

[45] L.F. Alday, *Large Spin Perturbation Theory for Conformal Field Theories*, Phys. Rev. Lett. 119 (2017) 111601 [arXiv:1611.01500] [INSPIRE].

[46] L.F. Alday, J. Henriksson and M. van Loon, *Taming the $\epsilon$-expansion with large spin perturbation theory*, JHEP 07 (2018) 131 [arXiv:1712.02314] [INSPIRE].

[47] X. Zhou, *Recursion Relations in Witten Diagrams and Conformal Partial Waves*, JHEP 05 (2019) 006 [arXiv:1812.01006] [INSPIRE].

[48] E. Lauria, M. Meineri and E. Trevisani, *Spinning operators and defects in conformal field theory*, JHEP 08 (2019) 066 [arXiv:1807.02522] [INSPIRE].

[49] P. Liendo and C. Meneghelli, *Bootstrap equations for $N = 4$ SYM with defects*, JHEP 01 (2017) 122 [arXiv:1608.05126] [INSPIRE].