On the modular operator of mutli-component regions in chiral CFT

TL;DR

This paper presents a novel method to compute the modular flow for multi-component regions in chiral conformal field theory, utilizing locality, analyticity, and the KMS condition to transform the problem into a Riemann-Hilbert problem.

Contribution

It introduces a new approach based on locality and analyticity to determine the modular flow in multi-region chiral CFT, connecting it to Riemann-Hilbert problems.

Findings

The method successfully transforms the modular flow problem into a Riemann-Hilbert problem.

Examples demonstrate the applicability of the approach to specific cases.

The approach offers a new perspective on modular operators in chiral CFT.

Abstract

We introduce a new approach to find the Tomita-Takesaki modular flow for multi-component regions in general chiral conformal field theory. Our method is based on locality and analyticity of primary fields as well as the so-called Kubo-Martin-Schwinger (KMS) condition. These features can be used to transform the problem to a Riemann-Hilbert problem on a covering of the complex plane cut along the regions, which is equivalent to an integral equation for the matrix elements of the modular Hamiltonian. Examples are considered.

Tables

Table 1. Table 1: Periodicity properties, with M = e − i ​ π ​ τ ​ e − 2 ​ π ​ i ​ u 𝑀 superscript 𝑒 𝑖 𝜋 𝜏 superscript 𝑒 2 𝜋 𝑖 𝑢 M=e^{-i\pi\tau}e^{-2\pi iu} .

	${\vartheta }_1$	${\vartheta }_2$	${\vartheta }_3$	${\vartheta }_4$
${\vartheta }_i\left(u+1\right)/{\vartheta }_i\left(u\right)$	$-1$	$-1$	$1$	$1$
${\vartheta }_i\left(u+\tau \right)/{\vartheta }_i\left(u\right)$	$-M$	$M$	$M$	$-M$

Equations

$$\mathrm{Tr}(\varphi (x){\rho }_A^{it}\varphi (y){\rho }_A^{1-it})=⟨{\Omega }_0|\varphi (x){\Delta }^{it}\varphi (y){\Omega }_0⟩.$$

$$F(s,x,y)=\{\begin{array}{ll}⟨{\Omega }_0|\varphi (x)[1-e^s(1-\Delta {)}^{-1}{]}^{-1}\varphi (y){\Omega }_0⟩& \text{if ℑ(x)<0,}\\ ⟨{\Omega }_0|\varphi (y)[1-e^s(1-{\Delta }^{-1}{)}^{-1}{]}^{-1}\varphi (x){\Omega }_0⟩& \text{if ℑ(x)>0.}\end{array}$$

$$(1-e^s)F(s,x-i0,y)-F(s,x+i0,y)=⟨{\Omega }_0|[\varphi (x),\varphi (y)]{\Omega }_0⟩.$$

$$⟨{\Omega }_0|\varphi (x)(\ln \Delta )\varphi (y){\Omega }_0⟩={\int }_0^{\infty }\mathrm{d}sF(s,x,y).$$

$${\phi }_{X,Y}(t-i)={\phi }_{Y,X}(-t).$$

$$(1-e^{-\beta E_k}){\phi }_{\mathbf{p},\mathbf{k}}(0)=e^{-\beta E_k}{\delta }^3(\mathbf{k}-\mathbf{p}).$$

$$SX|\Omega ⟩=X^{\ast }|\Omega ⟩.$$

$$t\mapsto {\phi }_{X,Y}(t)=\omega (X{\sigma }^t(Y))\equiv ⟨\Omega |X{\Delta }^{it}Y\Omega ⟩$$

$${\phi }_{X,Y}(t-i)=\omega ({\sigma }^t(Y)X).$$

$${\rho }_{\omega }=\sum _{j=1}^n{p}_j|j⟩⟨j|,\omega (X)={\mathrm{Tr}}_{{\mathbb{C}}^n}(X{\rho }_{\omega })(X\in \mathfrak{M}).$$

$${\Delta }^{\frac{1}{2}}={\rho }_{\omega }^{\frac{1}{2}}\otimes {\rho }_{\omega }^{-\frac{1}{2}}.$$

$${\alpha }^g(X)={\mathrm{Ad}}_{U(g)}(X)\equiv U(g)XU(g{)}^{\ast }.$$

$$A=\bigcup _{j=1}^p(a_j,b_j)\subset \mathbb{R}\text{or S below,}$$

$$\mathfrak{A}(A)=\underset{i=1}{\overset{p}{\bigvee }}\mathfrak{A}(I_i),$$

$${\Delta }^{it}=U(g_t),g_t(z)=\frac{a(b-z)e^{-2\pi t}+b(z-a)}{(b-z)e^{-2\pi t}+(z-a)},$$

$$[L_n,L_m]=(n-m)L_{n+m}+\frac{\kappa }{12}n(n^2-1){\delta }_{n,-m},[L_n,\kappa ]=0.$$

$$\Vert (1+L_0{)}^k{L}_n\Psi \Vert \le \sqrt{c/2}(|n|+1{)}^{k+3/2}\Vert (1+L_0{)}^{k+1}\Psi \Vert$$

$$T(z)=-\frac{1}{2\pi }\sum _{n=-\infty }^{\infty }{L}_n{z}^{-n-2}.$$

$$T(f)={\int }_{\text{SS}}T(z)f(z)\mathrm{d}z:=-\frac{1}{2\pi }\sum _{n=-\infty }^{\infty }\left({\int }_{\text{SS}}{z}^{-n-2}f(z)\mathrm{d}z\right)L_n$$

$$\Gamma f(z)=-z^2\overline{f(z)},$$

$${\mathfrak{A}}_{\mathrm{Vir}}(I)=\{e^{iT(f)}\mid f\in C^{\infty }(I),\Gamma f=f{\}}^{\prime \prime }$$

$$\varphi (z)=\frac{1}{\sqrt{2\pi }}\sum _{n\in \mathbb{Z}}{\varphi }_n{z}^{-n-h}.$$

$$\Vert (1+L_0{)}^k\varphi (f)\Psi \Vert \le C\Vert (1+L_0{)}^{k+h-1}\Psi \Vert \sum _{n\in \mathbb{Z}}|{\hat{f}}_{-n-h}|(1+|n|{)}^{h-\frac{1}{2}}\le C_{\Psi }\underset{j\le h}{\sup }\Vert f^{(j)}{\Vert }_{L^{\infty }},$$

$$i[L_m,\varphi (z)]=z^{m+1}\frac{\mathrm{d}}{\mathrm{d}z}\varphi (z)+h(\rho )z^m\varphi (z)\text{or}U(g)\varphi (z)U(g{)}^{\ast }=[g^{\prime }(z){]}^h\varphi (g(z)),$$

$$⟨{\Omega }_0|\varphi (x)\varphi (y)|{\Omega }_0⟩=\frac{e^{-i\pi h}}{2\pi (x-y{)}^{2h}}.$$

$$j(z)=\frac{1}{\sqrt{2\pi }}\sum _{n\in \mathbb{Z}}{j}_n{z}^{-n-1}$$

$$W(f)=e^{ij(f)},\text{f∈CΓ∞(SS),}$$

$$W(f)W(g)=e^{iC(f,g)/2}W(f+g),W(f{)}^{\ast }=W(-f),C(f,g)=\frac{1}{2}{\int }_{\text{SS}}(g^{\prime }f-f^{\prime }g)\mathrm{d}z.$$

$${\mathfrak{A}}_{U(1)}(I)\equiv \{W(f)\mid f\in C_{\Gamma }^{\infty }(I){\}}^{\prime \prime }$$

$$\psi (z)=\frac{1}{\sqrt{2\pi }}\sum _{n\in \mathbb{Z}}{\psi }_{n-\frac{1}{2}}{z}^{-n}$$

Full text

1 Introduction

The reduced density matrix of a subsystem induces an intrinsic internal dynamics called the “modular flow”. The flow is non-trivial only for non-commuting observable algebras – i.e., in quantum theory – and depends on both the subsystem and the given state of the total system. It has been subject to much attention in theoretical physics in recent times because it is closely related to information theoretic concepts. As examples for some topics such as Bekenstein bounds, Quantum Focussing Conjecture, c-theorems, holography we mention [1, 2, 3, 4, 5]. In mathematics, the modular flow has played an important role in the study of operator algebras through the work of Connes, Takesaki and others, see [6] for an encyclopedic account.

It has been known almost from the beginning that the modular flow has a geometric nature in local quantum field theory when the subsystem is defined by a spacetime region of a simple shape such as an interval in chiral conformal field theory (CFT) [7, 8, 9]: it is the 1-parameter group of Möbius transformations leaving the interval fixed. For more complicated regions, important progress was made only much later in a pioneering work by Casini et al. [10], who were able to determine the flow for multi-component regions for a chiral half of free massless fermions in two dimensions. Recently in [11] they have generalized their method to the conformal theory of a chiral $U(1)$-current. Unfortunately, the method by [10, 11], as well as all other concrete methods known to the author, is based in an essential way on special properties of free quantum field theories. The purpose of this paper is to develop methods that could give a handle on the problem in general chiral CFTs, i.e. the left-moving half of a CFT on a (compactified) lightray in 1+1 dimensional Minkowski spacetime, and to make some of the constructions in the literature rigorous by our alternative method.

Consider a (possibly mixed) state in the chiral CFT, described by a density matrix $\rho$. Typical states of interest are the vacuum $\rho=|\Omega_{0}\rangle\langle\Omega_{0}|$, or a thermal state $\rho=e^{-\beta L_{0}}/\operatorname{Tr}e^{-\beta L_{0}}$. Given a union $A=\cup_{j}(a_{j},b_{j})$ of intervals of the (compactified) lightray, we can consider its reduced density matrix $\operatorname{Tr}_{A^{\prime}}\rho=\rho_{A}$, where $A^{\prime}$ is the complement of $A$. For the purposes of this discussion, we restrict to the vacuum state, although in the main part, thermal states will play a major role as well. If $\phi(x)$ is a primary field localized at $x\in A$, the “modular flow” is the Heisenberg time evolution $\rho_{A}^{it}\phi(x)\rho^{-it}_{A}$. The object $\rho_{A}$ is not actually well-defined in quantum field theory, but the modular flow is. Below we will use the rigorous framework of Tomita-Takesaki theory in our construction, but for pedagogical purposes, we here pretend that $\rho_{A}$ exists. Formally, the Hilbert space $\mathcal{H}$ splits as $\mathcal{H}_{A}\otimes\mathcal{H}_{A^{\prime}}$ and if $\rho$ is pure, then $\rho_{A}$ is formally a density matrix on $\mathcal{H}_{A}$. Its – equally formal – logarithm $H_{A}=\operatorname{ln}\rho_{A}$ is called the modular Hamiltonian in the physics literature.

In mathematical terms, the quantity which is well defined is the operator $\Delta=\rho_{A}\otimes\rho_{A^{\prime}}^{-1}$. For $x\in A$, we can then also write $\rho_{A}^{it}\phi(x)\rho^{-it}_{A}=\Delta^{it}\phi(x)\Delta^{-it}$ and $\operatorname{ln}\Delta=H_{A}\otimes 1_{A^{\prime}}-1_{A}\otimes H_{A^{\prime}}$. Furthermore, one can write

[TABLE]

and since the conformal primaries generate the full Hilbert space (mathematically, the Reeh-Schlieder theorem), we see that we knowledge of this quantity for all primaries $\phi$ suffices, in principle, to determine all matrix elements of $\Delta^{it}$, hence the operator itself, hence the flow. Alternatively, to know the generator of the flow, it suffices to know $\langle\Omega_{0}|\phi(x)(\operatorname{ln}\Delta)\phi(y)\Omega_{0}\rangle$. It is those types of quantities which we will study in this paper.

Our main trick is the following observation and it variants. For $s>0$ and fixed $y\in A$, define a function of $x$ on the complex plane cut along the intervals $A$,

[TABLE]

Then not only do the usual properties of CFTs imply that this function is holomorphic on the mutliply cut plane, but we also know its jumps across the cuts, given by the functional equation

[TABLE]

The commutator on the right side is given by a sum of $\delta$-functions and their derivatives by locality. We also prove certain further general properties of this function such as the degree of divergences as $x$ approaches $y$ or any boundary of a cut which depend on the conformal dimension of $\phi$. Using this functional equation and a standard contour argument appearing frequently in the study of Riemann-Hilbert type problems, we then obtain a linear integral equation for $F$ of Cauchy-type, see cor. 2 (in the case of bosonic fields). The desired matrix elements of the modular Hamiltonian are related by the integral

[TABLE]

A variant of this method also works for fermionic fields and for thermal states where the corresponding function $F$ lives on a torus cut along $A$ and satisfies a corresponding integral equation, see cor. 1.

The basis of our method is in some sense an old trick in quantum statistical mechanics. Consider a statistical operator $\rho$. The expectation functional acting an observable $X$ is $\omega(X)=\operatorname{Tr}(X\rho)$ and the modular flow acting on an observable $X$ is, by definition, $\sigma^{t}(X)=\rho^{it}X\rho^{-it}$. For observables $X,Y$, consider the function $\varphi_{X,Y}(t)=\omega(X\sigma^{t}(Y))=\operatorname{Tr}(X\rho^{it}Y\rho^{1-it})$. Since $\rho$ is a positive operator, one expects this function to be analytic inside the strip $\{t\in\mathbb{C}|-1<\Im(t)<0\}$. The values at the two boundaries of the strip are evidently related by the functional equation

[TABLE]

This functional equation is called the “KMS-condition.” Its fundamental importance was first understood in [12].

Note that on the right side, $X,Y$ appear in opposite order. Thus, if we have information about their commutator, we can sometimes get a closed equation for $\varphi_{X,Y}(t)$ or related quantities. This is for instance the case for the ideal quantum Bose or Fermi gas, where $\rho=Z^{-1}e^{-\beta H},H=\sum E_{k}a^{*}_{\bf k}a_{\bf k}$, the modular flow is just the Heisenberg time evolution with ‘time’ parameter $-\beta t$, and where one takes $Y=a_{\bf k},X=a_{\bf p}^{*}$. Using the commutators and $\varphi_{\bf p,k}(t)=e^{-it\beta E_{k}}\varphi_{\bf p,k}(0)$, the KMS condition is thereby equivalent in the case of bosons to the condition

[TABLE]

In this way, one can easily derive the standard Bose-Einstein (or Fermi-Dirac-) formula for the 2-point function. In our problem, these ideas are modified and applied to the reduced density matrix $\rho_{A}$ of a CFT.

This paper is organized as follows. In secs. 2, 3, we review basic notions from operator algebras, Tomita-Takesaki theory, and the operator algebraic approach to CFT (conformal nets) in order to make the paper self-contained. In secs. 4, 5 we introduce our method and study several examples. We conclude in sec. 6. Some conventions for elliptic functions are described in the appendix.

Notations and conventions: Gothic letters ${\mathfrak{A}},{\mathfrak{M}},\dots$ denote $*$-algebras, usually v. Neumann algebras. Calligraphic letters $\mathcal{H},\mathcal{K},\dots$ denote linear spaces, always assumed to be separable. The inverse temperature $\beta$ and modular parameter $\tau$ are related by $-2\pi i\tau=\beta$. The branches of $\operatorname{ln}z$ and $z^{\alpha}$ are taken along the negative real axis. ${\mathbb{S}}=\{z\in\mathbb{C}\mid|z|=1\}$ denotes the unit circle, ${\mathbb{D}}^{\pm}$ its interior/exterior.

Note added in proof: After this preprint was submitted, it was pointed out to us by the authors of [39] that one of our calculations related to thermal states contained an error, creating a tension between some of our results and those by [39], see also [38]. We are grateful to these authors for making us aware of this issue, which has been fixed in the current version.

2 Review of modular theory

2.1 Modular flow

For the convenience of the unfamiliar reader we review the basic elements of modular ($=$ Tomita-Takesaki-) theory; detailed references are [6, 13]. Connections to quantum information theory are described in [14]. An exposition directed towards a theoretical physics audience is [15].

The notion of modular flow is embedded into the theory of v. Neumann algebras. Such an algebra, ${\mathfrak{M}}$, can be defined as a complex linear space of bounded operators on some Hilbert space111We always assume that $\mathcal{H}$ is separable. $\mathcal{H}$ that is closed under taking products, adjoints (denoted by $*$). Such limits are understood in the so called “weak” topology, i.e. convergence of matrix elements. It is common to denote by ${\mathfrak{M}}^{\prime}$ the commutant, defined as the set of all bounded operators on $\mathcal{H}$ commuting with all operators in ${\mathfrak{M}}$.

To define the objects of main interest of the theory, one has to assume that ${\mathfrak{M}}$ is in “standard form”, meaing: $\mathcal{H}$ contains a “cyclic and separating” vector for ${\mathfrak{M}}$, that is, a unit vector $|\Omega\rangle$ such that the set consisting of $X|\Omega\rangle$, $X\in{\mathfrak{M}}$ is a dense subspace of $\mathcal{H}$, and such that $X|\Omega\rangle=0$ implies $X=0$ for any $X\in{\mathfrak{M}}$. The point is that one can then consistently define the anti-linear Tomita operator $S$ on the domain ${\mathcal{D}}(S)=\{X|\Omega\rangle\mid X\in{\mathfrak{M}}\}$ by the formula

[TABLE]

The cyclic property is needed in order that $S$ is densely defined, whereas without the separating property the definition would not be self-consistent. One can show that $S$ is a closable operator. This technical property guarantees that $S$ has a polar decomposition. It is customarily denoted by $S=J\Delta^{\frac{1}{2}}$, where $J$ anti-linear and unitary and $\Delta$ self-adjoint and non-negative. Tomita-Takesaki theory is about the interplay between the operators $\Delta,J$ and the algebras ${\mathfrak{M}},{\mathfrak{M}}^{\prime}$. The basic theorem is:

1. (i)

$J$ exchanges ${\mathfrak{M}}$ with the commutant in the sense that $J{\mathfrak{M}}J={\mathfrak{M}}^{\prime}$. Furthermore, $J^{2}=1,J\Delta J=\Delta^{-1}$. 2. (ii)

The modular flow $\sigma^{t}(X)=\Delta^{it}X\Delta^{-it}$ leaves ${\mathfrak{M}}$ and ${\mathfrak{M}}^{\prime}$ invariant for all $t\in\mathbb{R}$. 3. (iii)

From the vector $|\Omega\rangle$, one can define the state functional $\omega(X)=\langle\Omega|X\Omega\rangle,\quad\omega:{\mathfrak{M}}\to\mathbb{C}.$ It is positive and normalized (meaning $\omega(X^{*}X)\geq 0\,\,\forall X\in{\mathfrak{M}},\omega(1)=1$), and invariant under the modular flow in the sense that $\omega\circ\sigma^{t}=\omega$ for all $t\in\mathbb{R}$. The KMS-condition holds: for all $X,Y\in{\mathfrak{M}}$, the bounded function

[TABLE]

has an analytic continuation to the strip $\{z\in\mathbb{C}\mid-1<\Im z<0\}$ with the property that its boundary value for $\Im z\to-1^{+}$ exists and is equal to

[TABLE]

A partial converse to (iii) is: If $\omega^{\prime}$ is a normal (i.e. continuous in the weak∗-topology) positive linear functional on ${\mathfrak{M}}$, then it has a unique vector representative $|\Omega^{\prime}\rangle$ in the natural cone ${{\cal P}}^{\sharp}:=\{Xj(X)|\Omega\rangle\mid X\in{\mathfrak{M}}\}$, where $j(X)=JXJ$; in other words $\omega^{\prime}(X)=\langle\Omega^{\prime}|X\Omega^{\prime}\rangle$ for all $X\in{\mathfrak{M}}$.

The objects $J,\Delta,{{\cal P}}^{\sharp}$ depend on the algebra ${\mathfrak{M}}$ and the state $|\Omega\rangle$.

Example 0: Even tough Tomita-Takesaki theory is most interesting in the case of infinite dimensional v. Neumann algebras of types II, III, it helps with intuition to have in mind the finite dimensional case, i.e. the type In’ (algebra of $n$ by $n$ matrices). In this case, ${\mathfrak{M}}=M_{n}(\mathbb{C})\otimes 1_{n}$, which acts on the Hilbert space $\mathcal{H}=\mathbb{C}^{n}\otimes\mathbb{C}^{n}$. Evidently, the commutant is ${\mathfrak{M}}^{\prime}=1_{n}\otimes M_{n}(\mathbb{C})$. A vector $|\Omega\rangle$ in this Hilbert space is cyclic and separating if $|\Omega\rangle=\sum_{j=1}^{n}\sqrt{p_{j}}|j\rangle\otimes|j\rangle$ in some ON basis $\{|j\rangle\}$ and iff all $p_{j}>0$, $\sum_{j=1}^{n}p_{j}=1$. The state functional $\omega$ can be written in this example in terms of the “reduced density matrix”

[TABLE]

In fact, any positive normalized state functional $\omega^{\prime}$ arises from a unique reduced density matrix $\rho_{\omega^{\prime}}$ in this way. It is easy to go through the definition of $\Delta,J$ via $S$ giving for instance that

[TABLE]

Therefore, the modular flow is $\sigma^{t}(X)=\rho^{it}_{\omega}X\rho^{-it}_{\omega}$. The “modular Hamiltonian” is defined as the self-adjoint operator $\operatorname{ln}\Delta$. In our example, therefore, $\operatorname{ln}\Delta=\operatorname{ln}\rho_{\omega}\otimes 1_{n}-1_{n}\otimes\operatorname{ln}\rho_{\omega},$ where the first term belongs to ${\mathfrak{M}}$ and the second to ${\mathfrak{M}}^{\prime}$. It is important to stress that the split of $\operatorname{ln}\Delta$ into a part from ${\mathfrak{M}}$ and one from ${\mathfrak{M}}^{\prime}$ is impossible for general v. Neumann algebras, in particular for the type III1-factors appearing in quantum field theories222The possibility of making the split implies that $\sigma^{t}$ is inner, i.e. can be written as $\sigma^{t}(X)=U(t)XU(t)^{*}$ for unitaries $U(t)$ in ${\mathfrak{M}}$. One characterization of type $III$ v. Neumann algebras is that $\sigma^{t}$ precisely cannot be inner for any normal state $\omega$.. Therefore, apart from trivial cases, the object $\operatorname{ln}\rho_{\omega}$, hence the reduced density operator $\rho_{\omega}$ itself, does not exist. On the other and, $\operatorname{ln}\Delta$ and $\omega$ always exist. We will make sure to work with these well-defined objects in our setting.

Sometimes, a state $\omega$ is only given as an abstract (weakly continuous) expectation functional on an abstract333We mean a $C^{*}$-algebra with a preferred “folium” of normal states, see [6]. In particular, it is not assumed that ${\mathfrak{M}}$ is a priori represented by bounded linear operators on some Hilbert space. v. Neumann algebra ${\mathfrak{M}}$. Then one can perform the basic but very important GNS construction in order to obtain a Hilbert space in which the state is represented by a vector.

The starting point of this construction is the simple observation that the algebra ${\mathfrak{M}}$ itself, as a linear space, always forms a representation $\pi$ by left multiplication, i.e. $\pi(X)Y\equiv XY$. To equip this representation with a Hilbert space structure, it is natural to define $\langle X|Y\rangle=\omega(X^{*}Y)$, but this will in general lead to non-zero vectors with vanishing norm, unless $\omega$ is separating. Introduce ${\mathfrak{J}}_{\omega}=\{X\in{\mathfrak{M}}\mid\omega(X^{*}X)=0\}$. By the Cauchy-Schwarz inequality, $|\omega(X^{*}Y)|\leq\omega(X^{*}X)^{1/2}\omega(Y^{*}Y)^{1/2}$, we have ${\mathfrak{J}}_{\omega}=\{X\in{\mathfrak{M}}\mid\forall Y\in{\mathfrak{M}},\omega(Y^{*}X)=0\}$, so it is a closed linear subspace and a left ideal of ${\mathfrak{M}}$ containing precisely the null vectors. We can then define $\mathcal{H}_{\omega}={\mathfrak{M}}/{\mathfrak{J}}_{\omega}$ and complete it in the induced inner product. The left representation induces a representation on $\mathcal{H}_{\omega}$ which is called $\pi_{\omega}$. It is the desired GNS-representation. The vector $|\Omega_{\omega}\rangle\in\mathcal{H}_{\omega}$ representing $\omega$ is simply the equivalence class of the unit operator, $1$. It is by construction “cyclic” in the sense that the set $\pi_{\omega}({\mathfrak{M}})|\Omega_{\omega}\rangle$ is dense in $\mathcal{H}_{\omega}$. The vector is standard if $\omega(X^{*}X)$ implies $X=0$ (meaning ${\mathfrak{J}}_{\omega}=\{0\}$), in which case we say that it is faithful.

3 Review of chiral CFTs

3.1 Conformal nets on the real line (lightray)

One way to formalize the structure of chiral conformal quantum field theories (CFTs) is via nets of operator algebras. A chiral conformal field theory is associated with one lightray. It is given abstractly by an assignment of an algebra of operators ${\mathfrak{A}}(I)$ with each open interval $I=(a,b)\subset\mathbb{R}$ of the this lightray.

This assignment is called a conformal net if it obeys the following rules (see [18] for a general introdution to algebraic quantum field theory and e.g. [19, 20] for conformal nets):

1. a1)

(Isotony) The algebras ${\mathfrak{A}}(I)$ are v. Neumann algebras acting on a common Hilbert space $\mathcal{H}$. If $I\subset J$ are intervals, then ${\mathfrak{A}}(I)\subset{\mathfrak{A}}(J)$. 2. a2)

(Causality) Setting $I^{\prime}=\mathbb{R}\setminus[a,b]$ if $I=(a,b)$, we have ${\mathfrak{A}}(I^{\prime})\subset{\mathfrak{A}}(I)^{\prime}$, i.e. observables from disjoint intervals commute. 3. a3)

(Covariance) On $\mathcal{H}$, there is a unitary representation $g\mapsto U(g)$ of the group ${\rm SL}(2,\mathbb{R})/\{\pm 1\}$. If we let elements $g=\left(\begin{matrix}a&b\\ c&d\end{matrix}\right)$ of this group act locally on $\mathbb{R}$ by fractional transformations $g(x)=\frac{ax+b}{dx+c}$, then it is assumed that $U(g){\mathfrak{A}}(I)U(g)^{*}={\mathfrak{A}}(g(I))$ for all intervals $I$ and $g\in{\rm SL}(2,\mathbb{R})$ such that $g(x)$ is well defined for all $x\in I$. We also use the notation

[TABLE] 4. a4)

(Spectrum) The representation $g\mapsto U(g)$ is strongly continuous. The infinitesimal generator $P$ of translations ${\rm tra}_{t}(x)=x+t$, i.e. $P=-i\frac{{\rm d}}{{\rm d}t}U({\rm tra}_{t})|_{t=0}$ has non-negative spectrum. 5. a5)

(Vacuum) There is a unique (unit) vector $|\Omega_{0}\rangle\in\mathcal{H}$ such that $U(g)|\Omega_{0}\rangle=|\Omega_{0}\rangle$. The corresponding state functional will be called $\omega_{0}(X)=\langle\Omega_{0}|X\Omega_{0}\rangle$ throughout. The vacuum should be cyclic for $\bigvee_{I}{\mathfrak{A}}(I)$, the v. Neumann algebra generated by all intervals.

The algebra of observables associated with the union of $p$ open intervals with disjoint closures,

[TABLE]

where each $I_{j}=(a_{j},b_{j})$ is an interval of $\mathbb{R}$ (or arc of the circle ${\mathbb{S}}$ below), is defined to be

[TABLE]

where the symbol $\vee$ means the v. Neumann algebra that is generated by the algebras for the individual arcs/intervals.

3.2 Conformal nets on the circle (compactified lightray)

If we want to insist on a global action of the Möbuis group ${\rm SL}(2,\mathbb{R})/\{\pm 1\}$ on the net, we must pass from the light ray to a compactified lightray, i.e. the circle. The compactification proceeds via the Caley transformation $C:{\mathbb{S}}\setminus\{+1\}\to\mathbb{R},C(x)=-i(x+1)/(x-1)$, and under this transformation intervals get mapped to arcs of the circle. The Caley transform intertwines the action of ${\rm SL}(2,\mathbb{R})/\{\pm 1\}$ on the lightray with the action $z\mapsto g(z)=\frac{\alpha z+\beta}{\bar{\beta}z+\bar{\alpha}}$ of ${\mathrm{SU}(1,1)}/\{\pm 1\}$ on the circle, where $g$ now corresponds to the matrix $\left(\begin{matrix}\alpha&\beta\\ \bar{\beta}&\bar{\alpha}\end{matrix}\right)\in\rm{SU}(1,1)$ under the standard isomorphism between the groups ${\rm SL}(2,\mathbb{R})$ and ${\mathrm{SU}(1,1)}$.

The axioms for a conformal field theory, i.e. net of operator algebras, on the circle are completely analogous to those for the lightray. In the circle picture, it is more standard and natural to use the generators of ${\mathrm{SU}(1,1)}$ called $L_{0},L_{\pm 1}$, where $L_{0}$ is the generator of rotations $z\to e^{it}z,t\in\mathbb{R}$. The requirement a4) is equivalent to the requirement that $L_{0}$ has non-negative spectrum. From a net on the circle, we may via the Caley transform always get a net on the lightray such that $P$ has non-negative spectrum, but not necessarily vice versa since the point at infinity is missing from the lightray. For the rest of the paper, we will assume the axioms on the circle. In sec. 5, we will also need:

1. a6)

(Finite trace) $\operatorname{Tr}e^{-\beta L_{0}}<\infty$ for $\beta>0$.

The above axioms (including the trace condition just mentioned) have a number of well-known consequences which are of interest for this paper:

For each interval ${\mathfrak{A}}(I)^{\prime}={\mathfrak{A}}(I^{\prime})$ (Haag duality [9]). 2. 2.

For each interval, the linear subspace ${\mathfrak{A}}(I)|\Omega_{0}\rangle$ is dense in $\mathcal{H}$ (Reeh-Schlieder theorem). As a consequence, the vector $|\Omega_{0}\rangle$ is cyclic and separating for each local algebra ${\mathfrak{A}}(I)$, and we can apply Tomita-Takesaki theory to the pair $({\mathfrak{A}}(I),|\Omega_{0}\rangle)$. 3. 3.

The modular operator $\Delta$ associated with an open arc $I=(a,b)$ acts geometrically in the sense that

[TABLE]

(Hislop-Longo-theorem [9]). 4. 4.

Each algebra ${\mathfrak{A}}(I)$ has in its central decomposition only hyperfinite type III1 factors [21, 22]. 5. 5.

$\overline{\bigcup_{I}{\mathfrak{A}}(I)}={\mathfrak{B}}(\mathcal{H})$ (irreducibility [23]).

Most of these axioms and results have a more or less obvious counterpart for graded local, i.e. “fermionic”, theories, see e.g. [24].

Example 4: (Virasoro-net) The Virasoro algebra is the Lie-algebra with generators $\{L_{n},\kappa\}_{n\in\mathbb{Z}}$ obeying

[TABLE]

A positive energy representation on a Hilbert space $\mathcal{H}$ is a representation such that (i) $L^{*}_{n}=L_{-n}$ (unitarity), (ii) $L_{0}$ is diagonalizable with non-negative eigenvalues, and (iii) the central element is represented by $\kappa=c1$. From now, we assume a positive energy representation. We assume that $\mathcal{H}$ contains a vacuum vector $|\Omega_{0}\rangle$ which is annihilated by $L_{-1},L_{0},L_{1}$, ($\mathfrak{sl}(2,\mathbb{R})$-invariance) and which is a highest weight vector (of weight 0), i.e. $L_{n}|\Omega_{0}\rangle=0$ for all $n>0$. One has the bound [20, 25, 26, 27]

[TABLE]

for $|\Psi\rangle\in{{\cal V}}\equiv\bigcap_{k\geq 0}{{\cal D}}(L_{0}^{k})\subset\mathcal{H}$ and any natural number $k$.

One next defines from the Virasoro algebra the stress tensor on the unit circle $\SS$, identified with points $z=e^{2\pi iu},u\in\mathbb{R}$ in $\mathbb{C}$. The stress tensor is an operator valued distribution on $\mathcal{H}$ defined in the sense of distributions by the series

[TABLE]

More precisely, for a test function $f\in C^{\infty}(\SS)$ on the circle, it follows from (17) that the corresponding smeared field

[TABLE]

is an operator defined e.g. on the dense invariant domain ${{\cal V}}=\bigcap_{k\geq 0}{{\cal D}}(L_{0}^{k})\subset\mathcal{H}$ (which can be shown to be a common core for the operators $T(f)$) and the assignment $f\mapsto T(f)|\psi\rangle$ is continuous in the topologies on $C^{\infty}(\SS)$ and $\mathcal{H}$ for any vector in this domain. Letting $\Gamma$ be the anti-linear involution

[TABLE]

the smeared stress tensor is a self-adjoint operator on ${{\cal D}}(L_{0})$ for $f$ obeying the reality condition $\Gamma f=f$, and one has $T(f)^{*}=T(\Gamma f)$ in general. It can be shown that the operators $e^{iT(f)}$ for real $f$ form a unitary projective representation of the (covering of the) group of orientation preserving diffeomorphisms (whose generators are the vector fields $f(z){\rm d}/{\rm d}z$) on the circle. The Virasoro net is then defined by

[TABLE]

where the double prime means the v. Neumann closure. The generators of this algebra hence correspond to diffeomorphisms acting trivially outside the arc $I\subset\SS$.

3.3 Pointlike fields

The standard setup of CFT commonly used in the physics literature is based on the use of pointlike fields rather than nets of algebras of bounded operators. Here we will sketch the connection. In fact, the full mathematical details of this connection are not understood in general, although in many important classes of examples, see [20].

On the circle, one typically postulates the existence of local fields having the “mode expansions”

[TABLE]

$h>0$ is called the conformal dimension of the field. The field is typically “energy bounded” i.e. that the modes $\phi_{n},n\in{\mathbb{Z}}$ of the field are linear operators on $\mathcal{H}_{0}$, which satisfy:

Assumption 1**.**

The local fields have a mode expansion (22) such that:

an energy bound of the type $\|(1+L_{0})^{k}\phi_{n}\Psi\|\leq C(1+|n|)^{k+h-\frac{1}{2}}\|(1+L_{0})^{k+h-1}\Psi\|$ for all $n\in{\mathbb{N}}_{0},|\Psi\rangle\in\mathcal{H}_{0}$, and for some $k\geq 0$, satisfying 2. 2.

the commutation relations $[L_{m},\phi_{n}]=((h-1)m-n)\phi_{n+m}$ for $|m|\leq 1$ where $L_{-1},L_{0},L_{1}$ are the generators of the action of ${\mathrm{SU}(1,1)}$ on $\mathcal{H}_{0}$ and where $h\in\mathbb{R}$ is called the conformal spin, satisfying 3. 3.

if $|\Omega_{0}\rangle\in\mathcal{H}_{0}$ is the vacuum vector, then $\phi_{n}|\Omega_{0}\rangle=0$ for $n>-h$, and satisfying 4. 4.

$\phi_{n}^{*}=\phi_{-n}$* for a self-adjoint local field.* 5. 5.

The fields $\phi(f),{\rm supp}(f)\subset I$ should be affiliated with ${\mathfrak{A}}(I)$, i.e. there exists a sequence $B_{n}$ such that $\lim_{n}B_{n}|\Psi\rangle=\phi(f)|\Psi\rangle$ for all $\Psi\in{{\cal V}}=\cap_{k}{{\cal D}}(L_{0}^{k})$.

These properties imply that the smeared fields are operator valued tempered distributions on the domain ${{\cal V}}=\cap_{k}{{\cal D}}(L_{0}^{k})$: Let $|\Psi\rangle\in{{\cal V}}$. Then 1) gives, with $\phi(f):=\int_{\SS}\phi(z)f(z){\rm d}z$,

[TABLE]

because $|\widehat{f}_{-n-h}|$ goes to zero for $|n|\to\infty$ faster than any inverse power. Thus, $\phi(f)$ is an operator valued distribution on the dense invariant domain ${{\cal V}}$, which is in fact a common core for the operators $\phi(f)$. By the same type of estimate the properties 1) and 2) imply furthermore that $\phi(z)|\Omega_{0}\rangle$ can be analytically continued to a $\mathcal{H}$-valued holomorphic function on ${\mathbb{D}}^{+}$ with vector valued distributional boundary value on ${\mathbb{S}}$. It follows from the commutation relations 2) that $\mathcal{H}$ carries a strongly continuous unitary representation $U$ of ${\mathrm{SU}(1,1)}$ generated by $L_{0},L_{\pm 1}$, and this representation satisfies transformation law

[TABLE]

where $g\in\widetilde{{\mathrm{SU}(1,1)}}$ is in the covering group of the Möbius group, and $g(z)$ its action on points $z$ of the circle. For integer $h\in{\mathbb{N}}_{0}$, we get a representation of ${\mathrm{SU}(1,1)}/\{\pm 1\}$. The restriction of $U$ to the invariant subspace ${\rm span}\{\phi_{n}|\Omega_{0}\rangle=0\mid n\leq-h\}$ is a discrete series representation (see e.g. IX, para. 3 of [28]). It also follows that primary fields can and will be normalized so that

[TABLE]

In particular, we see that the field can be local only if the dimension $h$ is a natural number. Fermionic fields are not local but satisfy a graded locality. In that case $h\in\tfrac{1}{2}{\mathbb{N}}_{0}$.

Example 5: (Stress tensor) The stress tensor $T(z)$ affiliated with the Virasoro net of central charge $c>0$ is a pointlike field of dimension $h=2$ satisfying the above assumptions.

Example 6: ($U(1)$-current, see e.g. [25, 24]) The net of the free $U(1)$ current on the circle can be defined e.g. starting from the Lie-algebra generated by a central element $1$ and the “modes” $j_{n},n\in\mathbb{Z}$ defined by $[j_{n},j_{m}]=in\delta_{n,-m}1$ with *-operation $j_{n}^{*}=j_{-n}$. The Hilbert space $\mathcal{H}$ is the closure of the linear span of $j_{n_{1}}\dots j_{n_{k}}|\Omega_{0}\rangle,n_{1}\leq\dots\leq n_{k}\leq-1$ on which the action of the $J_{m}$’s is obtained via the commutation relations and the condition $j_{n}|\Omega_{0}\rangle=0$ for $n>-1$. One sets $L_{n}=\frac{1}{2}\sum_{m\in{\mathbb{Z}}}:j_{n-m}j_{m}:$, where here and in the following, the normal ordering sign $:\ ,\ :$ means that modes with index $m>-1$ (or $>h$ if the field has dimension $h$) are always put to the right of the modes with index $n-m\leq-1$. The $L_{n}$’s the satisfy a Virasoro algebra of central charge $c=1$.

It can be checked that the corresponding current

[TABLE]

satisfies the above assumptions with $h=1$ and is hence an operator valued distribution satisfying $j(z)^{*}=z^{2}J(z)$ and $[j(z),j(w)]=i\delta^{\prime}(z-w)$. For any test-function $f$, the smeared operator $j(f):=\int_{\SS}j(z)f(z){\rm d}z$ has a dense set of analytic vectors (a space of such vectors is spanned by the eigenvectors of $L_{0}$), and hence is essentially self-adjoint by Nelson’s analytic vector theorem. Hence, we can unambiguously define the Weyl operators

[TABLE]

(here $\Gamma f(z)=-\overline{f(z)}$ and $C^{\infty}_{\Gamma}(\SS)$ is the set of invariant elements under $\Gamma$), satisfying the Weyl relations

[TABLE]

The corresponding net of v. Neumann algebras is defined by

[TABLE]

where $I\subset\SS$ is an open arc of the circle or a union thereof, and the double prime means the weak closure (double commutant). The local, unbounded field operators $j(f),{\mathrm{supp}\,}f\subset I$ are not contained in- but are affiliated with these algebras.

Example 7: (Free Fermi net, see e.g. [33, 24]) This net is constructed starting from the Clifford algebra generated by a central element $1$ and the “modes” $\psi_{n},n\in\mathbb{Z}+\frac{1}{2}$ subject to the relations $\psi_{n}\psi_{m}+\psi_{m}\psi_{n}=\delta_{n,-m}1,\psi_{n}^{*}=\psi_{-n}$. The vacuum Hilbert space $\mathcal{H}_{\rm NS}$ is the closure of the linear span of the vectors $\psi_{n_{1}}\cdots\psi_{n_{k}}|\Omega_{\rm NS}\rangle,n_{1}<n_{2}<\dots<0$. A $*$-representation is defined setting $\psi_{n}|\Omega_{\rm NS}\rangle=0$ for all $n\geq 0$ using the relations to define the action of an arbitrary $\psi_{n}$. The state $|\Omega_{\rm NS}\rangle$ is in this context called the “Neveu-Schwarz-vacuum”. The operators $L_{n}$ are defined by $L_{n}=\sum_{m\in{\mathbb{Z}}+\frac{1}{2}}m:\psi_{-m+n}\psi_{m}:$ which generate an action of the Virasoro algebra (in particular of the Lie algebra of ${\mathrm{SU}(1,1)}$ generated by $L_{n},n=-1,0,1$), at central charge $c=\frac{1}{2}$. The corresponding field

[TABLE]

is hence an operator valued distribution. It satisfies $\psi(z)^{*}=z\psi(z)$ and $\psi(z)^{*}\psi(w)+\psi(w)\psi(z)^{*}=\delta(z-w)1$. For any test-function $f$, the smeared operator $\psi(f):=\int_{\SS}\psi(z)f(z){\rm d}z$ is in fact a bounded operator satisfying the canonical anti-commutation relations

[TABLE]

The corresponding net of v. Neumann algebras is defined as the CAR-algebra [30]

[TABLE]

The net of local observables is not a local net but a graded local net, see e.g. [24]. There is another representation of the same net ${\mathfrak{A}}_{\rm Fermi}$, called the “Ramond” representation. It is given by the integer moded expansion

[TABLE]

where the modes satisfy the same relations as before. The Hilbert space $\mathcal{H}_{\rm R}$ is constructed as the linear span of the vectors $\psi_{n_{1}}\cdots\psi_{n_{k}}|\Omega_{\rm R}\rangle,n_{1}<n_{2}<\dots\leq 0$ setting $\psi_{n}|\Omega_{\rm R}\rangle=0$ for all $n>0$ using the relations to define the action of an arbitrary $\psi_{n}$. The Virasoro generators in the Ramond representation $\mathcal{H}_{\rm R}$ are $L_{n}=\sum_{m\in{\mathbb{Z}}}m:\psi_{-m+n}\psi_{m}:$.

4 Modular operators for conformal nets on ${\mathbb{S}}$

In this section we give a first prescription for computing modular operators of chiral conformal nets on ${\mathbb{S}}$ satisfying some natural extra conditions. It is related naturally to the matrix elements $\langle\Omega_{0}|\phi(x)\Delta^{it}\phi(y)|\Omega_{0}\rangle$, but leaves in general certain ambiguities that preclude so far their explicit calculation. This difficulty can be overcome to a certain extent by our second method, presented in sec. 5, more directly related to the matrix element $\langle\Omega|\phi(x)(\operatorname{ln}\Delta)\phi(y)|\Omega\rangle$. Since the material here will form the basis of our discussion in sec. 5, and since the arguments are also of independent interest, we nevertheless present this approach first.

4.1 General results

Quite generally, if ${\mathfrak{M}}$ is a v. Neumann algebra in standard form with cyclic and separating vector $\Omega$, then if $X,Y\in{\mathfrak{M}}$, the Fourier transform

[TABLE]

is well-defined in the sense of a tempered distribution in the variable $s$ – in fact for $Y=X^{*}$, $\Gamma_{X,X^{*}}(s){\rm d}s$ is a positive Radon measure on $\mathbb{R}$, see sec. 5.3 of [13]. By the Fourier inversion theorem, the operator $\Delta^{it}$ is hence fully characterized provided we know $\Gamma_{X,Y}(s)$ for all $X,Y\in{\mathfrak{M}}$ and all $s$, i.e. as a distribution in $s$.

Using the KMS condition (9) after shifting the integration contour from the real axis $\mathbb{R}$ to the line $\mathbb{R}-i$ parallel to the real axis immediately gives

[TABLE]

On the other hand, for $X\in{\mathfrak{M}},Y\in{\mathfrak{M}}^{\prime}$ or vice versa, we get

[TABLE]

using that the modular flow $\sigma^{t}(X)=\Delta^{it}X\Delta^{-it}$ preserves ${\mathfrak{M}},{\mathfrak{M}}^{\prime}$.

We now want to describe how the extra structure of chiral conformal field theory can help to characterize $\Gamma_{X,Y}(s)$. The case we want to consider is the v. Neumann algebra ${\mathfrak{M}}=\mathfrak{A}(A)$, associated with a region $A$ consisting of $p$ open arcs. The vector under consideration is the vacuum, $|\Omega\rangle=|\Omega_{0}\rangle$. It seems that the information is most easily retrieved if instead of bounded operators $X,Y$, we work with point-like unbounded field operators as described in the previous section. We define for a generic primary field $\phi$:

[TABLE]

Here, $x,y\in{\mathbb{S}}$ are to be smeared with test functions on $A$ or $A^{\prime}$, and $\Delta\equiv\Delta_{A}$ is the modular operator in the vacuum state for the multi-interval/arc $A$. The quantity $\Gamma$ should be considered as analogous to (34).

Example 8: For one arc, $A=(a,b)$, the modular flow of a local primary field of dimension $h$ is given by the Hislop-Longo theorem [8, 9] as

[TABLE]

Therefore $\Gamma$ can be found from (37) and (25), giving for $x,y\in(a,b)$

[TABLE]

The example suggests that the behavior of $\Gamma(s,x,y)$ near a boundary point $q_{i}$ of a multi-interval could be $(x-q_{i})^{-h}$. This is supported by the following lemma, formulated in the circle picture.

Lemma 1**.**

Under the assumptions on the CFT given in the previous subsections:

If $f\in C^{\infty}_{0}(\mathbb{R})$ and $\Gamma(f;x,y)\equiv\int\Gamma(s;x,y)f(s){\rm d}s$, then $G(f;x,y)$ is smooth in $x,y\in{\mathbb{S}}$ away from the $2p$ end-points of the $p$ intervals $I_{j}$. Moreover, for any $0<\varepsilon\leq\frac{1}{2}$,

[TABLE]

for some constant $C$ only depending on the end-points. Here, $\sigma=+1$ if both $x,y$ are in $A$, $\sigma=-1$ if $x,y$ are in $A^{\prime}={\mathbb{S}}\setminus\bar{A}$. 2. 2.

We have the KMS condition

[TABLE]

in the sense of distributions for $x,y\in A$. When $x\in A,y\in A^{\prime}$ (or vice versa), we have instead

[TABLE]

Proof.

1. We will compare the quantity $\Gamma$ of an arbitrary multi-arc $A$ to that corresponding to a single arc. First we assume $x,y\in A$. Let $I$ be the largest arc contained in $A$ that is symmetric around $x$, and $J$ the largest interval contained in $A$ symmetric around $y$. If, for example, $c_{j}$ resp. $c_{k}$ are to the left of $x$ resp. $y$, then $I=(c_{j},c_{j}^{-1}x^{2})$ resp. $J=(c_{k},c_{k}^{-1}y^{2})$. We let $\Delta_{A},\Delta_{J},\Delta_{I}$ be the modular operators for the corresponding local algebras. Then we use the well-known operator inequality $\Delta_{I}^{\alpha}\geq\Delta_{A}^{\alpha}$ for $0\leq\alpha\leq 1$. This follows from the fact that $\mathfrak{A}(A)\supset\mathfrak{A}(I)$, which is exploited as follows. Quite generally, let ${\mathfrak{M}}_{i}$ be two v. Neumann algebras on the same Hilbert space with common cyclic and separating vector $|\Omega\rangle$. We let $S_{i}$ be the Tomita operators for ${\mathfrak{M}}_{i}$ with polar decompositions $S_{i}=J_{i}\Delta_{i}^{1/2}$. Note that, if ${\mathfrak{M}}_{2}\subset{\mathfrak{M}}_{1}$, then ${\mathcal{D}}(S_{2})\subset{\mathcal{D}}(S_{1})$. The set ${\mathcal{D}}(S_{1})$ is a Hilbert space called $\mathcal{H}_{1}$ with respect to the inner product (graph norm)

[TABLE]

Letting $I:\mathcal{H}_{1}\to{\mathcal{D}}(S_{1})$ be the identification map, one shows that $I^{-1}{\mathcal{D}}(S_{2})$ is a closed subspace $\mathcal{H}_{2}\subset\mathcal{H}_{1}$ with associated orthogonal projection $P_{2}$. The operators $V_{j}=I^{-1}(1+e^{u}\Delta_{j})^{-1/2}$ are isometries from $\mathcal{H}$ to $\mathcal{H}_{j}$ ($j=1,2$) and their adjoints are $V_{j}^{*}=(1+e^{u}\Delta_{j})^{1/2}IP_{j}$ (with $P_{1}=1$). There follow the relations

[TABLE]

which can already be found in [21].

We multiply the first relation from the right with $X\in{\mathfrak{M}}_{2}$ and from the left with $Y^{*}\in{\mathfrak{M}}_{2}$ and take the expectation value in the state $|\Omega\rangle$. Then we use the second equation and obtain

[TABLE]

The fact that the right side is manifestly non-negative for $X=Y$ implies $(1+e^{u}\Delta_{2})^{-1}\leq(1+e^{u}\Delta_{1})^{-1}$, and that, combined with the operator identity

[TABLE]

for $0<\alpha<1$ gives the claim. Therefore $\Delta^{-\alpha/2}_{I}\Delta^{\alpha/2}_{A}(\Delta^{-\alpha/2}_{I}\Delta^{\alpha/2}_{A})^{*}=\Delta^{-\alpha/2}_{I}\Delta^{\alpha}_{A}\Delta^{-\alpha/2}_{I}\leq 1$ implying $\|\Delta_{A}^{\alpha/2}\Delta_{I}^{-\alpha/2}\|\leq 1$, and similarly for $J$. By the functional calculus, if ${\rm d}E(\lambda)$ is the spectral resolution of $\operatorname{ln}\Delta_{A}$:

[TABLE]

and then by the Cauchy-Schwarz inequality

[TABLE]

Since $\phi(x)|\Omega_{0}\rangle$ can be analytically continued to a $\mathcal{H}_{0}$-valued holomorphic function inside the unit disk ${\mathbb{D}}^{+}$ (by the mode expansion of $\phi$, see the previous subsection), the Hislop-Longo theorem applied to the modular flow $\sigma_{I}^{t}$ of $I$ can similarly be continued to imaginary flow time parameter $t=-i\alpha$ (see eq. (38)). In combination with (25) we thereby obtain

[TABLE]

and similarly for $y$. Since this holds for any pair of end points,

[TABLE]

Now we split the testfunction $f=f_{+}+f_{-}$, where the testfunction $f_{-}$ has support in $(-\infty,c]$, and the testfunction $f_{+}$ has support in $[-c,\infty)$. Then, for the contribution from $f_{-}$, we choose $\alpha=1-\varepsilon$, whereas for the contribution from $f_{+}$, we choose $\alpha=\varepsilon$. As a consequence, we find that (40) holds.

This shows in particular boundedness in $x,y$ of $\Gamma(f,x,y)$ away from the endpoints. A similar estimation can be made for descendant fields, i.e. the derivatives of the fields $\phi(x),\phi(y)$, and this shows smoothness. This finishes the proof of 1) when $x,y\in A$.

To cover the other case, we note that the modular operator for ${\mathfrak{A}}(A)^{\prime}\supset{\mathfrak{A}}(A^{\prime})$ is related to that of ${\mathfrak{M}}_{1}$ by $\Delta_{1}^{\prime}=\Delta_{1}^{-1}$. Then the other case of 1) follow by the same argument, namely, $x,y\in A^{\prime}$ we take $I^{\prime}$ be the largest interval contained in $A^{\prime}$ that is symmetric around $x$, and $J^{\prime}$ the largest interval contained in $A^{\prime}$ symmetric around $y$ and proceed in the same way as before.

2. Eq. (41) follows from the KMS condition (9) for bounded operators $X,Y\in\pi_{0}({\mathfrak{A}}(A))^{\prime\prime}$ because we are assuming about the local fields that for test-functions $f,g$ supported in $A$ there exist $X_{n},Y_{n}\in\pi_{0}({\mathfrak{A}}(A))^{\prime\prime}$ with the property that $\lim_{n}X_{n}|\Omega_{0}\rangle=\phi(f)|\Omega_{0}\rangle$ and $\lim_{n}Y_{n}|\Omega_{0}\rangle=\phi(g)|\Omega_{0}\rangle$ in the strong topology. A similar remark applies to (42). ∎

We next make analytic continuations in the variables $x,y$. This gives us the following: For fixed $x\in A$, and fixed test-function $f(s)$, the function $y\mapsto\Gamma(f;x,y)$ has an analytic extension to a holomorphic function of $y$ inside the unit disk ${\mathbb{D}}^{+}=\{z\in\mathbb{C}\mid|z|<1\}$, by the mode expansion of $\phi$. Similarly, for fixed $y\in A$, the function $x\mapsto\Gamma(s;x,y)$ has an analytic extension to $x$ outside the unit disk ${\mathbb{D}}^{-}=\{z\in\mathbb{C}\mid|z|>1\}$. How to extend to $y\in{\mathbb{D}}^{-}$ or $x\in{\mathbb{D}}^{+}$? The idea is that $\Gamma(f,y,x)$ is analytic in this domain, so we try to paste $\Gamma(s;x,y)$ and $\Gamma(f,y,x)$ together across the boundary of the disk and hope that we get an analytic function that way. This will turn out to be the case on account of the KMS condition which relates the two quantities in precisely the right way.

To this end, we define the following auxiliary quantity for fixed $x\in A,s\in\mathbb{R}$:

[TABLE]

We similarly define, for fixed $y\in A,s\in\mathbb{R}$:

[TABLE]

Note that $H$ implicitly also depends on the choice of $y,s$ and $K$ on the choice of $x,s$, but we suppress this since we are, for the moment, only interested in the dependence of $H$ on $x$ and of $K$ on $y$. Furthermore, note that $H$ and $K$ are, a priori defined only as holomorphic functions on the union ${\mathbb{D}}^{+}\cup{\mathbb{D}}^{-}=\mathbb{C}\setminus{\mathbb{S}}$, i.e. the complex plane minus the circle. On the circle, we define the boundary values from the inside resp. the outside of the disk ($\pm$):

[TABLE]

Then the “jump conditions” (41), (42) imply that

[TABLE]

Thus, both $H$ resp. $K$ are solutions to a Riemann-Hilbert-problem (across the contour ${\mathbb{S}}$). These problems are essentially completely understood, see e.g. [31, 32]. The number and type of solutions depends in general on the specification of the behavior of $H$ resp. $K$ near the boundary points $\{q_{j}\}$ of the multi-arc $A=\cup_{j}(a_{j},b_{j})$ and at infinity, see e.g. [31] (para. 79, pp 230) or [32] (para. 42, pp 420). In the case at hand, this behavior is restricted by (40) and by the mode expansions of the fields.

A minor technical complication arises at this stage due to the fact that, as functions of $s$, both $H$ resp. $K$ are only defined in the distributional sense, and our bound (40) likewise also involves a test-function in $s$. This complication would make the direct application of the results in [32, 31] somewhat cumbersome, so we give an explicit analysis of the implications imposed by the Riemann-Hilbert problem in the case at hand taking into account this complication.

First we define the shorthands

[TABLE]

and

[TABLE]

Notice that $Z_{+}(x(1\mp\varepsilon))=Z_{-}(x(1\mp\varepsilon))\mp i/2$ when $\varepsilon\to 0^{+}$, and that $Z_{+}$ has branch cuts on $A$, while $Z_{-}$ has branch cuts on $A^{\prime}={\mathbb{S}}\setminus\overline{A}$. Then we define

[TABLE]

We have:

Lemma 2**.**

$\widetilde{K}(y)$* is a polynomial in $y$ of degree at most $2(p-1)h$, $\widetilde{H}(x)$ is a polynomial in $x$ of degree at most $2(p-1)h$.*

Proof.

We first consider $\widetilde{K}(y)$. The function $\operatorname{ln}\frac{\Pi_{a}(y)}{\Pi_{b}(y)}=2\pi Z_{+}(y)$ jumps by $+2\pi i$ as $y$ crosses $A$ from the inside of the unit disk to the outside. That jump compensates precisely the jump (41) implied by the KMS condition, so that $\widetilde{K}(y)$ is continuous across $A$. Similarly, (42) implies that $\widetilde{K}(y)$ is continuous across the complement $A^{\prime}$ as a function of $y$, too. Therefore, by the edge-of-the-wedge theorem, $\widetilde{K}$ is an analytic function of $y$ in the entire complex plane minus the boundary points of the intervals. Since $\widetilde{K}$ is a tempered distribution on ${\mathbb{S}}$ which away from the boundary points is the boundary value of an analytic function (from either inside or outside the unit disk), it follows that $\widetilde{K}$ cannot have any essential singularities at the boundary points $y\in\{q_{j}\}$. Indeed, at any given boundary point, say $a_{j}$, since we have a tempered distribution, there exists a natural number $N$, such that, if we multiply $\widetilde{K}(y)$ by $(y-a_{j})^{N}$, we get a continuous function near $a_{j}$ in $y\in{\mathbb{S}}$. By the edge of the wedge theorem, this function, being a boundary value from both inside and outside the disk, must be holomorphic near $a_{j}$. Actually, by (40) we know the factor $(\Pi_{a}(y)\Pi_{b}(y))^{h}$ cancels the potential blow up near $a_{j}$ when $y\in A$, and therefore we can actually choose $N=0$. Since this argument can be repeated for any other boundary point, we learn that the function $\widetilde{K}(y)$ is analytic in $y$ throughout the entire complex plane. A similar statement holds for $x$ and $y$ interchanged and with $\widetilde{K}$ and $\widetilde{H}$ interchanged.

We now establish a bound on the modulus of $\widetilde{K}(y)$ for $|y|\to\infty$. We learn from the mode expansions and properties of the fields that $y^{-2ph+2h}\widetilde{K}(y)$ remains bounded. Thus, we conclude that $|\widetilde{K}(y)|\lesssim|y|^{2hp-2h}$ throughout the entire complex plane for some new constant possibly depending on $x\in A$ and on the test-function $f$. Therefore, $\widetilde{K}(y)$ must for fixed $x\in A$ be a polynomial in $y$ of degree at most $2ph-2h$. We may repeat the same argument with the roles of $x$ and $y$ and of $\widetilde{K}$ and $\widetilde{H}$ reversed, and this finishes the proof. ∎

The next lemma is a straightforward consequence of the preceding two lemmas.

Lemma 3**.**

As a distribution on $(x,y)\in{\mathbb{S}}\times{\mathbb{S}}$ (and $s\in\mathbb{R}$), we have (for $Z_{-}$, see (57)):

[TABLE]

where $\varepsilon^{2h}e^{(\frac{1}{2}-\varepsilon)|s|-s/2}c_{mn}(s)\in L^{1}(\mathbb{R},{\rm d}s)$ for each $\frac{1}{2}\geq\varepsilon>0$, with uniformly bounded $L^{1}$-norm in $\varepsilon$, and where $q_{n}$ are polynomials of degree $n$.

Proof.

We consider the distributional boundary values for $x\to A$ from within ${\mathbb{D}}^{-}$ and for $y\to A$ from within ${\mathbb{D}}^{+}$, respectively in the following expression

[TABLE]

This boundary value prescription coincides with that for $\Gamma(s,x,y)$ and thus the right side is well-defined as a distribution in on $A\times A$ (after smearing in $s$ against a testfunction $f(s)$), by elementary results on products of distributions that are boundary values of analytic functions. By lemma 2, $\widetilde{\Gamma}(f,x,y)$ is a polynomial both in $x$ and in $y$. The inequality (40) and the definition of $Z_{+}$ (57) gives the upper bound

[TABLE]

on this polynomial for all $x,y\in A$. Since the coefficients, $a_{mn}(f)$, of this polynomial can be reconstructed by interpolation from the values $\widetilde{\Gamma}(f,x_{\alpha},y_{\beta})$ for $2h(p-1)$ interpolation points $x_{\alpha}$ and $2h(p-1)$ interpolation points $y_{\beta}$ from $A$, this upper bound also holds for $a_{mn}(f)$. By the well-known duality between the Banach spaces $L^{1}(\mathbb{R})$ and $L^{\infty}(\mathbb{R})$ we can interpret this as saying that $a_{mn}(s)e^{(\frac{1}{2}-\varepsilon)|s|+\frac{1}{2}s}$ is a function in $L^{1}(\mathbb{R})$ with norm bounded from above by $C\varepsilon^{-2h}$. Now set $c_{nm}(s)=a_{nm}(s)e^{s}$ and use the relationship $Z_{+}(x(1\mp\varepsilon))=Z_{-}(x(1\mp\varepsilon))\mp i/2$ when $\varepsilon\to 0^{+}$. Then the proposition follows after expressing $\Gamma$ in terms of $\widetilde{\Gamma}$. ∎

If we wish, we can at this stage take an inverse Fourier transform of $G$ in $s$ to get a general expression for $\langle\Omega_{0}|\phi(x)\Delta^{it}\phi(y)|\Omega_{0}\rangle$. We set

[TABLE]

and we may take the polynomials $q$ in (59) as monomials, again for simplicity of notation. Then we immediately get:

Proposition 1**.**

1) As a distribution in $(x,y)\in A\times A$ (with boundary value prescription $(x,y)\in{\mathbb{D}}^{-}\times{\mathbb{D}}^{+}\to{\mathbb{S}}\times{\mathbb{S}}$ understood)

[TABLE]

where $Z\equiv Z_{-}$ is defined in (57), where $\widehat{c}_{mn}(t)$ is analytic in the strip $\{t\in\mathbb{C}\mid-1<\Im(t)<0\}$. There, it satisfies a bound

[TABLE]

and for real $t$ satisfies the property (in the distributional sense as a boundary value)

[TABLE]

2) We must have:

[TABLE]

in the distributional sense (for $x,y\in A\subset{\mathbb{S}}$), where the bi-variate polynomial $Q$ is as in (67).

Proof.

1. The formula (63) follows directly from lemma 3. In particular, the claimed analyticity and bound (64) follow from the corresponding bounds on $c_{mn}(s)$. The formula (65) follows from the KMS-condition. 2) For $t=0$ we evidently have $\Delta^{it}=1$. This condition gives a non-trivial constraint on the functions $\widehat{c}_{nm}$. Introduce the quantity $Z$ as in (57) and

[TABLE]

Note that $Q(x,y)$ is a polynomial in $x,y$ of degree $2(p-1)$ in each variable. We then get 2). ∎

Remark 1**.**

The domain of analyticity of $\widehat{c}_{mn}$ is large enough to permit us to take the limit $|y|\to\infty$ or $|x|\to 0$. The constraint then confirms that $\widehat{c}_{mn}(t)=0$ when $m,n>2h(p-1)$.

As we will see, in certain special cases eq. (66) and the properties given in proposition 1 suffice to determine $\widehat{c}_{mn}$ uniquely. For instance we will see in subsec. 4.2 that for a free fermion, the information we have obtained uniquely fixes the modular flow. For the $U(1)$-current, the proposition is however already less restrictive, although we are still able to get some results in subsec. 4.3. This is mainly because the polynomial $Q$ is of increasing degree and thus contains more free parameters for fields of higher dimension. Also for this reason, we will introduce in sec. 5 another method.

4.2 Example: Modular flow of free Fermi field in vacuum (NS)-state

As an application of these general results, we find the action of the modular flow of a multi-arc $A$ for the net ${\mathfrak{A}}_{\rm Fermi}$ in the vacuum state (Neveu-Schwarz sector), see example 7. Even though the free Fermi net is not local but graded local (the free Fermi field $\psi$ has $h=\frac{1}{2}$) we can easily adapt, in this simple case, our arguments leading to proposition 1 to fields obeying Fermi-statistics, i.e. fields of dimension $h\in\frac{1}{2}{\mathbb{N}}$. The main change appears in (42), where there is now a pre-factor $-1$ on the right side when $\phi=\psi$ obeys Fermi-statistics. This change propagates to eqs. (52) and (53), where there now appears a pre-factor $-1$ on the second line on the right sides in both equations. Following through this sign change one sees that proposition 1 still holds if we replace (65) by $\widehat{c}_{mn}(t-i)=-\widehat{c}_{nm}(-t)$.

To determine these functions, we may, in this simple case, test the relation (66) with $p$ points (and with $h=\frac{1}{2}$). We pick $\zeta,\eta\in\mathbb{R}$ not equal, and we let $x_{l},y_{l}\in A\subset{\mathbb{S}},k,l=1,\dots,p$ be the pre-images of $\zeta=Z(x_{l})\neq\eta=Z(y_{k})$, where $Z=Z_{-}$ is the function defined by (57), and where $A$ is the union of $p$ open disjoint arcs in ${\mathbb{S}}$ as in (13). Testing the constraint (66) with these points we get for the free Fermi field $\psi$

[TABLE]

We note that $v_{m}(x_{k})=(x_{k})^{m}$ and $v_{n}(y_{l})=(y_{l})^{n}$ where $m,n=0,\dots,p-1$ are $p\times p$ Vandermonde matrices whose determinants

[TABLE]

do not vanish since all the points $x_{l},l=1,\dots,p$ are from disjoint intervals in $A$. Thus, the Vandermonde matrices in (68) may be inverted and therefore $\widehat{c}_{mn}(t)$ is uniquely determined. However, rather than finding the coefficients $Q(x,y)=\sum_{m,n=0}^{p-1}Q_{nm}x^{n}y^{m}$ from (67) and inserting the inverses of the Vandermondians directly, we may observe that one solution to the constraint (68) is of the form

[TABLE]

and this must hence be the unique solution. It is a good check that this solution is also consistent with the general properties of proposition 1 (for $h=\frac{1}{2}$). Substituting the solution into proposition 1 (for $h=\frac{1}{2}$) then gives:

Theorem 1**.**

For the free massless real Fermi field on ${\mathbb{S}}$ and a multi-arc $A=\cup_{j=1}^{p}(a_{j},b_{j})\subset{\mathbb{S}}$, the associated modular flow of the Neveu-Schwarz state is

[TABLE]

where $x,y\in A$, with the usual boundary value prescription ($y$ approached from the within ${\mathbb{D}}^{+}$, $x$ approached from within ${\mathbb{D}}^{-}$) understood.

Since the action of the modular flow $\sigma^{t}$ is of second quantized form on the vacuum Hilbert space, it follows that modular flow is uniquely determined by (71). We now obtain the generator if the flow, thereby making contact with the original analysis due to [10] based on eigenfunctions of the Cauchy kernel.

First, we transform our result from the circle to the lightray via the Caley transformation $C:{\mathbb{S}}\setminus\{+1\}\to\mathbb{R},C(x)=-i(x+1)/(x-1)$. The lightray fields are then related to the circle fields by $\psi_{\mathbb{S}}(x)=\sqrt{C^{\prime}(x)}\psi_{\mathbb{R}}(C(x))$. In terms of the lightray fields, eq. (71) is seen to retain its form, where the arcs of the circle $(a_{j},b_{j})$ become intervals $(C(a_{j}),C(b_{j}))$ of the lightray. By abuse of notation, we can thus work with (71) and pretend that all quantities, such as $x,a_{j},b_{j},\psi,A$ (see (13)) refer to the lightray. Next, we go back from (71) to the Fourier transform (37) using (105). This gives us for $x,y\in A\subset\mathbb{R}$,

[TABLE]

In the last line we have substituted the functions $U_{s}^{k}(x)=(-\Pi_{a}(x)\Pi_{b}(x))^{-\frac{1}{2}}q_{k}(x)e^{isZ_{-}(x)}$ [compare (63)] with the choice $q_{k}(x)=N_{k}^{-1}\prod_{i\neq k}(x-a_{i})$ for the polynomials, where $N_{k}$ is the constant

[TABLE]

given in [10], and we have used the identity $\sum_{k=1}^{p}q_{k}(x)q_{k}(y)=Q(x,y)$ taken from [11] (eq. 2.55). For a function $f\in C^{\infty}_{0}(A,\mathbb{C})$ of compact support on our multi-interval $A$ (13), we next let $G_{A}$ be the restriction of the Cauchy kernel $G(x,y)=(1/2\pi i)(x-y-i0)^{-1}$ to $A$, defining an operator on $L^{2}(A)$ and we let

[TABLE]

which will be identified as the modular hamiltonian on the 1-particle space momentarily. As shown by [11] (eq. 2.11), we have $G_{A}U_{k,s}=(1+e^{-s})^{-1}u_{k,s}$, and the functions $U_{s}^{k}$ in fact give a spectral resolution of the operator $G_{A}$. We can consequently write, with $\psi(f)=\int_{\mathbb{R}}\psi(x)f(x){\rm d}x,f,g\in C^{\infty}_{0}(A,\mathbb{C})$:

[TABLE]

with the usual $L^{2}$-inner product on the right side. The kernel $k_{A}(x,y)$ has been computed in [11] (eqs. 2.72, 2.76), with the result ($Z=Z_{-}$)

[TABLE]

where the sum over $k$ is over all pre-images $y_{j}$ of $Z(y)$ not equal to $y$ itself. Below we will also see that $k_{A}=\operatorname{ln}(G_{A}^{-1}-1)$. 444In the literature [40], such formulas have previously been proven for finite-dimensional fermion algebras.

In order to re-interpret this result on Fock space, it is convenient to give a slightly different, but fully equivalent description of the theory $\{{\mathfrak{A}}_{\rm Fermi}(I)\}$ on the lightray. The $n$-point functions on the lightray are of “quasi-free” form in the sense of [30]

[TABLE]

where the sum is over all perfect matchings in the group of permutations on $n$ elements, and where $G$ is the operator defined by the Cauchy kernel. This operator is a projection which in momentum space corresponds to the multiplication with the characteristic function on $\mathbb{R}_{+}$, i.e. $\widehat{Gf}(k)=1_{(0,\infty)}(k)\widehat{f}(k)$. As shown in [30], this leads to an alternative but equivalent description of $\mathcal{H}$ as the fermionic Fock-space ${{\cal H}}=\oplus_{n}\wedge^{n}{{\cal K}}$ with 1-particle space ${{\cal K}}=\{f\in L^{2}(\mathbb{R})\mid\hat{f}(k)=0,\forall k\leq 0\}$ of square integrable functions $f(x)$ whose Fourier transform $\hat{f}(k)$ is non-zero only for $k\geq 0$. In terms of this Fock-space, the representation of the light ray fields can be written as

[TABLE]

where $a^{*}(g),g\in{{\cal K}}$ are smeared creation operators defined as $a^{*}(g)|\Psi\rangle=|g\wedge\Psi\rangle$ on any $n$-particle state $|\Psi\rangle=|\Psi_{1}\wedge\cdots\wedge\Psi_{n}\rangle\in\wedge^{n}{{\cal K}}\subset\mathcal{H}$.

The 1-particle version of the Reeh-Schlieder theorem implies that it is consistent to introduce on the dense domain ${{\cal D}}(h_{A})=\{Gf\mid f\in C^{\infty}_{0}(A,\mathbb{C})\}\subset{{\cal K}}$ the “1-particle” modular Hamiltonian $h_{A}$ as

[TABLE]

and as a consequence of (90), we can then write the modular flow of ${\mathfrak{A}}_{\rm Fermi}(A)$ for a multi-interval $A=\cup_{j=1}^{p}I_{j}$ in second quantized form as

[TABLE]

In view of (75) the final answer may also be (formally) rewritten as

[TABLE]

where $H_{A}=\tfrac{1}{2}\int_{A\times A}k_{A}(x,y)\psi(x)\psi(y){\rm d}x{\rm d}y$ and $k_{A}$ the kernel of the operator on the right side of (76). Our result for the modular flow is thereby seen to be equivalent to the result for the modular flow found previously by [10]. Our arguments therefore in particular provide a rigorous proof of the result by [10]. A rigorous proof in the case that $A$ is a symmetric $p$-interval was previously given by [33, 34] using different methods. They also give a (slightly) corrected way to write this equation in exponentiated form, see eq. 4.3 of [33].

4.3 Example: Modular flow of $U(1)$-current on ${\mathbb{S}}$

The conformal net ${\mathfrak{A}}_{U(1)}$ for the free $U(1)$ current algebra on the circle was defined in example 6. Via the Caley transform $C:{\mathbb{S}}\setminus\{+1\}\to\mathbb{R},C(x)=-i(x+1)/(x-1)$, one obtains a corresponding net indexed by open intervals $I\subset\mathbb{R}$ of the real line (lightray) or a union thereof. The circle and lightray currents are related by $j_{\mathbb{S}}(x)=C^{\prime}(x)j_{\mathbb{R}}(C(x))$ and thus $j_{\mathbb{R}}(x)^{*}=j_{\mathbb{R}}(x)$. The corresponding lightray Weyl operators satisfy the same relations as on the circle. The two-point function on the lightray is

[TABLE]

and thus takes the same form as on the circle (25) up to the precise form of the boundary value prescription.

We would next like to understand better the modular flow of the net ${\mathfrak{A}}_{U(1)}$ of the free $U(1)$ current algebra. In so far as proposition 1 is concerned, the discussion is actually identical for any bosonic field $\phi$ of dimension $d=1$. First we note that for local fields of conformal dimension $h=1,2,3,\dots$ the method used for the free massless Fermi field to determine $\widehat{c}_{mn}$ is inapplicable since the analog of the Vandermonde matrices, $V_{l}^{n}=(x_{l})^{n},l=1,\dots,p,n=1,\dots,2h(p-1)$, that now appear in the analogue of (68) for general $h$ are no longer square matrices and hence not invertible as the sum over $n,m$ would now go up to $2h(p-1)$ according to (66).

But we can obtain a weaker result for $h=1$ which will follow instantly from the following two lemmas. The first lemma is taken from [11].

Lemma 4**.**

Let $x_{l}\in A\subset{\mathbb{S}},l=1,\dots,p$ ($A$ the union of $p$ open disjoint arcs as in (13)) be the pre-images of $\zeta=Z_{-}(x_{l})$ as in (57). Then

[TABLE]

for all natural numbers $j$ in the range $0\leq j\leq 2p-2$ and all $\zeta\in\mathbb{R}$, where

[TABLE]

From this result, one gets:

Lemma 5**.**

Let $x_{l},y_{l}\in A\subset{\mathbb{S}},k,l=1,\dots,p$ be the pre-images of $\zeta=Z_{-}(x_{l})\neq\eta=Z_{-}(y_{k})$ as in (57). Then

[TABLE]

Proof.

Using the notation introduced in the previous proof (with $Z=Z_{-}$), we have:

[TABLE]

using the previous lemma in the last step. Now it follows from the definition of $Q$ (67) that

[TABLE]

which is also equal to $-\delta_{kl}\prod_{j\neq l}(a_{l}-a_{j})\prod_{i=1}^{p}(a_{l}-b_{i})$. Inserting this identity into (86) completes the proof. ∎

Now let $\phi$ be a bosonic field of dimension $d=1$. For fixed $\zeta$, consider the pre-images $x_{l}\in A,l=1,\dots,p$ of $\zeta=Z_{-}(x_{l})$ inside the $p$ open disjoint arcs (13). We can view the $x_{l}=x_{l}(\zeta)$ as functions of $\zeta$ and form the operator-valued distribution on $\mathbb{R}$ given by

[TABLE]

formally corresponding to the “transformation law” of a primary field of dimension 1. Our first result on the modular flow is

Theorem 2**.**

We have for any dimension 1 fields $\phi$ (e.g. the $U(1)$-current)

[TABLE]

in the sense of distributions in $\eta,\zeta\in\mathbb{R}$.

Proof.

First we apply eq. (25) to $\phi$ and we set $x=x_{k},y=y_{l}$ and sum over $k,l=1,\dots,p$. Then we get from the previous lemma:

[TABLE]

Next we set $x=x_{k},y=y_{l}$ in proposition 1 and sum over $k,l=1,\dots,p$. Then we get using the notation introduced in the previous proofs

[TABLE]

applying lemma 4 in the last step. Now let $f(t)=\sum_{m,n=0}^{2(p-1)}\widehat{c}_{mn}(t)K_{m}K_{n}$. Comparing (91) with (90), we conclude that $f(t-i0)=p[2\sinh(\pi t-i0)]^{-2}$ for real $t$, and hence for all $t$ in the strip $0>\Im t>-1$ by the edge-of-the-wedge theorem, completing the proof. ∎

Remark 2**.**

As is well-known, the $U(1)$ current can be represented on the Fock space of two independent free real Fermion fields $\psi_{1},\psi_{2}$ by $j(x)=i:\psi_{1}\psi_{2}:(x).$ In operator algebraic terms, the $U(1)$-net is a subnet of two copies of the free Fermi net, see [24]. If one could show that there was a unit norm vacuum preserving conditional expectation value from the Fermi algebra for region $A$ to the current algebra of $A$ (as follows from the work by [24] when $A$ is an interval by Haag duality), then the modular flow on the current algebra ${\mathfrak{A}}_{U(1)}(A)$ would be that induced by the flow for the Fermi net (thm. 1), by Takesaki’s theorem, see e.g. sec. 5 of [35]. At present, however, we do not know that such a conditional expectation exists for multi-component regions $A$.

As a test, Thm. 1, eq. (71), can then be applied to $\widetilde{j}(\zeta)$ defined as in (88). Alternatively, we my compute the left side of (71) using the modular flow of the free Fermi field(s) given explicitly in [10, 33]. One sees after a computation that either results are consistent with thm. 2. In the next section, we will again discuss the modular Hamiltonian for the $U(1)$-current from a different perspective.

5 Thermal states

It is possible to analyse the modular flow of a thermal state in a similar manner as for the vacuum state. However, we find it useful to use a variation of the method described in the previous section which, in essence, corresponds to replacing the matrix elements of the function $\Delta^{it}$ with functions closely related to the resolvent, $(\Delta-\lambda)^{-1}$. Such matrix elements will have certain analogous jump properties as the functions $K,H$ introduced above, but their exact form depends on the statistics (i.e. conformal dimension $h$) of the field $\phi(x)$. The cases of fermionic555In the case of fermionic fields, we have graded locality as in the example of the free Fermi field. and bosonic fields are treated separately in subsecs. 5.1, 5.3, respectively.

We will use the parametrization $x=e^{2\pi iu}$ of the circle. Under this map $A$ consists of intervals $\cup_{i=1}^{p}(a_{i},b_{i})\subset(0,1)$. $A^{\prime}$ is as before the interior of the complement. We define, by a slight abuse of notations,

[TABLE]

A Gibbs state is given by the usual formula

[TABLE]

where the trace is taken in the vacuum (i.e. the defining) representation $(\pi_{0},\mathcal{H}_{0})$ of the net. We shall mostly take $\beta$ to be real and positive, and occasionally use $\tau=i\beta/2\pi$, which is the periodicity of the correlation functions in imaginary direction in the coordinate $u$. The general case can be obtained usually by analytic continuation in the end. Then it follows immediately that $\omega_{\beta}$ is a $\beta$-KMS-state on ${\mathfrak{A}}={\mathfrak{B}}(\mathcal{H}_{0})$ relative to the 1-parameter automorphism group of rotations of the circle, i.e. translations in $u$. By the Reeh-Schlieder theorem, the GNS-vector $|\Omega_{\beta}\rangle$ corresponding to $\omega_{\beta}$ is (cyclic and) separating for ${\mathfrak{A}}(A)=\vee_{i=1}^{p}{\mathfrak{A}}(I_{i})$, and we can define a corresponding modular operator $\Delta\equiv\Delta_{\beta,A}$ as in the vacuum situation.

5.1 Fermionic fields

We begin by introducing a variant of the construction in the previous section involving resolvents. For $|\Omega\rangle\equiv|\Omega_{\beta}\rangle$ and $\Delta\equiv\Delta_{\beta,A}$, we set666Here and in the following, we write $\phi(u)$ for $\pi_{\omega}(\phi(u))$, where $\pi_{\omega}$ is the GNS representation of $\omega\equiv\omega_{\beta}$.

[TABLE]

The resolvents in these expressions are well defined if $\xi\in(\tfrac{1}{2},\infty)\cup(-\tfrac{1}{2},-\infty)$ in view of $\Delta>0$, and the analytic continuations in $u$ are justified by the KMS property for the Gibbs state, because $\phi(u)|\Omega\rangle$ is a vector-valued holomorphic function on the strip $\{u\in\mathbb{C}\mid 0<\Im(u)<\beta/2\pi\}$. $\phi$ is assumed to be a hermitian field of conformal dimension $h$ satisfying our assumptions 1. The relation to the modular Hamiltonian follows from the formula

[TABLE]

which trivially holds for positive real numbers $\Delta$ and for positive self adjoint operators in view of the spectral theorem. This immediately gives

[TABLE]

where $u$ has a small negative imaginary part. Thus we should try to find $F$. Our first lemma is the crucial tool expressing the analyticity/jump properties across the real $u$-axis.

Lemma 6**.**

Let $v\in A$ be fixed. If $u\in A^{\prime}$, then

[TABLE]

If $u\in A$, then

[TABLE]

in the sense of distributions. Here $\{\phi(u),\phi(v)\}=\phi(u)\phi(v)+\phi(v)\phi(u)$ is the anti-commutator.

Proof.

We may formally write $F(\xi,u,v)=\int_{\mathbb{R}}[(e^{-s}+1)^{-1}+\xi-\tfrac{1}{2}]^{-1}H(s,u,v){\rm d}s$, with $H(s,u,v)$ defined by analogy with (53), with “$-$” in the second line for fermionic fields. Then the statement of the lemma the formally follows from the properties analogous to (55). This argument is formal because $H(s)$ is defined a priori only after smearing with test functions falling off exponentially, see lemma 1. However, we may also argue directly in the following manner. First we let $u,v\in A$ and compute, with $\lambda=\tfrac{1}{2}-\xi$

[TABLE]

On the other hand, we also have

[TABLE]

Here, we used the relations $\phi(u)|\Omega\rangle=\phi(u)^{*}|\Omega\rangle=S\phi(u)|\Omega\rangle=J\Delta^{\tfrac{1}{2}}\phi(u)\Omega\rangle$, the anti-unitary property of $J$, and the property $J\Delta^{-1}J=\Delta$, which follow from Tomita-Takesaki theory. Adding up these relations imply

[TABLE]

which is equivalent to the statement of the lemma when $u,v\in A$.

In the other case, when $u\in A^{\prime},v\in A$, we use the formula (116). On this formula, we act from the left with $\phi(v)$ and from the left with $\phi(u)$ and take the expectation value in $|\Omega$. Then we obtain a formula for $F$. Since $u\in A^{\prime},v\in A$, we have $[\Delta^{-it}\phi(v)\Delta^{it},\phi(u)]=0$ by Tomita-Takesaki theory, and using this formula to commute the operators inside the expectation values, we get the claim of the lemma when $u\in A^{\prime},v\in A$.

∎

As in the previous section, the lemma shows that $F(\xi,u,v)$ defines a function of $u$ for fixed $v\in A$ that is analytic in the cut strip $\{u\in\mathbb{C}\mid|\Im(u)|<\beta/2\pi,u\notin A\}$. Furthermore, by the KMS condition, it can be checked that $H(s,u-i(\beta-0),v)=-H(s,u+i0,v)$, which is also $=-H(s,u+1+i0,v)$ by construction. Thus, $F(\xi,u,v)$ has the same periodicity as the 2-point function $\langle\Omega|\phi(u)\phi(v)\Omega\rangle$, that is

[TABLE]

This allows us to define $F$ as a function of $u$ on the entire complex plane cut by $A+\mathbb{Z}+\tau\mathbb{Z}$. The limits from below the real axis define hermitian distributional kernels $F(\xi,u-i0,v)$ that are of positive/negative type for $\xi>\tfrac{1}{2}$ resp. $\xi<-\tfrac{1}{2}$, which follows from $\Delta>0$; similarly for the limit from above the real axis. Similar statements hold when fixing $u\in A$ and viewing $F(\xi,u,v)$ as a function of $v$.

As in the vacuum case studied in sec. 4.1, we would next like to have a result like (40) of lemma 1 about the potential singularities of $F$ at the end points of the intervals. Unfortunately, the proof strategy of lemma 1 does not hold in the present case since we have no analogue of the Hislop-Longo theorem for thermal states.

Instead, we will prove first a result comparing the modular operator of a thermal state $\omega_{\beta}$ for the full algebra ${\mathfrak{A}}={\mathfrak{B}}(\mathcal{H}_{0})$ to the modular operator for the partial algebra ${\mathfrak{A}}(A)$. The point is that the former corresponds to rotations of the circle (i.e. translations of the coordinate $u$) and is thus known. The idea is more precisely to apply eq. (46) to the case ${\mathfrak{M}}_{1}={\mathfrak{A}}$, ${\mathfrak{M}}_{2}={\mathfrak{A}}(A)$ (viewed as operator algebras on the GNS-Hilbert space $\mathcal{H}_{\beta}\equiv\mathcal{H}$ of the thermal state $\omega_{\beta}\equiv\omega$), so we can define the modular operators $\Delta_{i}$ of ${\mathfrak{M}}_{i}$ on $\mathcal{H}$ with respect to the cyclic and separating vector $|\Omega\rangle\equiv|\Omega_{\omega}\rangle$.

Lemma 7**.**

Let $s\in\mathbb{R}$, let $u,v\in A=\cup_{i=1}^{p}(a_{i},b_{i})$. Then

[TABLE]

with implicit constant depending on $s$ and the endpoints $\{q_{j}\}$ of the intervals.

Proof.

Let $X,Y\in{\mathfrak{M}}_{2}$. From (46), we get with the notations introduced around (43)

[TABLE]

For $y>0$, we can write

[TABLE]

Therefore, by the spectral calculus applied to $y=e^{s}\Delta_{1}$,

[TABLE]

The key idea is now the following. Suppose that, for $|t|<t_{0}$ and some $t_{0}>0$, we knew that $\sigma_{1}^{t}(X)$ is in ${\mathfrak{M}}_{2}$, so $\Delta_{1}^{it}X|\Omega\rangle=\sigma_{1}^{t}(X)|\Omega\rangle$ is in the domain ${\mathcal{D}}(S_{2})$, so $I^{-1}\Delta_{1}^{it}X|\Omega\rangle$ is in $\mathcal{H}_{2}$, so $(1-P_{2})I^{-1}\Delta_{1}^{it}X|\Omega\rangle=0$. Then we can effectively restrict the range in the integral to $|t|\geq t_{0}$ and drop the $i0$-prescription, and this is the moral reason for the existence of our bound.

An even better estimate is obtained if instead we choose a, say even, real-valued smooth function $h(t)\geq 0$ such that $h(t)=0$ for $|t|<\frac{1}{2}$ and $h(t)=1$ for $|t|\geq 1$, say, and write

[TABLE]

Now we take a test function $f$ compactly supported inside $A$ and we denote the distance of the support of $f$ to the boundary of $A$ by $\delta={\rm dist}(\partial A,{\rm supp}(f))$; of course $\delta>0$. Furthermore, we let $\{X_{n}\}$ be a sequence in ${\mathfrak{M}}_{2}$ converging strongly to $\phi(f)$. Such a sequence exists since we assume that the local fields are affiliated. Then (107) also holds for $X=\phi(f)$. We next wish to use the mode expansion (22) inside (107). Because $\omega\equiv\omega_{\beta}$ is a $\beta$-KMS state on ${\mathfrak{M}}_{1}$ with respect to the automorphic actions of rotations $\alpha_{t}=\alpha_{g_{t}},g_{t}(z)=e^{-i\beta t}z$, the modular flow $\Delta_{1}^{it}$ corresponds to rotations in the sense that $\Delta^{it}_{1}\phi(z)\Delta_{1}^{-it}=e^{-iht\beta}\phi(e^{-i\beta t}z)$ for $z\in{\mathbb{S}}$. Then it is clear that $\Delta^{it}_{1}\phi(f)\Delta_{1}^{-it}$ will remain affiliated with ${\mathfrak{M}}_{2}$ as long as $|t|<\delta/\beta$, so (107) holds with $X=\phi(f)$ and $t_{0}=\delta/\beta$.

In terms of modes

[TABLE]

In particular, since the modular flow preserves the state, we must have $\langle\phi_{m}\Omega|\phi_{n}\Omega\rangle=\delta_{n,m}\|\phi_{n}\Omega\|^{2}=\delta_{n,m}\omega(\phi_{n}^{*}\phi_{n})$. From (iv) of assumption 1 it follows that $\phi_{n}^{*}\phi_{n}\lesssim(1+|n|)^{2h-1}(1+L_{0})^{2k}$ for some $k\geq 0$, and therefore we must have

[TABLE]

for all $n$. Since $\omega\equiv\omega_{\beta}$ is a Gibbs state, we conclude

[TABLE]

Using the definition of $I$ as well as $1-P_{2}\leq 1$, we get for any complex constants $c_{n}$,

[TABLE]

Taking the norm squared of (107) then gives:

[TABLE]

where $f_{n}=\int_{0}^{1}e^{2\pi inu}f(u){\rm d}u$ the Fourier components of $f(u)$, and where $\varphi_{t_{0}}(s)$ is some smooth function which can be chosen to satisfy for large $|s|$ a bound of the form $|\varphi_{t_{0}}(s)|\leq O[\,(1+t_{0}|s|)^{-N}]$ for as large an $N$ as we wish. The Fourier coefficients are trivially bounded by the $L^{1}$-norm of $f$. The bound (112) then gives, altogether

[TABLE]

with implicit constants depending on $s$. We can likewise finde a sequence $\{Y_{n}\}$ be a sequence in ${\mathfrak{M}}_{2}$ converging strongly to $\phi(g)$, and thereby obtain a similar result as (113) replacing $\phi(f)$ by $\phi(g)$. Combining these two results now with (46) and using the Cauchy-Schwarz inequality on the right side of that equation gives

[TABLE]

Letting $f,g$ tend to delta-distributions centered at $u,v$ then by defintion, $\delta_{f}\to{\rm dist}(u,\partial A)\lesssim\sum|u-a_{j}||u-b_{j}|$ and likewise for $g$, and the $L^{1}$ norms remain bounded. This gives the claim of the lemma. ∎

From here, we can get:

Lemma 8**.**

Let $u,v\in A=\cup_{i=1}^{p}(a_{i},b_{i})$. Then

[TABLE]

with implicit constant depending on $\beta,s$ and the endpoints $\{q_{j}\}$ of the intervals.

Proof.

For definiteness, consider $F(\xi,u-i0,v)$. From (105), we get for $\Delta>0$ the identity

[TABLE]

We use this for $\Delta=\Delta_{1}$, the modular operator for the full algebra ${\mathfrak{A}}$. Then $\Delta_{1}^{it}$ generates translations by $-\beta t/2\pi$ of the coordinate $u$. Now we sandwich the above identity between $\langle\Omega|\phi(u-i0)$ and $\phi(v+i0)|\Omega\rangle$. Then if $F_{1}$ is defined as $F$ but with $\Delta_{1}$, we get

[TABLE]

The second term is uniformly bounded in $u,v$ since the integrand is the boundary value of an analytic function that his holomorphic for $-\varepsilon<\Im(t)<0$ with algebraic singularities. The first term behaves as $\sim(u-v)^{-2h}$, since this is the UV-behavior of a thermal 2-point function of a field of conformal dimension $h$, as one may also prove rigorously by decomposing the field into modes and using (iv) of assumption 1. The proof now follows by writing $|F|\equiv|F_{2}|\leq|F_{2}-F_{1}|+|F_{1}|$ and using the preceeding lemma 7 on $|F_{2}-F_{1}|$. We remark that the proof actually shows $F(\xi,u-i0,v)\sim\frac{e^{-i\pi h}}{2\pi}\frac{1}{\xi+\tfrac{1}{2}}(u-v-i0)^{-2h}$ for $u$ near $v$ if the field $\phi$ is in standard normalization (25).

∎

Our assumptions 1 and the fermionic nature of $\phi$ (i.e. $h\in\tfrac{1}{2}{\mathbb{N}}$) imply the operator expansion in anti-commutator form,

[TABLE]

as one can see for instance using the relation between our operator algebraic formalism and Vertex Operator Algebras, see [20]. Here the $O_{n}$ are hermitian bosonic primary fields of conformal dimension $2h-1-n$, and $\delta^{(n)}(u)=\frac{{\rm d}^{n}}{{\rm d}u^{n}}\delta(u)$. We let $\langle O_{n}\rangle=\langle\Omega|O_{n}(u)\Omega\rangle$, which is independent of $u$ for our rotation invariant thermal state $|\Omega\rangle$, and use this relation in lemma 6, to obtain a new version of the jump condition involving the thermal expectation values $\langle O_{n}\rangle$.

We summarize the properties of $F$ (94):

Theorem 3**.**

Let $\phi$ be a fermionic field of conformal dimension $h$ in the normalization (25).

• •

For fixed $v\in A$ and $\xi\in(\tfrac{1}{2},\infty)\cup(-\tfrac{1}{2},-\infty)$, $F(\xi,u,v)$ is a holomorphic function of $u$ on the periodically cut complex plane $\mathbb{C}\setminus(A+\mathbb{Z}+\tau\mathbb{Z})$ with the periodicity $F(\xi,u+1,v)=F(\xi,u+\tau,v)=-F(\xi,u,v)$.

• •

Across the cuts, $u\in A$, $F$ satisfies the jump condition

[TABLE]

• •

Near the end-points $q_{i}$ of the intervals, $|F|\lesssim|q_{i}-u|^{-h}$, as a function of $u$.

• •

Near $v$,

[TABLE]

as a function of $u$.

• •

Analogous properties hold true for $u\leftrightarrow v$.

Now we let

[TABLE]

also equal to the thermal 2-point function of the free Fermi field (see below). Our conventions for the Jacobi $\vartheta$ functions are summarized in the appendix. Then the combination $G(u-u^{\prime})F(\xi,u^{\prime},v)$ is doubly periodic in $u^{\prime}$ for any fixed $v\in A$. Thus, integrating around a contour $\gamma_{\square}$ surrounding the fundamental parallelogram as in fig. 1, we get zero,

[TABLE]

Since $G(u-v)\sim\frac{1}{2\pi i}(u-v)^{-1}$ for $u\sim v$, and $u=v$ is the only singularity of $G$ inside the fundamental parallelogram, we can deform the contour $\gamma_{\square}$ to a contour tightly surrounding the intervals as shown in fig. 1, for $\Im(u)\to 0^{-}$ and $\Re(u)\in A$. Then we collect the residue, use the jump condition in the theorem, as well as $\sum_{n=0}^{2h-1}\langle O_{n}\rangle\,G^{(n)}(u-v-i0)=\langle\Omega|\phi(u)\phi(v)\Omega\rangle$.This gives:

Corollary 1**.**

Let $\phi$ be a fermionic primary field of dimension $h$. $F$ defined as in (94) satisfies (in the distributional sense for $u,v\in A$):

[TABLE]

Here, $\operatorname{\int\!\!\!\!\!--\,}$ denotes a regulated version of the integral described in the following remark.

Remark/Definition: The meaning of the regularized integral operation $\operatorname{\int\!\!\!\!\!--\,}_{A}$ is somewhat subtle because because $F^{\pm}(\xi,u,v)=F(\xi,u\pm i0,v)$ are distributions of $u\in\mathbb{R}$ for fixed $v\in A$. As such, there is no obvious meaning to the integral in cor. 1. To define $\operatorname{\int\!\!\!\!\!--\,}_{A}$, we must remember the origin of the formula from an integration along a contour fitting tightly around $A$, see fig. 1. We consider tempered distributions $\widetilde{F}^{\pm}(\xi,u,v)$ of $u\in\mathbb{R}$ with the property that

1. (a)

$\widetilde{F}^{+}(\xi,u,v)-\widetilde{F}^{-}(\xi,u,v)=F^{+}(\xi,u,v)-F^{-}(\xi,u,v)$, 2. (b)

$\widetilde{F}^{\pm}(\xi,u,v)|_{A^{\prime}}=0$, as a distribution tested with functions of $u$ having compact support in $A^{\prime}$ (thus away from any boundary point), and 3. (c)

$\widetilde{F}^{\pm}(\xi,u,v)|_{A}=0$ satisfies the jump condition (119) of thm. 3 as a distribution tested with functions of $u$ having compact support in $A$ (thus away from any boundary point).

We then define the regulated integral by replacing $F^{-}$ with $\widetilde{F}^{-}$, now integrated over all $u$, which is now well-defined. Cor. 1 holds with this prescription, since (a), (b), (c) are all that is used about $F$ in the vicinity of the real axis in the above proof.

The existence of $\widetilde{F}^{\pm}$ can be seen as follows. In $A$, it is defined to be equal to $F^{\pm}$. In $A^{\prime}$, it is defined to zero. This defines distributions for all $u\in\mathbb{R}$, except for the boundary points $q_{i}$. We wish to define $\widetilde{F}^{\pm}$ by a suitable distributional extension. At the boundary points the unextended distributions have “scaling degree” at most $-h$ (by thm. 3), i.e. finite. Therefore, by standard results (see e.g. [36]), an extension exists and is unique up to addition of $\sum_{q_{j}}\sum_{k\leq h-1}{}^{\pm}A_{j}^{k}\delta^{(k)}(u-q_{j})$. The free parameters must be adjusted in such a way that (c) holds for all $u\in\mathbb{R}$, not just inside $A$ away from the boundary points, and such that $(a)$ still holds. This uniquely determines $\widetilde{F}^{\pm}$, and hence our prescription $\operatorname{\int\!\!\!\!\!--\,}_{A}$.777 We remark that in the expressions considered below in sec. 5.4, our regulated integral has the same meaning as that explained in sec. 4 of [11].

The corollary imposes on $F$ an integral equation at which we will look at in more detail in the following example. The example also suggests that corollary 1 can be used to find the modular operator in the case of a general fermionic primary field of dimension $h$.

5.2 Example: Modular Hamiltonian for thermal state of free Fermi field

Now we analyze what we can learn from thm. 3 and cor. 1 for the free Fermi field, where $h=\tfrac{1}{2}$ and $\{\psi(u),\psi(v)\}=\delta(u-v)1$. For this theory, we can actually prescribe thermal states (93) using either the Neveu-Schwarz or the Ramond representations. For the Neveu-Schwarz (vacuum) sector, the representation is the half-integer moded expansion (30), now denoted as $\psi_{\rm NS}(z)$. For the Ramond-sector, the representation corresponds to an integer moded expansion given in (33) now denoted as $\psi_{\rm R}(z)$. The Gibbs states (93) in the Neveu-Schwarz/Ramond representation correspond, accordingly, to the thermal state vectors $|\Omega_{{\rm NS},\beta}\rangle$ and $|\Omega_{{\rm R},\beta}\rangle$, respectively. The thermal 2-point functions are [37]:

[TABLE]

where the subscript ${\rm X}\in\{{\rm R,NS}\}$ indicates the choice of boundary conditions (sector). For NS, the definition was already given above in (125), whereas for R, we set

[TABLE]

Next we define $F_{\rm X}$ as in (94). By thm. 3 we know its properties, with the only trivial difference that $F_{\rm R}(s,u+1,v)=F_{\rm R}(\xi,u,v)$ in the Ramond sector, where the subscript ${\rm X}\in\{{\rm R,NS}\}$ indicates the choice of boundary conditions. These properties uniquely determine $F_{\rm X}$ in our case (cf. table 1):

[TABLE]

assuming for instance $-\varepsilon<\Im(u)<0$. Here $|A|=\sum_{j}(b_{j}-a_{j})$ is the total length of the intervals,

[TABLE]

and

[TABLE]

etc.

By the well-known periodicity properties and pole/zero structure of the $\vartheta$-functions, it is easy to check that the above formula indeed satisfies the the properties known to hold by thm. 3. For instance, the role of the second line is to give the required jump across $A$, and the role of the pole at $u=v$ in the first line gives the delta-function in (116), also consistent with (141). To show that $F$ as given is the only solution to the properties in thm. 3, we multiply any solution $F$ with the inverse of the second line. This will cancel the jump of $F$ except for the delta function in (116). Thus, the multiplied $F$ can have a first order pole at $u=v$ mod $\mathbb{Z}$ with residue $1/(\xi-\tfrac{1}{2})$, by (116). It cannot have poles, however, at the end-points $q_{i}$ of the intervals, since thm. 3 tells us that $|F|\lesssim|q_{i}-u|^{-1/2}$. Thus, the multiplied function $F$ is meromorphic with a single pole in the fundamental parallelogram, with the same residue and periodicity properties as the first line. There can only be one such function, by well-known results on elliptic functions.

Knowing $F$, we get the kernel of the modular hamiltonian $\langle\Omega_{{\rm X},\beta}|\psi_{\rm X}(u)(\operatorname{ln}\Delta)\psi_{\rm X}(v)\Omega_{{\rm X},\beta}\rangle$ by integrating up (96). Thus, we have, in principle, found this kernel. The (non-trivial) integration over $\xi$ in (96) has recently performed by [39, 38], who have independently arrived at an analogous formula as for $F$ for a closely related quantity, by a different method based on resolvents. It is instructive to see more precisely how this method is related to ours, as it may also shed light on how to proceed in the case of fields with general conformal dimension $h$.

For this purpose, let us set $G_{{\rm X}}(u-v)$ be the thermal 2-point function given by (124), with ${{\rm X}}\in\{{\rm NS,R}\}$. By thm. 3, $F_{{\rm X}}$ has the same periodicity as $G_{{\rm X}}$, so the combination $G_{{\rm X}}(u-v^{\prime})F_{{\rm X}}(\xi,u,v)$ is doubly periodic in $u$ for any fixed $v\in A$. Thus, integrating again around a contour $\gamma_{\square}$ surrounding the fundamental parallelogram as in fig. 1, we get zero as in (122), with $G=G_{X}$ now in that equation. Such an identity also derived by [38] via a different argument using the special properties of the free Fermi field. By contrast, we have so far only used general properties and thm. 3, so our argument works for any CFT.

Now, we deform again the contour $\gamma_{\square}$ to a contour tightly surrounding the cuts $A$ inside the fundamental parallelogram, see fig. 1, use the jump properties and bounds on $F$ as given in thm. 3 and collect the residue. While these steps can still be performed in any CFT according to thm. 3, the free Fermi case is especially simple because the poles of $F$ near the end-points $q_{i}$ of the intervals are of the order $|u-q_{i}|^{\tfrac{1}{2}}$, and hence integrable. Thereby, we immediately get the integral equation

[TABLE]

without any need for a regulator in the integral as in cor. 1. Dropping the subscript “X” for ease of notation and indicating by $G_{A}\equiv G|_{A}$ the restriction of the kernels (124) to $A\times A$, we may view $G_{A}$ as an operator $L^{2}(A)\to L^{2}(A)$. In fact, the anti-commutation relations and the positive nature of the Hilbert space inner product in (124) imply the operator inequalities $0\leq G_{A}\leq 1$. Now we define $F^{\pm}_{A}(\xi)$ the operator defined by the restriction $F(\xi,u\pm i0,v)$ to $A$ in kernel notation, i.e. when $u,v\in A$. We can then also write the above integral equation in operator notation as $G_{A}=(G_{A}+\xi-\tfrac{1}{2})F^{-}_{A}(\xi)$ and solve it as

[TABLE]

when $\xi\in(\tfrac{1}{2},\infty)\cup(-\tfrac{1}{2},-\infty)$. In view of (96), we should then further calculate, in operator notation

[TABLE]

Thus, as a kernel, the modular hamiltonian is given by

[TABLE]

Formally, this may also be expressed as saying $H_{A}=\tfrac{1}{2}\int_{A\times A}k_{A}(x,y)\psi(x)\psi(y){\rm d}x{\rm d}y$ and $k_{A}=\operatorname{ln}(G_{A}^{-1}-1)$ is the hamiltonian appearing in eq. (81). This is in accord with a well known formula [41] for the modular Hamiltonian of the free Fermi field in terms the restricted 2-point function $G_{A}$, which can be derived using the Fock-space structure of theory. In the present case, however, our derivation was rather different and can be paralleled for fermionic fields in general CFTs which are not necessarily equivalent to free field theories. As derived here, our formula holds for thermal states in either the NS or the R sector, choosing for $G$ either one of the 2-point functions (124). Thus, apart from confirming our method, this relation to integral equations might also be useful in the case of more general CFTs where one does not have a priori relations like (132).

5.3 Bosonic fields

We now repeat a similar analysis for a bosonic field $\phi$ of dimension $h\in{\mathbb{N}}$, satisfying our standing assumptions 1. In fact, most results in sec. 5.1 for fermions hold just as well for bosons, with nearly identical proof. A difference, however, arises in the jump condition (98), which involves the anti-commutator of $\phi$, whereas for bosons, we would like to have the commutator. This means that we should work with a different definition of $F$, (94). In the bosonic case, we instead define using the shorthands $|\Omega\rangle\equiv|\Omega_{\beta}\rangle$ and $\Delta\equiv\Delta_{\beta,A}$

[TABLE]

Since $\Delta>0$, $F$ is well defined, in the distributional sense if $s>0$, and we refrain from using a different symbol for $F$ even though it is a different quantity compared to sec. 5.1. The connection to the modular Hamiltonian now follows from the elementary formula

[TABLE]

giving

[TABLE]

Thus, we should determine $F$. Again, the key result is a lemma expressing a jump condition.

Lemma 9**.**

Let $v\in A$ be fixed. If $u\in A^{\prime}$, then

[TABLE]

If $u\in A$, then

[TABLE]

in the sense of distributions. Here $[\phi(u),\phi(v)]=\phi(u)\phi(v)-\phi(v)\phi(u)$ is the commutator.

Proof.

Analogous to the proof of lemma 6. ∎

Our assumptions 1 and the bosonic nature of $\phi$ imply the operator expansion

[TABLE]

which may be used to eliminate the commutator in the jump condition (137). Certain obvious changes also apply to the periodicity of $F$, which is now (102) without the “$-$” sign. In the proof of lemma 8 we now use the formula

[TABLE]

instead of (116), leading now to $F(s,u-i0,v)\sim\frac{e^{-i\pi h}}{2\pi}(u-v-i0)^{-2h}$ for $u$ near $v$. Altogether, this gives the following variant of thm. 3 for bosonic fields:

Theorem 4**.**

• •

For fixed $v\in A$ and $s>0$, $F(s,u,v)$ is a holomorphic function of $u$ on the periodically cut complex plane $\mathbb{C}\setminus(A+\mathbb{Z}+\tau\mathbb{Z})$ with the periodicity $F(s,u+1,v)=F(s,u+\tau,v)=F(s,u,v)$.

• •

Across the cuts, $u\in A$, $F$ satisfies the jump condition

[TABLE]

• •

Near the end-points $q_{i}$ of the intervals, $|F|\lesssim|q_{i}-u|^{-h}$, as a function of $u$.

• •

Near $v$,

[TABLE]

as a function of $u$.

• •

Analogous properties hold true for $u\leftrightarrow v$.

Remark 3**.**

We will see below that the theorem, or alternatively the following cor. 2, may be used to find $F$, and thereby the modular Hamiltonian in view of (135). See sec. 5.4 for examples.

In the case of bosonic fields, it is apparently not straightforward to get an analog of the integral equation in cor. 1. The problem is that $G$, the analogue of (125) should have a simple pole at $u=v$ and be such that $G(u-u^{\prime})F(s,u^{\prime},v)$ is doubly periodic in $u^{\prime}$ for fixed $v\in A$. Since, by contrast to fermionic fields, $F(s,u,v)$ is itself doubly periodic in $u$, the simplest choice – in the absence of any other structural properties of $F$, would be to choose $G$ as doubly periodic, too. However, as is well known there is no such meromorphic function with only one simple pole in the fundamental parallelogram. We will, however, be able to make a similar construction below in the case of the $U(1)$ current, where $F$ has an additional property.

In the general case, we can nevertheless still find an integral equation in the vacuum case, where we can simply set $G(x-y)=(2\pi i)^{-1}(x-y)^{-1}$. Here, we go back to the real line picture in which the theory is living on $\mathbb{R}$ parameterized by $x,y$, and we correspondingly write $F(s,x,y)$, etc. Thm. 4 still applies to the vacuum case: By inspection of the proof, we can take in the end the limit $\beta\to\infty$. A similar type of argument as for cor. 1 now leads to (using $\langle O_{n}\rangle=0$ except when $O_{n}=1$ in the vacuum state):

Corollary 2**.**

$F$* defined as in (133) (for a bosonic primary field $\phi$ of dimension $h$ in standard nomalization (25)) satisfies for the vacuum state $|\Omega\rangle\equiv|\Omega_{0}\rangle$:*

[TABLE]

for $x,y\in A\subset\mathbb{R}$ (lightray picture).

Proof.

For the proof, we first work in the circle picture where $x,y\in\SS$. For our assumptions about the fields, we have for $z\in{\mathbb{D}}^{+}=\{z:|z|<1\}$ that $\phi(z)|\Omega_{0}\rangle=\sum_{n\geq h}z^{n-h}\phi_{n}|\Omega_{0}\rangle$ which is holomorphic in ${\mathbb{D}}^{+}$. Likewise, for $z\in{\mathbb{D}}^{-}=\{z:|z|>1\}$ that $\langle\Omega_{0}|\phi(z)=\sum_{n\geq h}z^{-n-h}\langle\phi_{n}\Omega_{0}|$ which is holomorphic in ${\mathbb{D}}^{-}$ and goes as $|z|^{-2h}$ for $|z|\to\infty$. This implies that, for fixed $x\in A$, $F(s,x,y)$ is a holomorphic function in $\mathbb{C}\setminus A$ decaying like $|y|^{-2h}$ for $|y|\to\infty$. Now we consider $y\in{\mathbb{D}}^{+}$, a small contour $\gamma_{0}$ around $y^{\prime}=0$, and the identity

[TABLE]

We then move the contour $\gamma_{0}$ across the cuts $A$ as in fig. 2, deform it to a very large circle $\gamma_{\infty}$, use the decay of $F(s,x,y^{\prime})$ as well as the jump conditions (140) across the cuts $A$ and collect the residue. The statement then follows after transforming back to the lightray picture. ∎

5.4 Example: Modular Hamiltonian of free $U(1)$-current

The $U(1)$ current, $j$, was introduced in example 6 above. It has dimension $h=1$ and commutator $[j(u),j(v)]=i\delta^{\prime}(u-v)$. The thermal 2-point function is

[TABLE]

Here, $\eta_{1},\eta_{2}$ are the constants appearing in connection with the Weierstrass $\wp$-function, see appendix.

Guessing a – hopefully unique – answer for $F$ just from the properties given in thm. 4 as in the case of the free Fermi field does not seem as straightforward for the current, so we proceed by the more deductive method of integral equations. We have already mentioned that for thermal states, we have not found a general way to obtain an analogue of cor. 1 for bosons. However, for the case of the U(1)-current, there is an additional structural property which helps.

Consider the Weierstrass $\zeta$-function satisfying $\zeta^{\prime}=-\wp$, see appendix. It has a single simple pole at $u=0$ in the fundamental parallelogram but is only quasi-periodic, $\zeta(u+1)=\zeta(u)+\eta_{1},\zeta(u+\tau)=\zeta(u;\tau)+\eta_{2}$. As a consequence, the function $\zeta(u)-\eta_{1}u$ is periodic under $u\to u+1$, and changes by $\eta_{2}-\eta_{1}\tau=-2\pi i$ under $u\to u+\tau$. The combination $\{\zeta(v^{\prime}-v)-(v^{\prime}-v)\eta_{1}\}F(\xi,u,v^{\prime})$ is also not doubly periodic in $v^{\prime}$ for any fixed $u\in A$. However, integrating around a contour $\gamma_{\square}$ surrounding the fundamental parallelogram as in fig. 1,

[TABLE]

still gives zero because of the special property of $F$ that, for $u\in A$,

[TABLE]

which follows from the mode expansion of $j(v)=\sum_{n=-\infty}^{\infty}j_{n}e^{-2\pi inv}$ and the fact that $j_{0}$ is the charge operator, which annihilates any state in the vacuum sector. Hence also $j_{0}|\Omega\rangle=0$ for our thermal state, since it can be viewed as a statistical operator on the vacuum Hilbert space.

If we now evaluate this contour integral as before taking the properties of $F$ in thm. 4 into account, then we get

[TABLE]

Here, $\operatorname{\int\!\!\!\!\!--\,}$ denotes the regulated integral described in and below cor. 1.

We now consider the connection with formulae in the literature for the kernel of the modular Hamiltonian. First, we decompose

[TABLE]

with $\frac{i}{2}C(u-v)=\delta^{\prime}(u-v)$ symmetric (imaginary) and $S(u-v)=-\frac{1}{2\pi}{\rm P.V.}\left(\wp(u-v)-\eta_{1}\right)$ the symmetric (real) part of $G$ as in (144) (P.V. means Cauchy principal value). Since the correlation function satisfies $\langle\Omega|j(f)^{*}j(f)\Omega\rangle\geq 0$, we get $G(\overline{f},f)\geq 0$ for any test-function $f\in C^{\infty}_{0}(A,\mathbb{C})$. Decomposing $f=g+ih$ into a real and imaginary part, we get

[TABLE]

for all $g,h\in C^{\infty}_{0}(A,\mathbb{R})$. We then use the positive definite $\mathbb{R}$-bi-linear form $S$ to define on $C^{\infty}_{0}(A,\mathbb{R})$ an inner product, which we then extend to a hermitian inner product by complex anti-linearity in the first entry. Let this hermitian positive definite sesqui-linear form be called $(\ ,\ )_{S}$. Its completion defines a Hilbert space, ${\mathcal{K}}_{A}$, which is contained in the Sobolev space $W_{0}^{1/2,2}(A)$. It is the 1-particle space of the GNS-Fock space built on $|\Omega\rangle\equiv|\Omega_{\beta}\rangle$. It follows from the definitions that ${\mathcal{K}}_{A}\owns f\mapsto\phi(f)|\Omega\rangle\in{\mathcal{H}}$ is an isometry. Furthermore, since $[1-e^{s}(1-\Delta^{\pm 1})^{-1}]^{-1}$ is a bounded, self-adjoint operator on ${\mathcal{H}}$ for $s>0$, it follows from (133) that the kernel $F^{\pm}(s,x,y)=F(s,x\pm i0,y)$ extends to a bounded quadratic form on ${\mathcal{K}}_{A}$, and hence can be identified with a bounded linear operator on ${\mathcal{K}}_{A}$.

By Riesz’ theorem, there is a self-adjoint operator on this Hilbert space, $\Sigma_{A}$, such that $\frac{i}{2}C(f,g)=(f,\Sigma_{A}g)_{S}$ for all $f,g\in{\mathcal{K}}_{A}$. This operator satisfies $|\Sigma_{A}|\leq 1$, so $|\Sigma_{A}^{-1}|\geq 1$ in view of (149). It follows that for $s>0$, $(e^{-s}-\tfrac{1}{2}+\tfrac{1}{2}\Sigma^{-1}_{A})^{-1}$ exits and is a bounded, self-adjoint operator on ${\mathcal{K}}_{A}$. By construction, the operator $\Sigma_{A}$ is expressible as $(i/2)S_{A}^{-1}C_{A}$, with $S_{A},C_{A},G_{A}$ the operators defined by the restriction of the kernels to $A$. Alternatively,

[TABLE]

Using $C=\delta^{\prime}$ and following the same argument as in sec. 4 of [11] (where the vacuum state was considered), we find that

[TABLE]

in the case of one interval $A=(a,b)$, leading to

[TABLE]

Note that the terms in $[\dots]$ do not depend on $u$, and therefore could be added in the integral (145) in view of (146). Therefore, we can write our integral equation (147) in operator form also as:

[TABLE]

where $F^{+}_{A}(s)$ is the bounded operator on ${\mathcal{K}}_{A}$ corresponding to the kernel $F(s,u,v+i0)$. Since we already know that the operator on the right side can be inverted for $s>0$, it follows that

[TABLE]

is a solution to our integral equation in operator notation. Integrating this operator identity over $s$ as in eq. (135), we obtain

[TABLE]

and thereby in view of eq. (135)

[TABLE]

This formula is sometimes quoted in the literature [40] for the matrix elements of the modular Hamiltonian of a free boson field, and should be thought of as analogous to (132) for the fermion field. A similar derivation also goes through if $A=\cup_{i}(a_{i},b_{i})$ consists of an arbitrary number of intervals. By contrast to the results in the literature, our derivation is based on the rigorously derived integral equations, rather than a formal analog with finite-dimensional quantum systems.

We now sketch how to find an explicit expression for the matrix element (156) of the modular Hamiltonian in the case of one interval. To calculate the logarithm of operators, we need the spectral decomposition of the self-adjoint operator $\Sigma_{A}^{-1}$, which we write as

[TABLE]

Here, $V_{s}$ is a generalized888 More precisely, $V_{s}$ is a linear form on a domain ${\mathcal{D}}\subset{\mathcal{K}}_{A}$ defined by $g\mapsto\operatorname{\int\!\!\!\!\!--\,}_{A}V_{s}(u)g(u){\rm d}u$.

left-eigenfunction of $\Sigma_{A}^{-1}$ with eigenvalue $V_{s}\Sigma_{A}^{-1}=-\coth(s)V_{s}$, and $U_{s}$ satisfies $\operatorname{\int\!\!\!\!\!--\,}_{a}^{b}V_{s}(u)U_{s^{\prime}}(u){\rm d}u=\delta(s-s^{\prime})$. Paralleling the steps in sec. 4 of [11], one can see that ($q_{i}$ ranging over the endpoints $a,b$ of the interval)

[TABLE]

with boundary value prescription $u+i0$ understood, is such a generalized eigenfunction. This eigenfunction is constructed in such a way that is has a multiplicative jump by $e^{s}$ across the interval $(a,b)$ (from the last term), and such that it is doubly periodic in $u$ (cf. table 1). Furthermore, near any one $q_{i}$ of the endpoints $a,b$, it is bounded by $\lesssim|u-q_{i}|^{-1}$ in accordance with thm. 4. The constants $c_{q_{i}}(s)$ are at first adjusted so that

[TABLE]

Then, we have again the relation (145) with $V_{s}$ in place of $F$, and evaluating this contour integral by deforming the contour around the fundamental parallelogram to a tight contour around $(a,b)$ as in fig. 1 then gives the desired eigenvalue equation. That (159) impose no actual loss of generality is demonstrated in the end.

Property (159) is evidently equivalent to

[TABLE]

where

[TABLE]

The constants $c_{q_{i}}(s)$ can further be adjusted so that $-i{\rm sgn}(s)C^{-1}_{A}\overline{V}_{s}=U_{s}$ is a right-eigenvector of $\Sigma_{A}^{-1}$, satisfying the desired normalization $\operatorname{\int\!\!\!\!\!--\,}_{a}^{b}V_{s}(u)U_{s^{\prime}}(u){\rm d}u=\delta(s-s^{\prime})$, see sec. 3 of [11] for further discussion of this point. The $\delta$-function in this expression can only come from the the contribution of the integrals near the boundary points $q_{i}$, which gives a practical way of evaluating these constraints.

Using $\operatorname{ln}[(\coth(s)-1)/(\coth(s)+1)]=-2s$ and $\Sigma_{A}^{-1}U_{s}=-\coth(s)U_{s}$, it immediately follows that

[TABLE]

Therefore in view of (156) and the commutation relation $[j(u),j(v)]=i\delta^{\prime}(u-v)$, the final answer may also be (formally) rewritten as

[TABLE]

where $H_{A}=\int_{A\times A}k_{A}(u,v)j(u)j(v){\rm d}u{\rm d}v$ and $k_{A}(u,v)$ the kernel of the operator on the right side of (162). We refrain here form analyzing in detail the remaining integrals.

This construction can be generalized to the case of $p$ intervals $A=\cup_{j=1}^{p}(a_{j},b_{j})$ in the following way. The ansatz for the generalized eigenfunction is now

[TABLE]

where $|A|=\sum_{j}(b_{j}-a_{j})$. On the $2p$ coefficients $c_{q_{i}}(s)$ we impose (159) and the $p$ constraints

[TABLE]

Of these, only $p-1$ are independent because (159) and the periodicity of $V_{s}(u)$ can be used to show that $\sum_{j}\operatorname{\int\!\!\!\!\!--\,}_{a_{j}}^{b_{j}}V_{s}(u)\,{\rm d}u=0$ by integrating $V_{s}(u)$ around a contour as in fig. 1. A priori, (165), (165) are again a restriction on the possible set of eigenfunctions. But in the end one shows that, in fact, all eigenfunctions of $\Sigma^{-1}_{A}$ satisfy this constraint. They are equivalent to

[TABLE]

where

[TABLE]

This leaves us with $p$ linearly independent left-eigenfunctions, $V^{k}_{s}(u),k=1,\dots,p$. The constants $c_{q_{i}}(s)$ are further be adjusted so that $-i{\rm sgn}(s)C^{-1}_{A}\overline{V_{s}^{k}}=U_{s}^{k}$ or equivalently

[TABLE]

is a right-eigenvector of $\Sigma_{A}^{-1}$, satisfying the desired ortho-normalization $\operatorname{\int\!\!\!\!\!--\,}_{A}V_{s}^{k}(u)U_{s^{\prime}}^{k^{\prime}}(u){\rm d}u=\delta_{kk^{\prime}}\delta(s-s^{\prime})$.

Following an argument given in sec. 5 of [11], let us finally explain why conditions (159), (165), initially only imposed for convenience in order to derive the desired eigenfunctions, impose no actual loss of generality. By expanding the elliptic functions near the end-points of the intervals and using the regularization prescription implicit in $\operatorname{\int\!\!\!\!\!--\,}$ described in the remark after cor. 1, that the limits of the functions $\lim_{s\to 0}U_{s}^{k}(u)=\chi_{k}(u)$ yields a set of functions whose span is equal to the set of indicator functions $\{1_{(a_{j},b_{j})}(u)\}_{j=1,\dots,p}$. Therefore, the orthogonality relation implies that the $p$ eigenfunctions $\{V_{s}^{k}(u)\}_{k=1,\dots,p}$ are already complete, because any eigenfunction must satisfy (165).

Then a similar analysis as in the case of one interval, expanding the bounded operator $F^{+}(s)$ on ${\mathcal{K}}_{A}$ in the generalized basis $U^{k}_{s},V^{k}_{s}$ of $\Sigma_{A}$, gives again (163) with $H_{A}=\int_{A\times A}k_{A}(u,v)j(u)j(v){\rm d}u{\rm d}v$, where

[TABLE]

now. Again, we refrain here form analyzing in detail the remaining integrals.

6 Conclusion

In this work we have studied the modular flows of multi-component regions in chiral CFTs using methods from complex analysis. Our main tool was the KMS-condition built into modular theory, which we have combined with input from CFT such as locality and analyticity. The main general results concern matrix elements of the modular operator of the general type $\langle\Omega|\phi(x)f(\Delta)\phi(y)\Omega\rangle$, where $f$ are various functions such as $\operatorname{ln}$, resolvents, complex powers, and $\phi$ a primary field. The states $|\Omega\rangle$ studied in this work were the vacuum and Gibbs (thermal) states. In some cases, our results for the matrix elements can be expressed in terms of integral equations of Cauchy-type [32, 31], see e.g. cor. 1, 2.

Solving these types of equations was possible in a number of examples, such as free fermions and the $U(1)$-current. To keep the paper at a reasonable length, we have not carried out in this work all the resulting integrations, which would be needed if one is interested in more explicit answers. This would be possible. In the examples considered, we were also able to explain the relationship with known results in the literature [10, 11, 39, 38], and also with certain general, but formal (in the QFT context), methods known specifically for such free field theories, see [40, 41]. In this sense, our methods also serve to make these formal approaches rigorous.

Our analysis is, in principle, not limited to such free CFTs, and in fact e.g., cor. 1, 2, refer to general primary fields in an arbitrary chiral CFT, subject only to certain standard assumptions. It would be interesting to study the implications of these general results further, which we leave for future work [42].

Acknowledgements: Part of this work was carried out while I was visiting IHES, Paris (March 2019), and Instituto di Atomico Balseiro, Bariloche (March 2018). It is a pleasure to thank these institutions for hospitality and for financially supporting those visits. I am grateful to the Max-Planck Society for supporting the collaboration between MPI-MiS and Leipzig U., grant Proj. Bez. M.FE.A.MATN0003. Discussions with H. Casini, P. Fries, M. Huerta, F. Otto, and S. Del Vecchio are gratefully acknowledged.

Appendix A Conventions for elliptic functions

Our conventions for the $\vartheta$-functions are:

[TABLE]

The $\wp$-function is given by

[TABLE]

The Weierstrass $\zeta$-function is a function defined so that $\zeta^{\prime}(u,\tau)=-\wp(u;\tau)$. Explicitly,

[TABLE]

It has the periodicity $\zeta(u+1)=\zeta(u)+\eta_{1},\zeta(u+\tau)=\zeta(u)+\eta_{2}$.