# Logarithmic correlation functions for critical dense polymers on the cylinder

Alexi Morin-Duchesne, Jesper Lykke Jacobsen

► To cite this version:

Alexi Morin-Duchesne, Jesper Lykke Jacobsen. Logarithmic correlation functions for critical dense polymers on the cylinder. SciPost Physics, 2019, 7 (3), pp.040. 10.21468/SciPostPhys.7.3.040. hal-02279088

HAL Id: hal-02279088

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

Submitted on 23 Aug 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
Logarithmic correlation functions
for critical dense polymers on the cylinder

Alexi Morin-Duchesne¹* and Jesper Lykke Jacobsen²,³,⁴†

1 Institut de Recherche en Mathématique et Physique,
Université catholique de Louvain, Louvain-la-Neuve, B-1348, Belgium

2 Laboratoire de Physique de l'Ecole Normale Supérieure, ENS, Université PSL,
CNRS, Sorbonne Université, Université de Paris, Paris, France

3 Sorbonne Université, École Normale Supérieure,
CNRS, Laboratoire de Physique (LPENS), 75005 Paris, France

4 Institut de Physique Théorique, Université Paris Saclay,
CEA, CNRS, F-91191 Gif-sur-Yvette, France

★ alexi.morin-duchesne@uclouvain.be, † jesper.jacobsen@ens.fr

Abstract

We compute lattice correlation functions for the model of critical dense polymers on a semi-infinite cylinder of perimeter $n$. In the lattice loop model, contractible loops have a vanishing fugacity whereas non-contractible loops have a fugacity $\alpha \in (0, \infty)$. These correlators are defined as ratios $Z(x)/Z_0$ of partition functions, where $Z_0$ is a reference partition function wherein only simple half-arcs are attached to the boundary of the cylinder. For $Z(x)$, the boundary of the cylinder is also decorated with simple half-arcs, but it also has two special positions 1 and $x$ where the boundary condition is different. We investigate two such kinds of boundary conditions: (i) there is a single node at each of these points where a long arc is attached, and (ii) there are pairs of adjacent nodes at these points where two long arcs are attached. We find explicit expressions for these correlators for finite $n$ using the representation of the enlarged periodic Temperley-Lieb algebra in the XX spin chain. The resulting asymptotics as $n \to \infty$ are expressed as simple integrals that depend on the scaling parameter $\tau = \frac{x-1}{n} \in (0, 1)$. For small $\tau$, the leading behaviours are proportional to $\tau^{1/4}$, $\tau^{1/4} \log \tau$, $\log \tau$ and $\log^2 \tau$. We interpret the lattice results in terms of ratios of conformal correlation functions. We assume that the corresponding boundary changing fields are highest weight states in irreducible, Kac or staggered Virasoro modules, with central charge $c = -2$ and conformal dimensions $\Delta = -\frac{1}{8}$ or $\Delta = 0$. With these assumptions, we obtain differential equations of order two and three satisfied by the conformal correlation functions, solve these equations in terms of hypergeometric functions, and find a perfect agreement with the lattice results. We use the lattice results to compute structure constants and ratios thereof which appear in the operator product expansions of the boundary condition changing fields. The fusion of these fields is found to be non-abelian.

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

Copyright Alexei Morin-Duchesne and J. L. Jacobsen.
This work is licensed under the Creative Commons
Attribution 4.0 International License.
Published by the SciPost Foundation.

Received 01-08-2019
Accepted 16-09-2019
Published 30-09-2019
doi:10.21468/SciPostPhys.7.3.040

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

Check for
updates
Contents

<table><tr><td>1</td><td>Introduction</td><td>2</td></tr><tr><td>2</td><td>Dense polymers with periodic boundary conditions</td><td>4</td></tr><tr><td>2.1</td><td>Dense polymers and two-point functions</td><td>4</td></tr><tr><td>2.2</td><td>The enlarged periodic Temperley-Lieb algebra</td><td>5</td></tr><tr><td>2.3</td><td>XX Hamiltonian</td><td>8</td></tr><tr><td>3</td><td>Lattice correlators for single entry points</td><td>8</td></tr><tr><td>3.1</td><td>Refined partition functions</td><td>8</td></tr><tr><td>3.2</td><td>Closed-form expressions in the limit m → ∞</td><td>9</td></tr><tr><td>3.3</td><td>Asymptotic behaviour</td><td>11</td></tr><tr><td>4</td><td>Lattice correlators for double entry points</td><td>13</td></tr><tr><td>4.1</td><td>Refined partition functions</td><td>13</td></tr><tr><td>4.2</td><td>Closed-form expressions in the limit m → ∞</td><td>15</td></tr><tr><td>4.3</td><td>Asymptotic behaviour</td><td>17</td></tr><tr><td>5</td><td>Conformal correlators for single entry points</td><td>19</td></tr><tr><td>5.1</td><td>Preliminaries</td><td>19</td></tr><tr><td>5.2</td><td>Differential equation for the four-point function</td><td>20</td></tr><tr><td>5.3</td><td>Solving for the unknowns and verifying smoothness</td><td>23</td></tr><tr><td>6</td><td>Conformal correlators for double entry points</td><td>25</td></tr><tr><td>6.1</td><td>Preliminaries</td><td>25</td></tr><tr><td>6.2</td><td>Representations of the Virasoro algebra at c = −2</td><td>26</td></tr><tr><td>6.3</td><td>Differential equation for the four-point function</td><td>27</td></tr><tr><td>6.4</td><td>Correlation function for the Kac module</td><td>28</td></tr><tr><td>6.5</td><td>Correlation function for the staggered module</td><td>29</td></tr><tr><td>7</td><td>Operator product expansion and structure constants</td><td>31</td></tr><tr><td>7.1</td><td>Analysis for the fields φ</td><td>31</td></tr><tr><td>7.2</td><td>Analysis for the fields μ</td><td>33</td></tr><tr><td>7.3</td><td>Analysis for the field ω</td><td>34</td></tr><tr><td>8</td><td>Discussion and conclusion</td><td>34</td></tr><tr><td></td><td>References</td><td>37</td></tr></table>

1 Introduction

A synthetic presentation of the study of critical phenomena, taking into account some of the main lessons learned during the first half-century following Onsager's celebrated exact solution of the two-dimensional Ising model [1], would run roughly as follows: A physical system standing at a second-order phase transition possesses a continuum limit which is scale invariant [2] and usually also conformally invariant [3, 4]. It is characterised by a set of critical exponents and universal amplitude ratios that define a given universality class, characteristic
of all systems with a given set of symmetries. The exponents are related to the eigenvalues of the dilation operator, and they characterise the algebraic decay of the correlation functions. This operator can be defined either in lattice models (as the transfer matrix or Hamiltonian) or in the quantum field theory describing the continuum limit, so a paramount goal of the theoretician is to diagonalise it. In two dimensions, this task can be efficiently accomplished by identifying suitably tractable lattice models whose spin chain descriptions can be solved exactly by integrability techniques [5–7], or by studying the properties of suitable highest-weight representations of infinite-dimensional conformal algebras [3, 4, 8–10]. The compatibility of the conclusions obtained from these two approaches has been witnessed in thousands of model studies.

As appealing as this outline may appear, a significant proviso was brought forward in the early 1990's [11, 12]: What happens if the dilation operator is not diagonalisable after all? After a lingering start, the importance of this question has become increasingly clear in the last two decades [13, 14] and has lead directly to the detailed study of non semi-simple representations of certain algebras, both in the lattice model [15–17] and conformal field theory (CFT) [18–22] contexts. The latter define what has become known as logarithmic conformal field theory (LCFT) [23], since the correlation functions exhibit both power-law and logarithmic factors. Quite remarkably, the interplay between the lattice models and their continuum limit has turned out to be as tight as ever before, with essentially the same indecomposable structures appearing in either [24].

A necessary condition for this logarithmic behaviour is provided by non-unitarity. While there are various ways of providing this ingredient, a very natural and physically relevant context arises in the study of models formulated in terms of extended degrees of freedom, such as loops and clusters. In this paper, we shall focus on so-called critical dense polymers [25], of which the basic manifestation is a single completely-packed closed curve that fills up the whole lattice. It is the first member, $\mathcal{LM}(1, 2)$, of the family of logarithmic minimal models [26], with its central charge and conformal weights given by

$$c = -2, \quad \Delta_{r,s} = \frac{(2r-s)^2 - 1}{8}. \qquad (1.1)$$

Theories with $c = -2$ are the most well understood logarithmic conformal field theories [27–31]. Other cognate models described by logarithmic CFTs include self-avoiding walks and critical percolation. A common feature of such lattice models is the appearance of cellular algebras [32] of the Temperley-Lieb type.

Logarithms in correlation functions were previously found in various lattice models, including the abelian sandpile model [33–35], critical dense polymers [36–38], critical percolation [39–41] and the Q-state Potts model [42–44]. In many cases, the results were obtained using conformal arguments and verified numerically on a computer. Our goal is to provide new examples of such logarithmic behaviour where the correlators are computed both rigorously from the lattice using the toolbox of integrability, and using the arguments of conformal invariance extended to the logarithmic theories [45].

In a previous paper [38], we have defined several types of two-point boundary correlation functions of critical dense polymers. We established their exact finite-size expressions on a semi-infinite strip of width $n$ and compared the corresponding asymptotic expansions with the field-theoretical predictions. These correlators describe boundaries with defect points, allowing them to be connected by the loops in various ways. Some of the correlators turned out to exhibit logarithmic features, while others did not. The main goal of the present paper is to extend parts of this study to periodic boundary conditions, where now the correlation functions are defined on a semi-infinite cylinder of circumference $n$. An interesting by-product of this modification is that we can now allow for non-contractible loops with a generic weight $\alpha$,
while contractible loops are forbidden (they have the dense polymer weight $\beta = 0$). This is also interesting from the CFT perspective, as the relevant conformal correlation functions involve boundary fields, as in our previous work, but also a bulk field $\psi_{\alpha}(z, z^*)$, which is responsible for assigning the weight $\alpha$ to the non-contractible loops.

In the usual CFT setting, there is a close connection between conformally invariant boundary conditions and highest-weight representations in the bulk theory [46, 47]. Moreover, a number of strong results (see, e.g., [48, 49] in the loop model context) have been obtained from the principle of modular invariance, that is, by comparing the equivalent results of a closed evolution operator (i.e. subject to periodic boundary conditions) acting between two given boundary states with those of an open evolution operator (with given boundary conditions) acting in periodic imaginary time. Our setup similarly replaces the open evolution operator of [38] with a closed one, but we keep the defect points at the boundary. The correct interpretation of our results thus remains within the boundary version of LCFT. The relation between bulk and boundary theories is much more involved in the LCFT setting [50, 51], and in particular it is known that in bulk LCFT primary operators can mix into Jordan cells of rank higher than two [52, 53].

The outline of the paper is as follows. We start out, in Section 2, by recalling the definition of the model of critical dense polymers and its connection with the Temperley-Lieb algebra and the XX spin chain. Due to our setup, we shall need the enlarged periodic version of the Temperley-Lieb algebra, and shall consider the spin chain with periodic boundary conditions. In Sections 3 and 4, we define two types of lattice correlation functions, where each insertion point on the boundary involves respectively a single node or a pair of nearest-neighbour nodes. We find exact expressions for each of these correlators, compute the leading large-$n$ asymptotics in the form of integral formulae, and extract the limiting cases of small and large distances. The same correlation functions are then discussed in the context of LCFT in Sections 5 and 6. Each lattice correlator is understood as a ratio of conformal correlation functions. We use conformal invariance to obtain differential equations for these conformal correlators and find that the solutions to these equations precisely reproduce the exact results of Sections 3 and 4. In Section 7, we study the operator product expansions of the boundary conformal fields and use the lattice results to compute some of the structure constants. The fusion of the corresponding boundary conformal fields distinguishes between operators that mark the start and end of arcs attached to the boundary, and is found to be non-abelian. The paper ends in Section 8 with a discussion of our results and with concluding comments.

# 2 Dense polymers with periodic boundary conditions

## 2.1 Dense polymers and two-point functions

We study the model of critical dense polymers on the $m \times n$ cylinder. We draw it in the plane, as in Figure 1, and choose $n$ to be even. A configuration of the model is the choice of a tile [2] or [3] for each face of the lattice. The boundary conditions are periodic in the horizontal direction, meaning that the left and right ends of the rectangle are identified in the planar representation. The top of the cylinder is decorated exclusively with simple arcs. Labeling the nodes from 1 to $n$, these arcs connect the points in the pairs (1, 2), (3, 4), ..., ($n-1$, $n$). In contrast, the bottom of the cylinder is decorated with a collection of defects and arcs. Their organisation depends on the correlation function that we are studying and is detailed below. A loop configuration $\sigma$ has the weight $w_{\sigma} = \alpha^{n_{\alpha}} \delta_{n_{\beta},0}$ where $n_{\alpha}$ and $n_{\beta}$ are the numbers of non-contractible and contractible loops in the configuration. A collection of arcs that connects
<figure><img src="image_4.png" /><figcaption>Figure 1: Loop configurations on the 10×8 cylinder, with the boundary conditions corresponding to $Z_0$, $Z_{\Pi}(x=6)$ and $Z_{\Pi\Pi}(x=7)$.</figcaption></figure>

two defects of the boundary has weight one. The partition function is defined as

$$Z = \sum_{\sigma} w_{\sigma}. \qquad (2.1)$$

We denote by $Z_0$ the partition function for the model where the lower segment is identical to the top segment of the cylinder and is therefore exclusively decorated with simple arcs. We consider $\alpha \in (0, \infty)$ for which $Z_0 \neq 0$. Likewise, we denote by $Z_{\Pi}(x)$ the partition function wherein the bottom of the cylinder has simple arcs, and two defects inserted in positions 1 and $x$. Finally, we denote by $Z_{\Pi\Pi}(x)$ the partition function where the bottom of the cylinder has two pairs of adjacent defects in the positions (1, 2) and $(x, x+1)$, and simple arcs occupying the other nodes. Examples of configurations for the three partition functions are given in Figure 1. The corresponding two-point correlation functions are then defined as follows:

$$C_{\Pi}(x) = \lim_{m \to \infty} \frac{Z_{\Pi}(x)}{Z_0}, \qquad C_{\Pi\Pi}(x) = \lim_{m \to \infty} \frac{Z_{\Pi\Pi}(x)}{Z_0}. \qquad (2.2)$$

## 2.2 The enlarged periodic Temperley-Lieb algebra

**Definition of the algebra.** This periodic Temperley-Lieb algebra [54–58] is a unital associative algebra that is used to describe many classes of statistical models on periodic geometries. The terminology and conventions vary, and here we work with the enlarged periodic Temperley-Lieb algebra $\mathcal{EPTL}_n(\alpha, \beta)$ defined in [59]:

$$\mathcal{EPTL}_n(\alpha, \beta) = \langle I, \Omega, \Omega^{-1}, e_j; j = 1, \dots, n \rangle. \qquad (2.3)$$

Each element of the algebra is associated to a connectivity diagram drawn inside a rectangle that has periodic boundary conditions in the horizontal direction. For the generators, this identification is

$$I = \begin{tikzpicture}[baseline=(current bounding box.center)]
    
    \node[draw, fill=blue!20, minimum size=20pt] (top) at (0,0) {};
    \node[draw, fill=blue!20, minimum size=20pt] (bottom) at (0,2) {};
    \node at (top) {$\cdots$};
    \node at (bottom) {$\cdots$};
    \node at (0,0) {1};
    \node at (0,2) {n};
    \node at (1,0) {2};
    \node at (1,2) {3};
    \end{tikzpicture}, \quad e_j = \begin{tikzpicture}[baseline=(current bounding box.center)]
    
    \node[draw, fill=blue!20, minimum size=20pt] (top) at (0,0) {};
    \node[draw, fill=blue!20, minimum size=20pt] (bottom) at (0,2) {};
    \node at (top) {$\cdots$};
    \node at (bottom) {$\cdots$};
    \node at (0,0) {1};
    \node at (0,2) {$j$};
    \node at (1,0) {$j+1$};
    \node at (1,2) {$\cdots$};
    \end{tikzpicture}, \quad e_n = \begin{tikzpicture}[baseline=(current bounding box.center)]
    
    \node[draw, fill=blue!20, minimum size=20pt] (top) at (0,0) {};
    \node[draw, fill=blue!20, minimum size=20pt] (bottom) at (0,2) {};
    \node at (top) {$\cdots$};
    \node at (bottom) {$\cdots$};
    \node at (0,0) {1};
    \node at (0,2) {n};
    \node at (1,0) {2};
    \node at (1,2) {3};
    \end{tikzpicture}, \qquad (2.4a)$$

$$\Omega = \begin{tikzpicture}[baseline=(current bounding box.center)]
    
    \node[draw, fill=blue!20, minimum size=20pt] (top) at (0,0) {};
    \node[draw, fill=blue!20, minimum size=20pt] (bottom) at (0,2) {};
    \node at (top) {$\cdots$};
    \node at (bottom) {$\cdots$};
    \node at (0,0) {1};
    \node at (0,2) {n};
    \node at (1,0) {2};
    \node at (1,2) {3};
    \end{tikzpicture}, \quad \Omega^{-1} = \begin{tikzpicture}[baseline=(current bounding box.center)]
    
    \node[draw, fill=blue!20, minimum size=20pt] (top) at (0,0) {};
    \node[draw, fill=blue!20, minimum size=20pt] (bottom) at (0,2) {};
    \node at (top) {$\cdots$};
    \node at (bottom) {$\cdots$};
    \node at (0,0) {1};
    \node at (0,2) {n};
    \node at (1,0) {2};
    \node at (1,2) {3};
    \end{tikzpicture}. \qquad (2.4b)$$

The defining relations of the algebra are

$$e_j^2 = \beta e_j, \quad e_j e_{j+1} e_j = e_j, \quad e_i e_j = e_j e_i \quad (|i-j| > 1), \qquad (2.5a)$$

$$\Omega\Omega^{-1} = \Omega^{-1}\Omega = I, \quad \Omega e_i\Omega^{-1} = e_{i-1}, \quad \Omega^n e_n = e_n\Omega^n, \quad (\Omega^{\pm 1} e_n)^{n-1} = \Omega^{\pm n}(\Omega^{\pm 1} e_n), \qquad (2.5b)$$
the overlaps in the link state basis are encoded in the Gram matrix. To illustrate, in the bases (2.9), the Gram matrices for $n = 2$ and $n = 4$ are

$$ \begin{pmatrix} 0 & \alpha \\ \alpha & 0 \end{pmatrix} \quad \text{and} \quad \begin{pmatrix} 0 & 0 & 0 & \alpha & 0 & 0 \\ 0 & 0 & \alpha & 0 & \alpha & \alpha^2 \\ 0 & \alpha & 0 & 0 & \alpha^2 & \alpha \\ \alpha & 0 & 0 & 0 & 0 & 0 \\ 0 & \alpha & \alpha^2 & 0 & 0 & \alpha \\ 0 & \alpha^2 & \alpha & 0 & \alpha & 0 \end{pmatrix}. \qquad (2.11) $$

**The spin-chain representation X_n.** The representation X_n of EPTL_n(α, 0) is defined on the vector space (C²)⊗n. The generators e₁, ..., eₙ₋₁ are represented by the matrices

$$ X_n(e_j) = \mathbb{I}_{j-1} \otimes \begin{pmatrix} 0 & 0 & 0 & 0 \\ 0 & i & 1 & 0 \\ 0 & 1 & i^{-1} & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} \otimes \mathbb{I}_{n-j-1}, \quad j = 1, \dots, n-1, \qquad (2.12) $$

where $\mathbb{I}_j$ is the identity matrix of size $2^j$. The representatives for $e_n$ and $\Omega$ depend on a twist angle $\phi$,

$$ X_n(e_n) = t \left[ \begin{pmatrix} 0 & 0 & 0 & 0 \\ 0 & i & e^{i\phi} & 0 \\ 0 & e^{-i\phi} & i^{-1} & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix} \otimes \mathbb{I}_{n-2} \right] t^{-1}, \quad X_n(\Omega) = t e^{i\phi\sigma_1^z/2}, \qquad (2.13) $$

where t is the translation operator:

$$ t|v_1\rangle \otimes |v_2\rangle \otimes \cdots \otimes |v_n\rangle = |v_2\rangle \otimes \cdots \otimes |v_n\rangle \otimes |v_1\rangle. \qquad (2.14) $$

One can check that the defining relations (2.5) and (2.6) are satisfied, with the weight $\alpha$ of the non-contractible loops related to the twist angle via the relation

$$ \alpha = 2 \cos \left( \frac{\phi}{2} \right). \qquad (2.15) $$

For $\alpha \in (0, 2]$, $\phi$ is real and in the range $[0, \pi)$, whereas for $\alpha \in (2, \infty)$, $\phi$ is purely imaginary in the range $(i0, i\infty)$.

**Homomorphism and overlaps.** There exists a homomorphism between the standard module $W_{n,0}$ and the spin-chain module $X_n$ of EPTL_n(α, 0). It is defined via the following local maps:

$$ |\underline{\alpha}\rangle = \omega |\uparrow\downarrow\rangle + \omega^{-1} |\downarrow\uparrow\rangle, \quad |\underline{\beta}\rangle = \omega e^{i\phi/2} |\downarrow\uparrow\rangle + \omega^{-1} e^{-i\phi/2} |\uparrow\downarrow\rangle, \quad \omega = e^{i\pi/4}. \qquad (2.16) $$

For link states with more than one arc, the local map is applied multiplicatively to each arc. For instance:

$$ |\underline{\alpha\beta}\rangle = e^{-i\phi/2} |\uparrow\downarrow\uparrow\downarrow\rangle + \omega^{-2} e^{-i\phi/2} |\uparrow\downarrow\downarrow\uparrow\rangle + \omega^2 e^{i\phi/2} |\downarrow\uparrow\uparrow\downarrow\rangle + e^{i\phi/2} |\downarrow\uparrow\downarrow\uparrow\rangle. \qquad (2.17) $$

It is well known that this is indeed a homomorphism, namely for each $v$ in $W_{n,0}$ we have

$$ X_n(e_j)|v\rangle = |e_j v\rangle, \quad X_n(\Omega^{\pm 1})|v\rangle = |\Omega^{\pm 1} v\rangle, \qquad (2.18) $$

where $\alpha$ and $\phi$ are related as in (2.15). The dual states are then defined as $\langle v| = |v\rangle^t|_{\phi \to -\phi}$. The spin-chain overlap between $v$ and $w$ then equals the Gram overlap between $v$ and $w$:

$$ \langle v|w\rangle = v \cdot w, \quad v, w \in W_{n,0}. \qquad (2.19) $$
## 2.3 XX Hamiltonian

The periodic XX Hamiltonian with the twist $\phi$ is defined as

$$H = X_n(\mathcal{H}) = -\left(\sum_{j=1}^{n-1} \sigma_j^+ \sigma_{j+1}^- + \sigma_j^- \sigma_{j+1}^+\right) - e^{i\phi} \sigma_n^+ \sigma_1^- - e^{-i\phi} \sigma_n^- \sigma_1^+. \quad (2.20)$$

It is hermitian and therefore has real eigenvalues. In terms of the fermions

$$c_j = (-1)^{j-1} \left( \prod_{k=1}^{j-1} \sigma_k^z \right) \sigma_j^-, \quad c_j^\dagger = (-1)^{j-1} \left( \prod_{k=1}^{j-1} \sigma_k^z \right) \sigma_j^+, \quad (2.21)$$

it is expressed as

$$H = -\left(\sum_{j=1}^{n-1} c_j^\dagger c_{j+1} + c_{j+1}^\dagger c_j\right) - e^{i\left(\frac{\pi}{2}(n+2S^z+2)+\phi\right)} c_n^\dagger c_1 - e^{-i\left(\frac{\pi}{2}(n+2S^z+2)+\phi\right)} c_1^\dagger c_n, \quad (2.22)$$

where $S^z = \frac{1}{2} \sum_{j=1}^n \sigma_j^z$ is the magnetisation. Applying a Fourier transform allows us to put $H$ in diagonal form:

$$H = - \sum_{k=1}^{n} 2 \cos(\theta_k) \eta_k^\dagger \eta_k, \quad (2.23)$$

where

$$\eta_k = \frac{1}{\sqrt{n}} \sum_{j=1}^{n} e^{i j \theta_k} c_j, \quad \eta_k^\dagger = \frac{1}{\sqrt{n}} \sum_{j=1}^{n} e^{-i j \theta_k} c_j^\dagger, \quad \theta_k = \begin{cases} \frac{2\pi k - \phi}{n} & \frac{n}{2} + S^z \text{ odd,} \\ \frac{2\pi(k - \frac{1}{2}) - \phi}{n} & \frac{n}{2} + S^z \text{ even.} \end{cases} \quad (2.24)$$

A full set of fermionic operators is obtained by taking $k$ in the set $\{1, \dots, n\}$. In what follows, it is however convenient to extend the definition of $\eta_k$ to integer values of $k$ that are negative, using the periodicity properties $\eta_{k+n} = \eta_k$ and $\eta_{k+n}^\dagger = \eta_k^\dagger$. Then, for $\phi \in [0, \pi)$ and $\phi \in (i0, i\infty)$, the groundstate of $H$ lies in the magnetisation sector $S^z = 0$ and is given by

$$|w_0\rangle = \begin{cases} \eta_{(2-n)/4}^\dagger \cdots \eta_{(n-2)/4}^\dagger |0\rangle & \frac{n}{2} \text{ odd,} \\ \eta_{(4-n)/4}^\dagger \cdots \eta_{n/4}^\dagger |0\rangle & \frac{n}{2} \text{ even,} \end{cases} \quad |0\rangle = |\downarrow \cdots \downarrow\rangle. \quad (2.25)$$

This state is also the groundstate of the transfer matrix $T(u)$ for $0 \le u \le \frac{\pi}{2}$. The corresponding eigenvalue $\Lambda_0$ of $T = T(\frac{\pi}{4})$ is [59, 60]

$$\Lambda_0 = \frac{\cos(\frac{\phi}{2})}{2^{n-1}} \prod_{j=1}^{n/2} (1 + \tan x_j) \prod_{j=n/2+1}^{n} (1 - \tan x_j), \quad x_j = \frac{\pi(j - \frac{1}{2}) - \frac{\phi}{2}}{n}. \quad (2.26)$$

# 3 Lattice correlators for single entry points

## 3.1 Refined partition functions

For each loop configuration that contributes to $Z_{\Pi}(x)$, the two defects can connect either via the back or the front of the cylinder. We denote these two possible connections as $\mathcal{L}$ and $\mathcal{A}$ and define the refined partition functions $Z_{\mathcal{L}}(x)$ and $Z_{\mathcal{A}}(x)$. The partition function $Z_{\Pi}(x)$ decomposes as

$$Z_{\Pi}(x) = Z_{\mathcal{L}}(x) + Z_{\mathcal{A}}(x). \quad (3.1)$$
Let us also define two more partition functions: $Z_{\text{C}}(x)$ and $Z_{\text{R}}(x)$. They are defined in a similar way to $Z_{\text{Π}}(x)$, but with the two defects on the lower boundary attached together to become a long arc, in the two possible ways. This means that, in the middle panel of Figure 1, the boundary condition at the lower end of the cylinder is replaced with the following states flipped vertically:

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

Clearly, we have

$$ \begin{pmatrix} Z_{\text{C}}(x) \\ Z_{\text{R}}(x) \end{pmatrix} = \begin{pmatrix} 0 & \alpha \\ \alpha & 0 \end{pmatrix} \begin{pmatrix} Z_{\text{C}}(x) \\ Z_{\text{R}}(x) \end{pmatrix}. \qquad (3.3) $$

The matrix on the right side is the Gram matrix for $W_{2,0}$ given in (2.11). We also note that the state attached to the top of the cylinder in Figure 1 is none other than $v_2^a$. Its translation by one node is $v_n^b$.

We define the refined correlation functions corresponding to C and R:

$$ C_{\text{C}}(x) = \lim_{m \to \infty} \frac{Z_{\text{C}}(x)}{Z_0}, \qquad C_{\text{R}}(x) = \lim_{m \to \infty} \frac{Z_{\text{R}}(x)}{Z_0}. \qquad (3.4) $$

In constrast to $C_{\text{Π}}(x)$, these two functions have a well-defined $\alpha \to 0$ limit. Indeed, the numerators and denominators are polynomials in $\alpha$ with the lowest-degree term proportional to $\alpha$. In the limit $\alpha \to 0$, these correlation functions give ratios of partition functions where the configurations with non-zero weights have a single non-contractible loop, and no contractible loops.

The partition functions $Z_{\text{C}}(x)$ and $Z_{\text{R}}(x)$ are expressed in terms of Gram overlaps as

$$ Z_{\text{C}}(x) = 2^{mn/2} (v_x^a \cdot T^m v_2^a), \qquad Z_{\text{R}}(x) = 2^{mn/2} (v_x^b \cdot T^m v_2^b), \qquad (3.5) $$

where the factor of $2^{mn/2}$ ensures that the tiles have weight 1 and not $1/\sqrt{2}$, as they do in (2.7) for $u = \pi/4$. Because of the invariance under translation, we can equivalently write

$$ Z_{\text{C}}(x) = 2^{mn/2} (v_{n+2-x}^a \cdot T^m v_n^b). \qquad (3.6) $$

These partition functions are then expressed as overlaps in the spin chain:

$$ Z_{\text{C}}(x) = 2^{mn/2} \langle v_x^a | T^m | v_2^a \rangle, \qquad Z_{\text{R}}(x) = 2^{mn/2} \langle v_{n+2-x}^a | T^m | v_n^b \rangle, \qquad (3.7) $$

where $T = X_n(T)$ is the transfer matrix.

## 3.2 Closed-form expressions in the limit $m \to \infty$

In the limit $m \to \infty$, the leading contribution to the overlaps comes from the groundstate,

$$ Z_{\text{C}}(x) \approx 2^{mn/2} \Lambda_0^m \langle v_x^a | w_0 \rangle \langle w_0 | v_2^a \rangle, \qquad Z_{\text{R}}(x) \approx 2^{mn/2} \Lambda_0^m \langle v_{n+2-x}^a | w_0 \rangle \langle w_0 | v_n^b \rangle, \qquad (3.8) $$

where $\approx$ indicates an equality up to terms that are exponentially small in $m$ compared to the leading term. Because $|w_0\rangle$ is invariant under translation, we know that $\langle w_0 | v_2^a \rangle = \langle w_0 | v_n^b \rangle$. The partition function $Z_0$ is equal to $Z_{\text{C}}(2)$, and as a result we have

$$ C_{\text{C}}(x) = \frac{\langle v_x^a | w_0 \rangle}{\langle v_2^a | w_0 \rangle}, \qquad C_{\text{R}}(x) = \frac{\langle v_{n+2-x}^a | w_0 \rangle}{\langle v_2^a | w_0 \rangle} = C_{\text{C}}(n+2-x). \qquad (3.9) $$
We therefore focus on computing $C_{\phi}(x)$. Following (2.16), the state $v_x^a$ is represented in the spin chain by

$$ \langle v_x^a | = \langle 0 | a_{n-1} a_{n-3} \cdots a_{x+1} a_{x-2} a_{x-4} \cdots a_2 (\omega c_1 + (-1)^{\frac{x-2}{2}} \omega^{-1} c_x), \quad (3.10) $$

where

$$ a_j = \omega c_j + \omega^{-1} c_{j+1}, \quad \omega = e^{i\pi/4}. \quad (3.11) $$

The sign $(-1)^{\frac{x-2}{2}}$ of the last factor in (3.10) comes from anticommuting $c_x$ with the factors $a_2, a_4, \dots, a_{x-2}$. We have the following identity,

$$ (\omega c_1 + (-1)^{\frac{x-2}{2}} \omega^{-1} c_x) = \sum_{\ell=1}^{x-1} i^{\ell-1} a_{\ell}, \quad (3.12) $$

and the following commutation relation:

$$ \{a_j, \eta_k^\dagger\} = 2n^{-1/2} \cos\left(\frac{\theta_k}{2} + \frac{\pi}{4}\right) e^{-i(j+1/2)\theta_k}. \quad (3.13) $$

Using Wick's theorem, we find

$$ \langle v_x^a | w_0 \rangle = \sum_{\ell=1}^{x-1} i^{\ell-1} \det M^{(\ell)}, \quad (3.14a) $$

where

$$ M_{jk}^{(\ell)} = \begin{cases} \{a_\ell, \eta_k^\dagger\} & j=1, \\ \{a_{2j-2}, \eta_k^\dagger\} & j=2, \dots, \frac{x}{2}, \\ \{a_{2j-1}, \eta_k^\dagger\} & j=\frac{x+2}{2}, \dots, \frac{n}{2}, \end{cases} \quad k \in K = \begin{cases} \{\frac{4-n}{4}, \dots, \frac{n}{4}\} & \frac{n}{2} \text{ even,} \\ \{\frac{2-n}{4}, \dots, \frac{n-2}{4}\} & \frac{n}{2} \text{ odd.} \end{cases} \quad (3.14b) $$

For $\ell$ even, the determinant is zero, as the first row is identical to the row with label $j = \frac{\ell+2}{2}$. The factors appearing in (3.13) that are constant or depend only on $k$ can be extracted from the denominator, and we find

$$ C_{\phi}(x) = \frac{\langle v_x^a | w_0 \rangle}{\langle v_2^a | w_0 \rangle} = \sum_{\ell=1,3,\dots,x-1} (-1)^{(\ell-1)/2} \frac{\det \tilde{M}^{(\ell)}}{\det \tilde{M}^{(1)}}, \quad (3.15a) $$

$$ \tilde{M}_{jk}^{(\ell)} = \begin{cases} e^{-i(\ell+1/2)\theta_k} & j=1, \\ e^{-i(2j-3/2)\theta_k} & j=2, \dots, \frac{x}{2}, \\ e^{-i(2j-1/2)\theta_k} & j=\frac{x+2}{2}, \dots, \frac{n}{2}, \end{cases} \quad k \in K. \quad (3.15b) $$

The inverse of the matrix $\tilde{M}^{(1)}$ has the entries

$$ (\tilde{M}^{(1)})^{-1}_{kj} = \frac{2}{n} e^{i\theta_k(2j-1/2)}, \quad j=1, \dots, \frac{n}{2}, \quad k \in K. \quad (3.16) $$

The next step is to multiply $\tilde{M}^{(\ell)}$ by $(\tilde{M}^{(1)})^{-1}$ from the right:

$$ \begin{bmatrix} \tilde{M}^{(\ell)} & (\tilde{M}^{(1)})^{-1} \end{bmatrix}_{jk} = \begin{cases} \delta_{i,(\ell+1)/2} & j=1, \\ -\frac{2(-1)^{j+k} e^{i\phi(2j-2k-1)/n}}{n \sin(\frac{\pi}{n}(2j-2k-1))} & 2 \le j \le \frac{x}{2}, \quad k=1, \dots, \frac{n}{2}. \\ \delta_{j,k} & \frac{x+2}{2} \le j \le \frac{n}{2}, \end{cases} \quad (3.17) $$
This holds for both parities of $\frac{n}{2}$. The presence of the Kronecker-$\delta$ functions stems from the fact that the rows of $\tilde{M}^{(\ell)}$ and $\tilde{M}^{(1)}$ with labels $j = 1$ and $j = \frac{x+2}{2}, \dots, \frac{n}{2}$ are identical. The determinant of $\tilde{M}^{(\ell)}(\tilde{M}^{(1)})^{-1}$ thus reduces to the determinant of a matrix of size $\frac{x-2}{2}$. Except for the $\sin(\frac{\pi}{n}(2j-2k-1))$ in the denominator, all the other factors can be extracted from the determinant. After simplification, we find

$$C_{\phi}(x) = \left(\frac{-4i}{n}\right)^{\frac{x-2}{2}} e^{-\frac{i\phi(x+2)}{2n}} e^{\frac{i\pi(x-1)(x+2)}{2n}} \sum_{\ell=1}^{x/2} (-1)^{\ell-1} e^{\frac{2i\ell(\phi-\pi)}{n}} \det N^{(\ell)} \quad (3.18a)$$

$$N_{jk}^{(\ell)} = \begin{cases} (w_j - z_k)^{-1} & 1 \le k \le \ell - 1, \\ (w_j - z_{k+1})^{-1} & \ell \le k \le \frac{x-2}{2}, \end{cases} \quad j = 1, \dots, \frac{x-2}{2}, \quad (3.18b)$$

where $w_j = e^{4i\pi(j-\frac{1}{2})/n}$ and $z_k = e^{4i\pi k/n}$. The determinants are evaluated using Cauchy's identity:

$$\det_{j,k=1}^{a} \left( \frac{1}{w_j - z_k} \right) = \frac{\prod_{1 \le j < k \le a} (w_k - w_j)(z_j - z_k)}{\prod_{j,k=1}^{a} (w_j - z_k)}. \quad (3.19)$$

This gives us a closed-form expression for $C_{\phi}(x)$ in terms of a complicated sum, wherein the summand is written as a ratio of trigonometric functions. Many factors are independent of $\ell$ and can be factored from the sum. After simplification, we find

$$C_{\phi}(x) = \left(\frac{-2}{n}\right)^{\frac{x-2}{2}} \frac{\prod_{1 \le j < k \le (x-2)/2} \sin\left(\frac{2\pi}{n}(k-j)\right) \prod_{1 \le j < k \le x/2} \sin\left(\frac{2\pi}{n}(j-k)\right)}{\prod_{j=1}^{(x-2)/2} \prod_{k=1}^{x/2} \sin\left(\frac{2\pi}{n}(j-k + \frac{1}{2})\right)} \times \sum_{\ell=1}^{x/2} \cos\left(\frac{\phi}{n}\left(2\ell - \frac{x+2}{2}\right)\right) \frac{\prod_{j=1}^{(x-2)/2} \sin\left(\frac{2\pi}{n}\left(j - \ell + \frac{1}{2}\right)\right)}{\prod_{\substack{j=1 \\ j \ne \ell}}^{x/2} \sin\left(\frac{2\pi}{n}(\ell - j)\right)}. \quad (3.20)$$

### 3.3 Asymptotic behaviour

We compute the asymptotic behaviour of (3.20) as $n \to \infty$ with the ratio of $x$ and $n$ fixed. We achieve this by setting $\tau = \frac{x-1}{n}$ and sending $n$ to infinity with $\tau$ fixed in the range $(0, 1)$. As a first step, we write all the sine functions in terms of the cardinal sine function $s[x] = \frac{\sin(x)}{x}$. The products of integers that appear are then rewritten in terms of the Gamma and Barnes functions:

$$\Gamma(z+1) = z \Gamma(z), \quad G(z+1) = \Gamma(z)G(z). \quad (3.21)$$

The result simplifies to

$$C_{\phi}(x) = \frac{G(\frac{x}{2})G(\frac{x+2}{2})G^2(\frac{1}{2})}{G^2(\frac{x+1}{2})} \prod_{j=0}^{(x-2)/2} \frac{s[\frac{2\pi}{n}(j+\frac{1}{2})]}{s[\frac{2\pi}{n}j]} \prod_{j=0}^{(x-2)/2} \prod_{k=0}^{(x-2)/2} \frac{s[\frac{2\pi}{n}(j-k)]}{s[\frac{2\pi}{n}(j-k-\frac{1}{2})]} \times \sum_{\ell=0}^{(x-2)/2} \frac{\cos(\frac{\phi}{n}(2\ell - \frac{x-2}{2})) \Gamma(\ell + \frac{1}{2}) \Gamma(\frac{x-1}{2} - \ell)}{s[\frac{2\pi}{n}(\ell + \frac{1}{2})] \Gamma(\ell + 1) \Gamma(\frac{x}{2} - \ell)} \prod_{j=0}^{(x-2)/2} \frac{s[\frac{2\pi}{n}(j - \ell - \frac{1}{2})]}{s[\frac{2\pi}{n}(j - \ell)]}, \quad (3.22)$$
where we use the convention $s[0] = 1$. The large-$z$ asymptotics of $G(z)$ is

$$ \log G(z) = \left( \frac{(z-1)^2}{2} - \frac{1}{12} \right) \log(z-1) - \frac{3(z-1)^2}{4} + \frac{z-1}{2} \log(2\pi) + \frac{1}{12} - \log A + \mathcal{O}(z^{-1}), \quad (3.23) $$

where $A$ is the Glaisher-Kinkelin constant. This yields

$$ \frac{G(\frac{n\tau+1}{2})G(\frac{n\tau+3}{2})}{G^2(\frac{n\tau+2}{2})} = \left(\frac{n\tau}{2}\right)^{1/4} + \mathcal{O}(n^{-3/4}). \quad (3.24) $$

To compute the asymptotics of the products in (3.22), we use the $1/n$ expansions

$$ \sum_{j=0}^{(n\tau-1)/2} \log s\left[\frac{2\pi(j-a)}{n} - b\right] = n \int_0^{\tau/2} dx \log s\left[2\pi x - b\right] + \left(a + \frac{1}{2}\right) \left(\log s[b] - \log s[\pi\tau - b]\right) + \mathcal{O}(n^{-1}), \quad (3.25a) $$

$$ \sum_{j,k=0}^{(n\tau-1)/2} \log s\left[\frac{2\pi(j-k-a)}{n}\right] = n^2 \int_0^{\tau/2} dx dy \log s\left[2\pi(x-y)\right] + \left(a^2 - \frac{1}{6}\right) \log s[\pi\tau] + \mathcal{O}(n^{-1}), \quad (3.25b) $$

which are obtained from the Euler-Maclaurin formula. This yields

$$ \prod_{j=0}^{(x-2)/2} \frac{s\left[\frac{2\pi}{n}\left(j + \frac{1}{2}\right)\right]}{s\left[\frac{2\pi}{n}j\right]} = s[\pi\tau]^{1/2} + \mathcal{O}(n^{-1}), \quad \prod_{j,k=0}^{(x-2)/2} \frac{s\left[\frac{2\pi}{n}\left(j - k\right)\right]}{s\left[\frac{2\pi}{n}\left(j - k - \frac{1}{2}\right)\right]} = s[\pi\tau]^{-1/4} + \mathcal{O}(n^{-1}). \quad (3.26) $$

To evaluate the asymptotic behaviour of the sum in (3.22), we first change the summation index to $z = \frac{\ell}{n}$, so that the sum goes from 0 to $\tau$ with increments of $\frac{1}{n}$. We have the following asymptotics

$$ \prod_{j=0}^{(x-2)/2} \frac{s\left[\frac{2\pi}{n}\left(j - \frac{1}{2}\right) - 2\pi z\right]}{s\left[\frac{2\pi j}{n} - 2\pi z\right]} = \frac{s\left[2\pi z\right]^{1/2}}{s\left[\pi(\tau - 2z)\right]^{1/2}} + \mathcal{O}(n^{-1}). \quad (3.27a) $$

$$ \frac{\Gamma(nz + \frac{1}{2})\Gamma(\frac{n\tau}{2} - nz)}{\Gamma(nz + 1)\Gamma(\frac{n\tau+1}{2} - nz)} = \frac{1}{n z^{1/2}(\frac{\tau}{2} - z)^{1/2}} + \mathcal{O}(n^{-2}). \quad (3.27b) $$

The first is obtained from (3.25a) and the second from the known asymptotic expansion of $\log \Gamma(z)$. As $n \to \infty$, the sum has a well-defined limit in terms of an integral. Putting all the results together yields

$$ C_{\phi}(x)|_{x=n\tau+1} = n^{1/4} \sin(\pi\tau)^{1/4} (2\pi)^{3/4} G^2(\frac{1}{2}) \int_0^{\tau/2} dz \frac{\cos\left(\phi\left(2z - \frac{\tau}{2}\right)\right)}{\sin(2\pi z)^{1/2} \sin\left(\pi(\tau - 2z)\right)^{1/2}} + \mathcal{O}(n^{-3/4}). \quad (3.28) $$

This is our final formula for the asymptotics for generic values of $\tau$ in $(0, 1)$. The results are displayed in Figure 2.

We analyse its behaviour in the limiting cases $\tau \to 0^+$ and $1^-$. The former is easily evaluated to

$$ C_{\phi}(x)|_{x=n\tau+1} \xrightarrow{n \gg 1, \tau \to 0^+} \left(\frac{n\tau}{2}\right)^{1/4} \pi G^2(\frac{1}{2}). \quad (3.29) $$
The latter requires closer scrutiny, as the integral in (3.28) diverges for $\tau = 1$. Let us denote this integral by $\mathcal{J}$. To obtain its asymptotics for $\tau \to 1^-$, we first use the reflection symmetry of the integrand under $z \to \frac{\tau}{2} - z$:

$$ \mathcal{J} = \int_{0}^{\tau/2} dz \frac{\cos\left(\phi\left(2z - \frac{\tau}{2}\right)\right)}{\sin\left(2\pi z\right)^{1/2} \sin\left(\pi\left(\tau - 2z\right)\right)^{1/2}} = 2 \int_{0}^{\tau/4} dz \frac{\cos\left(\phi\left(2z - \frac{\tau}{2}\right)\right)}{\sin\left(2\pi z\right)^{1/2} \sin\left(2\pi\left(z + \frac{1-\tau}{2}\right)\right)^{1/2}} \quad (3.30) $$

We subtract and add a counter-term in the form of an integral wherein the numerator and denominator are simpler:

$$ \mathcal{J} = 2 \int_{0}^{\tau/4} dz \left[ \frac{\cos\left(\phi\left(2z - \frac{\tau}{2}\right)\right)}{\sin\left(2\pi z\right)^{1/2} \sin\left(2\pi\left(z + \frac{1-\tau}{2}\right)\right)^{1/2}} - \frac{\cos\frac{\phi}{2}}{2\pi z^{1/2}\left(z + \frac{1-\tau}{2}\right)^{1/2}} \right] \\ + \int_{0}^{\tau/4} dz \frac{\cos\frac{\phi}{2}}{\pi z^{1/2}\left(z + \frac{1-\tau}{2}\right)^{1/2}} \quad (3.31) $$

The first integral has a well-defined limit for $\tau \to 1$, whereas the second one is evaluated explicitly:

$$ \mathcal{J} \approx 2 \int_{0}^{1/4} dz \left[ \frac{\cos\left(\phi\left(2z - \frac{1}{2}\right)\right)}{\sin\left(2\pi z\right)} - \frac{\cos\frac{\phi}{2}}{2\pi z} \right] + \frac{2\cos\left(\frac{\phi}{2}\right)}{\pi} \log\left(\frac{\sqrt{1-\frac{\tau}{2}} + \sqrt{\frac{\tau}{2}}}{\sqrt{1-\tau}}\right), \quad (3.32) $$

where $\approx$ indicates an equality up to terms of order $(1 - \tau)^1$. The logarithmic divergence as $\tau \to 1^-$ is explicit in the last term. This yields

$$ C_{\phi}(x) \Big|_{x=n\tau+1} \xrightarrow{n \gg 1, \tau \to 1^{-}} n^{1/4}(1-\tau)^{1/4} 2^{-1/4} G^2(\frac{1}{2}) (K - \alpha \log(\frac{1-\tau}{2})), \quad (3.33) $$

where

$$ K = 2\pi \int_{0}^{1/2} dy \left( \frac{\cos\left(\phi\left(y - \frac{1}{2}\right)\right)}{\sin\pi y} - \frac{\cos\left(\frac{\phi}{2}\right)}{\pi y} \right). \quad (3.34) $$

Recalling from (3.9) that $C_{\phi}(x) = C_{\phi}(n + 2 - x)$, the final results for the asymptotics of these correlation functions for $1 \ll x \ll n$ are

$$ C_{\phi}(x) \xrightarrow{1 \ll x \ll n} \left(\frac{x}{2}\right)^{1/4} \pi G^2\left(\frac{1}{2}\right), \qquad C_{\phi}(x) \xrightarrow{1 \ll x \ll n} \left(\frac{x}{2}\right)^{1/4} G^2\left(\frac{1}{2}\right) \left(\alpha \log n + K - \alpha \log\left(\frac{x}{2}\right)\right). \quad (3.35) $$

The first has a pure power-law behaviour, namely it is proportional to $x^{-2\Delta}$ with $\Delta = -\frac{1}{8}$. The second has the same power-law exponent as the first, but with an added logarithmic dependence upon the position $x$. The conformal interpretation of these results is discussed in Section 5.

# 4 Lattice correlators for double entry points

## 4.1 Refined partition functions

For each loop configuration that contributes to $Z_{\Pi\Pi\Pi}(x)$, the four defects can connect in six possible ways, according to the six link states of $W_{4,0}$, see (2.9). We thus define six refined partition functions, one for each connectivity of the defects:

$$ Z_{\Pi\Pi\Pi}(x) = Z_{\phi\phi}(x) + Z_{\phi\phi}(x) + Z_{\phi\phi}(x) + Z_{\phi\phi}(x) + Z_{\phi\phi}(x) + Z_{\phi\phi}(x). \quad (4.1) $$
The partition functions in the numerators in (4.5) are expressed in terms of spin-chain overlaps as

$$Z_{\text{CFT}}(x) = 2^{mn/2} \langle v_x^d | T^m | v_2^a \rangle, \quad Z_{\text{FT}}(x) = 2^{mn/2} \langle v_x^f | T^m | v_2^a \rangle. \quad (4.6)$$

## 4.2 Closed-form expressions in the limit $m \to \infty$

In the limit $m \to \infty$, the leading contribution for the overlaps is

$$Z_{\text{CFT}}(x) \approx 2^{mn/2} \Lambda_0^m \langle v_x^d | w_0 \rangle \langle w_0 | v_2^a \rangle, \quad Z_{\text{FT}}(x) \approx 2^{mn/2} \Lambda_0^m \langle v_x^f | w_0 \rangle \langle w_0 | v_2^a \rangle, \quad (4.7)$$

and therefore

$$C_{\text{CFT}}(x) = \frac{\langle v_x^d | w_0 \rangle}{\langle v_2^a | w_0 \rangle}, \quad C_{\text{FT}}(x) = \frac{\langle v_x^f | w_0 \rangle}{\langle v_2^a | w_0 \rangle}. \quad (4.8)$$

The states $v_x^d$ and $v_x^f$ are represented in the spin chain by

$$\langle v_x^d | = \langle 0 | a_{n-1} a_{n-3} \cdots a_{x+2} a_{x-2} a_{x-4} \cdots a_3 (\omega c_2 + (-1)^{\frac{x-3}{2}} \omega^{-1} c_x)(\omega c_1 + (-1)^{\frac{x-1}{2}} \omega^{-1} c_{x+1}), \quad (4.9a)$$

$$\langle v_x^f | = \langle 0 | a_{n-1} a_{n-3} \cdots a_{x+2} a_{x-2} a_{x-4} \cdots a_3 (\omega c_2 + (-1)^{\frac{x-3}{2}} \omega^{-1} c_x) \times (\omega^{-1} e^{i\phi/2} c_1 + (-1)^{\frac{x-1}{2}} \omega e^{-i\phi/2} c_{x+1}). \quad (4.9b)$$

As a result, we have

$$C_{\text{CFT}}(x) = \omega^2 f_{1,2} + (-1)^{\frac{x-3}{2}} f_{1,x} + (-1)^{\frac{x-3}{2}} f_{2,x+1} + \omega^{-2} f_{x,x+1}, \quad (4.10a)$$

$$C_{\text{FT}}(x) = e^{i\phi/2} f_{1,2} + e^{i\phi/2} \omega^{-2} (-1)^{\frac{x-3}{2}} f_{1,x} + e^{-i\phi/2} \omega^2 (-1)^{\frac{x-3}{2}} f_{2,x+1} + e^{-i\phi/2} f_{x,x+1}, \quad (4.10b)$$

where

$$f_{a,b} = \frac{\det P^{(a,b)}}{\det M^{(1)}}, \quad P_{ij}^{(a,b)} = \begin{cases} \{c_a, \eta_k^\dagger\} & j=1, \\ \{c_b, \eta_k^\dagger\} & j=2, \\ \{a_{2j-3}, \eta_k^\dagger\} & j=3, \dots, \frac{x+1}{2}, \\ \{a_{2j-1}, \eta_k^\dagger\} & j=\frac{x+3}{2}, \dots, \frac{n}{2}, \end{cases} \quad k \in K. \quad (4.11)$$

We use (3.13), remove from the determinant some of the factors that depend on $k$, and find

$$f_{a,b} = \frac{\det \tilde{P}^{(a,b)}}{\det \tilde{M}^{(1)}}, \quad \tilde{P}_{jk}^{(a,b)} = \begin{cases} \frac{e^{-i\theta_k}}{2\cos(\frac{\theta_k}{2} + \frac{\pi}{4})} & j=1, \\ \frac{e^{-ib\theta_k}}{2\cos(\frac{\theta_k}{2} + \frac{\pi}{4})} & j=2, \\ e^{-i\theta_k}(2j-5/2) & j=3, \dots, \frac{x+1}{2}, \\ e^{-i\theta_k}(2j-1/2) & j=\frac{x+3}{2}, \dots, \frac{n}{2}, \end{cases} \quad k \in K. \quad (4.12)$$

All the rows of $P^{(a,b)}$ except for the first two appear in $M^{(1)}$. Using (3.16), we find

$$Q_{j\ell}^{(a,b)} = [\tilde{P}^{(a,b)}(\tilde{M}^{(1)})^{-1}]_{j\ell} = \begin{cases} \frac{1}{n} \sum_{k \in K} \frac{e^{i\theta_k(2\ell - a - 1/2)}}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})} & j=1, \\ \frac{1}{n} \sum_{k \in K} \frac{e^{i\theta_k(2\ell - b - 1/2)}}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})} & j=2, \\ \delta_{j,k+1} & j=3, \dots, \frac{x+1}{2}, \\ \delta_{j,k} & j=\frac{x+3}{2}, \dots, \frac{n}{2}. \end{cases} \quad (4.13)$$
The determinant of $Q^{(a,b)}$ thus reduces to the determinant of a 2 × 2 matrix:

$$ 
\begin{aligned}
f_{a,b} &= (-1)^{\frac{x-3}{2}} \det \begin{pmatrix} Q_{1,1}^{(a,b)} & Q_{1,(x+1)/2}^{(a,b)} \\ Q_{1,2}^{(a,b)} & Q_{2,(x+1)/2}^{(a,b)} \end{pmatrix} \\
&= \frac{(-1)^{\frac{x-3}{2}}}{n^2} \sum_{k,l \in K} e^{-i(a\theta_k+b\theta_l)} \frac{e^{\frac{3}{2}i\theta_k+i\theta_l(x+\frac{1}{2})} - e^{\frac{3}{2}i\theta_l+i\theta_k(x+\frac{1}{2})}}{\cos(\frac{\theta_k}{2}+\frac{\pi}{4})\cos(\frac{\theta_l}{2}+\frac{\pi}{4})}.
\end{aligned} 
\quad (4.14) $$

After simplifications, this yields

$$ C_{\mathcal{G}}(x) = \frac{2}{n^2} \sum_{k,l \in K} \left( \frac{\cos(\frac{\theta_l + \theta_k}{2}) - \cos((x - \frac{1}{2})\theta_l - (x - \frac{3}{2})\theta_k)}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})\cos(\frac{\theta_l}{2} + \frac{\pi}{4})} + (-1)^{\frac{x-1}{2}} \frac{\sin((x - \frac{3}{2})\theta_l + \frac{\theta_k}{2}) - \sin((x - \frac{1}{2})\theta_l - \frac{\theta_k}{2})}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})\cos(\frac{\theta_l}{2} + \frac{\pi}{4})} \right), \quad (4.15a) $$

$$ C_{\mathcal{G}}(x) = \frac{2}{n^2} \sum_{k,l \in K} \left( \frac{\sin(\frac{\theta_l + \theta_k}{2} + \frac{\phi}{2}) - \sin((x - \frac{1}{2})\theta_l - (x - \frac{3}{2})\theta_k + \frac{\phi}{2})}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})\cos(\frac{\theta_l}{2} + \frac{\pi}{4})} + (-1)^{\frac{x-3}{2}} \frac{\cos((x - \frac{3}{2})\theta_l + \frac{\theta_k}{2} + \frac{\phi}{2}) - \cos((x - \frac{1}{2})\theta_l - \frac{\theta_k}{2} + \frac{\phi}{2})}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})\cos(\frac{\theta_l}{2} + \frac{\pi}{4})} \right). \quad (4.15b) $$

Many of these double sums can be reduced to single sums. For instance,

$$ 
\begin{aligned}
\sum_{k,l \in K} \frac{\cos(\frac{\theta_l + \theta_k}{2})}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})\cos(\frac{\theta_l}{2} + \frac{\pi}{4})} &= \sum_{k,l \in K} \frac{\cos(\frac{\theta_l}{2} + \frac{\pi}{4})\sin(\frac{\theta_k}{2} + \frac{\pi}{4}) + \cos(\frac{\theta_k}{2} + \frac{\pi}{4})\sin(\frac{\theta_l}{2} + \frac{\pi}{4})}{\cos(\frac{\theta_k}{2} + \frac{\pi}{4})\cos(\frac{\theta_l}{2} + \frac{\pi}{4})} \\
&= \sum_{k,l \in K} \tan(\frac{\theta_k}{2} + \frac{\pi}{4}) + \tan(\frac{\theta_l}{2} + \frac{\pi}{4}) = n \sum_{k \in K} \tan(\frac{\theta_k}{2} + \frac{\pi}{4}). \quad (4.16)
\end{aligned}
 $$

The other terms are simplified using similar ideas, and after some algebra we find

$$ C_{\mathcal{G}}(x) = \frac{4}{n} \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right) \sin^2\left((x-1)\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right)\right), \quad (4.17a) $$

$$ 
\begin{aligned}
C_{\mathcal{G}}(x) = & \frac{4 \sin(\frac{\phi}{2})}{n} \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right) \sin^2\left((x-1)\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right)\right) \\
& + \frac{2 \cos(\frac{\phi}{2})}{n} \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right) \sin\left(2(x-1)\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right)\right) \\
& - \frac{2 \cos(\frac{\phi}{2})}{n^2} \left[ \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right) \cos\left(2(x-1)\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right)\right) \right]^2 \\
& - \frac{2 \cos(\frac{\phi}{2})}{n^2} \left[ \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right) \sin\left(2(x-1)\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right)\right) \right]^2 \\
& + \frac{2 \cos(\frac{\phi}{2})}{n^2} \left[ \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right) \right]^2 - \frac{\cos(\frac{\phi}{2})}{2}.
\end{aligned} 
\quad (4.17b) $$
In the second expression, some of the double sums could not be simplified to single sums and were instead rewritten as the square of single sums. We simplify these expressions further using the identity

$$ \frac{\sin(\frac{\ell \xi}{n})}{\cos(\frac{\xi}{n})} = \sum_{t=0}^{\ell-1} (-1)^t \sin(\frac{\xi}{n}(\ell - 1 - 2t)), \quad \ell \text{ even.} \qquad (4.18) $$

After some manipulations, we find

$$ C_{\mathcal{D}}(x) = \frac{4}{n} \sum_{t=0}^{(x-3)/2} \frac{\cos\left(\frac{\phi}{n}(2t+1)\right)}{\sin\left(\frac{\pi}{n}(2t+1)\right)}, \qquad (4.19a) $$

$$ C_{\mathcal{D}}(x) = \left( \frac{4\sin(\frac{\phi}{2})}{n} + \frac{8Y(n)\cos(\frac{\phi}{2})}{n^2} \right) \sum_{t=0}^{(x-3)/2} \frac{\cos\left(\frac{\phi}{n}(2t+1)\right)}{\sin\left(\frac{\pi}{n}(2t+1)\right)} \\ - \frac{8\cos(\frac{\phi}{2})}{n^2} \left( \left[ \sum_{t=0}^{(x-3)/2} \frac{\cos\left(\frac{\phi}{n}(2t+1)\right)}{\sin\left(\frac{\pi}{n}(2t+1)\right)} \right]^2 + \left[ \sum_{t=0}^{(x-3)/2} \frac{\sin\left(\frac{\phi}{n}(2t+1)\right)}{\sin\left(\frac{\pi}{n}(2t+1)\right)} \right]^2 \right), \qquad (4.19b) $$

where

$$ Y(n) = \sum_{k=0}^{(n-2)/2} \cot\left(\frac{\pi}{n}\left(k + \frac{1}{2}\right) + \frac{\phi}{2n}\right). \qquad (4.20) $$

## 4.3 Asymptotic behaviour

We compute the asymptotics for $\tau = \frac{x-1}{n}$ with $\tau \in (0, 1)$ and $n \to \infty$. We find

$$ \frac{1}{n} \sum_{t=0}^{(x-3)/2} \frac{\cos\left(\frac{\phi}{n}(2t+1)\right)}{\sin\left(\frac{\pi}{n}(2t+1)\right)} = \frac{1}{2\pi}\left(\log(n\tau) + \log 2 + \gamma\right) + \mathcal{I}_1(\tau) + \mathcal{O}(n^{-1}), \qquad (4.21a) $$

$$ \frac{1}{n} \sum_{t=0}^{(x-3)/2} \frac{\sin\left(\frac{\phi}{n}(2t+1)\right)}{\sin\left(\frac{\pi}{n}(2t+1)\right)} = \mathcal{I}_2(\tau) + \mathcal{O}(n^{-1}), \qquad (4.21b) $$

$$ \frac{Y(n)}{n} = \frac{1}{\pi} \left( \log n - \log \pi - \psi\left(\frac{1}{2} + \frac{\phi}{2\pi}\right) + \mathcal{O}(n^{-1}) \right), \qquad (4.21c) $$

where $\psi(z) = \Gamma'(z)/\Gamma(z)$ is the digamma function, $\gamma$ is the Euler-Mascheroni constant, and

$$ \mathcal{I}_1(\tau) = \frac{1}{2} \int_0^\tau dt \left( \frac{\cos(\phi t)}{\sin(\pi t)} - \frac{1}{\pi t} \right), \quad \mathcal{I}_2(\tau) = \frac{1}{2} \int_0^\tau dt \left( \frac{\sin(\phi t)}{\sin(\pi t)} \right). \qquad (4.22) $$

From (4.19a), we obtain

$$ C_{\mathcal{D}}(x)|_{x=n\tau+1} = \frac{2}{\pi} \left( \log(n\tau) + \log 2 + \gamma + 2\pi \mathcal{I}_1(\tau) + \mathcal{O}(n^{-1}) \right). \qquad (4.23) $$

The functions $\mathcal{I}_1(\tau)$ and $\mathcal{I}_2(\tau)$ have the following asymptotic behaviour for $\tau \to 0^+$:

$$ \mathcal{I}_1(\tau) = \frac{\pi \tau^2}{8} \left[ \frac{1}{3} - \left(\frac{\phi}{\pi}\right)^2 \right] + \mathcal{O}(\tau^4), \quad \mathcal{I}_2(\tau) = \frac{\phi \tau}{2\pi} + \mathcal{O}(\tau^3). \qquad (4.24) $$
For $\tau \to 1^-$, we have

$$ \mathcal{I}_1(\tau) = -\frac{\cos \phi}{2\pi} \log(1-\tau) + \cos^2(\frac{\phi}{2}) \int_0^{1/2} dt \left( \frac{\cos(\phi t)}{\sin(\pi t)} - \frac{1}{\pi t} \right) + \frac{\sin \phi}{2} \int_0^{1/2} dt \frac{\sin(\phi t)}{\sin(\pi t)} \\ - \frac{\cos^2(\frac{\phi}{2}) \log 2}{\pi} + \frac{1}{2\pi} (1-\tau)(1-\phi \sin \phi) + \mathcal{O}((1-\tau)^2). \qquad (4.25) $$

The asymptotic behaviour for $C_{\mathcal{I}}(x)$ is then given by

$$ C_{\mathcal{I}}(x) \Big|_{x=n\tau+1} \xrightarrow{n \gg 1, \tau \to 0^+} \frac{2}{\pi} (\log(n\tau) + \log 2 + \gamma), \qquad (4.26a) $$

$$ C_{\mathcal{I}}(x) \Big|_{x=n\tau+1} \xrightarrow{n \gg 1, \tau \to 1^-} \frac{2}{\pi} \left( \log n - \cos \phi \log(1-\tau) - \cos \phi \log 2 + \gamma \right) + \hat{K}, \qquad (4.26b) $$

where

$$ \hat{K} = (2 \cos(\frac{\phi}{2}))^2 \int_0^{1/2} dt \left( \frac{\cos(\phi t)}{\sin(\pi t)} - \frac{1}{\pi t} \right) + 2 \sin \phi \int_0^{1/2} dt \frac{\sin(\phi t)}{\sin(\pi t)}. \qquad (4.27) $$

Likewise from (4.19a), we obtain

$$ C_{\mathcal{I}}(x) \Big|_{x=n\tau+1} = \frac{2}{\pi^2} \left[ \cos(\frac{\phi}{2})(\log n)^2 - \cos(\frac{\phi}{2})\left( (2\pi \mathcal{I}_2(\tau))^2 + (\gamma + 2\pi \mathcal{I}_1(\tau) + \log 2 + \log \tau)^2 \right) \right. \\ \left. + (\gamma + 2\pi \mathcal{I}_1(\tau) + \log 2 + \log(n\tau))(\pi \sin(\frac{\phi}{2}) - 2 \cos(\frac{\phi}{2})(\log \pi + \psi(\frac{1}{2} + \frac{\phi}{2\pi}))) \right] + \mathcal{O}(n^{-1}). \qquad (4.28) $$

In the limit $\tau \to 0^+$, this yields

$$ C_{\mathcal{I}}(x) \Big|_{x=n\tau+1} \xrightarrow{n \gg 1, \tau \to 0^+} \frac{2}{\pi^2} \left[ \cos(\frac{\phi}{2})(\log n)^2 - \cos(\frac{\phi}{2})\left( (\gamma + \log 2 + \log \tau)^2 \right) \right. \\ \left. + (\gamma + \log 2 + \log(n\tau))(\pi \sin(\frac{\phi}{2}) - 2 \cos(\frac{\phi}{2})(\log \pi + \psi(\frac{1}{2} + \frac{\phi}{2\pi}))) \right]. \qquad (4.29) $$

The results are plotted in the right panel of Figure 2. As expected from (4.4b), $C_{\mathcal{I}}(x)$ is symmetric under the transformation $\tau \to 1-\tau$. This is not immediately obvious from (4.28), but can be shown using the identities

$$ \mathcal{I}_1(1-\tau) = -\frac{\cos \phi}{2\pi} \left( \log \tau - \log 2 + 2\pi \mathcal{I}_1(\tau) \right) - \mathcal{I}_2(\tau) \sin \phi - \frac{\log(\frac{1-\tau}{2})}{2\pi} \\ - \frac{\cos^2(\frac{\phi}{2})}{\pi} \left( \gamma + \log \pi + 2 \log 2 + \psi(\frac{1}{2} + \frac{\phi}{2\pi}) \right) + \frac{\sin \phi}{4}, \qquad (4.30a) $$

$$ \mathcal{I}_2(1-\tau) = -\frac{\sin \phi}{2\pi} \left( \gamma + \log 2 + \log \pi + \log \tau + \psi(\frac{1}{2} + \frac{\phi}{2\pi}) + 2\pi \mathcal{I}_1(\tau) \right) + \mathcal{I}_2(\tau) \cos \phi + \frac{1}{2} \sin^2(\frac{\phi}{2}). \qquad (4.30b) $$

Recalling from (4.4) that $C_{\mathcal{I}}(x) = C_{\mathcal{I}}(n+2-x)$, the final results for the asymptotics for
<figure><img src="image_6.png" /><figcaption>Figure 2: The correlation functions $C_{\mathcal{G}}(x)$ and $C_{\mathcal{G}^c}(x)$ for $\phi = 1$. The points are the exact values for $n = 200$. In the left panel, the blue, orange and purple solid curves are drawn from (4.23), (4.26a) and (4.26b). In the right panel, the blue and orange curves are drawn from (4.28) and (4.29).</figcaption></figure>

1 ≪ x ≪ n are

$$C_{\mathcal{G}}(x) \xrightarrow{1 \ll x \ll n} \frac{2}{\pi}(\log x + \log 2 + \gamma), \qquad (4.31a)$$

$$C_{\mathcal{G}^c}(x) \xrightarrow{1 \ll x \ll n} \frac{2}{\pi}(2 \cos^2(\frac{\phi}{2}) \log n - \cos \phi \log x) + \hat{K}, \qquad (4.31b)$$

$$C_{\mathcal{G}^c}(x) \xrightarrow{1 \ll x \ll n} \frac{2}{\pi^2}(\gamma + \log 2 + \log x)(\pi \sin(\frac{\phi}{2}) - 2 \cos(\frac{\phi}{2})(\log n + \log \pi + \psi(\frac{1}{2} + \frac{\phi}{2\pi}))) \qquad (4.31c)$$

$$-\frac{2 \cos(\frac{\phi}{2})}{\pi^2}(\gamma + \log 2 + \log x)^2.$$

These lattice correlators have a purely logarithmic behavior. The power-law behavior is thus absent, consistent with the conformal weight $\Delta = 0$. We discuss the conformal interpretation of these results in Section 6.

# 5 Conformal correlators for single entry points

## 5.1 Preliminaries

In this section, we provide an interpretation of the lattice results of Section 3 in terms of logarithmic conformal field theory. We denote by $\mathbb{V}$ the semi-infinite cylinder drawn in the plane. We claim that, in the conformal description, the lattice correlation function $C_{\mathcal{G}}(x)$ is a
ratio of correlation functions on V:

$$C_{\phi}(x_{21}) = \lim_{x_3 \to i\infty} \frac{\langle \phi_{\phi}(x_1) \phi_{\phi}(x_2) \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathcal{V}}}{\langle \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathcal{V}}}, \qquad (5.1)$$

where $x_{21} = x_2 - x_1$.

The fields $\phi_{\phi}(x)$ and $\phi_{\psi}(x)$ are chiral primary fields of conformal weight $-\frac{1}{8}$ that live on the boundary. As their labels suggest, their insertion on a boundary respectively marks the start or end of a long boundary arc. In our derivation below, we assume that these fields are highest weight states $|\phi\rangle$ in irreducible representations of the Virasoro algebra with $c = -2$ and $\Delta = \Delta_{1,2} = -\frac{1}{8}$. By contrast, $\psi_{\alpha}(z, z^*)$ is a field that lives in the bulk and is therefore not chiral. It changes the weight of the loops encircling the point $z$ from 0 to $\alpha$. Its action is therefore non-local. The derivation below uses the method of images [46] and works with the assumption that correlators involving $\psi_{\alpha}(z, z^*)$ on the upper-half plane are equal to correlators on the full plane with this field replaced by $\psi_{\alpha}(z)\psi_{\alpha}(z^*)$, where $\psi_{\alpha}(z)$ is a primary field. The calculations below strongly support this claim and will allow us to extract the value of the conformal weight $\Delta$ of $\psi_{\alpha}(z)$ as a function of $\alpha$, see (5.36). The transformation laws of the fields under conformal maps are

$$\phi(z) \mapsto \phi(y) = \left(\frac{dy}{dz}\right)^{1/8} \phi(z), \qquad \psi_{\alpha}(z) \mapsto \psi_{\alpha}(y) = \left(\frac{dy}{dz}\right)^{-\Delta} \psi_{\alpha}(z), \qquad (5.2)$$

and the two-point functions are

$$\langle \phi(z_1) \phi(z_2) \rangle = \tilde{\kappa}^{\phi\phi} (z_2 - z_1)^{1/4}, \qquad \langle \psi_{\alpha}(z_1) \psi_{\alpha}(z_2) \rangle = \frac{\kappa^{\psi\psi}}{(z_2 - z_1)^{2\Delta}}, \qquad (5.3)$$

where $\tilde{\kappa}^{\phi\phi}$ and $\kappa^{\psi\psi}$ are constants. The fields $\phi_{\phi}(x)$ and $\phi_{\psi}(x)$ transform identically under conformal transformations. As we shall see in Section 7.1, the distinction between these two fields lies in their fusion rules and in the values of the structure constants. For convenience, their labels are omitted in this section.

The conformal map from V to the upper-half plane H is

$$z(x) = \frac{\sin\left(\frac{\pi}{n}(x - x_1)\right)}{\sin\left(\frac{\pi}{n}(x_2 - x_1)\right)} \frac{\sin\left(\frac{\pi x_2}{n}\right)}{\sin\left(\frac{\pi x}{n}\right)}. \qquad (5.4)$$

It is illustrated in Figure 3. The points $x_1$ and $x_2$ are sent to 0 and 1, whereas the point at $x = i\infty$ is mapped to a finite point in the upper-half plane:

$$\lim_{x \to i\infty} z(x) = e^{i\pi x_1 / n} \frac{\sin\left(\frac{\pi x_2}{n}\right)}{\sin\left(\frac{\pi}{n}(x_2 - x_1)\right)}. \qquad (5.5)$$

The derivative of this transformation is

$$\frac{dz(x)}{dx} = \frac{\pi}{n} \frac{\sin\left(\frac{\pi x_1}{n}\right) \sin\left(\frac{\pi x_2}{n}\right)}{\sin\left(\frac{\pi}{n}(x_2 - x_1)\right) \sin\left(\frac{\pi x}{n}\right)^2}. \qquad (5.6)$$

## 5.2 Differential equation for the four-point function

We consider the correlation functions appearing in (5.1), but defined on H. Using the method of images, the correlation functions on H are equal to correlation functions in the full complex plane:

$$\langle \phi(z_1) \phi(z_2) \psi_{\alpha}(z, z^*) \rangle_{\mathcal{H}} = \langle \phi(z_1) \phi(z_2) \psi_{\alpha}(z) \psi_{\alpha}(z^*) \rangle_{\mathcal{C}}, \qquad \langle \psi_{\alpha}(z, z^*) \rangle_{\mathcal{H}} = \langle \psi_{\alpha}(z) \psi_{\alpha}(z^*) \rangle_{\mathcal{C}}. \qquad (5.7)$$
<figure><img src="image_7.png" /><figcaption>Figure 3: The map $z(x)$ from the semi-infinite cylinder $V$ to the upper half-plane $H$.</figcaption></figure>

We proceed to compute the four-point function

$$G = \langle \phi(z_1)\phi(z_2)\psi_\alpha(z_3)\psi_\alpha(z_4) \rangle_C. \qquad (5.8)$$

We use the known techniques of conformal field theory to compute four-point functions of primary fields in the case where one of them has a null descendant [3, 4]. The conformal Ward identities are

$$\left(\sum_{i=1}^{4} \frac{\partial}{\partial z_i}\right) G = 0, \quad \left(\sum_{i=1}^{4} z_i \frac{\partial}{\partial z_i} + \Delta_i\right) G = 0, \quad \left(\sum_{i=1}^{4} z_i^2 \frac{\partial}{\partial z_i} + 2\Delta_i z_i\right) G = 0, \qquad (5.9)$$

and imply that $G$ is of the form

$$G = \frac{(z_2 - z_1)^{1/4}}{(z_4 - z_3)^{2\Delta}} \tilde{G}(\eta), \quad \eta = \frac{(z_1 - z_2)(z_3 - z_4)}{(z_1 - z_3)(z_2 - z_4)}, \qquad (5.10)$$

where $\eta$ is the cross-ratio. Because $|\phi\rangle$ has the null descendant $(2L_{-1}^2 - L_{-2})|\phi\rangle$ at level two, $G$ satisfies the second-order partial differential equation

$$\left[ 2 \frac{\partial^2}{\partial z_1^2} - \sum_{i=2}^{4} \left( \frac{1}{z_1 - z_i} \frac{\partial}{\partial z_i} + \frac{\Delta_i}{(z_1 - z_i)^2} \right) \right] G = 0. \qquad (5.11)$$

This differential equation is rewritten as an ordinary differential equation for $\tilde{G}(\eta)$:

$$\eta(1 - \eta)\tilde{G}''(\eta) + (1 - \frac{3\eta}{2})\tilde{G}'(\eta) - \frac{\Delta}{2} \frac{\eta}{1 - \eta} \tilde{G}(\eta) = 0. \qquad (5.12)$$

By setting

$$\tilde{G}(\eta) = (1 - \eta)^{\rho/2} F(\eta), \quad \rho = \frac{1}{2}(1 - \sqrt{1 + 8\Delta}), \qquad (5.13)$$

we find that $F(\eta)$ satisfies the hypergeometric differential equation

$$z(1-z)F''(z) + (c - (1 + a + b)z)F'(z) - abF(z) = 0 \qquad (5.14)$$

for

$$a = \rho, \quad b = \frac{1}{2}, \quad c = 1. \qquad (5.15)$$

For c generic, the two solutions are $_2F_1(a, b; c|z)$ and $z^{1-c} {}_2F_1(1+a-c, 1+b-c; 2-c|z)$. For c = 1, these two solutions coincide and a second independent solution [61] is given by

$$H(a, b|z) = {}_2F_1(a, b; 1|z) \log z + \sum_{k=0}^{\infty} \frac{(a)_k (b)_k z^k}{(k!)^2} (\psi(a+k) + \psi(b+k) - 2\psi(1+k)), \qquad (5.16)$$
where $\psi(x)$ is the digamma function. The general solution is a linear combination of these two solutions.

Before continuing, let us immediately use the transformation laws (5.2) of the fields to obtain the form of the lattice correlator on the cylinder. This computation is achieved by first considering the insertion points of the fields $\psi_\alpha(x_3, x_4)$ on $\mathbb{V}$ to be in finite positions $x_3$ and $x_4$, instead of $\pm i\infty$. This yields

$$ \langle \phi(x_1)\phi(x_2)\psi_\alpha(x_3, x_4) \rangle_{\mathbb{V}} = \left(\frac{n}{\pi}\sin\left(\frac{\pi}{n}(x_2 - x_1)\right)\right)^{1/4} \left(\frac{n}{\pi}\sin\left(\frac{\pi}{n}(x_4 - x_3)\right)\right)^{-2\Delta} (1-\eta)^{\rho/2} F(\eta), \quad (5.17a) $$

$$ \langle \psi_\alpha(x_3, x_4) \rangle_{\mathbb{V}} = \left( \frac{n}{\pi} \sin\left(\frac{\pi}{n}(x_4 - x_3)\right) \right)^{-2\Delta} \kappa^{\psi\psi}. \quad (5.17b) $$

For $x_4 = x_3^*$ and $x_3 \to i\infty$, the expression $\sin(\frac{\pi}{n}(x_4 - x_3))$ diverges. However, the ratio of these correlation functions, as in (5.1), has a well-defined limit. In the same $x_3, x_4$ limit, the cross-ratio is given by

$$ \eta = 1 - e^{2\pi i \tau}, \qquad \tau = \frac{x_2 - x_1}{n}. \quad (5.18) $$

The parameter $\tau$ is the same as the one that was defined in Section 3.3. It lies in the range $(0, 1)$. As $\tau$ explores this range, $\eta$ takes complex values and draws a counterclockwise circle of unit radius around the point $\eta = 1$, starting and ending at the origin. We therefore have

$$ C_{\mathbb{V}}(x_{21}) = \frac{1}{\kappa^{\psi\psi}} \tilde{G}(1 - e^{2\pi i \tau}). \quad (5.19) $$

The hypergeometric function ${}_2F_1(a, b; c|\eta)$ has a branch cut on the real $\eta$-axis for $\eta > 1$, which the circle crosses at $\eta = 2$, for $\tau = \frac{1}{2}$. One way to obtain a smooth function for $C_{\mathbb{V}}(x_{21})$ is to search for two solutions to the differential equation, one for $\tau \in (0, \frac{1}{2})$ and another for $\tau \in (\frac{1}{2}, 1)$, and to impose that the function be smooth at $\tau = \frac{1}{2}$.

Before proceeding with this plan, we note that, instead of (5.13), one can choose to set

$$ \tilde{G}(\eta) = \tilde{F}\left(\frac{\eta}{4(\eta - 1)}\right). \quad (5.20) $$

Then the function $\tilde{F}(z)$ satisfies (5.14) with

$$ a = \frac{\rho}{2}, \qquad b = \frac{1-\rho}{2}, \qquad c = 1. \quad (5.21) $$

The general solution is again given by a linear combination of ${}_2F_1(a, b; 1|z)$ and $H(a, b|z)$ for the appropriate values of $a$ and $b$. This basis of solutions is more convenient because, for $z_1 = 0, z_2 = 1, z_3 = z(i\infty)$ and $z_4 = z(-i\infty)$, the expression

$$ \frac{\eta}{4(\eta - 1)} = \sin^2(\pi\tau) \quad (5.22) $$

is real for $\tau \in (0, 1)$. In this basis, the argument $z$ of the two solutions remains real and bounded in $(0, 1]$. It makes contact with the endpoint of the branch cut at $\tau = \frac{1}{2}$. The derivative of ${}_2F_1(\frac{\rho}{2}, \frac{1-\rho}{2}; 1|\sin^2(\pi\tau))$ with respect to $\tau$ is discontinuous at $\tau = \frac{1}{2}$, and likewise for $H(\frac{\rho}{2}, \frac{1-\rho}{2}|\sin^2(\pi\tau))$, so we search for separate solutions on the two sub-intervals. We write the general solution as

$$ C_{\mathbb{V}}(x_{21}) = (n \sin(\pi \tau))^{1/4} \times \begin{cases} A_1 {}_2F_1(\frac{\rho}{2}, \frac{1-\rho}{2}; 1|\sin^2(\pi \tau)) + A_2 H(\frac{\rho}{2}, \frac{1-\rho}{2}|\sin^2(\pi \tau)) & \tau < \frac{1}{2}, \\ B_1 {}_2F_1(\frac{\rho}{2}, \frac{1-\rho}{2}; 1|\sin^2(\pi \tau)) + B_2 H(\frac{\rho}{2}, \frac{1-\rho}{2}|\sin^2(\pi \tau)) & \tau > \frac{1}{2}, \end{cases} \quad (5.23) $$

where $A_1, A_2, B_1, B_2$ are constants to be determined.
## 5.3 Solving for the unknowns and verifying smoothness

Near $z = 0$, the functions ${}_2F_1(a, b; 1|z)$ and $H(a, b|z)$ behave as

$${}_2F_1(a, b; 1|z) = 1 + \mathcal{O}(z), \quad H(a, b|z) = \log(z) + (\psi(a) + \psi(b) + 2\gamma) + \mathcal{O}(z). \quad (5.24)$$

We fix $A_1, A_2, B_1, B_2$ by imposing that the expression (5.23) for $C_\phi(x_{21})$ has the same asymptotic behaviour, for $\tau \to 0^+$ and $\tau \to 1^-$, as the result we obtained from the lattice, see (3.29) and (3.33). We find

$$A_1 = \frac{\pi^{3/4}}{2^{1/4}} G^2(\frac{1}{2}), \quad A_2 = 0, \quad B_2 = -\frac{\alpha}{2} \frac{G^2(\frac{1}{2})}{(2\pi)^{1/4}}, \quad (5.25a)$$

$$B_1 = \frac{G^2(\frac{1}{2})}{(2\pi)^{1/4}} \left[ K + \alpha \left( \log 2\pi + \gamma + \frac{1}{2} \left( \psi(\frac{\rho}{2}) + \psi(\frac{1-\rho}{2}) \right) \right) \right]. \quad (5.25b)$$

The only remaining unknown is the conformal weight $\Delta$. It is fixed by imposing that the solution is continuous at $\tau = \frac{1}{2}$:

$$(A_1 - B_1) {}_2F_1(\frac{\rho}{2}, \frac{1-\rho}{2}; 1|1) = B_2 H(\frac{\rho}{2}, \frac{1-\rho}{2}|1). \quad (5.26)$$

The left side is evaluated using the relation

$${}_2F_1(a, b; c|1) = \frac{\Gamma(c)\Gamma(c-a-b)}{\Gamma(c-a)\Gamma(c-b)}, \quad \text{Re}(c-a-b) > 0. \quad (5.27)$$

To evaluate the right side, we use the identity

$$H(a, b|z) = \lim_{c \to 1} \left( \frac{\partial}{\partial a} + \frac{\partial}{\partial b} + 2 \frac{\partial}{\partial c} + \psi(a) + \psi(b) - 2\psi(c) \right) {}_2F_1(a, b; c|z). \quad (5.28)$$

Combining this with (5.27) yields

$$H(a, b|1) = \frac{\Gamma(1-a-b)}{\Gamma(1-a)\Gamma(1-b)} (\psi(a) + \psi(b) - \psi(1-a) - \psi(1-b)). \quad (5.29)$$

The equality (5.26) then simplifies to

$$K = \pi - \alpha \left( \log 2\pi + \gamma + \frac{1}{2} \left( \psi\left(\frac{2-\rho}{2}\right) + \psi\left(\frac{1+\rho}{2}\right) \right) \right). \quad (5.30)$$

We compare this with the integral expression (3.34) for $K$, which we now rewrite in terms of digamma functions:

$$\begin{aligned} K &= 2\pi i \int_0^{1/2} dy \frac{\cos\left(\phi\left(y - \frac{1}{2}\right)\right) - \cos\left(\frac{\phi}{2}\right)}{\sin\pi y} + 2\pi \cos\left(\frac{\phi}{2}\right) \int_0^{1/2} dy \left(\frac{1}{\sin(\pi y)} - \frac{1}{\pi y}\right) \quad (5.31) \\ &= e^{i\phi/2} L(\phi) + e^{-i\phi/2} L(-\phi) + \alpha(2\log 2 - \log \pi), \end{aligned}$$

where

$$\begin{aligned} L(\phi) &= 2\pi i \int_0^{1/2} dy \frac{e^{-i\phi y} - 1}{e^{i\pi y} - e^{-i\pi y}} = - \int_{e^{i0}}^{e^{-i\pi}} \frac{dt}{t^{1/2}} \frac{t^{\frac{\phi}{2\pi}} - 1}{1 - t} = \int_{-1}^{1} \frac{dt}{t^{1/2}} \frac{t^{\frac{\phi}{2\pi}} - 1}{1 - t} \\ &= \int_0^1 \frac{dt}{t^{1/2}} \frac{t^{\frac{\phi}{2\pi}} - 1}{1 - t} + \int_0^1 \frac{dt}{t^{1/2} e^{-i\pi/2}} \frac{t^{\frac{\phi}{2\pi}} e^{-i\frac{\phi}{2}} - 1}{1 + t} \\ &= \psi(\frac{1}{2}) - \psi(\frac{\phi}{2\pi} + \frac{1}{2}) - \frac{i\pi}{2} + i e^{-i\phi/2} (\psi(\frac{\phi}{4\pi} + \frac{3}{4}) - \psi(\frac{\phi}{4\pi} + \frac{1}{4})). \end{aligned} \quad (5.32)$$
At the last step, we used the following integral formula for the digamma function:

$$ \psi(z) = -\gamma + \int_{0}^{1} dt \frac{1 - t^{z-1}}{1 - t}, \qquad \mathrm{Re}(z) > 0. \qquad (5.33) $$

The corresponding expression for K is simplified using properties of $\psi(z)$:

$$ \psi(1-z) - \psi(z) = \pi \cot \pi z, \qquad \psi(2z) = \log 2 + \frac{1}{2} \left( \psi(z) + \psi(z + \frac{1}{2}) \right), \qquad \psi(\frac{1}{2}) = -2 \log 2 - \gamma. \qquad (5.34) $$

We find

$$ K = \pi - \alpha \left( \log 2\pi + \gamma + \frac{1}{2} \left( \psi\left(\frac{3}{4} + \frac{\phi}{4\pi}\right) + \psi\left(\frac{3}{4} - \frac{\phi}{4\pi}\right) \right) \right). \qquad (5.35) $$

Comparing with (5.30), we see that the equality is satisfied for $\rho = \frac{1}{2} \pm \frac{\phi}{2\pi}$, both of which lead to the following value for the conformal weight:

$$ \Delta = \frac{(\frac{\phi}{\pi})^2 - 1}{8}. \qquad (5.36) $$

This is precisely the value expected from Coulomb gas arguments. Indeed, the groundstate of the transfer matrix of the six-vertex model with periodic twisted boundary condition in the zero-magnetization sector has precisely this conformal dimension [62, 63]. We also note that, with these values of $\rho$, the expression for $B_1$ in (5.25b) simplifies to

$$ B_1 = -\frac{\pi^{3/4}}{2^{1/4}} G^2(\frac{1}{2}) = -A_1. \qquad (5.37) $$

All the unknowns have been solved for. It remains to check that the expression (5.23) is a smooth function of $\tau$ at $\tau = \frac{1}{2}$, namely

$$ \lim_{\tau \to (1/2)^-} \frac{d^k C_{\phi}(x_{21})}{d\tau^k} = \lim_{\tau \to (1/2)^+} \frac{d^k C_{\phi}(x_{21})}{d\tau^k}, \qquad k = 1, 2, \dots \qquad (5.38) $$

It in fact suffices to check that the equality holds for $k = 1$. Indeed, the function $C_{\phi}(x_{21})$ satisfies the same second-order differential equation in each of the two sub-intervals. The higher-order derivatives at $\tau = \frac{1}{2}$ are then uniquely fixed by $C_{\phi}(x_{21})$ and $\frac{d}{d\tau}C_{\phi}(x_{21})$ at this point. The continuity of $C_{\phi}(x_{21})$ and of its first derivative at $\tau = \frac{1}{2}$ therefore implies that all the higher-order derivatives are continuous and that the function is smooth at this point.

The equality (5.38) for $k = 1$ boils down to

$$ \lim_{\tau \to 1/2} \frac{d}{d\tau} \sum_{k=0}^{\infty} \frac{(a)_k (b)_k \sin^{2k}(\pi\tau)}{(k!)^2} (\psi(a+k) + \psi(b+k) - 2\psi(1+k)) = 0. \qquad (5.39) $$

This identity is readily verified by applying the derivative to each term in the sum, rewriting the result in terms of a differential operator acting on a hypergeometric function similarly to (5.28), and finally evaluating the limit $\tau \to \frac{1}{2}$ using (5.27).

We have therefore produced a prediction for $C_{\phi}(x_{21})$ using conformal invariance. By plotting the corresponding curve alongside the exact data for $n = 200$, we find that the expression (5.23) with the constants fixed as in (5.25a) and (5.37) precisely reproduces the solid curve in Figure 2. It is indeed non-trivial that the solution obtained from conformal invariance, which is defined separately on the two sub-intervals, is equal to the integral expression (3.28), which is defined with a unique expression on the full interval.
To prove that the two expressions are equal, we check that the function in (3.28) satisfies the differential equation predicted by conformal field theory. Rewriting (5.12) as a second-order differential equation in $\tau$ and setting $\Delta$ to its value (5.36), we find that the following identity should hold:

$$ \left( \frac{d^2}{d\tau^2} + \pi \cot(\pi\tau) \frac{d}{d\tau} + 2\pi^2\Delta \right) \mathcal{J} = 0, \qquad (5.40) $$

where $\mathcal{J}$ is defined in (3.30). With a simple rescaling of the integration variable, it is expressed as

$$ \mathcal{J} = \frac{1}{2} \int_{0}^{1} dz \frac{\tau \cos\left(\frac{\phi}{2}(1-2z)\right)}{\sin(\pi\tau z)^{1/2} \sin(\pi\tau(1-z))^{1/2}}. \qquad (5.41) $$

The differential operator in (5.40) can then be applied to the integrand. After a bit of work, we find that (5.40) translates to

$$ \frac{1}{4} \int_{0}^{1} dz \frac{d}{dz} \left[ \frac{\pi \cos\left(\frac{\phi}{2}(1-2z)\right) \sin\left(\pi\tau(1-z)\right)^{1/2} \sin(\pi\tau z)^{1/2}}{\sin(\pi\tau)} \left( \frac{(1-z)^2}{\sin(\pi\tau(1-z))^2} - \frac{z^2}{\sin(\pi\tau z)^2} \right) \right. \\ \left. - \frac{2\phi z(1-z) \sin\left(\frac{\phi}{2}(1-2z)\right)}{\sin(\pi\tau(1-z))^{1/2} \sin(\pi\tau z)^{1/2}} \right] = 0. \qquad (5.42) $$

The expression in the bracket vanishes at $z = 0$ and $z = 1$, so the equality indeed holds.

As a result, the lattice expression for $C_{\phi}(x_{21})$ is a linear combination of the two hypergeometric solutions on each of the sub-intervals, as in (5.23). It is a continuous function of $\tau$ and, by construction of the conformal solution in terms of the lattice data, the constants $A_1$, $A_2$, $B_1$ and $B_2$ are equal to those given in (5.25a) and (5.37). This concludes the proof that the lattice and conformal expressions for $C_{\phi}(x_{21})$ are equal.

# 6 Conformal correlators for double entry points

## 6.1 Preliminaries

In the conformal picture, the lattice correlation functions $C_{\phi}(x)$ and $C_{\phi\phi}(x)$ are ratios of conformal correlation functions on the semi-infinite cylinder:

$$ C_{\phi}(x_{21}) = \lim_{x_3 \to i\infty} \frac{\langle \mu_{\phi}(x_1) \mu_{\phi}(x_2) \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathbb{V}}}{\langle \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathbb{V}}}, \qquad (6.1a) $$

$$ C_{\phi\phi}(x_{21}) = \lim_{x_3 \to i\infty} \frac{\langle \omega_{\phi}(x_1) \omega_{\phi}(x_2) \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathbb{V}}}{\langle \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathbb{V}}}. \qquad (6.1b) $$

The evidence given below will support the following claims for the conformal interpretation:

(a) The fields $\mu_{\phi}(z)$ and $\mu_{\phi\phi}(z)$ form a pair of highest weight states of dimension $\Delta = \Delta_{1,3} = 0$. One of these fields belongs to a Kac module and the other belongs to a staggered module.

(b) The field $\omega_{\phi}(z)$ is a logarithmic highest weight state in a staggered module, also with $\Delta = 0$.
The structure and specifics of these representations are described in Section 6.2.

The transformation law for $\omega(z)$ under a conformal transformation $z \mapsto y(z)$ involves its primary partner $\varphi(z)$:

$$ \omega(z) \mapsto \omega(y) = \left( \omega(z) - \lambda \varphi(z) \log \left( \frac{dy}{dz} \right) \right). \qquad (6.2) $$

In CFT derivations, the value of $\lambda$ can be adjusted by changing the normalisation of the field $\omega(z)$, so it is common to fix $\lambda = 1$. In the current context, however, the field $\omega(z)$ arises from the lattice with a natural normalisation, so we leave $\lambda$ free. We will in fact compute its value below, see (6.23). The invariance under the global conformal transformations fixes [12] the two-point functions to

$$ \langle \omega(z_1)\omega(z_2) \rangle_C = \kappa^{\omega\omega} - 2\lambda\kappa^{\omega}\log z_{21}, \quad \langle \omega(z_1)\varphi(z_2) \rangle_C = \kappa^{\omega}, \quad \langle \varphi(z_1)\varphi(z_2) \rangle_C = 0. \qquad (6.3) $$

Likewise, the three-point functions involving two fields $\psi_\alpha$ and one field of the pair $(\varphi, \omega)$ are fixed by conformal invariance to

$$ \langle \varphi(z_1) \psi_\alpha(z_2) \psi_\alpha(z_3) \rangle_C = \frac{\kappa^{\omega\psi\psi}}{z_{32}^{2\Delta}}, \quad \langle \omega(z_1) \psi_\alpha(z_2) \psi_\alpha(z_3) \rangle_C = \frac{\kappa^{\omega\psi\psi} - \lambda\kappa^{\omega\psi\psi} \log\left(\frac{z_{21}z_{31}}{z_{32}}\right)}{z_{32}^{2\Delta}}. \qquad (6.4) $$

In the context of loop models, the boundary field $\varphi(z)$ has a simple interpretation: It inserts a simple half-arc $\bar{\omega}$ at the point $z$ on the boundary. The Dirichlet boundary conditions that we are considering consist of a macroscopic collection of these simple arcs. The insertion of $\varphi(z)$ on such a boundary is thus trivial. This implies that, in any correlator involving other fields, if $\varphi(z)$ appears, it can be simply removed, with the corresponding correlator independent of the position $z$. For instance, from (6.3), we have $\langle \omega(z) \rangle_C = \kappa^\omega$. For generic value of $\beta$, the one-point function $\langle \varphi(z) \rangle_\bar{\omega}$ then corresponds to the partition function of the loop model on the cylinder, with all loops (contractible and non-contractible) having fugacity $\beta$. This partition function vanishes for $\beta = 0$, implying that the one-point function of $\varphi(z)$ equals zero at this value.

## 6.2 Representations of the Virasoro algebra at $c = -2$

This subsection discusses the structure of five modules over the Virasoro algebra with $c = -2$ and $\Delta = 0$. The Loewy diagrams of these modules are given in Figure 4. To start, let us consider the module $M$ defined from the free action of the Virasoro generators $L_m$, $m \le 0$, on two highest-weight states $|\varphi\rangle$ and $|\omega\rangle$ satisfying

$$ L_0|\varphi\rangle = 0, \quad L_0|\omega\rangle = \lambda|\varphi\rangle, \quad L_m|\varphi\rangle = L_m|\omega\rangle = 0, \quad m > 0. \qquad (6.5) $$

The state $|\varphi\rangle$ and its descendants generate a submodule of $M$ isomorphic to the Verma module of highest weight $\Delta = 0$, which we denote by $V$. This submodule is reducible. In particular, it has the singular vector $L_{-1}|\varphi\rangle$ at level 1 and the singular vector $(L_1^2 - 2L_2)L_{-1}|\varphi\rangle$ at level 3. This last state is a descendant of $L_{-1}|\varphi\rangle$, which is consistent with the known structure for this Verma module for $c = -2$ in the form of an imbricated chain [9, 10]. By quotienting $V$ by the singular state at level 1, we obtain the irreducible module $I$. Likewise, by quotienting $V$ by the singular state at level 3, we obtain the so-called Kac module $K$.

We now investigate the descendants of the logarithmic state $|\omega\rangle$ in $M$. Clearly, if $\lambda = 0$, $M$ splits as the direct sum of two copies of $V$. For $\lambda \ne 0$, we observe that the state $L_{-1}|\omega\rangle$ is not singular because $L_1L_{-1}|\omega\rangle = 2\lambda|\varphi\rangle$. It is however sub-singular: if one quotients $M$ by the
<figure><img src="image_8.png" /><figcaption>Figure 4: The Loewy diagrams for the modules V, I, K, M and S.</figcaption></figure>

submodule generated by $|\varphi\rangle$, then the state $L_{-1}|\omega\rangle$ becomes singular. One can check that the state

$$ (L_{-1}^3 - 2L_{-2}L_{-1} + L_{-3}L_0)|\omega\rangle = (L_{-1}^2 - 2L_{-2})L_{-1}|\omega\rangle + \lambda L_{-3}|\varphi\rangle \quad (6.6) $$

is also subsingular: by acting on it with any element of the Virasoro algebra that has a positive multi-index, one obtains a linear combination of states belonging to the tower of states generated by $L_{-1}|\varphi\rangle$. In Rohsiepe's original paper [18], these free modules (like M) were referred to as Jordan-Verma modules. Their structure is known [18, 22], and for the case we are considering, it is as depicted in the fourth panel of Figure 4.

The staggered module S is obtained by quotienting M by the submodule generated by the singular state at level 3 given in (6.6). The three composition factors of S are irreducible modules of weights $\Delta = 0, 0, 1$, and $L_0$ has rank-two Jordan cells tying the states from the two irreducible factors with $\Delta = 0$. We also note that for $\lambda = 0$, the module S decomposes as the direct sum $I \oplus K$. This is used in our calculations below, where we obtain a differential equation for a correlation function involving the field $|\omega\rangle$ in S, and obtain the corresponding result for K by setting $\lambda = 0$.

## 6.3 Differential equation for the four-point function

The correlators on the right sides of (6.1) involve one non-chiral field, which we rewrite as the product of two chiral fields using the method of images. We proceed to compute the correlator

$$ G = \langle \omega_1(z_1) \omega_2(z_2) \psi_\alpha(z_3, z_4) \rangle_{\mathbb{H}} = \langle \omega_1(z_1) \omega_2(z_2) \psi_\alpha(z_3) \psi_\alpha(z_4) \rangle_{\mathbb{C}}. \quad (6.7) $$

We work under the hypothesis that $\omega_1(z)$ and $\omega_2(z)$ are potentially different fields, and in particular that the corresponding constants $\lambda_1$ and $\lambda_2$ dictating their behaviour under conformal transformations, as in (6.2), may not be equal. The correlator $G$ satisfies the logarithmic version [45] of the Ward identities, namely

$$ \left(\sum_{i=1}^{4} \frac{\partial}{\partial z_i}\right) G = 0, \quad \left(\sum_{i=1}^{4} z_i \frac{\partial}{\partial z_i} + \Delta_i + \hat{\delta}_i\right) G = 0, \quad \left(\sum_{i=1}^{4} z_i^2 \frac{\partial}{\partial z_i} + 2\Delta_i z_i + 2z_i \hat{\delta}_i\right) G = 0, \quad (6.8) $$

with $\Delta_1 = \Delta_2 = 0$ and $\Delta_3 = \Delta_4 = \Delta$ given in (5.36). In its action on $G$, the operator $\hat{\delta}_i$ replaces $\omega_i$ by its primary partner $\varphi_i$ and gives zero otherwise. These identities imply that $G$ has the following form:

$$ G = \frac{1}{z_{43}^{2\Delta}} \left[ F(\eta) - \lambda_1 \kappa_2^{\omega\psi\psi} \log\left(\frac{z_{41}z_{31}}{z_{43}}\right) - \lambda_2 \kappa_1^{\omega\psi\psi} \log\left(\frac{z_{42}z_{32}}{z_{43}}\right) + \lambda_1 \lambda_2 \kappa^{\psi\psi} \log\left(\frac{z_{41}z_{31}}{z_{43}}\right) \log\left(\frac{z_{42}z_{32}}{z_{43}}\right) \right], \quad (6.9) $$
where $\eta$ is the cross-ratio, defined in (5.10), and $\kappa_1^{\omega\psi\psi}$ and $\kappa_2^{\omega\psi\psi}$ are the constants appearing in the three-point function (6.4) for $\omega_1(z)$ and $\omega_2(z)$. Mapping this result to the cylinder, we find

$$ \langle \omega_1(x_1)\omega_2(x_2)\psi_\alpha(x_3, x_4) \rangle_V = \left(\frac{n}{\pi} s_{43}\right)^{-2\Delta} \left[ \begin{aligned} & F(\eta) - \lambda_1 \kappa_2^{\omega\psi\psi} \log\left(\frac{n s_{41}s_{31}}{\pi s_{43}}\right) - \lambda_2 \kappa_1^{\omega\psi\psi} \log\left(\frac{n s_{42}s_{32}}{\pi s_{43}}\right) \\ & + \lambda_1 \lambda_2 \kappa^{\omega\psi\psi} \log\left(\frac{n s_{41}s_{31}}{\pi s_{43}}\right) \log\left(\frac{n s_{42}s_{32}}{\pi s_{43}}\right) \end{aligned} \right], \quad (6.10) $$

where we use the compact notation

$$ s_{ij} = \sin\left(\frac{\pi x_{ij}}{n}\right). \quad (6.11) $$

For $x_4 = x_3^*$ and $x_3 \to i\infty$, the ratios $s_{41}s_{31}/s_{43}$ and $s_{42}s_{32}/s_{43}$ tend to $i/2$, and as a result we have

$$ \lim_{x_3 \to i\infty} \frac{\langle \omega_1(x_1)\omega_2(x_2)\psi_\alpha(x_3, x_3^*) \rangle_V}{\langle \psi_\alpha(x_3, x_3^*) \rangle_V} = \frac{1}{\kappa^{\omega\psi\psi}} \left[ F(\eta) - (\lambda_1 \kappa_2^{\omega\psi\psi} + \lambda_2 \kappa_1^{\omega\psi\psi}) \log\left(\frac{in}{2\pi}\right) + \lambda_1 \lambda_2 \kappa^{\omega\psi\psi} \log^2\left(\frac{in}{2\pi}\right) \right]. \quad (6.12) $$

In the staggered module, the state $|\omega\rangle$ has the null descendant (6.6) at level 3, which results in the following differential equation for $G$:

$$ (\mathcal{L}_{-1}^3 - 2\mathcal{L}_{-2}\mathcal{L}_{-1} + \mathcal{L}_{-3}\mathcal{L}_0)G = 0, \quad (6.13) $$

where

$$ \mathcal{L}_{-n} = \sum_{i=2}^{4} \frac{(n-1)(\Delta_i + \hat{\delta}_i)}{(z_i - z_1)^n} - \frac{1}{(z_i - z_1)^{n-1}} \frac{\partial}{\partial z_i}. \quad (6.14) $$

The partial differential equation (6.13) simplifies to an ordinary differential equation satisfied by $F(\eta)$:

$$ \eta^3(1-\eta)^2 F'''(\eta) + 2(1-2\eta)(1-\eta)\eta^2 F''(\eta) - 2\eta^2(1-(1-\Delta)\eta)F'(\eta) = (2-\eta)\lambda_1\lambda_2\kappa^{\omega\psi\psi}. \quad (6.15) $$

## 6.4 Correlation function for the Kac module

We solve the differential equation (6.15) for the case $\lambda_1\lambda_2 = 0$. This means that either $\lambda_1$ or $\lambda_2$ vanishes. The analysis below is independent of which one is set to zero. We choose $\lambda_1 = 0$, for which $\omega_1$ becomes a highest state $\mu$ of a Kac module $K$, and we write $\kappa_1^{\omega\psi\psi} = \kappa_1^{\mu\psi\psi}$. Comparing (4.23) and (6.12), we see that the terms proportional to $\log n$ are equal for

$$ \lambda_2 \frac{\kappa_1^{\mu\psi\psi}}{\kappa^{\omega\psi\psi}} = -\frac{2}{\pi}, \quad (6.16) $$

implying that $\lambda_2 \neq 0$. We note that, for $\tau \to 0$, the term proportional to $\log \tau$ in (4.26a) has the overall constant $2/\pi$, which is identical to the universal value predicted in [42, equation (76)]. The result obtained here is however more general, as the correlation function $C_{\psi_\alpha}(x_{21})$ involves the field $\psi_\alpha(z, z^*)$ as well.
We use the notation $F_0(\eta) = F(\eta)|_{\lambda_1\lambda_2=0}$ for the solution to (6.15) in the homogeneous case. The differential equation is in fact a second-order homogeneous differential equation for $F'_0(\eta)$. Setting $\Delta$ to its value given in (5.36), the general solution is

$$F'_0(\eta) = A_1 \frac{(1-\eta)^{-\frac{1}{2}(1-\frac{\phi}{\pi})}}{\eta} + A_2 \frac{(1-\eta)^{-\frac{1}{2}(1+\frac{\phi}{\pi})}}{\eta}, \quad (6.17)$$

where $A_1$ and $A_2$ are constants to be determined. According to (6.1a) and (6.12), we have

$$C(x_{21}) = \frac{1}{\kappa^{\psi\psi}} F_0(\eta), \quad \eta = 1 - e^{2\pi i \tau}, \quad (6.18)$$

and therefore

$$\frac{\partial C(x_{21})}{\partial \tau} = \frac{1}{\kappa^{\psi\psi}} \frac{\partial F_0}{\partial \eta} \frac{\partial \eta}{\partial \tau} = \frac{\pi}{\kappa^{\psi\psi}} \left( \frac{A_1 e^{i\phi\tau} + A_2 e^{-i\phi\tau}}{\sin \pi\tau} \right). \quad (6.19)$$

Comparing this with the lattice result, we see that the derivative of the expression (4.23) with respect to $\tau$ is easily evaluated because $\tau$ only appears in the upper integration limit:

$$\frac{\partial C(x_{21})}{\partial \tau} = \frac{2 \cos(\phi \tau)}{\sin(\pi \tau)}. \quad (6.20)$$

We find that (6.19) and (6.20) coincide for

$$A_1 = A_2 = \frac{\kappa^{\psi\psi}}{\pi}. \quad (6.21)$$

This consistency between the lattice and conformal results supports our claim that the fields $\mu_{\bar{\nu}}$ and $\mu_{\bar{\tau}}$ form a pair of highest-weight states of dimension $\Delta = 0$, one of which belongs to a Kac module K and the other to a staggered module S.

## 6.5 Correlation function for the staggered module

For the correlator (6.1b), the four-point correlator contains two copies of the same logarithmic field: $\omega_1 = \omega_2 = \omega_{\bar{\nu}}$. We therefore write $\lambda_1 = \lambda_2 = \lambda$ and $\kappa_1^{\omega\psi\psi} = \kappa_2^{\omega\psi\psi} = \kappa^{\omega\psi\psi}$. For $\lambda_1\lambda_2 = \lambda^2 \neq 0$, the differential equation (6.15) has an extra inhomogeneous term. The complete solution is equal to the two-parameter solution of the homogeneous equation plus a second function that takes into account the inhomogeneous term:<fnref n="1" />

$$F'(\eta) = F'_0(\eta) + \frac{2\lambda^2 \kappa^{\psi\psi}}{\frac{\phi}{\pi}(1-\eta)^2} \left[ 2 \left( \frac{2F_1(\frac{1}{2} + \frac{\phi}{2\pi}, 1, \frac{5}{2} + \frac{\phi}{2\pi} \middle| \frac{1}{1-\eta})}{3 + \frac{\phi}{\pi}} - \frac{2F_1(\frac{1}{2} - \frac{\phi}{2\pi}, 1, \frac{5}{2} - \frac{\phi}{2\pi} \middle| \frac{1}{1-\eta})}{3 - \frac{\phi}{\pi}} \right) \right. \\ \left. - \eta \left( \frac{2F_1(\frac{1}{2} + \frac{\phi}{2\pi}, 1, \frac{3}{2} + \frac{\phi}{2\pi} \middle| \frac{1}{1-\eta})}{1 + \frac{\phi}{\pi}} - \frac{2F_1(\frac{1}{2} - \frac{\phi}{2\pi}, 1, \frac{3}{2} - \frac{\phi}{2\pi} \middle| \frac{1}{1-\eta})}{1 - \frac{\phi}{\pi}} \right) \right], \quad (6.22)$$

where the constants $A_1$ and $A_2$ of $F'_0(\eta)$ in (6.17) remain to be fixed.

We fix the remaining unknowns using the lattice expression. Comparing (4.28) with (6.12), we see that the $\log^2(n)$ and $\log(n)$ terms are identical for

$$\lambda = \frac{(2 \cos \frac{\phi}{2})^{1/2}}{\pi}, \quad \frac{\kappa^{\omega\psi\psi}}{\kappa^{\psi\psi}} = \frac{1}{\pi(2 \cos \frac{\phi}{2})^{1/2}} \left( 2 \cos(\frac{\phi}{2}) \left( \frac{i\pi}{2} - \log 2 + \psi(\frac{1}{2} + \frac{\phi}{2\pi}) \right) - \pi \sin(\frac{\phi}{2}) \right). \quad (6.23)$$

<footnote n="1">¹This solution to the differential equation (6.15) was obtained using Mathematica.</footnote>
It may seem odd that $\lambda$ depends on $\phi$, as this parameter dictates the behaviour of the boundary field $\omega_{\tau\tau}(z)$, which a priori one expects not to depend on $\alpha$. This will be discussed further in Section 7. We also note that the structure constant $\kappa^{\omega\psi\psi}$ originates from the correlator, given in (6.4), of a bulk field and two boundary fields. So it is perhaps not surprising that the ratio $\kappa^{\omega\psi\psi}/\kappa^{\psi\psi}$ is complex. The remaining constant factors in (4.28), which are independent of $n$ and $\tau$, are accounted for in the conformal solution by the integration constant that appears when one integrates (6.22) to obtain $F(\eta)$.

To fix the constants $A_1$ and $A_2$, we compare the asymptotic behaviour of the lattice and conformal expressions for $F'(\eta)$ in the neighbourhood of $\eta = 1$, corresponding to $\tau \to 0^+$. For the lattice expression, (4.29) yields

$$ \frac{\partial C_{\tau\tau}}{\partial \tau} \xrightarrow{n \gg 1, \tau \to 0^+} \frac{2}{\pi^2 \tau} \left[ \pi \sin\left(\frac{\phi}{2}\right) - 2 \cos\left(\frac{\phi}{2}\right) \left( \gamma + \log 2 + \log \pi + \psi\left(\frac{1}{2} + \frac{\phi}{2\pi}\right) + \log \tau \right) \right]. \quad (6.24) $$

For the conformal expression, we extract the behaviour of the hypergeometric functions in (6.22) at $z = 1$ by using the relations

$$ {}_2F_1(a, b, a+b|z) = -\frac{\Gamma(a+b)}{\Gamma(a)\Gamma(b)}H(a, b|1-z), \quad (6.25a) $$

$$ {}_2F_1(a, b, a+b+1|z) = -\frac{\Gamma(a+b+1)}{\Gamma(a+1)\Gamma(b+1)}(1-z)\frac{\partial H(a, b|1-z)}{\partial z}, \quad (6.25b) $$

and the series expansion (5.16) of $H(a, b|z)$. This yields

$$ \begin{aligned} \frac{\partial C_{\tau\tau}}{\partial \tau} &= \frac{1}{\kappa^{\psi\psi}} \frac{\partial \eta}{\partial \tau} \frac{\partial F}{\partial \eta} \\ &= -\frac{1}{\tau} \left( A_1 + A_2 + 2\lambda^2 \left( \gamma + \log \pi + \log 2 + \psi\left(\frac{1}{2} + \frac{\phi}{2\pi}\right) + \log \tau + \frac{i\pi}{2} + \frac{\pi}{2} \tan\left(\frac{\phi}{2\pi}\right) \right) \right) \\ &\quad - 2i\phi \left( A_1 - A_2 + \pi\lambda^2 \tan\left(\frac{\phi}{2}\right) \right). \end{aligned} \quad (6.26) $$

With $\lambda$ evaluated to its value in (6.23), we find that (6.24) and (6.26) are equal for

$$ \frac{A_1}{\kappa^{\psi\psi}} = \frac{e^{-i\phi/2}}{i\pi}, \quad \frac{A_2}{\kappa^{\psi\psi}} = \frac{e^{i\phi/2}}{i\pi}. \quad (6.27) $$

The final conformal expression for $\frac{\partial C_{\tau\tau}}{\partial \tau}$ is

$$ \begin{aligned} \frac{\partial C_{\tau\tau}}{\partial \tau} = & \frac{8 \cos\left(\frac{\phi}{2}\right)}{i\phi(1 - e^{2\pi i\tau})^2} \left[ 2 \left( \frac{2F_1\left(\frac{1}{2} + \frac{\phi}{2\pi}, 1, \frac{5}{2} + \frac{\phi}{2\pi} \middle| e^{-2\pi i\tau}\right)}{3 + \frac{\phi}{\pi}} - \frac{2F_1\left(\frac{1}{2} - \frac{\phi}{2\pi}, 1, \frac{5}{2} - \frac{\phi}{2\pi} \middle| e^{-2\pi i\tau}\right)}{3 - \frac{\phi}{\pi}} \right) \right. \\ & \left. - (1 - e^{2\pi i\tau}) \left( \frac{2F_1\left(\frac{1}{2} + \frac{\phi}{2\pi}, 1, \frac{3}{2} + \frac{\phi}{2\pi} \middle| e^{-2\pi i\tau}\right)}{1 + \frac{\phi}{\pi}} - \frac{2F_1\left(\frac{1}{2} - \frac{\phi}{2\pi}, 1, \frac{3}{2} - \frac{\phi}{2\pi} \middle| e^{-2\pi i\tau}\right)}{1 - \frac{\phi}{\pi}} \right) \right] \\ & + \frac{2i \cos\left(\phi(\tau - \frac{1}{2})\right)}{\sin(\pi\tau)}. \end{aligned} \quad (6.28) $$

It is in fact not obvious that this function is real for $\tau \in (0, 1)$, but indeed it is. Plotting this function alongside the lattice expression for $\frac{\partial C_{\tau\tau}}{\partial \tau}$, we find that the two functions precisely coincide, on the full interval. To prove the equality, we show that the lattice expression also satisfies the differential equation (6.15). In terms of $\tau$, this equation reads

$$ \left( \frac{\partial^2}{\partial \tau^2} + 2\pi \cot(\pi\tau) \frac{\partial}{\partial \tau} + 8\pi^2 \Delta \right) \frac{\partial C_{\tau\tau}}{\partial \tau} = 2\pi^3 \lambda^2 \frac{\cos(\pi\tau)}{\sin^3(\pi\tau)}. \quad (6.29) $$
To check this, we take the derivative of (4.28) and simplify it to

$$ \frac{\partial C_{\tau\tau}}{\partial \tau} = -8 \cos\left(\frac{\phi}{2}\right) \int_{0}^{\tau} dt \frac{\sin\left(\frac{t\phi}{2}\right) \sin\left((\tau - \frac{t}{2})\phi\right)}{\sin(\pi\tau) \sin(\pi t)} + 2 \sin\left(\frac{\phi}{2}\right) \frac{\cos(\phi\tau)}{\sin(\pi\tau)} \\ - \frac{4 \cos\left(\frac{\phi}{2}\right)}{\pi} \frac{\cos(\phi\tau)}{\sin(\pi\tau)} \left(\gamma + 2 \log 2 + \log \tan\left(\frac{\pi\tau}{2}\right) + \psi\left(\frac{1}{2} + \frac{\phi}{2\pi}\right)\right). \quad (6.30) $$

Applying the operator that appears on the left-hand side of (6.29) to the terms proportional to $\frac{\cos(\phi\tau)}{\sin(\pi\tau)}$ gives zero, as this is a solution of the homogeneous differential equation discussed in Section 6.4. For the remaining terms, after simplifications we find

$$ \left( \frac{\partial^2}{\partial \tau^2} + 2\pi \cot(\pi\tau) \frac{\partial}{\partial \tau} + 8\pi^2 \Delta \right) \left[ \frac{4 \cos\left(\frac{\phi}{2}\right)}{\pi} \frac{\cos(\phi\tau)}{\sin(\pi\tau)} \log \tan\left(\frac{\pi\tau}{2}\right) \right] \\ = 8\pi \cos\left(\frac{\phi}{2}\right) \left( \frac{\cos(\pi\tau) \sin\left(\frac{\phi\tau}{2}\right)^2}{\sin^3(\pi\tau)} - \frac{\phi}{\pi} \frac{\sin(\phi\tau)}{\sin^2(\pi\tau)} \right) - 2\pi^3 \lambda^2 \frac{\cos(\pi\tau)}{\sin^3(\pi\tau)} \quad (6.31) $$

and

$$ \int_0^1 dt \left( \frac{\partial^2}{\partial \tau^2} + 2\pi \cot(\pi\tau) \frac{\partial}{\partial \tau} + 8\pi^2 \Delta \right) \left[ \frac{\tau \sin\left(\frac{\phi\tau t}{2}\right) \sin\left(\phi\tau(1-\frac{t}{2})\right)}{\sin(\pi\tau t) \sin(\pi\tau)} \right] \\ = \int_0^1 dt \frac{d}{dt} \left[ \frac{t}{\sin(\pi\tau) \sin(\pi\tau t)} \left( \phi \sin(\phi\tau) - (1-\frac{t}{2})\phi \sin(\phi\tau(1-t)) \right) \right. \\ \left. - \pi t \cot(\pi\tau t) \sin\left(\frac{\phi\tau t}{2}\right) \sin\left(\phi\tau(1-\frac{t}{2})\right) \right] \\ = \phi \frac{\sin(\phi\tau)}{\sin^2(\pi\tau)} - \pi \frac{\cos(\pi\tau) \sin^2\left(\frac{\phi\tau}{2}\right)}{\sin^3(\pi\tau)}. \quad (6.32) $$

Combining these results, we confirm that the lattice expression (6.30) indeed satisfies (6.29). The lattice and conformal solutions therefore satisfy the same differential equation. Recalling that the constants $A_1$ and $A_2$ were fixed so that the lattice and conformal expressions have the same asymptotics, we conclude that the two expressions coincide.

# 7 Operator product expansion and structure constants

## 7.1 Analysis for the fields $\phi(z)$

In this section, we use the lattice expressions for the correlators $C_{\tau\tau}(x)$ and $C_{\tau\tau}(x)$ to obtain conformal structure constants and ratios thereof. We claim that the fields $\phi_{\tau}(z)$ and $\phi_{\tau}(z)$ satisfy the following operator product expansions (OPEs):

$$ \phi_{\tau}(z_1)\phi_{\tau}(z_2) = z_{21}^{1/4} \kappa^{\phi\phi} \varphi(z_2) + \dots, \quad (7.1a) $$

$$ \phi_{\tau}(z_1)\phi_{\tau}(z_2) = z_{21}^{1/4} \left( \kappa^{\phi\phi\omega} (\omega_{\tau\tau}(z_2) + \lambda \varphi(z_2) \log z_{21}) + \kappa^{\phi\phi} \varphi(z_2) \right) + \dots. \quad (7.1b) $$

From (6.3), this yields the following two-point functions:

$$ \langle \phi_{\tau}(z_1) \phi_{\tau}(z_2) \rangle_C = 0, \quad \langle \phi_{\tau}(z_1) \phi_{\tau}(z_2) \rangle_C = z_{21}^{1/4} \kappa^{\phi\phi\omega} \kappa^{\omega}. \quad (7.2) $$
We therefore see that the constant $\tilde{\kappa}^{\phi\phi}$ appearing in (5.3) vanishes for the first correlator, but equals $\kappa^{\phi\phi\omega}\kappa^{\omega}$ for the second correlator.

With these OPEs, we can easily obtain the behaviour of the four-point functions studied in Section 5.2 in the regime where $z_1$ approaches $z_2$. Indeed, in this regime, the four-point functions reduce to two- and three-point functions of the form (6.4), and we find

$$ \langle \phi_{\tau}(z_1) \phi_{\tau}(z_2) \psi_{\alpha}(z_3) \psi_{\alpha}(z_4) \rangle_C \xrightarrow{z_1 \to z_2} \frac{z_{21}^{1/4}}{z_{43}^{2\Delta}} \kappa^{\phi\phi} \kappa^{\psi\psi}, \quad (7.3a) $$

$$ \langle \phi_{\tau}(z_1) \phi_{\tau}(z_2) \psi_{\alpha}(z_3) \psi_{\alpha}(z_4) \rangle_C \xrightarrow{z_1 \to z_2} \frac{z_{21}^{1/4}}{z_{43}^{2\Delta}} \left[ \kappa_{\tau\tau}^{\phi\phi\omega} \kappa^{\omega\psi\psi} + \kappa_{\tau\tau}^{\phi\phi\omega} \kappa^{\psi\psi\lambda} \log\left(\frac{z_{21}z_{43}}{z_{32}z_{42}}\right) + \kappa_{\tau\tau}^{\phi\phi} \kappa^{\psi\psi} \right]. \quad (7.3b) $$

The corresponding correlators on the cylinder are obtained by applying the transformation laws (5.2). The regime $z_1 \to z_2$ then corresponds to $x_2 - x_1 \ll n$, or equivalently to $\tau \to 0^+$. In this regime, we also have $s_{21} \simeq \pi\tau$ and $\frac{x_1}{n} \simeq \frac{x_2}{n}$, where $\simeq$ indicates an equality up to corrections of order $\tau$. We obtain expressions for the correlators $C_{\tau\tau}(x_{21})$ and $C_{\tau\tau}(x_{21})$ by dividing by $\langle \psi_{\alpha}(x_3) \psi_{\alpha}(x_4) \rangle_{\tau}$, setting $x_4 = x_3^*$ and sending $x_3$ to $\pm i\infty$. This yields

$$ \begin{aligned} C_{\tau\tau}(x_{21}) &\xrightarrow{\tau \to 0^+} (n\tau)^{1/4} \kappa_{\tau\tau}^{\phi\phi}, \\ C_{\tau\tau}(x_{21}) &\xrightarrow{\tau \to 0^+} (n\tau)^{1/4} \left[ \kappa_{\tau\tau}^{\phi\phi\omega} \frac{\kappa^{\omega\psi\psi}}{\kappa^{\psi\psi}} + \kappa_{\tau\tau}^{\phi\phi\omega} \lambda \log(-2\pi i\tau) + \kappa_{\tau\tau}^{\phi\phi} \right]. \end{aligned} \quad (7.4) $$

We compare these expressions with (3.29) and (3.33), assuming that the structure constants are independent of $n$ and $\tau$. We find

$$ \kappa_{\tau\tau}^{\phi\phi} = \frac{\pi G^2(\frac{1}{2})}{2^{1/4}}, \qquad \kappa_{\tau\tau}^{\phi\phi\omega} = -\frac{(2 \cos(\frac{\phi}{2}))^{1/2} \pi G^2(\frac{1}{2})}{2^{1/4}}, \qquad \kappa_{\tau\tau}^{\phi\phi} = \frac{2 \cos(\frac{\phi}{2}) G^2(\frac{1}{2})}{2^{1/4}} (-\gamma + \log 2), \quad (7.5) $$

where we used (6.23) for the ratio $\kappa^{\omega\psi\psi}/\kappa^{\psi\psi}$.

We have therefore computed the structure constants for the fields $\phi_{\tau}(z)$ and $\phi_{\tau}(z)$ in their OPEs (7.1). We see that the fusion of these fields is non-abelian. We expect that this is a general feature of logarithmic CFT. If the fusion product of two fields involves a logarithmic field, then interchanging the order of the two fields induces a non-trivial change of the value of the correlator. The (weak) non-locality of the two fields endows them with a non-trivial monodromy, which can for instance be compared with that of disorder operators in the Ising model or in more general parafermionic field theories [64]. The unusual distinction made here between defect insertions that mark the start and end of boundary arcs goes beyond the introduction of a mere phase. We note that the same behaviour was observed for the abelian sandpile model in [33], wherein two fields $\mu_{D,N}(z)$ and $\mu_{N,D}(z)$ that mark the transition between dissipative and non-dissipative sites on the boundary satisfy different fusion rules according to the order in which they are fused.

Finally, we note that, similarly to the value of $\lambda$ in (6.23), the structure constants $\kappa_{\tau\tau}^{\phi\phi\omega}$ and $\kappa_{\tau\tau}^{\phi\phi}$ depend on $\alpha$. This is somewhat unexpected, as one might expect the boundary fields $\phi_{\tau}(z)$, $\phi_{\tau}(z)$ and $\omega_{\tau}(z)$ not to depend on $\alpha$. We believe that the resolution to this conundrum lies in the non-local nature of the field $\psi_{\alpha}(z)$. Its action on the lattice is indeed non-local, as it modifies the weight of all loops that encircle its insertion point, independently of how distant these are. We believe that, in the conformal interpretation, the non-locality of this field means that it modifies some conformal properties of the other fields of the theory, namely their structure constants and transformation laws, even if these are inserted a large distance away.
## 7.2 Analysis for the fields $\mu(z)$

The results of Section 6.4 are derived with the assumption that $\mu_{\mathcal{Z}}(z)$ is a field that belongs to a Kac module, whereas $\mu_{\mathcal{Z}}(z)$ belongs to a staggered module. (We recall that these identifications could equivalently have been made the other way around). Here, we claim that the OPEs for these fields are of the form

$$ \mu_{\mathcal{Z}}(z_1)\mu_{\mathcal{Z}}(z_2) = \kappa^{\mu\mu} \varphi(z_2) + \dots, \qquad (7.6a) $$

$$ \mu_{\mathcal{Z}}(z_1)\mu_{\mathcal{Z}}(z_2) = \kappa^{\mu\mu\tilde{\omega}}(\tilde{\omega}_{\mathcal{Z}}(z_2) + \lambda_2 \varphi(z_2) \log z_{21}) + \kappa^{\mu\mu} \varphi(z_2) + \dots. \qquad (7.6b) $$

We believe that the field $\tilde{\omega}_{\mathcal{Z}}(z)$ is a logarithmic field with weight $\Delta = 0$ in a rank-two Jordan cell. It is potentially different from the field $\omega_{\mathcal{Z}}(z)$, although the evidence below cannot rule out that they are in fact identical. Similarly to $\mu_{\mathcal{Z}}(z)$, $\tilde{\omega}_{\mathcal{Z}}(z)$ transforms conformally according to (6.2), with the constant $\lambda$ replaced by $\lambda_2$. Its three-point correlator with two copies of the field $\psi_\alpha(z)$ is

$$ \langle \tilde{\omega}_{\mathcal{Z}}(z_1) \psi_\alpha(z_2) \psi_\alpha(z_3) \rangle_{\mathcal{C}} = \frac{\kappa^{\tilde{\omega}\psi\psi} - \lambda_2 \kappa^{\psi\psi} \log\left(\frac{z_{21}z_{31}}{z_{32}}\right)}{z_{32}^{2\Delta}}, \qquad (7.7) $$

where $\kappa^{\tilde{\omega}\psi\psi}$ is a constant.

In the regime $z_1 \to z_2$, we have the following four-point functions:

$$ \langle \mu_{\mathcal{Z}}(z_1) \mu_{\mathcal{Z}}(z_2) \psi_\alpha(z_3) \psi_\alpha(z_4) \rangle_{\mathcal{C}} \xrightarrow{z_1 \to z_2} \frac{1}{z_{43}^{2\Delta}} \kappa^{\mu\mu} \kappa^{\psi\psi}, \qquad (7.8a) $$

$$ \langle \mu_{\mathcal{Z}}(z_1) \mu_{\mathcal{Z}}(z_2) \psi_{\alpha}(z_3) \psi_{\alpha}(z_4) \rangle_{\mathcal{C}} \xrightarrow{z_1 \to z_2} \frac{1}{z_{43}^{2\Delta}} \left[ \kappa_{\mathcal{Z}}^{\mu\mu\tilde{\omega}} \kappa^{\tilde{\omega}\psi\psi} + \lambda_2 \kappa_{\mathcal{Z}}^{\mu\mu\tilde{\omega}} \kappa^{\psi\psi} \log\left(\frac{z_{21}z_{43}}{z_{32}z_{42}}\right) + \kappa_{\mathcal{Z}}^{\mu\mu} \kappa^{\psi\psi} \right]. \qquad (7.8b) $$

We map these correlators on the cylinder, divide by $\langle \psi_\alpha(x_3)\psi_\alpha(x_4) \rangle_\mathcal{V}$, take the limit $x_3 \to i\infty$ with $x_4 = x_3^*$, and obtain

$$ C_{\mathcal{Z}}(x_{21}) \xrightarrow{\tau \to 0^+} \kappa_{\mathcal{Z}}^{\mu\mu} - \lambda_2 \frac{\kappa_1^{\mu\psi\psi}}{\kappa^{\psi\psi}} \log(n\tau), \qquad (7.9a) $$

$$ C_{\mathcal{Z}}(x_{21}) \xrightarrow{\tau \to 0^+} \kappa_{\mathcal{Z}}^{\mu\mu\tilde{\omega}} \frac{\kappa^{\tilde{\omega}\psi\psi}}{\kappa^{\psi\psi}} + \lambda_2 \kappa_{\mathcal{Z}}^{\mu\mu\tilde{\omega}} \log(-2\pi i \tau) + \kappa_{\mathcal{Z}}^{\mu\mu} - \lambda_2 \frac{\kappa_1^{\mu\psi\psi}}{\kappa^{\psi\psi}} \log(n\tau). \qquad (7.9b) $$

In each of these expressions, the last term originates from the logarithmic transformation law for the field $\mu_{\mathcal{Z}}(z)$, and $\kappa_1^{\mu\psi\psi}$ is the constant that appears in the correlator $\langle \mu_{\mathcal{Z}}(z_1)\psi_\alpha(z_2)\psi_\alpha(z_3) \rangle_{\mathcal{C}}$. Comparing these expressions with (4.31a) and (4.31b), we equate the different terms, assuming that the structure constants are independent of $n$ and $\tau$. We find

$$ \kappa_{\mathcal{Z}}^{\mu\mu} = \frac{2}{\pi}(\log 2 + \gamma), \qquad \lambda_2 \kappa_{\mathcal{Z}}^{\mu\mu\tilde{\omega}} = -\frac{(2 \cos(\frac{\phi}{2}))^2}{\pi}, \qquad (7.10a) $$

$$ \kappa_{\mathcal{Z}}^{\mu\mu\tilde{\omega}} \frac{\kappa^{\tilde{\omega}\psi\psi}}{\kappa^{\psi\psi}} + \kappa_{\mathcal{Z}}^{\mu\mu} = \frac{2}{\pi} \left( \gamma + \log 2 + 2 \cos^2\left(\frac{\phi}{2}\right) \left( \log \pi - \frac{i\pi}{2} \right) \right) + \hat{K}, \qquad (7.10b) $$

where we used (6.16).

The lattice derivations of Section 4 are therefore insufficient to fix individually all the structure constants appearing in (7.6). In comparison, all the constants for the OPEs of the fields
$\phi_{\tau}(z)$ and $\phi_{\tau}(z)$ were obtained in (7.5). This was possible because the four-point correla- tor involving the field $\omega_{\tau}(z)$ was computed from the lattice in Section 4. Presumably, by computing

$$ \lim_{x_3 \to i\infty} \frac{\langle \tilde{\omega}_{\tau\tau}(x_1) \tilde{\omega}_{\tau\tau}(x_2) \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathcal{V}}}{\langle \psi_{\alpha}(x_3, x_3^*) \rangle_{\mathcal{V}}} \quad (7.11) $$

from the lattice, one would be able to compute $\lambda_2$ and $\kappa^{\tilde{\omega}\psi\psi}/\kappa^{\psi\psi}$ directly and then solve for the remaining unknowns, namely $\kappa^{\mu\mu\tilde{\omega}}$, $\kappa^{\mu\mu}$ and $\kappa_1^{\mu\psi\psi}/\kappa^{\psi\psi}$.

Thus, in the present state of affairs, these results do not allow us to determine whether $\omega_{\tau}(z)$ and $\tilde{\omega}_{\tau}(z)$ are the same field. If they turn out not to be, one may have to consider the possibility that the LCFT for the model of critical dense polymers involves an infinite number of boundary logarithmic fields of conformal dimension $\Delta = 0$, one for each link state of the form $\tau$, $\tau\tau$, $\tau\tau\tau$, etc.

## 7.3 Analysis for the fields $\omega(z)$

Following the notation of [33], we write the OPE of the field $\omega_{\tau}(z)$ with itself as

$$ \omega_{\tau}(z_1)\omega_{\tau}(z_2) = (a - 2\lambda \log z_{21})\omega_{\tau}(z_2) + (b + a\lambda \log z_{21} - \lambda^2 \log^2 z_{21})\varphi(z_2) + \dots \quad (7.12) $$

where $a$ and $b$ are constants. The two-point function is then given by (6.3) with $\kappa^{\omega\omega} = a\kappa^{\omega}$. We use this OPE to obtain the four-point function involving two copies of $\omega_{\tau}(z)$ and two copies of $\psi_{\alpha}(z)$, in the regime $z_1 \to z_2$:

$$ \langle \omega_{\tau}(z_1) \omega_{\tau}(z_2) \psi_{\alpha}(z_3) \psi_{\alpha}(z_4) \rangle_C \xrightarrow{z_1 \to z_2} \frac{1}{z_{43}^{2\Delta}} \left[ (a - 2\lambda \log z_{21}) (\kappa^{\omega\psi\psi} - \lambda \kappa^{\psi\psi} \log \left( \frac{z_{32}z_{42}}{z_{43}} \right)) \right. \\ \left. + (b + a\lambda \log z_{21} - \lambda^2 \log^2 z_{21}) \kappa^{\psi\psi} \right]. \quad (7.13) $$

We obtain the same correlation function on $\mathcal{V}$ using the conformal mapping. We divide by $\langle \psi_{\alpha}(x_3) \psi_{\alpha}(x_4) \rangle_{\mathcal{V}}$, set $x_4 = x_3^*$, send $x_3 \to i\infty$ and find

$$ C_{\tau\tau}(x_{21}) \xrightarrow{\tau \to 0^+} (a - 2\lambda \log(n\tau))(\frac{\kappa^{\omega\psi\psi}}{\kappa^{\psi\psi}} + \lambda \log(-2\pi i\tau)) + b + \lambda^2 \log^2(n\tau). \quad (7.14) $$

We compare this result with (4.29) and equate separately the terms in $\log^2 n$, $\log^2 \tau$, $\log n$, $\log \tau$ and the constant term. This confirms the values obtained in (6.23) and yields

$$ a = -\frac{2}{\pi} \left( 2 \cos\left(\frac{\phi}{2}\right) \right)^{1/2} (\gamma + \log 2), \quad b = -\frac{2 \cos\left(\frac{\phi}{2}\right)}{\pi^2} (\gamma + \log 2)^2 = -\frac{a^2}{4}. \quad (7.15) $$

Interestingly, comparing with the results of [33], we see that the values of $a$ and $b$ are different, but that the relation $b = -\frac{a^2}{4}$ is the same, leading us to wonder whether this is a universal feature of logarithmic CFTs at $c = -2$.

## 8 Discussion and conclusion

In this paper, we computed correlation functions for the model of critical dense polymers on a semi-infinite cylinder. These were obtained from an exact lattice derivation using the XX spin-chain representation of the enlarged periodic Temperley-Lieb algebra. The asymptotic behaviour as the system size $n$ grows to infinity was obtained in terms of integral formulae
involving the scaling parameter $\tau = \frac{x-1}{n} \in (0, 1)$. For small $\tau$, the leading behaviour of these correlators are proportional to $\tau^{1/4}$, $\tau^{1/4} \log \tau$, $\log \tau$ and $\log^2 \tau$. This logarithmic dependence of the correlation functions upon the position of the fields is the defining feature of logarithmic conformal field theories. This behaviour involving the square of the logarithm of the distance is rather uncommon, but has previously been observed in the Potts model for certain Fortuin-Kasteleyn probabilities [42].

In the conformal interpretation, these lattice observables were understood to be ratios of conformal correlation functions: a four-point function divided by a two-point function. Using conformal invariance, we derived differential equations for the four-point functions and solved these equations. We found a perfect agreement with the lattice results, with the logarithmic behaviour arising in two ways: (i) from degenerate solutions to the hypergeometric differential equation, and (ii) from the logarithmic generalisation of the differential equations for highest-weight fields in rank-two Jordan cells that have null descendants.

Admittedly, our conformal derivation is really a hybrid approach, as the constants that arose in the solutions of the differential equations were fixed using the asymptotics of the lattice results. It should in fact be possible to fix these constants directly using the conformal bootstrap [3, 65], and thus to avoid any input coming from the lattice. It would certainly be interesting to understand how this method applies to the case at hand.

We also used the lattice results to compute the structure constants that appear in the operator product expansions of the boundary fields of the model. An intriguing feature that we found pertains to the role played by the field $\psi_\alpha(z)$. Indeed, the fields $\phi_{\tau^-}(z)$, $\phi_{\tau^+}(z)$, $\mu_{\tau^-}(z)$, $\mu_{\tau^+}(z)$ and $\omega_{\tau^-}(z)$ are boundary fields that are expected not to depend on $\alpha$, but the calculations of Section 7 show that the structure constants appearing in their OPEs do in fact depend on $\alpha$. These are of the form $\alpha^\ell$ times constants, with $\ell \in \frac{1}{2}\mathbb{Z}$. We interpret this unexpected dependence on $\alpha$ as a consequence of the non-locality of the field $\psi_\alpha(z)$, which appears to modify the conformal behaviour of the other fields that are present in the theory. It is currently uncertain whether this feature extends to other values of $\beta$, and in particular if it occurs only for values of $\beta$ where the boundary CFT has non-trivial indecomposable modules and Jordan cells. There was also another logarithmic field that appeared in the OPEs: $\omega_{\tau^-}(z)$. Our calculation of the structure constants did not allow us to determine whether it coincides with $\omega_{\tau^+}(z)$, and left open the possibility that the boundary CFT for critical dense polymers has an infinite number of such logarithmic fields with dimension $\Delta = 0$. (It certainly does have an infinity of rank-two Jordan cells, and some of the corresponding logarithmic couplings have been computed in [16].)

We believe there is more to be learned about the field $\psi_\alpha(z)$. Notably, our lattice results have only allowed us to put constraints on ratios of structure constants involving the field $\psi_\alpha(z)$, see (6.23) and (7.10b). This can be traced back to the fact that, in the ratios of conformal correlation functions that we considered, the field $\psi_\alpha(z)$ appears in both the numerator and denominator. For $\beta = 0$, we believe that the fusion of this field with itself is of the form $\psi_\alpha \times \psi_\alpha = \varphi + \bar{\omega} + \dots$ where $\bar{\omega}$ is a logarithmic bulk field with dimension $\Delta = 0$. Indeed, the two-point function $\langle \psi_\alpha(z_1)\psi_\alpha(z_2) \rangle_C$ is non-vanishing, and because $\langle \varphi(z) \rangle_C = 0$, the fusion of $\psi_\alpha$ with itself must include a second field with a vanishing conformal dimension. There is, in fact, a simple geometrical interpretation for the fields $\psi_\alpha(z)$ and $\bar{\omega}(z)$. When two fields $\psi_\alpha(z_1)$ and $\psi_\alpha(z_2)$ are inserted, the loops that encircle only $z_1$ or only $z_2$ have a fugacity $\alpha$, whereas the loops that encircle neither or both have a fugacity $\beta$. On the lattice, one can bring $z_1$ close to $z_2$ until these points are one lattice spacing apart. In this case, there is a unique loop that passes between the two insertion points, and for $\beta = 0$, this is the only loop in the configuration. This loop is space-filling and its fugacity $\alpha$ appears as an overall prefactor. This is consistent with a field $\bar{\omega}$ whose conformal dimension is independent of $\alpha$.

There also remains some light to be shed upon the fields $\mu_{\tau^-}(z)$ and $\mu_{\tau^+}(z)$. Although
the consistency of the lattice and conformal results is quite convincing, we struggle to find a proper justification for treating these two fields in an asymmetric way, with one of them belonging to a Kac module and the other one to a staggered module. Moreover, the structure constant $\kappa_1^{\mu\psi\psi}$ that appears in (7.9) belongs to a three-point function $\langle \mu_{\Sigma}(z_1)\psi_\alpha(z_2)\psi_\alpha(z_3) \rangle_C$ for which we have no lattice interpretation. One possibility is that our identification of $\mu_{\Sigma}(z)$ and $\mu_{\Sigma\Sigma}(z)$ as primary fields is incorrect, and that the correct identification should be made in terms of specific linear combinations of link states inserted in the four marked nodes, instead of only $\Sigma$. One attempt at resolving this issue would be to compute the five-point correlator $\langle \mu_{\Sigma\Sigma}(z_1)\phi_\tau(z_2)\phi_\tau(z_3)\psi_\alpha(z_4)\psi_\alpha(z_5) \rangle_C$ from the lattice and investigate its conformal behaviour.

It is by now well-known that the six-vertex model at $q=i$ and the model of critical dense polymers are intimately related. The model of critical dense polymers has the central charge $c = -2$, whereas the six-vertex model has the central charge $c = 1$. On the torus, the two models share the same modular invariant partition function [60,66], and many questions that arise in one model have a natural interpretation in the other. For instance, assigning a fugacity $\alpha$ to non-contractible loops in the loop model corresponds to inserting the field $\psi_\alpha(z, z^*)$ at $i\infty$. In the vertex model, this corresponds to inserting two electric operators of opposite charges, one at $i\infty$ and another one on the boundary of the cylinder, with a twist line connecting the two. In the related Coulomb gas formalism [62, 66], these electric operators correspond to vertex operators. Their scaling dimensions are to be measured relative to the corresponding groundstate in which non-contractible and contractible loops get the same weight, namely $\alpha = \beta = 0$, and this results in the value (5.36) for $\Delta$. Thus, the insertion of the vertex operators in the six-vertex model changes the central charge from $c = 1$ to $c = -2$ and shifts all the scaling dimensions. The indecomposable structures of the representations, which eventually lead to the logarithmic behaviour of the correlation functions of critical dense polymers, are key features of the conformal description with $c = -2$.

The boundary fields that we studied also have an elegant geometric interpretation in terms of the Q-state Potts model and its Fortuin-Kasteleyn high-temperature expansion. The loop fugacity is related to the number of states in the Potts model via the relation $\beta = \sqrt{Q}$. The Potts spins occupy one half of the lattice (for instance the odd sublattice) and the closed loops draw the contours of the Fortuin-Kasteleyn clusters [67]. On the boundary, free boundary conditions for the Potts spins correspond to having a segment of the boundary decorated with simple half-arcs, $\sigma\sigma...\sigma$, where the red circles mark the positions of the Potts spins. One can also choose a boundary where the Potts spins are fixed to one of the Q values. The Potts model has $S_Q$ symmetry, so in computing the partition function on a domain with this boundary, one can equivalently impose that the boundary spins take the same value, which can be any of the Q spin values, and then divide by Q. In the loop model, the corresponding boundary condition also consists of simple arcs, but with the spins shifted with respect to the half-arcs: $\sigma\sigma...\sigma$. This is usually referred to as *wired boundary conditions*. As a result, the fields $\phi_\tau(z)$ and $\phi_\alpha(z)$ mark transitions between free and wired (or free and fixed, up to a factor of Q) boundary conditions for the Potts spins. Such boundary operators that change the boundary condition from free to wired were previously studied in [68, 69], and found to have the conformal weight $\Delta_{1,2}$. Likewise, $\omega_\tau(z)$ marks a transition between two adjacent large connected clusters of boundary Potts spins. Finally, the fields $\mu_{\Sigma\Sigma}(z)$ and $\mu_{\Sigma\tau}(z)$ correspond to having two spins on the boundary that are in the same cluster, but where all the other adjacent spins are free. Their conformal weight is $\Delta_{1,3}$. These five boundary operators thus appear to be refined versions of the similar operators that insert $d$ adjacent straight defects in the boundary (for $d=1,2$), and whose conformal weights are known [70] to equal $\Delta_{1,d+1}$. The logarithmic nature of the fields $\mu$ and $\omega$ is special to the value $\beta = \sqrt{Q} = 0$, as in this
case there is a *collision* [13] of the values of the conformal weights: $\Delta_{1,3} = \Delta_{1,1}$.

To conclude, we expect that it will be possible to generalise our results to other values of $\beta$ using the algebraic Bethe ansatz, and similarly to cases where the boundary condition at the edge of the cylinder has blobbed half-arcs [71, 72]. We believe that both cases are accessible using the methods of logarithmic conformal field theory. For other values of $\beta$ of the form $2 \cos(\frac{\pi(p'-p)}{p'})$ with $p, p'$ coprime integers, the nature of the logarithmic fields will differ compared to what occurs for $\beta = 0$, as the coincidences of the conformal weights are different. The blobbed case is in particular interesting because it gives access to weights $\Delta_{r,r\pm s}$ for $r \in \mathbb{R}$ [69], thus offering more possibilities for the formation of indecomposable structures.

## Acknowledgments

AMD is an FNRS Postdoctoral Researcher under the project CR28075116. He acknowledges support from the EOS-contract O013018F. JLJ is grateful for support from the European Research Council under the Advanced Grant NuQFT. The authors thank the staff and members of the School of Mathematics and Statistics at the University of Melbourne for their kind hospitality during a simultaneous visit in 2018. The authors thank Yacine Ikhlef, Philippe Ruelle and David Ridout for helpful discussions, and Gilles Parez for proofreading the manuscript.

## References

[1] L. Onsager, *Crystal statistics. I. A two-dimensional model with an order-disorder transition*, Phys. Rev. **65**, 117 (1944), doi:10.1103/PhysRev.65.117.

[2] K. Wilson, *The renormalization group and the $\epsilon$ expansion*, Phys. Rep. **12**, 75 (1974), doi:10.1016/0370-1573(74)90023-4.

[3] A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, *Infinite conformal symmetry in two-dimensional quantum field theory*, Nucl. Phys. B **241**, 333 (1984), doi:10.1016/0550-3213(84)90052-X.

[4] P. Di Francesco, P. Mathieu and D. Sénéchal, *Conformal field theory*, Springer New York (1997), doi:10.1007/978-1-4612-2256-9.

[5] E. K. Sklyanin, L. A. Takhtadzhyan and L. D. Faddeev, *Quantum inverse problem method*. I, Theor. Math. Phys. **40**, 688 (1979), doi:10.1007/BF01018718.

[6] R.J. Baxter. *Exactly solved models in statistical mechanics*, Academic Press (1982).

[7] V. E. Korepin, N. M. Bogoliubov and A. G. Izergin, *Quantum inverse scattering method and correlation functions*, Cambridge University Press (1993), doi:10.1017/CBO9780511628832.

[8] D. Friedan, Z. Qiu and S. Shenker, *Conformal invariance, unitarity and two dimensional critical exponents*, in Mathematical Sciences Research Institute Publications, Springer US (1985), doi:10.1007/978-1-4613-9550-8_21.

[9] V. G. Kac, A. K. Raina and N. Rozhkovskaya, *Bombay lectures on highest weight representations of infinite dimensional lie algebras*, World Scientific, Singapore (2013), doi:10.1142/8882.

[10] K. Iohara and Y. Koga, *Representation theory of the Virasoro algebra*, Springer London (2011), doi:10.1007/978-0-85729-160-8.
[11] L. Rozansky and H. Saleur, *Quantum field theory for the multi-variable Alexander-Conway polynomial*, Nucl. Phys. B **376**, 461 (1992), doi:10.1016/0550-3213(92)90118-U.

[12] V. Gurarie, *Logarithmic operators in conformal field theory*, Nucl. Phys. B **410**, 535 (1993), doi:10.1016/0550-3213(93)90528-W.

[13] J. Cardy, *Logarithmic correlations in quenched random magnets and polymers* (1999), arXiv:cond-mat/9911024.

[14] V. Gurarie and A. W. W. Ludwig, *Conformal algebras of two-dimensional disordered systems*, J. Phys. A: Math. Gen. **35**, L377 (2002), doi:10.1088/0305-4470/35/27/101.

[15] J. Dubail, J. L. Jacobsen and H. Saleur, *Conformal field theory at central charge $c = 0$: A measure of the indecomposability (b) parameters*, Nucl. Phys. B **834**, 399 (2010), doi:10.1016/j.nuclphysb.2010.02.016.

[16] R. Vasseur, J. L. Jacobsen and H. Saleur, *Indecomposability parameters in chiral logarithmic conformal field theory*, Nucl. Phys. B **851**, 314 (2011), doi:10.1016/j.nuclphysb.2011.05.018.

[17] D. Ridout and Y. Saint-Aubin, *Standard modules, induction and the structure of the Temperley-Lieb algebra*, Adv. Theor. Math. Phys. **18**, 957 (2014), doi:10.4310/ATMP2014.v18.n5.a1.

[18] F. Rohsiepe, *On reducible but indecomposable representations of the Virasoro algebra* (1996), arXiv:hep-th/9611160.

[19] M. R. Gaberdie and H. G. Kausch, *Indecomposable fusion products*, Nucl. Phys. B **477**, 293 (1996), doi:10.1016/0550-3213(96)00364-1.

[20] P. Mathieu and D. Ridout, *From percolation to logarithmic conformal field theory*, Phys. Lett. B **657**, 120 (2007), doi:10.1016/j.physletb.2007.10.007.

[21] P. Mathieu and D. Ridout, *Logarithmic $M(2,p)$ minimal models, their logarithmic couplings, and duality*, Nucl. Phys. B **801**, 268 (2008), doi:10.1016/j.nuclphysb.2008.02.017.

[22] K. Kytölä and D. Ridout, *On staggered indecomposable Virasoro modules*, J. Math. Phys. **50**, 123503 (2009), doi:10.1063/1.3191682.

[23] A. Gainutdinov, D. Ridout and I. Runkel, *Logarithmic conformal field theory*, J. Phys. A: Math. Theor. **46**, 490301 (2013), doi:10.1088/1751-8113/46/49/490301.

[24] A. M. Gainutdinov, J. L. Jacobsen, H. Saleur and R. Vasseur, *A physical approach to the classification of indecomposable Virasoro representations from the blob algebra*, Nucl. Phys. B **873**, 614 (2013), doi:10.1016/j.nuclphysb.2013.04.017.

[25] P. A. Pearce and J. Rasmussen, *Solvable critical dense polymers*, J. Stat. Mech. P02015 (2007), doi:10.1088/1742-5468/2007/02/P02015.

[26] P. A. Pearce, J. Rasmussen and J.-B. Zuber, *Logarithmic minimal models*, J. Stat. Mech. P11017 (2006), doi:10.1088/1742-5468/2006/11/P11017.

[27] H. G. Kausch, *Curiosities at $c=-2$* (1995), arXiv:hep-th/9510149.

[28] M. R. Gaberdiel and H. G. Kausch, *A local logarithmic conformal field theory*, Nucl. Phys. B **538**, 631 (1999), doi:10.1016/S0550-3213(98)00701-9.
[29] H. G. Kausch, *Symplectic fermions*, Nucl. Phys. B **583**, 513 (2000), doi:10.1016/S0550-3213(00)00295-9.

[30] S. Kawai and J. F. Wheater, *Modular transformation and boundary states in logarithmic conformal field theory*, Phys. Lett. B **508**, 203 (2001), doi:10.1016/S0370-2693(01)00503-2.

[31] A. Bredthauer and M. Flohr, *Boundary states in $c=2$ logarithmic conformal field theory*, Nucl. Phys. B **639**, 450 (2002), doi:10.1016/S0550-3213(02)00466-2.

[32] J. J. Graham and G. I. Lehrer, *Cellular algebras*, Invent. Math. **123**, 1 (1996), doi:10.1007/BF01232365.

[33] G. Piroux and P. Ruelle, *Pre-logarithmic and logarithmic fields in a sandpile model*, J. Stat. Mech.: Theor. Exp. P10005 (2004), doi:10.1088/1742-5468/2004/10/P10005.

[34] M. Jeng, G. Piroux and P. Ruelle, *Height variables in the Abelian sandpile model: scaling fields and correlations*, J. Stat. Mech. P10015 (2006), doi:10.1088/1742-5468/2006/10/P10015.

[35] P. Ruelle, *Logarithmic conformal invariance in the abelian sandpile model*, J. Phys. A: Math. Theor. **46**, 494014 (2013), doi:10.1088/1751-8113/46/49/494014.

[36] E. V. Ivashkevich, *Correlation functions of dense polymers and $c=-2$ conformal field theory*, J. Phys. A: Math. Gen. **32**, 1691 (1999), doi:10.1088/0305-4470/32/9/015.

[37] E. V. Ivashkevich and C.-K. Hu, *Exact multileg correlation functions for the dense phase of branching polymers in two dimensions*, Phys. Rev. E **71**, 015104 (2005), doi:10.1103/PhysRevE.71.015104.

[38] A. Morin-Duchesne and J. L. Jacobsen, *Two-point boundary correlation functions of dense loop models*, SciPost Phys. **4**, 034 (2018), doi:10.21468/SciPostPhys.4.6.034.

[39] R. Vasseur, J. L. Jacobsen and H. Saleur, *Logarithmic observables in critical percolation*, J. Stat. Mech. L07001 (2012), doi:10.1088/1742-5468/2012/07/L07001.

[40] G. Gori and J. Viti, *Four-point boundary connectivities in critical two-dimensional percolation from conformal invariance*, J. High Energ. Phys. **12**, 131 (2018), doi:10.1007/JHEP12(2018)131.

[41] X. Tan, R. Couvreur, Y. Deng and J. L. Jacobsen, *Observation of nonscalar and logarithmic correlations in 2D and 3D percolation*, Phys. Rev. E **99**, 050103 (2019), doi:10.1103/PhysRevE.99.050103.

[42] R. Vasseur and J. L. Jacobsen, *Operator content of the critical Potts model in $d$ dimensions and logarithmic correlations*, Nucl. Phys. B **880**, 435 (2014), doi:10.1016/j.nuclphysb.2014.01.013.

[43] R. Couvreur, J. L. Jacobsen and R. Vasseur, *Non-scalar operators for the Potts model in arbitrary dimension*, J. Phys. A: Math. Theor. **50**, 474001 (2017), doi:10.1088/1751-8121/aa7f32.

[44] G. Gori and J. Viti, *Exact logarithmic four-point functions in the critical two-dimensional Ising model*, Phys. Rev. Lett. **119**, 191601 (2017), doi:10.1103/PhysRevLett.119.191601.

[45] M. A. I. Flohr, *Bits and pieces in logarithmic conformal field theory*, Int. J. Mod. Phys. A **18**, 4497 (2003), doi:10.1142/S0217751X03016859.
[46] J. L. Cardy, Boundary conditions, fusion rules and the Verlinde formula, Nucl. Phys. B **324**, 581 (1989), doi:10.1016/0550-3213(89)90521-X.

[47] J. Cardy, Boundary conformal field theory, in Encyclopedia of Mathematical Physics, Elsevier (2006), arXiv:hep-th/0411189.

[48] J. Dubail, J. L. Jacobsen and H. Saleur, Conformal two-boundary loop model on the annulus, Nucl. Phys. B **813**, 430 (2009), doi:10.1016/j.nuclphysb.2008.12.023.

[49] J. de Gier, J. L. Jacobsen and A. Ponsaing, Finite-size corrections for universal boundary entropy in bond percolation, SciPost Phys. **1**, 012 (2016), doi:10.21468/SciPostPhys.1.2.012.

[50] A. M. Gainutdinov, N. Read and H. Saleur, Continuum limit and symmetries of the periodic spin chain, Nucl. Phys. B **871**, 245 (2013), doi:10.1016/j.nuclphysb.2013.01.018.

[51] A. M. Gainutdinov, N. Read and H. Saleur, Bimodule structure in the periodic $g\ell(1|1)$ spin chain, Nucl. Phys. B **871**, 289 (2013), doi:10.1016/j.nuclphysb.2013.02.017.

[52] A. Morin-Duchesne and Y. Saint-Aubin, Jordan cells of periodic loop models, J. Phys. A: Math. Theor. **46**, 494013 (2013), doi:10.1088/1751-8113/46/49/494013.

[53] A. M. Gainutdinov, N. Read, H. Saleur and R. Vasseur, The periodic $s\ell(2|1)$ alternating spin chain and its continuum limit as a bulk logarithmic conformal field theory at $c = 0$, J. High Energ. Phys. **05**, 114 (2015), doi:10.1007/JHEP05(2015)114.

[54] D. Levy, Algebraic structure of translation-invariant spin-1/2 xxz and q-Potts quantum chains, Phys. Rev. Lett. **67**, 1971 (1991), doi:10.1103/PhysRevLett.67.1971.

[55] P. Martin and H. Saleur, On an algebraic approach to higher dimensional statistical mechanics, Commun. Math. Phys. **158**, 155 (1993), doi:10.1007/BF02097236.

[56] J.J. Graham and G.I. Lehrer, The representation theory of affine Temperley-Lieb algebras, Enseign. Math. **44**, 173 (1998).

[57] R.M. Green, On representations of affine Temperley-Lieb algebras, Can. Math. Soc. Conf. Proc. **24**, 245 (1998).

[58] K. Erdmann and R. M. Green, On representations of affine Temperley-Lieb algebras, II (1998), arXiv:math/9811017.

[59] P. A. Pearce, J. Rasmussen and S. P. Villani, Solvable critical dense polymers on the cylinder, J. Stat. Mech. P02010 (2010), doi:10.1088/1742-5468/2010/02/P02010.

[60] A. Morin-Duchesne, P. A. Pearce and J. Rasmussen, Modular invariant partition function of critical dense polymers, Nucl. Phys. B **874**, 312 (2013), doi:10.1016/j.nuclphysb.2013.05.016.

[61] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, NIST Handbook of Mathematical Functions, Cambridge University Press (2010), doi:10.1111/j.1751-5823.2011.00134_18.x.

[62] B. Nienhuis, Critical behavior of two-dimensional spin models and charge asymmetry in the Coulomb gas, J. Stat. Phys. **34**, 731 (1984), doi:10.1007/BF01009437.
[63] V. Pasquier and H. Saleur, *Common structures between finite systems and conformal field theories through quantum groups*, Nucl. Phys. B **330**, 523 (1990), doi:10.1016/0550-3213(90)90122-T.

[64] A. B. Zamolodchikov and V. A. Fateev, *Disorder fields in two-dimensional conformal quantum field theory and N = 2 extended supersymmetry*. Sov. Phys. JETP **63**, 913 (1986).

[65] Vl. S. Dotsenko and V. A. Fateev, *Conformal algebra and multipoint correlation functions in 2D statistical models*, Nucl. Phys. B **240**, 312 (1984), doi:10.1016/0550-3213(84)90269-4.

[66] P. Di Francesco, H. Saleur and J. B. Zuber, *Relations between the Coulomb gas picture and conformal invariance of two-dimensional critical models*, J. Stat. Phys. **49**, 57 (1987), doi:10.1007/BF01009954.

[67] R. J. Baxter, S. B. Kelland and F. Y. Wu, *Equivalence of the Potts model or Whitney polynomial with an ice-type model*, J. Phys. A: Math. Gen. **9**, 397 (1976), doi:10.1088/0305-4470/9/3/009.

[68] J. L. Cardy, *Critical percolation in finite geometries*, J. Phys. A: Math. Gen. **25**, L201 (1992), doi:10.1088/0305-4470/25/4/009.

[69] J. L. Jacobsen and H. Saleur, *Conformal boundary loop models*, Nucl. Phys. B **788**, 137 (2008), doi:10.1016/j.nuclphysb.2007.06.029.

[70] H. Saleur and M. Bauer, *On some relations between local height probabilities and conformal invariance*, Nucl. Phys. B **320**, 591 (1989), doi:10.1016/0550-3213(89)90014-X.

[71] P. Martin and H. Saleur, *The blob algebra and the periodic Temperley-Lieb algebra*, Lett. Math. Phys. **30**, 189 (1994), doi:10.1007/BF00805852.

[72] P. P. Martin and D. Woodcock, *On the structure of the blob algebra*, J. Algebra **225**, 957 (2000), doi:10.1006/jabr.1999.7948.