## Bootstrapping bulk locality Part II: Interacting functionals

Nat Levine a,b and Miguel F. Paulos a

a Laboratoire de Physique, b Institut Philippe Meyer, ´ cole Normale Sup´rieure, E e Universit´ PSL, CNRS, Sorbonne Universit´, Universit´ Paris Cit´, e e e e 24 rue Lhomond, F-75005 Paris, France

nat.levine@phys.ens.fr , miguel.paulos@ens.fr

## Abstract

Locality of bulk operators in AdS imposes stringent constraints on their description in terms of the boundary CFT. These constraints are encoded as sum rules on the bulkto-boundary expansion coefficients. In this paper, we construct families of sum rules that are (i) complete and (ii) 'dual' to sparse CFT spectra. The sum rules trivialise the reconstruction of bulk operators in strongly interacting QFTs in AdS space and allow us to write down explicit, exact, interacting solutions to the locality problem. Technically, we characterise 'completeness' of a set of sum rules by constructing Schauder bases for a certain space of real-analytic functions. In turn, this allows us to prove a Paley-Weiner type theorem characterising the space of sum rules. Remarkably, with control over this space, it is possible to write down closed-form 'designer sum rules', dual to a chosen spectrum of CFT operators meeting certain criteria. We discuss the consequences of our results for both analytic and numerical bootstrap applications.

## Contents

| 1 Introduction   | 1 Introduction                                           | 1 Introduction                                                                | 1 Introduction                                                                |   3 |
|------------------|----------------------------------------------------------|-------------------------------------------------------------------------------|-------------------------------------------------------------------------------|-----|
| 2                | Locality bootstrap and hyperfunctions                    | Locality bootstrap and hyperfunctions                                         | Locality bootstrap and hyperfunctions                                         |  11 |
|                  | 2.1                                                      | Form factors as hyperfunctions . . . . . . . .                                | Form factors as hyperfunctions . . . . . . . .                                |  11 |
|                  | 2.2                                                      | Functionals and locality conditions . . . . . . . . . . .                     | Functionals and locality conditions . . . . . . . . . . .                     |  14 |
|                  | 2.3                                                      | Generalized Free Field functionals . . . . . . . . .                          | Generalized Free Field functionals . . . . . . . . .                          |  17 |
| 3                | A general construction: interacting functionals          | A general construction: interacting functionals                               | A general construction: interacting functionals                               |  19 |
|                  | 3.1 A useful L 2 space . . . . . . . . . . . . . . . . . | 3.1 A useful L 2 space . . . . . . . . . . . . . . . . .                      | 3.1 A useful L 2 space . . . . . . . . . . . . . . . . .                      |  21 |
|                  |                                                          | 3.1.1                                                                         | Embeddings . . . . . . . .                                                    |  21 |
|                  |                                                          | 3.1.2                                                                         | Riesz bases . . . . . . . .                                                   |  22 |
|                  |                                                          | 3.1.3                                                                         | How to use L 2 . . . . . . .                                                  |  23 |
|                  | 3.2                                                      | Completeness of functionals using L 2 . . . . . . . . . . . . . . . . . .     | Completeness of functionals using L 2 . . . . . . . . . . . . . . . . . .     |  25 |
|                  |                                                          | 3.2.1 GFF functionals . . . . . . . .                                         | 3.2.1 GFF functionals . . . . . . . .                                         |  25 |
|                  |                                                          | 3.2.2                                                                         | Interacting functionals using                                                 |  27 |
|                  | 3.3                                                      | On general values of ˜ α and 'subtractions' . . . .                           | On general values of ˜ α and 'subtractions' . . . .                           |  29 |
|                  | 3.4                                                      | Completeness and extremality .                                                | .                                                                             |  31 |
|                  |                                                          | 3.4.1                                                                         | Completeness . . . . . . .                                                    |  31 |
|                  |                                                          | 3.4.2                                                                         | Extremal solutions . . . .                                                    |  32 |
| 4                | The space of functional actions                          | The space of functional actions                                               | The space of functional actions                                               |  33 |
|                  | 4.1                                                      | An intrinsic definition .                                                     | . . . . .                                                                     |  34 |
|                  | 4.2                                                      | Exact interacting functional actions: a product formula . . . . . . . . . . . | Exact interacting functional actions: a product formula . . . . . . . . . . . |  36 |
| 5                | Application: a case study in extremality                 | Application: a case study in extremality                                      | Application: a case study in extremality                                      |  38 |
|                  | 5.1                                                      | Preamble: general expectations .                                              | Preamble: general expectations .                                              |  38 |
|                  | 5.2                                                      | An extremal family . .                                                        | . . . . . .                                                                   |  41 |
|                  | 5.3                                                      | Interacting basis: perturbative computation . . .                             | Interacting basis: perturbative computation . . .                             |  42 |
| 6                | Discussion and outlook                                   | Discussion and outlook                                                        | Discussion and outlook                                                        |  47 |
|                  | 6.1                                                      | Extremal solutions vs bases . . . . . . . .                                   | Extremal solutions vs bases . . . . . . . .                                   |  47 |
|                  | 6.2                                                      | Implications for numerics . . .                                               | . .                                                                           |  49 |
|                  | 6.3                                                      | Outlook . . . . . . . . . . . . . . . . . . .                                 | Outlook . . . . . . . . . . . . . . . . . . .                                 |  50 |

| B   | Riesz bases                                       | Riesz bases                                                   |   53 |
|-----|---------------------------------------------------|---------------------------------------------------------------|------|
|     | B.1                                               | GFF Riesz bases . . . . . . . . . . . . . . . . . . . .       |   54 |
|     | B.2                                               | Interacting Riesz bases . . . . . . . . . . . . . . . . .     |   60 |
| C   | Analyticity and bounds on interacting functionals | Analyticity and bounds on interacting functionals             |   63 |
|     | C.1                                               | Functional actions and analyticity . . . . . . . . . . .      |   63 |
|     | C.2                                               | Analyticity of functional kernels . . . . . . . . . . . .     |   67 |
|     | C.3                                               | Polynomial bound at z = ∞ . . . . . . . . . . . . . .         |   72 |
|     | C.4                                               | Polynomial bound at z = 1 . . . . . . . . . . . . . . .       |   72 |
| D   | On general values of ˜ α and subtractions         | On general values of ˜ α and subtractions                     |   74 |
|     | D.1                                               | Almost all values of ˜ α . . . . . . . . . . . . . . . . .    |   74 |
|     | D.2                                               | An identity for interacting functionals . . . . . . . . .     |   77 |
|     |                                                   | D.2.1 How many times can the Casimir be applied?              |   77 |
|     |                                                   | D.2.2 Deriving the identity . . . . . . . . . . . . . .       |   77 |
|     | D.3                                               | All values of ˜ α . . . . . . . . . . . . . . . . . . . . . . |   78 |

## 1 Introduction

A special class of d -dimensional conformal field theories (CFTs) admit a dual description as local theories in a ( d +1)-dimensional bulk Anti de Sitter (AdS) spacetime. Although the main case of interest is when this theory has a dynamical metric, much can be learned by considering the simpler case of UV-complete QFTs in a non-fluctuating AdS space. Such theories support well-defined, exactly localised operators which are mutually local, i.e. commuting at spacelike separations (unlike the gravitational case). Expressing these local bulk fields in terms of the boundary CFT data is an old problem studied by several authors [1-14].

In our previous work [15], we considered this problem from a modern bootstrap perspective. Our starting point was the state-operator correspondence in AdS, which guarantees that these operators admit a representation in terms of boundary CFT operators, schematically

$$\Psi ^ { \text{AdS} } \sim \sum _ { \Delta } \mu _ { \Delta } \, \mathcal { O } _ { \Delta } ^ { \text{CFT} } \.$$

The task of reconstructing a bulk field Ψ AdS means producing an admissible choice of the dynamical information in this boundary operator expansion (BOE): the coefficients µ ∆ . To

see that most choices of µ ∆ would not be admissible, one may inspect the correlator of the putative bulk field against two boundary operators O 1 , O 2 . After stripping off universal kinematic prefactors, the resulting 2-point AdS form factor F z ( ) depends on a single conformal cross-ratio, denoted z , and admits a BOE expansion:

$$\langle \Psi \, | \, \mathcal { O } _ { 1 } \, \mathcal { O } _ { 2 } \rangle \sim F ( z ) = \sum _ { \Delta } \mu _ { \Delta } \lambda _ { \Delta } ^ { 1 2 } \, G _ { \Delta } ( z ) \,,$$

with further details given below.

Microcausality in the bulk implies that correlation functions can only become singular when local insertions enter each other's lightcones. In the present setup, this only happens for z ≤ 0, so the only singularity of F should be there. On the other hand, the BOE allows the existence of another, unphysical branch cut at z ≥ 1 (at Lorentzian but spacelikeseperated kinematics), which is present in each summand appearing above. Our locality constraint is the vanishing of this acausal branch cut [12]:

$$\mathcal { I } _ { z \geq 1 } F ( z ) = 0 \.$$

Since each block G ∆ has a discontinuity at z ≥ 1, this clearly places constraints on the coefficients µ ∆ , λ ∆ . However, characterizing these constraints is not immediately easy since the discontinuity appearing in this equation does not commute with the BOE expansion.

We resolved this issue in [15] by smearing the form factor against specially designed analytic functional kernels , along a contour that wraps the hypothetical z ≥ 1 cut. After this smearing, one can commute the locality condition with the BOE to produce sum rules on the BOE data. Concretely, we constructed a basis of kernels { f n } that lead to sum rules satisfying two desirable properties:

## Property 1): Completeness

$$\mathcal { I } _ { z \geq 1 } F ( z ) = 0 \quad \Longleftrightarrow \quad \sum _ { \Delta } \mu _ { \Delta } \lambda _ { \Delta } ^ { 1 2 } \ \theta _ { n } ( \Delta ) = 0 \quad \text{for } \ n \geq 1 \.$$

## Property 2): Duality

For some { ∆ n } ∞ n =1 , θ n (∆ m ) = δ nm for m,n ≥ 1 .

Here θ n (∆) denotes the action of a functional f n on the block G ∆ . Let us emphasise that these conditions draw up a duality between a basis of functionals and some distinguished spectrum,

$$\{ f _ { n } \} \quad \longleftrightarrow \quad \{ \Delta _ { n } \} \.$$

In [15], we constructed functionals satisfying these properties, dual to the specific choice of a Generalized Free Field-type spectrum { ∆ = 2˜ + 2 n α n } . In this paper, we will do it for a very general class of interacting spectra.

Before detailing our new results, let us explain the significance of the above properties. Completeness is clearly desirable, guaranteeing that the sum rules fully capture the locality constraint. However, it is still a difficult task to actually solve the infinite set of sum rules for the BOE data. Duality strikes at the core of this issue: it means the sum rules are trivialised when the BOE spectrum { ∆ is taken to contain the particular distinguished one } { ∆ n } that is 'dual' to the basis of functionals. As a result, the sum rules, and hence the locality problem, admit certain canonical solutions. For example, making the ansatz for the BOE spectrum { ∆ } = { ∆ n } ∪ { ∆ 0 } consisting of the dual spectrum, plus one extra block of dimension ∆ , the sum rules are 0 diagonalised and admit an exact solution for the BOE data c n ≡ µ ∆ n λ 12 ∆ n :

$$\mathcal { I } _ { z \geq 1 } \left [ G _ { \Delta _ { 0 } } + \sum _ { n = 1 } ^ { \infty } c _ { n } \, G _ { \Delta _ { n } } \right ] = 0 \quad \underset { \theta _ { n } ( \Delta _ { m } ) = \delta _ { n m } } { \text{Quality} } \quad c _ { n } = - \theta _ { n } ( \Delta _ { 0 } ) \.$$

For the choice ∆ n = 2˜ + 2 α n of [15], the form factors corresponding to these solutions were called 'local blocks', and have a simple interpretation in terms of weakly coupled QFTs in AdS space. For instance, with ˜ = α ∆ 0 2 = ∆ , they correspond to the form factors ϕ ⟨ Φ 2 | ϕϕ ⟩ in the theory of a free scalar Φ in AdS with mass m 2 = ∆ (∆ ϕ ϕ -d ).

We will refer to these special solutions as extremal solutions . These solutions share many features with the familiar extremal solutions that saturate unitarity bounds in other problems, such as the crossing equation [16-20]. Namely, (i) they are 'minimal' solutions, in the sense that no blocks can be removed, (ii) they possess dual bases of functionals and (iii) they can be explicitly bootstrapped using those functionals. Indeed, understanding extremal solutions of bootstrap problems is an important motivation for us: it is equivalent to understanding the bootstrap bounds that they saturate, and would open the door for their analytic understanding. In the case of the 1d crossing equation for 4-point correlators, free extremal solutions have been studied analytically [17-19], while interacting ones have only been studied numerically [20]. For this simpler locality problem, we can do better: for a very large set of choices of interacting spectra { ∆ n } , we will analytically construct sets of sum rules satisfying properties 1) and 2) and the corresponding extremal solutions (1.5).

Our main results concern the locality problem for d &lt; 4 and are stated as follows:

Assumptions. Take a ' sparse ' spectrum, ∆ = 2˜ + 2 n α n + γ n ( n ∈ N &gt; 0 ), meaning:

- · γ n = n →∞ O n ( -ϵ ) , ϵ &gt; 0 ,
- · γ n is an analytic function at large n ∈ C (for small arg( n )).

## Results.

- · Then there exist functionals θ n satisfying properties 1) and 2),
- · The functionals form a Schauder basis for a certain topological vector space,
- · The functionals' actions on blocks, θ n (∆), can be written down explicitly:

̸

$$\theta _ { n } ( \Delta ) = \prod _ { \substack { m = 1 \\ m \neq n } } ^ { \infty } \left ( \frac { \lambda _ { \Delta } - \lambda _ { \Delta _ { m } } } { \lambda _ { \Delta _ { n } } - \lambda _ { \Delta _ { m } } } \right ) \,, \quad \lambda _ { \Delta } = \Delta ( \Delta - d ) \,.$$

Let us now comment on the relevance of these results.

Firstly, as discussed above, they immediately produce huge classes of interacting, highly non-trivial solutions to the locality problem. The most elementary such solution has the BOE coefficients given by equation (1.5). Using the third result above, the resulting local form factor takes the explicit form

̸

$$F _ { \Delta _ { 0 } } ^ { \{ \Delta _ { n } \} } = G _ { \Delta _ { 0 } } - \sum _ { n = 1 } ^ { \infty } \left [ \, \prod _ { \substack { m = 1 \\ m \neq n } } ^ { \infty } \left ( \frac { \lambda _ { \Delta _ { 0 } } - \lambda _ { \Delta _ { m } } } { \lambda _ { \Delta _ { n } } - \lambda _ { \Delta _ { m } } } \right ) \, \right ] G _ { \Delta _ { n } } \.$$

These non-perturbative, interacting solutions are valid for any desired choice of { ∆ n } satisfying the assumptions above. The products and sums are absolutely convergent and are readily evaluated numerically: for an example of this technology in action, the reader may look ahead to plots of these form factors in Figure 6.

Secondly, for general d , these results can be taken as another piece of evidence that Generalized Free Fields and their deformations arise from local AdS Lagrangians [21]. Indeed, for a CFT correlator whose (scalar) spectrum is close to GFF, and without additional operators, our results allow us to readily construct associated solutions to the locality problem in AdS. Moreover, for the d = 1 case (corresponding to AdS /CFT ), our results 2 1 indicate that arbitrary 'extremal' solutions ⟨ ϕϕϕϕ ⟩ of the 1d crossing equation admit local AdS fields Ψ, such that the form factors ⟨ Ψ | ϕϕ ⟩ are manifestly local and computable. This is because those extremal solutions (saturating unitarity bounds) have sparse spectra of the type assumed above [20]. 1 This means that, in order to improve numerical bootstrap bounds, one must do more than simply imposing locality constraints on top of the standard

1 At least with respect to sparseness. Analyticity properties also hold in perturbation theory and in certain (flat space) limits, but have not be established in full generality.

semi-definite programming for the crossing equation. We will discuss this matter further in Section 6.

Finally, while the second result in the box above sounds abstract and technical, we would like to emphasise the conceptual importance of the functionals forming a complete basis. For the first time, we will make precise the notion of completeness of functionals for the bootstrap; to which space they belong; and the deep connection between complete bases of functionals and extremal solutions of the bootstrap (i.e. solutions uniquely determined by bootstrap equations, with sparsest possible OPEs, and saturating bounds). Our results make concrete a correspondence:

Extremal solutions ←→ Bases of a topological vector space

Establishing the existence of the interacting bases will lead us through some unfamiliar but rich mathematics. While our construction will be detailed for the locality problem, we will argue that the same logic is also likely to apply to the crossing problem for CFTs. More generally, we believe that the formalism and constructions of this work are the right mathematical language for thinking about analytic functionals - which have found important applications not only in physical, but also mathematical, problems (see [22] and refs. therein).

The remainder of this introduction is a bird's eye view of the ideas leading to the above results and an outline of the paper's structure.

## Key ideas of the paper

Our first task is to understand carefully what it means for a set of functionals to be complete. A familiar notion of completeness arises in Hilbert spaces, which possess orthonormal bases, with respect to which arbitrary elements can be decomposed. Given such an orthonormal basis | e n ⟩ , an element is shown to be zero by checking that its expansion coefficients in that basis are all zero,

$$| F \rangle = 0 \quad \Leftrightarrow \quad \theta _ { e _ { n } } [ F ] \coloneqq \langle e _ { n } | F \rangle = 0 \quad \forall n \.$$

We emphasise that e n has the interpretation as an element both of the Hilbert space and of its dual, i.e. as a bounded linear functional acting on the Hilbert space. Can we apply the same logic to our problem? We need to check whether a certain discontinuity of an analytic function vanishes (eq. (1.3)). As mentioned above, we will obtain our sum rules by integrating the unphysical cut of the form factor against a functional kernel:

$$\theta _ { n } [ F ] = \int _ { 1 } ^ { \infty } \mathrm d z \, f _ { n } ( z ) \, \mathcal { I } _ { z } F \.$$

This expression would suggest defining some suitable Hilbert space of functions on [1 , ∞ ] for which f n would form a basis. Making this idea concrete will require overcoming several

difficulties and the introducing some new concepts.

The first such difficulty deals with understanding what kind of object is I z F and in what space it lives. We propose that it is both natural and convenient to view this object as a hyperfunction : the discontinuity of an analytic function across a branch cut. 2 Hyperfunctions are a very large class of distributions: they contain derivatives of delta-functions of arbitrary orders, but also non-Schwartzian distributions such as the discontinuity of an essential singularity,

$$\mathcal { I } _ { z } \, \frac { 1 } { z ^ { n } } \sim \partial _ { z } ^ { n } \, \delta ( z ) \, \quad \mathcal { I } _ { z } \, e ^ { 1 / z } \, \sim \, ^ { \omega } \partial _ { z } ^ { \infty } \, \delta ( z )$$

The locality condition is then simply the requirement that it I z ≥ 1 F is the zero hyperfunction. Choosing to work in this space means that we can rely on elegant tools of complex analysis to understand the constraints of locality. It also makes only minimal assumptions about the properties of form factors (essentially absolute convergence of the BOE), meaning the scope of our work extends to cases where we may in future want to allow non-trivial discontinuities (say when reconstructing non local bulk operators). 3

We are thus led to consider a space of hyperfunctions HF with a branch cut along an interval V . This is a topological vector space, so it has a notion of convergence, but does not have an inner-product or norm, and in particular it is not a Hilbert space. However, it does have a dual space , which is the space A of real-analytic functions on V . There is then a dual pairing between f ∈ A and F ∈ HF :

$$\theta _ { f } [ F ] = ( f, F ) \coloneqq \oint _ { \Gamma _ { V } } \frac { \mathrm d z } { 2 \pi i } \ f ( z ) \, F ( z )$$

$$F ) \coloneqq \oint _ { \Gamma _ { V } } \frac { \aleph \epsilon } { 2 \pi i } \ f ( z ) \, F ( z )$$

where Γ V is a contour encircling the interval V , inside the shared domain of analyticity of f and F . As we will discuss, defining the functionals in this way makes the swapping property of [24] automatic, since the pairing ( · · , ) is continuous.

The space A is again not a Hilbert space, but it still makes sense to ask for complete sets of analytic functions f n . This is a set such that any f ∈ A is a limit of linear combinations of the f n . Integrating against such a complete set then provides a characterisation of whether a hyperfunction vanishes (generalising (1.7)):

$$\mathcal { I } _ { z } F = 0 \ \in \ \mathcal { H F } \quad \Leftrightarrow \quad ( f, F ) = 0 \ \forall f \in \mathcal { A } \quad \Leftrightarrow \quad ( f _ { n }, F ) = 0 \ \forall n \. \ \ ( 1. 1 1 )$$

Having understood the proper meaning of completeness, i.e. property 1), we are still left with the task of constructing not only complete sets of f n , but ones that are dual to particular

2 Formally, hyperfunctions on an interval V are defined as functions analytic on a complex neighbourhood of V but singular at V , up to equivalence by adding functions that are non-singular at V . In this sense, the hyperfunction 'only cares' about the discontinuity of F .

3 For CFT 4-point correlators, it was argued in [23] that the discontinuities are tempered distributions, with the dual space of linear functionals being the Schwartz space of test functions. We comment further on the relation of their work with ours in footnote 12 and Section 3.4.1.

spectra of blocks, i.e. satisfying property 2). We have found a natural way to obtain both these properties. The key idea is that, while I z F is taken to be a hyperfunction, it may decomposed into the discontinuities I z G ∆ of individual blocks thanks to the BOE (1.1). Each of these is actually a function, and they are naturally associated to a Hilbert space: L 2 functions on [1 , ∞ ). This L 2 space has a continuous embedding as a dense subspace of HF . In turn, the analytic functional kernels are also naturally identified with a dense subspace of L 2 . Thus we naturally arrive at a hierarchy of spaces called a Gelfand triple (or rigged Hilbert space )

$$\mathcal { A } = \left \{ _ { \text{functional kernels} } \text{real-analytic} \right \} \, \subset \, L ^ { 2 } \, \subset \, \mathcal { H F } = \left \{ \text{hyperfunctions} \right \} \,.$$

The L 2 space will buy us technical power: it is a Hilbert space with well-behaved notions of basis and decomposition. In particular, one can already guess that a useful approach may be first to construct bases of L 2 consisting of analytic functions, and then argue that they form the desired complete sets in A , leading to property 1).

To achieve property 2), we have to introduce one last idea: we should work not with orthonormal bases for L 2 , but with the more general Riesz bases [25]. A Riesz basis can be seen as a structure-preserving 'skewing' of an orthonormal basis: it is the image of an orthonormal basis under a bounded, invertible, linear map. While Riesz basis elements { p n } are not necessarily orthonormal, they automatically possess another, dual Riesz basis { h n } such that

$$\langle h _ { n } \, | \, p _ { m } \rangle = \delta _ { n m }$$

This suggests choosing sets ∆ n such that the discontinuities p n = I z G ∆ n form a Riesz basis, and defining the functional kernels to be the dual Riesz basis, f n = h n . In this way we have the desired duality:

$$\theta _ { f _ { n } } ( \Delta _ { m } ) = \langle f _ { n } \, | \, \mathcal { I } _ { z } G _ { \Delta _ { m } } \rangle = \delta _ { n m }$$

and the completion of our program.

As well as understanding which choices of ∆ m form Riesz bases for L 2 , we must establish which of those can be mapped back to the spaces A and HF . In general this is a complicated mathematical problem. Our strategy will be first to reformulate the functionals dual to the 'free' spectrum ∆ n = 2˜ + 2 α n (constructed in our previous work [15]) in this language, showing that they satisfy all the required properties to give a complete set of functional kernels in A . After that, we will show that the deformation ∆ n → 2˜ + 2 α n + γ n , with γ n satisfying the assumptions above, does not modify this conclusion. In particular the decay of γ n is necessary to ensure that we still have a Riesz basis, while the fact that the γ n admits an analytic continuation for complex n is required to establish that the dual Riesz basis elements h n (which are a priori just L 2 functions) are actually analytic.

Finally, we will study the space PW of functional actions θ (∆), defined 4 as

$$\mathcal { P W } = \left \{ \begin{array} { c } \theta \, \colon \mathbb { C } \longrightarrow \mathbb { C } \\ \theta ( \Delta ) = \langle f \, | \, \mathcal { I } _ { z } G _ { \Delta } \rangle & \text{for } f \in \mathcal { A } \end{array} \right \} \,.$$

These functions of the complex variable ∆ are physically very significant: they encode all of the valid sum rules ∑ ∆ µ ∆ λ 12 ∆ θ (∆) = 0 characterising the locality problem. By studying the transform ⟨ f | I z G ∆ ⟩ , we will show that PW is also a dense subspace of a Hilbert space, for which the θ f n (∆) provide a (Schauder) basis. Most importantly, we provide an independent definition of PW as a space of entire functions with prescribed behaviour near infinity. Such entire functions given by an explicit product of their zeros, leading to the third result above.

Outline. The contents of the paper are organised as follows.

- · In Section 2, we set out the problem of finding complete sets of sum rules for locality in terms of the space of hyperfunctions and its dual space of analytic functionals.
- · Section 3 explains how it is possible to construct Schauder bases for these spaces by using an auxiliary Hilbert space that sits between them. The outcome is a large class of complete sets of interacting functionals and associated sum rules.
- · In Section 4, we consider the space of functional actions, and show that it inherits the same Hilbert space structure of the space of functional kernels. We give an independent characterization of this space as a set of entire functions with certain exponential type. With this definition, we write down explicit expressions for the functional actions forming our interacting bases.
- · In Section 5, we verify our construction in an explicit example. First we set out the general relationship between the large-∆ behaviour of the functional actions and the asymptotics of the functional kernels. Then we start from a known family of extremal 4-point correlators solving the 1d crossing equation. By constructing a basis of interacting functionals dual to the spectrum appearing in these correlators' OPEs, we reconstruct local form factors for bulk operators consistent with these correlators.
- · In Section 6, we discuss the implications of our results and future directions of study.
- · Several technical appendices complement the paper. Appendix A reviews the results of our paper [15] and some useful asymptotic formulae. The most non-trivial technical work is: proving in Appendix B that our sparse sets of discontinuities of blocks form Riesz bases; and proving in Appendix C that the dual functional kernels are realanalytic functions of the appropriate class. Appendix D addresses the case of the most general values of ˜ and the finite subtraction procedure relevant for that case. α Appendix E is a technical complement to Section 5.

This subject being rather mathematical and technical, we have attempted to tread a fine line (barring inevitable missteps) between rigour and accessibility - with a physicist's penchant for the latter.

4 The definition given here is, for simplicitly, somewhat schematic. See Section 4 for precise definitions.

## 2 Locality bootstrap and hyperfunctions

## 2.1 Form factors as hyperfunctions

Let us review the properties of form factors and show there is a natural mapping of these objects into a certain space of hyperfunctions. Crucially, this mapping will be defined in a manner compatible with the BOE decomposition (1.1). As in [15], we are considering form factors ⟨ Ψ | O 1 O ⟩ 2 for bulk scalar fields Ψ, so that all the primary operators O ∆ appearing in the BOE are scalars. In this paper, we will set:

$$\Delta _ { 1 } = \Delta _ { 2 } \quad \text{ and } \quad d < 4$$

with d the spacetime dimension of the boundary CFT. These assumptions simplify some technicalities, roughly corresponding to avoiding needing certain subtractions; but we do not expect them to be decisive for our general approach.

We will call form factor any function F z ( ) satisfying certain properties described below, which imply in particular that F z ( ) is analytic on the cut plane C \ (( -∞ , 0] ∪ [1 , ∞ )). A local form factor is one that is actually analytic on C \ -∞ ( , 0], i.e. with no discontinuity at [1 , ∞ ) - solving the locality constraint (1.3).

The first property is that F z ( ) should admit a BOE expansion, 5

$$F ( z ) = \sum _ { \Delta } \mu _ { \Delta } \lambda _ { \Delta } ^ { 1 2 } G _ { \Delta } ( z ) = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } ( z )$$

where, for convenience, we have changed normalisation by setting

$$c _ { \Delta } = \mathcal { N } _ { \Delta } ^ { - 1 } \ \mu _ { \Delta } \, \lambda _ { \Delta } ^ { 1 2 } \, \quad \ \ g _ { \Delta } ( z ) = \mathcal { N } _ { \Delta } \ G _ { \Delta } ( z )$$

with

$$\mathcal { N } _ { \Delta } \coloneqq \frac { \Gamma ( \frac { \Delta } { 2 } ) ^ { 2 } } { \pi \, \Gamma ( \Delta + 1 - \frac { d } { 2 } ) } \,, \quad \ \ G _ { \Delta } ( z ) = \, z ^ { \frac { \Delta } { 2 } } \, _ { 2 } F _ { 1 } ( \frac { \Delta } { 2 }, \frac { \Delta } { 2 } ; \Delta + 1 - \frac { d } { 2 } ; z ) \. \quad \ ( 2. 4 )$$

This is an important physical consequence of the state-operator correspondence in AdS (see discussion in [15] and refs. there).

The second property is to do with the sense in which the BOE converges, and it implies that F z ( ) satisfies certain boundedness properties near the branch points at z = 1 and ∞ . Concretely, we will ask that

$$| F ( z ) | \leq \sum _ { \Delta } | c _ { \Delta } | | g _ { \Delta } ( z ) | \sim \begin{cases} \begin{array} { c c } | z | ^ { \alpha _ { F } } & \text{as $z\to\infty$} \\ | z - 1 | ^ { - \alpha _ { F } } & \text{as $z\to1$} \end{array} \end{cases}$$

5 The sum over states is implicitly restricted to unitary ∆ ∈ { 0 } ∪ [ d -2 2 , ∞ ).

for some fixed number α F . 6 This specific polynomial bound is physically motivated by considering the bulk QFT's UV completion: see [26]. Note that these constraints imply that F z ( ) is analytic on C \ ( ( -∞ , 0] ∪ [1 , ∞ ) ) .

Let us now examine the discontinuities of these functions on the cut [1 , ∞ ). It is natural to think of these discontinuities as elements of a space of hyperfunctions . A hyperfunction on an interval is defined to be an equivalence class of functions which are analytic in a complex neighbourhood of the interval and which share the same discontinuity (see, e.g., the textbook [27]). Clearly, given one representative in the equivalence class, an equally good one is obtained by adding any analytic function. Choosing once and for all a representative, let us define the subtracted form factor ,

$$F ( z ) \quad \rightarrow \quad F ^ { ( \tilde { \alpha } ) } ( z ) = \oint _ { \Gamma } \frac { \mathrm d w } { 2 \pi i } \, \frac { F ( w ) \, w ^ { - 1 - \tilde { \alpha } } } { w - z } \, \quad \ z \in \mathbb { C } \, \sum [ 1, \infty )$$

where Γ is a contour wrapping the unphysical cut [1 , ∞ ), as shown in Figure 1. F (˜) α is an analytic function on C \ [1 , ∞ ) that vanishes if and only if the original form factor is local

$$F ^ { ( \tilde { \alpha } ) } \equiv 0 \quad \Longleftrightarrow \quad \mathcal { I } _ { z \geq 1 } F = 0 \.$$

This is because its only singularities are at [1 , ∞ ), given by

$$\mathcal { I } _ { z } F ^ { ( \tilde { \alpha } ) } = z ^ { - 1 - \tilde { \alpha } } \mathcal { I } _ { z } F \, \quad 1 \leq z < \infty \.$$

In order to produce this hyperfunction capturing only the unphysical cut, we had to integrate in between the two cuts [1 , ∞ ) and ( -∞ , 0] that touch at ∞ (see Figure 1). To make the integral (2.6) well-defined, we required a subtraction by a large enough power of w :

$$\tilde { \alpha } > \alpha _ { F } - 1 \.$$

We note that the largez behaviour of the subtracted form factor is now:

$$| F ^ { ( \tilde { \alpha } ) } ( z ) | \stackrel { z \to \infty } { \sim } | z | ^ { - 1 - \tilde { \alpha } + \alpha _ { F } } = | z | ^ { - \epsilon } \, \quad \epsilon > 0$$

whereas behaviour near z = 1 is the same as F :

$$| F ^ { ( \tilde { \alpha } ) } ( z ) | \stackrel { z \to 1 } { \sim } | z - 1 | ^ { - \alpha _ { F } } = | z - 1 | ^ { - \tilde { \alpha } - 1 + \epsilon } \, \quad \epsilon > 0 \.$$

These considerations motivate the following definition of a space H F of hyperfunctions

6 In fact for our purposes a weaker assumption would be sufficient: see equation (3.39) and the discussion in Section 3.4.1.

Figure 1: The integration contour Γ used to define the subtracted form factor F (˜) α in equation (2.6). The contour goes through ∞ , passing between the two cuts meeting there (the physical one in blue and the unphysical one in red).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-52479b60-dbb2-4c1b-a76c-1fe8ec109901.jpg)

with properties similar to those of F (˜) 7 α :

$$\mathcal { H F } = \left \{ H ( z ) \text{ analytic on } \mathbb { C } \, \langle \, [ 1, \infty ) \, \left | \, \begin{array} { c } | H ( z ) | ^ { z \to \infty } \end{array} | z | ^ { - \epsilon }, \, | H ( z ) | ^ { z \to 1 } \, | z - 1 | ^ { - \tilde { a } - 1 + \epsilon } \right \}. \ ( 2. 1 2 )$$

We have seen that for any form factor F , its subtraction F (˜) α is an element of this space (provided ˜ α &gt; α F -1). Our task will be to understand under which circumstances F (˜) α coincides with the zero element of H F , since this is precisely the locality condition (2.7) for the form factor.

An important property of the map from form factors to HF is that it is nicely compatible with the BOE. In particular, the BOE descends to a convergent decomposition of the subtracted form factor in HF :

$$F ^ { ( \tilde { \alpha } ) } = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } ^ { ( \tilde { \alpha } ) }$$

where g (˜) α ∆ are the hyperfunctions obtained by applying the subtraction (2.6) to g ∆ . The above sum converges in the space of hyperfunctions with its natural compact-open topology - meaning it converges uniformly on compact regions away from the cut [1 , ∞ ). 8

7 It may seem odd that the space HF depends on the parameter ˜, which may be chosen freely as long α as ˜ α &gt; α F -1. This subtlety is addressed below in Section 3.4.1.

8 Formally, because of the bounds at the endpoints of the cut, we are using a slightly finer topology than the standard compact-open one: we add some extra open sets so that H n → 0 in HF also implies ∫ ∞ | d z || z | -1 | H n ( z ) | → n →∞ 0 and ∫ 1 | d z || z - | | 1 ˜ α H n ( z ) | → n →∞ 0 on arbitrary integration contours touching the endpoints ∞ and 1. In the same vein, one may understand the definition of the space (2.12) more formally as requiring ∫ ∞ | d z || z | -1 | H z ( ) | &lt; ∞ and ∫ 1 | d z || z - | | 1 ˜ α H n ( z ) | &lt; ∞ . One can check that the space defined in this way is complete. The dual space A defined in (2.17) may also be formalised in a compatible way.

To see that (2.13) converges in this way, we will use the fact that a pointwise limit of analytic functions converges uniformly on compact sets inside the shared domain of analyticity of the functions and their limit. 9 Hence it is sufficient to show that (2.13) converges pointwise away from [1 , ∞ ). To prove this, let us apply the subtraction transform (2.6) to the BOE (2.2): if we can commute the integral with the sum over ∆, then we are done. This is indeed possible by splitting the integral into two regions: one near ∞ and one elsewhere covered by a compact domain. For the region near ∞ , the bound in (2.5) allows to swap the sum and integral by Dominated Convergence. On the other region covered by a compact domain, the BOE (2.2) converges uniformly by the fact mentioned above, so the sum and integral can be swapped.

Before moving on let us make a quick observation. The formulation (2.7) of locality as the vanishing of F (˜) α is closely related to the dispersion relation of [15] since we can write

$$F ^ { ( \tilde { \alpha } ) } ( z ) = z ^ { - 1 - \tilde { \alpha } } \left ( F ( z ) - \oint _ { \Gamma _ { [ - \infty, 0 ] } } \frac { \mathrm d w } { 2 \pi i } \frac { ( z / w ) ^ { 1 + \tilde { \alpha } } } { w - z } \, F ( w ) \right ) \,.$$

where the dispersion relation is the vanishing of the RHS. Applying this relation to a single block gives a relation between subtracted blocks g (˜) α ∆ in (2.13) and the 'local blocks' of [15]:

$$g _ { \Delta } ^ { ( \tilde { \alpha } ) } ( z ) = \mathcal { N } _ { \Delta } \, z ^ { - 1 - \tilde { \alpha } } \left ( G _ { \Delta } ( z ) - \mathcal { L } _ { \Delta } ^ { \tilde { \alpha } } ( z ) \right ) \,.$$

## 2.2 Functionals and locality conditions

Now that we have mapped form factors into hyperfunctions, we will describe how such objects can be probed by the action of linear functionals. As discussed in the introduction, HF is a very broad class of distributions. Like for other classes of distributions, they may be defined as the continuous linear functionals acting on a space of test functions: 10

$$\mathcal { H F } = \mathcal { A } ^ { \prime }$$

where the test functions A are the real-analytic functions on the interval. In our case, the hyperfunctions have particular bounded behaviour at the endpoints 1 and ∞ (cf. (2.10)

9 A precise statement (a result of Osgood [28]; see the review in [29]) is that, if g ∆ are holomorphic on a certain domain and the series F = ∑ ∆ c ∆ ∆ g converges pointwise, then F is holomorphic on a dense, open subset Ω of the domain, and the convergence is uniform in compact subsets of Ω. If F is holomorphic on (at least) the same domain as g ∆ , then Ω is the whole domain.

10 The identification of hyperfunctions as the dual to real-analytic functions on a compact real domain is due to the Grothendieck-K¨the-Silva duality. o The compact-open topology that we use on HF is then the natural strong-dual topology inherited from this duality.

Figure 2: The integration contour (in black) defining the action of functionals on subtracted form factors. The shared domain of analyticity is shown in blue, the branch cut of the form factor in red and possible non-analyticities of the functional kernel in green.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-1003162c-ff9a-4efb-8a81-5efb2ea11929.jpg)

and (2.11)). Hence the test function space can be taken to be:

$$\mathcal { A } = \left \{ f ( z ) \, \text{real-analytic on } ( 1, \infty ) \, \left | \begin{array} { c } \text{analytic on finite angle at endpoints} \\ \left | f ( z ) \right | \underset { z \rightarrow 1 } { \sim } \left | z \right | ^ { - 1 } \\ f ( z ) \underset { z \rightarrow 1 } { \sim } \left ( \text{analytic} \right ) + \left | z - 1 \right | ^ { \tilde { \alpha } } \end{array} \right \} \,. \quad ( 2. 1 7 )$$

A typical domain of analyticity for an element of A is shown in blue in Figure 2. For example, f ( z ) may have branch cuts running along ( -∞ , 1), as well as other singularities away from, but arbitrarily close to, the interval [1 , ∞ ].

The definition (2.16) means there is a continuous dual pairing between A and HF , defined by

$$( f, G ) \coloneqq \oint _ { \Gamma _ { [ 1, \infty ] } } \frac { \mathrm d z } { 2 \pi i } \, f ( z ) \, G ( z ) \, \quad \ \ f \in \mathcal { A } \, \ \ G \in \mathcal { H F }$$

where Γ [1 , ∞ ] is a contour wrapping the cut [1 , ∞ ], inside the shared domain of analyticity of F and f . Since f is permitted non-analyticity at the end points, the contour must in general be 'pinched' to go through the points 1 and ∞ , approaching these endpoints at a finite angle from the cut [1 , ∞ ] (see Figure 2). 11

11 The analytic term allowed in the functionals at z = 1 is understood as acting on a contour pushed to the left at z = 1 to avoid the singularity of the hyperfunction.

We thus identify A as a natural class of functional kernels acting on HF as 12

$$\theta _ { f } [ H ] \coloneqq ( f, H ) \quad \ \ f \in \mathcal { A } \, \ \ H \in \mathcal { H F } \.$$

↦

Since the pairing ( · · , ) is continuous, then these functionals are automatically continuous. Formally, what is happening here is that there is a canonical embedding A → A ′′ = HF ′ , which maps f → θ f . In this natural mathematical language, the swapping property of [24] becomes automatic: the functionals (2.19) are compatible with the BOE:

$$\theta _ { f } \left [ F ^ { ( \tilde { \alpha } ) } \right ] = \sum _ { \Delta } c _ { \Delta } \, \theta _ { f } \left [ g _ { \Delta } ^ { ( \tilde { \alpha } ) } \right ] \,.$$

This follows because the functionals are continuous and the BOE (2.13) converges in the topology of HF .

Let us summarise what we have found so far. Firstly, we described a mapping of form factors into the space H F , such that

$$F = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } \ \text{local} \quad \Longleftrightarrow \ \ F ^ { ( \tilde { \alpha } ) } = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } ^ { ( \tilde { \alpha } ) } = 0 \ \in \mathcal { H F }$$

Secondly, we showed there is a dual space of objects which act on hyperfunctions continuously. These objects are analytic functions on the interval (1 , ∞ ) with certain boundedness properties at the endpoints. This means in particular that

$$H = 0 \ \in \mathcal { H F } \ \iff \ \theta _ { f } [ H ] = 0 \ \text{ for all } f \in \mathcal { A } \.$$

Combining both of these, we see that:

$$F = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } \text{ local} \quad \Longleftrightarrow \quad \sum _ { \Delta } c _ { \Delta } \, \theta _ { f } \left [ g _ { \Delta } ^ { ( \tilde { \alpha } ) } \right ] = 0 \quad \text{for all } \, f \in \mathcal { A } \quad \quad ( 2. 2 3 )$$

At this point let us recall our two initial goals defined in the introduction: to find a complete set of sum rules for locality satisfying duality with respect to a particular spectrum. In the language of hyperfunctions, this becomes the following task:

Task: Produce a set of functional kernels { f n } ∞ n =1 in A such that:

- 1) { f n } are complete in A .
- 2) f n are dual to a set of blocks, ( f m , g (˜) α ∆ n ) = δ mn , for some { ∆ n } . .

12 The analytic functionals of the type used in [17,19,30] clearly fall into the present class of functionals A . Hence they are also in the larger class of smooth functional kernels considered in [23]. Note that the standard derivative functionals are also in our class: ∂ n z H ( 1 2 ) = D H n [ ] := ( n !( z -1 2 ) - -n 1 , H ).

The word 'complete' now has precise meaning: { f n } is a complete set in the topological vector space A if any element of A is a limit of linear combinations of the f n . 13 Together with (2.23), this implies that the associated sum rules are complete in the sense of Property 1) from the introduction:

$$F = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } \text{ local} \quad \Longleftrightarrow \quad \sum _ { \Delta } c _ { \Delta } \, \theta _ { n } ( \Delta ) = 0 \quad \text{for all } n \geq 1$$

where we have defined

$$\theta _ { n } ( \Delta ) \coloneqq \theta _ { f _ { n } } \left [ g _ { \Delta } ^ { ( \tilde { \alpha } ) } \right ] = ( f _ { n }, \, g _ { \Delta } ^ { ( \tilde { \alpha } ) } ) \.$$

In particular, the second element of this task allows us then to achieve Property 2) as stated in the introduciton.

## 2.3 Generalized Free Field functionals

Let us now rephrase the results of our previous work [15] in the language of the preceding discussion. In that work we produced a set of analytic functions ˆ ( f n z ), given explicitly in Appendix A. They are simple polynomials of degree n in 1 /z without constant term, and thus clearly belong to A . They were found by demanding that they satisfy the duality property with respect to the sparse spectrum

$$\hat { \Delta } _ { n } = 2 \tilde { \alpha } + 2 n$$

i.e. a Generalized Free Field (GFF) spectrum of double-traces (with a hat denoting GFF quantities). To see that { ˆ ( f n z ) } are complete is, in this case, rather simple, since one can just appeal to the fact that, for any value of ˜, α { zf ˆ ( n z ) } are a complete set of polynomials (one of each order). It is a mathematical fact that such functions are dense in the space of real analyic functions, and it then follows that { f n } are dense in A .

Suppose however that we were not aware of this mathematical fact. How can we convince ourselves that the GFF functional kernels are complete? We will now set out a proof of this property that will generalise to the case of functionals dual to interacting spectra ∆ n -whose kernels will no longer be polynomials. The idea is to trade the completeness condition in our task for something equivalent:

13 Let us make a technical comment about the topology of A . There are several natural topologies on this space, but since we are only using A as functionals acting on HF , it suffices for us to equip it with the weak topology . This means that convergence f n → f in A just means that the action on any hyperfunction converges: ( f n , H ) → ( f, H ) for any H ∈ HF .

Task': Produce a set of functionals { f n } ∞ n =1 in A such that:

- 1') Any H ∈ HF decomposes as H = ∑ n ( f , H n ) g (˜) α ∆ n ,
- 2) f n are dual to a set of blocks, ( f m , g (˜) α ∆ n ) = δ mn , for some { ∆ n } .

This immediately implies the completeness condition 1) above, as follows: for any f ∈ A , acting with it on the equation in 1'), we have

$$( f, H ) = \sum _ { n } ( f _ { n }, H ) \, ( f, g _ { n } ^ { ( \tilde { \alpha } ) } ) = \left ( \sum _ { n } ( f, g _ { n } ^ { ( \tilde { \alpha } ) } ) f _ { n }, \, H \right ) \quad \forall H \in \mathcal { H F }, \, f \in \mathcal { A } \. \quad \ ( 2. 2 7 )$$

This means that (in the weak topology)

$$f = \sum _ { n } ( f, g _ { n } ^ { ( \tilde { \alpha } ) } ) \, f _ { n } \quad \forall f \in \mathcal { A } \.$$

Hence f n are dense in A . In fact, we learn something even more than this: { f n } and { g (˜) α ∆ n } are Schauder bases for A and HF respectively. This means that all elements have unique decompositions in terms of the basis elements.

Proof of Property 1') for GFF functionals. To show that the GFF functionals (A.7) indeed satisfy property 1'), consider the 'master functional' of [15],

$$f _ { w } ( z ) = \frac { 1 } { z - w } \.$$

This is a 1-parameter family of functionals that completely characterises locality:

$$F \text{ local } \iff \ \theta _ { f _ { w } } [ F ^ { ( \tilde { \alpha } ) } ] = 0 \.$$

That is a simple consequence of Cauchy integral formula

$$\theta _ { f _ { w } } [ H ] = \oint _ { \gamma } \frac { \mathrm d z } { 2 \pi i } \, \frac { 1 } { z - w } \, H ( z ) = H ( w ) \, \quad H \in \mathcal { H F }$$

where γ can be taken as any contour wrapping [1 , ∞ ], and with w on its outside.

Now we recall [15] that the master functional admits an expansion into the GFF functionals (A.7):

$$\frac { 1 } { z - w } = w ^ { - \tilde { \alpha } - 1 } \sum _ { n = 1 } ^ { \infty } \hat { f } _ { n } ( z ) \, g _ { \hat { \Delta } _ { n } } ( w ) \.$$

One can check that this is true, at least as a formal Taylor series for to each order O w ( N ). Moreover, as we shall now establish, it actually converges as an analytic function of z and as a hyperfunction of w : i.e. for z in some neighbourhood of [1 , ∞ ] and for w away from [1 , ∞ ]. This follows by the estimate given in Appendix A:

$$\hat { f } _ { n } ( z ) \, g _ { \hat { \Delta } _ { n } } ( w ) \stackrel { n \to \infty } { \sim } \left ( \frac { r ( w ) } { r ( z ) } \right ) ^ { \tilde { \alpha } + n } \times ( \dots )$$

where we only show the dependence on n and where r z ( ) maps C \ [1 , ∞ ] to the disk | r | &lt; 1,

$$r ( z ) = \frac { 1 - \sqrt { 1 - z } } { 1 + \sqrt { 1 - z } } \.$$

It follows that the sum converges absolutely for | r w ( ) | &lt; r z | ( ) | ≤ 1, as required.

We then observe that

$$w ^ { - \tilde { \alpha } - 1 } \, g _ { \hat { \Delta } _ { n } } ( w ) = g _ { \hat { \Delta } _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \.$$

This follows, for example by acting with equation (2.32) on g (˜) α ˆ ∆ n ( z ) and using the duality ( f m , g (˜) α n ) = δ mn . Hence we may re-write (2.32) as

$$\frac { 1 } { z - w } = \sum _ { n = 1 } ^ { \infty } \hat { f } _ { n } ( z ) \, g _ { \hat { \Delta } _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \.$$

Now substituting this decomposition into (2.31), one finds

$$H ( w ) = \sum _ { n = 1 } ^ { \infty } ( \hat { f } _ { n }, H ) \, g _ { \hat { \Delta } _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \quad \text{ for all } \ H \in \mathcal { H } \mathcal { F }$$

where commuting the series with the integral was possible because the sum converges as an analytic function of w away from [1 , ∞ ]. Note that our results imply that (2.37) converges in HF , proving Property 1').

## 3 A general construction: interacting functionals

In this section we will formally establish the existence of a large set of solutions of the conditions 1') and 2), leading to a large set of Schauder bases. A practical method to construct the corresponding kernels will later be given in Section 5, and an explicit expression for the functionals' actions on blocks will be given in Section 4.

Figure 3: The 'wedge' domain of analyticity W Λ ,δ (in blue) assumed for the functions α n ( ) and β n ( ) that interpolate the anomalous dimensions γ n according to equation (3.2). The dashed line is the real axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-5d7ccfa5-73ff-42c0-a47d-be74049fe5ad.jpg)

Concretely, we will produce functionals dual to a large class of possible { ∆ n } ,

$$\Delta _ { n } = \hat { \Delta } _ { n } + \gamma _ { n } \, \quad ( \hat { \Delta } _ { n } = 2 \tilde { \alpha } + 2 n )$$

which may be seen as 'small' deformations of GFF spectra. In detail, we assume that γ n takes the form

$$\gamma _ { n } = \alpha ( n ) + ( - 1 ) ^ { n } \beta ( n ) \,, \quad \alpha ( n ), \beta ( n ) \underset { n \to \infty } { = } O ( n ^ { - \epsilon } ) \,, \quad \epsilon > 0$$

with α n ( ) and β n ( ) analytic in n in some complex 'wedge' (see Figure 3), 14

$$W _ { \Lambda, \delta } \coloneqq \left \{ z \in \mathbb { C } \ | \ | \arg \left ( n - \Lambda \right ) | < \delta \ \right \} \,, \quad \ \Lambda > 0 \,, \delta > 0 \.$$

Note that it is sufficient to have Λ quite large (and δ quite small). We will call such spectra ' sparse ' and, for any sparse spectrum, we will argue the existence of a dual basis of functionals satisfying Properties 1') and 2) above.

14 For the case d = 1, we need to assume that the oscillating term is slightly more suppressed, β n ( ) ∼ n -1 2 -ϵ where ϵ ≥ 0.

## 3.1 A useful L 2 space

## 3.1.1 Embeddings

In order to make functionals dual to interacting spectra (3.1), one may think of 'rotating' from the GFF functionals (A.7) by a 'small' amount. Since A and HF are topological vector spaces without norms, it is not clear how to meaningfully measure such a rotation. To do this, we will use the trick of introducing an intermediate L 2 space:

$$\mathcal { A } \ \to \ \ L _ { \mu } ^ { 2 } \ \to \ \mathcal { H F } \,,$$

for an appropriate choice of inner product defined by:

$$\langle f | g \rangle _ { \mu } \coloneqq \int _ { 1 } ^ { \infty } \mathrm d z \, \mu ( z ) \, f ( z ) \, g ( z ) \, \quad \ \ f, g \in L _ { \mu } ^ { 2 } \.$$

Such a hierarchy of spaces (a Hilbert space in between a space A and its dual) is called a Gelfand triple or rigged Hilbert space . As we will see below, in practice, the L 2 space will buy us an elegant and user-friendly shortcut for constructing Schauder bases of functionals in A satisfying 1') and 2). First, we will construct bases for L 2 itself - a simpler and more familiar task - and the desired Schauder bases will follow under general conditions.

We need to choose the L 2 measure µ in (3.5) and describe the embeddings in (3.4). A natural choice for µ follows from the observation that the blocks g ∆ are eigenfunctions of a Casimir operator (see, e.g., [31]):

$$\left ( C _ { z } - \lambda _ { \Delta } \right ) g _ { \Delta } = 0 \, \quad \ \lambda _ { \Delta } = \Delta ( \Delta - d )$$

with

$$C _ { z } \coloneqq 4 ( 1 - z ) ^ { 1 - d / 2 } z ^ { 1 + d / 2 } \partial _ { z } [ ( 1 - z ) ^ { d / 2 } z ^ { 1 - d / 2 } \partial _ { z } ] \.$$

The Casimir operator is self-adjoint (up to boundary terms) on L 2 if we choose the particular measure

$$\mu ( z ) = \frac { 1 } { z ^ { 2 } } \left ( \frac { z } { z - 1 } \right ) ^ { 1 - d / 2 } \,,$$

so let us choose this one. With respect to this measure, the discontinuities of blocks are square-integrable: I z g ∆ ∈ L 2 µ . Henceforth, we will omit the subscript µ , but let us emphasise that there is only one L 2 norm used in this paper, defined with the measure (3.8).

Next, we need to describe the arrows appearing in the diagram (3.4), i.e. the embeddings of A into L 2 and of L 2 in HF . We would like the L 2 functions I z g ∆ to correspond to the subtracted blocks g (˜) α ∆ in the hyperfunction space, so let us define the following ˜-dependent α

embeddings: 15

↦

$$i ^ { ( \tilde { \alpha } ) } \colon & \mathcal { A } \to L ^ { 2 } \, & j ^ { ( \tilde { \alpha } ) } \colon & L ^ { 2 } \to \mathcal { H } \mathcal { F } \, \\ & f ( z ) \mapsto \pi ^ { - 1 } \mu ( z ) ^ { - 1 } z ^ { - \tilde { \alpha } - 1 } f ( z ) & \mathcal { I } _ { z } [ j ^ { ( \tilde { \alpha } ) } ( G ) ] ( z ) = z ^ { - \tilde { \alpha } - 1 } G ( z ) \. & \quad ( 3. 1 0 )$$

When unambiguous, we will use the streamlined notation

$$f ^ { ( \tilde { \alpha } ) } = i ^ { ( \tilde { \alpha } ) } ( f ) \, \quad f \in \mathcal { A } \, \quad \ \ \ G ^ { ( \tilde { \alpha } ) } = j ^ { ( \tilde { \alpha } ) } ( G ) \, \quad G \in L ^ { 2 } \. \quad \ \ ( 3. 1 1 )$$

The hyperfunction defining j (˜) α ( G ) may be written down explicitly as

$$j ^ { ( \tilde { \alpha } ) } ( G ) ( z ) = \int _ { 1 } ^ { \infty } \frac { \mathrm d w } { \pi } \, \frac { G ( w ) } { w ^ { \tilde { \alpha } + 1 } ( w - z ) } = \left \langle \frac { 1 } { \pi \mu ( w ) w ^ { \tilde { \alpha } + 1 } ( w - z ) } \Big | G ( w ) \right \rangle \,.$$

Importantly, we have the desired property j (˜) α ( I z g ∆ ) = g (˜) α ∆ .

The mappings i (˜) α and j (˜) α have several important properties:

- · The maps are compatible in the sense

$$( f, j ^ { ( \tilde { \alpha } ) } ( G ) ) = \langle i ^ { ( \tilde { \alpha } ) } ( f ) \, | \, G \rangle \, \quad \ \ \ f \in \mathcal { A } \, \ \ G \in L ^ { 2 } \. \quad \ \ ( 3. 1 3 )$$

- · Both i (˜) α and j (˜) α are continuous , i.e. they preserve convergence of sequences of functions. For i (˜) α this is trivial while for j (˜) α it follows from (3.12) and the CauchySchwarz inequality as long as ˜ α &gt; -1 2 .
- · Both i (˜) α and j (˜) α are injective (i.e. one-to-one).
- · Both i (˜) α and j (˜) α have dense image . 16

Altogether, these properties allow us to re-express the dual pairing of arbitrary elements of A and HF in terms of L 2 inner products, over which we will have better technical control.

## 3.1.2 Riesz bases

Now that we have set up the hierarchy of spaces, we would like to profit from the fact that L 2 is a Hilbert space by introducing a basis. Given the requirement 1') of our task, it seems natural to look for bases { p n } made out of discontinuities of blocks (i.e. the L 2 members

15 To make these maps well-defined (with good enough behaviour at z = 1), we require that ˜ α &gt; d 4 -1. Note that one can formally construct functionals for other values of ˜ but the statements of completeness α in the spaces A and HF as defined here apply only in this case.

16 For i (˜) α , this is related to the statement that the subset of L 2 functions suppressed by an extra power at ∞ is still dense. For j (˜) α , it follows from the statement (2.37) that any hyperfunction can be represented as a sum of hyperfunctions whose discontinuity is in L 2 .

which map to g (˜) α ∆ under j (˜) α ). A clue that this is the right track is that these discontinuities are given by:

$$\mathcal { I } _ { z \geq 1 } \, g _ { \Delta } ( z ) = \frac { \Gamma \left ( \frac { \Delta } { 2 } \right ) } { \Gamma \left ( \frac { \Delta } { 2 } + 1 - \frac { d } { 2 } \right ) } \, \left ( \frac { z - 1 } { z } \right ) ^ { 1 - \frac { d } { 2 } } J _ { \frac { \Delta - 2 } { 2 } } ^ { ( 1 - \frac { d } { 2 }, 0 ) } ( \frac { 2 - z } { z } ) \,.$$

For the particular case of ∆ n = 2 n even integers (i.e. ∆ n = 2˜ +2 α n with ˜ = 0), the Jacobi α functions reduce to Jacobi polynomials and the discontinuities form an orthonormal basis up to an overall normalisation factor. More generally, we will need to relax the requirement of orthonormality and ask instead that these discontinuities form instead what is known as a Riesz basis (see, e.g., the textbook [25]).

A Riesz basis is defined as the image { p n } of an orthonormal basis { e n } under a bounded, invertible map T :

$$p _ { n } = T ( e _ { n } ) \, \quad T \colon L ^ { 2 } \to L ^ { 2 } \text{ bounded, invertible} \.$$

It follows that any Riesz basis automatically comes with a dual Riesz basis , consisting of elements h n := T -1 ( e n ). These dual bases form a biorthogonal system:

$$\langle h _ { n } | p _ { m } \rangle = \delta _ { m n } \.$$

Crucially, any element in the Hilbert space has a unique, unconditionally convergent decomposition with respect to a Riesz basis:

$$f = \sum _ { n } \langle f | h _ { n } \rangle \, p _ { n } = \sum _ { n } \langle f | p _ { n } \rangle \, h _ { n } \, \quad \forall f \in L ^ { 2 } \.$$

Riesz bases are exactly the tool we need to achieve our task. On the one hand, by acting with j (˜) α , we may hope to lift completeness of the p n ∝ I z g ∆ n in L 2 into completeness in the space of hyperfunctions, i.e. property 1). On the other hand, the dual Riesz basis { h n } are prime candidates for the functional kernels, satisfying the duality property 2) by construction.

## 3.1.3 How to use L 2

Now that we have set up this L 2 space, let us explain how Riesz bases of L 2 can be used to construct complete bases of functionals solving our task. Suppose then that:

- (a) { p n := I z g ∆ n } forms a Riesz basis, and denote its dual Riesz basis { h n } . 17

17 In fact what we call p n , and their duals h n , are non-normalised Riesz bases - i.e. they become Riesz bases when normalised instead as p n → n 3 2 -d 2 p n and h n → n d 2 -3 2 h n (making their norms bounded) - see eq. (B.3) of Appendix B. In a slight abuse of language, we will still call them Riesz basis even when not appropriately normalised. This makes no difference to the properties that we will need.

$$( \mathbf b ) \ h _ { n } = i ^ { ( \tilde { \alpha } ) } ( f _ { n } ) \text{ for some } f _ { n } \in \mathcal { A } \.$$

A priori , from (a), the h n are only guaranteed to be L 2 functions so that (b) is a non-trivial requirement. It means that h n is real-analytic on (1 , ∞ ), bounded as | z | -˜ α at ∞ , and that h n ( z ) ∼ z → 1 ( z -1) 1 -d/ 2 [(analytic) + O z ( -1) ˜ α ]. In what follows, one of the major subtleties will be proving this assumption (b) to be satisfied for our class of 'sparse' spectra ∆ n . It will be satisfied exactly as stated here for ˜ α &lt; 1 2 ; for higher values of ˜, α we will instead need to consider particular finite linear combinations of the h n , but the conclusions will be unmodified (see Section 3.3 below).

Let us show that the resulting functionals f n ∈ A then solve the task set out above. Firstly, the duality property (3.16) of Riesz bases give

$$( f _ { m }, \, g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ) = \langle h _ { m } | p _ { n } \rangle = \delta _ { m n } \.$$

so that the functional kernels f n satisfy property 2). For property 1'), we may think of using the dual pair of Riesz bases to write down a partition of unity for the Hilbert space L 2 ,

$$\mathbf 1 = | p _ { n } \rangle \langle h _ { n } | \quad \Leftrightarrow \quad \delta ( z - w ) = \mu ( z ) \sum _ { n = 1 } ^ { \infty } p _ { n } ( w ) \, h _ { n } ( z ) \.$$

Of course such identities, familiar in quantum mechanics, only make rigorous sense when understood in the language of rigged Hilbert spaces. Since h n = i (˜) α ( f n ) for analytic functionals f n ∈ A in this special class, we can expect such a partition of unity to be valid when acting on hyperfunctions of the z variable. Thus, we may more rigorously rewrite the above distributional identity as the discontinuity of the following hyperfunction of z (cf. equation (2.36) in the GFF case):

$$\frac { 1 } { z - w } = \sum _ { n = 1 } ^ { \infty } g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, f _ { n } ( z ) \.$$

This is indeed a 'partition of unity', since the LHS acts as the identity operator on hyperfunctions: ( 1 z -w , H ( z )) z = H w ( ). As we will explain below, this decomposition is always correct, at least converging in some sense (its image under i (˜) α z converges in L 2 ). But for it to be useful, we must make an additional assumption that:

- (c) The sum (3.20) converges as an analytic function of z and a hyperfunction of w .

Then, copying the logic from the GFF case above, Property 1') follows by integrating (3.20) against H z ( ):

$$H ( w ) = \left ( \frac { 1 } { z - w }, H ( z ) \right ) _ { z } = \sum _ { n = 1 } ^ { \infty } ( f _ { n }, H ) \ g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \.$$

Let us come back to derive a version of equation (3.20) as an identity in L 2 . We would like to decompose the analytic function 1 z -w of z as a sum of the functions f n ( z ) ∈ A -but we do not know how to do that. Instead let us consider the L 2 element i (˜) α z ( 1 z -w ) , which we can decompose into the Riesz basis h n ( z ) = i (˜) α ( f n )( z ):

$$i ^ { ( \tilde { \alpha } ) } \left ( \frac { 1 } { z - w } \right ) = \sum _ { n } g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, i ^ { ( \tilde { \alpha } ) } ( f _ { n } ) ( z )$$

Indeed these are the correct w -dependent expansion coefficients: ⟨ i (˜) α z ( 1 z -w ) | p n ⟩ = g (˜) α ∆ n ( w ). The assumption (c) implies that one can act on this equation with ( i (˜) α z ) -1 , obtaining (3.20) converging as an analytic function of z .

Gathering all of the assumptions we used, we have established (for the range ˜ α &lt; 1 2 ) the following sufficient conditions in L 2 language to produce a basis of functionals f n ∈ A satisfying Properties 1') and 2):

## Criteria in L 2 language: (˜ α &lt; 1 2 )

- (a) A set of discontinuities of blocks { p n = I z g ∆ n } forms a Riesz basis for L 2 .
- (b) The dual Riesz basis { h n } are given by h n = i (˜) α ( f n ) for some f n ∈ A .
- (c) The decomposition (3.20) converges as an analytic function of z and a hyperfunction of w .

## 3.2 Completeness of functionals using L 2

We have derived a set of sufficient conditions (a)-(c) for the existence of a complete basis of analytic functionals dual to a given spectrum ∆ n . In this section, we will first show that these conditions are met by the GFF spectrum ∆ n = 2˜+2 , giving an alternate derivation α n of the completeness of the GFF functionals. Then we will turn to interacting sparse spectra ∆ n = 2˜ + 2 α n + γ n and demonstrate the existence of a dual basis of analytic functionals for any γ n satisfying the assumptions (3.2). The analysis in this section focuses on the case ˜ α &lt; 1 2 , but the conclusions are also valid in the general case, which is treated below in Section 3.3.

## 3.2.1 GFF functionals

We begin by showing that the conditions (a)-(c) are satisfied for the spectrum ∆ ˆ n = 2˜+2 . α n This will be instructive for proving the same thing for interacting functionals. As above, we use hats to denote GFF quantities (as opposed to interacting ones).

## (a) For any ˜ α ∈ -( 1 2 , 1 2 ) , the GFF double trace blocks ˆ p n := I z g 2˜+2 α n form a Riesz basis for L 2 .

As mentioned before, this is natural because we know that, for the value ˜ = 0, the α ˆ p n are related to an orthonormal basis of Jacobi polynomials (see (3.14)). Therefore, we expect that, for close enough values of ˜, α the basis property is preserved (although not orthonormality).

Proving this mathematical statement is rather technical, and is relegated to Appendix B.1. Here will only give a sketch of the argument. The task is to construct a bounded, invertible operator

↦

$$T \colon L ^ { 2 } & \to L ^ { 2 } \\ e _ { n } & \mapsto \hat { p } _ { n }$$

where e n is some orthonormal basis (there is no need to specify which, since any is equivalent). In this case, we already explicitly know of a dual family of L 2 functions, ˆ h n := i (˜) α ( ˆ ), f n obtained by embedding the GFF functional kernels (A.7) into L 2 , which satisfy

$$\langle \hat { h } _ { n } | \hat { p } _ { m } \rangle = \delta _ { m n } \.$$

This allows us to construct an inverse map by defining an operator

↦

$$S \colon L ^ { 2 } & \to L ^ { 2 } \\ e _ { n } & \mapsto \hat { h } _ { n } \.$$

Assuming that T and S extend to well-defined, bounded operators, (3.24) gives

$$S ^ { \dagger } \, T = 1 \.$$

This is almost enough to conclude that S † = T -1 , but formally for a Hilbert space it only implies that T is left-invertible. However, the right-invertibility of T (or the left-invertibility of S ∗ ) may be argued using the fact that the ˆ f n are polynomials, so that ˆ f (˜) α n can be written as finite combinations of some known orthonormal basis (see Appendix B.1).

The only remaining step is to show that (3.23) and (3.25) give well-defined, bounded linear operators T and S . This is shown in Appendix B.1 by computing and estimating the matrix elements ⟨ ˆ p n | ˆ p m ⟩ and ⟨ ˆ h n | ˆ h m ⟩ for large n, m . We then apply Schur's test - a classic mathematical tool giving a sufficient condition for such operators to be bounded.

## (b) Analyticity and bounds on functionals.

For the GFF case, we already know the dual Riesz basis ˆ h n : by construction, they must match the GFF functionals (A.7) from our previous work: ˆ h n = i (˜) α ( ˆ ). f n Hence property

## (b) is immediate.

## (c) Decomposition of the Cauchy kernel.

The decomposition (3.20) was already derived in the GFF case in eq. (2.32). As we explained above, it can also be derived as an L 2 decomposition with respect to the Riesz basis ˆ h n . The estimate (2.33) then applies, implying that this decomposition converges as an analytic function of z and a hyperfunction of w .

## 3.2.2 Interacting functionals using L 2

Now we are ready to do the same thing for interacting spectra, ∆ n = 2˜ + 2 α n + γ n , i.e. prove that the properties (a)-(c) are satisfied.

## (a) For any ˜ α ∈ ( -1 2 , 1 2 ) , the interacting set of blocks p n := I z g 2˜+2 α n + γ n form a Riesz basis for L 2 µ .

Our strategy to establish this is to study the mapping between the GFF blocks and the interacting ones:

↦

$$M \colon \ \hat { p } _ { n } \mapsto p _ { n } \.$$

In particular, if we can show that M defines a bounded, invertible operator, then so is the composition MT which now relates p n to an orthonormal basis (where T is the map in (3.23)).

Our approach is first to show that (1 -M ) is a compact operator so that, in particular, M is bounded. Essentially this follows from the fact that the γ n decay for large n so that the inner products of p n are sufficiently close to those of ˆ . p n Then, using the fact that the p n are distinct eigenfunctions of a Casimir operator, this implies that M is injective, and since (1 -M ) is compact this implies it is also invertible as desired. The detailed argument may be found in Appendix B.2.

## (b) Analyticity and bounds on functionals.

Unlike the GFF case, we do not have a stand-alone definition of the interacting functionals kernels h n : they are simply defined as the Riesz basis to dual p n . However, they may be expanded as an L 2 -convergent sum of the kernels ˆ h m (the latter being a Riesz basis):

$$h _ { n } = \sum _ { m } \langle h _ { n } | \hat { p } _ { m } \rangle \, \hat { h } _ { m } \.$$

We would like to prove that h n = i (˜) α ( f n ) for some f n ∈ A , i.e. that h n are real-analytic

and suitably bounded at ∞ and 1. This is carefully argued in Appendix C, making crucial use of the assumptions (3.2) of polynomial boundedness and analyticity of the anomalous dimensions. The argument consists of three steps:

- 1. The expansion coefficients above can be written in terms of the putative interacting functional actions:

$$b _ { m } \coloneqq \langle h _ { n } | \hat { p } _ { m } \rangle = ( f _ { n }, g _ { \hat { \Delta } _ { m } } ^ { ( \tilde { \alpha } ) } ) = \theta _ { n } ( \hat { \Delta } _ { m } )$$

We show that these coefficients b m admit a suitable analytic continuation into the complex m plane.

- 2. Next, using this analytic continuation, we prove that the interacting functional kernels given by sums ∑ m b m m ˆ h , are real-analytic functions. We do this via a SommerfeldWatson type trick.
- 3. Finally we argue that the decaying anomalous dimensions lead to the right boundedness properties at z = 1 and ∞ . Roughly this is because, at large ∆, interacting and GFF functionals behave similarly, and that limit is controlled by the z = 1 and z = ∞ regions (however, the interacting functionals are non-analytic at the endpoints 1 and ∞ , and may have a branch cut along ( -∞ , 1)).

Let us elaborate on the first step. We will be able to construct this analytic continuation explicitly thanks to the explicit formula derived in Section 4 which expresses the interacting functional actions θ n (∆) in terms of the sparse spectrum { ∆ n } . Essentially the reason why this is possible is analyticity of θ n ( ˆ ∆ ) only has to hold for sufficiently large m m . In this regime we have θ n ( ˆ ∆ ) m ∼ -γ m ˆ θ n ′ ( ˆ ∆ ), and the result now follows from our explicit m knowledge of the ˆ θ n and the assumptions made on the analytic properties of the γ m , cf. (3.2) and (3.3).

## (c) Decomposition of the Cauchy kernel.

The decomposition (3.20) can again be argued to converge as an analytic function of z and a hyperfunction of w , using the same estimate (2.33) as in the GFF case. Since the anomalous dimensions are assumed to decay and g ∆ ( w ) has exponential dependence r w ( ) ∆ we have

$$g _ { \Delta _ { n } } ( w ) \underset { n \to \infty } { \sim } g _ { \hat { \Delta } _ { n } } ( w ) \left ( 1 + O ( n ^ { - \epsilon } ) \right ) \,.$$

Similarly, one can show

$$f _ { n } ( z ) = \hat { f } _ { n } ( z ) \left ( 1 + O ( n ^ { - \epsilon } ) \right ) \,.$$

Essentially this is because for large n we have p n ∼ ˆ p n and hence h n ∼ ˆ h n . (More rigorously, this follows from eq. (5.12), or by substituting (C.38) into (C.40)). Hence the estimate

(2.33) carries over to the interacting case, thereby guaranteeing the absolute convergence of (3.20) for | r w ( ) | &lt; r z | ( ) | &lt; 1, as required.

## 3.3 On general values of ˜ α and 'subtractions'

The criteria (a)-(c) above focus on the particular case ˜ α ∈ ( -1 2 , 1 2 ). However, the special role of ˜ = α 1 2 is simply an artefact of our use of L 2 , which is, after all, only a convenient mathematical trick. In terms of the physically meaningful spaces A and HF , all of the conclusions above still hold for general values of ˜. α Moreover, they can still be derived using L 2 , but through a slightly more intricate relation interpreted as a 'subtraction' procedure.

In the following, we shall generalise the argument above, underlining the necessary assumptions as we use them. Casual readers may wish to skip this section. On the other hand, interested readers should note that there is a remaining boundary case ˜ α ∈ 1 2 + Z , which is treated in Appendix D.3.

For ˜ α &gt; 1 2 , clearly the blocks I z g ∆ n with ∆ n ∼ 2˜ + 2 α n cannot be a Riesz basis for L 2 , since we know that (a) they are a proper subset of a Riesz basis (by adding a few additional blocks):

$$\{ p _ { n } = \mathcal { I } _ { z } \, g _ { \Delta _ { n } ^ { \beta } } \} \, \quad \ \Delta _ { n } ^ { \beta } \coloneqq \begin{cases} \Delta _ { n } ^ { \text{add} } \, \quad & 1 \leq n \leq M \, \\ \Delta _ { n - M } ^ { \tilde { \alpha } } \, \quad & n > M \, \end{cases}$$

where the additional scaling dimensions ∆ add 1 , . . . , ∆ add n can be chosen arbitrarily. This new spectrum then has the large n behaviour ∆ β n ∼ 2 β +2 , where n

$$\beta = \tilde { \alpha } - M \in ( - \frac { 1 } { 2 }, \frac { 1 } { 2 } ) \, \quad M \in \mathbb { N } \.$$

Applying the argument of Section 3.1.3 with this Riesz basis, we reach the following equation before encountering any issues:

$$i _ { z } ^ { ( \tilde { \alpha } ) } \Big ( \frac { 1 } { z - w } \Big ) = \sum _ { n } g _ { \Delta _ { n } ^ { \beta } } ^ { ( \tilde { \alpha } ) } ( w ) \, h _ { n } ( z ) \.$$

In the ˜ α &lt; 1 2 case, we had assumed h n = i (˜) α ( f n ), but for ˜ α &gt; 1 2 that assumption will not be true. Instead, (b) certain 'subtracted' combinations of h n will be in the image of i (˜) α :

$$h _ { n + M } - \sum _ { k = 1 } ^ { M } c _ { n + M, k } ^ { \tilde { \alpha } } \, h _ { k } = i ^ { ( \tilde { \alpha } ) } ( f _ { n } ) \, \quad f _ { n } \in \mathcal { A }$$

at the price of 'losing' M functionals h , . . . , h 1 M . Due to the duality of the spaces A and H F , the loss of these functionals is reflected in the fact that the hyperfunctions g (˜) α ∆ β n ( w ) on the RHS of (3.34) are not independent when viewed as elements of HF . Intuitively, this is

because they are 'oversubtracted' by ˜, which is α M units greater than β . As a result, there are M relations between them (see Appendix D.1) involving precisely the same coefficients,

$$g _ { \Delta _ { k } ^ { \beta } } ^ { ( \tilde { \alpha } ) } ( w ) = - \sum _ { n > M } c _ { n k } ^ { \tilde { \alpha } } \, g _ { \Delta _ { n } ^ { \beta } } ^ { ( \tilde { \alpha } ) } ( w ) \, \quad k = 1, \dots, M \.$$

Using these relations to eliminate the lowest M blocks, the remaining expression only contains the original blocks (not the added ones) and the subtracted combinations of functionals (3.35):

$$i _ { z } ^ { ( \tilde { \alpha } ) } \Big { ( } \frac { 1 } { z - w } \Big { ) } = \sum _ { n } g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, i ^ { ( \tilde { \alpha } ) } ( f _ { n } ) ( z ) \.$$

Now we would like to apply ( i (˜) α ) -1 and pull it through the sum we obtain the same result (3.20) as above:

$$\frac { 1 } { z - w } = \sum _ { n = 1 } ^ { \infty } g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, f _ { n } ( z ) \.$$

We now need to assume that (c) this decomposition converges as an analytic function of z and a hyperfunction of w . As in the unsubtraced case, this implies that the functionals f n solve our task and form a complete Schauder basis for A .

Let us comment that the arbitrary choice of the M additional scaling dimensions ∆ add n drops out in the eventual functionals f n ∈ A . This must be the case, since f n are dual to the Schauder basis g (˜) α ∆ n , so the difference between functionals generated by making different choices is orthogonal to all basis elements, and hence vanishes.

Gathering up the assumptions, a set of sufficient conditions for a complete basis of analytic functionals for general values of ˜ is: α

## Criteria in L 2 language: (general ˜ α )

- (a) A set of discontinuities of blocks { p n = I z g ∆ n } can be completed to a Riesz basis for L 2 by adding a finite number M of blocks.
- (b) Finite combinations of M of the dual basis Riesz basis { h n } take the form i (˜) α ( f n ) for some f n ∈ A .
- (c) The decomposition (3.38) converges as an analytic function of z and a hyperfunction of w .

In Appendix D.1, we show that any 'sparse' spectrum ∆ n = 2˜+2 α n + γ n (with γ n satisfying the assumptions (3.2)) satisfies these conditions - generalising the arguments from the ˜ α &lt; 1 2 case above in Section 3.2.

## 3.4 Completeness and extremality

## 3.4.1 Completeness

The results of the previous sections show that certain spectra ∆ n = 2˜ + 2 α n + γ n possess complete bases of functionals dual to them. In this section we will discuss in more detail the meaning of completeness and what exactly can be learned from the sum rules associated to these functionals.

We have established the following statement:

## Locality of form factors:

Suppose ∑ ∆ | c ∆ || g ∆ ( z ) | satisfies the assumed polynomial bound (2.5) for some α F . Then:

$$\Big | \quad F ( z ) = \sum _ { \Delta } c _ { \Delta } \, g _ { \Delta } ( z ) \ \text{ local} \quad \Longleftrightarrow \quad \sum _ { \Delta } c _ { \Delta } \, \theta _ { n } ( \Delta ) = 0 \ \text{ for } n \geq 1$$

where the basis of functionals on the RHS can have any value ˜ α &gt; α F -1.

It may seem puzzling that the choice of ˜ above is almost arbitrary. α For instance, when the choice of basis has γ n = 0, then increasing ˜ by one unit seems to remove a constraint, as α the resulting functionals are now dual to a proper subset of the basis elements. How can this be?

The answer lies in the assumptions of polynomial boundedness. While it is true that, as we increase ˜, we get a less restrictive set of sum rules, this just reflects the fact that the sum α rules progressively allow for local solutions with harder and harder polynomial behaviour. Since the correct behaviour is being imposed as an external assumption on the allowed c ∆ , then there is no contradiction.

Relaxing this a priori requirement on the c ∆ can be at least partially achieved by working with the largest possible set of functionals. To see this, note that the polynomial boundedness of the form factor F roughly correlates with the growth of the c ∆ for large ∆. Since θ n (∆) ∼ ∆ -2˜+ α d 2 -3 (cf. Appendix A), we see that lower values of ˜ automatically rule α out unwanted kinds of polynomial growth of F . To avoid as many extraneous solutions as possible from the start, one ought to choose ˜ not too far from α α F . In practice, a natural sweet spot is to choose ˜ = α α F . The reason is that, if we believe some bulk operator to have a local form factor with polynomial behaviour given by α F , then there are certainly local form factors of the same operator acted on by the bulk AdS Laplacian, whose growth will now be α F + 1. To pick up the first but not the second in our solutions to locality, we should choose ˜ α &gt; α F -1 and ˜ α ≤ α F . To speed up convergence of the sum rules, the maximum value ˜ = α α F is a good choice.

Hyperfunctions vs. tempered distributions. This perspective of choosing the largest possible set of functionals suggests the following question. Given that we know discontinuities of form factors actually live in the smaller space of tempered distributions 18 rather than hyperfunctions, we could have chosen this space from the start and worked with functionals in its dual space of Schwartz test functions, which is larger than A . This means that there exist additional sum rules, distinct from ours, that could have been written down. What do such sum rules buy us? Clearly they do not change anything regarding locality: the space A is dense in the Schwartz space, so our bases of functionals are still complete in that larger space. The new sum rules would simply be more wild combinations of our functionals, converging only in a generalised sense, but not in A .

Instead, such exotic sum rules should be again understood as replacing some a priori requirements on the c ∆ . In particular, they would constrain the c ∆ such that ∑ ∆ | c ∆ || g ∆ | is polynomially bounded - and otherwise the new sum rules would not converge. Note that this is not the case for our own sum rules: strictly speaking, our sum rules require only polynomial boundedness of F near the points z = 1 and z = ∞ . This is evident from the fact that the functional actions are defined by contour integrals touching the cut z ∈ (1 , ∞ ) only at these endpoints. A sufficient but not necessary condition for achieving this is indeed to request polynomial boundedness of ∑ ∆ | c ∆ || g ∆ | as we have in (2.5). But this can be relaxed, e.g., to

$$\left | \sum _ { \Delta > \Delta ^ { * } } c _ { \Delta } \, g _ { \Delta } ( z ) \right | \lesssim \begin{cases} \begin{array} { c c } | z | ^ { \alpha _ { F } } & \text{as $z\to\infty$} \\ | z - 1 | ^ { - \alpha _ { F } } & \text{as $z\to 1$} \end{array} \end{cases} \quad \text{for all $\Delta^{*}$} \,.$$

Thus our sum rules are actually compatible with a wilder set of c ∆ , e.g. those which a priori could have led to form factors with essential singularities on the cut z ∈ (1 , ∞ ). The extra functionals gained by working in the space of Schwartz functions would rule out such esoteric possibilities. By working with hyperfunctions with their duals, we instead choose to rule these out by directly constraining the allowed c ∆ .

## 3.4.2 Extremal solutions

Moving on, let us now note that our locality sum rules have particularly simple solutions which we call ' extremal '. For each complete basis of functionals, such solutions can be labelled by a choice of fixed BOE data which acts as inputs. These are then completed by the sparse set { g ∆ n } of blocks dual to our functionals. In the simplest case of a single input block g ∆ , we write the ansatz for a tentative local form factor.

$$\mathcal { L } _ { \Delta } ^ { \{ \Delta _ { n } \} } = g _ { \Delta } + \sum _ { n = 1 } ^ { \infty } c _ { n } \, g _ { \Delta _ { n } } \.$$

18 This follows from the analysis of [23] for the crossing equation, since our assumed polynomial bound implies the 'slow-growth' condition there.

Applying the functionals θ n to the above construction, the sum rules are 'diagonalised' by duality, leading to an immediate solution for the coefficients:

$$c _ { n } = - \theta _ { n } ( \Delta ) \.$$

The resulting extremal solutions are thus interacting generalisations of the 'local block' solutions of our previous work [15], to which they reduce for the special case ∆ n = ˆ ∆ . n In Section 4 below we will give an explicit expression for θ n (∆) in terms of { ∆ n } - see equation (4.20).

We note that the above local blocks should correspond to an alternate formulation of the locality constraint of the form

$$F \text{ is local } \quad \Leftrightarrow \quad \sum c _ { \Delta } \left [ g _ { \Delta } - \mathcal { L } _ { \Delta } ^ { \{ \Delta _ { n } \} } \right ] = 0$$

$$\Leftrightarrow \quad \sum _ { \Delta } c _ { \Delta } \left [ g _ { \Delta } - \mathcal { L } _ { \Delta } ^ { \psi \psi ^ { n j } } \right ] = 0$$

That is, for local form factors the usual conformal block expansion may be swapped by a manifestly local one. Again, this was established in [15] for the special GFF case, but we expect it to be generally valid for any set { ∆ n } satisfying our assumptions. To prove it we would have to show that the 'master functional' [32] for these bases is a valid functional.

## 4 The space of functional actions

Until now, we have been concerned with the space A of functional kernels, and have constructed bases for it. This section proposes to consider instead the space of functional actions ,

$$\theta ( \Delta ) = \langle h _ { \theta } \, | \, \mathcal { I } _ { z } g _ { \Delta } \rangle = \int _ { 1 } ^ { \infty } \mathrm d z \, \mu ( z ) \, h _ { \theta } ( z ) \, \mathcal { I } _ { z } g _ { \Delta } ( z ) \.$$

These are the functions of the complex variable ∆ defined by integrating functional kernels h θ against the blocks g ∆ . If one can understand the space of such functional actions θ (∆), then one understands all sum rules, ∑ ∆ c ∆ θ (∆) = 0, characterising the locality constraints.

In this section it will be convenient to define:

$$P _ { \lambda _ { \Delta } } ( z ) \coloneqq \mathcal { I } _ { z } \, g _ { \Delta } ( z ) = \frac { \Gamma ( \frac { \Delta } { 2 } ) } { \Gamma ( \frac { \Delta + 2 - d } { 2 } ) } \, \left ( \frac { z - 1 } { z } \right ) ^ { 1 - \frac { d } { 2 } } \, J _ { \frac { \Delta - 2 } { 2 } } ^ { ( 1 - \frac { d } { 2 }, 0 ) } ( \frac { 2 - z } { z } ) \.$$

This is indeed a function of the Casimir eigenvalue λ ∆ = ∆(∆ -d ), since it is invariant under ∆ ↔ d -∆. In the remainder of this section we will abuse notation and equate θ (∆) ≡ θ λ ( ∆ ).

Next, we introduce an operator that formalises the mapping between kernels and

functional actions. We set:

$$\theta \coloneqq K [ h _ { \theta } ] \, \quad K [ h ] ( \lambda ) \coloneqq \int _ { 1 } ^ { \infty } \mathrm d z \, \mu ( z ) \, h ( z ) \, P _ { \lambda } ( z ) \.$$

The image of this transform is typically called a Paley-Wiener space:

$$\mathcal { P W } ^ { ( \tilde { \alpha } ) } \coloneqq K [ i ^ { ( \tilde { \alpha } ) } ( \mathcal { A } ) ] \quad \subset \quad \mathcal { P W } \coloneqq K [ L ^ { 2 } ] \.$$

We recall that that i (˜) α ( A ) is the subspace of L 2 -functions that are real-analytic and sufficiently bounded at the endpoints - so as to have finite action on hyperfunctions and produce valid sum rules. Its image PW (˜) α is therefore the space of actual functional actions. But, as above in Section 3, it will be useful to consider the larger space PW .

The operator K : L 2 → PW is a bijection: i.e. to each functional action θ λ ( ) ∈ PW , there corresponds precisely one functional kernel h θ ∈ L 2 . Surjectivity follows by the definition of PW . For injectivity, we note that

$$K [ h ] ( \lambda ) \coloneqq \langle h \, | \, P _ { \lambda } \rangle = 0 \ \forall \lambda \ \implies \ \langle h \, | \, P _ { \lambda _ { n } } \rangle = 0 \ \forall n \geq 1 \ \implies \ \ h = 0 \, \quad ( 4. 5 )$$

for any set P λ n forming a Riesz basis of L 2 . (We have shown that there are many.) The Hilbert space structure on L 2 then defines a Hilbert space structure on PW with inner product

$$\langle \theta | \theta ^ { \prime } \rangle \coloneqq \langle h _ { \theta } | h _ { \theta ^ { \prime } } \rangle$$

so that K is a Hilbert space isomorphism. With this definition, since our analytic functionals i (˜) α ( A ) are dense in L 2 , then their image PW (˜) α is clearly a dense subspace of PW .

## 4.1 An intrinsic definition

We would like to give an intrinsic definition of the space PW , i.e. one that does not rely on the mapping K . To see how this may be achieved, let us examine more closely the properties of the functions θ λ ( ). First, we note that θ λ ( ) are entire functions (analytic functions with no singularities in C ): this follows since the integral in (4.3), which is controlled by L 2 , converges locally uniformly over λ ∈ C . Moreover, it follows from the asymptotic form (given in detail in (A.16))

$$P _ { \lambda _ { \Delta } } ( z ) \underset { | \Delta | \gg 1 } { \sim } \sin \left [ \frac { \Delta } { 2 i } \log ( r ( z ) ) \right ] \times ( \dots )$$

that θ λ ( ) is bounded in all directions of the complex plane as

$$| \theta ( \lambda _ { \Delta } ) | \leq C \frac { \exp \left ( \frac { \pi } { 2 } | \Delta | \right ) } { | \Delta \, \text{Im} ( \Delta ) | ^ { 1 / 2 } } \.$$

In particular, it is an entire function of λ of order 1 / 2 and type π/ 2, i.e.

$$| \theta ( \lambda ) | \leq D \exp \left ( \frac { \pi } { 2 } \sqrt { | \lambda | } \right ) \quad \text{ for all } \ \lambda \in \mathbb { C } \,.$$

For positive real λ , the functional actions are actually much more constrained. To see this, first let us set

$$\check { \lambda } _ { n } \coloneqq \hat { \lambda } _ { n } | _ { \tilde { \alpha } = 0 } = 2 n ( 2 n - d ) \.$$

It is easily checked from the relation of P ˇ λ n to Jacobi polynomials that

$$\langle P _ { \tilde { \lambda } _ { n } } | P _ { \tilde { \lambda } _ { m } } \rangle = \delta _ { n, m } | | P _ { \tilde { \lambda } _ { n } } | | ^ { 2 } \, \quad | | P _ { \tilde { \lambda } _ { n } } | | ^ { 2 } = \frac { 1 } { \sqrt { \frac { d ^ { 2 } } { 4 } + \check { \lambda } _ { n } } \left ( \frac { \check { \Delta } _ { n } } { 2 } \right ) _ { 1 - \frac { d } { 2 } } ^ { 2 } } \.$$

Finiteness of the L 2 norm of h θ now implies

$$\sum _ { n = 1 } ^ { \infty } \frac { | \theta ( \check { \lambda } _ { n } ) | ^ { 2 } } { | | P _ { \check { \lambda } _ { n } } | | ^ { 2 } } < \infty \, \quad$$

which roughly translates as

$$| \theta ( \lambda ) | = O ( \lambda ^ { \frac { d } { 4 } - 1 - \epsilon } ) \quad \text{ for } \lambda > 0 \quad \text{ (where } \epsilon > 0 \text{)} \.$$

We are ready to write down a tentative definition of the space PW : it is the space of entire functions satisfying conditions (4.8),(4.12) equipped with the inner product

$$\langle \theta | \theta ^ { \prime } \rangle \coloneqq \sum _ { n = 1 } ^ { \infty } \frac { \theta ( \check { \lambda } _ { n } ) \theta ^ { \prime } ( \check { \lambda } _ { n } ) } { | | P _ { \check { \lambda } _ { n } } | | ^ { 2 } }$$

To prove that this is indeed our original space we must show that to any element θ there corresponds an L 2 function h θ such that θ = K h [ θ ]. Starting from θ λ ( ) let us then set:

$$h _ { \theta } ( z ) = \sum _ { n = 1 } ^ { \infty } \theta ( \check { \lambda } _ { n } ) \frac { P _ { \check { \lambda } _ { n } } ( z ) } { | | P _ { \check { \lambda } _ { n } } | | ^ { 2 } } \.$$

This sum converges in L 2 thanks to condition (4.12). Furthermore, with this definition, we

have ⟨ θ θ | ′ ⟩ = ⟨ h θ | h θ ′ ⟩ as required. It remains to show that the prescription (4.15) satisfies K h [ θ ]( λ ) = θ λ ( ). To see this, we first compute:

$$K [ h _ { \theta } ] ( \lambda ) = \sum _ { n = 1 } ^ { \infty } \theta ( \check { \lambda } _ { n } ) K _ { \check { \lambda } _ { n } } ( \lambda )$$

with

$$K _ { \tilde { \lambda } _ { n } } ( \lambda ) \coloneqq \int _ { 1 } ^ { \infty } \mathrm d z \, \mu ( z ) \, \frac { P _ { \lambda _ { n } } ( z ) P _ { \lambda } ( z ) } { | | P _ { \tilde { \lambda } _ { n } } | | ^ { 2 } } = \frac { 4 } { \pi } \, ( - 1 ) ^ { n } \frac { \sin [ \frac { \pi } { 2 } \Delta ] } { \lambda - \check { \lambda } _ { n } } \, \frac { \left ( \frac { \tilde { \Delta } _ { n } } { 2 } \right ) _ { 1 - \frac { d } { 2 } } } { \left ( \frac { \Delta } { 2 } \right ) _ { 1 - \frac { d } { 2 } } } \.$$

To see that this indeed is the same as θ λ ( ), consider the identity

$$\theta ( \lambda ) = \oint _ { \lambda } \frac { d \lambda ^ { \prime } } { 2 \pi i } \, \frac { \theta ( \lambda ^ { \prime } ) } { \lambda ^ { \prime } - \lambda } \, \frac { \sin [ \frac { \pi } { 2 } \Delta ] } { \sin [ \frac { \pi } { 2 } \Delta ^ { \prime } ] } \, \frac { \left ( \frac { \Delta ^ { \prime } } { 2 } \right ) _ { 1 - \frac { d } { 2 } } } { \left ( \frac { \Delta } { 2 } \right ) _ { 1 - \frac { d } { 2 } } } \.$$

Note that the factors inside the integrand should be really thought of as functions of λ, λ ′ by setting e.g. ∆ = d 2 + √ ( d 2 ) 2 + λ . In particular, it can be checked that they respect invariance under ∆ ↔ -d ∆. Deforming the contour, we can drop the arcs thanks to assumptions (4.8) and (4.12), leading to a sum over poles, which precisely reproduces the expression (4.16). Thus we indeed conclude θ = K h [ θ ].

This completes the construction of the space PW . Since PW is isomorphic to L 2 , Riesz bases of the latter lift into Riesz bases of the former. In particular, the Riesz bases h n that we constructed for L 2 (dual to a sparse set of blocks p ∆ n ) are mapped to Riesz bases θ n := K h [ n ] for PW , and any functional action can be decomposed into such a basis:

$$\theta ( \lambda ) = \sum _ { n = 1 } ^ { \infty } \theta ( \lambda _ { n } ) \, \theta _ { n } ( \lambda ) \.$$

## 4.2 Exact interacting functional actions: a product formula

The interacting functional actions live in the PW (˜) α space. The classic factorisation theorem of Hadamard tells us that entire functions of order less than 1 are entirely determined by their zeros. This is in fact valid for any θ ∈ PW , and in particular the subspace PW (˜) α . Consider then an interacting functional basis element. Since we know the positions of its zeros, then the following remarkable result follows almost directly:

Consider a basis { θ m } dual to { ∆ n } , i.e. satisfying θ m ( λ n ) = δ mn . They are given by

̸

$$\theta _ { n } ( \lambda ) = \prod _ { \substack { m = 1 \\ m \neq n } } ^ { \infty } \left ( \frac { \lambda - \lambda _ { m } } { \lambda _ { n } - \lambda _ { m } } \right ) \,.$$

Before we continue, there is a small argument necessary to actually establish this formula, since we must prove that the set λ n exhausts all of the zeros of the functional actions. To show this, let us first check that the formula above holds for the GFF bases. Starting from expression (A.12), we can rewrite it as (with now ˆ λ n = 2˜ + 2 α n ):

$$\hat { \theta } _ { n } ( \lambda _ { \Delta } ) = \lim _ { \Delta ^ { \prime } \to \Delta _ { n } } \left ( \frac { \lambda _ { \Delta ^ { \prime } } - \hat { \lambda } _ { n } } { \lambda _ { \Delta } - \hat { \lambda } _ { n } } \right ) \, \frac { \Gamma ( \frac { 2 + 2 \tilde { \alpha } - \Delta ^ { \prime } } { 2 } ) \Gamma ( \frac { 2 + 2 \tilde { \alpha } + d - \Delta ^ { \prime } } { 2 } ) } { \Gamma ( \frac { 2 + 2 \tilde { \alpha } - \Delta } { 2 } ) \Gamma ( \frac { 2 + 2 \tilde { \alpha } + d - \Delta } { 2 } ) }$$

and (4.20) now follows from the identity Γ(1 + z ) = ∏ ∞ m =1 (1+ 1 n ) z 1+ z n . We can now show for the interacting bases that there are still no other zeros. Firstly, for any basis dual to a spectrum with γ n = ∆ n -ˆ ∆ decaying, the largen λ behaviour of θ n ( λ ) and ˆ θ n ( λ ) must be the same, as argued in eq. (5.5) below. Let us then write θ n ( λ ) as a polynomial factor R N ( λ ) capturing the possible extra zeros times the infinite product in (4.20). Since λ n ∼ ˆ λ n for large n , we can write

$$\frac { \theta _ { n } ( \lambda ) } { \hat { \theta } _ { n } ( \lambda ) } \sim \frac { R _ { n } ^ { N } ( \lambda ) } { \hat { \theta } _ { n } ( \lambda _ { n } ) } \prod _ { m = 1 } ^ { N - 1 } \left [ \left ( \frac { \lambda - \lambda _ { m } } { \lambda - \hat { \lambda } _ { m } } \right ) \left ( \frac { \lambda _ { n } - \hat { \lambda } _ { m } } { \lambda _ { n } - \lambda _ { m } } \right ) \right ]$$

Requiring the same large λ behaviour now fixes R N n ( λ ) ∼ 1.

Incidentally, this expression is also a useful way to approximate the infinite product in practice. One finds to subleading order in N :

$$R _ { n } ^ { N } ( \lambda ) \sim 1 + ( \lambda - \hat { \lambda } _ { n } ) \sum _ { m = N } ^ { \infty } \frac { \lambda _ { m } - \hat { \lambda } _ { m } } { ( \lambda - \hat { \lambda } _ { m } ) ( \hat { \lambda } _ { n } - \hat { \lambda } _ { m } ) } \.$$

This makes clear that for fixed λ , ˆ λ n and γ n = O n ( -ϵ ), the correction term to R N ∼ 1 is O N ( - -2 ϵ ). In applications, we know γ n as an analytic function of n and it is possible to do better. For instance, with γ n ∼ g/n 2 we get

$$\frac { R _ { n } ^ { N } ( \lambda ) - 1 } { g } \sim \frac { 2 \log \left ( 1 - \frac { \dot { \lambda } _ { n } } { \lambda _ { N } } \right ) } { \hat { \lambda } _ { n } } - \frac { 2 \log \left ( 1 - \frac { \lambda } { \lambda _ { N } } \right ) } { \lambda }$$

More generally, one may use the exact formula (C.15) derived in Appendix C.1, which expresses R N n essentially as the exponential of an integral.

We conclude this section with two comments.

Firstly, the explicit formula (4.20) means, as promised in Section 3.4, that we can now write down the interacting, extremal solutions (3.40) corresponding to any particular interacting basis in closed form:

̸

$$\mathcal { L } _ { \Delta } ^ { ( \Delta _ { n } ) } ( z ) = g _ { \Delta } ( z ) - \sum _ { n = 1 } ^ { \infty } \left [ \, \prod _ { \substack { m = 1 \\ m \neq n } } ^ { \infty } \left ( \frac { \lambda _ { \Delta } - \lambda _ { m } } { \lambda _ { n } - \lambda _ { m } } \right ) \, \right ] g _ { \Delta _ { n } } ( z ) \.$$

Moreover, any linear combination of such solutions (possibly with different choices of basis) is also a solution to locality.

Secondly, as we mentioned in Section 3.2.2 (item (b)), having this explicit expression for the interacting functional actions θ n (∆) crucially allows us to study the sequence θ n ( ˆ ∆ ) m appearing in the decomposition f n = ∑ m θ n ( ˆ ∆ ) ˆ m f m . In particular, in Appendix C.1,we use it to establish that this sequence has appropriate analyticity properties in m , which are then used to show that f n ∈ A .

## 5 Application: a case study in extremality

In this section we will apply our formalism to a concrete example. Specifically, we will construct a functional basis, and associated extremal solutions, which diagonalise the bulk locality constraints for a certain interacting CFT spectrum. We will do this in multiple ways, as a consistency check on our various results.

Before doing so, it is useful to have a general idea of what interacting functionals and their actions look like, which we turn to next.

## 5.1 Preamble: general expectations

We have established the existence of large families of interacting functionals, in terms of their kernels f n ( z ) and their actions θ n (∆) on blocks. We would like to understand how the large ∆ behaviour of θ n (∆) correlates with properties of the functional kernel. The functionals act on blocks as

$$\theta _ { n } ( \Delta ) = ( f _ { n }, \, g _ { \Delta } ^ { ( \tilde { \alpha } ) } ) = \int _ { 1 } ^ { \infty } \frac { \text{d} z } { \pi z ^ { 1 + \tilde { \alpha } } } \, f _ { n } ( z ) \, \mathcal { I } _ { z } g _ { \Delta } ( z ) \.$$

For large ∆, the integral is dominated by contributions around z = 1 and ∞ . In the neighbourhood of these points we can deform the contour to write

$$\theta _ { n } ( \Delta ) & = \int _ { - \infty } \frac { \mathrm d z } { \pi } \mathcal { I } _ { z } \left [ f _ { n } ( z ) e ^ { i \pi \frac { \Delta - 2 \tilde { a } } { 2 } } \right ] \, ( - z ) ^ { - 1 - \tilde { a } } \tilde { g } _ { \Delta } ( z ) - \int ^ { 1 } \frac { \mathrm d z } { \pi } \mathcal { I } _ { z } [ z ^ { - 1 - \tilde { a } } f _ { n } ( z ) ] \, g _ { \Delta } + \dots \, \\ \dots \dots \, \vdots \, \dots \, \dots \, \dots \, \dots \, \vdots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \end{array}$$

where the ellipsis represents contributions away from those points and ˜ g ∆ is defined in equation (A.5). Note that, in the above, the contour has been deformed so that the integrals shown run along the real axis along z &lt; 0 and z &lt; 1 respectively in the neighbourhoods of ∞ and 1. Let us now assume

$$f _ { n } ( z ) \underset { z \to 1 } { \sim } d _ { 1 } \, ( z - 1 ) ^ { \tilde { \alpha } + \frac { n _ { 1 } } { 2 } } + \text{analytic} \, \\ f _ { n } ( z ) \underset { z \to \infty } { \sim } - \frac { 1 } { z } + \frac { d _ { \infty } } { z ^ { 1 + \frac { n _ { \infty } } { 2 } } }$$

where we need η 1 , ∞ ≥ 0 in order that f n belong to the space A . Using the asymptotic expressions given in appendix (A.17) we can perform the integrals and obtain an expression for the functional action at large ∆. Defining the constants,

$$r _ { \infty } & \coloneqq \frac { 2 } { \pi } \left ( 1 + \tilde { \alpha } \right ) _ { \frac { \eta _ { \infty } } { 2 } } ^ { 2 } \sin [ \frac { \pi } { 2 } \eta _ { \infty } ] \, \\ r _ { 1 } & \coloneqq \frac { 2 } { \pi } \sin [ \frac { \pi } { 2 } \left ( 2 \tilde { \alpha } + \eta _ { 1 } \right ) ] \left ( 1 + \tilde { \alpha } \right ) _ { \frac { \eta _ { 1 } + 2 - d } { 2 } } \left ( 1 + \tilde { \alpha } \right ) _ { \frac { \eta _ { 1 } + 2 - d } { 2 } } \,,$$

then we obtain:

$$\theta _ { n } ( \Delta ) \underset { \Delta \to \infty } { \sim } \frac { \Gamma ( 1 + \tilde { \alpha } ) ^ { 2 } } { \pi ^ { 2 } ( \frac { \Delta } { 2 } ) ^ { 2 \tilde { \alpha } - \frac { 4 } { 2 } + 3 } } } \, \left ( 1 - \frac { d _ { \infty } r _ { \infty } } { ( \frac { \Delta } { 2 } ) ^ { \eta _ { \infty } } } \right ) \\ \times 2 \sin [ \frac { \pi } { 4 } ( \Delta - 2 \tilde { \alpha } - \gamma ^ { ( + ) } ( \Delta ) ] \cos [ \frac { \pi } { 4 } ( \Delta - 2 \tilde { \alpha } - \gamma ^ { ( - ) } ( \Delta ) ] \,, \ \ ( 5. 5 )$$

with

$$\frac { \gamma ^ { ( + ) } ( \Delta ) + \gamma ^ { ( - ) } ( \Delta ) } { 2 } = - \frac { d _ { \infty } r _ { \infty } } { ( \frac { \Delta } { 2 } ) ^ { \eta _ { \infty } } } \, \quad \frac { \gamma ^ { ( + ) } ( \Delta ) - \gamma ^ { ( - ) } ( \Delta ) } { 2 } \ = \frac { d _ { 1 } r _ { 1 } } { ( \frac { \Delta } { 2 } ) ^ { \eta _ { 1 } - \frac { d - 2 } { 2 } } } \. \quad \ ( 5. 6 )$$

This result shows that the behaviour of the functional kernel near z = 1 , ∞ directly correlates with the positions of the zeros at large ∆. In particular the zeros occur at ∆ = 2˜ + 2 n α n + γ n with

$$\gamma _ { n } \underset { n \to \infty } { \sim } - \frac { d _ { \infty } r _ { \infty } } { n ^ { \eta _ { \infty } } } + ( - 1 ) ^ { n } \, \frac { d _ { 1 } r _ { 1 } } { n ^ { \eta _ { 1 } + 1 - \frac { d } { 2 } } } \,.$$

If anomalous dimensions are to decay at large n then necessarily:

$$\eta _ { \infty } > 0 \,, \quad \eta _ { 1 } > \frac { d - 2 } { 2 } \,.$$

For d &lt; 2 we note that the earlier condition η 1 ≥ 0 dominates so that, for any basis, the alternating-sign piece of γ n must decay at least as n -1+ d/ 2 - i.e. n -1 2 / in the d = 1 case. This is in agreement with the analysis of Appendix C.4, and is the reason for this slightly stronger assumption in the d = 1 case (see footnote 14).

We now recall that the interacting functionals admit the following decomposition in terms of the GFF basis:

$$f _ { n } = \sum _ { m = 1 } ^ { \infty } b _ { n m } \, \hat { f } _ { m } \,, \quad \theta _ { n } ( \Delta ) = \sum _ { m = 1 } ^ { \infty } b _ { n m } \, \hat { \theta } _ { m } ( \Delta ) \.$$

We would like to understand how this can be consistent with the results above. The coefficients b nm in these decompositions will be fixed in terms of the interacting functional's zeros ∆ n . We will see that inputting the behaviour (5.7) of the zeros reproduces the right power-law behaviour (5.3) of the interacting functional kernels by summing up the GFF ones. In passing, we note that the direct analysis here is distinct, but morally equivalent to, that of Appendices C.3 and C.4, which instead passes through an integral representation for the functionals.

Let us set γ m as in (5.7). We want to check that the sum over GFF kernels produces the correct behaviour at z = 1 , ∞ given in (5.3). Since none of the individual GFF kernels has satisfies it, this behaviour must arise from the infinite sum, and it is sufficient to focus on its tail. We now note that demanding the right duality properties for large ∆ gives

$$\theta _ { n } ( \hat { \Delta } _ { m } + \gamma _ { m } ) = 0 \quad \Rightarrow \quad b _ { n m } = - \gamma _ { m } \hat { \theta } _ { n } ^ { \prime } ( \hat { \Delta } _ { m } ) + O ( \gamma _ { m } ^ { 2 } )$$

with

$$\theta _ { n } ^ { \prime } ( \hat { \Delta } _ { m } ) \underset { m \to \infty } { \sim } ( - 1 ) ^ { m } O ( m ^ { - 2 \tilde { \alpha } + \frac { d } { 2 } - 3 } ) \.$$

Zooming in on z → ∞ and z → 1 regions while doing the sum over large m with m/ √ z and m z √ -1 fixed respectively, we can use the asymptotic expressions (A.18) to find

$$& \sum _ { m } b _ { n m } \, \hat { f } _ { m } ( z ) \supsetneq \hat { f } _ { n } ( z ) \times \frac { d _ { \infty } } { z ^ { \frac { \eta _ { 1 } } { 2 } } } \, \\ & \sum _ { m } b _ { n m } \, \hat { f } _ { m } ( z ) \supsetneq \hat { f } _ { n } ( 1 ) \times 1 + d _ { 1 } ( z - 1 ) ^ { \tilde { \alpha } + \frac { \eta _ { 1 } } { 2 } }$$

precisely as expected.

After these preliminaries, let us move on to the analysis of a concrete set of bases relevant

for an actual CFT computation.

## 5.2 An extremal family

Our example originates from considering a family of 1d CFT 4-point correlators of an elementary field ϕ with scaling dimension ∆ ϕ . In the OPE of ϕ with itself, after the identity operator, we denote the lowest-dimension operator as ϕ 2 and its dimension as ∆ , which is 0 chosen as an input. Below we set

$$g \coloneqq \gamma _ { 0 } = \Delta _ { 0 } - 2 \Delta _ { \phi }$$

and consider g ≥ 0 for definiteness. Asking that the OPE coefficient λ ϕϕϕ 2 is maximised then determines the correlator to take the form [19,30]

$$z ^ { 2 \Delta _ { \phi } } \left \langle \phi ( \infty ) \phi ( 1 ) \phi ( z ) \phi ( 0 ) \right \rangle = 1 + \sum _ { n = 0 } ^ { \infty } a _ { n } \, G _ { \Delta _ { n } } ^ { \text{CFT} } \,,$$

where G CFT ∆ = z ∆ 2 F 1 (∆ ∆ 2∆ , , , z ) and

$$\Delta _ { n } & = 2 \Delta _ { \phi } + 2 n + \gamma _ { n } \quad \text{ with } \quad \gamma _ { n } = O ( n ^ { - 2 } ) \, \\ a _ { n } & = a _ { n } ^ { g g f } ( 1 + \delta _ { n } ) \quad \text{ with } \quad \delta _ { n } = O ( n ^ { - 2 } ) \.$$

The OPE contains only a discrete set of operators, denoted ( ϕ 2 ) n , and both their OPE coefficients a n ≡ ( λ ϕϕ ϕ ( 2 ) n ) 2 and their dimensions ∆ n are uniquely determined by the choice of ∆ 0 - and known numerically to high precision [20,33]. For g ≪ 1, both δ n and γ n are O g ( ) and hence the correlator approaches that of a bosonic Generalized Free Field (GFF). For finite g , the above corresponds to a family of deformations of the GFF that preserve the density of states in the ϕ × ϕ OPE. This family cannot be extended beyond g = 1, since ∆ = 1+2∆ 0 ϕ represents an absolute upper bound on the dimension of the first operator in the OPE [17]. As we approach g = 1, the correlator approaches that of a Generalised Free Fermion, which saturates this bound. In detail, as g → 1 we have ∆ n → 1 + 2∆ ϕ +2 n (the fermionic solution spectrum), but non-uniformly in 1 -g , i.e. such that the asymptotics (5.15) remain valid for all g &lt; 1.

For small g , the above CFT data matches that of a free scalar field Φ in AdS 2 (dual to the generalized free field ϕ ) with an additional quartic interaction ∝ g Φ . 4 Thus, in that regime at least, it makes sense to say that the above CFT data arises from a local bulk dual. 19 For instance, to a few orders in perturbation theory in g , it is possible to consider the form factor ⟨ Φ 2 | ϕϕ ⟩ . After stripping off the universal prefactors, this form factor is

19 Note that, for finite coupling g , this data is not described by a Φ 4 theory in AdS . 2 Specifically, the mismatch will occur at O g ( 4 ) where in Φ 4 theory we would have new states ∼ ϕ 4 appearing in the ϕ × ϕ OPE.

written:

$$F ^ { \Phi ^ { 2 } } ( z ) = \sum _ { n = 0 } ^ { \infty } c _ { n } \, g _ { \Delta _ { n } } ( z ) \,, \quad c _ { n } \coloneqq \lambda _ { \phi \phi ( \phi ^ { 2 } ) _ { n } } \, \mu _ { ( \phi ^ { 2 } ) _ { n } } ^ { \Phi ^ { 2 } } \mathcal { N } _ { \Delta _ { n } } ^ { - 1 }$$

where for convenience we continue to use the normalisation (2.3) of the locality blocks. In our previous work [15] we considered this form factor, using the GFF functional bases to solve for the coefficients c n to leading order in the coupling g .

Given that the non-perturbative (albeit numerical) CFT data is available, a natural question is whether this local Φ 2 operator can still be constructed at finite coupling g . Indeed, our construction of interacting functionals above answers this question affirmatively. Identifying ∆ = ∆ 0 as the fixed 'input' and ∆ n ≥ 1 as the 'sparse' spectrum dual to a basis of functionals, such form factors must be proportional to our extremal 'local block' soltuions:

$$F ^ { \Phi ^ { 2 } } = c _ { 0 } \, \mathcal { L } _ { \Delta _ { 0 } } ^ { \{ \Delta _ { n } \} } \.$$

These interacting local block solutions L { ∆ n } ∆ have been constructed explicitly in (3.40) and (4.25). They are guaranteed to be local as long as the dimensions ∆ n satisfy the required analyticity assumptions. There is much evidence for this: it is true in perturbation theory in the coupling g and, even for g = O (1), the γ n admit a systematic expansion in powers of 1 /λ ˆ n . In what follows we will take this as a working assumption. 20

The local blocks, and hence the form factor above, are explicitly known, so the easiest way to evaluate them is by directly using the explicit formula (4.25). Nonetheless, we would also like to use this example to check our whole formalism. In the following subsections, we will construct the interacting functional kernels, functional actions, and extremal solutions corresponding to this example - first perturbatively for small g and then numerically for finite g .

## 5.3 Interacting basis: perturbative computation

Let us first describe what happens in perturbation theory for g ≪ 1. To understand how locality constrains the form factor (5.16), the simplest approach is to use the GFF sum rules with ˜ = ∆ . Writing α ϕ c n = c (0) n + g c (1) n and expanding in g , we find (setting ∆ ˆ n ≡ 2∆ +2 ϕ n and assuming γ n = O g ( )):

$$\sum _ { m = 0 } ^ { \infty } c _ { m } \, \hat { \theta } _ { n } ( \Delta _ { m } ) = 0 \quad \Longleftrightarrow \quad \begin{cases} c _ { 0 } ^ { ( 0 ) } \, \hat { \theta } _ { n } ( \hat { \Delta } _ { 0 } ) + c _ { n } ^ { ( 0 ) } = 0, \\ g \, c _ { n } ^ { ( 1 ) } + \sum _ { m = 0 } ^ { \infty } c _ { m } ^ { ( 0 ) } \, \gamma _ { m } \, \hat { \theta } _ { n } ^ { \prime } ( \hat { \Delta } _ { m } ) = 0 \, \end{cases}$$

20 We are also assuming that F Φ 2 is sufficiently polynomially bounded, so that L { ∆ n } ∆ 0 is the unique such solution up to proportionality. Indeed, this is physically reasonable since it is the case in perturbation theory in g .

which uniquely fixes the BOE data to this order.

Now let us see how to derive the same results by constructing an interacting functional basis satisfying

$$\theta _ { n } ( \Delta _ { m } ) = \delta _ { n m } \.$$

These functionals may always be expanded in the GFF basis:

$$\theta _ { n } = \sum _ { m = 1 } ^ { \infty } b _ { n m } \, \hat { \theta } _ { m } \.$$

Imposing the duality relation (5.19) to O g ( ) we find

$$\theta _ { n } & = \hat { \theta } _ { n } - \sum _ { m = 1 } ^ { \infty } \gamma _ { m } \, \hat { \theta } _ { n } ^ { \, \prime } ( \hat { \Delta } _ { m } ) \, \hat { \theta } _ { m } + \dots \\ \Longrightarrow \quad f _ { n } & = \hat { f } _ { n } - \sum _ { m = 1 } ^ { \infty } \gamma _ { m } \, \hat { \theta } _ { n } ^ { \, \prime } ( \hat { \Delta } _ { m } ) \, \hat { f } _ { n } + \dots \.$$

We can also check whether this expression for the interacting functionals is consistent with our general formula (4.20). In perturbation theory that formula gives

̸

$$\theta _ { n } ( \lambda ) & = \hat { \theta } _ { n } ( \lambda ) \left [ 1 - ( \lambda _ { n } - \hat { \lambda } _ { n } ) \theta _ { n } ^ { \prime } ( \hat { \lambda } _ { n } ) \right ] + \sum _ { \substack { m = 1 \\ m \neq n } } ^ { \infty } ( \lambda _ { m } - \hat { \lambda } _ { m } ) \frac { ( \lambda - \hat { \lambda } _ { n } ) \hat { \theta } _ { n } ( \lambda ) } { ( \lambda - \hat { \lambda } _ { m } ) ( \hat { \lambda } _ { n } - \hat { \lambda } _ { m } ) } \\ & = \hat { \theta } _ { n } ( \lambda ) - \sum _ { m = 1 } ^ { \infty } ( \lambda _ { m } - \hat { \lambda } _ { m } ) \hat { \theta } _ { n } ^ { \prime } ( \hat { \lambda } _ { m } ) \theta _ { m } ( \lambda )$$

which indeed matches the above.

Continuing, the interacting basis sum rules applied to the interacting solution are trivialized:

$$\sum _ { \Delta } c _ { n } \, \theta _ { n } ( \Delta _ { n } ) = 0 \quad \Leftrightarrow \quad c _ { n } = - \theta _ { n } ( \Delta _ { 0 } ) \ c _ { 0 } \.$$

Plugging in (5.21) and expanding in g this gives:

$$c _ { n } ^ { ( 0 ) } + g \, c _ { n } ^ { ( 1 ) } = - c _ { 0 } \, \theta _ { n } ( \hat { \Delta } _ { 0 } ) + c _ { 0 } \sum _ { m = 0 } ^ { \infty } \gamma _ { m } \, \hat { \theta } _ { n } ^ { \, \prime } ( \hat { \Delta } _ { m } ) \, \hat { \theta } _ { m } ( \hat { \Delta } _ { 0 } )$$

which indeed agrees with the previous computation (5.18).

Further details can be worked out for special values of ∆ ϕ . For instance, for ∆ ϕ = 1 the

Figure 4: Functional actions. The figure shows a comparison between a GFF functional action and an interacting one, the latter corresponding to an extremal spectrum with γ 0 = 0 99 . (rescaled by a positive function of ∆ for clarity). The dashed lines at the bottom of the plot in green and red mark the functionals' zeros, i.e. the dimensions of the spectra of the corresponding solutions to crossing. For large ∆, as expected, the actions become identical.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-68d9ab8e-5139-4548-8eb9-5b8b07fecb60.jpg)

anomalous dimensions are given analytically by (see e.g. [19]):

$$\gamma _ { n } = \frac { \gamma _ { 0 } } { ( n + 1 ) ( 2 n + 1 ) } = \frac { 2 \gamma _ { 0 } } { \hat { \Delta } _ { n } ( \hat { \Delta } _ { n } - 1 ) } \.$$

In this case, as explained in appendix E, it is possible to perform the sum in (5.21) explicitly, finding

$$\frac { f _ { 1 } ( z ) - \hat { f } _ { 1 } ( z ) } { \gamma _ { 0 } } = \frac { z ^ { \frac { 3 } { 2 } } \sinh ^ { - 1 } ( \sqrt { z - 1 } ) } { 3 \sqrt { z - 1 } } - \frac { ( 2 4 z \log ( 2 ) + 6 ( 2 z + 1 ) \log ( z ) + 7 ) } { 7 7 } \. \quad ( 5. 2 6 )$$

As is manifest from this explicit expression, upon adding interactions, the functional kernels have now developed a branch cut running along ( -∞ , 1) - consistently with the general expectations of Section 5.1 and Appendix C.

## 5.4 Interacting bases: numerical computation

We will now discuss a method for a non-perturbative but numerical construction of the local Φ operator. This method provides a bottom-up consistency check for the above closed form 2 expressions for interacting functionals. The idea is as follows. Starting from the CFT data ∆ , we will construct interacting functionals as the largen N limit of combinations of GFF

functionals:

$$\theta _ { n } ^ { N } = \sum _ { m = 1 } ^ { N } b _ { n m } ^ { N } \, \hat { \theta } _ { m } \.$$

The duality conditions:

$$\theta _ { n } ^ { N } ( \Delta _ { m } ) = \delta _ { n m } \quad \Longleftrightarrow \quad \sum _ { m = 1 } ^ { N } b _ { n m } ^ { N } \, \hat { \theta } _ { m } ( \Delta _ { n } ) = \delta _ { n m } \.$$

provide finitely many linear constraints on the coefficients b N nm , which may be solved. As N increases, this leads to systematically better approximations to the exact interacting functional θ n . This limit is guaranteed to converge in the sense b N nm → b nm . 21 This follows from the special properties of Riesz bases, and would not be true for a general functional basis (i.e. derivatives).

In fact, it is possible to improve the convergence in N by subtracting off the following perturbative-like 'tail'

$$\theta _ { n } ^ { N } \rightarrow \theta _ { n } ^ { N } + T _ { n } ^ { N } \,, \quad T _ { n } ^ { N } \coloneqq \sum _ { m = N + 1 } ^ { \infty } \gamma _ { m } \, \theta _ { n } ^ { N \prime } ( \hat { \Delta } _ { m } ) \, \hat { \theta } _ { m } \.$$

The tail T N n approximates the contributions of the operators of large dimension. The addition of this tail introduces a small error in the duality conditions that were previously solved. To fix this, one may follow an iterative procedure, writing

$$\theta _ { n } ^ { N, k } = \sum _ { m = 1 } ^ { N } b _ { n m } ^ { N, k } \, \hat { \theta } _ { m } + T _ { n } ^ { N, k - 1 }$$

and solving for b N,k nm iteratively. In practice, two or three iterations suffice to find negligible changes in the coefficients. Once we have the coefficients b nm , then we can obtain not only the functional actions but also the interacting kernels f n = ∑ m b nm ˆ f m .

The results of this procedure are shown in Figures 4 and 5. The first figure plots the interacting functional action θ n (∆) while the second plots the functional kernels f n ( z ), for various choices of the coupling g (i.e. ∆ ). As an important cross-check, we find that: 0

$$\lim _ { N \to \infty } b _ { n m } ^ { N } = \theta _ { n } ( \hat { \Delta } _ { m } )$$

21 This can be seen by noting that the functional kernels f N n corresponding to (5.27) are the first N elements of a Riesz basis dual to { p n ≤ N , ˆ p n&gt;N } . By the analysis of Appendix B.2, the latter are the image of { ˆ p n } under a bounded, invertible operator M N = 1 -K N , with K N defined in (B.43). This operator converges in operator norm M N → M where M p (ˆ ) = n p n , and hence so does M -1 N → M -1 . Then we find b N m n = ⟨ ˆ f n | M -1 N | ˆ p m ⟩ → b nm = ⟨ ˆ f n | M -1 | ˆ p m ⟩ as required.

Figure 5: Interacting functional kernels with n = 1 (rescaled by the GFF functional kernel, which is just a power of z ). The figure shows the kernels for various choices of g , with the inset displaying a better view of the kernel for g = 0 2. . The black dashed curve corresponds to the analytic result given in equation (5.26). As g → 1, the curves approach z -1 2 (the red dashed curve), so that f (˜) α 1 → ˆ f (˜+ α 1 2 ) 1 (with ˜ = ∆ α ϕ ). The plots were done for the particular choice ∆ ϕ = 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-47494166-a707-480c-9540-91af69c31970.jpg)

with the RHS computed using the exact formula (4.20) (or in practice the approximations (4.22) and (4.23), which ultimately correspond to the same approximation scheme used above). Incidentally, the plots demonstrate that, as the coupling is increased, the interacting functionals approach the GFF basis with ˜ = ∆ α ϕ +1 2 as they should, since the CFT data / is approaching the Generalized Free Fermion solution with spectrum ∆ n = 2(∆ +1 2)+2 ϕ / n .

With the functional actions well under control, we can compute the BOE data. For instance, the construction of the interacting functionals above gives

$$c _ { n } = - c _ { 0 } \, \theta _ { n } ^ { N } ( \Delta _ { 0 } ) + O ( \gamma _ { N } ^ { 2 } ) \quad \text{ for } \ 1 \leq n \leq N$$

where quadratic accuracy follows from our inclusion of the tail. It is also possible to predict the BOE data to arbitrarily high n in terms of the lower n data, simply by using the highern GFF functional sum rules to get

$$c _ { n } = - \sum _ { m = 0 } ^ { N } \frac { \hat { \theta } _ { n } ( \Delta _ { m } ) } { \hat { \theta } _ { n } ( \Delta _ { m } ) } + O ( \gamma _ { N } ) \quad \text{ for all } \ n > N$$

This in turn can be improved to O γ ( 2 N ) by making the replacement

$$\hat { \theta } _ { n } \rightarrow \hat { \theta } _ { n } - \sum _ { m = N + 1 } ^ { \infty } \gamma _ { m } \, \hat { \theta } _ { n } ^ { \prime } ( \hat { \Delta } _ { m } ) \, \hat { \theta } _ { m }$$

Figure 6: Form factors. In the figure we show the interacting form factors corresponding to various values of the coupling g . As the coupling increases the form factor interpolates between z ∆ + ϕ 1 2 (the red dashed sloping line) and z ∆ ϕ (the horizontal dashed curve). The form factors were computed by representing them as sums of local blocks, with the corresponding coefficients µλ computed as discussed in the main text. The result is local by construction, and in particular it shows no singularity for z ≥ 1. As a consistency check we show in the inset that in the Euclidean region ( z ≤ 1) this representation agrees with a direct compution of the form factor as sum of blocks. As we include more terms in that sum (shown as successive dashed curves), the result nicely approaches the local block result in the entire Euclidean region.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-befa2ab4-2d0d-4c9d-a4ab-cd6cfb94bf91.jpg)

everywhere in the above formula. The outcome of this procedure is to improve our estimates of the BOE data, such that the error in all quantities is now O γ ( 2 N ).

With the BOE data in hand, we can reconstruct the local form factor. As noted above, this form factor can be identified with the interacting local block L { ∆ n } ∆ 0 . We plot it in Figure 6. The result contains rather rich physics: for instance, by scaling the coupling g → 1 and z →∞ at the same time, one probes the flat space RG flow for the bulk QFT (which happens well above the AdS scale), ending in a free fermion CFT in the IR [33].

## 6 Discussion and outlook

## 6.1 Extremal solutions vs bases

In this paper we described a natural language for the locality bootstrap problem. Candidate form factors can be mapped to hyperfunctions in a suitable space HF , with dual space A whose elements are functional kernels. This language makes precise the question of

completeness of a set of functionals { f n } , which had to be tackled in previous works [17, 19,30] from a more 'bottom-up' perspective: completeness simply means that the set { f n } is complete in A , i.e. that any element of A is a limit of linear combinations of the f n . We constructed a large family of such complete sets by asking for the f n to form Schauder bases for A dual to distinguished bases in HF made up of blocks with specific sets of scaling dimensions { ∆ n } .

Such complete bases are associated with canonical solutions to the locality problem, 'local blocks' L { ∆ n } ∆ given in equation (4.25). These are local form factors which can be thought of as 'extremal', in the sense that their BOE contains the minimal amount of blocks compatible with locality. 22 We have thus established a link between bases and extremal solutions . We expect this link to be quite general, and hold across a variety of distinct bootstrap problems. It seems therefore that characterising bases for the relevant spaces is a useful path for an understanding of bootstrap equations.

In our construction, an intermediate step involved understanding Riesz bases for L 2 . This is a classical problem in functional analysis. For example, for the space L 2 ( -π, π ), the exponential functions { e int } n ∈ Z form an orthonormal basis, but one may consider different 'spectra' { λ n } and ask whether { e iλ n t } n ∈ Z still forms a Riesz basis. The solutions of that problem have been classified [34] (in Russian; see the recent review [35]), with a sufficient condition being sup n | λ n -n | &lt; 1 4 (a famous result of Kadets [36]). In our case, we dealt instead with the sets {I z g ∆ n ( z ) } and established sufficient conditions on the spectrum ∆ n for them to form a Riesz basis for L 2 µ (1 , ∞ ):

$$| \Delta _ { n } - ( 2 \tilde { \alpha } + 2 n ) | \stackrel { n \to \infty } { \lesssim } n ^ { - \epsilon } \, \quad | \tilde { \alpha } | < \frac { 1 } { 2 }, \, \epsilon > 0 \.$$

Importantly, in our setting, in order to for these to uplift to Schauder bases for the spaces A and HF , it was also necessary that the sequence { ∆ n } has certain analyticity properties, adding a significant further constraint.

Let us mention an interesting set of bases that do not satisfy our assumptions on { ∆ n } , but nevertheless lead to functionals satisfying our properties 1) and 2) in the introduction. These 'Zettl bases' [37,38] rely on the fact that the block discontinuities I z g ∆ are eigenfunctions of the Casimir operator (3.7). This operator can be made Hermitian on L 2 µ by choosing appropriate boundary conditions at z = 1 , ∞ , in such a way that special sets {I z g ∆ n } remain allowed eigenfunctions. The net result is a continuous family of sequences { ∆ n } (labeled by ∆ 1 and computable numerically) such that {I z g ∆ n } form orthonormal bases - thus with dual functional kernels trivially given by the blocks themselves. The reason we did not find these basis in our construction is that such sequences have ∆ n -2 n = O (1 / log( n )) (not satisfying our assumptions), which does not seem to be immediately physically relevant for bootstrap applications.

22 For a fixed polynomial growth at z = ∞ .

## 6.2 Implications for numerics

The results of this paper have interesting implications for the numerical bootstrap program in d = 1 (and probably for higher d as well, but we focus the following discussion on that case). One question of interest is: in demanding that the CFT data supports the existence of local bulk operators, does one get stronger bounds? In other words, does combining crossing with locality lead to stronger constraints?

To get a handle on this question, first we note that, to understand the bootstrap bounds, it is sufficient to focus on the construction of the extremal solutions of crossing that saturate them. 23 Thus, to answer the question, one only has to consider extremal solutions of crossing, and ask whether it is possible to solve the locality constraints using them as inputs.

In Section 5, we saw for a particular family of extremal solutions that it was indeed possible to solve the locality constraints, and we explicitly reconstructed a local form factor associated to it, corresponding to some operator Φ 2 in the bulk of AdS. In fact, we expect this to be possible for all extremal solutions of crossing, which always consist of finite sets of fixed 'input' CFT data (specifying which bound is to be saturated), completed by an infinite collection of 'double-trace' operators, which are dynamically determined by crossing and optimality. There are both analytical arguments and ample numerical evidence that the scaling dimensions of such operators must always asymptote to either ∆ n = 2∆ ϕ +2 n or ∆ n = 1 + 2∆ ϕ +2 n and meet our other assumptions on 'sparse' spectra: being analytic and decaying as some power for large n . The results of this paper then imply that we can always choose a basis of locality functionals that will be dual to such spectra, and hence reconstruct a local form factor, which will in practice always be a finite sum of the 'local block' solutions for that basis.

Thus we find that generic extremal solutions are compatible with the existence of local bulk operators. Intuitively what has happened is that, although bulk locality introduces infinitely many new constraints, it also introduces infinitely many new pieces of data to be fixed: the BOE coefficients. Non-trivially, these two infinities line up on the nose, thanks to the fact that the bases we were able to construct for locality are precisely adapted to those spectra that saturate CFT bounds.

The good news is that we can still learn something from locality and obtain additional bounds on top of crossing, by including some additional input. In fact, since we already have a tight matching of infinite constraints and degrees of freedom, any finite mismatch will suffice to create the required tension. The simplest thing is just to lower the parameter ˜. α In particular, lowering this value by one unit corresponds to adding one extra locality constraint. This can be motivated if, in some case, one already has a notion of what the correct value should be. In a numerical setup, bounds would thus have to become tighter. For instance, for the family of examples we considered in Section 5, choosing ˜ = ∆ α ϕ -1 would lead to the absence of a solution (contingent on fixing an overall normalisation of the

23 A detailed characterization and study of such extremal solutions will be given elsewhere. For the purposes of this section, it is sufficient to think of them as the solutions of crossing saturating a bound on the CFT data.

bulk field, e.g. fixing one of the c n coefficients). The upper bound on the OPE coefficient a 0 in (5.14) would thus have to decrease, as the extremal solution under consideration is now ruled out. Perhaps more interestingly, special choices of bulk operators such as the stresstensor or conserved currents already come with added constraints - e.g. those developed for the stress-tensor in [26]. This makes sense, since such operators know about the bulk Lagrangian and cannot be compatible with arbitrary boundary data.

## 6.3 Outlook

To conclude, let us briefly mention a few of the interesting research directions in this subject.

Firstly, it would be very interesting to apply our machinery to other bootstrap problems, particularly the crossing equation. In particular, one of our most significant results is the characterisation of the space of functional actions in Section 4 and the explicit product formula (4.20) expressing the functional actions in terms of their zeros. We believe that this approach may similarly give explicit formulae for the functional actions in the crossing case - perhaps yielding greater analytic control over the interacting extremal solutions there [20]. In general, working with the space of functional actions is appealing since it bypasses the difficult construction of functional kernels and the laborious computation of their action. In fact, similar logic already had success in the past, although it was not formulated in our language [39-42].

Secondly, we sketched in [15] the eventual generalisations of our work for bulk theories with gauge and gravitational symmetries. An appropriate choice of dressing will enforce specific non-localities in certain form factors, which will generally be distributions:

$$F ^ { ( \tilde { \alpha } ) } = H ^ { \text{dressing} } \.$$

Our formulation of functionals as Schauder bases immediately provides a way to convert this 'source' term to source terms in the locality sum rules. Decomposing

$$H ^ { \text{dressing} } = \sum _ { n } ( f _ { n }, H ^ { \text{dressing} } ) \, g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) }$$

the sum rules become

$$\sum _ { \Delta } c _ { \Delta } \, \theta _ { n } ( \Delta ) = ( f _ { n }, H ^ { \text{dressing} } ) \.$$

For instance it would be very interesting to see if this prescription reproduces bulk radial Wilson-line dressings for charged operators.

Finally, our work has been concerned with characterizing locality around the AdS vacuum. However, in gravitational theories with a large N limit, the bulk Hilbert space splits into various superselection sectors occupied, e.g., by black holes states. It would be interesting to try to formulate locality constraints in this context and make contact with

the proposals of [43,44].

## Acknowledgements

This work was co-funded by the European Union (ERC, FUNBOOTS, project number 101043588). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Research Council. Neither the European Union nor the granting authority can be held responsible for them. NL was supported by the Institut Philippe Meyer at the Ecole Normale Sup´rieure in Paris. ´ e We acknowledge discussions with K. Ghosh, P. Kravchuk, M. Meineri and S. Rychkov.

## A Compendium of useful formulae

For completeness and convenience, we recall in this section the results of our previous work [15] as well as various useful results. The locality constraint is equivalent to the following sum rules on the BOE data:

$$\sum _ { \Delta } c _ { \Delta } \hat { \theta } _ { n } ( \Delta ) = 0 \, \quad n \geq 1$$

with c ∆ defined in (2.3). The functional actions are given as

$$\hat { \theta } _ { n } ( \Delta ) & = ( \hat { f } _ { n }, g _ { \Delta } ^ { ( \tilde { \alpha } ) } ) = \int _ { 1 } ^ { \infty } \frac { \mathrm d z } { \pi } \, \hat { f } _ { n } ( z ) \frac { \mathcal { I } _ { z } g _ { \Delta } ( z ) } { z ^ { 1 + \tilde { \alpha } } } \\ & = \langle \hat { h } _ { n } | \mathcal { I } _ { z } g _ { \Delta } \rangle = \int _ { 1 } ^ { \infty } \mathrm d z \, \mu ( z ) \, \hat { h } _ { n } ( z ) \, \mathcal { I } _ { z } g _ { \Delta } ( z )$$

with the measure

$$\mu ( z ) = \frac { 1 } { z ^ { 2 } } \left ( \frac { z } { z - 1 } \right ) ^ { 1 - \frac { d } { 2 } } \,.$$

To derive the functional action, a useful representation is obtained by deforming the contour of integration to wrap the cut running along the negative real axis,

$$\theta _ { n } ( \Delta ) = \sin [ \frac { \pi } { 2 } ( \Delta - 2 \tilde { \alpha } ) ] \int _ { - \infty } ^ { 0 } \frac { \mathrm d z } { \pi } f _ { n } ( z ) ( - z ) ^ { - 1 - \tilde { \alpha } } \tilde { g } _ { \Delta } ( z )$$

with

$$\tilde { g } _ { \Delta } ( z ) = \mathcal { N } _ { \Delta } \, \left ( \frac { z } { z - 1 } \right ) ^ { \frac { \Delta } { 2 } } \, _ { 2 } F _ { 1 } \left ( \frac { \Delta } { 2 }, \frac { \Delta + 2 - d } { 2 }, \Delta + 1 - \frac { d } { 2 } ; \frac { z } { z - 1 } \right ) \.$$

The ˆ h n ∈ L 2 are given by

$$\hat { h } _ { n } ( z ) = i ^ { \tilde { \alpha } } ( \hat { f } _ { n } ) = \frac { \hat { f } _ { n } ( z ) } { \pi \mu ( z ) z ^ { 1 + \tilde { \alpha } } } \,,$$

and the functional kernels f n are polynomials of degree n in 1 /z without constant term, given explicitly by

$$\hat { f } _ { n } ( z ) = \pi \frac { ( - 1 ) ^ { n - 1 } } { ( n - 1 )! } ( \hat { \Delta } _ { n } - \frac { d } { 2 } ) \Gamma ( 2 \tilde { \alpha } + n + 1 - \frac { d } { 2 } ) \\ \times \frac { 1 } { z } \, _ { 3 } \tilde { F } _ { 2 } ( 1, 1 - n, 2 \tilde { \alpha } + n + 1 - \frac { d } { 2 } ; 1 + \tilde { \alpha }, 1 + \tilde { \alpha } ; \frac { 1 } { z } ) \.$$

Note that ˆ h n solves an inhomogeneous Casimir equation [15],

$$\left [ C _ { z } - \lambda _ { \hat { \Delta } _ { n } } \right ] \hat { h } _ { n } ( z ) = \frac { A _ { n } } { \mu ( z ) z ^ { 1 + \hat { \alpha } } } \,, \quad A _ { n } = - \frac { 4 \, \Gamma ( \hat { \Delta } _ { n } + 1 - \frac { d } { 2 } ) } { \Gamma ( \tilde { \alpha } ) ^ { 2 } \Gamma ( n ) \left ( 1 + \frac { d } { 2 } - \hat { \Delta } _ { n } \right ) _ { n - 1 } } \,, \quad \quad ( A. 8 )$$

with

$$C _ { z } \coloneqq 4 ( 1 - z ) ^ { 1 - d / 2 } z ^ { 1 + d / 2 } \partial _ { z } [ ( 1 - z ) ^ { d / 2 } z ^ { 1 - d / 2 } \partial _ { z } ] \.$$

The resulting functional actions are:

$$\hat { \theta } _ { n } ( \Delta ) = \frac { 2 } { \pi } \, \frac { \left ( \frac { \hat { \Delta } _ { n } - 2 \tilde { \alpha } } { 2 } \right ) _ { 1 - \frac { d } { 2 } + 2 \tilde { \alpha } } } { \left ( \frac { \Delta - 2 \tilde { \alpha } } { 2 } \right ) _ { 1 - \frac { d } { 2 } + 2 \tilde { \alpha } } } \frac { ( 2 \hat { \Delta } _ { n } - d ) \, \sin \left [ \frac { \pi } { 2 } ( \Delta - \hat { \Delta } _ { n } ) \right ] } { ( \hat { \Delta } _ { n } + \Delta - d ) ( \Delta - \hat { \Delta } _ { n } ) }$$

with ∆ ˆ n = 2˜ + 2 α n . In particular they satisfy the duality conditions

$$\hat { \theta } _ { n } ( \hat { \Delta } _ { m } ) = \delta _ { n, m } \, \quad n, m \geq 1 \.$$

They can be rewritten in a way that makes their invariance under ∆ ↔ d -∆ manifest:

$$\hat { \theta } _ { n } ( \Delta ) = \frac { 2 \, ( - 1 ) ^ { n } \, ( 2 \hat { \Delta } _ { n } - d ) } { ( \Delta + \hat { \Delta } _ { n } - d ) ( \Delta - \hat { \Delta } _ { n } ) } \, \frac { \left ( \frac { \hat { \Delta } _ { n } - 2 \tilde { \alpha } } { 2 } \right ) _ { 1 - \frac { d } { 2 } + 2 \tilde { \alpha } } } { \Gamma \left ( \frac { 2 + 2 \tilde { \alpha } - \Delta } { 2 } \right ) \Gamma \left ( \frac { 2 + 2 \tilde { \alpha } - d + \Delta } { 2 } \right ) } \.$$

The actions satisfy the asymptotics:

$$| \hat { \theta } _ { n } ( \Delta ) | \underset { n \to \infty } { = } O ( \Delta ^ { - 2 \tilde { \alpha } + \frac { d } { 2 } - 3 } ) \\ | \hat { \theta } _ { n } ( \Delta ) | \underset { n \to \infty } { = } O ( n ^ { 2 \tilde { \alpha } - \frac { d } { 2 } } )$$

The conformal blocks and the functional kernels have the following asymptotic forms for z ∈ C \ (1 , ∞ ).

$$g _ { \Delta } ( z ) \underset { \Delta \to \infty } { \sim } \frac { 2 ^ { 3 } _ { \frac { 3 } { 2 } } \Delta ^ { \frac { d - 3 } { 2 } } } { \sqrt { \pi } } \, \frac { [ 1 - r ( z ) ] ^ { 1 - \frac { d } { 2 } } } { \sqrt { 1 + r ( z ) } } \, r ( z ) ^ { \frac { \Delta } { 2 } } \, \\ \frac { \hat { f } _ { n } ( z ) } { z ^ { 1 + \tilde { \alpha } } } \underset { n \to \infty } { \sim } \frac { d r } { d z } \, \frac { \sqrt { \pi } } { 2 ^ { \frac { 3 } { 2 } } \left ( \hat { \Delta } _ { n } \right ) ^ { \frac { d - 3 } { 2 } } } \, \frac { \sqrt { 1 + r ( z ) } } { [ 1 - r ( z ) ] ^ { 1 - \frac { d } { 2 } } } \, r ( z ) ^ { - 1 - \tilde { \alpha } - n }$$

where

$$r ( z ) = \frac { 1 - \sqrt { 1 - z } } { 1 + \sqrt { 1 - z } } \.$$

We also have with θ = 1 i log ( 1+ i √ z -1 1 -i √ z -1 ) :

$$\mathcal { I } _ { z \geq 1 } \, g _ { \Delta } ( z ) \underset { \Delta \to \infty } { \sim } \frac { \left ( \frac { \Delta } { 2 } \right ) ^ { \frac { d - 3 } { 2 } } } { \sqrt { \pi } } \, \frac { [ \sin ( \theta / 2 ) ] ^ { \frac { 1 - d } { 2 } } } { \sqrt { \cos ( \theta / 2 ) } } \, \sin \left [ \frac { \Delta } { 2 } \theta + \frac { 1 } { 4 } [ \pi ( d - 1 ) - d \theta ] \right ] \.$$

Other useful asymptotic expressions are given by:

$$g _ { \Delta } ( z ) \underset { \overset { z \to 1 } { \overset { 1 } { \zeta } \to \infty } } { \sim } \frac { 2 } { \pi } \frac { \Gamma ( \frac { \Delta } { 2 } ) } { \Gamma ( \frac { \Delta + 2 - d } { 2 } ) } \left ( 1 - z \right ) ^ { \frac { 2 - d } { 4 } } K _ { \frac { d - 2 } { 2 } } ( \Delta \sqrt { 1 - z } ) \,, \\ \tilde { g } _ { \Delta } ( z ) \underset { \overset { z \to - \infty } { \overset { 1 } { \zeta } \to \infty } } { \sim } \frac { 2 } { \pi } \frac { \Gamma ( \frac { \Delta } { 2 } ) } { \Gamma ( \frac { \Delta + 2 - d } { 2 } ) } K _ { 0 } ( \Delta / \sqrt { - z } ) \,,$$

and similarly for the kernels

$$\hat { f } _ { n } ( z ) \underset { z \to 1 } { \sim } \frac { 2 \pi ( - 1 ) ^ { n - 1 } n ^ { 2 \tilde { \alpha } - \frac { d } { 2 } } } { \Gamma ( \tilde { \alpha } ) ^ { 2 } } + 2 \pi n ^ { 2 - \frac { d } { 2 } } ( z - 1 ) ^ { \frac { d - 2 } { 4 } } J _ { 1 - \frac { d } { 2 } } ( 2 n \sqrt { z - 1 } ) \,, \\ \hat { f } _ { n } ( z ) \underset { z \sim n ^ { 2 } } { \sim } \frac { 2 \pi } { z } \, ( - 1 ) ^ { n - 1 } n ^ { 2 \tilde { \alpha } + 2 - \frac { d } { 2 } } \, _ { 1 } \tilde { F } _ { 2 } ( 1, 1 + \tilde { \alpha }, 1 + \tilde { \alpha }, - n ^ { 2 } / z )$$

All asymptotic expressions can be derived from the fact that the blocks and functionals satisfy a Casimir equation, cf. (3.6) and (A.8).

## B Riesz bases

This appendix contains the proofs that the GFF functionals and the interacting functionals form Riesz bases for L 2 . Before beginning, we point out that a sequence of functions being a Riesz basis depends on the choice of normalisation: their norms need to be bounded. Thus

in this section it will be convenient to use the normalisation

$$p _ { \Delta } ( z ) \coloneqq \sqrt { \Delta - \frac { d } { 2 } } \, \left ( \frac { \Delta } { 2 } \right ) _ { 1 - \frac { d } { 2 } } \, \mathcal { I } _ { z } \, g _ { \Delta } ( z )$$

and define

$$p _ { n } \coloneqq p _ { \Delta _ { n } } \,, \quad \hat { p } _ { n } \coloneqq p _ { \hat { \Delta } _ { n } } \,, \quad \hat { \Delta } _ { n } = 2 \tilde { \alpha } + 2 n.$$

where for p n the particular implied sequence of ∆ n will be clear from context. Similarly we have to modify the normalisation of the dual elements h n . Overall there is a mismatch in normalisation conventions with the rest of the paper, which can be summarised as:

$$p _ { n } = \nu _ { n } \, p _ { n } ^ { \text{above} } \, \quad h _ { n } = \nu _ { n } ^ { - 1 } \, h _ { n } ^ { \text{above} } \quad \text{ with } \ \nu _ { n } = \sqrt { \Delta _ { n } - \frac { d } { 2 } } \, \left ( \frac { \Delta _ { n } } { 2 } \right ) _ { 1 - \frac { d } { 2 } } \.$$

## B.1 GFF Riesz bases

In this section we will prove property (a) for the GFF 'double trace' spectrum in Section 3.2, i.e. that the set { ˆ p n } forms a Riesz basis for L 2 . To achieve this, we just need to show there is a bounded, invertible, linear map:

↦

$$T \colon L ^ { 2 } & \to L ^ { 2 } \\ e _ { n } & \mapsto \hat { p } _ { n }$$

where e n is some orthonormal basis (any choice is equivalent). We remind the reader that our L 2 space is always defined with the non-trivial measure

$$\langle f | g \rangle \coloneqq \int _ { 1 } ^ { \infty } \frac { \mathrm d z } { z ^ { 2 } } \left ( \frac { z } { z - 1 } \right ) ^ { 1 - d / 2 } f ( z ) \, g ( z ) \.$$

Invertibility. Our approach for demonstrating invertibility of the putative map (B.4) is to explicitly construct an inverse. This is easy since we already explicitly know that the elements ˆ h n that are dual to the GFF blocks ˆ p n (as reviewed in Section A):

$$\langle \hat { h } _ { m } | \hat { p } _ { n } \rangle = \delta _ { m n } \, \quad m, n \geq 1 \.$$

Hence we may define a new map S as

$$S \colon L ^ { 2 } \rightarrow L ^ { 2 }$$

↦

$$e _ { n } & \mapsto \hat { h } _ { n } \.$$

Assuming we can prove that T and S define bounded operators, then S † = T -1 , since the duality relation (B.6) gives

$$S ^ { \dagger } T e _ { n } = e _ { n } \, \quad n \geq 1 \.$$

More precisely, this only shows that S † is a left-inverse of T . To conclude the proof of invertibility we need to show that S † is one-to-one (or equivalently that T is surjective). The trick will be that the GFF kernels ˆ h n are (up to some measure factors) polynomials in 1 /z , so can be related to known bases by taking finite linear combinations.

Let us see this in detail. We must show that if S g † = 0 for g ∈ L 2 , then g = 0. For any n ≥ 1, we have

$$\langle \hat { h } _ { n } | g \rangle = \langle e _ { n } | S ^ { \dagger } g \rangle = 0 \.$$

But the GFF functional kernels take the special form

$$\hat { h } _ { n } ( z ) = \nu _ { n } ^ { - 1 } \, z ^ { - \tilde { \alpha } } \, \left ( \frac { z } { z - 1 } \right ) ^ { \frac { d } { 2 } - 1 } \, z \hat { f } _ { n } ( z )$$

where zf ˆ ( n z ) are polynomials of degree ( n -1) in 1 /z . Hence it follows that

$$\langle z ^ { - \tilde { \alpha } } \, ( \frac { z } { z - 1 } ) ^ { \frac { d } { 2 } - 1 } \, j _ { n } ( 1 / z ) | g \rangle = 0$$

for any polynomials j n . Now let us take j n (1 /z ) ∝ J (1 -d/ , 2 2˜) α n -1 ( 2 z -1) Jacobi polynomials. These form a complete orthonormal basis for a different space, ˜ L 2 :

$$\delta _ { n, m } = \langle j _ { n } | j _ { m } \rangle _ { \bar { \mu } } = \int _ { 1 } ^ { \infty } \mathrm d z \, \tilde { \mu } ( z ) \, j _ { n } ( 1 / z ) j _ { m } ( 1 / z ) \, \quad \tilde { \mu } ( z ) = \frac { 1 } { z ^ { 2 + 2 \tilde { a } } } \left ( \frac { z } { z - 1 } \right ) ^ { \frac { d } { 2 } - 1 } \. \quad \text{(B.13)}$$

We now reinterpret equation (B.12) as

$$\langle j _ { n } | \tilde { g } \rangle _ { \tilde { \mu } } = 0 \,, \quad \tilde { g } ( z ) = g ( z ) \, z ^ { \tilde { \alpha } } \, \left ( \frac { z } { z - 1 } \right ) ^ { 1 - \frac { d } { 2 } } \,.$$

Since j n are complete in ˜ L 2 then g = ˜ = 0 as long as ˜ g g ∈ ˜ L 2 , but this is true since || ˜ g || 2 ˜ µ = || g || 2 as is easily checked. We conclude that S ∗ is 1-1.

↦

Schur's test. We want to show that the operators S and T defined above are well-defined, bounded operators. Note that in general if M : e n → q n for e n an orthonormal basis and q n some elements of L 2 , then M M , ∗ and M M ∗ are of course well-defined on finite linear combinations of the e n , and they satisfy

$$\| M x \| ^ { 2 } = \langle x | M ^ { \dagger } M x \rangle \leq \| x \| \, \| M ^ { \dagger } M x \| \leq \| M ^ { \dagger } M \| \, \| x \| ^ { 2 }$$

where we have used the Cauchy-Schwarz inequality and defined ∥ M M † ∥ = sup x ∥ M Mx † ∥ as usual (with x running over finite combinations of e n ).

We wish to extend M to a linear operator on the full L 2 space by the definition M (lim x n ) = lim M x ( n ), and similarly for M M † . To show this is possible, it is sufficient to apply Schur's test to the matrix elements M 2 mn := ⟨ e m | M M e † | n ⟩ = ⟨ q m | q n ⟩ . This classic result states that, if there exist positive sequences β n , α n and numbers λ , µ such that

$$\sum _ { n } | M _ { m n } ^ { 2 } | \, \beta _ { n } \leq \lambda \, \alpha _ { m } \, \quad \sum _ { m } | M _ { m n } ^ { 2 } | \, \alpha _ { m } \leq \mu \, \beta _ { n }$$

then M M † is bounded with ∥ M M † ∥ ≤ √ λµ , and hence is a well-defined map on all of L 2 . It follows from (B.15) that both M and M † are also well-defined and bounded.

Below we will apply this strategy to prove first boundedness of S and then T .

## Boundedness of S . In this case we have

$$\langle e _ { n } | S ^ { * } S | e _ { m } \rangle = \langle \hat { p } _ { n } | \hat { p } _ { m } \rangle$$

so we must study the latter matrix of inner products. In doing so it will become clear why we chose the particular L 2 measure (B.5). The functions ˆ p n are eigenfunctions of the conformal quadratic Casimir operator defined by (A.9),

$$\left ( C _ { z } - \lambda _ { \hat { \Delta } _ { n } } \right ) \hat { p } _ { n } ( z ) = 0 \.$$

This L 2 measure (B.5) is precisely the one for which C z is 'almost self-adjoint', i.e. selfadjoint up to boundary terms:

$$\langle f \, | \, C _ { z } g \rangle - \langle C _ { z } f \, | \, g \rangle = - 4 \left [ ( z - 1 ) ^ { d / 2 } z ^ { 1 - d / 2 } \, \left ( f \, \partial _ { z } g - g \, \partial _ { z } f \right ) \right ] _ { 1 } ^ { \infty } \,,$$

by simply integrating by parts. Taking ( f, g ) = (ˆ p , p n ˆ m ) and using (B.18) leads to a standard trick for evaluating the desired matrix elements:

̸

$$\langle \hat { p } _ { n } | \hat { p } _ { m } \rangle = \begin{cases} \frac { [ \psi _ { 1 } ( 1 - \hat { \frac { \Delta } { 2 } } _ { n } ) - \psi _ { 1 } ( 1 - \frac { d - \hat { \Delta } _ { n } } { 2 } ) ] \sin ^ { 2 } [ \frac { \pi } { 2 } \hat { \Delta } _ { n } ] } } { \pi ^ { 2 } } \sim 1 + O ( \frac { 1 } { n } ) & n = m \, \\ \frac { \sqrt { \hat { \Delta } _ { n } - d / 2 } \sqrt { \hat { \Delta } _ { m - d / 2 } } } { \hat { \Delta } _ { m + \hat { \Delta } _ { n } - d } } \frac { 4 \sin [ \frac { \pi } { 2 } \hat { \Delta } _ { n } ] \sin [ \frac { \pi } { 2 } \hat { \Delta } _ { m } ] } { \pi ^ { 2 } } & n \neq m \, \\ \times \frac { [ \psi _ { 0 } ( 1 - \hat { \Delta } _ { n } } _ { 2 } ) + \psi _ { 0 } ( 1 - \frac { d - \hat { \Delta } _ { n } } { 2 } ) - \psi _ { 0 } ( 1 - \frac { d - \hat { \Delta } _ { m } } { 2 } ) - \psi _ { 0 } ( 1 - \frac { d - \hat { \Delta } _ { m } } { 2 } ) } } { \hat { \Delta } _ { m } - \hat { \Delta } _ { n } } \end{cases}$$

where we recall that that ∆ ˆ n = 2˜ + 2 α n , and ψ p denotes the polygamma function of order p .

Let us apply Schur's test (B.16) with β n = α n = 1. Approximating ∑ m → n ∫ ∞ 1 /n dx

with m = xn , we find

$$\sum _ { m } | \langle \hat { p } _ { n } | \hat { p } _ { m } \rangle | \sim \int _ { 1 / n } ^ { \infty } d x \frac { \sqrt { x } \log x } { 1 - x ^ { 2 } } \lesssim 1 \.$$

Hence the test is satisfied and S is bounded.

Boundedness of T . To derive boundedness of T we must work a bit harder, essentially because the relevant matrix elements ⟨ ˆ h m | ˆ h n ⟩ are difficult to compute in closed form. Firstly, let us define

$$d _ { n } ( z ) \equiv ( - 1 ) ^ { n + 1 } \left [ \hat { h } _ { n } ( z ) - \hat { p } _ { n } ( z ) \right ] \.$$

We claim based on many examples that this quantity and its derivative have definite sign:

$$d _ { n } ( z ) \geq 0 \,, \quad \partial _ { z } d _ { n } ( z ) \geq 0 \quad \text{ for } \ z > 1 \.$$

We also claim that, for large n and fixed z ,

$$d _ { n } ( z ) \underset { n \to \infty } { \sim } \frac { ( - 1 ) ^ { n } A _ { n } } { \hat { \Delta } _ { n } ( \hat { \Delta } _ { n } - d ) \nu _ { n } \mu ( z ) z ^ { 1 + \tilde { \alpha } } } \sim \frac { \sqrt { 2 } } { \Gamma ( \tilde { \alpha } ) ^ { 2 } } \, n ^ { 2 \tilde { \alpha } - \frac { 3 } { 2 } } \, \frac { 1 } { \mu ( z ) z ^ { 1 + \tilde { \alpha } } } \.$$

Indeed, while ˆ p n satisfies the homogeneous Casimir equation (B.18), we know ˆ h n solves an inhomogeneous Casimir equation quoted in (A.8), which in the present normalisation for ˆ h n gives:

$$\left [ C _ { z } - \lambda _ { \hat { \Delta } _ { n } } \right ] \hat { h } _ { n } ( z ) = \nu _ { n } ^ { - 1 } \frac { A _ { n } } { \mu ( z ) z ^ { 1 + \tilde { \alpha } } }$$

with

$$\nu _ { n } ^ { - 1 } A _ { n } \underset { n \to \infty } { \sim } \frac { ( - 1 ) ^ { n } 4 \sqrt { 2 } } { \Gamma ( \tilde { \alpha } ) ^ { 2 } } \, n ^ { 2 \tilde { \alpha } + \frac { 1 } { 2 } } \.$$

At large n and fixed z , one can easily check that d n as given in (B.24) correctly reproduces the RHS of (B.25), up to corrections that become important when z ∼ n 2 . Such corrections must of course be present, since the exact ˆ h n ( z ) decay as z -˜ α for large z , cf. the second asymptotic expression A.18

To estimate the matrix elements ⟨ ˆ h n | ˆ h m ⟩ for -1 2 &lt; α &lt; ˜ 1 2 , we will first assume ˜ α &gt; 0 and then analytically continue to ˜ α &lt; 0. Using relation (A.8) and integrating by parts, the

off-diagonal matrix elements are given by

̸

$$\langle \hat { h } _ { n } | \hat { h } _ { m } \rangle = \frac { \nu _ { n } ^ { - 1 } A _ { n } \int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + d } } \hat { h } _ { m } ( z ) - \nu _ { m } ^ { - 1 } A _ { m } \int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + d } } \hat { h } _ { n } ( z ) } { ( \hat { \Delta } _ { m } - \hat { \Delta } _ { n } ) ( \hat { \Delta } _ { m } + \hat { \Delta } _ { n } - d ) } \, \quad n \ne m \.$$

For ˜ α &gt; 0 there are no boundary terms (unlike (B.19) above). The integrals in the numerator can actually be done exactly:

$$\int _ { 1 } ^ { \infty } \frac { \text{d} z } { z ^ { 1 + \tilde { \alpha } } } \hat { h } _ { n } ( z ) & = - \, \frac { A _ { n } } { 4 \nu _ { n } } \, \Gamma ( 2 - \frac { d } { 2 } ) \, \Gamma ( \tilde { \alpha } ) ^ { 2 } \Gamma ( 2 \tilde { \alpha } ) \\ & \quad \times \, _ { 4 } \tilde { F } _ { 3 } \left ( 1, 1 - n, 2 \tilde { \alpha }, 2 \tilde { \alpha } + n + 1 - \frac { d } { 2 } ; \tilde { \alpha } + 1, \tilde { \alpha } + 1, 2 \tilde { \alpha } + 2 - \frac { d } { 2 } ; 1 \right ) \.$$

Numerically, we find that for large n this scales as n 1 2 -2˜ α . We confirm the same scaling by a different argument. First note that:

$$\int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + \tilde { \alpha } } } \, \hat { h } _ { n } ( z ) = ( - 1 ) ^ { n + 1 } \int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + \tilde { \alpha } } } \, d _ { n } ( z )$$

since we claim that

$$\int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + \tilde { \alpha } } } \, \hat { p } _ { n } ( z ) = 0 \quad \text{ for all } n \geq 1 \.$$

This can be seen by (i) noting that this integral converges for ˜ α &gt; 0, (ii) deforming the contour to wrap the left cut ( -∞ , 0], and (iii) noting that z - -1 ˜ α g ˆ ∆ n ( z ) = g (˜) α ˆ ∆ n ( z ) has no left cut (see equation (2.35)). 24

Using the approximation (B.24) for d n ( z ) we have

$$\int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + \tilde { \alpha } } } \, d _ { n } ( z ) \sim \frac { \sqrt { 2 } } { \Gamma ( \tilde { \alpha } ) ^ { 2 } } \, n ^ { 2 \tilde { \alpha } - \frac { 3 } { 2 } } \int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 2 \tilde { \alpha } } } \, ( \frac { z - 1 } { z } ) ^ { 1 - d / 2 }$$

This is divergent in the large z region, but since the approximation breaks down for z ∼ n 2 we expect that the integral gets regulated at that scale, yielding

$$\int _ { 1 } ^ { \infty } \frac { d z } { z ^ { 1 + \alpha } } d _ { n } ( z ) \, \in \, n ^ { 2 \tilde { \alpha } - \frac { 3 } { 2 } } \times ( n ^ { 2 ( 1 - 2 \tilde { \alpha } ) } + \mathcal { O } ( 1 ) ) = n ^ { \frac { 1 } { 2 } - 2 \tilde { \alpha } } + \mathcal { O } ( n ^ { 2 \tilde { \alpha } - \frac { 3 } { 2 } } )$$

24 Eq. (B.30) may also be understood as the action of the first GFF functional with ˜ α → ˜ α -1 on ˆ . p n One is able to consider such functionals as long as the integral (B.30) converges. By duality, it must then vanish.

Overall we get

̸

$$\langle \hat { h } _ { n } | \hat { h } _ { m } \rangle _ { n, m \to \infty } = c \, ( - 1 ) ^ { n + m } \sqrt { n m } \, \frac { \left [ \left ( \frac { n } { m } \right ) ^ { 2 \tilde { \alpha } } - \left ( \frac { m } { n } \right ) ^ { 2 \tilde { \alpha } } \right ] } { m ^ { 2 } - n ^ { 2 } } \,, \quad n \neq m$$

where c is some numerical constant dependent on ˜. α

We must now concern ourselves with the case n = m . First let us write

$$\langle \hat { h } _ { n } | \hat { h } _ { m } \rangle = \langle d _ { n } | d _ { m } \rangle - ( - 1 ) ^ { n } \langle \hat { p } _ { m } | d _ { n } \rangle - ( - 1 ) ^ { n } \langle \hat { p } _ { n } | d _ { m } \rangle + \langle \hat { p } _ { n } | \hat { p } _ { m } \rangle$$

Combining previous results we find that ⟨ d n | d m ⟩ is suppressed for large n, m . Importantly this cannot change for n = m because of the monotonicity properties (B.23). 25 Therefore we have

$$\langle \hat { h } _ { n } | \hat { h } _ { n } \rangle ^ { d } = \langle d _ { n } | d _ { n } \rangle ^ { d } - 2 ( - 1 ) ^ { n } \langle \hat { p } _ { n } | d _ { n } \rangle ^ { d } + \langle \hat { h } _ { n } | \hat { h } _ { n } \rangle ^ { d } = 1 + O ( 1 / n )$$

We have checked that both (B.33) and (B.35) hold up to numerical scrutiny. Thus we have

̸

$$\langle \hat { h } _ { m } | \hat { h } _ { n } \rangle = \begin{cases} \ 1 + \mathcal { O } ( \frac { 1 } { n } ) &, \ \ n = m \\ c \sqrt { n m } \, \frac { \left [ \left ( \frac { n } { m } \right ) ^ { 2 \delta } - \left ( \frac { m } { n } \right ) ^ { 2 \delta } \right ] } { m ^ { 2 } - n ^ { 2 } } &, \ \ n \neq m \end{cases}$$

We now analytically continue this result to negative value ˜ α &lt; 0 (and remark that it passes numerical checks also in that case).

Let us apply Schur's test with the ansatz β n = α n = n δ . The sums ∑ m |⟨ ˆ h m | ˆ h n ⟩| m δ and ∑ n |⟨ ˆ h m | ˆ h n ⟩| n δ are finite only when

$$- 3 / 2 + 2 | \tilde { \alpha } | < \delta < \frac { 1 } { 2 } - 2 | \tilde { \alpha } | \.$$

Approximating ∑ m → n ∫ ∞ 1 /n dx where m = xn , we get

$$\sum _ { m } | \langle \hat { h } _ { m } | \hat { h } _ { n } \rangle | \, m ^ { \delta } \sim n ^ { \delta } \int _ { 1 / n } ^ { \infty } d x \frac { x ^ { 1 / 2 + \alpha } ( x ^ { 2 \tilde { \alpha } } - x ^ { - 2 \tilde { \alpha } } ) } { 1 - x ^ { 2 } } \sim n ^ { \delta }$$

25 This is also consistent with the estimate

$$\langle d _ { n } | d _ { m } \rangle \, \in \, n ^ { 2 \tilde { \alpha } - \frac { 3 } { 2 } } m ^ { 2 \tilde { \alpha } - \frac { 3 } { 2 } } \int _ { 1 } ^ { n m } \frac { d z } { z ^ { 2 \tilde { \alpha } } } ( \frac { z - 1 } { z } ) ^ { 1 - d / 2 } \, \in \, \frac { 1 } { \sqrt { n m } } ( 1 + \mathcal { O } ( n ^ { 2 \tilde { \alpha } - 1 } m ^ { 2 \tilde { \alpha } - 1 } ) ) \.$$

and Schur's test is passed (by symmetry between m,n ):

$$\sum _ { m } | \langle \hat { h } _ { m } | \hat { h } _ { n } \rangle | \, m ^ { \delta } \lesssim n ^ { \delta } \, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \sum _ { n } | \langle \hat { h } _ { m } | \hat { h } _ { n } \rangle | \, n ^ { \delta } \lesssim m ^ { \delta } \.$$

An appropriate value of δ satisfying (B.37) can be chosen if and only if | ˜ α &lt; | 1 2 .

We conclude that both T, S are bounded and by our previous argument T is invertible, thus establishing that { ˆ p n } , { ˆ h n } are biorthogonal Riesz bases for L 2 .

## B.2 Interacting Riesz bases

Now let us turn to general sparse spectra ∆ n = 2˜ + 2 α n + γ n and prove that p n ∝ I z g ∆ n form a Riesz basis for L 2 as long as | ˜ α &lt; | 1 2 .

We recall that the GFF spectrum and corresponding basis are denoted ∆ ˆ n = 2˜ + 2 α n and ˆ . It suffices to demonstrate that the map, p n

↦

$$M \colon \ \hat { p } _ { n } \mapsto p _ { n }$$

↦

defines a bounded, invertible operator from L 2 → L 2 , as it may then be composed with T of the previous subsection to obtain a new bounded invertible MT that maps e n → p n , thus establishing p n is Riesz. Below we will study the operator

↦

$$K = 1 - M \colon \ \hat { p } _ { n } \mapsto \hat { p } _ { n } - p _ { n } \.$$

We will begin by proving that K is a compact operator and hence bounded, implying that M is also bounded. Then we will prove that M is invertible.

Compactness of K . We claim that the operator K is compact (meaning it is a bounded operator and K U ( ) is compact for any bounded set U ⊂ L 2 ). The anomalous dimensions have been assumed to decay as a power for large n ,

$$| \gamma _ { n } | \stackrel { n \to \infty } { \lesssim } n ^ { - \varepsilon } \,,$$

but they can be large for finite values of n . For compactness, it is sufficient to show that K is a limit of finite-rank operators in the operator-norm. Defining the finite-rank operators

$$K _ { N } ( p _ { n } ^ { 0 } ) = \Theta ( n \leq N ) \left ( \hat { p } _ { n } - p _ { n } \right ) \,,$$

(with Θ( x ) the Heaviside function for the predicate x ), we will show that ∥ K N -K ∥ → 0 as N →∞ .

We will apply Schur's test to the matrix elements of ( K N -K ) ( † K N -K ),

$$\Delta K _ { m n } ^ { 2 } \coloneqq \Theta ( m, n > N ) \left \langle p _ { m } - \hat { p } _ { m } | p _ { n } - \hat { p } _ { n } \right \rangle \,.$$

Since γ n is assumed to decay then p n -ˆ p n ∼ ∂ p n ˆ n + . . . . Note that here it is understood

$$\partial _ { n } \hat { p } _ { n } \equiv \partial _ { n } \left [ \nu _ { n } ^ { - 1 } \mathcal { I } _ { z } g _ { \Delta } | _ { \Delta = 2 \tilde { \alpha } + 2 n } \right ] \,.$$

Hence

$$\Delta K _ { m n } ^ { 2 } \sim \Theta ( m, n > N ) \, \gamma _ { n } \gamma _ { m } \partial _ { m } \partial _ { n } \langle \hat { p } _ { m } | \hat { p } _ { n } \rangle \.$$

Differentiating the matrix elements (B.20) (which hold even for m,n non-integer) and using the bound (B.42), we find

$$| \Delta K _ { m n } ^ { 2 } | \lesssim \Theta ( m, n > N ) \left ( m n \right ) ^ { - \varepsilon } \langle \hat { p } _ { m } | \hat { p } _ { n } \rangle$$

Schur's test for sequences m , n δ δ gives

$$\sum _ { m > N } | \Delta K _ { n m } ^ { 2 } | m ^ { \delta } \sim n ^ { \delta - 2 \varepsilon } \int _ { N / n } ^ { \infty } \mathrm d x \frac { x ^ { \frac { 1 } { 2 } - \varepsilon + \delta } } { 1 - x ^ { 2 } } \lesssim n ^ { \delta } N ^ { - 2 \varepsilon }$$

for e.g. δ = . This implies that ( ε K N -K ) ( † K N -K ) is bounded with norm decaying like ( N -2 ε ) , and hence

$$\| K _ { N } - K \| = \sqrt { \| ( K _ { N } - K ) ^ { \dagger } ( K _ { N } - K ) \| } \lesssim N ^ { - \varepsilon } \rightarrow 0$$

as required. Hence K is compact and therefore bounded, and trivially M = 1 -K is also bounded.

Invertibility of M . To prove that M is invertible we will argue by contradiction. Suppose that M is not invertible, so that it is either not injective or not surjective. Since K is compact then M is what is known as a Fredholm operator , which are injective if and only if surjective. Hence M must be neither injective and nor surjective.

In particular, since M is not injective, there is a non-zero element ∑ n c n ˆ p n such that

$$M \left ( \sum _ { n } c _ { n } \hat { p } _ { n } \right ) = \sum _ { n } c _ { n } p _ { n } = 0 \.$$

Let us define :

̸

$$\tilde { g } _ { \Delta } \coloneqq \sum _ { n } \frac { c _ { n } } { \lambda _ { \Delta _ { n } } - \lambda _ { \Delta } } p _ { n } \quad \text{ for } \Delta \neq \Delta _ { n } \quad \left ( \lambda _ { \Delta } = \Delta ( \Delta - d ) \right ) \,.$$

This sum clearly converges in L 2 since the coefficients are more suppressed than c n . It also converges pointwise on (1 , ∞ ) since

$$\sum _ { n } | c _ { n } | ^ { 2 } < \infty \, \Longrightarrow \, | c _ { n } | \lesssim n ^ { - 1 / 2 } \,, \quad \ \lambda _ { n } \sim n ^ { 2 } \,,$$

and | p n ( z ) | &lt; const for fixed z (by the estimate (A.16)). Differentiating term-by-term, we see that (B.51) satisfies the same Casimir equation as I z g ∆ ( z ):

$$\left ( C _ { 2 } - \lambda _ { \Delta } \right ) \tilde { g } _ { \Delta } = \sum _ { n } c _ { n } \, p _ { n } = 0 \.$$

Moreover, both I z g ∆ ( z ) and ˜ g ∆ are invariant under ∆ → d -∆. There is only one such eigenfunction (up to normalisation), so we must have

$$\tilde { g } _ { \Delta } ( z ) = h _ { \Delta } \, \mathcal { I } _ { z } \, g _ { \Delta } ( z ) \.$$

The proportionality coefficient h ∆ may be determined by evaluating ˜ g ∆ ( z ) log z ∣ ∣ ∣ z = ∞ :

$$h _ { \Delta } = \frac { 1 } { \pi } \, \Gamma ( 1 - \frac { \Delta } { 2 } ) \, \Gamma ( 1 - \frac { d - \Delta } { 2 } ) \, \sum _ { n } \frac { c _ { n } \, \sqrt { \Delta _ { n } - \frac { d } { 2 } } \, \sin \left ( \frac { \pi } { 2 } \Delta _ { n } \right ) } { \lambda _ { \Delta _ { n } } - \lambda _ { \Delta } } \.$$

̸

Since the sum converges locally uniformly over ∆ = ∆ n (the summand behaving as n -2 ), the function h ∆ is meromorphic in ∆, with a countable set of zeros and poles. Hence for almost all ∆, we have

$$\mathcal { I } _ { z } \, g _ { \Delta } = h _ { \Delta } ^ { - 1 } \, \tilde { g } _ { \Delta } = h _ { \Delta } ^ { - 1 } \, M \left ( \sum _ { n } \frac { c _ { n } } { \lambda _ { \Delta _ { n } } - \lambda _ { \Delta } } \hat { p } _ { n } \right ) \in \text{Im} ( M )$$

Therefore the closure contains them all: I z g ∆ ∈ Im( M ) for all ∆. In particular Im( M ) contains the Riesz basis ˆ p n ∝ I z g ˆ ∆ n , so it is the whole space: Im( M ) = L 2 . Since M is Fredholm then its image is closed, and so Im( M ) = L 2 . Thus M is surjective, but since it is Fredholm, this implies it is also injective, giving the desired contradiction.

## C Analyticity and bounds on interacting functionals

The interacting functional kernels h n were defined in Section 3.2.2 as L 2 functions. In this appendix, we shall prove that, under our assumptions on the sparse spectrum, they are actually h n = i (˜) α ( f n ) for some f n ∈ A . In other words they are real-analytic and sufficiently polynomially bounded near 1 and ∞ .

Our strategy will be to use the following expansion:

$$h _ { m } = \sum _ { n } \theta _ { m } ( \hat { \Delta } _ { n } ) \, \hat { h } _ { n } \,, \quad \quad \theta _ { m } ( \Delta ) = \langle h _ { m } \, | \, \mathcal { I } _ { z } g _ { \Delta } \rangle \.$$

This is just the canonical L 2 -convergent decomposition of the L 2 elements h m with respect to the Riesz basis of free functional kernels ˆ h n . In this sum the ˆ h n are analytic - in fact they are essentially polynomials. The basic problem is that this representation of h m is too 'coarse': in particular the sum above cannot even be commuted with differentiation. We must therefore produce with a different representation that will make manifest the analyticity of h m . The key will be to replace the sum over integer n with a complex integral.

We will proceed as follows:

- · In Section C.1, we will show that the coefficients θ m ( ˆ ∆ ) admit an analytic extension n with respect to n that is polynomially bounded in a certain complex domain. This is non-trivial as θ m (∆) actually grows exponentially for complex ∆.
- · In Section C.2, we will use this fact to rewrite the sum in (C.1) as a integral over complex n ` la Sommerfeld-Watson. After deforming the contour, this representation a establishes analyticity of h m in a finite domain containing z ∈ (1 , ∞ ).
- · Finally in Sections C.3 and C.4, we will prove that this representation also establishes that the kernels h m are sufficiently bounded at z = ∞ and z = 1 so that indeed h m = i (˜) α ( f n ).

## C.1 Functional actions and analyticity

Here we will show that the coefficients θ m ( ˆ ∆ ) in (C.1) admit an analytic extension in n n of the form

$$\theta _ { m } ( \hat { \Delta } _ { n } ) = a ( n ) + ( - 1 ) ^ { n } c ( n ) \quad \text{ for } \ n > \Lambda$$

with a n , c n ( ) ( ) analytic and polynomially bounded in a 'wedge'-shaped region extending to infinity around the positive-real axis (see Figure 3),

$$W _ { \Lambda, \delta } = \left \{ n \in \mathbb { C } \ \Big | \ | \arg \left ( n - \Lambda \right ) | < \delta \right \} \,.$$

As we will see in Section C.2, this will be sufficient for establishing analyticity of h m . We will need to use our assumption (3.2) that the sparse spectrum { ∆ n } admits precisely this type of analytic extension on such a wedge. Equivalently, for the Casimir eigenvalues λ n := λ ∆ n , by assumption there are two analytic functions λ (+) ( n , λ ) ( -) ( n ) on the same region such that

$$\lambda _ { n } = \left ( \frac { 1 + ( - 1 ) ^ { n } } { 2 } \right ) \lambda ^ { ( + ) } ( n ) + \left ( \frac { 1 - ( - 1 ) ^ { n } } { 2 } \right ) \lambda ^ { ( - ) } ( n ) \,, \quad n > \Lambda \.$$

We note that this is definitely true for the GFF spectrum,

$$\hat { \lambda } _ { n } = \hat { \lambda } ( n ) \,, \quad \hat { \lambda } ( n ) \coloneqq \hat { \Delta } _ { n } ( \hat { \Delta } _ { n } - d ) \.$$

To prove (C.2), one needs to use the analyticity of the interacting spectrum ∆ n to learn about the analyticity of the interacting functional actions. The explicit formula (4.20) is the perfect tool to do this: it expresses the interacting functionals actions in terms of the dual interacting spectrum:

̸

$$\theta _ { m } ( \Delta ) = \prod _ { \substack { n = 1 \\ n \neq m } } ^ { \infty } \left ( \frac { \lambda _ { \Delta } - \lambda _ { n } } { \lambda _ { m } - \lambda _ { n } } \right ) \,, \quad \lambda _ { \Delta } = \Delta ( \Delta - d ) \,, \quad \lambda _ { n } \equiv \lambda _ { \Delta _ { n } }$$

Note that validity of this formula does not require that the kernels h m be analytic.

Furthermore, we will need to use the assumption (3.2) that the anomalous dimensions ∆ n -ˆ ∆ decay like a power on this wedge. We will write this as n

$$\lambda ^ { ( \pm ) } ( n ) & \underset { | n | \to \infty } { \sim } \hat { \lambda } ( n ) + c ^ { \pm } \hat { \lambda } ( n ) ^ { \frac { 1 - \eta _ { \pm } } { 2 } } + \dots \\ & = 4 n ^ { 2 } + 4 ( 2 \tilde { \alpha } - \frac { d } { 2 } ) n + c ^ { \pm } ( 2 n ) ^ { 1 - \eta _ { \pm } } + \dots$$

for constants c ± and η ± &gt; 0. These statements imply that there exists N such that for integer n ,

$$\hat { \lambda } _ { n } < \lambda ^ { ( + ) } ( n ) < \hat { \lambda } _ { n + 1 } \quad \forall n \geq N \quad \text{ or } \quad \hat { \lambda } _ { n - 1 } < \lambda ^ { ( + ) } ( n ) < \hat { \lambda } _ { n } \quad \forall n \geq N \,, \quad \text{ (C.8)}$$

and similarly for λ ( -) ( n ). That is, the eigenvalues are eventually always above or always below the GFF ones. This property will help us control the overall sign of the functional action.

Let us therefore set n &gt; N &gt; max( p, Λ) and write the functional action in the following form:

$$\theta _ { p } ( \hat { \Delta } _ { n } ) = \Pi _ { p \text{-ly} } ( n ) \, \Pi ^ { ( + ) } ( n ) \, \Pi ^ { ( - ) } ( n )$$

with

̸

$$\Pi _ { p \circ l y } ( n ) = \prod _ { \substack { m = 1 \\ m \neq p } } ^ { N - 1 } \left ( \frac { \hat { \lambda } _ { n } - \lambda _ { m } } { \lambda _ { p } - \lambda _ { m } } \right ) \,, \quad \Pi ^ { ( \pm ) } ( n ) = \prod _ { \substack { m = N \\ m \text{ even/odd} } } ^ { \infty } \left ( \frac { \hat { \lambda } _ { n } - \lambda ^ { ( \pm ) } ( m ) } { \lambda _ { p } - \lambda ^ { ( \pm ) } ( m ) } \right ) \,. \quad \text{(C.10)}$$

Since Π poly ( n ) is polynomial in n we may focus on analysing the infinite products. The analysis is similar for ( ± ), so for definiteness we consider here just the (+) case. We write

$$\Pi ^ { ( + ) } ( n ) = S ^ { ( + ) } ( n ) \, \exp \left \{ \sum _ { \substack { m = N \\ m \text{ even} } } ^ { \infty } \log \left | \frac { \hat { \lambda } _ { n } - \lambda ^ { ( + ) } ( m ) } { \lambda _ { p } - \lambda ^ { ( + ) } ( m ) } \right | \right \} \, \quad \, \quad \, \text{(C.11)}$$

with S (+) ( n ) an overall sign to account for the absolute values, as we will discuss below. The argument of the exponential can be written as:

$$- \frac { 1 } { 2 } \int _ { N - \delta } ^ { \infty } \mathrm d m \, \text{Im} \left [ \cot ( \frac { \pi } { 2 } m ) \right ] \, \text{Re} \left [ \log \left ( \frac { \hat { \lambda } _ { n } - \lambda ^ { ( + ) } ( m ) } { \lambda _ { p } - \lambda ^ { ( + ) } ( m ) } \right ) \right ] \,.$$

This can be massaged into:

$$& - \frac { 1 } { 2 } \int _ { N - \delta } ^ { \infty } \, \text{dm} \, \text{Im} \left [ \cot ( \frac { \pi } { 2 } m ) \, \log \left ( \frac { \hat { \lambda } _ { n } - \lambda ^ { ( + ) } ( m ) } { \lambda _ { p } - \lambda ^ { ( + ) } ( m ) } \right ) \right ] + \frac { \pi } { 2 } \int _ { N - \delta } ^ { \hat { n } } \text{dm} \, \text{Re} \left [ \cot ( \frac { \pi } { 2 } m ) \right ] \\ & = - \frac { 1 } { 2 } \oint _ { < } \frac { \text{dm} } { 2 i } \, \cot ( \frac { \pi } { 2 } m ) \, \log \left ( \frac { \hat { \lambda } _ { n } - \lambda ^ { ( + ) } ( m ) } { \lambda _ { p } - \lambda ^ { ( + ) } ( m ) } \right ) + \log \left ( \frac { \sin [ \frac { \pi } { 2 } \hat { n } ] } { \sin [ \frac { \pi } { 2 } ( N - \delta ) ] } \right ) \.$$

Here we define ˆ n (+) implicitly by the equation

$$\lambda ^ { ( + ) } ( \hat { n } ^ { ( + ) } ) = \hat$$

In fact, in deriving (C.13), we had to assume that ∂ m λ (+) &gt; 0 throughout the integration range. Thanks to (C.7) this may be assured by taking N large enough. It also implies that n → λ n ( ) is invertible, so (C.14) gives a valid definition of ˆ. n As for the symbol &lt; , it means that the integration contour has been deformed to a wedge in the complex m plane with apex at N -δ . That this can be done follows from the assumed analyticity properties of λ (+) m .

Performing a similar analysis for Π ( -) we find:

$$\theta _ { p } ( \hat { \Delta } _ { n } ) = ( - 1 )$$

with

$$E ^ { ( + ) } ( n ) & = - \frac { 1 } { 2 } \int _$$

Note that we have written ˆ λ n → ˆ ( λ n ). This representation makes clear that both E ( ± ) ( n ) are now also analytically defined in a certain domain in complex n space containing the real line. This domain contains a wedge of slope δ , i.e. the same slope appearing in (C.3) since for large n we have ˆ ( λ n ) ∼ λ ( ± ) ( n ).

We will now show that both these factors grow at most logarithmically with ˆ ( λ n ) (and hence with n ). To see this we first split the integration region between large | m | and | m | ∼ N , with the split performed at a scale M such that λ ± ( m ) is well approximated by its asymptotic form (C.7). Clearly in the small | m | region both integrals are bounded by the logarithm of ˆ ( λ n ). In the large | m | region we have analytic control over the asymptotics of λ ± ( m ) and we find that region contributes

$$\cdot \, \cdot \, \cdot \, & \quad \quad \quad$$

Thus we conclude that

$$P ( n ) \coloneqq \frac { \Pi _ { _ { \text{poly} } }$$

is analytic and polynomially bounded in a complex wedge in n space.

Finally, we need to analyse the various prefactors. Define

$$s _ { \pm } = \text{sign} \, c _ { \pm } \,.$$

We then have

$$( - 1 ) ^ { N + 1 } S ^ { ( + ) } ( n ) S ^ { ( - ) } ( n ) = \begin{cases} \begin{array} { c } - s _ { + } \,, & \text{for $n$ even} \\ + s _ { - } \,, & \text{for $n$ odd} \end{array} \end{cases}$$

As for the remaining factors, note that setting δn ˆ ( ± ) := ˆ n ( ± ) -n , we have

$$\sin \left [ \frac { \pi } { 2 } \hat { n } ^ { ( + ) } \right ] \cos \left [ \frac { \pi } { 2 } \hat { n } ^ { ( - ) } \right ] = \begin{cases} \begin{array} { c c } \sin \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( + ) } \right ] \cos \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( - ) } \right ] & \text{for $n$ even} \end{array}, \\ \cos \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( + ) } \right ] \sin \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( - ) } \right ] & \text{for $n$ odd} \end{cases},$$

and for large n ,

$$\delta \hat { n } ^ { \pm } & \sim - \frac { c _ { \pm } } { n ^ { 1 + \eta _ { \pm } } }$$

which implies that these prefactors are also analytic and polynomially bounded on a wedge in complex n space, after splitting between even and odd. Hence we conclude:

$$\theta _ { p } ( \hat { \Delta } _ { n } ) = \begin{cases} \ - 2 s _ { + } \sin \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( + ) } \right ] \cos \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( - ) } \right ] P ( n ) & \text{for $n$ even} \, \\ + 2 s _ { - } \sin \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( - ) } \right ] \cos \left [ \frac { \pi } { 2 } \delta \hat { n } ^ { ( + ) } \right ] P ( n ) & \text{for $n$ odd} \, \end{cases}$$

thus establishing the desired result.

## C.2 Analyticity of functional kernels

Each of the interacting functionals can be expanded as an L 2 -convergent sum of the GFF functionals ˆ h n (see eq. (3.28))

$$h ( z ) = \sum _ { n } b _ { n } \, \hat { h } _ { n } ( z )$$

We showed in Section C.1 that our assumptions on the sparse spectrum ∆ n imply that these coefficients satisfy

$$b _ { n } = a ( n ) + ( - 1 ) ^ { n } c ( n ) \, \quad n \in \mathbb { Z }$$

where a n ( ) and c n ( ) are polynomially bounded, analytic functions on some 'wedge' domain (C.3) near infinity. In this section we will show that this means that that h z ( ) is then realanalytic on (1 , ∞ ). The strategy will be to convert the infinite sum into a contour integral using a Sommerfeld-Watson type trick.

Firstly, let us set

$$\tilde { b } _ { n } \coloneqq b _ { n } \times ( \hat { \Delta } _ { n } - \frac { d } { 2 } ) \,,$$

Figure 7: The analytic structure of the integral representation (C.30) for interacting functional kernels in the auxiliary x -plane. K x,z ( ) is analytic in x away from the red branch cuts (whose positions depend on z , shown here for z = 3 2 ). g x ( ) is analytic outside the 'twoleaf' domain in blue (shown here for δ = π 7 ). The eventual integration contour used in (C.40) is shown here in black.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-ac274980-ea27-475b-b4c3-f1a5f8709bd0.jpg)

and introduce a generating function for the ˆ h n :

$$K ( x, z ) & \coloneqq \sum _ { n = 1 } ^ { \infty } x ^ { n - 1 } \, \frac { \hat { h } _ { n } ( z ) } { \hat { \Delta } _ { n } - \frac { d } { 2 } } \. \\ & = z ^ { - \tilde { \alpha } } \left ( \frac { z } { z - 1 } \right ) ^ { \frac { d } { 2 } - 1 } ( 1 + x ) ^ { - 2 - 2 \tilde { \alpha } + \frac { d } { 2 } } \, \Gamma ( 2 + 2 \tilde { \alpha } - \frac { d } { 2 } ) \\ & \quad \times _ { 3 } \tilde { F } _ { 2 } ( 1, 1 + \tilde { \alpha } - \frac { d } { 4 }, \frac { 3 } { 2 } + \tilde { \alpha } - \frac { d } { 4 } ; 1 + \tilde { \alpha }, 1 + \tilde { \alpha } ; \frac { 4 x } { ( 1 + x ) ^ { 2 } z } ) \.$$

For any z ∈ (1 , ∞ ), this is an analytic function of x with branch points at

$$x = - 1 \, \quad x = t + i \sqrt { 1 - t ^ { 2 } } \, \quad x = t - i \sqrt { 1 - t ^ { 2 } } \quad ( t \equiv \frac { 2 - z } { z } ) \. \quad \text{(C.28)}$$

The three branch cuts may be chosen to all meet at 1, running along the negative real axis ( -∞ -, 1) and along the unit circle from 1 to t + i √ 1 -t 2 and to t -i √ 1 -t 2 (in red in Figure 7). If any z ∈ (1 , ∞ ) gains a small enough imaginary part, the branch points t ± i √ 1 -t 2 move slightly off the unit circle and K x,z ( ) remains an analytic function of x away from the cuts.

We may now write the expansion (C.24) as

$$h ( z ) = \sum _ { n = 1 } ^ { \infty } \tilde { b } _ { n } \, \oint _ { 0 } \frac { d x } { 2 \pi i } \, \frac { 1 } { x ^ { n } } \, K ( x, z ) = \sum _ { n = 1 } ^ { \infty } \tilde { b } _ { n } \, \oint _ { \Gamma } \frac { d x } { 2 \pi i } \, \frac { 1 } { x ^ { n } } \, K ( x, z ) \.$$

In the first expression the contour tightly encloses x = 0, while in the second it is deformed towards the edge of the unit disk. We claim that we can exchange the order of summation and integration in (C.29) to obtain

$$h ( z ) = \oint _ { \Gamma } \frac { \mathrm d x } { 2 \pi i } \, g ( x ) \, K ( x, z ) \, \quad \ g ( x ) \coloneqq \sum _ { n = 1 } ^ { \infty } \tilde { b } _ { n } \, x ^ { - n } \.$$

First, we shall establish the domain of analyticity of g . Clearly this series converges absolutely to an analytic function for | x | &gt; 1, since ˜ b n is polynomially bounded. We will show that it admits an analytic extension to a strictly larger region: the non-analyticities are confined to a smaller 'two-leaf' region (in blue in Figure 7), which only touches the unit disc at ± 1 . It is sufficient to treat the tail of the sum, for which the analytic expression (C.25) for ˜ b n applies:

$$\tilde { g } ( x ) \coloneqq g ( x ) - \sum _ { n \leq \Lambda } \tilde { b } _ { n } \, x ^ { - n } = \sum _ { n > \Lambda } \left ( a ( n ) + ( - 1 ) ^ { n } c ( n ) \right ) x ^ { - n } \.$$

Following the trick of Sommerfeld and Watson, we may recast the sum over n as an integral:

$$\tilde { g } ( x ) = \oint _ { \Gamma _ { ( \Lambda, \infty ) } } \frac { \mathrm d n } { 2 \pi i } \, \frac { a ( n ) ( - x ) ^ { - n } + c ( n ) \, x ^ { - n } } { \sin \pi n } \.$$

To begin with, the integration contour wraps the segment n &gt; Λ of the real axis, but we will now deform it away from the real axis. To check that the integral will converge on such a contour, let us estimate the integrand for large complex n with Re( n ) &gt; 0: focusing first on the second term in (C.32), we note that

$$\left | \frac { x ^ { - n } } { \sin \pi n } \right | ^ { \, | n | \to \infty } \sim \exp \left ( - \, \pi \left | \text{Im} ( n ) \right | - \arg ( x ) \text{Im} ( n ) - \log | x | \text{Re} ( n ) \right ) \quad \text{(C.33)}$$

For | x | &gt; 1, the final term in the exponent is negative, and the first two terms combine to be non-positive. The same analysis applies (with x → -x ) to the other term. Hence the integrand of (C.32) is exponentially suppressed for n large with Re( n ) &gt; 0. One may then deform the contour to the wedge-shaped one (denoted &lt; ) inside the domain of analyticity W Λ ,δ for a n ( ), c n ( ) (see Figure 8)

$$\tilde { g } ( x ) = \oint _ { < } \frac { \mathrm d n } { 2 \pi i } \, \frac { a ( n ) ( - x ) ^ { n } + c ( n ) \, x ^ { - n } } { \sin \pi n } \.$$

Figure 8: The integration contour denoted &lt; in the n -plane in the expression (C.34). The domain of analyticity W Λ ,δ of the integrand is shown in blue, with its poles shown in red.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-531d5fed-e627-4f8f-9906-44768a5ec4b5.jpg)

This expression is now manifestly analytic for a larger region in x : using the same estimate (C.33) with arg n ∼ δ , the second term of the integrand is exponentially suppressed outside of a leaf-shaped region,

$$| x | < e ^ { - \tan \delta ( \pi - | \arg x | ) } \, \quad \ \ | \arg ( x ) | \leq \pi \.$$

Replacing x →-x , we see that the first term of the integrand is exponentially suppressed outside a mirrored leaf-shaped region, and the whole integrand is exponentially suppressed outside of a two-leaf region (the blue region in Figure 7),

$$R _ { \text{two-leaf} } \coloneqq \left \{ | x | < e ^ { - \tan \delta | \arg x | } \right \} \cup \left \{ | x | < e ^ { - \tan \delta ( \pi - | \arg x | ) } \right \} \,.$$

(Note that, if the angle δ of the wedge on which a n , c n ( ) ( ) are analytic tends to zero, the two-leaf region becomes the unit disc.) Hence it follows that the integral (C.34) converges locally uniformly over x ∈ C \ R two -leaf , so ˜( g x ), and hence g x ( ), is analytic on this region.

Now let us return to the question of interchanging the sum and integral in (C.29) to obtain (C.30). The integration contour Γ is pinched at x = -1 between the branch cut singularities of K x,z ( ) and the singularities of g in R two -leaf , as shown in Figure 7. We will prove at the same time that the integral in (C.30) converges, and that the sum and integral in (C.29) swap to give (C.30). First, these statements are clearly true away from x = -1 since the integrand is analytic and the sum converges locally uniformly.

Near the problematic point x = -1, we shall argue by dominated convergence. One can

check that

$$K ( x, z ) \underset { x \rightarrow - 1 } { \sim } \underset { \text{analytic} } { \sim } O ( 1 ) + \underbrace { O \left ( | x + 1 | ^ { - 2 \tilde { \alpha } + d / 2 } \right ) } _ { \text{non-analytic} } \.$$

Meanwhile, using the estimate (5.21) to get the leading behaviour

$$\tilde { b } _ { n } \sim - \left ( \hat { \Delta } _ { n } - \frac { d } { 2 } \right ) \gamma _ { n } \, \hat { \theta } ^ { \prime } ( \hat { \Delta } _ { n } ) = O ( n ^ { d / 2 - 2 - 2 \tilde { \alpha } - \epsilon } )$$

one finds as x approches the unit circle that

$$g ( x ) = \sum _ { n = 1 } ^ { \infty } \tilde { b } _ { n } \, x ^ { - n } \underset { | x | \to 1 } { \sim } \, \underset { \text{analytic} } { \underbrace { O ( 1 ) } { \, } + } \, \underbrace { O \left [ ( | x | - 1 ) ^ { 1 - d / 2 + 2 \tilde { \alpha } + \epsilon } \right ] } _ { \text{non-analytic} } \,.$$

For terms in the integral (C.30) involving the analytic part of either K or g , the contour can be deformed respectively to the left or the right and we are done. The term involving only the non-analytic pieces of K and g , which really has to pass through x = -1, is dominated by the integrable function | x +1 . |

We conclude that h z ( ) may represented as

$$h ( z ) = \oint _ { \Gamma _ { R _ { \text{two-leaf} } } } \frac { \text{d} x } { 2 \pi i } \, g ( x ) \, K ( x, z ) \.$$

with the contour wrapping the region R two -leaf as in Figure 7.

Note now that any z ∈ (1 , ∞ ) has a complex neighbourhood over which the integrand is analytic in z and the integral over x converges uniformly. The latter is true because the cuts along the unit circle exit the point -1 with arguments ± π -1 2 arg z (for example, see Figure 9a). For any δ &gt; 0, the two-leaf region departs from x = -1 at an angle strictly bounded away from the unit circle. Hence z can always be given a finite imaginary part without the cuts and the two-leaf region overlapping, meaning the integration contour does not get pinched anywhere other than x = -1.

Hence h z ( ) is real-analytic on (1 , ∞ ). Note that the radius of convergence of h z ( ) does generically shrink to zero at the endpoints z = 1 and z = ∞ , with the branch points x = t ± √ 1 -t 2 colliding and pinching the integration contour at x = 1 and x = -1 respectively (see Figure 9). As a result, we typically expect the interacting functionals h z ( ) to have a branch cut running along ( -∞ , 1). In the next sections we will verify that h z ( ) satisfies the correct properties at z = 1 and z = ∞ , namely appropriate polynomial boundedness, so that it lies in the image of A inside L 2 .

## C.3 Polynomial bound at z = ∞

Using the integral representation (C.40) for an interacting kernel h = h n , we shall prove the polynomial bound

$$| h _ { n } ( z ) | \underset { z \rightarrow \infty } { \leq } | z | ^ { - \tilde { \alpha } } \,.$$

This is the limit when the branch points at x = t ± i √ 1 -t 2 pinch the contour at x = -1, as shown in Figure 9a. First, we note that the parts of the contour away from x = -1 indeed scale as K x,z ( ) ∼ z -˜ α in this limit. The same applies for the terms involving the analytic part of g x ( ) in (C.39) or the analytic part of K x,z ( ) in (C.37), for which the contour can be deformed away from x = -1.

Turning to the part of the contour near x = -1 involving only the non-analytic parts, we find the bound

$$\left | \int _ { x \sim - 1 } \mathrm d x \, g ( x ) \, K ( x, z ) \right | \supsetneq \int _ { x \sim - 1 } | z | ^ { - \tilde { \alpha } } \int _ { x \sim - 1 } | \mathrm d x | \, | x + 1 | ^ { - 1 + \epsilon } \, \left | _ { 3 } \tilde { F } _ { 2 } ( y / z ) \right | \,, \ \ y \equiv \frac { 4 x } { ( 1 + x ) ^ { 2 } } \,, \ \text{(C.4.2)}$$

where 3 ˜ F 2 function is the one appearing in (C.27). Changing the integration variable x → y , properly accounting for the Jacobian factor, and expanding at large y , one obtains the bound

$$\left | \int _ { x \sim - 1 } \mathrm d x \, g ( x ) \, K ( x, z ) \right | \supsetneq \left | z \right | ^ { - \tilde { \alpha } } \int _ { z } ^ { \infty } \mathrm d y \, | y | ^ { - 1 - \epsilon / 2 } \left | _ { 3 } \tilde { F } _ { 2 } ( y / z ) \right | \sim | z | ^ { - \tilde { \alpha } - \frac { \epsilon } { 2 } }$$

giving the desired bound (C.41).

## C.4 Polynomial bound at z = 1

Finally, let us study the behaviour at z = 1. In order for the interacting kernel h ∈ L 2 to belong to the subspace i (˜) α ( A ) the definition (2.17) requires a bound on its behaviour as z → 1:

$$( z - 1 ) ^ { d / 2 - 1 } \, h _ { n } ( z ) \underset { z \to 1 } { \sim } \, \text{(analytic)} + O ( | z - 1 | ^ { \tilde { \alpha } } ) \.$$

We will prove this using the integral representation (C.40). As z → 1, the cuts along the unit circle pinch the contour at x = 1 as shown in Figure 9b. The parts of the contour away from x = 1 will lead to purely analytic contributions to ( z -1) d/ 2 -1 h z ( ) (with the prefactor coming from (C.27)). The same applies for the analytic part of g x ( ) in (C.39), for which the integration contour can be deformed inwards, away from x = 1.

The non-analytic terms that we need to bound in (C.44) correspond to the part of the contour passing through x = 1 containing only the non-analytic part of g x ( ). This

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-3870de0f-c7e7-4c47-9f5b-6a8537018486.jpg)

- (a) Pinching at x = -1 when z →∞

Figure 9: Close-ups of the branch cuts (red) pinching the integration contour (black) against the 'two-leaf' region (blue) at x = -1 and x = 1 in the respective limits z →∞ and z → 1. Here, as above, we denote t ≡ 2 -z z .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00572v2-figure-55328b09-71e6-4d7d-84ee-35447c3a7d89.jpg)

contribution is bounded by

$$\left | ( z - 1 ) ^ { d / 2 - 1 } \int _ { x \sim 1 } d x \, g ( x ) \, K ( x, z ) \right | \underset { | x - 1 | \sim | z - 1 | ^ { 1 / 2 } } { \lesssim } \int _ { \substack { | x - 1 | \sim | z - 1 | ^ { 1 / 2 } } } \, | d x | \, | x - 1 | ^ { - d / 2 + 2 \tilde { a } + \epsilon } \, \left | _ { 3 } \tilde { F } _ { 2 } ( y / z ) \right | \,, \, \text{(C.45)}$$

where y and 3 ˜ ( F 2 y ) are defined as in (C.42). Again changing integration variable x → y and including the relevant Jacobian, we obtain the bound

$$\left | ( z - 1 ) ^ { d / 2 - 1 } \int _ { x \sim 1 } \mathrm d x \, g ( x ) \, K ( x, z ) \right | & \lesssim \int _ { \substack { z \rightarrow 1 \\ | y - 1 | \sim | z - 1 | } } | \mathrm d y | \, | y - 1 | ^ { - d / 4 + \tilde { \alpha } + \epsilon / 2 } \, ( 1 + | y - 1 | ^ { ( d - 3 ) / 2 } ) \\ & \lesssim | z - 1 | ^ { \tilde { \alpha } + ( 4 - d ) / 4 + \epsilon / 2 } + | z - 1 | ^ { \tilde { \alpha } + ( d - 2 ) / 4 + \epsilon / 2 } \, \quad \, ( C. 4 6 )$$

Since d ≤ 4 then the first term is bounded by | z - | 1 ˜ α as required. The same applies for the second term as long as d ≥ 2. For the remaining case d = 1, it is necessary to make a slightly stronger boundedness assumption on the sparse spectrum ∆ n = 2˜ + 2 α n + γ n :

$$\gamma _ { n } = \alpha ( n ) + ( - 1 ) ^ { n } \beta ( n ) \quad \text{ with } \quad \alpha ( n ) \sim n ^ { - \epsilon } \, \quad \beta ( n ) \sim n ^ { - \frac { 1 } { 2 } - \epsilon }$$

instead of β n ( ) ∼ n -ϵ in (3.2).

## D On general values of ˜ α and subtractions

In Section 3, we used an intermediate L 2 space to generate complete Schauder bases of analytic functionals f n ∈ A , dual to particular sets of blocks g ∆ n . We stated there that such a basis exists dual to any sparse set, ∆ n = 2˜ + 2 α n + γ n , under assumptions of largen analyticity and power-law decay of γ n . We gave an argument for this fact in the main text in the special case ˜ α ∈ -( 1 2 , 1 2 ), and this appendix will treat the case of general values of ˜. α

## D.1 Almost all values of ˜ α

First let us consider the case ˜ α / ∈ 1 2 + N , i.e. all values related to ( -1 2 , 1 2 ) by an integer shift (these are almost all values of ˜). α

In Section 3.3, we explained the extra subtleties that arise when ˜ α &gt; 1 2 . Riesz bases for L 2 can still be used to produce Schauder bases for A and HF , but the embedding of L 2 becomes more non-trivial, and the bases are subject to a finite subtraction procedure reviewed in Section 3.3. There, we derived a set of sufficient conditions (a)-(c) in L 2 language to produce a Schauder basis f n ∈ A . In this section, we will show that these conditions are satisfied by any sparse spectrum satisfying our assumptions.

## (a) The set of discontinuities of blocks { p n = I z g ∆ n } can be completed to a Riesz basis by adding finitely many blocks.

As explained in eq. (3.32), if ˜ = α β + M for β ∈ ( -1 2 , 1 2 ) and M ∈ N , then one can just add M arbitrary blocks to the sparse spectrum ∆ n = ˜ + 2 α n + γ n to obtain an extended spectrum ∆ β n ∼ 2 β +2 n + γ n . Then the analysis of Section 3.2 applies, so the corresponding blocks { p n = I z g ∆ β n } form a Riesz basis.

## (b) Finite combinations of M of the dual Riesz basis { h n } take the form i (˜) α ( f n ) for some f n ∈ A .

We already know from Appendix C that each h n is analytic on (1 , ∞ ) and satisfies the same bounds (2.17) as the GFF functionals at z = 1 and z = ∞ . At z = 1, this is immediately sufficient. At z = ∞ , we encounter the need for additional subtractions for ˜ α &gt; 1 2 - even for the GFF functionals - since h n behave as a z -β , which is M powers less suppressed than the desired functionals i (˜) α ( f n ) ∼ z -α .

Therefore one should take the combinations of M of the h n that are sufficiently bounded at z = ∞ . More precisely, the combinations of interacting functionals h n that are in the

image of i (˜) α are:

$$h _ { n } - \sum _ { k = 1 } ^ { M } c _ { n k } ^ { \tilde { \alpha } } \, h _ { k } \sum _ { \substack { \in i ^ { ( \tilde { \alpha } ) } ( \mathcal { A } ) \, \quad \ n > M \\ \ell } } \,, \quad \, n > M$$

̸

$$c _ { n k } ^ { \tilde { \alpha } } \coloneqq \left ( \, \prod _ { \substack { l \neq k \\ 0 \leq l \leq M } } \frac { \lambda _ { \Delta _ { n } ^ { \beta } } - \lambda _ { \Delta _ { l } ^ { \beta } } } { \lambda _ { \Delta _ { k } ^ { \beta } } - \lambda _ { \Delta _ { l } ^ { \beta } } } \right ) \frac { \theta _ { n } ( \Delta _ { 0 } ^ { \beta } ) } { \theta _ { k } ( \Delta _ { 0 } ^ { \beta } ) }$$

where θ n (∆) := ⟨ h n |I z g ∆ ⟩ . This expression appears to depend on another (arbitrary) variable ∆ β 0 , but in fact it is independent of all the arbitrary parameters ∆ β 0 , . . . , ∆ . This β M is because any dependence on these arbitrary parameters would lead to multiple subtracted functionals (D.1) dual to the same subtracted blocks g ˜ α ∆ n , which in turn form a Schauder basis of HF . Hence the difference between them would be orthogonal to all basis elements, so must vanish.

To see that the combination (D.1) is suppressed by z -M at ∞ relative to h n , let us consider its action against blocks. Using the identity (D.15) derived below, one finds

$$\langle h _ { n } - \sum _ { k = 1 } ^ { M } c _ { n k } ^ { \tilde { \alpha } } \, h _ { k } \, | \, \mathcal { I } _ { z } \, g _ { \Delta } \rangle = \left ( \prod _ { 1 \leq k \leq M } \frac { \lambda _ { \Delta _ { n } ^ { \beta } } - \lambda _ { \Delta _ { k } ^ { \beta } } } { \lambda _ { \Delta } - \lambda _ { \Delta _ { k } ^ { \beta } } } \right ) \, \theta _ { n } ( \Delta ) \, \quad n > M \, \quad \text{(D.3)}$$

i.e. the functional action defined by this combination is suppressed by ( λ ∆ ) -M ∼ ∆ -2 M relative to θ n (∆). The large ∆ behaviour of the functional action is tied to the large z behaviour of the kernel with the identification z ∼ ∆ , signaling the desired 2 z -M suppression.

## (c) Decomposition of Cauchy kernel.

As explained in Section 3.3, the loss of M functionals is related by duality to the appearance of M relations (3.36) between the subtracted blocks in the space of hyperfunctions. Let us explain the origin of these relations.

We know that the blocks p n = I z g ∆ β n form a Riesz basis for L 2 so, for example, any block can be uniquely decomposed into them:

$$\mathcal { I } _ { z } \, g _ { \Delta } = \sum _ { n = 1 } ^ { \infty } \theta _ { n } ( \Delta ) \, \mathcal { I } _ { z } \, g _ { \Delta _ { n } ^ { \beta } } \, \in L ^ { 2 } \.$$

Since the embedding j (˜) α : L 2 → HF is continuous, then we obtain a convergent decomposition in HF :

$$g _ { \Delta } ^ { ( \tilde { \alpha } ) } = \sum _ { n = 1 } ^ { \infty } \theta _ { n } ( \Delta ) \, g _ { \Delta _ { n } ^ { \beta } } ^ { ( \tilde { \alpha } ) } \, \in \, \mathcal { H F } \,.$$

However, we claim that, when ˜ α &gt; 1 2 , the blocks on the RHS are not independent as elements of HF . Rather, there are relations among them - morally because they have been 'oversubtracted' by ˜, which is α M units greater than β . (In other words, for ˜ α &gt; 1 2 , the map j (˜) α does not have a continuous inverse.)

One can eliminate M blocks by applying the Casimir operator to the decomposition (D.5). We define the 'subtracted Casimir', acting on hyperfunctions, by

$$C _ { z } ^ { ( \tilde { a } ) } H \coloneqq z ^ { - \tilde { \alpha } - 1 } C _ { z } ( z ^ { \tilde { \alpha } + 1 } H ) \, \quad C _ { z } \coloneqq 4 ( 1 - z ) ^ { 1 - d / 2 } z ^ { 1 + d / 2 } \partial _ { z } [ ( 1 - z ) ^ { d / 2 } z ^ { 1 - d / 2 } \partial _ { z } ]$$

defined equivalently by either acting on the analytic function H away from the cuts, or on its discontinuity I z H itself. Note that in general the result C (˜) α z H may not be a member of the same space HF since its behaviour at ∞ may be worse, but this is not a concern for individual blocks since they are eigenfunctions of the subtracted Casimir:

$$\left ( C _ { z } ^ { \left ( \tilde { \alpha } \right ) } - \lambda _ { \Delta } \right ) g _ { \Delta } ^ { \left ( \tilde { \alpha } \right ) } = 0 \.$$

̸

One now needs to analyse whether the action of the Casimir preserves the convergence of the sum (D.5) in HF . In Section D.2.1 below, we show that it can be applied precisely M times. Hence, denoting ∆ ≡ ∆ in (D.5) and applying the combination 0 ∏ l = k 0 ≤ ≤ l M ( C (˜) α z -λ ∆ β l ) for 0 &lt; k ≤ M , we obtain the desired relations expressing the lowest M blocks in terms of the higher ones:

$$g ^ { ( \tilde { \alpha } ) } _ { \Delta ^ { \beta } _ { k } } = - \sum _ { n > M } c ^ { \tilde { \alpha } } _ { n k } \, g ^ { ( \tilde { \alpha } ) } _ { \Delta ^ { \beta } _ { n } } \, \quad k = 1, \dots, M$$

where c ˜ α nk are the same coefficients defined in (D.2).

As a result, we may eliminate all of the 'additional' blocks and, as argued in Section 3.3, one obtains the following decomposition of the Cauchy kernel into the physical blocks ∆ n :

$$\frac { 1 } { z - w } = \sum _ { n = 1 } ^ { \infty } g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, f _ { n } ( z )$$

converging as a hyperfunction in w . We further need to prove that this expression converges as an analytic function of z . This essentially follows from the fact the dimensions ∆ n decay to the GFF ones at large n . One may argue similarly to (3.30),(3.31) in the unsubtracted case that this sum converges in the same way as in the GFF case:

$$g _ { \Delta _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, f _ { n } ( z ) \sim g _ { \hat { \Delta } _ { n } } ^ { ( \tilde { \alpha } ) } ( w ) \, \hat { f } _ { n } ( z ) \, \left ( 1 + O ( n ^ { - \epsilon } ) \right ) \,.$$

The subtracted GFF functionals ˆ f n coincide with the analytic continuation of the GFF functionals (A.7) to ˜ α &gt; 1 2 (see discussion in Appendix D.3 below). Hence they satisfy the estimate (2.33), so the sum (D.10) converges absolutely, as required, for z in a neighbourhood

of the cut [1 , ∞ ] and w away from the cut.

## D.2 An identity for interacting functionals

## D.2.1 How many times can the Casimir be applied?

We will derive an identity for interacting functional actions by applying the Casimir to the infinite sum (D.5). First we need understand whether it will preserve the sum's convergence, and how many times it can be applied.

To begin with, we note that the sum (D.4) in L 2 converges with the bounds O (log z ) and [ O (( z -1) 1 -d/ 2 ) + O (1)] at the respective endpoints. Hence it follows that the sum (D.5) in HF converges with behaviour at least as good as O z ( - - -˜ α 1 ϵ ) and [ O (( z -1) 1 -d/ 2 ) + O (1)] (in the sense discussed in footnote 8).

Applying the subtracted Casimir (D.6) to hyperfunctions worsens their z → ∞ and z → 1 behaviour by one power, so applying it P times will give a sum converging with behaviour at least as good as O z ( P - - -˜ α 1 ϵ ) and [ O (( z -1) -P +1 -d/ 2 ) + O (( z -1) -P )]. In order to have convergence in the space HF , we require that these exponents satisfy:

$$P - \tilde { \alpha } - 1 - \epsilon < 0 \quad \text{ and } \quad \max \left ( P, P - 1 + \frac { d } { 2 } \right ) < \tilde { \alpha } + 1 \. \quad \text{ (D.11)}$$

In the worse case d = 3, this requirement becomes

$$P < \tilde { \alpha } + \frac { 1 } { 2 } \.$$

By definition, the largest integer less than ˜ + α 1 2 is M . Hence the Casimir can be applied to (D.5) at most M times.

## D.2.2 Deriving the identity

Let us derive an identity for the interacting functionals, which shows that the subtracted functionals (D.1) have an action (D.3) on blocks that is more suppressed for large ∆.

(˜) α

̸

On one hand, setting ∆ ≡ ∆ and applying the combination 0 ∏ l = k 0 ≤ ≤ l M ( C z -λ ∆ β l ) to (D.5) leads to the relation (3.36) between Riesz basis elements. Substituting this back into (D.5) gives the Schauder basis decomposition in HF ,

̸

$$g _ { \Delta } ^ { ( \tilde { \alpha } ) } = \sum _ { n > M } \left [ \theta _ { n } ( \Delta ) - \sum _ { 0 < k \leq M } \left [ \, \prod _ { 0 \leq l \leq M } \frac { \lambda _ { \Delta _ { n } ^ { \beta } } - \lambda _ { \Delta _ { l } ^ { \beta } } } { \lambda _ { \Delta _ { k } ^ { \beta } } - \lambda _ { \Delta _ { l } ^ { \beta } } } \, \right ] \frac { \theta _ { n } ( \Delta _ { 0 } ^ { \beta } ) } { \theta _ { k } ( \Delta _ { 0 } ^ { \beta } ) } \, \theta _ { k } ( \Delta ) \right ] g _ { \Delta _ { n } ^ { \beta } } ^ { ( \tilde { \alpha } ) }$$

On the other hand, applying the combination ∏ 1 ≤ ≤ l M ( C (˜) α z -λ ∆ β l ) and setting ∆ = ∆ k

gives another decomposition of g (˜) α ∆ into the same basis,

$$g _ { \Delta } ^ { ( \tilde { \alpha } ) } = \sum _ { n > M } \left [ \, \prod _ { 1 \leq k \leq M } \frac { \lambda _ { \Delta _ { n } ^ { \beta } } - \lambda _ { \Delta _ { k } ^ { \beta } } } { \lambda _ { \Delta } - \lambda _ { \Delta _ { k } ^ { \beta } } } \, \right ] \theta _ { n } ( \Delta ) \, g _ { \Delta _ { n } ^ { \beta } } ^ { ( \tilde { \alpha } ) }$$

Equating the two decompositions (D.13) and (D.14) and noting, e.g. by applying the functionals, that their coefficients must coincide, we obtain the non-trivial identity (for n &gt; M ):

̸

$$\theta _ { n } ( \Delta ) - \sum _ { 0 < k \leq M } \left [ \, \prod _ { \substack { l \neq k \\ 0 \leq l \leq M } } \frac { \lambda _ { \Delta _ { n } ^ { \beta } } - \lambda _ { \Delta _ { l } ^ { \beta } } } { \lambda _ { \Delta _ { k } ^ { \beta } } - \lambda _ { \Delta _ { l } ^ { \beta } } } \, \right ] \frac { \theta _ { n } ( \Delta _ { 0 } ^ { \beta } ) } { \theta _ { k } ( \Delta _ { 0 } ^ { \beta } ) } \, \theta _ { k } ( \Delta ) = \left [ \, \prod _ { 0 < k \leq M } \frac { \lambda _ { \Delta _ { n } ^ { \beta } } - \lambda _ { \Delta _ { k } ^ { \beta } } } { \lambda _ { \Delta } - \lambda _ { \Delta _ { k } ^ { \beta } } } \, \right ] \theta _ { n } ( \Delta ) \, \cdot \, \text{ (DD.15)}$$

One can check that the explicit product formula (4.20) for θ n (∆) indeed solves this identity: this provides a non-trivial cross-check of that formula and the results of Section 4.

## D.3 All values of ˜ α

For all values ˜ α / ∈ 1 2 + Z , we have constructed complete Schauder bases of functionals { f ˜ α n } for the spaces A ˜ α , dual to the spectra ∆ ˜ α n := 2˜ +2 α n + γ n . (In this section, we shall indicate explicitly the dependence on ˜.) α

Let us now define functionals for the final case ˜ α ∈ 1 2 + Z by simply taking a pointwise limit of the functional kernels we have already constructed:

$$f _ { n } ^ { 1 / 2 + M } \coloneqq \lim _ { \tilde { \alpha } \to ( 1 / 2 + M ) ^ { + } } f _ { n } ^ { \tilde { \alpha } } \, \quad M \in \mathbb { Z }$$

In the rest of this section, we will argue that (D.16) leads to well-defined functionals f 1 2+ / M n that are: (i) elements of A 1 2+ / M , (ii) complete in A 1 2+ / M and (iii) dual to the spectrum ∆ 1 2+ / M n .

For clarity of presentation, we will set M = 0 and treat here the ˜ = α 1 2 case, but the others work in the same way.

Elements of A 1 2 / . First, it is easy to argue that f 1 2 / n ( z ) satisfy the appropriate bounds at the endpoints z = ∞ and z = 1 by starting from the polynomial bounds (C.41) and (C.44) derived above for ˜ α / ∈ 1 2 + Z . The subleading terms depend on ˜ in a way that is α uniform as ˜ α → 1 + 2 : in other words, none of the subleading contributions blow up and spoil the estimates, so they still apply in this limit. To show that f 1 2 / n ∈ A 1 2 / , all that remains is then to show real-analyticity.

To do this, let us take α, β / ∈ 1 2 + Z with α &gt; β . Since A ⊂ A α β , and { f β n } is a complete

basis for A β , then f α n ∈ A α may be decomposed with respect to this basis:

$$f _ { n } ^ { \alpha } = \sum _ { m } ( f _ { n } ^ { \alpha }, g _ { \Delta _ { m } ^ { \beta } } ^ { ( \beta ) } ) \, f _ { m } ^ { \beta } \.$$

The coefficients ( f α n , g ( β ) ∆ β m ) appearing here are not simple in terms of our functional actions θ ˜ α n (∆) ≡ ( f α n , g ( α ) ∆ ). Instead, let us consider decomposing the elements z β -α f α n :

$$z ^ { \beta - \alpha } f _ { n } ^ { \alpha } = \sum _ { m } ( z ^ { \beta - \alpha } f _ { n } ^ { \alpha }, g _ { \Delta _ { m } ^ { \beta } } ^ { ( \beta ) } ) \, f _ { m } ^ { \beta } = \sum _ { m } \theta _ { n } ^ { \alpha } ( \Delta _ { m } ^ { \beta } ) \, f _ { m } ^ { \beta } \.$$

We claim that this sum converges in A β (with the weak topology), uniformly over α in a neighbourhood of 1 2 . To show this, it is sufficient to note that the largem behaviour of the expansion coefficients θ α n (∆ β m ) is uniform over ˜. α

Hence, we may take the limit α → 1 + 2 in (D.18) with β &lt; 1 2 to obtain

$$z ^ { \beta - 1 / 2 } \, f _ { n } ^ { 1 / 2 } = \sum _ { m } \theta _ { n } ^ { 1 / 2 } ( \Delta _ { m } ^ { \beta } ) \, f _ { m } ^ { \beta } \quad \text{in} \, \mathcal { A } ^ { \beta } \.$$

Since the space A β is complete, then z β -1 2 / f 1 2 / n is an element and is hence real-analytic. But this means f 1 2 / n is also real-analytic on (1 , ∞ ). We have therefore shown that the pointwise limit (D.16) converges and is an element of A 1 2 / .

Completeness. On the other hand, to show that { f 1 2 / m } are complete in A 1 2 / , we would instead like to take the limit β → 1 + 2 in (D.18). For this, we note that (D.18) can be derived from the Cauchy kernel decomposition (3.38), which we repeat here:

$$\frac { 1 } { z - w } = \sum _ { m } g _ { \Delta _ { m } } ^ { ( \beta ) } ( w ) \, f _ { m } ^ { \beta } ( z ) \.$$

This converges as an analytic function of z and a hyperfunction of w . Moreover, this convergence, controlled by the estimate (2.33), is completely uniform over β near 1 2 . In particular, it converges in the weak topology for the space A β with respect to the z variable, uniformly over β . Acting on this hyperfunction of w with the functional z β -α f α n gives the decomposition (D.18), and tells us that its weak-topology convergence is uniform over β .

We may therefore take the limit β → 1 + 2 to obtain

$$z ^ { 1 / 2 - \alpha } f _ { n } ^ { \alpha } = \sum _ { m } \theta _ { n } ^ { \alpha } ( \Delta _ { m } ^ { 1 / 2 } ) \, f _ { m } ^ { 1 / 2 } \quad \text{ in } \mathcal { A } ^ { 1 / 2 } \.$$

Since { f α n } are complete in A α , the objects on the LHS are complete in a space z 1 2 / -α A α with slightly softer behaviour at z = ∞ . Thus we learn that the space generated by { f 1 2 / m }

includes at least that one. Taking α arbitrarily close to 1 2 , that space becomes A 1 2 / and we learn that { f 1 2 / m } are complete.

Duality. Now that we have established that { f 1 2 / n } are a complete set in A , their duality to the desired spectrum is immediate. For example, starting with α &gt; 1 2 , we have

$$\delta _ { m n } = ( f _ { m } ^ { \tilde { \alpha } }, \, g _ { \Delta _ { n } ^ { \tilde { \alpha } } } ^ { ( \tilde { \alpha } ) } ) \equiv \oint _ { \Gamma } \frac { \text{d} z } { 2 \pi i } \, f _ { m } ^ { \tilde { \alpha } } ( z ) \, g _ { \Delta _ { n } ^ { \tilde { \alpha } } } ^ { ( \tilde { \alpha } ) } ( z ) \.$$

This integral clearly converges uniformly over ˜, so taking the limit α ˜ α → 1 + 2 through the integral gives

$$( f _ { m } ^ { 1 / 2 }, \, g _ { \Delta _ { n } ^ { 1 / 2 } } ^ { ( 1 / 2 ) } ) = \delta _ { m n } \.$$

Continuity across thresholds. Finally, we note that this whole arugment could have been made instead with the lower limit ˜ α → ( 1 2 + M ) -in (D.16). The resulting functionals would be identical, since they are dual to the same spectrum. 26 Therefore, while the relation to the L 2 space is subject to a subtraction procedure across the values ˜ α ∈ 1 2 + Z , everything physical (i.e. involving the physical spaces A and HF ) is perfectly continuous across these thresholds - with the upper and lower limits coinciding.

## E Perturbative computation of interacting bases for special values of ∆ ϕ

This short appendix contains technical details for the analytic computation of interacting functionals in the particular example (5.25) with ∆ ϕ = 1. We begin by setting

$$\delta h _ { n } ( z ) \coloneqq h _ { n } ( z ) - \hat { h } _ { n } ( z ) = \frac { f _ { n } ( z ) - \hat { f } _ { n } ( z ) } { \pi \mu ( z ) z ^ { 1 + \Delta _ { \phi } } } = - \sum _ { m = 1 } ^ { \infty } \gamma _ { m } \, \hat { \theta } _ { n } ^ { \prime } ( \hat { \Delta } _ { m } ) \, \hat { h } _ { m } ( z ) \,,$$

and recall that γ n = 2 γ 0 ˆ ∆ ( ˆ n ∆ n -1) . Acting with the Casimir operator on δh n and using (A.8) the above relation leads to two new sums. One of them is the same as the original but with γ n set to a constant, and is thus proportional to ∂ h ˜ α n α | ˜=∆ ϕ . The second sum can be performed explicitly, as the z -dependence of the summands factors out. We thus get

$$C _ { z } \, \delta h _ { n } ( z ) = \frac { \gamma _ { 0 } } { 2 } \, \partial _ { \tilde { \alpha } } \hat { h } _ { n } | _ { \tilde { \alpha } = \Delta _ { \phi } } - \frac { B _ { n } } { \mu ( z ) \, z ^ { 1 + \Delta _ { \phi } } }$$

26 Explicitly, the difference between the upper and lower limits would be a functional that vanishes on all blocks g (1 / 2+ M ) ∆ 1 / 2+ M n . These blocks will be dense in HF 1 / 2+ M (by an argument similar to that for A above), forcing a functional vanishing at all of them to be the zero functional.

where

$$B _ { n } = \sum _ { m = 0 } ^ { \infty } \gamma _ { m } \, A _ { m } \, \theta _ { n } ^ { \prime } ( \hat { \Delta } _ { m } ) \.$$

We can now solve for δh n ( z ) by making an ansatz:

$$\delta h _ { n } ( z ) = a _ { 1 } \, \partial _ { \tilde { a } } \hat { h } _ { n } | _ { \tilde { a } = \Delta _ { \phi } } + a _ { 2 } \, \hat { h } _ { n } + a _ { 3 } \, g _ { \Delta = 0 } ( z ) + a _ { 4 } \, g _ { \Delta = 1 } ( z ) + a _ { 5 } \, m ( z ) + a _ { 6 } \, n ( z ) \.$$

Here m z ( ) and n z ( ) satisfy

$$C _ { z } \, m ( z ) = \frac { 1 } { \mu ( z ) z ^ { 1 + \Delta _ { \phi } } }, \quad C _ { z } \, n ( z ) = \frac { 1 } { \mu ( z ) z ^ { 1 + \Delta _ { \phi } } } \log ( z ) \.$$

and can be solved for explicitly for ∆ ϕ = 1. Imposing (E.2), as well as δh n being sufficiently well behaved for large z , uniquely determines all a i .

## References

- [1] I. Bena, On the construction of local fields in the bulk of AdS(5) and other spaces , Phys. Rev. D 62 (2000) 066007, [ hep-th/9905186 ].
- [2] A. Hamilton, D. N. Kabat, G. Lifschytz, and D. A. Lowe, Local bulk operators in AdS/CFT: A Boundary view of horizons and locality , Phys. Rev. D 73 (2006) 086003, [ hep-th/0506118 ].
- [3] A. Hamilton, D. N. Kabat, G. Lifschytz, and D. A. Lowe, Holographic representation of local bulk operators , Phys. Rev. D 74 (2006) 066009, [ hep-th/0606141 ].
- [4] A. Hamilton, D. N. Kabat, G. Lifschytz, and D. A. Lowe, Local bulk operators in AdS/CFT: A Holographic description of the black hole interior , Phys. Rev. D 75 (2007) 106001, [ hep-th/0612053 ]. [Erratum: Phys.Rev.D 75, 129902 (2007)].
- [5] A. Hamilton, D. N. Kabat, G. Lifschytz, and D. A. Lowe, Local bulk operators in AdS/CFT and the fate of the BTZ singularity , AMS/IP Stud. Adv. Math. 44 (2008) 85-100, [ arXiv:0710.4334 ].
- [6] D. Kabat, G. Lifschytz, and D. A. Lowe, Constructing local bulk observables in interacting AdS/CFT , Phys. Rev. D 83 (2011) 106009, [ arXiv:1102.2910 ].
- [7] D. Kabat, G. Lifschytz, S. Roy, and D. Sarkar, Holographic representation of bulk fields with spin in AdS/CFT , Phys. Rev. D 86 (2012) 026004, [ arXiv:1204.0126 ].
- [8] D. Kabat and G. Lifschytz, CFT representation of interacting bulk gauge fields in AdS , Phys. Rev. D 87 (2013), no. 8 086004, [ arXiv:1212.3788 ].

- [9] D. Kabat and G. Lifschytz, Decoding the hologram: Scalar fields interacting with gravity , Phys. Rev. D 89 (2014), no. 6 066010, [ arXiv:1311.3020 ].
- [10] D. Kabat and G. Lifschytz, Finite N and the failure of bulk locality: Black holes in AdS/CFT JHEP , 09 (2014) 077, [ arXiv:1405.6394 ].
- [11] D. Kabat and G. Lifschytz, Bulk equations of motion from CFT correlators , JHEP 09 (2015) 059, [ arXiv:1505.03755 ].
- [12] D. Kabat and G. Lifschytz, Locality, bulk equations of motion and the conformal bootstrap , JHEP 10 (2016) 091, [ arXiv:1603.06800 ].
- [13] V. F. Foit, D. Kabat, and G. Lifschytz, Bulk reconstruction for spinor fields in AdS/CFT JHEP , 02 (2020) 129, [ arXiv:1912.00952 ].
- [14] D. Kabat and G. Lifschytz, Dressing bulk fields in AdS 3 , JHEP 10 (2020) 189, [ arXiv:2008.01198 ].
- [15] N. Levine and M. F. Paulos, Bootstrapping bulk locality. Part I: Sum rules for AdS form factors , arXiv:2305.07078 .
- [16] S. El-Showk and M. F. Paulos, Bootstrapping conformal field theories with the extremal functional method , Phys. Rev. Lett. 111 (Dec, 2013) 241601, [ arXiv:1211.2810 ].
- [17] D. Mazac, Analytic bounds and emergence of AdS 2 physics from the conformal bootstrap , JHEP 04 (2017) 146, [ arXiv:1611.10060 ].
- [18] D. Maz´ aˇ, L. Rastelli, and X. Zhou, c An analytic approach to BCFT d , JHEP 12 (2019) 004, [ arXiv:1812.09314 ].
- [19] 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 ].
- [20] M. F. Paulos and B. Zan, A functional approach to the numerical conformal bootstrap , JHEP 09 (2020) 006, [ arXiv:1904.03193 ].
- [21] I. Heemskerk, J. Penedones, J. Polchinski, and J. Sully, Holography from Conformal Field Theory , JHEP 10 (2009) 079, [ arXiv:0907.0151 ].
- [22] T. Hartman, D. Maz´ aˇ, and L. Rastelli, c Sphere Packing and Quantum Gravity , JHEP 12 (2019) 048, [ arXiv:1905.01319 ].
- [23] P. Kravchuk, J. Qiao, and S. Rychkov, Distributions in CFT. Part I. Cross-ratio space , JHEP 05 (2020) 137, [ arXiv:2001.08778 ].
- [24] J. Qiao and S. Rychkov, Cut-touching linear functionals in the conformal bootstrap , JHEP 06 (2017) 076, [ arXiv:1705.01357 ].

| [25]   | O. Christensen et al. , An introduction to frames and Riesz bases , vol. 7. Springer, 2003.                                                                                                                                      |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [26]   | M. Meineri, J. Penedones, and T. Spirig, Renormalization group flows in AdS and the bootstrap program , arXiv:2305.11209 .                                                                                                       |
| [27]   | U. Graf, Introduction to hyperfunctions and their integral transforms: an applied and computational approach . Springer Science & Business Media, 2010.                                                                          |
| [28]   | W. F. Osgood, Note on the functions defined by infinite series whose terms are analytic functions of a complex variable; with corresponding theorems for definite integrals , The Annals of Mathematics 3 (1901), no. 1/4 25-34. |
| [29]   | S. G. Krantz, On limits of sequences of holomorphic functions , The Rocky Mountain Journal of Mathematics (2013) 273-283.                                                                                                        |
| [30]   | 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 ].                                                                                      |
| [31]   | F. A. Dolan and H. Osborn, Conformal partial waves and the operator product expansion , Nucl. Phys. B678 (2004) 491-507, [ hep-th/0309180 ].                                                                                     |
| [32]   | M. F. Paulos, Dispersion relations and exact bounds on CFT correlators , JHEP 08 (2021) 166, [ arXiv:2012.10454 ].                                                                                                               |
| [33]   | B. Girault, N. Suchel, K. Ghosh, and M. F. Paulos, CFT 1 /QFT 2 : RG flows, massless scattering and bulk reconstruction , To appear.                                                                                             |
| [34]   | B. S. Pavlov, Basicity of an exponential system and muckenhoupt's condition , in Doklady Akademii Nauk , vol. 247 (1), pp. 37-40, Russian Academy of Sciences, 1979.                                                             |
| [35]   | G. Semmler, Complete interpolating sequences, the discrete muckenhoupt condition, and conformal mapping , Ann. Acad. Sci. Fenn. Math 35 (2010) 23-46.                                                                            |
| [36]   | M. I. Kadets, The exact value of the paley-wiener constant , in Dokl. Akad. Nauk SSSR , vol. 155 (6), pp. 1253-1254, 1964.                                                                                                       |
| [37]   | L. L. Littlejohn and A. Zettl, The legendre equation and its self-adjoint operators , Electronic Journal of Differential Equations 2011 (2011), no. 69 1-33.                                                                     |
| [38]   | A. Zettl, Sturm-liouville theory . No. 121. American Mathematical Soc., 2005.                                                                                                                                                    |
| [39]   | B. Mukhametzhanov and A. Zhiboedov, Modular invariance, tauberian theorems and microcanonical entropy , JHEP 10 (2019) 261, [ arXiv:1904.06359 ].                                                                                |
| [40]   | S. Pal, Bound on asymptotics of magnitude of three point coefficients in 2D CFT , JHEP 01 (2020) 023, [ arXiv:1906.11223 ].                                                                                                      |
| [41]   | S. Ganguly and S. Pal, Bounds on the density of states and the spectral gap in CFT 2 Phys. Rev. D 101 (2020), no. 10 106022, [ arXiv:1905.12636 ].                                                                               |

- [42] S. Pal and Z. Sun, Tauberian-Cardy formula with spin , JHEP 01 (2020) 135, [ arXiv:1910.07727 ].
- [43] E. Bahiru, A. Belin, K. Papadodimas, G. Sarosi, and N. Vardian, State-dressed local operators in AdS/CFT , arXiv:2209.06845 .
- [44] S. Leutheusser and H. Liu, Subalgebra-subregion duality: emergence of space and time in holography , arXiv:2212.13266 .
Nat Levine Miguel F. Paulos
Bootstrapping bulk locality. Part II: Interacting functionals
Locality of bulk operators in AdS imposes stringent constraints on their
description in terms of the boundary CFT. These constraints are encoded as sum
rules on the bulk-to-boundary expansion coefficients. In this paper, we
construct families of sum rules that are (i) complete and (ii) `dual' to sparse
CFT spectra. The sum rules trivialise the reconstruction of bulk operators in
strongly interacting QFTs in AdS space and allow us to write down explicit,
exact, interacting solutions to the locality problem. Technically, we
characterise `completeness' of a set of sum rules by constructing Schauder
bases for a certain space of real-analytic functions. In turn, this allows us
to prove a Paley-Weiner type theorem characterising the space of sum rules.
Remarkably, with control over this space, it is possible to write down
closed-form `designer sum rules', dual to a chosen spectrum of CFT operators
meeting certain criteria. We discuss the consequences of our results for both
analytic and numerical bootstrap applications.