# A functional approach to the numerical conformal bootstrap

Bernardo Zan, Miguel F. Paulos

► To cite this version:

Bernardo Zan, Miguel F. Paulos. A functional approach to the numerical conformal bootstrap. Journal of High Energy Physics, 2020, 2020 (09), pp.006. 10.1007/JHEP09(2020)006. hal-02410945

HAL Id: hal-02410945

https://hal.science/hal-02410945v1

Submitted on 21 Aug 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A functional approach to the numerical conformal bootstrap

Miguel F. Paulos<sup>a</sup> and Bernardo Zan<sup>b,a</sup>

E-mail: miguel.paulos@lpt.ens.fr, bernardo.zan@epfl.ch

**ABSTRACT:** We apply recently constructed functional bases to the numerical conformal bootstrap for 1D CFTs. We argue and show that numerical results in this basis converge much faster than the traditional derivative basis. In particular, truncations of the crossing equation with even a handful of components can lead to extremely accurate results, in opposition to hundreds of components in the usual approach. We explain how this is a consequence of the functional basis correctly capturing the asymptotics of boundsaturating extremal solutions to crossing. We discuss how these methods can and should be implemented in higher dimensional applications.

**KEYWORDS:** Conformal Field Theory, Field Theories in Lower Dimensions, Nonperturbative Effects

ARXIV ePRINT: 1904.03193

<footnote><sup>a</sup>Laboratoire de Physique de l'École Normale Supérieure, PSL University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, 24 rue Lhomond, 75231 Paris Cedex 05, France</footnote>

<footnote><sup>b</sup>Theoretical Particle Physics Laboratory (LPTP), Institute of Physics, EPFL, Lausanne, Switzerland</footnote>
# Contents

1 Introduction 1

2 Setup 3

2.1 Crossing symmetry 3

2.2 Free functional bases 4

3 Bounds and the choice of basis 6

3.1 Bounds and extremality 7

3.2 Arguments for the new bases 9

4 Applications 13

4.1 Bootstrapping bosons with fermions 14

4.2 OPE maximization, or $\phi^4$ flow 17

5 Outlook — higher dimensions 20

A Functional kernels 23

B Computation of the functional action for special cases 24

B.1 The functional kernels 24

B.2 Computation of beta functionals 25

B.3 Computation of the alpha functionals 26

## 1 Introduction

More than ten years ago, [1] introduced a simple and effective method to extract information from crossing equations in conformal field theories. In its simplest guise,<fnref n="1" /> the equation takes the form

$$ \sum_{O \in \phi \times \phi} \lambda_{\phi\phi O}^2 F_{\Delta_O}(z|\Delta_{\phi}) = 0., \qquad (1.1) $$

with kinematically determined functions $F_\Delta$ given below. This equation should be thought of as expressing an infinite set of sum rules, labeled by the cross-ratio $z$, on the OPE data $\lambda_{\phi\phi O}$ between primary operators $\phi, \phi, O$. The basic idea of [1] is as follows.<fnref n="2" /> Firstly we abstract from any particular CFT realization and think of the above as a general constraint

<footnote n="1">¹In this example the equation expresses crossing symmetry for a correlator of four identical fields $\phi$ on a line, when decomposed into SL(2, $\mathbb{R}$) blocks.</footnote>

<footnote n="2">²For a recent review and an extensive list of references see [2].</footnote>
on all possible consistent OPEs. Secondly we act on the crossing equation with suitably chosen linear functionals $\omega$ to get simpler sum rules:

$$ \sum_{\Delta \in S} a_{\Delta} \omega(\Delta | \Delta_{\phi}) = 0, \quad \omega[F_{\Delta}(z | \Delta_{\phi})] = \omega(\Delta | \Delta_{\phi}), \qquad (1.2) $$

where for unitary solutions we have $a_{\Delta_O} \equiv \lambda_{\phi\phi O}^2 \ge 0$. We can now make various assumptions on the set $S$ of allowed quantum numbers, on the external dimension $\Delta_{\phi}$, or even the allowed values of $a_{\Delta}$ for specific $\Delta$. Depending on these assumptions, it may be possible to find $\omega(\Delta|\Delta_{\phi})$ non-negative in $S$. The sum rule is then impossible to satisfy for positive $a_{\Delta}$, which means our original assumptions cannot hold for any unitary CFT. In this way, it is possible to derive bounds on CFT data such as scaling dimensions and OPE coefficients.

In practice, the search for such functionals is done numerically inside a finite dimensional linear vector space. This amounts to truncating the crossing constraints to a more manageable finite subset. We want the size of this space to be as large as possible to get the strongest bounds, but we are constrained in doing this by the available computational resources. Clearly we would like to choose this space in a way which gives us the most stringent constraints possible for a fixed dimensionality. In the original reference [1] this choice corresponds to taking $\omega$ to be a finite sum of derivatives with respect to the cross-ratio $z$ at a specific point. For its simplicity and effectiveness, this has remained the default choice up to now.<fnref n="3" />

Very recently, new bases of functionals were constructed in [4] (see also [5–9] for related works) and found to lead to exact, optimal, bounds on CFT data. The goal of this note is to explore how these bases can be used to obtain better numerical bounds. As we will see, one of the key points is that in these bases each individual functional automatically encodes the expected asymptotics and spectrum of bound-saturating (extremal [10, 11]) solutions to crossing. This leads to a decoupling between dynamics at small and large scaling dimensions, which means the truncated bootstrap equations not only fully capture the interesting physics, they may also be effectively completed to the full set of constraints once the truncation size is large enough.

In this note we will consider two simple applications to test our proposal. In the first we will bootstrap the generalized free boson solution using a functional basis associated to free fermions. This somewhat perverse choice was chosen as a worse-case scenario for our approach, and yet we will still find that it is a significant improvement over the derivative basis. For our second application we will find a family of universal bounds on OPE coefficients in 1D. The bound is saturated by a family of solutions to crossing which neatly interpolates between the generalized free boson and fermionic solutions. Near the boson point it is described by (boundary correlators of) a free scalar in AdS$_2$ with a $\phi^4$ interaction. We are able to easily find the approximate exact bound, obtaining several digits accuracy using a handful of components.

The outline of this note is as follows. In the next section we review the functional bases obtained in [4]. Section 3 reviews the link between bounds and extremal solutions,

<footnote n="3">³See [3] for some alternative explorations.</footnote>
and explains why the new basis is expected to lead to highly improved numerical bounds. In section 4 we do the two numerical applications explained above, which confirm these expectations. We conclude with an outlook on higher dimensional applications.

This note is complemented by two appendices. The reader interested in using the functionals will find the explicit form of the associated integral kernels in appendix A. Furthermore for integer or half-integer external dimensions, the kernels simplify and we explain how to compute the functional actions explicitly in appendix B.

# 2 Setup

## 2.1 Crossing symmetry

We are interested in studying the crossing equation for 1D CFT correlators of identical operators $\phi$:

$$ \langle \phi(x_1)\phi(x_2)\phi(x_3)\phi(x_4) \rangle = \frac{\mathcal{G}(z)}{x_{13}^{2\Delta_{\phi}} x_{24}^{2\Delta_{\phi}}}, \quad x_{ij} = x_i - x_j, \quad z = \frac{x_{12}x_{34}}{x_{13}x_{24}}. \quad (2.1) $$

For $z \in (0, 1)$ the correlator can be expressed in terms of SL(2, $\mathbb{R}$) blocks as,

$$ \mathcal{G}(z) = \sum_{O \in \phi \times \phi} a_{\Delta_O} \frac{G_{\Delta_O}(z)}{z^{2\Delta_O}}, \quad a_{\Delta_O} = \lambda_{\phi\phi O}^2, \quad (2.2) $$

with the conformal block $G_{\Delta}(z)$ given by

$$ G_{\Delta}(z) := z^{\Delta} {}_2F_1(\Delta, \Delta, 2\Delta, z). \quad (2.3) $$

Crossing symmetry is the statement $\mathcal{G}(z) = \mathcal{G}(1-z)$. For a generic CFT this leads to the equation

$$ \sum_{\Delta} a_{\Delta} F_{\Delta}(z|\Delta_{\phi}) = 0, \quad (2.4) $$

with $F_{\Delta}(z|\Delta_{\phi}) = z^{-2\Delta_{\phi}}G_{\Delta}(z) - (z \leftrightarrow 1-z)$. We will often simplify notation by omitting the dependence on $\Delta_{\phi}$. For unitary CFTs the sum ranges over all scaling dimensions $\Delta \ge 0$. We will study this equation as a mathematical set of constraints on the $a_{\Delta}$ which hold for any and all (unitary) CFTs.

Two simple and important solutions to these constraints are associated with generalized free fields in $D=1$, which describe boundary correlators of free fields in AdS$_2$. Introduce:

$$ a_{\Delta}^{\text{free}} = \frac{2 \Gamma(\Delta)^2 \Gamma(\Delta + 2\Delta_{\phi} - 1)}{\Gamma(2\Delta_{\phi})^2 \Gamma(2\Delta - 1) \Gamma(\Delta - 2\Delta_{\phi} + 1)}. \quad (2.5) $$

Then we have

$$ F_0(z|\Delta_{\phi}) + \sum_{n=0}^{+\infty} a_{\Delta_n^{B,F}}^{\text{free}} F_{\Delta_n^{B,F}}(z|\Delta_{\phi}) = 0 \quad (2.6) $$

with $\Delta_n^B = 2\Delta_{\phi} + 2n$ and $\Delta_n^F = 1 + 2\Delta_{\phi} + 2n$. For $z \in (0, 1)$ the corresponding correlators are

$$ \mathcal{G}^{B,F}(z) = \pm 1 + z^{-2\Delta_{\phi}} + (1-z)^{-2\Delta_{\phi}} \quad (2.7) $$

with the + and - signs for (B)osons and (F)ermions respectively.
## 2.2 Free functional bases

We will now review the analytic functionals constructed in [4]. There are two functional bases which are in a sense dual to the simple solutions to crossing of the previous subsection.<fnref n="4" /> We will denote a generic functional by $\omega$ and functional actions will be abbreviated as

$$ \omega(\Delta) \equiv \omega[F_{\Delta}(z|\Delta_{\phi})]. \qquad (2.8) $$

Each basis is made up of infinite sets of functionals $\alpha_n, \beta_n$ with $n$ a non-negative integer, and satisfying the duality conditions

$$ \begin{aligned} \alpha_n^F(\Delta_m^F) &= \delta_{nm}, & \partial_\Delta \alpha_n^F(\Delta_m^F) &= 0, \\ \beta_n^F(\Delta_m^F) &= 0, & \partial_\Delta \beta_n^F(\Delta_m^F) &= \delta_{nm}, \end{aligned} \qquad (2.9) $$

for the fermionic (F) basis, and

$$ \begin{aligned} \alpha_n^B(\Delta_m^B) &= \delta_{nm}, & \partial_\Delta \alpha_n^B(\Delta_m^B) &= -\delta_{0m}c_n, \\ \beta_n^B(\Delta_m^B) &= 0, & \partial_\Delta \beta_n^B(\Delta_m^B) &= \delta_{nm} - \delta_{0m}d_n, \end{aligned} \qquad (2.10) $$

for the bosonic (B) basis, with $\beta_0^B \equiv 0$. The coefficients appearing here are given by

$$ c_n = \frac{1}{2} \partial_n d_n, \qquad d_n = \frac{(\Delta_{\phi})_n^4 (4\Delta_{\phi} - 1)_{2n}}{(n!)^2 (2\Delta_{\phi})_n^2 (4\Delta_{\phi} + 2n - 1)_{2n}}. \qquad (2.11) $$

The fact that $c_n, d_n$ are non-zero is related to the fact the generalized free boson solution admits a relevant deformation that does not introduce new states. Thinking of the free boson as being the holographic dual to a massive scalar field $\phi$ in AdS$_2$, this deformation is accounted for by the $\phi^4$ interaction [4].

The action of the functionals is defined by the formula:

$$ \omega[F_{\Delta}] = \frac{1}{2} \int_{\frac{1}{2}}^{\frac{1}{2}+i\infty} dz f_{\omega}(z) F_{\Delta}(z) + \int_{\frac{1}{2}}^{1} dz g_{\omega}(z) F_{\Delta}(z) \qquad (2.12) $$

with $g_{\omega}(z) = \pm(1-z)^{2\Delta_{\phi}-2}f_{\omega}(\frac{1}{1-z})$ and the $+(-)$ sign for boson (fermion) functionals respectively. The general functional kernels $f_{\omega}$ are given very explicitly in appendix A. They are chosen such that not only the duality relations hold, but also that the functionals $\omega$ are compatible with the crossing equation, in the sense that:

$$ F_0 + \sum_{\Delta} a_{\Delta} F_{\Delta} = 0 \quad \Rightarrow \quad \omega \left[ F_0 + \sum_{\Delta} a_{\Delta} F_{\Delta} \right] = \omega(0) + \sum_{\Delta} a_{\Delta} \omega(\Delta). \qquad (2.13) $$

The infinite sum appearing in the crossing equation means this swapping condition is not trivial, and it in fact constrains the behaviour of $f(z)$ near $z = 0, 1$ and $\infty$ [12].

<footnote n="4">⁴Note we use the word “basis” loosely — what we mean here is that the functionals fully encode the constraints implied by the $z = \bar{z}$ SL(2, $\mathbb{R}$) crossing equation under the assumption of Regge boundedness.</footnote>
One of the main results of [4] is that the crossing equation is equivalent to the set of sum rules arising from applying a complete functional basis:

$$ \sum_{\Delta} a_{\Delta} F_{\Delta}(z) = 0 \quad \text{for all} \quad z \in (0, 1) $$

is equivalent to

$$ \sum_{\Delta} a_{\Delta} \omega(\Delta) = 0 \quad \text{for all} \quad \omega \in \{\alpha_n, \beta_n, \ n \in \mathbb{N}_{\ge 0}\}. $$

The equivalence holds whether we use the bosonic or the fermionic basis of functionals. Because of this equivalence between the crossing equation and the sum rules involving the functional actions $\alpha_n(\Delta)$, $\beta_n(\Delta)$ it is clearly of utmost importance to obtain expressions for these, or at least to have the means for their speedy numerical evaluation. In the derivative basis, this would be analogous to being able to compute arbitrary derivatives of conformal blocks efficiently. In our case, a useful result is that for each basis, for $\Delta > \Delta_m$ a contour rotation argument leads to the expression

$$ \omega_m(\Delta) = [1 \pm \cos \pi(\Delta - 2\Delta_{\phi})] \int_0^1 dz g_{\omega_m}(z) \frac{G_{\Delta}(z)}{z^{2\Delta_{\phi}}} \quad (2.14) $$

where again the $+(-)$ sign corresponds to the bosonic (fermionic) functionals. This expression is essentially what lies at the origin of the duality conditions, thanks to the oscillating prefactor. The limited validity of this formula comes about due to singularities in $g(z)$ as $z \to 0$. However, if one manages to evaluate the integral above exactly, we can of course extend the result to all $\Delta$ by analytic continuation.

With the explicit forms of the $f_{\omega}(z)$ given in appendix A, it is relatively straightforward to evaluate numerically the functional actions for any $\Delta$, using either (2.12) or (2.14), and this is all we need for bootstrap applications.<fnref n="5" /> However, the kernels take on particularly simple expressions when $\Delta_{\phi}$ is a half-integer, for the fermion basis, or an integer, for the bosonic basis. With these simplified kernels we are able compute the functional actions analytically, and this is discussed in appendix B. These analytic expressions can be evaluated much more rapidly than by doing the integrals numerically. It would be very useful to be able to determine similar analytic results for any $\Delta_{\phi}$, although this is strictly speaking not necessary for numerical applications.<fnref n="6" />

**Examples.** Since the discussion above is perhaps unfamiliar, it may be useful to consider a couple of concrete examples. For general $\Delta_{\phi}$, $f(z)$ are sums of $_3F_2$ hypergeometric functions. In the special case $\Delta_{\phi} = 1$ we have instead, for instance, the simple expression:

$$ f_{\alpha_0^B}(z) = \frac{2 \left( (z-1)z + (z-2)(z+1)(2z-1) \coth^{-1}(1-2z) + 1 \right)}{\pi^2 (z-1)z} \quad (2.15) $$

The functional action can now easily be computed numerically by plugging in this expression into (2.12) with $g_{\alpha_0^B}(z) = (1-z)^{2\Delta_{\phi}-2} f_{\alpha_0^B} \left( \frac{1}{1-z} \right)$. Equivalently, we can use the

<footnote n="5">⁵In our experience, it is faster to numerically evaluate (2.14) when available, rather than (2.12).</footnote>

<footnote n="6">⁶Along these lines the results of [13] seem to be very promising at least for the bosonic functional actions.</footnote>
methods described in appendix B to get an analytic expression. Either way the motivated reader can easily check that $\alpha_0^B(0) = -2$, $\alpha_0^B(2+2n) = \delta_{n,0}$, with the general shape of $\alpha_0^B(\Delta)$ as in figure 1. This is perfectly consistent with the fact that for $\Delta_\phi = 1$ the generalized free boson solution takes the form:

$$F_0(z) + 2F_2(z) + \sum_{n \ge 1}^{\infty} a_n F_{2+2n}(z) = 0. \qquad (2.16)$$

Similarly, for $\Delta_\phi = 1/2$ the fermionic basis also simplifies. We have for instance:

$$g_{\beta_0^F}(z) = \frac{(1-z)(2z^2+z+2)}{\pi^2 z^2}. \qquad (2.17)$$

Instead of using (2.12), in this case we will do the integral (2.14) directly, finding the analytic result<fnref n="7" />

$$\begin{aligned} \beta_0^F(\Delta) &= \frac{2 \sin^2\left(\frac{\pi \Delta}{2}\right)}{\pi^2} \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2} \left[ \frac{1}{(\Delta-2)(\Delta+1)} \right. \\ & \qquad \left. + \Gamma(\Delta)^3 \left( {}_3\tilde{F}_2(\Delta, \Delta, \Delta; 2\Delta, \Delta+1; 1) - 2\Delta {}_3\tilde{F}_2(\Delta, \Delta, \Delta+1; 2\Delta, \Delta+2; 1) \right) \right] = \\ &= \frac{2 \sin^2\left(\frac{\pi \Delta}{2}\right)}{\pi^2} \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2} \left\{ \frac{1}{(\Delta-2)(\Delta+1)} + \right. \\ & \qquad \left. + \left(\Delta(\Delta-1) + \frac{1}{2}\right) \left[ \psi^{(1)}\left(\frac{\Delta}{2}\right) - \psi^{(1)}\left(\frac{\Delta+1}{2}\right) \right] - 2 \right\} \quad (2.18) \end{aligned}$$

where ${}_3\tilde{F}_2$ stands for the regularized hypergeometric function. One can now explicitly check that the properties (2.9) are satisfied: the term in parenthesis is finite for every value of $\Delta = \Delta_m^F = 2+2m$ for $m > 0$, while it has a single pole for $\Delta = \Delta_0^F = 2$. However, the prefactor has double zeroes at $\Delta = \Delta_m^F = 2+2m$ for every $m \ge 0$, and we have indeed that $\beta_n^F(\Delta_m^F) = 0$.

### 3 Bounds and the choice of basis

The goal of this section is to review how bounds on CFT data are associated to extremal functionals and sparse solutions to crossing, and how this perspective explains why our proposed bases are ideally suited for numerics. In the next subsection we remind the reader of the basic idea of [1], and the link between optimal bounds and extremal solutions to the crossing equation [10, 11]. In subsection 3.2 we will argue for our new basis of functionals. To do this we will use the link between extremal solutions and bounds, studying the former to understand the latter. In particular we will show that the functional bases described in the previous section strongly constrain the form of such solutions, leading to fast convergence of numerical bootstrap computations.

<footnote n="7">⁷This particular functional and corresponding functional action first appeared in [6].</footnote>
## 3.1 Bounds and extremality

Consider a set of solutions to crossing of the form

$$F_0(z) + a_0 F_{\Delta_0} + \sum_{\Delta \ge \Delta_g} a_\Delta F_\Delta(z) = 0. \quad (3.1)$$

We would like to determine an upper bound on $a_0$ [14]. To do this we apply a linear functional $\omega$ to this equation, setting without loss of generality $\omega(\Delta_0) = 1$. We will consider $\omega$ to lie inside some linear vector space $W_N$ of finite dimension. We get

$$a_0 = -\omega(0) - \sum_{\Delta \ge \Delta_g} a_\Delta \omega(\Delta). \quad (3.2)$$

Hence for any solution to the crossing equation (3.1) we have the bound

$$a_0 \le -\omega^{\text{opt},N}(0) \equiv \min_{\omega \in W_N} -\omega(0) \quad \text{such that} \quad \omega(\Delta) \ge 0 \quad \text{for all} \quad \Delta \ge \Delta_g. \quad (3.3)$$

If $\hat{\omega}_n$ are a basis for $W_N$, then we can define the dual space to be:

$$W_N^* = \operatorname{span} \{\vec{F}_\Delta(z), \Delta \ge 0\}, \quad \vec{F}_\Delta := (\hat{\omega}_1(\Delta), \dots, \hat{\omega}_N(\Delta)). \quad (3.4)$$

Then we can also determine a bound by solving a problem in this space:

$$a_0 \le a_0^{\text{opt},N} \equiv \max a_0 \quad \text{such that} \quad \vec{F}_0 + a_0 \vec{F}_{\Delta_0} + \sum_{\Delta \ge \Delta_g} a_\Delta \vec{F}_\Delta = 0, \quad (3.5)$$

for some $a_\Delta \ge 0$, since increasing $N$ adds more constraints and hence lowers $a_0$. These kinds of problems are known as semi-infinite linear programs [15], and the theory of linear convex optimisation gives us an important result:

$$-\omega^{\text{opt},N}(0) = a_0^{\text{opt},N}. \quad (3.6)$$

In other words the functional problem and the dual problem are equivalent.

This deserves some explanations. In the dual space formulation, we always find an explicit solution to the truncated form of the crossing equation (i.e. that involving the $\vec{F}_\Delta$), since we are satisfying the constraints throughout the maximization process. Meanwhile in the original formulation of the problem one constructs functionals which are supposed to rule out such solutions. How can these problems then be equivalent? The reason is that for any fixed $N$, at optimality there is indeed a solution to crossing but a very special one. The statement is that there exists a (generically unique) set of $a_n \ge 0$, $\Delta_n$ with $n = 1, \dots, N-1$ such that not only

$$\vec{F}_0 + a_0^{\text{opt},N} \vec{F}_{\Delta_0} + \sum_{n=1}^{N-1} a_n \vec{F}_{\Delta_n} = 0, \quad (3.7)$$

but furthermore

$$\omega^{\text{opt},N}(\Delta_n) = 0, \quad \partial_\Delta \omega^{\text{opt},N}(\Delta_n) = \delta_{\Delta_n, \Delta_g} \quad \text{if} \quad a_n > 0. \quad (3.8)$$
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Schematic form of the optimal functional. Such a functional provides a valid upper bound on the OPE coefficient at $\Delta_0$ for any choice of $\Delta_g$ down to where the functional first becomes negative.</figcaption></figure>

Note that these conditions are clearly necessary to make sure that $\omega(\Delta)$ is positive for $\Delta \ge \Delta_g$, see figure 1. The conditions (3.7) together (3.8) are the so-called Karush-Kuhn-Tucker second order optimality conditions. They are necessary conditions for optimality, with sufficiency being guaranteed as long as the functional $\omega$ is positive above $\Delta_g$.

For the kinds of problems we are interested in the conformal bootstrap, the fact that the set of allowed $\Delta$ is continuous typically leads to very sparse (extremal) solutions. For the 1D bootstrap problems considered in this note and elsewhere in the literature, one always finds (experimentally) that only $\lceil \frac{N-1}{2} \rceil$ operators appear in the extremal solution (beyond the identity and that at $\Delta_0$).<fnref n="8" />

Thus, extremal functionals go hand in hand with extremal solutions to crossing. What is perhaps less obvious is that at extremality we get not only the functional we were looking for, but actually a full basis for $W_N$. To see this very explicitly, let us consider the case where none of the $\Delta_n$ sit at the gap, set $N$ to be an odd number and consider the case where the extremal solution has only $(N-1)/2$ non-zero $a_n$. In this case the $N$ constraints from the crossing equation neatly match up with the unknown $(N+1)/2$ OPE coefficients (i.e the previously mentioned $(N-1)/2$ non-zero $a_n$ and $a_0$) and $(N-1)/2$ scaling dimensions of the extremal solution. We can then choose:

$$W_N^* = \operatorname{span} \left\{ \vec{F}_{\Delta_0}; \vec{F}_{\Delta_n}, \partial_{\Delta} \vec{F}_{\Delta_n}, \quad n = 1, \dots, \frac{N-1}{2} \right\}. \qquad (3.9)$$

We now construct a basis for the space of functionals $W_N$ by setting

$$W_N = \operatorname{span} \left\{ \alpha_0; \alpha_n, \beta_n, \quad n = 1, \dots, \frac{N-1}{2} \right\} \qquad (3.10)$$

with functionals $\alpha, \beta$: satisfying the duality conditions:

$$\begin{align}
\alpha_n(\Delta_m) &= \delta_{nm}, & \partial_\Delta \alpha_n(\Delta_m) &= -\delta_{0m}C_n, \tag{3.11} \\
\beta_n(\Delta_m) &= 0, & \partial_\Delta \beta_n(\Delta_m) &= \delta_{nm} - \delta_{0m}D_n. \nonumber
\end{align}$$

<footnote n="8">⁸A similar statement apparently holds for the modular bootstrap when ignoring spin [16]. The most general situation, which holds for instance in higher dimensions, is analysed in [10].</footnote>
with $\beta_0 \equiv 0$ and some constants $C_n, D_n$. These constants, as well as the functionals, can be obtained completely explicitly by putting the basis vectors of $W_N$ into a matrix and computing its inverse [10], leading to expressions in terms of the original set of $\hat{\omega}_i$. In particular, we see that our original OPE maximization functional $\omega^{\text{opt},N}$ is nothing but $\alpha_0$, as it satisfies (3.8). Essentially the $\alpha$ functionals provide OPE bounds while $\beta$ functionals provide bounds on dimensions of operators, subject to certain assumptions. These bounds must exist, as they encode the uniqueness of the extremal solution under such assumptions.

To summarize, the message of this section is that optimal bounds are associated to extremely simple solutions to crossing containing as few operators as possible, which saturate these bounds. Such solutions in turn determine a pair of bases for $W$ and $W^*$. Hence, bounds on the CFT data follow from the study of these extremal solutions and their associated functionals.<fnref n="9" />

## 3.2 Arguments for the new bases

From now on we shift our perspective from determining bounds to constructing extremal solutions to the truncated crossing equations:

$$ 
\begin{align} 
& \vec{F}_0 + a_0 \vec{F}_{\Delta_0} + \sum_{m=1}^{\frac{N-1}{2}} a_m \vec{F}_{\Delta_m} = 0 \tag{3.12} \\ 
& \Leftrightarrow \hat{\omega}_n(0) + a_0 \hat{\omega}_n(\Delta_0) + \sum_{m=1}^{\frac{N-1}{2}} a_m \hat{\omega}_n(\Delta_m) = 0, \quad n = 1, \dots, N. \nonumber 
\end{align} 
$$

This is because as we have just reviewed, understanding and constructing these solutions automatically leads to the desired bounds.

We are particularly interested in how the truncated solution approaches the exact one. As we dial up $N$, the number of operators in the solution increases, and the CFT data of those operators already in the solution changes, eventually approaching their exact values. Depending on the choice of basis elements $\hat{\omega}_1, \dots, \hat{\omega}_N$, the rate at which we approach the exact extremal solution, and with it the exact bounds, can be very different as we shall see.

The original choice of [1] set $\hat{\omega}_n = \partial_z^{2n-1}|_{z=\frac{1}{2}}$, and this has been the default approach ever since. We will argue that a better choice is to use the functionals $\alpha_n^{B,F}, \beta_n^{B,F}$. Whether bosonic or fermionic functionals should be used depends on the expected asymptotics of the extremal solution. Although the full set of bosonic functionals is equivalent to the fermionic one, for finite $N$ we could choose a mixed basis combining elements from the two. For the sake of the arguments in this section we will stick to just using the bosonic basis, and set

$$ W_N = \operatorname{span} \left\{ \alpha_0^B; \alpha_n^B, \beta_n^B, n = 1, \dots, \frac{N-1}{2} \right\}, \quad (3.13) $$

with $N$ odd.

<footnote n="9">⁹Gliozzi [17] was the first to propose that sparseness (what we call extremality) can be used to bootstrap interesting solutions to crossing, even in the absence of unitarity or positivity.</footnote>
To make our case we will focus on the OPE maximization problem of an operator with dimension $\Delta_0 < 2\Delta_\phi + 1$, the same problem which will be considered numerically in the next section. We will start by choosing $\Delta_0 = 2\Delta_\phi$ and then discuss perturbations around this point up to second order. This problem was already considered and solved in [4], but here our focus will be on understanding the dependence of the solution on the truncation size. At the end we will give a more general argument that justifies why this basis should be generally useful for any bootstrap problem.

**Perturbative arguments — zeroth order.** Consider first the case where $\Delta_0 = 2\Delta_\phi$. We must search for a functional $\omega^{\text{opt},N}$ in $W_N$ which will give us the best possible bound. In this case it is not hard to see that the optimal functional is $\alpha_0^B$, independently of $N$. Indeed it can be checked that $\alpha_0^B(\Delta) \ge 0$ for $\Delta \ge \Delta_0$ (and that $\alpha_0^B(\Delta_0) = 1$). Hence $\alpha_0^B$ definitely provides a valid bound. That this bound is optimal can be checked by computing

$$-\alpha_0^B(0) = a_{\Delta_0^B}^{\text{free}}, \qquad (3.14)$$

which establishes optimality since the generalized free boson solution (2.6) saturates it. This is not too surprising, since the exact functional bases were constructed precisely so as to be dual to the exact generalized free solutions. The extremal solution in this case will consist of $\frac{N+1}{2}$ vectors with dimensions $\Delta_0^B, \dots, \Delta_{\frac{N-1}{2}}^B$, with

$$\vec{F}_0 + a_{\Delta_0^B}^{\text{free}} \vec{F}_{\Delta_0} + \sum_{n=1}^{\frac{N-1}{2}} a_{\Delta_n^B}^{\text{free}} \vec{F}_{\Delta_n^B} = 0. \qquad (3.15)$$

That this equation holds is a simple consequence of the duality conditions (2.10). In other words, the truncated extremal solution is obtained by simply considering the first few operators in the exact generalized free boson solution, which appear with their *exact* scaling dimensions and OPE coefficients. As we add more functionals by increasing $N$ the solution to crossing systematically approaches the exact one, one operator at a time, in the best way possible. If we are only interested in obtaining the bound on $a_{\Delta_0}$, a truncation with $N=1$ will. Of course in the derivative basis the correct result would only be achieved with $N = \infty$.

**Perturbative arguments — first order.** For a more interesting result, consider instead maximizing the OPE coefficient of an operator whose dimension $\Delta_0$ is close to but not quite $2\Delta_\phi$. We expect that the extremal solution (and accordingly, the optimal bound) can be obtained by perturbing around the generalized free boson. We will do our analysis by expanding in $g \equiv \Delta_0 - 2\Delta_\phi$:

$$\Delta_m \sim \sum_{k=0}^{2} g^k \Delta_m^{(k)}, \qquad a_{\Delta_m} \sim \sum_{k=0}^{2} g^k a_m^{(k)}, \qquad (3.16)$$

where the zeroth order values are the free ones. Expanding to leading order the truncated crossing equations (3.12) and using the duality relations (2.10) we find

$$a_m^{(1)} = a_0^{(0)} c_m, \qquad a_m^{(0)} \Delta_m^{(1)} = a_0^{(0)} d_m, \qquad (3.17)$$
which fixes the extremal solution completely, including the new optimal bound

$$a_0^{\text{opt}} = a_0^{(0)} (1 + g c_0) + O(g^2). \qquad (3.18)$$

These results are completely independent of the truncation size, and are therefore exact for any $N$. So, just as before, the only role that the truncation size has is to systematically add more operators which appear with their correct dimensions and OPE coefficients (up to $O(g^2)$ corrections). In particular, we again find that the exact OPE bound is captured, to leading order in $g$, even with $N = 1$.

On the other hand, given this extremal solution to crossing we can also reconstruct the new extremal functionals. Imposing the duality relations on the modified solution we can find for instance<fnref n="10" />

$$\alpha_0 = (1 + g c_0) \left[ \alpha_0^B - g \sum_{k=1}^{\frac{N-1}{2}} \Delta_k^{(1)} \partial_{\Delta}^2 \alpha_0^B (\Delta_k^{(0)}) \beta_k^B \right] + O(g^2). \qquad (3.19)$$

Notice that this result now depends on the explicit truncation, and hence this is actually not the correct exact functional for any finite $N$. The reason for this is that this functional still has double zeros at $\Delta_m = 2\Delta_\phi + 2m$ for $m > N/2$ (i.e. above the truncation level), instead of the corrected $O(g)$ ones. In spite of this, this functional does compute the correct new bound to $O(g^2)$, $-\alpha_0(0) = a_0^{(0)} + ga_0^{(1)}$, independently of the truncation size.

**Perturbative arguments — second order.** We now go to quadratic order in $g$. In this case the crossing equations lead to:

$$a_m^{(2)} = a_0^{(1)} c_m - \frac{1}{2} \sum_{n=0}^{\frac{N-1}{2}} a_n^{(0)} (\Delta_n^{(1)})^2 \partial_{\Delta}^2 \alpha_m^B (\Delta_n^{(0)}), \qquad (3.20)$$

$$a_m^{(0)} \Delta_m^{(2)} + a_m^{(1)} \Delta_m^{(1)} = a_0^{(1)} d_m - \frac{1}{2} \sum_{n=0}^{\frac{N-1}{2}} a_n^{(0)} (\Delta_n^{(1)})^2 \partial_{\Delta}^2 \beta_m^B (\Delta_n^{(0)}).$$

The results now depend explicitly on $N$, a dependence which is directly inherited from that of the functionals at the previous order. To proceed we need to know something about the functional actions and the anomalous dimensions. Using the asymptotics derived in [4] (quoted below in equation (3.23)) as well as the large $n$ behaviour of $d_n$ from (2.11), we have:

$$\Delta_n^{(1)} \underset{n \gg 1}{\sim} 1/n^2, \qquad a_n^{(0)} \partial_{\Delta}^2 \omega_m^B (\Delta_n^{(0)}) \underset{n \gg 1}{\sim} 1/n^3, \qquad (3.21)$$

with $\omega = \alpha, \beta$, and holding $m$ fixed in the last equality. We conclude that both scaling dimensions and OPE coefficients are expected to converge to their exact ($N = \infty$) values at a rate $O(N^{-6})$.

<footnote n="10">¹⁰Note that in this process we are not free to choose the $C_n, D_n$, rather their precise values emerge automatically after imposing the other duality relations.</footnote>
**General argument.** To conclude we will do a more general analysis that does not rely on perturbation theory. As we will see in section 4, at large dimension the operator spectrum and the OPE coefficients of an extremal solution often asymptote to those of a generalized free solution. We show here that this implies that if we have solved the functional equations (numerically) with N components, and N is large enough so that the free asymptotics have kicked in, we can effectively complete the solution for arbitrarily high dimensional operators.

The full constraints of crossing symmetry on the extremal solution are written in the following way:

$$ \begin{gathered} \hat{\omega}_n(0) + \sum_{m=0}^{M} a_m \hat{\omega}_n(\Delta_m) = - \sum_{m=M+1}^{+\infty} a_m \hat{\omega}_n(\Delta_m), \quad 0 \le n \le N \tag{3.22} \\ \sum_{m=M+1} a_m \hat{\omega}_{N+n}(\Delta_m) = -\hat{\omega}_{N+n}(0) - \sum_{m=0}^{M} a_m \hat{\omega}_{N+n}(\Delta), \quad 0 < n. \end{gathered} $$

Truncating the solution means that we set to zero the righthand side of the first set of equations, while ignoring the constraints arising from the second set. We have seen that for extremal solutions we have $M = \frac{N-1}{2}$, so we make that choice, and use for the functional basis the one introduced in equation (3.13).

Let us suppose that the exact solution has a spectrum which asymptotes to the generalized free boson for high enough scaling dimension, that is, the high energy spectrum consists of “double-trace” operators whose anomalous dimensions $\gamma_m \equiv \Delta_m - \Delta_m^B$ fall off with some power of $m$. For definiteness let us say that $\gamma_m$ is sufficiently small for $m > M^*$. The bounds derived in [4] imply that the corresponding OPE coefficients must also approach the generalized free ones. In this case, as long as $M = N/2 > M^*$, the righthand side of the first set equations is indeed small for all $n < N$. This follows from the asymptotic behaviour, determined in the same reference:

$$ \begin{gathered} \beta_n(\Delta) \underset{\Delta \to \infty}{\sim} \sin\left[\frac{\pi}{2}(\Delta - \Delta_n^B)\right]^2 \left(\frac{a_{\Delta_n}^{\text{free}}}{a_{\Delta}^{\text{free}}}\right) \frac{4\Delta(\Delta_n^B)^2}{\Delta^4 - (\Delta_n^B)^4}, \quad \text{with } \Delta/\Delta_n \text{ fixed,} \tag{3.23} \\ \beta_n(\Delta) \underset{\Delta \to \infty}{\sim} \sin\left[\frac{\pi}{2}(\Delta - \Delta_n^B)\right]^2 \left(\frac{a_{\Delta_n}^{\text{free}}}{a_{\Delta}^{\text{free}}}\right) \frac{4(\Delta_n^B)^2}{\Delta^3}, \quad \text{with } \Delta_n \text{ fixed,} \end{gathered} $$

with similar expressions for the $\alpha_n$. In particular the effect of tails on the functional equations for fixed $n$ falls off as $\sim \gamma_n^2/N^2$. In the OPE maximization problem, which will also be analysed numerically in the next section, $\gamma_n \sim n^{-2}$ which leads to an overall expected $1/N^6$ convergence rate, which matches what we found in perturbation theory. More generally, we expect a $1/N^2$ fall-off behaviour until $N$ is large enough that we reach the scale where anomalous dimensions are sufficiently small, followed then by the faster $1/N^6$ fall-off.
Let us now discuss the second set of equations which we ignored. Again for large enough $N$ that we reach the weakly coupled regime, these equations become

$$ 
\begin{align} 
a_n &= -\alpha_n(0) - \sum_{m=0}^{\frac{N-1}{2}} a_m \alpha_n(\Delta_m) + O(\gamma_N^2), & n &> \frac{N-1}{2} \tag{3.24} \\ 
a_n \gamma_n &= -\sum_{m=0}^{\frac{N-1}{2}} a_m \beta_n(\Delta_m) + O(\gamma_N^2), & n &> \frac{N-1}{2} 
\end{align} 
$$

In fact, we may freely drop terms in the sums above for which $m$ is above the strongly coupled scale. These equations tell us that at high scaling dimension the CFT data is determined entirely in terms of the strongly coupled region where operators have order one anomalous dimensions.

Overall we see that using our proposed basis, the truncation procedure effectively converges as soon as $N$ is large enough that we reach the weakly coupled regime of small anomalous dimensions. This means that, in practice, we can use numeric algorithms to solve for the low energy spectrum, and then use equations (3.24) to complete the solution up to arbitrarily large scaling dimensions. The reason this is possible is because all functionals in the basis automatically encode the asymptotics of the generalized free spectrum. In particular any finite linear combination of functionals will have double zeros at the positions of the free boson spectrum for sufficiently high dimension, and so truncation preserves the asymptotics. Of course, it is clear that should the extremal solution asymptote to a generalized free fermion we would do well then to use the associated functional basis.

## 4 Applications

In this section we will test our proposed functional basis in concrete numerical applications.

In the first application, we will bootstrap the generalized free boson solution using the fermionic functional basis. This is interesting for several reasons. Firstly, we have an exact solution to compare with. Secondly, the fermionic functionals are adapted to solutions to crossing that asymptote to a spectrum of the form $\Delta_n = 1 + 2\Delta_\phi + 2n$, i.e. not the generalized free boson's $\Delta_n = 2\Delta_\phi + 2n$. From the fermionic functionals perspective the free boson is as strongly coupled as it is possible to get, so many of our motivating arguments in the previous section aren't applicable. It is interesting then to see how our basis performs in this worse case scenario.

In our second application we will maximize the OPE coefficient of the first operator above the identity, whose dimension $\Delta_0$ we take to continuously interpolate between the bosonic value $2\Delta_\phi$ and the fermionic value $1 + 2\Delta_\phi$. We find that right up to $\Delta_0 = 1 + 2\Delta_\phi$, where there is a discontinuity, the bound is saturated by a solution whose spectrum asymptotes to the generalized free boson's, and as such we will use the bosonic functional basis throughout. This is a concrete numerical test of the analytic arguments of the previous section. For $\Delta_0$ perturbatively close to $2\Delta_\phi$ the extremal solution matches with what would be obtained by adding a $\phi^4$ interaction to a free scalar field in AdS$_2$, and as such we call this family of solutions the $\phi^4$ flow.
All computations in this section were performed with JuliBootS [18]. As an input, we need to provide to this program the functional actions $\alpha_n(\Delta)$, $\beta_n(\Delta)$, analogously to $\omega^\partial(\Delta) \equiv \partial_z^{2n+1} z^{-2\Delta_\phi} G_\Delta(z)|_{z=1/2}$ for the derivative basis. Unlike the latter, for the former we do not yet have a way of obtaining efficient rational function representations, and hence we simply computed a table of functional actions evaluated at several values of $\Delta$ and fed it into the linear programming algorithm. To improve the accuracy of the overall result, interpolation was then used to iteratively generate new points on the table in the vicinity of functional zeros until convergence was observed. Further details are available upon request.

Explicit forms of the functional kernels $f_\omega(z)$ for both the bosonic and the fermionic functionals for any external dimension $\Delta_\phi$ are given in appendix A. We have chosen to focus on $\Delta_\phi = 1/2$ for the fermionic basis and $\Delta_\phi = 1$ for the bosonic one, where these $f_\omega$ take a simpler form. In this case we are able to compute the functional actions analytically, see appendix B, making the numerical evaluation of the tables faster. We emphasize that apart from this there is nothing special about these particular external dimensions: for a generic $\Delta_\phi$ it may take longer to compute the value of the functional actions, but once this is done the linear programming part of the algorithm proceeds in the same way. In particular, all nice properties of the functional bases, such as fast convergence, are generic.<fnref n="11" />

## 4.1 Bootstrapping bosons with fermions

As we have reviewed in the previous section, in 1D CFTs the OPE maximization problem for an operator with dimension $2\Delta_\phi$ is solved by the generalized free boson, with the exact functional providing this bound given by $\alpha_0^B$.<fnref n="12" /> We now attempt to reproduce this result using the fermionic functional basis and compare with the same computation as done with the derivative basis. In detail the two basis used are:

$$ 
\begin{align} 
W_N^\omega &\equiv \operatorname{span} \left\{ \alpha_n^F, \beta_n^F, \quad n = 0, \dots, \frac{N-2}{2} \right\}, \tag{4.1} \\ 
W_N^\partial &\equiv \operatorname{span} \left\{ \partial_z^{1+2n} \bigg|_{z=\frac{1}{2}}, \quad n = 0, \dots, N-1 \right\}. 
\end{align} 
$$

We focus on $\Delta_\phi = 1/2$, for which the functional actions can be computed more efficiently, as shown in appendix B. Our results are shown in figures 2 and 3. We find the functional basis outperforms the derivative basis in two different ways. Firstly, as expected for both bases the extremal solution gains a new operator every time that $N$ is increased by two units (meaning we have two additional functional components), however with functionals this operator appears already with a small initial error, which is in fact bounded in absolute value by $\sim 0.6$. No such bound is observed with the derivative basis, with this initial error seemingly growing linearly with $N$. After their initial appearance operators systematically approach their correct generalized free values as $N$ is increased further. In the derivative basis this takes longer than with functionals, since not only the initial error is larger, but also the rate at which the correct value is approached is worse. Experimentally

<footnote n="11">¹¹As a sanity check, we have explicitly checked this up to $N = 10$ for some non half-integer external dimension in the fermionic basis.</footnote>

<footnote n="12">¹²To be precise this assumes a certain minimum gap in the spectrum. This gap can be chosen to be the scaling dimension, always below $2\Delta_\phi$, where $\alpha_0^B$ first becomes negative.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Spectrum of the solution saturating an OPE maximization bound for $\Delta_0 = 2\Delta_\phi$, with $\Delta_\phi = 1/2$, using the fermionic functional basis (top) and derivatives (bottom). As $N$ is increased by two units a new operator appears, while previous operators approach their correct values shown as dashed lines. With functionals the initial error is smaller and convergence is faster.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Spectrum of the solution saturating OPE maximization bound at $\Delta_0 = 2\Delta_\phi$ with $\Delta_\phi = 1/2$. $\delta\Delta_n \equiv |\Delta_n - \Delta_n^B|$. Left: using the fermionic functional basis, as $N$ is increased by two units a new operator appears with initial bounded absolute error in dimension, which converges to about $\sim 0.6$ at large $N$. As $N$ increases further this error then rapidly decreases. Right: convergence of scaling dimension $\Delta_1$. In red the derivative basis, in blue the functional basis. The lines represent fits to a $1/N$ and $1/N^2$ behaviour respectively.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** The single bin functional $\cup$, here shown for $\Delta_\phi = \frac{1}{2}$. Contributions to the OPE in the region where the functional is negative must cancel those where it is positive. Given our gap assumptions on the extremal solution, this implies an operator must appear between $2\Delta_\phi + 1$ and $2\Delta_\phi + 3$.</figcaption></figure>

we find convergence rates of order $N^{-2}$ and $N^{-1}$ for the functional and derivative bases respectively. The former is consistent with our analysis of the previous section setting $\gamma_n \sim 1$.

The reason why the initial error of operators is bounded when using the functional basis can be understood from the results of [4]. There it was shown that by combining functionals it is possible to find lower bounds on the OPE density. This forces operators to appear in certain bins in scaling dimension space. These bounds can be obtained by combining a finite number of functionals, and hence hold even in our truncated numerical setting. To see this explicitly, in figure 4 we plot the functional combination $\cup_1 := \beta_1^F - c\beta_0^F$, where $c$ is a suitably chosen coefficient. For $N \ge 4$, this functional lies inside $W_N$. Any solution to crossing, and in particular the extremal solution for the same $N$, has to be compatible
<figure><img src="image_6.png" /><figcaption>**Figure 5.** OPE maximization at $\Delta_{\phi} = 1$ with $N = 19$ components, although no visible changes are seen beyond $N = 5$. Main plot: bound on the OPE coefficient of the first operator as a function of its dimension $\Delta_0$. Inset: dimension of the second operator $\Delta_1$ at extremality as a function of $\Delta_0$. In both plots, the dashed lines represent analytic perturbation theory computations up to cubic order in $g \equiv \Delta_1 - 2\Delta_{\phi}$, as described in the main text.</figcaption></figure>

with the constraints imposed by this functional. This constraint is that contributions from the OPE density in regions where $\cup_1$ is positive have to be cancelled by those where it is negative. Since we are imposing a gap in the spectrum up to $2\Delta_{\phi}$, the constraint from $\cup_1$ here implies a non-zero OPE density inside the region $(1 + 2\Delta_{\phi}, 3 + 2\Delta_{\phi})$. At extremality this will be achieved by a single operator inside this region, and so the error on the dimension of this operator is bounded by one unit (since the correct value is $2 + 2\Delta_{\phi}$). This argument generalizes to higher values of $n$.

## 4.2 OPE maximization, or $\phi^4$ flow

For our next application, we will derive bounds on the OPE coefficient of an operator whose dimension $\Delta_0$ is varied between $2\Delta_{\phi}$ and $2\Delta_{\phi} + 1$. As mentioned earlier, we will set $\Delta_{\phi} = 1$ for simplicity and use the bosonic functional basis since as it turns out, the extremal solution saturating the bound has a spectrum which rapidly asymptotes to that of a free boson. That is, we will now do the numerical version of what we examined in section 3.2. The basis of functionals was given in (3.13). Note that the existence of the $\alpha_0^{B,F}$ functionals imply that the exact bound will be saturated by the bosonic solution when $\Delta_0 = 2\Delta_{\phi}$ and by the fermionic solution when $\Delta_0 = 1 + 2\Delta_{\phi}$. Hence we expect that using the bosonic basis should be less optimal as we approach $\Delta_0 = 1 + 2\Delta_{\phi}$ (but still better than derivatives, since in the worse case scenario we go back to a similar situation to that of the previous subsection).
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Extremal spectrum as a function of the lowest dimension operator with $N = 19$ components. The dashed lines represent the generalized free boson spectrum, while the red dots stand for the free fermion spectrum.</figcaption></figure>

Our results are summarized in figures 5 and 6. In the main plot in figure 5 we show the upper bound on the OPE coefficient. The bound is essentially unchanged beyond truncations with $N = 5$ with even $N = 3$ already providing an excellent approximation to the bound except very close to $\Delta_0 = 1 + 2\Delta_\phi = 3$. The results shown have $N = 19$. In the inset of the same figure we show the scaling dimension $\Delta_1$ of the subleading scaling operator (since $\Delta_0$ is fixed by hand) of the extremal solution saturating the bound.<fnref n="13" /> The numerical results are in excellent agreement with the analytic results of reference [4], where CFT data for the $\phi^4$ deformation was determined up to two loops. Our numerics can be thought of as providing a non-perturbative completion of these results into the strongly coupled regime.

When $\Delta_0 = 2\Delta_\phi$ or $\Delta_0 = 1 + 2\Delta_\phi$ the OPE bound is saturated by the free values 2 and $2\Delta_\phi$ respectively (which for $\Delta_\phi = 1$ are the same), as expected, and $\Delta_1$ also correctly interpolates between $2\Delta_\phi + 2 = 4$ and $1 + 2\Delta_\phi + 2 = 5$. In between those values, we have a non-trivial bound and an interacting solution to crossing. Figure 6 shows a more complete picture of the spectrum obtained from the same truncation with $N = 19$. We see that most of the spectrum is essentially identical to that of the free boson for almost all values of $\Delta_0$, which explains why using the bosonic functional basis is such a good idea here. It is only when $\Delta_0$ approaches $1 + 2\Delta_\phi$ that we see rapid rearrangements in the spectrum. For sufficiently large $\Delta$ the spectrum essentially jumps nearly discontinuously from that of the free boson to the free fermion one. This is consistent with the fact that there are no UV-

<footnote n="13">¹³A different way to think of this curve is as an upper bound on the dimension of the first subleading scalar in a solution to crossing, where the leading scalar is assumed to have dimension $\Delta_0$. In other words, the extremal solution that saturates the OPE bound can also be found by doing gap maximization on the subleading scalar.</footnote>
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Bound on $a_0$ for $\Delta_0 = 2\Delta_\phi + \frac{1}{2} = 2.5$ as a function of the number of components. In black the bound obtained using functionals, in red the one with derivatives. The dashed line corresponds to a prolongation of the $N = 19$ functional result. The inset allows us to better see how fast the functional basis result converges. In particular these results show that in this case the bound using the functional basis with $N = 3$ significantly outperforms the derivative basis with $N = 90$.</figcaption></figure>

preserving deformations of the generalized free fermion that do not introduce new states, as follows from the duality conditions (2.9). Presumably irrelevant contact interactions in AdS$_2$ should account for the spectrum close to the free fermion point.

We observe that as we increase $N$, new operators appear with a value of the scaling dimension which is essentially already their correct final value. That is, after an operator appears, increasing $N$ does not significantly change its OPE or scaling dimension. In this sense, we find quite remarkably that the solution to the numerical problem converges. This is in line with the expectations of section 3.2: once we have successfully bootstrapped the region of the spectrum where scaling dimensions are very different from the free ones, there is very little interest in pursuing the computation, since the OPE data of higher dimension operators may be determined by equations (3.24). We have checked that those equations are indeed satisfied to high accuracy.

It is illuminating to compare the relative merits of the functional and derivative bases in the current setup. In figure 7 we compare the bounds on $a_0$ obtained using the derivative and the functional bases for the case $\Delta_0 = 2\Delta_\phi + \frac{1}{2} = 2.5$, which is halfway between the cases $\Delta_0 = 2\Delta_\phi = 2$ and $\Delta_0 = 2\Delta_\phi + 1 = 3$ (where the theories saturating the bounds are the free boson and fermion respectively) and hence “maximally interacting” in some sense. We see that the functional basis outperforms the derivative basis dramatically: using just three functional components in this case outperforms (at least) 90 derivative components. In figure 8 we show more detailed convergence properties of the functional
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Convergence of OPE maximization. Shown are the relative errors in the determination of $\Delta_1$ (as compared to its value for $N = 19$), for two distinct values of $\Delta_0$, namely close to the bosonic and fermionic points. Although the initial error is smaller in the former (since we are using the bosonic basis), the rate of convergence is the same in either case and compatible with a $1/N^6$ behaviour.</figcaption></figure>

basis. Although it is hard to estimate the correct fall-off power from our data (we have considered only functionals up to $N = 19$), the results are consistent with a very fast convergence rate $\sim N^{-6}$, as argued in section 3.2. For the derivative basis, there is no special difference between setting $\Delta_0 = 2\Delta_\phi$ or some other value, and so the convergence properties look the same as those of the previous section. In particular, one again finds that increasing $N$ successive operators appear in with larger and larger initial errors, and converge then with a power law falloff $\sim 1/N$.

## 5 Outlook — higher dimensions

The main conclusion of this note is that the subset of constraints arising from the $z = \bar{z}$ section of the crossing equation are efficiently captured by the functional bases introduced in [4]. In particular, we have argued and shown that for 1D CFTs these functionals provide what is perhaps the best possible choice of basis with which to compute numerical bounds. Our results make even more exciting the possibility that there are analytic functionals in higher dimensions, analogous to the ones in 1D, which would allow us to distinguish between different spin channels. If such functionals can be constructed, even discounting their analytic implications, they will surely revolutionize what is possible to do in the numerical bootstrap. Since such functionals are not yet available, we would like to comment on what can be done right away using the functionals we do have.

Upon restriction to the $z = \bar{z}$ line, higher dimension solutions to crossing equation take the form.

$$ \sum_{\Delta, \ell} a_{\Delta, \ell} F_{\Delta, \ell}(z, z|\Delta_{\phi}) = 0. \qquad (5.1) $$
Hence there are now many kinds of functional actions, depending on spin:

$$ \omega(\Delta) \rightarrow \omega(\Delta, \ell). \qquad (5.2) $$

However, it is straightforward to show that the functional actions still encode generalized free asymptotics in every spin channel, i.e. we have relations such as

$$ \alpha_n^B(\Delta_m^B, \ell) = \delta_{n,m}, \quad \text{for } m \ge n \qquad (5.3) $$

independently of spin.<fnref n="14" />

When bootstrapping higher dimensional CFTs, the class of functionals typically used amounts to derivatives in the two cross-ratios $z, \bar{z}$, which may be taken along directions parallel and transverse to the $z = \bar{z}$ line. So, one obvious thing to do is to replace purely parallel derivatives by the functional basis proposed here. More generally one could imagine Taylor expanding the $F_{\Delta, \ell}(z, \bar{z})$ along the $z = \bar{z}$ line and acting with our functionals on each coefficient function. However, there are good reasons to believe that by simply adding our functionals to the existing approach will already lead to significant and in some cases dramatic improvements.

The first basic reason is that even the set of $z = \bar{z}$ constraints in higher dimensions are non-trivial, especially for moderate and large values of $\Delta_\phi$. Consider for instance the typical bootstrap problem, which is a bound on the leading scalar dimension, while allowing for higher spin operators with dimensions above or equal to the unitarity bound $d-2+\ell$. Since all higher D CFTs are also $D=1$ CFTs, this bound must be equal to or stronger than the 1D result $1+2\Delta_\phi$ for $1+2\Delta_\phi \le d$.<fnref n="15" /> The point now is that even mustering sufficient computation firepower to obtain a bound as strong as the $D=1$ result takes some work. As a simple example for instance, in $D=4$ (cf. [19], figure 2) with $\Delta_\phi = 3/2$, the usual derivative bound only becomes competitive with the $D=1$ result after using at least 21 derivative components, whereas the same result would be obtained with a single functional ($\beta_0^F$).

More interestingly, higher dimensional bounds can in fact sometimes be saturated by 1D solutions, where our functionals can then provide the optimal bound immediately. A trivial example is that if we take the free theory solution in any spacetime dimension, from a 1D perspective this is nothing but the generalized free boson, for which we already have associated extremal functionals. In particular, this predicts that the free theory solution maximizes the OPE bound on the scalar operator of dimension $2\Delta_\phi = d-2$. Hence, in the neighbourhood of the free points, using the functional basis above should not be too bad, since using a truncation containing a single functional $\alpha_0^B$ (i.e. the analog of using a single derivative) already yields the exact optimal bound at that point.

Two more interesting examples occur in $D=2$. There the bound on the leading scalar dimension looks schematically as in figure 9.<fnref n="16" /> The first example is at the point

<footnote n="14">¹⁴This is a consequence that the derivation of the functional action (2.14), with the characteristic sine squared oscillations, actually depends on very little details of the conformal blocks.</footnote>

<footnote n="15">¹⁵When we set a gap from a 1D perspective we are gapping all spin channels in higher dimensions. The restriction then arises since we don't want to gap spin 2 operators above the unitarity bound.</footnote>

<footnote n="16">¹⁶The first 2D bound appeared in [20], and a wider plot where the second kink at (1,4) is visible has appeared recently in [21]. However, the existence of this kink was communicated to us by S. El-Showk a long time ago, who seems to have been the first to notice it.</footnote>
<figure><img src="image_10.png" /><figcaption>**Figure 9.** Schematic bound on leading scalar dimension $\Delta_{\epsilon}$ in 2D CFTs. Three notable points $(\Delta_{\phi}, \Delta_{\epsilon})$ are the 2D Ising model $\langle\sigma\sigma\sigma\sigma\rangle$ correlator at (1/8, 1), the large $m$ limit of the minimal model $\phi_{1,2}$ four-point correlator at (1/2, 2) and the 2D Ising $\langle\epsilon\epsilon\epsilon\epsilon\rangle$ correlator at (1, 4). The dashed lines correspond to the 1D generalized free fermion ($\Delta_{\epsilon} = 1+2\Delta_{\phi}$) and boson ($\Delta_{\epsilon} = 2\Delta_{\phi}$) theories.</figcaption></figure>

($\Delta_{\phi}, \Delta_{\epsilon}$) = (1/2, 2), where the bound crosses the 1D generalized free fermion line. This point corresponds to the $m \to \infty$ limit of the minimal model correlator of four $\phi_{1,2}$ [22],$^{17}$ which for $z = \bar{z}$ is indeed the same as the generalized free fermion. Hence, at this point we already have the extremal functional available which determines the exact optimal bound, namely $\beta_0^F$. Establishing this with the derivative basis, even approximately, is highly non-trivial. The bound determined by $\beta_0^F$ alone is suboptimal everywhere below $\Delta_{\phi} = 1/2$ as it should, and above this value it gaps the spin-2 channel as well, so there is no inconsistency with higher-D results.

The second example is at $(\Delta_{\phi}, \Delta_{\epsilon}) = (1, 4)$ which is the four $\epsilon$ correlator for the 2D Ising model. This correlator restricted to the $z = \bar{z}$ line becomes a generalized free boson. From the 1D perspective the operator of dimension $2\Delta_{\phi} = 2$ is the 2D stress tensor, and $\epsilon$ is the operator with dimension $2+2\Delta_{\phi} = 4$. In particular the OPE bound on an operator of dimension $2\Delta_{\phi}$ is saturated by the free boson, which in higher $D$ terms translates into the statement that central-charge minimization will be saturated by the $\epsilon$ four point function, with corresponding functional $\alpha_0^B$.

In both cases, we see that in our basis truncations with even a single component already yield exact optimal bounds, so that transverse derivatives are only required to disentangle the spin structure of the extremal spectrum. It is reasonable to expect that in the vicinity of these points including our functionals in the basis should lead to great improvements in current numerical results. It would be very interesting to investigate this in a more systematic fashion.

<footnote>$^{17}$This is also the correlator of four energy operator in the 4-state Potts model at criticality, see for example [23].</footnote>
Acknowledgments

We would like to acknowledge A. Kaviraj, D. Mazáč, S. Rychkov, A. Vichi and especially E. Trevisani for discussions and comments at various points of this collaboration. BZ is supported by the National Centre of Competence in Research SwissMAP funded by the Swiss National Science Foundation.

A Functional kernels

In this section we show how to obtain the explicit kernels corresponding to functionals $\alpha_n, \beta_n$. For further details the reader should consult [4]. We start with the results:

$$f_{\beta_0^F}(z|\Delta_\phi) = -\kappa(\Delta_\phi) \frac{2z-1}{w^{3/2}} \left[ 3\tilde{F}_2\left(-\frac{1}{2}, \frac{3}{2}, 2\Delta_\phi + \frac{3}{2}; \Delta_\phi + 1, \Delta_\phi + 2; -\frac{1}{4w}\right) \right. \quad \text{(A.1a)} \\ \left. + \frac{9}{16w} 3\tilde{F}_2\left(\frac{1}{2}, \frac{5}{2}, 2\Delta_\phi + \frac{5}{2}; \Delta_\phi + 2, \Delta_\phi + 3; -\frac{1}{4w}\right) \right],$$

$$f_{\alpha_0^F}(z|\Delta_\phi) = \kappa(\Delta_\phi) \frac{2(z-2)(z+1)}{(2z-1)w^{3/2}} \left[ \begin{aligned} & 3\tilde{F}_2\left(-\frac{1}{2}, -\frac{1}{2}, 2\Delta_\phi + \frac{3}{2}; \Delta_\phi + 2, \Delta_\phi + 2; -\frac{1}{4w}\right) \quad \text{(A.1b)} \\ & + \frac{(2\Delta_\phi + 3)(2\Delta_\phi + 5)}{16w} 3\tilde{F}_2\left(\frac{1}{2}, \frac{1}{2}, 2\Delta_\phi + \frac{5}{2}; \Delta_\phi + 3, \Delta_\phi + 3; -\frac{1}{4w}\right) \\ & - \frac{3(4\Delta_\phi + 5)}{256w^2} 3\tilde{F}_2\left(\frac{3}{2}, \frac{3}{2}, 2\Delta_\phi + \frac{7}{2}; \Delta_\phi + 4, \Delta_\phi + 4; -\frac{1}{4w}\right) \end{aligned} \right].$$

where $w \equiv z(z-1)$ and we have added an extra argument to the kernel $f(z) \to f(z|\Delta_\phi)$ to keep track of which $\Delta_\phi$ it corresponds to. Here $3\tilde{F}_2$ stands for the regularized hypergeometric function, $w = z(z-1)$ and the normalization factor is given by

$$\kappa(\Delta_\phi) = \frac{\Gamma(4\Delta_\phi + 4)}{2^{8\Delta_\phi + 5} \Gamma(\Delta_\phi + 1)^2}. \quad \text{(A.2)}$$

Recall that $g_\omega(z) = -(1-z)^{2\Delta_\phi-2}f_\omega(\frac{1-z}{z})$ in the fermionic case (for bosons, the overall minus sign isn't there). We can now construct shifted functionals:

$$f_{s\omega_m^F}(z|\Delta_\phi) = \left[ \frac{1+z(z-1)}{z(z-1)} \right]^{2m} f_{\omega_0^F}(z|\Delta_\phi + m) \quad \text{(A.3)}$$

where $\omega$ can be $\alpha$ or $\beta$. The corresponding shifted functionals $s\beta_m^F$, $s\alpha_m^F$ satisfy the correct duality conditions for $\Delta_n > \Delta_m$, and hence we have e.g.

$$\beta_n^F = \sum_{m=0}^{n} q_n s\alpha_n^F + r_n s\beta_n^F \quad \text{(A.4)}$$

for some $q_n, r_n$. The constants may be determined by imposing the full duality relations (2.9). We note that the solution to this orthonormalisation step is not currently known in closed form for all $n$, i.e. it has to be done case by case, except for special values of $\Delta_\phi$. However, it is clear that for the purpose of the numerical computations we are
interested in any basis will do, and so we may as well work with the shifted functionals instead of the orthonormal ones.

Now let us turn to the bosonic functionals. We now define:

$$f_{s\omega_{m+1}^B}(z|\Delta_\phi) = -\frac{1}{z(z-1)} \left[ \frac{1+z(z-1)}{z(z-1)} \right]^{2m} f_{\omega_0^F}(z|\Delta_\phi + 3/2 + m). \quad (A.5)$$

That is, the bosonic kernels can be obtained from the fermionic ones. The kernels above correspond to shifted bosonic functionals analogous to the fermionic ones, and again an orthonormalisation procedure may be applied if so wished, by imposing the duality conditions (2.10). The above should be used for $m \ge 0$. The $\alpha_0^B$ functional must be treated separately, and is given by:

$$f_{\alpha_0^B}(z|\Delta_\phi) = -z(z-1) \left[ f_{\alpha_0^F}(z|\Delta_\phi - 3/2) + \frac{1}{2\Delta_\phi - 1} f_{\beta_0^F}(z|\Delta_\phi - 3/2) \right]. \quad (A.6)$$

## B Computation of the functional action for special cases

In this section we compute the action of the bosonic functionals for external dimension $\Delta_\phi = 1$ and the fermionic ones for $\Delta_\phi = 1/2$. We will give the explicit forms of these functionals and show how to evaluate them. The techniques shown here apply in the bosonic case with integer external dimension and fermionic case with half integer external dimension, where the functional kernels have a simpler form.

### B.1 The functional kernels

The functionals act on $F_\Delta(z)$ via the formula:

$$\omega_n[F_\Delta] \equiv \omega_n(\Delta) = \frac{1}{2} \int_{\frac{1}{2}}^{\frac{1}{2}+i\infty} dz f_{\omega_n}(z) F_\Delta(z) + \int_{\frac{1}{2}}^{1} dz g_{\omega_n}(z) F_\Delta(z) \quad (B.1)$$

where $g(z) = \pm(1-z)^{2\Delta_\phi-2}f(\frac{1}{1-z})$ and the $+(-)$ sign for the boson (fermion) basis respectively. This can be massaged into the alternative expression:

$$\omega_n(\Delta) = [1 \pm \cos \pi(\Delta - 2\Delta_\phi)] \int_0^1 dz g_{\omega_n}(z) \frac{G_\Delta(z)}{z^{2\Delta_\phi}}. \quad (B.2)$$

For the fermionic functionals at external dimension $\Delta_\phi = 1/2$, the functional kernels are

$$g_{\beta_n}^{-}(z) = \frac{2\Gamma(2+2n)^2}{\pi^2\Gamma(3+4n)} \left[ \frac{1}{z} P_{1+2n}\left(\frac{2-z}{z}\right) - P_{1+2n}(2z-1) \right],$$

$$g_{\alpha_n}^{-}(z) = \frac{2\Gamma(2+2n)^2}{\pi^2\Gamma(3+4n)} \left[ \frac{1}{z} DP_{1+2n}\left(\frac{2-z}{z}\right) - DP_{1+2n}(2z-1) - \frac{\Gamma(2+2n)^2}{\Gamma(4+4n)} \frac{G_{\Delta_n}(1-z)}{1-z} \right] \\ = \frac{1}{2} \partial_n g_{\beta_n}^{-}(z) - \frac{2}{\pi^2} \frac{\Gamma(2+2n)^4}{\Gamma(3+4n)\Gamma(4+4n)} \frac{G_{\Delta_n}(1-z)}{1-z}.$$
where $\Delta_n = 2n + 2$, $P_n(z)$ is the Legendre polynomial and the derivative with respect to its parameter is:

$$ 
\begin{aligned}
DP_n(y) &:= P_n(y) \log \left( \frac{1+y}{2} \right) + 2 \sum_{k=0}^{n-1} \left[ (-1)^{k+n} \frac{2k+1}{(n-k)(n+k+1)} P_k(y) \right] \\
&= \frac{\Gamma(2n+1)}{\Gamma(n+1)^2} \partial_n \left[ \frac{\Gamma(n+1)^2}{\Gamma(2n+1)} P_n(y) \right].
\end{aligned} 
\quad (B.3) $$

The motivated reader can check these expressions are perfectly compatible with those of appendix A. The bosonic functional kernels for external dimension are $\Delta_\phi = 1$ are instead

$$ 
\begin{aligned}
g_{\beta_m}^{+}(z) &= \frac{2}{\pi^2} \frac{\Gamma(2+2m)^2}{\Gamma(3+4m)} \left[ P_{2m+1}\left(\frac{2-z}{z}\right) + P_{2m+1}(2z-1) - P_1\left(\frac{2-z}{z}\right) - P_1(2z-1) \right], \\
g_{\alpha_m}^{+}(z) &= \frac{1}{2} \partial_m g_{\beta_m}^{+}(z) - \frac{2}{\pi^2} \frac{\Gamma(2+2m)^4}{\Gamma(3+4m)\Gamma(4+4m)} G_{\Delta_m}(1-z).
\end{aligned} 
\quad (B.4) $$

Numerically, the integral (B.2) converges only for $\Delta > \Delta_n$, and below that value we must use (B.1). Instead, we will now compute (B.2) exactly for $\Delta > \Delta_n$ and analytically continue the result to general $\Delta$.

## B.2 Computation of beta functionals

Let us start by showing how to compute the integral for the $\beta$ functionals. First of all we notice that $P_{1+2n}(\frac{2}{z}-1)$ is an eigenfunction of the Casimir $C_2 = z^2(1-z)\frac{d^2}{dz^2} - z^2\frac{d}{dz}$ with eigenvalue $c_2(\Delta_n) = \Delta_n(\Delta_n-1) = (2n+2)(2n+1)$. We use the fact that the measure $z^{-2}$ makes the Casimir operator self adjoint

$$ \int_0^1 \frac{dz}{z^2} f(z) (C_2 g(z)) = [(1-z)(f(z)g'(z) - f'(z)g(z))]_0^1 + \int_0^1 \frac{dz}{z^2} (C_2 f(z)) g(z). \quad (B.5) $$

Combining this with the fact that both $G_\Delta(z)$ and $P_{2n+1}(\frac{2}{z}-1)$ are eigenfunctions of the Casimir, we get

$$ \int_0^1 \frac{dz}{z^2} P_{1+2n} \left( \frac{2}{z} - 1 \right) G_\Delta(z) = \frac{1}{c_2(\Delta) - c_2(\Delta_n)} \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2}. \quad (B.6) $$

The l.h.s. converges for $\Delta > \Delta_n$ only, but we use this formula as the analytic continuation. There is a pole at $\Delta = \Delta_n$, which when combined with the prefactor in (B.2) is responsible for the single zero of $\beta_n$, see (2.9).

There is no such trick for the second part of the $\beta$ functional. Evaluating the Legendre polynomial for a given $n$ just leaves us to evaluate a finite number of integrals of the kind

$$ I_a(\Delta) = \int_0^1 \frac{dz}{z^2} z^a G_\Delta(z). \quad (B.7) $$

This can be evaluated for generic values of $a$ and turns out to be

$$ I_a(\Delta) = \frac{3F_2(\Delta, \Delta, a + \Delta - 1; 2\Delta, a + \Delta; 1)}{a + \Delta - 1}. \quad (B.8) $$
However, it can be convenient to compute $I_a(\Delta)$ recursively, since this is numerically faster. Using the fact that the conformal block is an eigenvalue of the Casimir operator

$$ 
\begin{aligned}
0 &= \int_0^1 \frac{dz}{z^2} z^a (C_2 - c_2(\Delta)) G_\Delta(z) \\
&= \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2} - \int_0^1 \frac{dz}{z^2} C_2(z^a) G_\Delta(z) - c_2(\Delta) I_a(\Delta) \\
&= \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2} + (c_2(a) - c_2(\Delta)) I_a(\Delta) - a^2 I_{a+1}(\Delta),
\end{aligned} 
\quad (B.9) $$

we get the recursion relation

$$ I_{a+1}(\Delta) = \frac{1}{a^2} \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2} + \frac{c_2(a) - c_2(\Delta)}{a^2} I_a(\Delta). \quad (B.10) $$

Now, we can find the value of $I_a(\Delta)$ for any $a \in \mathbb{Z}$ using that

$$ I_0(\Delta) = \frac{\Gamma(2\Delta)}{(\Delta - 1)\Delta\Gamma(\Delta)^2}, \quad (B.11) $$

$$ I_1(\Delta) = \frac{\Gamma(2\Delta)\Gamma(\Delta + 1)\left(\psi^{(1)}\left(\frac{\Delta}{2}\right) - \psi^{(1)}\left(\frac{\Delta+1}{2}\right)\right)}{2\Gamma(\Delta)^3}, \quad (B.12) $$

where for the second one we used equation (B.11) of [24]. By using this recursion relation there is no need to evaluate ${}_3F_2$ at unit argument, and this helps speeding up the numerics.

## B.3 Computation of the alpha functionals

The computation of the $\alpha$ functional is a bit more complicated. Using the derivative relation (B.3) between $DP_n$ and $P_n$ and the result (B.6) we get that

$$ 
\begin{aligned}
\int_0^1 \frac{dz}{z^2} DP_{2n+1} \left(\frac{2}{z}-1\right) G_\Delta(z) &= \\
&= \left[ 2(H_{2n+1}-H_{4n+2}) + \frac{2\Delta_n-1}{c_2(\Delta)-c_2(\Delta_n)} \right] \frac{1}{c_2(\Delta)-c_2(\Delta_n)} \frac{\Gamma(2\Delta)}{\Gamma(\Delta)^2},
\end{aligned} 
\quad (B.13) $$

where $H_n$ is the n-th Harmonic number. We can see that this integral has a double pole at $\Delta = \Delta_n$, and therefore $\alpha_n$ is finite at this point, see (2.9).

The term with $G_{\Delta_n}(1-z)$ can be computed in terms of ${}_4F_3$ with unit argument by using the result (3.38) of [25], which we report here for completeness:

$$ 
\begin{aligned}
\int_0^1 \frac{dz}{z^2} \left(\frac{z}{1-z}\right)^p G_\Delta(z) G_{\Delta_n}(1-z) &= \\
&= \frac{\Gamma(2\Delta_n)\Gamma(p+\Delta-1)^2\Gamma(-p-\Delta+\Delta_n+1)}{\Gamma(\Delta_n)^2\Gamma(p+\Delta+\Delta_n-1)} \\
&\quad \times {}_4F_3(\Delta, \Delta, p+\Delta-1, p+\Delta-1; 2\Delta, p+\Delta-\Delta_n, p+\Delta+\Delta_n-1; 1) + \\
&\quad + \frac{\Gamma(2\Delta)\Gamma(p+\Delta-\Delta_n-1)\Gamma(-p+\Delta_n+1)^2}{\Gamma(\Delta)^2\Gamma(-p+\Delta+\Delta_n+1)} \\
&\quad \times {}_4F_3(\Delta_n, \Delta_n, -p+\Delta_n+1, -p+\Delta_n+1; 2\Delta_n, -p-\Delta+\Delta_n+2, -p+\Delta+\Delta_n+1; 1).
\end{aligned} 
\quad (B.14) $$
Finally, we are left with expanding the $DP_{2n+1}(2z - 1)$ term and doing a finite number of integrals of the kind

$$J_a(\Delta) = \int_0^1 \frac{dz}{z^2} z^a \log(z) G_{\Delta}(z), \qquad (B.15)$$

which have the solution

$$J_a(\Delta) = - \frac{\mathfrak{F}_3(\Delta, \Delta, a + \Delta - 1, a + \Delta - 1; 2\Delta, a + \Delta, a + \Delta; 1)}{(a + \Delta - 1)^2}. \qquad (B.16)$$

When computing the action of the functionals the bottleneck is given by the evaluation of the $\mathfrak{F}_3$ at unit argument, and so we would like to compute the least possible amount of them. By noting that $J_a(\Delta) = \partial_a I_a(\Delta)$ and deriving formula (B.10), we get

$$J_{a+1}(\Delta) = - \frac{2 \Gamma(2\Delta)}{a^3 \Gamma(\Delta)^2} + \frac{a + 2c_2(\Delta)}{a^3} I_a(\Delta) + \frac{c_2(a) - c_2(\Delta)}{a^2} J_a(\Delta) \qquad (B.17)$$

so that every $J_a(\Delta)$ can be expressed in terms of $J_1(\Delta)$ and $J_0(\Delta) = - \frac{\Gamma(2\Delta)}{(\Delta-1)^2\Gamma(\Delta+1)^2}$.

**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] 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].

[2] D. Poland, S. Rychkov and A. Vichi, *The conformal bootstrap: theory, numerical techniques and applications*, Rev. Mod. Phys. 91 (2019) 015002 [arXiv:1805.04405] [INSPIRE].

[3] A. Castedo Echeverri, B. von Harling and M. Serone, *The effective bootstrap*, JHEP 09 (2016) 097 [arXiv:1606.02771] [INSPIRE].

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

[5] A. Kaviraj and M.F. Paulos, *The functional bootstrap for boundary CFT*, JHEP 04 (2020) 135 [arXiv:1812.04034] [INSPIRE].

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

[7] D. Mazáč, *A crossing-symmetric OPE inversion formula*, JHEP 06 (2019) 082 [arXiv:1812.02254] [INSPIRE].

[8] D. Mazac and M.F. Paulos, *The analytic functional bootstrap. Part I: 1D CFTs and 2D S-matrices*, JHEP 02 (2019) 162 [arXiv:1803.10233] [INSPIRE].

[9] D. Mazáč, L. Rastelli and X. Zhou, *An analytic approach to BCFT$_d$, JHEP 12 (2019) 004 [arXiv:1812.09314] [INSPIRE].

[10] S. El-Showk and M.F. Paulos, *Extremal bootstrapping: go with the flow*, JHEP 03 (2018) 148 [arXiv:1605.08087] [INSPIRE].
[11] 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].

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

[13] X. Zhou, Recursion relations in Witten diagrams and conformal partial waves, JHEP 05 (2019) 006 [arXiv:1812.01006] [INSPIRE].

[14] R. Rattazzi, S. Rychkov and A. Vichi, Central charge bounds in 4D conformal field theory, Phys. Rev. D 83 (2011) 046011 [arXiv:1009.2725] [INSPIRE].

[15] R. Reemtsen and J.-J. Rückmann eds., Semi-infinite programming, Springer, U.S.A. (1998).

[16] N. Afkhami-Jeddi, T. Hartman and A. Tajdini, Fast conformal bootstrap and constraints on 3d gravity, JHEP 05 (2019) 087 [arXiv:1903.06272] [INSPIRE].

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

[18] M.F. Paulos, JuliBootS: a hands-on guide to the conformal bootstrap, arXiv:1412.4127 [INSPIRE].

[19] D. Poland, D. Simmons-Duffin and A. Vichi, Carving out the space of 4D CFTs, JHEP 05 (2012) 110 [arXiv:1109.5176] [INSPIRE].

[20] V.S. Rychkov and A. Vichi, Universal constraints on conformal operator dimensions, Phys. Rev. D 80 (2009) 045006 [arXiv:0905.2211] [INSPIRE].

[21] C.N. Gowdigere, J. Santara and Sumedha, Conformal bootstrap signatures of the tricritical Ising universality class, Phys. Rev. D 101 (2020) 116020 [arXiv:1811.11442] [INSPIRE].

[22] P. Liendo, L. Rastelli and B.C. van Rees, The bootstrap program for boundary CFT$_d$, JHEP 07 (2013) 113 [arXiv:1210.4258] [INSPIRE].

[23] V. Gorbenko, S. Rychkov and B. Zan, Walking, weak first-order transitions, and complex CFTs II. Two-dimensional Potts model at $Q > 4$, SciPost Phys. 5 (2018) 050 [arXiv:1808.04380] [INSPIRE].

[24] M. Milgram, On hypergeometrics $3F2(1) - a$ review, arXiv:1011.4546.

[25] J. Liu, E. Perlmutter, V. Rosenhaus and D. Simmons-Duffin, $d$-dimensional SYK, AdS loops and $6j$ symbols, JHEP 03 (2019) 052 [arXiv:1808.00612] [INSPIRE].