Bootstrap equations for N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 4 SYM with defects

This paper focuses on the analysis of 4dN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 4 superconformal theories in the presence of a defect from the point of view of the conformal bootstrap. We will concentrate first on the case of codimension one, where the defect is a boundary that preserves half of the supersymmetry. After studying the constraints imposed by supersymmetry, we will obtain the Ward identities associated to two-point functions of 12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \frac{1}{2} $$\end{document} -BPS operators and write their solution as a superconformal block expansion. Due to a surprising connection between spacetime and R-symmetry conformal blocks, our results not only apply to 4dN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 4 superconformal theories with a boundary, but also to three more systems that have the same symmetry algebra: 4dN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 4 superconformal theories with a line defect, 3dN\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N} $$\end{document} = 4 superconformal theories with no defect, and OSP(4∗|4) superconformal quantum mechanics. The superconformal algebra implies that all these systems possess a closed subsector of operators in which the bootstrap equations become polynomial constraints on the CFT data. We derive these truncated equations and initiate the study of their solutions.


Introduction
The conformal bootstrap has seen a revival in recent years and given a plethora of results both numerical and analytical. It is a non-perturbative approach to quantum field theory in which symmetry principles play a fundamental role without relying on weak-coupling expansions. Some of the highlights that the modern reincarnation of the bootstrap has given include high precision estimates of critical exponents [1][2][3], universal behavior of CFTs at large spin [4,5], and the discovery of solvable subsectors in superconformal theories [6][7][8][9]. An important class of systems and the main focus of this work are CFTs in the presence of an extended operator or "defect". Conformal defects have numerous applications in condensed matter systems, holography, and formal aspects of superconformal field theory. From the bootstrap point of view most of the work has been on correlation functions of local operators in the presence of a defect. Apart from the rich and well understood case of 2d CFTs [10][11][12], two early papers in the higher dimensional case are [13,14] where two-point functions of CFTs near a boundary were studied. More recently, the authors of [15] presented a thorough analysis of two-point functions in the presence of a defect of any codimension.
The most interesting feature of defect CFTs is the presence of additional CFT data. In addition to the standard OPE coefficients of bulk operators and OPE coefficients of local operators on the defects, there are also bulk-to-defect couplings. 1 The latter control the convergent expansion of local operators in the bulk in terms of local operators on the defect. The increase in data is balanced by the fact that there are also more constraints. Two-point functions in the presence of the defect exemplify this feature. In figure 1 there is a pictorial representation of "crossing symmetry" for this system.
The channel on the left depicts the standard bulk OPE followed by taking the expectation value in the presence of the defect. Because the full conformal symmetry is partially broken, certain bulk local operators can have a non-zero one-point function. The CFT data for this channel are thus the standard three-point couplings times the one-point function coefficients. The channel on the right corresponds to expanding each bulk local operator in terms of local operators on the boundary. The resulting two-point function of boundary operators is then fixed by the defect conformal symmetry. This channel thus involves only the new bulk-to-defect couplings, a subset of which is given by the one-point function 1 There is also the interesting possibility of introducing an expansion of the defect itself in terms of local bulk operators as proposed in [16], and recently investigated in [17]. We will not consider such expansions in this paper.

JHEP01(2017)122
= Figure 1. Crossing symmetry for two-point functions in the presence of a defect. The channel on the left represents the standard bulk OPE between local operators. The channel on the right is the bulk-to-defect OPE in which each local operator can be written as a convergent sum of defect operators.
coefficients. It should be noticed that, given CFT data for the theory in the bulk and for the theory restricted to the boundary, the crossing equations in figure 1 do not imply by themselves that the full boundary CFT is consistent. They should be supplemented with the generalization of the case presented in figure 1, when an arbitrary local operator on the boundary is added. Note that even when the external operators are identical, the coefficients of the conformal block expansion do not exhibit any positivity property. The numerical bootstrap program was applied to defect CFTs in [18,19] using the method of [20], and in [21,22] using the method of [23].
In this work we will be interested in flat defect CFTs that also exhibit supersymmetry. In particular, N = 4 SYM with 1 2 -BPS boundary conditions. This set of boundary conditions was extensively studied by Gaiotto and Witten in [24]. 2 Even though this is an interesting system on its own, a particular motivation for our work is to present a set of bootstrap equations that might be tractable analytically. The analytical understanding of the bootstrap has improved significantly in recent years, however, much remains to be done. Two-point functions in the presence of a flat defect are a promising arena: the blocks depend on only one cross-ratio and are therefore relatively simple. The addition of supersymmetry makes the setup even more attractive. Supersymmetry is a powerful tool that gives good analytic control on certain quantities which are otherwise hard to study. If a non-trivial analytic solution to the bootstrap equations is within reach, the equations presented in this paper are a promising candidate for it.
The system we discuss possesses the three-dimensional superconformal symmetry OSP(4|4). As explained in [8,9], any three-dimensional N = 4 superconformal theory contains a closed subsector of operators whose correlators are described by a one-dimensional topological theory. It was also mentioned in [9] that such considerations extend to the case in which the three-dimensional theory lives on the boundary of a d = 4, N = 4 theory. The construction consists in truncating the system by restricting to the cohomology of a special supercharge inside osp(4|4). We will therefore refer to this subsector as the "cohomological sector". The approach championed in the original papers [6][7][8][9] is that the superconformal 2 A class of 1 2 -BPS interfaces in N = 4 SYM were constructed in [25]. In the context of integrability, there have been perturbative calculations of one-point functions in this setup, see [26][27][28][29] and references therein. Defect CFTs have also been studied from the point of view of holography, see [30][31][32][33] and references thereafter.

JHEP01(2017)122
bootstrap can be implemented as a two-step process. First, we solve for the truncated cohomological sector using analytical tools, and then we proceed to the harder task of studying non-protected quantities, maybe resorting to numerics. 3 In this paper we set the stage for this two-step program. We start by obtaining the full superconformal block expansion for two-point functions of 1 2 -BPS operators. This expansion, when restricted to the cohomological sector, implies an infinite set of polynomial equations relating bulk and boundary data. We present a preliminary study of solutions to this truncated system, leaving a more complete analysis for future work. We also explore the full bootstrap equations with the goal of eventually using modern numerical techniques to study their dynamics.
It should be pointed out that all the results we present for the codimension one 1 2 -BPS defect in a d = 4, N = 4 superconformal theory, automatically extend to the case of a codimension three 1 2 -BPS defect. The connection is spelled out in the main text and includes the superconformal blocks and the existence of a cohomological sector. Moreover, and somewhat surprisingly, we find that superconformal blocks for four-point functions in d = 3, N = 4 theories and OSP(4 * |4) superformal quantum mechanics with no defect, are essentially equal to the defect superblocks we determine. The latter superblocks have also not appeared in the literature yet, and we will therefore solve four different systems in one blow.
The outline of the paper is as follows. In section 2 we study one-and two-point functions using a novel superspace setup. For two-point functions of 1 2 -BPS operators we obtain the corresponding Ward identities, which capture in an elegant and compact form the constraints of superconformal invariance. In section 3 we present the complete solution of the Ward identities in the form of a superconformal block expansion for the correlator. In section 4 we initiate the study of solutions to the bootstrap equations concentrating mostly in the restricted cohomological sector. We conclude with section 5 and gather several technical details in the appendices.

Correlation functions in superspace
Let us start by introducing the superspace we will use to describe correlation functions in a d = 4, N = 4 superconformal theory in the presence of a flat 1 2 -BPS defect. It is particularly useful for correlators of 1 2 -BPS operators and will allow us to write the Ward identities in a very compact form. Similar superspaces have already been used in the literature to study 1 2 -BPS correlators in different superconformal setups [35][36][37].

Superspace setup
The N = 4 superconformal group in four dimensions is PSL(4|4). 4 The four-dimensional Minkowski space can be extended to a superspace with coordinates
We will now turn to the discussion of the symmetry group preserved by the presence of a 1 2 -BPS boundary, namely, the three-dimensional N = 4 superconformal group OSP(4|4).
The superalgebra osp(4|4). The even and odd parts of the superalgebra g = osp(4|4) are given by where the latter notation indicates that the fermionic generators transform in the trifundamental representation of the bosonic subalgebra g 0 . This super-algebra possess a Z 2 outer automorphism M that exchanges su(2) + with su(2) − and the fermionic generators accordingly. 5 Figure 2. On the left we have the configuration space for the spacetime coordinates and on the right for the R-symmetry coordinates. The subgroup of the conformal group unbroken in the presence of the defect is clear from the picture. On the left is the group of conformal transformations of the three-dimensional boundary SO (3,2). On the right is the product of rotations SO(3) in the space orthogonal to the line with the SL(2, R) conformal transformations on the line.

JHEP01(2017)122
Superspace coordinates and their transformations. The next step is to decompose the coordinate (2.1) of the four-dimensional superspace in boundary coordinates and distance coordinates as follows where the matrix A ψ is defined in (2.3) and the subscripts b and d stand for boundary and distance respectively. More explicitly When the fermions are zero, the geometric interpretation of these coordinates is given in figure 2.
It should be noticed that the parameter ψ, which parametrizes different embeddings of OSP(4|4) ⊂ PSL(4|4) affects only the relation between fermionic coordinates as From now on we will set ψ = 1. The general case can be easily recovered using (2.7). 6 The transformation properties of these coordinates under OSP(4|4) follow from (2.2), (2.4) and (2.7). More explicitly, super-translations and super-rotations are given by where B st = −Σ B and A ∈ GL(2|2). Special superconformal transformations mix boundary and distance coordinates and act as follows It is worth remarking that on the boundary these transformations reduce to Special superconformal transformations of the three-dimensional boundary theory, namely The element η given in (2.5) generates the superconformal inversion η•X = A 2 ψ (X −1 Σ)A −2 ψ . A remark concerning the action of R-symmetry on the y coordinates is in order. The R-symmetry preserved by the boundary conditions is o(4) su(2) ⊕ su (2). One su(2) acts linearly on the indices a, b, . . . and will be denoted by su(2) + , the second su(2) acts projectively on the R-symmetry boundary coordinate y b and will be denoted by su(2) − . This observation can be better understood by noticing that the superspace setup just described can be applied to the case of a 1 2 -BPS line defect in N = 4 as well. The only difference is that we reinterpret the coordinate y as spacetime coordinates and x as Rsymmetry coordinates, see figure 2. At the complexified level both the codimension one and codimension three cases preserve an OSP(4|4) embedded in the bulk symmetry as described by (2.4). At the level of real forms, in the codimension one case the symmetry is OSP(4|4, R) ⊂ PSU(4|2, 2), while in the codimension three case one has OSP(4 * |4) ⊂ PSU(2, 2|4). The fact that a 1 2 -BPS Wilson line in N = 4 SYM is symmetric under OSP(4 * |4) was first observed in [41], the SP (4) SO(5) ⊂ SO(6) is the symmetry that survives after fixing a non-null direction in the space of six scalars.
Superspace description of codimension two defects. Four dimensional N = 4 superconformal theories also admit codimension two superconformal defects. Let us present how this set up can be realized in superspace in the case in which the defect preserves half of the supersymmetries. As before one splits the superspace coordinate in coordinates on the surface and coordinates perpendicular to the surface X = X S + X ⊥ where The subgroup of the conformal group preserved by an 1 2 -BPS surface defect is P (SU(2|2) × SU(2|2)) ⊂ PSU(2, 2|4), see e.g. [42], corresponding to elements of the form with sdet āb cd = sdet a b c d = 1. The projective condition corresponding to the letter "P " is automatic in the representation (2.2). The action of the subgroup (2.14) on (2.13)

JHEP01(2017)122
follows from (2.2). For convenience we spell it out here. Super-translations and superrotation respectively act as Concerning special superconformal transformations, forc = 0 they take the form The transformations for c = 0 are obtained from the above by changing barred and unbarred quantities. Notice that on the surface, where x ⊥ =x ⊥ = 0, the coordinates x S ,x S transform independently. We remark that both boundary defect and surface defect preserve an su(2) ⊗ su(2) subgroup of the R-symmetry algebra su(4), but with different embeddings.

Correlation functions
We now turn our attention to correlation functions in the presence of a 1 2 -BPS boundary or interface. In the first part of this section we will classify which bulk operators can have a non-zero one-point function, and which boundary operators can have a non-zero two-point function with a 1 2 -BPS bulk operator. The second part of the section contains a derivation of the so-called superconformal Ward identities (WI) for two-point functions of 1 2 -BPS bulk operators. These identities follow from the requirement that the correlation function is free of spurious singularities in the R-symmetry coordinates, a condition that goes under the name of superspace analyticity.
Supermultiplets/superfields. Let us start by reviewing the relevant representations of superconformal algebras in three and four dimensions and set up the notation. A unitary highest weight representation χ blk of psu(2, 2|4) is usually identified using the Dynkin labels {∆, ( ,¯ ), [q, p, r]}, where ,¯ and [q, p, r] are spin and su(4) labels respectively. For the discussion of correlation functions it is convenient to reorganize these labels into a representation R blk of ps(u(2|2)⊕u(2|2)), see [36]. Similarly, a unitary highest weight representation χ bdy of osp(4|4) can be identified using the Dynkin labels {δ, s, (k + , k − )}. Where δ is the dilatation weight and we use conventions with s, k ± ∈ 1 2 Z ≥0 . For the discussion of correlation functions it is convenient to reorganize these labels into a representation R bdy of u(2|2). Notice that the bosonic subalgebra of this u(2|2) is su(2) + ⊕ su(2) Lorentz ⊕ u(1) − ⊕ u(1) D , where u(1) − corresponds to the Cartan generator of su(2) − and u(1) D corresponds to dilatations in spacetime. The list of the unitary irreducible highest weight representations is given in appendix A.1.

JHEP01(2017)122
The 1 2 -BPS operators in the bulk have Dynkin labels {p, (0, 0), [0, p, 0]} and will be denoted by W p (X). The 1 2 -BPS operators in the boundary come in two families (B, +) k and (B, −) k with Dynkin labels given by {k, 0, (0, k)} and {k, 0, (k, 0)} respectively. Their superspace description is not symmetrical. The (B, +) k supermultiplets correspond to the one-dimensional representation of the u(2|2) introduced above and are denoted by ϕ +,k (X b ). The (B, −) k supermultiplets correspond to the rank 2k totally symmetric representation of the same u(2|2), they are denoted by ϕ −,k (X b , V ), where V A is an auxiliary variable with V α fermionic and V a bosonic and ϕ −, The description of superfields just provided is actually oversimplified. What we did not specify are certain analyticity constraints on the superfields. These enforce the requirement that null states are actually zero (like the familiar ∂ µ ∂ µ ϕ = 0 for a massless scalar) and include the condition that only finite dimensional representation of R-symmetry, which is non-linearly realized in our setup, can appear. While such constraints are known in various examples, see e.g. [43], we feel that they have not been fully explored in the current literature. A complete classification is left to the future.
One-point functions of bulk operators. Since the full N = 4 superconformal symmetry of the four-dimensional theory is reduced to OSP(4|4) ⊂ PSL(4|4) by the presence of the boundary, certain bulk operators can have non-vanishing one-point function. This follows form the covariance property (2.20) that we will now analyze. Due to super-translation invariance (2.10) the one-point function can depend only on the distance coordinate X d , so we introduce the notation G O (X d ) := O(X) . For given X, one can restrict attention to g ∈ OSP(4|4) that leave X invariant, i.e. such that X = g • X = X. We call this stability group S X . It follows from (2.20) that In the second line we added the covariance properties with respect to the relevant GL(1) ⊂ GL(2|2). In general the two conditions (2.23) and (2.24) should be supplemented by the JHEP01(2017)122 requirement of superspace analyticity 7 discussed below. It is not hard to verify that the S X is isomorphic to OSP(2|2) × OSP(2|2), see appendix A.2.4. The first condition (2.23) implies that for the one-point function to be non-vanishing the representation R blk has to contain a singlet with respect to the decomposition osp(2|2)⊕osp(2|2) ⊂ ps(u(2|2)⊕u(2|2)). Such representations are classified in appendix A.3.2. The second condition (2.24) then fixes the one-point function uniquely up to a multiplicative constant. Let us illustrate this point in the simple example of the one-point function of a 1 2 -BPS operator W p (X). In this case the representation R blk is trivial and the condition (2.23) is automatically satisfied. The second condition implies that where the super-Pfaffian is defined in (A.18) and the constant a p , which is left undetermined by conformal symmetry, encodes dynamical information. Notice that the condition that p is even comes from the requirement of superspace analyticity. The case of a generic bulk operator is similar, but we will not give nor use the explicit expression for their onepoint function here. Using the criterion above together with the requirement of superspace analiticity we arrive at the following Classification: the only operators, apart from the identity, that can have non-zero onepoint function in the presence of an OSP(4|4) symmetric codimension one defect are • Long representations A ∆ [2m,2n,2m],(0,0) withm, m, n ∈ Z >0 .
Notice that the one-point function of a bulk operator is non-zero if and only if the top component of the supermultiplet, obtained by setting all the Grassmann variables to zero, has a non-zero one point function. The latter can be studied using the bosonic symmety, i.e. three-dimensional conformal and R-symmetry. The first symmetry implies that only scalar operators, i.e. representations χ blk with =¯ = 0 can have non-zero one-point function, the second that only representations [q, p, r] with q, p, r all even can have non-zero one-point function. This follows by looking at the decomposition of [p, q, r] in representations of o R (4). This condition is not sufficient, as can be seen in the case of the so-called semi-short multiplets C, see appendix A.1. It is gratifying that the same classification can be obtained independently by solving the superconformal Ward identities for bulk superconformal blocks. This is done in section 3. 7 An illustrative example is the case of the three-point function of 1 2 -BPS operators. The object where . . . 0 means that there is no defect, is superconformally covariant for any p1, p2, p3. The requirement of analyticity gives the condition

JHEP01(2017)122
Bulk-boundary two-point functions. While the one-point function of a bulk operator in the presence of a defect can be considered as the analogue of the two-point function without defect, the two-point function of one bulk and one boundary operator in the presence of a defect is the analogue of the three-point function in the CFT without defect. Generically, it is a linear combination of a finite number of superconformal structures. Given a bulk point X 1 and a point on the boundary X 2,b , it is convenient to construct the combination This combination is like a distance coordinate, see (2.7), and transforms covariantly with respect with the boundary point under special superconformal transformations. One can perform a super-translation to set X 2,b = 0 followed by a special superconformal transformation, see (2.11), to set (X 1,b , X 1,d ) = (0, X 1, 2 ). In this frame the definition (2.28) reduces to an identity. Let where The last equality, which follows from (2.20), allows to reconstruct in a simple way the bulk-boundary two-point function from G O O . We can now write equations analogous to (2.23) and (2.24) in the case of one bulk and one boundary operator: In the equation above S X 1 ,X 2,b denotes the stability group of one bulk point and one boundary point, specialized to S X 1, 2 ,0 by the choice of frame (2.29). This stability group is the subgroup of the stability group of one bulk point S X=X 1, 2 OSP(2|2) × OSP(2|2) that leaves the second boundary point X 2,b = 0 fixed. It is not hard to see that this is the diagonally embedded OSP(2|2). In summary, the elements of S X 1, 2 ,0 are super-rotations, see (2.10), such that A X 1, 2 A st = X 1, 2 . We conclude that superconformal structures for the two-point function (2.29) corresponds to osp(2|2) invariant states in the triple tensor product of the sl(2|2) representations R left blk,1 , R right blk,1 , R bdy,2 , see tables 2 and 1, regarded as reducible osp(2|2) representations. Among these, the superconformal structures that are not superspace analytic should be discarded as an allowed bulk-boundary twopoint function.
The 1 2 -BPS bulk operators correspond to R left/right blk being the trivial representation. Using the result of appendix A.3.2 and by looking at table 2 one immediately concludes that the only boundary operators that can appear in the boundary OPE of a 1 2 -BPS bulk operator are the ones given in (A.4), (A.5). The range of the R-symmetry labels (k + , k − ) is dictated by superspace analyticity. We conclude that

JHEP01(2017)122
Classification: the only boundary operators O, apart from the identity, that can have a non-zero two point function with a 1 2 -BPS bulk operator in the presence of a OSP(4|4) symmetric codimension one defect are Where we used the notation (A.4), (A.5). The precise range of representation labels for the operators appearing in the bulk-boundary OPE of a given 1 2 -BPS bulk operator W p (X) can in principle be derived by imposing the requirement of harmonic analyticity. Let us illustrate the derivation in the simplest case of 1 2 -BPS boundary operators. The general case will be analyzed in a different way using the superconformal block expansion of the two-point function of 1 2 -BPS bulk operators in section 3, see equations (3.24), (3.25). The explicit form of the simplest bulk-boundary two-point function is (2.32) The middle expression in (2.32) should be compared to (2.29), the second equality follows from the definition (2.28). The constants µ + p,k are not fixed by superconformal symmetry and encode dynamical information about the boundary conditions. Analyticity in the Rsymmetry coordinates implies that p − k ∈ 2 Z ≥0 . As already discussed, the choice of superspace is not symmetric in the descriptions of 1 2 -BPS boundary operators (B, +) and (B, −). In the latter case the relevant two point function is given by (2.33) We will not consider more general bulk-boundary two-point functions in this paper as they are not directly relevant for the study of the bootstrap equations for 1 2 -BPS operators. Nevertheless we stress that the criterion given above works in general. . See (A.1) for the notation. This classification corresponds to the supermultiplets of psu(2, 2|4) for which the inducing representations R left blk and R right blk are either atypical representations whose Young diagram is a single row of even length, or long representations of the form [0, 2n] γ . This follows from the analysis in appendix A.3.2 after a rotation by 90 degrees.

Correlation functions in the presence of a 1 2 -BPS line defect and other examples
In order to discuss the classification of boundary operators that can have a nonvanishing two-point function with a 1 2 -BPS bulk operator, it is necessary to first discuss some representation theory of OSP(4 * |4) superconformal quantum mechanics. The relevant representations are reviewed in appendix A.1. The local operators living on the line defect that can have a non-vanishing two-point function with a 1 2 -BPS bulk operator are listed in (A.8).
Correlators in d = 3, N = 4 superconformal theories. These considerations are relevant for the discussion in the next section. Using the same ideas as in section 2.2, one can determine the structure of boundary three-point functions. The stability group of three points on the boundary is OSP(2|2) where O(2) ⊂ OSP(2|2) is a subgroup of the conformal group. It follows that superconformal structures for three boundary operators are in one-to-one correspondence with OSP(2|2) invariant states in the triple tensor product R bdy,1 ⊗ R bdy,2 ⊗ R bdy,3 . It follows from this observation, consulting table 2, that the only operators that can appear in the OPE of (B, +) k 1 with (B, +) k 2 are Codimension two defect. In this case the one-point functions of bulk operators contain more than one covariant structure. Even in the simplest example of 1 2 -BPS operators one has where the perpendicular supercoordinates are defined in (2.13). The selection rule that determines which supermultiplets can have a non-zero one-point function is found by looking at the stability group for one bulk point, which in this case is SU(2|2). The class of operators that can have a non-vanishing one-point function is thus much larger than in the codimension one and three cases.

The super-displacement operator
In any defect CFT there is a distinguished boundary operator known as the displacement operator. We will denote it by D(x b ). It is associated to the breaking of the translation symmetry in the direction perpendicular to the defect. In the case of a CFT in d dimensions with a d − 1 dimensional defect, it is a scalar of dimension δ = d. The goal of this section

JHEP01(2017)122
is to determine which three-dimensional N = 4 superconformal multiplet contains the displacement operator. The main condition, apart from the fact that the supermultiplet should contain a conformal primary with the correct quantum numbers, is that it preserves supersymmetry. This is translated into the requirement that in other words, Q 3d D is a conformal descendant. This condition can be understood by recalling that The condition (2.35) is strong enough to leave only two possible super-displacement operators. They correspond to the 1 2 -BPS osp(4|4) supermultiplets (B, ±) 2 . By comparison recall that flavor currents sit in (B, ±) 1 supermultiplets. The structure of the (B, +) 2 supermultiplet can be found e.g. in equation (2.5) of [9] with r = 2. The case of (B, −) 2 is obtained by applying the mirror automorphism M defined below (2.6). The supermultiplets contain also other defect primaries corresponding to the bulk conserved currents that are broken by the defect. Recall that the bulk R-symmetry current, which is in the adjoint of su(4), decomposes as (1, 0) ⊕ (0, 1) ⊕ (1, 1) with respect to the relevant embedding of the defect R-symmetry su(2) ⊕ su(2) ⊂ su(4). The supermultiplets (B, ±) 2 contain the conformal primary associated to the broken R-symmetry which is a scalar of dimension δ = 3 and R-symmetry representation (1, 1). Finally, there are operators associated to the breaking of supersymmetry with dimension δ = 3 + 1 2 spin 1 2 and R-symmetry representation ( 1 2 , 1 2 ). As expected they are part of the (B, +) 2 and (B, −) 2 supermultiplets.

Ward identities
In the presence of a superconformal defect, the functional form of the two-point function of bulk operators is not fixed by the unbroken superconformal symmetry: there are certain combinations of coordinates that are invariant under its action. The first step is to determine such "cross ratios". The second step is to impose that the two-point function is free from certain superspace singularities. This requirement produces very powerful constraints known as superconformal Ward identities (WI). While these identities are new in this specific context, identical identities have appeared before, see e.g. [35,37] . This analogy suggests a deeper connection revealed in section 3.4.2.
Superconformal invariants made of two bulk points. The set of eigenvalues of the super-matrix

JHEP01(2017)122
It turns out that the three distinguished eigenvalues 9 of Y are convenient variables in which to write the superconformal Ward identities. A possible way to visualize the invariants above is to use superconformal transformations 10 to choose a frame in which this frame is still acted upon by GL(2|2). The only invariants under this action are the eigenvalues of the matrix X 2,d X −1 1,d . It is not hard to verify that The fact that Y is defined up to inversion is translated to the requirement of Bose symmetry when changing the first and the second point.
A useful Z 2 action. It follows from the definition of z, w 1 , w 2 that these variables are defined up to three involutions: 40) The first involution corresponds to the standard action of the symmetric group on the eigenvalues of a graded matrix. The second and third follow from the fact that the matrix Y is defined by (2.37) up to inversion. 11 There is an extra Z 2 that acts on these variables as Superconformal Ward Identities. It follows from superconformal invariance and the previous discussion that . While (2.43) transforms covariantly under superconformal transformations for any function F p 1 ,p 2 , in general it does not possess the correct analyticity properties. If one expands the eigenvalues of the supermatrix (2.37) in Grassmann coordinates, one encounters poles 9 The fact that the 4 × 4 matrix Z has only three distinguished eigenvalues can be verified by setting the fermionic coordinates to zero (this can be achieved with a superconformal transformation). 10 Explicitly, this can be achieved by performing a super-translation (2.10) followed by a special supercon- The fact that a graded-symmetric matrix B that satisfies the last equality exists is not obvious but true. 11 More concretely, z, w1, w2 are defined by the three equations

JHEP01(2017)122
when the R-symmetry cross ratios are equal to the spacetime one. The vanishing of the residue at these spurious poles translates into the equations (2.44) See [35,37] for more details. In particular, (2.44) imply that F p 1 ,p 2 (t, t, t) is a constant. These equations will play a crucial role in the derivation of superconformal blocks in the next section. It should be noticed that F p 1 ,p 2 is different from zero only when p 1 and p 2 have the same parity (both even or both odd). The structure of F p 1 ,p 2 depends only on M := min(p 1 , p 2 ) and is given by Invariance under (2.40) implies that A ab (z) = A ab (1/z) and the various A ab (z) are not independent. The number of independent functions depends on whether M is even or odd. If we rewrite M = 2n + ε with ε ∈ {0, 1} the number of independent functions is (n + 1)(n + ε + 1), which coincides with the number of R-symmetry channels as discussed in appendix A.3.1. It turns out that the general solution of the Ward identities is given in terms of two constants and n(n + ε) single variable functions. This is what we will now explain.
Solving the WI in two examples. Let us first see how the superconformal Ward identities can be solved explicitly in the two examples (p 1 , p 2 ) = (1, 1) and (p 1 , p 2 ) = (2, 2). The general solution is discussed below. The goal is to solve the equations (2.44) imposing the form (2.45). It is not hard to find that in the case (p 1 , p 2 ) = (1, 1) the WI fix the two-point function uniquely up to two constants This is not surprising as in this case the external operators correspond to N = 4 free fields. We will come back to a discussion of these constants in section 4.2.1. The solution in the (p 1 , p 2 ) = (2, 2) is more interesting. It is not hard to show that the most general solution to the WI in this case is given by where 48) Notice that, given a solution of the WI, like a superconformal block, the coefficients C ± and the function H(z) can be extracted unambiguously. In this case, all the dynamics is contained in the two constants C ± and the single-variable function H(z). It should be noticed that since the structure of F p 1 ,p 2 depends only on min(p 1 , p 2 ), the results (2.46), (2.47) can be applied to the cases of F 1,2n+1 and F 2,2n+2 with n > 0 as well.

JHEP01(2017)122
The general solution of the WI. It is not hard to convince oneself that any function of the form (2.46) subject to (2.40), can be rewritten with the help of the quantities introduced in (2.48) as 49) where min(p 1 , p 2 ) = 2n + ε with ε ∈ {0, 1}. Compare to (A. 24). Notice that the part multiplied by κ is odd under (2.42) while the remaining part is even. Using the rewriting (2.49) it is not hard to show that the general solution of the WI takes the form where C 1 , C 2 are constants and where the degree of F even/odd N in w 1 , w 2 is lowered by the presence of the factor (g 1 g 2 ) 2 . This function should be invariant under (2.40), so that it depends on 1 2 N (N − 1) functions of z individually invariant under z → z −1 . The functions G even/odd N are thus specified in terms of 1 2 N (N + 1) single variable functions, moreover they vanish when evaluated at z = w 1 = w 2 and its images under (2.40). Notice that in order to extract f n (z) and F (z, w 1 , w 2 ), it is convenient to first evaluate the quantity above at the special kinematical point w 1 = z so that g 1 = 0. We should remark that while (2.50) is the most general solution to the Ward identities, it is not clear whether more convenient parameterizations exist.

Superconformal blocks
In this section we present the superconformal blocks for the two-point function of 1 2 -BPS operators, obtained by solving the Ward identities derived in the previous section. Each superconformal block captures the contribution coming from the exchange of a particular superconformal multiplet in the relevant OPE, see figure 1. A superconformal block is therefore a sum of products of spacetime and R-symmetry blocks, whose relative coefficients are fixed by superconformal symmetry.
A straightforward way to proceed is to consider the most general linear combination of spacetime and R-symmetry blocks dictated by the operator content of the exchanged supermultiplet. The latter can be extracted, for example, using superconformal characters. Notice that out of the full content of the supermultiplets only a small fraction of spacetime and R-symmetry multiplets can appear due to the bosonic symmetry, e.g. conformal representations with non-zero spin are excluded. Once the ansatz is written, we can plug it in the WI and look for a solution for the unknown coefficients. In this section we systematically scan over all possible supermultiplets both in the boundary and bulk channels. If the WI have no solution, it indicates that the corresponding supermultiplet JHEP01(2017)122 cannot be exchanged in the given channel. On the other hand, if the supermultiplets can be exchanged, the WI can be solved and the solution is unique, up to normalization. The results of this section will confirm the superspace analysis of section 2, and further refine it by providing the range of representation labels for the exchanged operators.
An alternative way to proceed is by defining bulk and boundary superblocks as eigenfunctions of the Casimir operator, supplemented by certain boundary condition reflecting the OPE behavior of the correlator. In the boundary channel, the relevant operator is the osp(4|4) Casimir acting on one of the two points, in the bulk channel, the relevant operator is the psu(2, 2|4) Casimir acting on both points. The method we used here to determine the superblocks is to solve the Casimir equations only for the bosonic subalgebras and then impose the Ward identities. This method turned out to be sufficient in this case, but in more general situations it is not. 12 To conclude the section we will present an interesting connection between several systems that exhibit osp(4|4) symmetry. It turns out that from the superblocks of the 1 2 -BPS codimension one case, one can also obtain the superblocks relevant for the codimension three case, as well as the superblocks for four-point functions of 1 2 -BPS operators in d = 3 N = 4 and d = 1 OSP(4 * |4) theories with no defect.

Boundary channel
Let us write the boundary superblock as where w 2 ± := w 1 w ±1 2 , χ bdy is a representation of osp(4|4), (k + , k − ) is a representation of su(2) ⊕ su (2), and δ denotes the dimension of the three-dimensional operator. The precise range of the finite summation in (3.1) depends on the sp(4) ⊕ su(2) ⊕ su(2) content of the exchanged supermultiplet χ bdy and can be extracted using the characters of [47]. An efficient way to solve for the spacetime and R-symmetry blocks is to use the bosonic Casimir equation. The spacetime blocks f bdy δ (z) have already been obtained in the literature [14,18]. In four dimensions they take the form Notice that the boundary block is independent of the dimensions of the external operators p 1 , p 2 . The normalization has been chosen so that f bdy δ (z) ∼ z δ in the boundary OPE channel, corresponding to z close to zero. The R-symmetry boundary block solves two Casimir equations corresponding to the two su(2) R-symmetry factors and therefore they have a factorized form. The relevant Casimir equation for each factor iŝ 12 Two examples are N = 3 theories in 4d and N = 6 theories in 3d [45,46].

JHEP01(2017)122
as can be derived by acting with the su(2) Casimir on one of the two points in R-symmetry space. The solution with the right asymptotics is given by Above, k is a non-negative integer and we emphasize that, up to normalization, the Rsymmetry blocks h bdy k (w) are Gegenbauer polynomials with argument cos φ where w = e iφ . Their asymptotic behavior is given by h bdy k (w) ∼ w −k for w ∼ 0. Notice that only integer spin representations of su(2) appear, 13 as can be understood by looking at the branching ratios (A.23).
We will now proceed to fix the coefficients in (3.1) by solving the Ward identities (2.44). In the following we will only present the non-zero solutions. These solutions correspond to the osp(4|4) supermultiplets that can have a non-vanishing two-point function with a 1 2 -BPS bulk operator. As we will see, the results of this section fully agree with the analysis in section 2.2.

Boundary superblock for χ bdy = (B, ±)
The simplest superblocks correspond to short multiplets of the type (B, ±) being exchanged. The expansion reads Notice that we suppressed the dependence on z, w + , w − from the left-hand side. The superblocks F bdy (B,−) k are obtained from F bdy (B,+) k by the replacement w + ↔ w − .

Boundary superblock for χ bdy = (B, 1)
Next in line are the blocks for multiplets of type (B, 1), the most general ansatz consistent with the character is 13 In our conventions the dimension of the su(2) representation is 2k + 1.

JHEP01(2017)122
where we have suppressed the coordinate dependence to avoid cluttering and δ = k + + k − . By solving the Ward identities one finds . (3.8)

Boundary superblock for χ bdy = L δ
For the long blocks the solutions are a quite involved, but the procedure is the same as before. We start with the most general ansatz consistent with the content of the supermultiplet and the bosonic symmetries and fix the relative coefficients using (2.44). We have written the full solution in appendix C.1. It is interesting to note that the (B, 1) blocks presented above can be obtained as a special limit of the long block: . (3.9)

Bulk channel
Now we calculate the superblocks for the bulk channel. Recall that, unlike F p 1 ,p 2 (z, w 1 , w 2 ) in the boundary channel, here Ω − p 1 +p 2 2 F p 1 ,p 2 (z, w 1 , w 2 ) is the quantity to be expanded in bulk superblocks, where and ξ R and ξ are given in (B.9) and (B.10). The origin of this factor is clear by looking at (2.43). The general form of the superblock in the bulk OPE channel is given by where χ blk is a representation of psu(2, 2|4) for which the corresponding operator has nonzero one-point function, and R is a representation of su(4) with Dynkin label [2m, 2n, 2m], n, m,m ∈ Z ≥0 . In the following we will restrict our attention to the case in which the external operators have the same quantum numbers p 1 = p 2 so that only the case m =m

JHEP01(2017)122
is relevant. Again, the most efficient way for obtaining the bosonic spacetime f blk ∆ (z) and Rsymmetry blocks h blk R (w 1 , w 2 ) is by solving the Casimir equation. In this case the relevant operator is the Casimir of the bulk symmetry acting on both points in the two-point function. As in the boundary channel, the spacetime blocks are already given in the literature [14,18]: The normalization has been chosen so that f blk ∆ (z) ∼ (z − 1) ∆ in the bulk OPE channel that corresponds to z close to one. Similarly, the R-symmetry bulk blocks can be defined as polynomial solutions to certain second and fourth order differential equations coming from the Casimir, see (B.1) for the explicit form. The normalization is chosen so that in the bulk OPE limit w i ∼ 1 we have 13) and P m (x) are Legendre polynomials normalized as P m (1) = 1. It is also useful to ob- while they are invariant under (2.40). We found it convenient in our analysis of the Casimir equation to introduce the following variables which seem to play a special role in the mathematical literature [48]. 14 The advantage of the variables t i is that we managed to write the R = [0, 2n, 0] and R = [2m, 0, 2m] blocks in closed form, see (B.3). It turns out that this change of variables brings a pleasant surprise: the R-symmetry block h blk [2q,2k,2q] can be considered as the analytic continuation in the representation labels of the standard four-point function d = 3 bosonic block of a CFT without defects. 15 The precise correspondence is obtained by comparing the Casimir equations and is given in section 3.4.2 and appendix B. For now let us notice the relation between representation labels and cross-ratios It is worth mentioning that identical external operators in the defect block correspond to non-identical external operators in the d = 3 picture: ∆ 12 = ∆ 34 = − 1 2 where ∆ ij := ∆ i − ∆ j . It turns out that such a relation can be extended to relations between full superblocks and is explained with more detail in section 3.4.2.
As for the superblocks in the defect OPE channel, the basic idea is that the Ward identities (2.44) fix the coefficients c ∆,R (χ blk ) uniquely up to an overall normalization.

JHEP01(2017)122
We fix the normalization by c ∆ min ,R (χ blk ) = 1. If the representation χ blk is such that superconformal symmetry implies O χ blk = 0, see the classification in section 2.2, then the corresponding Ward identity has no solution. This provides a cross-check of our analysis.

Bulk superblock for
Apart from the identity superblock which is just 1, the next simplest superblock correspond to the 1 2 -BPS supermultiplet B [0,p,0] . The corresponding block is non-trivial only for p = 2n even. Looking at table I in [50] we see which Lorentz scalar conformal multiplets are in the B [0,2n,0] supermultiplets. One should further restrict to operators with zero U(1) Y charge, see below (2.5) for the definition. This leads to an ansatz for the superblock of the form Notice that we suppressed the dependence on z, w 1 , w 2 . The coefficients are fixed by the Ward identities (2.44) to be (3.18) As an example, in the case n = 1, which corresponds to the exchange of the stress-tensor supermultiplet in the bulk channel, the superblock takes the form This expression is remarkably simple. It is useful to check its properties. Firstly, it is invariant under the transformations (2.40) and picks up a sign under (2.42). Secondly, when specialized to w 1 = w 2 = z it reduces to a constant. In our normalizations this constant is one. For ∆ = 2n with n ∈ Z >0 the spacetime part of the blocks takes the form where Q n (ξ) and the quantity in square bracket are Laurent polynonials in ξ. .

JHEP01(2017)122
Cohomological sector from superconformal blocks. It is instructive to see how the cohomological sector of [8,9], which we will study with more detail in section 4, is singled out by restricting the correlation functions to certain values of the cross-ratios. Superconformal blocks possess the following remarkable properties where n, k ∈ Z >0 . This means that the complicated crossing equations of figure 1 reduce to identities involving a finite number of operators at a time. The relations (3.22) can be verified using the explicit expressions for the superblocks, but it should also be possible to prove them directly from their definition. The crucial input comes from the Ward identities (2.44). Indeed, they imply that The value of the constants corresponding to (3.22) can then be fixed by a careful analysis of the relevant OPE limit.

OPE summary
Here we give a summary of the selection rules obtained in the bulk and boundary OPE channels. For the boundary channel we have where ε ∈ {0, 1} and we define the set S n,ε = n Concerning the bulk channel, the OPE of 1 2 -BPS operators is well known, see [43,51,52], and is independent of the presence of the defect. To determine which operators contribute to the two-point function under study, one should further project to the exchanged operators that can have a non-zero one-point function. The latter operation is denote by Pr. Finally, it is convenient to split the operators appearing in the OPE into even and odd parts with respect to (2.42). In summary we obtain For simplicity, in (3.25) we restricted to the case of identical operators but the general case is immediately obtained. Summation over the values of δ and ∆ for the long representations L δ (k + ,k − ) and A ∆ R is understood.

One block to rule them all
In this section we elaborate on the relation between the polynomial R-symmetry blocks and spacetime bosonic blocks. Thanks to this analysis we will be able to write, with minimal effort, the superblocks for three additional systems: d = 4, N = 4 theories in the presence of a 1 2 -BPS line defect, 3d N = 4 superconformal theories, and OSP(4 * |4) superconformal quantum mechanics. In the latter two cases the theories have no defect, and the blocks we obtain are relevant for the expansion of four-point functions of 1 2 -BPS operators. All these systems exhibit the same osp(4|4) symmetry as 4d N = 4 theories in the presence of a flat codimension one defect.
It is clear from section 2 that the superblocks for the codimension three defects can be obtained from the codimension one case by a sort of analytic continuation in the representation labels. This is a special feature of the 4d N = 4 superconformal setup, and is therefore not surprising that the two-point function superblocks are related to each other. On the other hand, the connection to four-point function superblocks in theories without defect is unexpected. It also appears to be more general and certainly warrants further study.

Line defect superblocks
The superconformal blocks for 1 2 -BPS codimension three defects in 4d N = 4 superconformal theories can be obtained in a simple way from the codimension one case presented above. This fact is manifest by comparing the superspace setup for the two defects.
Superblock dictionary. The bulk channel superblocks are related as the defect channel superblocks are related as It was mentioned in the previous section that R-symmetry blocks in the bulk OPE channel of a two-point function in the presence of a codimension one defect are essentially analytic continuation in the representation labels of standard 3d conformal blocks, see (3.16). This observation is best understood by comparing their definition via the Casimir equations.
We will now sharpen and extend this surprising observation.
Superblock dictionary: line defect/3d N = 4 theory. Let us start by giving the dictionary between blocks where the cross-ratios are identified as (3.37) The superblocks G 3d,N =4 for the 3d N = 4 superconformal theory are relevant for the expantion of the four-point function The relations (3.34) follow from the spacetime and R-symmetry block identities together with the fact that both superblocks satisfy the same superconformal Ward identities, see (2.44) and [35,37]. Above, g 3d δ, are standard 3d blocks relevant for the expansion of a four-point function of scalar operators with ∆ 3d 12 = ∆ 3d 34 = − 1 2 where ∆ ij := ∆ i − ∆ j . Moreover, the representation label k on the left hand side is an su(2) labels, while on the right-hand side is part of the usp(4) Dynkin labels [0, k]. See appendix B for more details. It should be noticed that the first equation in (3.39) fills a gap in the literature. The conformal blocks for the line defect in the bulk channel were not determined in [15]. Now, thanks to (3.39), we can say that they are as well understood as 3d bosonic blocks.

JHEP01(2017)122
Superblock dictionary: boundary/1d OSP(4 * |4) quantum mechanics. In this case the dictionary between superblocks is given by At the level of bosonic blocks, the identity reduces to an analytic continuation in the representation labels of (3.39). In particular, the spacetime part is realized by the identification Remark. There is a relation between blocks for scalar four-point functions in d dimensions, and bulk channel blocks for scalar two-point functions in the same dimensionality d in the presence of a codimension two defect [15]. It is then natural to ask whether this statement can be extended to N = 4 superconformal blocks. An equality of blocks appears rather unlikely due to the fact that the relevant one-point functions contain multiple structures, see (2.34). However, we still expect that the superconformal WI for the twopoint function of bulk operators in the codimension two case take a similar form as the superconformal WI for the four-point function of 1 2 -BPS operators in 4d, N = 4 theories. Clarifying this point remains an interesting problem for the future.

Bootstrap equations
In this exploratory section we initiate the analysis of the dynamical constraints imposed by the bootstrap equations. As discussed in the introduction, supersymmetry allows to divide the implementation of the bootstrap in two steps. Section 4.1 concentrates on the cohomological sector while 4.2 considers the full bootstrap equations.

Microbootstrap
We will now present the cohomologically truncated bootstrap equations. 16 The space of (4.1) 16 The term microbootstrap, introduced in [53], is justified by the fact that, from the point of view of the bulk theory, it corresponds to a further reduction of the so called miniboostrap equations of [6].

JHEP01(2017)122
For each element in V p there is an associated superconformal primary for the supermultiplets B [0,p,0] . Since B [0,2,0] contains the stress tensor, which one requires to be unique, the space V 2 is one-dimensional. The space V 1 corresponds to massless representations and is expected to be part of a decoupled free theory. 17 Each V p is finite dimensional. The two-point function gives a non-degenerate pairing G : V × V → C that respects the grading in (4.1). The three-point function gives a trilinear map which is invariant under permutation of the three vectors and due to the superconformal selection rules is non-vanishing only if p 1 ∈ S p 2 ,p 3 := {|p 2 − p 3 |, |p 3 − p 3 | + 2, . . . , p 2 + p 3 } , compare to (2.25). The OPE is encoded in the bilinear map Above id corresponds to the identity operator. Notice that, as opposed to the full OPE of the 4d N = 4 superconformal theory, the cohomologically truncated OPE (4.3) contains only finitely many terms on the right hand side. The crossing equations, i.e. the requirement of associativity of the OPE, take the form As an example, let us consider the case in which all V p are one-dimensional. If we choose an orthonormal basis {e p } and define C p 1 ,p 2 ,p 3 := C(e p 1 , e p 2 , e p 3 ), the crossing equations read 18 Notice that the summation is finite. The simplest solution to these equations, which corresponds to G gauge = U(1), N = 4 SYM, is given by δ a+b−2n,c n! c! a n b n . (4.6) The relations (4.5) together with the condition above correspond to the definition of an infinite dimensional associative, commutative algebra with some extra structure. It thus makes sense to consider the subalgebra generated by the unique state in V 2 which is the stress-tensor supermultiplet. This subalgebra contains a unique element for each V 2n and the relevant bootstrap equations take the form (4.5) with p i even. The extremal three-point structures C 2n,2m,2(n+m) in this case can be fixed by general principles [54] to be C 2n,2m,2(n+n) = g n+m g n g m , g n := Γ(n + 1)Γ(n + 2 c 4d ) Γ(2 c 4d ) , where c 4d is the four-dimensional central charge, which for a Lagrangian theory is c 4d = 1 4 dim(G gauge ). Moreover, it follows from the conformal Ward identities for the stress tensor 17 Notice that it cannot really be decoupled if the stress tensor is unique. 18 Solutions to these equations were also explored in [53].
Let us now introduce the boundary data. The space of 1 2 -BPS operators in a 3d N = 4 superconformal theory splits into For each element in U ± k there is an associated superconformal primary for the supermultiplets (B, ±) k . Notice that we restricted to integer values of k as they are the only ones entering in the bulk-boundary OPE. The spaces U ± 1 correspond to flavor symmetry currents and the spaces U ± 2 to the displacement operator. We thus require that dim(U + 2 ) + dim(U − 2 ) = 1 in the case of a boundary and dim(U + 2 ) = dim(U − 2 ) = 1 in the case of an interface. The boundary two-point function gives non-degenerate pairings g ± : U ± × U ± → C that respect the grading. The one-point function of bulk operators gives a map a : V → C and the bulk boundary two-point function defines two maps µ ± : V ⊗ U ± → C. It follows from superconformal symmetry that they are non-vanishing only if p − k ∈ 2Z ≥0 . The bulk-boundary OPE coefficients µ ± : V → U ± , (4.10) can then be obtained as µ ± (u, v) = g ± (u, µ ± (v)). The bootstrap equations in figure 1 take the form a( C(v 1 , v 2 )) = g ± ( µ ± (v 1 ), µ ± (v 2 )) . In the case in which all the involved spaces are one-dimensional and choosing an orthonormal basis the equations above reduce to p∈Sp 1 ,p 2 (±1) p 1 +p 2 +p 2 C p 1 p 2 p a p = δ p 1 ,even δ p 2 ,even a p 1 a p 2 + min(p 1 ,p 2 ) where in the sum on the right hand side m ∈ {ε, ε + 2, . . . } and ε = 1 for p 1 , p 2 odd and ε = 2 for p 1 , p 2 even. One can study the equations (4.12) in the simplest example in which C p 1 p 2 p are the bulk OPE coefficients corresponding to free N = 4 Maxwell theory given above. Notice that considering this equation makes sense also when dim(U ± k ) > 1, as the image of a one-dimensional space under µ + or µ − is at most one dimensional. We also remark that the equations are invariant under µ ± pm → σ m µ ± pm where σ 2 m = 1 and that the two sets of equations map to each other under a p → (−1) The simplest solution is obtained by setting µ − p,m = 0 by hand. Then the solution of the remaining equations is given by

JHEP01(2017)122
One verifies that the displacement operator condition The result (4.13) can be confirmed by a relatively simple calculation using Wick contractions and the definition of composite operators in the free theory by normal ordering. It is interesting to observe that, after imposing µ ± p,2 = M ± a p p, it appears that (4.12) has a unique solution parametrized by (M + , M − ). It is easy to generate the solution on a computer but we could not find a closed form expression. Let us also notice that the trivial defect corresponds to a p = 0, µ ± = (±1) p 2 δ p,k , in this case M ± → ∞. In order to have a sufficient set of conditions for the mircrobootstrap one needs to add the crossing equations for correlators involving two bulks and one boundary operator. The corresponding equation takes the form (4.14) Remark. It is known that the three-point functions encoded in the map C entering (4.2), are independent of marginal deformations of the N = 4 four-dimensional theory. This is a rather unique property of 4d N = 4, which is canonically equipped with a conformal manifold of complex dimension one. The boundary data encoded in µ ± and a, on the other hand, are expected to vary under such marginal deformations. It should be possible to constraint their dependence on the marginal couplings by using superconformal perturbation theory.
Microbootstrap equations in the presence of a 1 2 -BPS line defect. The considerations above can be extended to the case of a 1 2 -BPS line defect. In this case correlation functions of 1 2 -BPS local operator can be related to certain quantities in 2d YM, see [55,56].
In these examples, it should be possible to verify the microbootstrap equations directly.

The full bootstrap
In the remaining of the section we will take a closer look at the complete bootstrap equations in two examples: the two-point function of massless representations and the two-point function of the super stress-tensor. where B [0,2,0] is the 4d N = 4 stress-tensor supermultiplet and (B, ±) 1 are 3d N = 4 multiplets corresponding to conserved currents. As in the non-supersymmetric case, see [18,21,57], the crossing equation for massless representations is solved by a finite number of blocks in both channels as

JHEP01(2017)122
In this case the full bootstrap equations are satisfied if the microbootstrap equations (4.12), that in this case read (µ ± 1,1 ) 2 = ±1 + C 112 a 2 , are. The parameter λ interpolates between the trivial interface, corresponding to λ = 0, and the boundary, i.e. the interface with an empty theory on one side, corresponding to λ = ±1. It should be possible to have a Lagrangian realization of all other values of the parameter λ as well. for the bulk channel and for the boundary channel. This time the OPE contains infinitely many operators. We will use the notation T = B [0,2,0] for the stress-tensor supermultiplet. Recall that (B, +) 2 and (B, −) 2 are the correct representations to be the super-displacement operator. The corresponding OPE coefficients µ ± 2,2 should be considered on the same footing as the central charge c 4d . In this case the microbootstrap equations read The odd part of the full crossing equations is then automatically satisfied once (4.20) holds. As in the previous example, the crossing equation (4.21) involves only a finite number of blocks. The fact that a solution of (4.21) exists relies on the identity between superblocks where Ω is defined in (3.10). Let us turn to the even part of the full crossing equation.
After rewriting bulk and boundary blocks in the form (2.47) and using (4.20) to cancel the C ± contributions from the two sides of the equation, we are left with where R := 1 2 (µ T ,(B,+) 2 ) 2 + (µ T ,(B,−) 2 ) 2 , we have defined

JHEP01(2017)122
and used the shorthand notation A ∆ = A ∆ [0,0,0] , L δ = L δ (0,0) . The reduced superblocks H blk/bdy are defined by writing F bdy and Ω 2 F blk in the form (2.47). Following this procedure one finds the following explicit expressions for the boundary channel contribution and for the bulk channel where β = 1 2 ∆ + 1. Notice that strictly above the unitarity bound we have β − 2 = and H blk B [2,0,2] are immediately obtained using the relations (3.9) and (3.21).
A toy model two-point function. We will now consider the superconformal block expansion of the toy model two-point function for T T given by compare to (2.43) and (2.47), where C ± and g are certain constants. This three parameter family of functions can be obtained by postulating that T T is a linear combination of T T and the square of the two-point function of massless fields corresponding to (4.16). By expanding this function in superblocks in the two channels we extract the CFT data in terms of the parameters C + , C − , g and The data associated to the bulk channel are and but implies a remarkable identity for the bulk superblock corresponding to the exchange of the identity operator F blk I = 1, namely Notice that this result works for any N = 4 theory in the case of a trivial interface. It would be interesting to learn a general lesson from this example on how the exchange of the identity operators in the bulk channel is reproduced upon summing infinitely many contributions from the boundary channel, along the lines of [4,5].

Conclusions
In this work we studied 4d N = 4 superconformal theories in the presence of defects from the point of view of the conformal bootstrap. We obtained the full superconformal block expansion for the two-point function of 1 2 -BPS operators in the codimension one case and thanks to the analytic continuations explained in section 3, these results also apply to 4d N = 4 superconformal theories in the presence of a line defect. Moreover, they also capture the block expansion of 3d N = 4 superconformal theories and 1d OSP(4 * |4) superconformal quantum mechanics.
In the cohomological sector we presented an infinite set of polynomial equations that relate defect and bulk data. Apart from the equations coupling bulk and boundary, one also needs to consider the truncated equations for pure bulk and pure boundary CFT data separately. While solutions to the bulk equations correspond to a commutative algebra with certain properties, the associative algebra that describes the boundary data is not commutative [9]. We observed that the cohomological bulk bootstrap equations admit a truncation to a subsector corresponding to the subalgebra generated by the super stresstensor, whose structure constants are fully determined in terms of the central charge. For this to be the case, we supplemented the equations with the knowledge of so-called extremal three-point couplings which had been obtained from rather general principles in [54].
Understanding the solutions of the truncated bootstrap equations is an important step in order to fully characterize the superconformal theory, and we leave a more thorough analysis of their solutions for the future. For example, it would very interesting to map the space of solutions to the boundary conditions studied by Gaiotto and Witten in [24]. This set of superconformal 1 2 -BPS boundary conditions in N = 4 SYM is extremely rich, as they can be engineered by gauging a flavor symmetry of a generic three-dimensional N = 4 SCFT. It should be pointed out that defect configurations that are S-dual correspond to the same solution of the bootstrap equations.
Another interesting question is to understand how deep is the connection between the four systems studied in this paper: 4d N = 4 superconformal theories with codimension one and three defects, 3d N = 4 superconformal theories and 1d OSP(4 * |4) quantum mechanics. If the solutions of the truncated equations can be mapped to each other, it will bring the kinematic relation uncovered in this paper to the realm of dynamics.

JHEP01(2017)122
Given a solution of the bootstrap equations in the cohomological sector, one would like to determine a, or many, solutions to the full bootstrap equations. This is a complicated problem, in particular due to the fact that long operators with unknown conformal dimension appear. The usual approach is to apply modern numerical techniques. The absence of positivity makes the original method of [58] unsuitable, however, there are alternative techniques in the literature. Apart from the already mentioned Gliozzi approach, a new promising method is the one developed in [59].
Another interesting venue is the study superconformal defects in superconformal theories in other dimensions. Promising candidates are codimension one and two 1 2 -BPS defects in 3d N = 4 theories [60,61] and codimension two and four defects in 6d (2, 0) theories, [62,63]. It would also be interesting to investigate whether certain bulk-boundary CFT data can be computed using localization methods. This idea has been applied successfully to extract the one-point function of the super stress tensor for the case in which the defect is a Wilson line, see [64] and references therein. For the case of a codimension one defect, the results of [65,66] may already contain the necessary information to extract some boundary CFT data.

A.1 Supermultiplets
We will now review aspects of representation theory of superconformal algebras relevant for the discussion in the main text. The representation theory of the 4d N = 4 superconformal algebra is summarized in where A * is either A 1 or A 2 depending on the value of . For representations of the 3d N = 4 superconformal algebra see table 2. We give the complete list of representations however, only a few of them will appear in the boundary OPE of a 1 2 -BPS bulk operator. We will therefore simplify the notation a bit following [47]:

Supermultiplet
Unitarity bound R bdy (0;0) δ are distinguished by the fact that the sl(2|2) inducing representation becomes a totally symmetric representation with respect to each of the inducing factors. These multiplets are referred to as conserved currents as they include higher-spin conserved currents; the s = 0 corresponds to the stress-tensor supermultiplet. The sl(2|2) central charge in our normalization is given by δ − k −σ , it is half the number of boxes in the Young tableaux.

JHEP01(2017)122
We also recall that at the unitarity bound, long representation decompose as where in the second line s > 0. These relations imply that certain boundary superblocks can be obtained as the residue of the pole in δ of the long block, see (3.9). We also remark that the decompositions (A.6) take place already at the level of inducing representation R bdy given in table 2, see [68].
Representations of superconformal quantum mechanics. The relevant representation theory of OSP(4 * |4) superconformal quantum mechanics can essentially be extracted from [69]. In the following we will describe three classes of representations that are relevant for the discussion in the main text. The d = 1 superconformal algebra osp(4 * |4) contains sl(2, R) ⊕ su(2) ⊕ usp(4) as bosonic subalgebra where sl(2, R) plays the role of conformal algebra in one dimension. The operators with smallest value ∆ 1d of the one dimensional dilatation generator form a finite dimensional irreducible representation of the R-symmetry su(2) ⊕ usp(4). We use Dynkin labels to characterize such representation {n, One should also notice that T 2 R = η R , T R T † R = 1, det(T R ) = −1.

A.2.2 Conformal factors
The conformal factors are given by (1 + X C) −1 , (1 + C X) , g = 1 0 The form for general g follows by group composition law.

A.2.3 Definitions
We recall the definitions of superdeterminant and super-Pfaffian

A.3.1 R-symmetry channels
For the boundary channel we use the su(4) → su(2) ⊕ su(2) branching ratios where Pr denotes a projection into representations of the form [2a, 2b, 2c]. Notice that the sum (A.24) splits into two sums, the blocks corresponding to each sum have definite Z 2 transformation properties. The number of multiplet is the same in bulk and boundary channel and it is equal to ( + 1)( + ε + 1). We also collect the defect channel R-symmetry OPE relevant for the line defect. In this case one has su(4) → sp (4) Table 3. Brief summary of the osp(2|2) representations appearing in the branching ratios considered here, see [70] for a complete description of these representations. Above 2j ∈ Z >0 and b = ±j.
A.3.2 sl(2|2) → osp(2|2) branching ratios In this appendix we determine which sl(2|2) representation listed in tables 2 and 1 contains a state (vector) which is invariant under osp(2|2) ⊂ sl(2|2) where the choice of embedding is specified in appendix A.2.4. It turns out that, apart from the case in which the sl(2|2) representation is trivial, this invariant state will not appear as the trivial representation of osp(2|2), but as part of the so called projective cover of the trivial representation denoted by P(0), see [70]. We refer to table 3 for the notation and a few facts about the relevant representations of osp(2|2) sl(2|1). Since we could not find the branching ratios in the literature we present them below. The branching rules can be determined from the knowledge of tensor products in sl(2|2) and osp(2|2) together with the fact that the operations of taking tensor product and branching ratios commute and the tensor product is distributive with respect to direct sums. Let us illustrate this procedure in a simple example. The fundamental representation of sl(2|2) decomposes as → {0, 1 2 }. This branching ratio specifies the embedding osp(2|2) ⊂ sl(2|2). Taking tensor products and branching ratios of the fundamental representation with itself gives  table 3 for details about the representations on the right hand side. Iterating this procedure, with extensive use of osp(2|2) sl(2|1) tensor product decompositions from [70], one arrives at the following rules (A.31) Some remarks on the notation are in order: the sums are integer spaced between the limits of summation, means that one should exclude from the sum terms that give contributions of the form {±j, j}, the latter are exactly the one that recombines into projective representations P, finally g ∈ Z ≥0 .
The only osp(2|2) representation appearing in the right hand side of the list above that contains an invariant state is P(0). Notice that where ∼ means that they have the same character. In both cases the trivial representation {0} appears. In the first case it appears twice, only one of them is a true osp(2|2) invariant state. In the second case the state associated to {0} by the character is not an invariant, see [70]. This important difference would have been missed by a character analysis. Looking at the branching ratios given above we conclude that, apart from the trivial representation of sl(2|2), only the totally antisymmetric representations in 2g+2 indices and the long representations [a, s] γ = [2g, 0] γ with g ∈ Z ≥0 contain an osp(2|2) invariant state.

B Comparison of conformal blocks
In this short appendix we present the details of the comparison between our conformal blocks and the results already present in the literature.  w i (w i − 1) 2 ∂ 2 w i + k(w 1 , w 2 )∂ w 1 + k(w 2 , w 1 )∂ w 2 h blk R = C R h blk R , (B.1)