Source author record

Jesper Lykke Jacobsen

Jesper Lykke Jacobsen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

72works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

72 published item(s)

preprint2022arXiv

Geometric Algebra and Algebraic Geometry of Loop and Potts Models

We uncover a connection between two seemingly separate subjects in integrable models: the representation theory of the affine Temperley-Lieb algebra, and the algebraic structure of solutions to the Bethe equations of the XXZ spin chain. We study the solution of Bethe equations analytically by computational algebraic geometry, and find that the solution space encodes rich information about the representation theory of Temperley-Lieb algebra. Using these connections, we compute the partition function of the completely-packed loop model and of the closely related random-cluster Potts model, on medium-size lattices with toroidal boundary conditions, by two quite different methods. We consider the partial thermodynamic limit of infinitely long tori and analyze the corresponding condensation curves of the zeros of the partition functions. Two components of these curves are obtained analytically in the full thermodynamic limit.

preprint2022arXiv

Global symmetry and conformal bootstrap in the two-dimensional $O(n)$ model

We define the two-dimensional $O(n)$ conformal field theory as a theory that includes the critical dilute and dense $O(n)$ models as special cases, and depends analytically on the central charge. For generic values of $n\in\mathbb{C}$, we write a conjecture for the decomposition of the spectrum into irreducible representations of $O(n)$. We then explain how to numerically bootstrap arbitrary four-point functions of primary fields in the presence of the global $O(n)$ symmetry. We determine the needed conformal blocks, including logarithmic blocks, including in singular cases. We argue that $O(n)$ representation theory provides upper bounds on the number of solutions of crossing symmetry for any given four-point function. We study some of the simplest correlation functions in detail, and determine a few fusion rules. We count the solutions of crossing symmetry for the $30$ simplest four-point functions. The number of solutions varies from $2$ to $6$, and saturates the bound from $O(n)$ representation theory in $21$ out of $30$ cases.

preprint2021arXiv

Nested Closed Paths in Two-Dimensional Percolation

For two-dimensional percolation on a domain with the topology of a disc, we introduce a nested-path operator (NP) and thus a continuous family of one-point functions $W_k \equiv \langle \mathcal{R} \cdot k^\ell \rangle $, where $\ell$ is the number of independent nested closed paths surrounding the center, $k$ is a path fugacity, and $\mathcal{R}$ projects on configurations having a cluster connecting the center to the boundary. At criticality, we observe a power-law scaling $W_k \sim L^{X_{\rm NP}}$, with $L$ the linear system size, and we determine the exponent $X_{\rm NP}$ as a function of $k$. On the basis of our numerical results, we conjecture an analytical formula, $X_{\rm NP} (k) = \frac{3}{4}ϕ^2 -\frac{5}{48}ϕ^2/ (ϕ^2-\frac{2}{3})$ where $k = 2 \cos(πϕ)$, which reproduces the exact results for $k=0,1$ and agrees with the high-precision estimate of $X_{\rm NP}$ for other $k$ values. In addition, we observe that $W_2(L)=1$ for site percolation on the triangular lattice with any size $L$, and we prove this identity for all self-matching lattices.

preprint2021arXiv

The action of the Virasoro algebra in quantum spin chains. I. The non-rational case

We investigate the action of discretized Virasoro generators, built out of generators of the lattice Temperley-Lieb algebra ("Koo-Saleur generators"[arXiv:hep-th/9312156]), in the critical XXZ quantum spin chain. We explore the structure of the continuum-limit Virasoro modules at generic central charge for the XXZ vertex model, paralleling [arXiv:2007.11539] for the loop model. We find again indecomposable modules, but this time not logarithmic ones. The limit of the Temperley-Lieb modules $\mathcal{W}_{j,1}$ for $j\neq 0$ contains pairs of "conjugate states" with conformal weights $(h_{r,s},h_{r,-s})$ and $(h_{r,-s},h_{r,s})$ that give rise to dual structures: Verma or co-Verma modules. The limit of $\mathcal{W}_{0,\mathfrak{q}^{\pm2}}$ contains diagonal fields $(h_{r,1},h_{r,1})$ and gives rise to either only Verma or only co-Verma modules, depending on the sign of the exponent in $\mathfrak{q}^{\pm 2}$. In order to obtain matrix elements of Koo-Saleur generators at large system size $N$ we use Bethe ansatz and Quantum Inverse Scattering methods, computing the form factors for relevant combinations of three neighbouring spin operators. Relations between form factors ensure that the above duality exists already at the lattice level. We also study in which sense Koo-Saleur generators converge to Virasoro generators. We consider convergence in the weak sense, investigating whether the commutator of limits is the same as the limit of the commutator? We find that it coincides only up to the central term. As a side result we compute the ground-state expectation value of two neighbouring Temperley-Lieb generators in the XXZ spin chain.

preprint2021arXiv

The action of the Virasoro algebra in the two-dimensional Potts and loop models at generic $Q$

The spectrum of conformal weights for the CFT describing the two-dimensional critical $Q$-state Potts model (or its close cousin, the dense loop model) has been known for more than 30 years. However, the exact nature of the corresponding $\hbox{Vir}\otimes\overline{\hbox{Vir}}$ representations has remained unknown up to now. Here, we solve the problem for generic values of $Q$. This is achieved by a mixture of different techniques: a careful study of "Koo--Saleur generators" [arXiv:hep-th/9312156], combined with measurements of four-point amplitudes, on the numerical side, and OPEs and the four-point amplitudes recently determined using the "interchiral conformal bootstrap" in [arXiv:2005.07258] on the analytical side. We find that null-descendants of diagonal fields having weights $(h_{r,1},h_{r,1})$ (with $r\in \mathbb{N}^*$) are truly zero, so these fields come with simple $\hbox{Vir}\otimes\overline{\hbox{Vir}}$ ("Kac") modules. Meanwhile, fields with weights $(h_{r,s},h_{r,-s})$ and $(h_{r,-s},h_{r,s})$ (with $r,s\in\mathbb{N}^*$) come in indecomposable but not fully reducible representations mixing four simple $\hbox{Vir}\otimes\overline{\hbox{Vir}}$ modules with a familiar "diamond" shape. The "top" and "bottom" fields in these diamonds have weights $(h_{r,-s},h_{r,-s})$, and form a two-dimensional Jordan cell for $L_0$ and $\bar{L}_0$. This establishes, among other things, that the Potts-model CFT is logarithmic for $Q$ generic. Unlike the case of non-generic (root of unity) values of $Q$, these indecomposable structures are not present in finite size, but we can nevertheless show from the numerical study of the lattice model how the rank-two Jordan cells build up in the infinite-size limit.

preprint2020arXiv

$N$-cluster correlations in four- and five-dimensional percolation

We study $N$-cluster correlation functions in four- and five-dimensional (4D, 5D) bond percolation by extensive Monte Carlo simulation. We reformulate the transfer Monte Carlo algorithm for percolation [Phys. Rev. E {\bf 72}, 016126 (2005)] using the disjoint-set data structure, and simulate a cylindrical geometry $L^{d-1}\times \infty$, with the linear size up to $L=512$ for 4D and $128$ for 5D. We determine with a high precision all possible $N$-cluster exponents, for $N \! =\!2$ and $3$, and the universal amplitude for a logarithmic correlation function. From the symmetric correlator with $N \! = \!2$, we obtain the correlation-length critical exponent as $1/ν\! =\! 1.4610(12)$ for 4D and $1/ν\! =\! 1.737 (2)$ for 5D, significantly improving over the existing results. Estimates for the other exponents and the universal logarithmic amplitude have not been reported before to our knowledge. Our work demonstrates the validity of logarithmic conformal field theory and adds to the growing knowledge for high-dimensional percolation.

preprint2020arXiv

Analytic continuation of Bethe energies and application to the thermodynamic limit of the $SL(2,\mathbb{C})$ non-compact spin chains

We consider the problem of analytically continuing energies computed with the Bethe ansatz, as posed by the study of non-compact integrable spin chains. By introducing an imaginary extensive twist in the Bethe equations, we show that one can expand the analytic continuation of energies in the scaling limit around another 'pseudo-vacuum' sitting at a negative number of Bethe roots, in the same way as around the usual pseudo-vacuum. We show that this method can be used to compute the energy levels of some states of the $SL(2,\mathbb{C})$ integrable spin chain in the infinite-volume limit, and as a proof of principle recover the ground-state value previously obtained in [1] (for the case of spins $s=0, \bar{s}=-1$) by extrapolating results in small sizes. These results represent, as far as we know, the first (partial) description of the spectrum of $SL(2,\mathbb{C})$ non-compact spin chains in the thermodynamic limit.

preprint2020arXiv

Critical percolation on the kagome hypergraph

We study the percolation critical surface of the kagome lattice in which each triangle is allowed an arbitrary connectivity. Using the method of critical polynomials, we find points along this critical surface to high precision. This kagome hypergraph contains many unsolved problems as special cases, including bond percolation on the kagome and $(3,12^2)$ lattices, and site percolation on the hexagonal, or honeycomb, lattice, as well as a single point for which there is an exact solution. We are able to compute enough points along the critical surface to find a very accurate fit, essentially a Taylor series about the exact point, that allows estimations of the critical point of any system that lies on the surface to precision rivaling Monte Carlo and traditional techniques of similar accuracy. We find also that this system sheds light on some of the surprising aspects of the method of critical polynomials, such as why it is so accurate for certain problems, like the kagome and $(3,12^2)$ lattices. The bond percolation critical points of these lattices can be found to 17 and 18 digits, respectively, because they are in close proximity, in a sense that can be made quantitative, to the exact point on the critical surface. We also discuss in detail a parallel implementation of the method which we use here for a few calculations.

preprint2020arXiv

Cylinder partition function of the 6-vertex model from algebraic geometry

We compute the exact partition function of the isotropic 6-vertex model on a cylinder geometry with free boundary conditions, for lattices of intermediate size, using Bethe ansatz and algebraic geometry. We perform the computations in both the open and closed channels. We also consider the partial thermodynamic limits, whereby in the open (closed) channel, the open (closed) direction is kept small while the other direction becomes large. We compute the zeros of the partition function in the two partial thermodynamic limits, and compare with the condensation curves.

preprint2020arXiv

Four-point geometrical correlation functions in the two-dimensional $Q$-state Potts model: connections with the RSOS models

The "bootstrap determination" of the geometrical correlation functions in the two-dimensional Potts model proposed in a paper [arXiv:1607.07224] was later shown in [arXiv:1809.02191] to be incorrect, the actual spectrum of the model being considerably more complex than initially conjectured. We provide in this paper a geometrical interpretation of the four-point functions built in [arXiv:1607.07224], and explain why the results obtained by these authors, albeit incorrect, appeared so close to those of their numerical simulations of the Potts model. Our strategy is based on a cluster expansion of correlation functions in RSOS minimal models, and a subsequent numerical and algebraic analysis of the corresponding $s$-channel spectrum, in full analogy with our early work on the Potts model [arXiv:1809.02191]. Remarkable properties of the lattice amplitudes are uncovered, which explain in particular the truncation of the spectrum of [arXiv:1809.02191] to the much simpler one of the RSOS models, and which will be used in a forthcoming paper to finally determine the geometric four-point functions of the Potts model itself.

preprint2020arXiv

Integrable boundary conditions in the antiferromagnetic Potts model

We present an exact mapping between the staggered six-vertex model and an integrable model constructed from the twisted affine $D_2^2$ Lie algebra. Using the known relations between the staggered six-vertex model and the antiferromagnetic Potts model, this mapping allows us to study the latter model using tools from integrability. We show that there is a simple interpretation of one of the known K-matrices of the $D_2^2$ model in terms of Temperley-Lieb algebra generators, and use this to present an integrable Hamiltonian that turns out to be in the same universality class as the antiferromagnetic Potts model with free boundary conditions. The intriguing degeneracies in the spectrum observed in related works (arXiv:nlin/0002050 and arXiv:1707.09260) are discussed.

preprint2020arXiv

Lattice regularisation of a non-compact boundary conformal field theory

Non-compact Conformal Field Theories (CFTs) are central to several aspects of string theory and condensed matter physics. They are characterised, in particular, by the appearance of a continuum of conformal dimensions. Surprisingly, such CFTs have been identified as the continuum limits of lattice models with a finite number of degrees of freedom per site. However, results have so far been restricted to the case of periodic boundary conditions, precluding the exploration via lattice models of aspects of non-compact boundary CFTs and the corresponding D-brane constructions. The present paper follows a series of previous works on a $\mathbb{Z}_2$-staggered XXZ spin chain, whose continuum limit is known to be a non-compact CFT related with the Euclidian black hole sigma model. By using the relationship of this spin chain with an integrable $D_2^2$ vertex model, we here identify integrable boundary conditions that lead to a continuous spectrum of boundary exponents, and thus correspond to non-compact branes. In the context of the Potts model on a square lattice, they correspond to wired boundary conditions at the physical antiferromagnetic critical point. The relations with the boundary parafermion theories are discussed as well. We are also able to identify a boundary renormalisation group flow from the non-compact boundary conditions to the previously studied compact ones.

preprint2019arXiv

Bond percolation thresholds on Archimedean lattices from critical polynomial roots

We present percolation thresholds calculated numerically with the eigenvalue formulation of the method of critical polynomials; developed in the last few years, it has already proven to be orders of magnitude more accurate than traditional techniques. Here we report the result of large parallel calculations to produce what we believe may become the reference values of bond percolation thresholds on the Archimedean lattices for years to come. For example, for the kagome lattice we find $p_{\rm c}=0.524\,404\,999\,167\,448\,20 (1)$, whereas the best estimate using standard techniques is $p_{\rm c}=0.524\,404\,99(2)$. We further provide strong evidence that there are two classes of lattices: one for which the first three scaling exponents characterizing the finite-size corrections to $p_{\rm c}$ are $Δ=6,7,8$, and another for which $Δ=4,6,8$. We discuss the open questions related to the method, such as the full scaling law, as well as its potential for determining critical points of other models.

preprint2019arXiv

Conformally invariant boundary conditions in the antiferromagnetic Potts model and the $SL(2,\mathbb{R})/U(1)$ sigma model

We initiate a study of the boundary version of the square-lattice $Q$-state Potts antiferromagnet, with $Q \in [0,4]$ real, motivated by the fact that the continuum limit of the corresponding bulk model is a non-compact CFT, closely related with the $SL(2,\mathbb{R})_k/U(1)$ Euclidian black-hole coset model. While various types of conformal boundary conditions (discrete and continuous branes) have been formally identified for the the $SL(2,\mathbb{R})_k/U(1)$ coset CFT, we are only able in this work to identify conformal boundary conditions (CBC) leading to a discrete boundary spectrum. The conformal boundary conditions (CBC) we find are of two types. The first is free boundary Potts spins, for which we confirm an old conjecture for the generating functions of conformal levels, and show them to be related to characters in a non-linear deformation of the $W_\infty$ algebra. The second type of CBC - which corresponds to restricting the values of the Potts spins to a subset of size $Q_1$, or its complement of size $Q-Q_1$, at alternating sites along the boundary - is new, and turns out to be conformal in the antiferromagnetic case only. Using algebraic and numerical techniques, we show that the corresponding spectrum generating functions produce all the characters of discrete representations for the coset CFT. The normalizability bounds of the associated discrete states in the coset CFT are found to have a simple interpretation in terms of boundary phase transitions in the lattice model. For $\sqrt{Q} = 2 \cos \fracπ{k}$, with $k \ge 3$ integer, we show also how our boundary conditions can be reformulated in terms of a RSOS height model. The spectrum generating functions are then identified with string functions of the compact $SU(2)_{k-2} / U(1)$ parafermion theory (with symmetry $Z_{k-2}$). The new alt conditions are needed to cover all the string functions.

preprint2016arXiv

Finite-size corrections for universal boundary entropy in bond percolation

We compute the boundary entropy for bond percolation on the square lattice in the presence of a boundary loop weight, and prove explicit and exact expressions on a strip and on a cylinder of size $L$. For the cylinder we provide a rigorous asymptotic analysis which allows for the computation of finite-size corrections to arbitrary order. For the strip we provide exact expressions that have been verified using high-precision numerical analysis. Our rigorous and exact results corroborate an argument based on conformal field theory, in particular concerning universal logarithmic corrections for the case of the strip due to the presence of corners in the geometry. We furthermore observe a crossover at a special value of the boundary loop weight.

preprint2016arXiv

On the growth constant for square-lattice self-avoiding walks

The growth constant for two-dimensional self-avoiding walks on the honeycomb lattice was conjectured by Nienhuis in 1982, and since that time the corresponding results for the square and triangular lattices have been sought. For the square lattice, a possible conjecture was advanced by one of us (AJG) more than 20 years ago, based on the six significant digit estimate available at the time. This estimate has improved by a further six digits over the intervening decades, and the conjectured value continued to agree with the increasingly precise estimates. We discuss the three most successful methods for estimating the growth constant, including the most recently developed Topological Transfer-Matrix method, due to another of us (JLJ). We show this to be the most computationally efficient of the three methods, and by parallelising the algorithm we have estimated the growth constant significantly more precisely, incidentally ruling out the conjecture, which fails in the 12th digit. Our new estimate of the growth constant is $$μ(\mathrm{square}) = 2.63815853032790\, (3).$$

preprint2016arXiv

Q-colourings of the triangular lattice: Exact exponents and conformal field theory

We revisit the problem of Q-colourings of the triangular lattice using a mapping onto an integrable spin-one model, which can be solved exactly using Bethe Ansatz techniques. In particular we focus on the low-energy excitations above the eigenlevel g_2, which was shown by Baxter to dominate the transfer matrix spectrum in the Fortuin-Kasteleyn (chromatic polynomial) representation for Q_0 <= Q <= 4, where Q_0 = 3.819671... We argue that g_2 and its scaling levels define a conformally invariant theory, the so-called regime IV, which provides the actual description of the (analytically continued) colouring problem within a much wider range, namely 2 < Q <= 4. The corresponding conformal field theory is identified and the exact critical exponents are derived. We discuss their implications for the phase diagram of the antiferromagnetic triangular-lattice Potts model at non-zero temperature. Finally, we relate our results to recent observations in the field of spin-one anyonic chains.

preprint2016arXiv

The continuum limit of $a_{N-1}^{(2)}$ spin chains

Building on our previous work for $a_2^{(2)}$ and $a_3^{(2)}$ we explore systematically the continuum limit of gapless $a_{N-1}^{(2)}$ vertex models and spin chains. We find the existence of three possible regimes. Regimes I and II for $a_{2n-1}^{(2)}$ are related with $a_{2n-1}^{(2)}$ Toda, and described by $n$ compact bosons. Regime I for $a_{2n}^{(2)}$ is related with $a_{2n}^{(2)}$ Toda and involves $n$ compact bosons, while regime II is related instead with $B^{(1)}(0,n)$ super Toda, and involves in addition a single Majorana fermion. The most interesting is regime III, where {\sl non-compact} degrees of freedom appear, generalising the emergence of the Euclidean black hole CFT in the $a_{2}^{(2)}$ case. For $a_{2n}^{(2)}$ we find a continuum limit made of $n$ compact and $n$ non-compact bosons, while for $a_{2n-1}^{(2)}$ we find $n$ compact and $n-1$ non-compact bosons. We also find deep relations between $a_{N-1}^{(2)}$ in regime III and the gauged WZW models $SO(N)/SO(N-1)$.

preprint2015arXiv

A new look at the collapse of two-dimensional polymers

We study the collapse of two-dimensional polymers, via an O($n$) model on the square lattice that allows for dilution, bending rigidity and short-range monomer attractions. This model contains two candidates for the theta point, $Θ_{\rm BN}$ and $Θ_{\rm DS}$, both exactly solvable. The relative stability of these points, and the question of which one describes the `generic' theta point, have been the source of a long-standing debate. Moreover, the analytically predicted exponents of $Θ_{\rm BN}$ have never been convincingly observed in numerical simulations. In the present paper, we shed a new light on this confusing situation. We show in particular that the continuum limit of $Θ_{\rm BN}$ is an unusual conformal field theory, made in fact of a simple dense polymer decorated with {\sl non-compact degrees of freedom}. This implies in particular that the critical exponents take continuous rather than discrete values, and that corrections to scaling lead to an unusual integral form. Furthermore, discrete states may emerge from the continuum, but the latter are only normalizable---and hence observable---for appropriate values of the model's parameters. We check these findings numerically. We also probe the non-compact degrees of freedom in various ways, and establish that they are related to fluctuations of the density of monomers. Finally, we construct a field theoretic model of the vicinity of $Θ_{\rm BN}$ and examine the flow along the multicritical line between $Θ_{\rm BN}$ and $Θ_{\rm DS}$.

preprint2015arXiv

Critical points of Potts and O($N$) models from eigenvalue identities in periodic Temperley-Lieb algebras

In previous work with Scullard, we defined a graph polynomial P_B(q,T) that gives access to the critical temperature T_c of the q-state Potts model on a general two-dimensional lattice L. It depends on a basis B, containing n x m unit cells of L, and the relevant root of P_B(q,T) was observed to converge quickly to T_c in the limit n,m to infinity. Moreover, in exactly solvable cases there is no finite-size dependence at all. We reformulate this method as an eigenvalue problem within the periodic Temperley-Lieb algebra. This corresponds to taking m to infinity first, so the bases B are semi-infinite cylinders of circumference n. The limit implies faster convergence in n, while maintaining the n-independence in exactly solvable cases. In this setup, T_c(n) is determined by equating the largest eigenvalues of two topologically distinct sectors of the transfer matrix. Crucially, these two sectors determine the same critical exponent in the continuum limit, and the observed fast convergence is thus corroborated by results of conformal field theory. We obtain similar results for the dense and dilute phases of the O(N) loop model, using now a transfer matrix within the dilute periodic Temperley-Lieb algebra. The eigenvalue formulation allows us to double the size n for which T_c(n) can be obtained, using the same computational effort. We study in details three significant cases: (i) bond percolation on the kagome lattice, up to n = 14; (ii) site percolation on the square lattice, to n = 21; and (iii) self-avoiding polygons on the square lattice, to n = 19. Convergence properties of T_c(n) and extrapolation schemes are studied in details for the first two cases. This leads to rather accurate values for the percolation thresholds: p_c = 0.524404999167439(4) for bond percolation on the kagome lattice, and p_c = 0.59274605079210(2) for site percolation on the square lattice.

preprint2015arXiv

Dilute oriented loop models

We study a model of dilute oriented loops on the square lattice, where each loop is compatible with a fixed, alternating orientation of the lattice edges. This implies that loop strands are not allowed to go straight at vertices, and results in an enhancement of the usual O(n) symmetry to U(n). The corresponding transfer matrix acts on a number of representations (standard modules) that grows exponentially with the system size. We derive their dimension and those of the centraliser by both combinatorial and algebraic techniques. A mapping onto a field theory permits us to identify the conformal field theory governing the critical range, $n \le 1$. We establish the phase diagram and the critical exponents of low-energy excitations. For generic n, there is a critical line in the universality class of the dilute O(2n) model, terminating in an SU(n+1) point. The case n=1 maps onto the critical line of the six-vertex model, along which exponents vary continuously.

preprint2015arXiv

High-precision percolation thresholds and Potts-model critical manifolds from graph polynomials

The critical curves of the q-state Potts model can be determined exactly for regular two-dimensional lattices G that are of the three-terminal type. Jacobsen and Scullard have defined a graph polynomial P_B(q,v) that gives access to the critical manifold for general lattices. It depends on a finite repeating part of the lattice, called the basis B, and its real roots in the temperature variable v = e^K - 1 provide increasingly accurate approximations to the critical manifolds upon increasing the size of B. These authors computed P_B(q,v) for large bases (up to 243 edges), obtaining determinations of the ferromagnetic critical point v_c > 0 for the (4,8^2), kagome, and (3,12^2) lattices to a precision (of the order 10^{-8}) slightly superior to that of the best available Monte Carlo simulations. In this paper we describe a more efficient transfer matrix approach to the computation of P_B(q,v) that relies on a formulation within the periodic Temperley-Lieb algebra. This makes possible computations for substantially larger bases (up to 882 edges), and the precision on v_c is hence taken to the range 10^{-13}. We further show that a large variety of regular lattices can be cast in a form suitable for this approach. This includes all Archimedean lattices, their duals and their medials. For all these lattices we tabulate high-precision estimates of the bond percolation thresholds p_c and Potts critical points v_c. We also trace and discuss the full Potts critical manifold in the (q,v) plane, paying special attention to the antiferromagnetic region v < 0. Finally, we adapt the technique to site percolation as well, and compute the polynomials P_B(p) for certain Archimedean and dual lattices (those having only cubic and quartic vertices), using very large bases (up to 243 vertices). This produces the site percolation thresholds p_c to a precision of the order 10^{-9}.

preprint2015arXiv

Potts-model critical manifolds revisited

We compute the critical polymials for the q-state Potts model on all Archimedean lattices, using a parallel implementation of the algorithm of (Jacobsen, J. Phys. A: Math. Theor. 47 135001) that gives us access to larger sizes than previously possible. The exact polynomials are computed for bases of size $6 \times 6$ unit cells, and the root in the temperature variable $v=e^K-1$ is determined numerically at $q=1$ for bases of size $8 \times 8$. This leads to improved results for bond percolation thresholds, and for the Potts-model critical manifolds in the real $(q,v)$ plane. In the two most favourable cases, we find now the kagome-lattice threshold to eleven digits and that of the $(3,12^2)$ lattice to thirteen. Our critical manifolds reveal many interesting features in the antiferromagnetic region of the Potts model, and determine accurately the extent of the Berker-Kadanoff phase for the lattices studied.

preprint2015arXiv

Recursive Percolation

We introduce a simple lattice model in which percolation is constructed on top of critical percolation clusters, and show that it can be repeated recursively any number $n$ of generations. In two dimensions, we determine the percolation thresholds up to $n=5$. The corresponding critical clusters become more and more compact as $n$ increases, and define universal scaling functions of the standard two-dimensional form and critical exponents that are distinct for any $n$. This family of exponents differs from any previously known universality class, and cannot be accommodated by existing analytical methods. We confirm that recursive percolation is well defined also in three dimensions.

preprint2015arXiv

Three-point functions in c <= 1 Liouville theory and conformal loop ensembles

The possibility of extending the Liouville Conformal Field Theory from values of the central charge $c \geq 25$ to $c \leq 1$ has been debated for many years in condensed matter physics as well as in string theory. It was only recently proven that such an extension -- involving a real spectrum of critical exponents as well as an analytic continuation of the DOZZ formula for three-point couplings -- does give rise to a consistent theory. We show in this Letter that this theory can be interpreted in terms of microscopic loop models. We introduce in particular a family of geometrical operators, and, using an efficient algorithm to compute three-point functions from the lattice, we show that their operator algebra corresponds exactly to that of vertex operators $V_{\hatα}$ in $c \leq 1$ Liouville. We interpret geometrically the limit $\hatα \to 0$ of $V_{\hatα}$ and explain why it is not the identity operator (despite having conformal weight $Δ=0$).

preprint2014arXiv

High-precision phase diagram of spin glasses from duality analysis with real-space renormalization and graph polynomials

We propose a duality analysis for obtaining the critical manifold of two-dimensional spin glasses. Our method is based on the computation of quenched free energies with periodic and twisted periodic boundary conditions on a finite basis. The precision can be systematically improved by increasing the size of the basis, leading to very fast convergence towards the thermodynamic limit. We apply the method to obtain the phase diagrams of the random-bond Ising model and $q$-state Potts gauge glasses. In the Ising case, the Nishimori point is found at $p_N = 0.10929 \pm 0.00002$, in agreement with and improving on the precision of existing numerical estimations. Similar precision is found throughout the high-temperature part of the phase diagram. Finite-size effects are larger in the low-temperature region, but our results are in qualitative agreement with the known features of the phase diagram. In particular we show analytically that the critical point in the ground state is located at finite $p_0$.

preprint2014arXiv

Non compact conformal field theory and the a_2^{(2)} (Izergin-Korepin) model in regime III

The so-called regime III of the a_2^{(2)} Izergin-Korepin 19-vertex model has defied understanding for many years. We show in this paper that its continuum limit involves in fact a non compact conformal field theory (the so-called Witten Euclidian black hole CFT), which leads to a continuous spectrum of critical exponents, as well as very strong corrections to scaling. Detailed numerical evidence based on the Bethe ansatz analysis is presented, involving in particular the observation of discrete states in the spectrum, in full agreement with the string theory prediction for the black hole CFT. Our results have important consequences for the physics of the O(n) model, which will be discussed elsewhere.

preprint2014arXiv

Non compact continuum limit of two coupled Potts models

We study two $Q$-state Potts models coupled by the product of their energy operators, in the regime $2 < Q \le 4$ where the coupling is relevant. A particular choice of weights on the square lattice is shown to be equivalent to the integrable $a_3^{(2)}$ vertex model. It corresponds to a selfdual system of two antiferromagnetic Potts models, coupled ferromagnetically. We derive the Bethe Ansatz equations and study them numerically for two arbitrary twist angles. The continuum limit is shown to involve two compact bosons and one non compact boson, with discrete states emerging from the continuum at appropriate twists. The non compact boson entails strong logarithmic corrections to the finite-size behaviour of the scaling levels, the understanding of which allows us to correct an earlier proposal for some of the critical exponents. In particular, we infer the full set of magnetic scaling dimensions (watermelon operators) of the Potts model.

preprint2014arXiv

Operator content of the critical Potts model in d dimensions and logarithmic correlations

Using the symmetric group $S_Q$ symmetry of the $Q$-state Potts model, we classify the (scalar) operator content of its underlying field theory in arbitrary dimension. In addition to the usual identity, energy and magnetization operators, we find fields that generalize the $N$-cluster operators well-known in two dimensions, together with their subleading counterparts. We give the explicit form of all these operators -- up to non-universal constants -- both on the lattice and in the continuum limit for the Landau theory. We compute exactly their two- and three-point correlation functions on an arbitrary graph in terms of simple probabilities, and give the general form of these correlation functions in the continuum limit at the critical point. Specializing to integer values of the parameter $Q$, we argue that the analytic continuation of the $S_Q$ symmetry yields logarithmic correlations at the critical point in arbitrary dimension, thus implying a mixing of some scaling fields by the scale transformation generator. All these logarithmic correlation functions are given a clear geometrical meaning, which can be checked in numerical simulations. Several physical examples are discussed, including bond percolation, spanning trees and forests, resistor networks and the Ising model. We also briefly address the generalization of our approach to the $O(n)$ model.

preprint2014arXiv

Phase diagram and the strong-coupling fixed point in the disordered O(n) loop model

We numerically study the phase diagram and critical properties of the two-dimensional disordered O(n) loop model by using the transfer matrix and the worm Monte Carlo methods. The renormalization group flow is extracted from the landscape of the effective central charge obtained by the transfer matrix method based on the Zamolodchikov's C-theorem. We find a line of random fixed points (FPs) for $n_c < n <1$, with $n_c \sim 0.5$, for which the central charge and critical exponents agree well with the results of the $1-n$ perturbative expansion. Furthermore, for $n> n_c$, we find a line of multicritical FPs at strong randomness. The FP at $n=1$ has $c=0.4612(4)$, which suggests that it belongs to the universality class of the Nishimori point in the $\pm J$ random-bond Ising model. For $n>2$, we find another critical line that connects the hard-hexagon FP in the pure model to a finite-randomness zero-temperature FP.

preprint2014arXiv

Universal entanglement crossover of coupled quantum wires

We consider the entanglement between two one-dimensional quantum wires (Luttinger Liquids) coupled by tunneling through a quantum impurity. The physics of the system involves a crossover between weak and strong coupling regimes characterized by an energy scale $T_B$, and methods of conformal field theory therefore cannot be applied. The evolution of the entanglement in this crossover has led to many numerical studies, but has remained little understood, analytically or even qualitatively. We argue in this Letter that the correct universal scaling form of the entanglement entropy $S$ (for an arbitrary interval of length $L$ containing the impurity) is $\partial S/\partial \ln L = f(L T_B)$. In the special case where the coupling to the impurity can be refermionized, we show how the universal function $f(L T_B)$ can be obtained analytically using recent results on form factors of twist fields and a defect massless-scattering formalism. Our results are carefully checked against numerical simulations.

preprint2013arXiv

A generalized Beraha conjecture for non-planar graphs

We study the partition function Z_{G(nk,k)}(Q,v) of the Q-state Potts model on the family of (non-planar) generalized Petersen graphs G(nk,k). We study its zeros in the plane (Q,v) for 1<= k <= 7. We also consider two specializations of Z_{G(nk,k)}, namely the chromatic polynomial P_{G(nk,k)}(Q) (corresponding to v=-1), and the flow polynomial Phi_{G(nk,k)}(Q) (corresponding to v=-Q). In these two cases, we study their zeros in the complex Q-plane for 1 <= k <= 7. We pay special attention to the accumulation loci of the corresponding zeros when n -> infinity. We observe that the Berker-Kadanoff phase that is present in two-dimensional Potts models, also exists for non-planar recursive graphs. Their qualitative features are the same; but the main difference is that the role played by the Beraha numbers for planar graphs is now played by the non-negative integers for non-planar graphs. At these integer values of Q, there are massive eigenvalue cancellations, in the same way as the eigenvalue cancellations that happen at the Beraha numbers for planar graphs.

preprint2013arXiv

A physical approach to the classification of indecomposable Virasoro representations from the blob algebra

In the context of Conformal Field Theory (CFT), many results can be obtained from the representation theory of the Virasoro algebra. While the interest in Logarithmic CFTs has been growing recently, the Virasoro representations corresponding to these quantum field theories remain dauntingly complicated, thus hindering our understanding of various critical phenomena. We extend in this paper the construction of Read and Saleur (2007), and uncover a deep relationship between the Virasoro algebra and a finite-dimensional algebra characterizing the properties of two-dimensional statistical models, the so-called blob algebra (a proper extension of the Temperley--Lieb algebra). This allows us to explore vast classes of Virasoro representations (projective, tilting, generalized staggered modules, etc.), and to conjecture a classification of all possible indecomposable Virasoro modules (with, in particular, L_0 Jordan cells of arbitrary rank) that may appear in a consistent physical Logarithmic CFT where Virasoro is the maximal local chiral algebra. As by-products, we solve and analyze algebraically quantum-group symmetric XXZ spin chains and sl(2|1) supersymmetric spin chains with extra spins at the boundary, together with the "mirror" spin chain introduced by Martin and Woodcock (2004).

preprint2013arXiv

Two-dimensional Potts antiferromagnets with a phase transition at arbitrarily large q

We exhibit infinite families of two-dimensional lattices (some of which are triangulations or quadrangulations of the plane) on which the q-state Potts antiferromagnet has a finite-temperature phase transition at arbitrarily large values of q. This unexpected result is proven rigorously by using a Peierls argument to measure the entropic advantage of sublattice long-range order. Additional numerical data are obtained using transfer matrices, Monte Carlo simulation, and a high-precision graph-theoretic method.

preprint2012arXiv

An integrable spin chain for the SL(2,R)/U(1) black hole sigma model

We introduce a spin chain based on finite-dimensional spin-1/2 SU(2) representations but with a non-hermitian `Hamiltonian' and show, using mostly analytical techniques, that it is described at low energies by the SL(2,R)/U(1) Euclidian black hole Conformal Field Theory. This identification goes beyond the appearance of a non-compact spectrum: we are also able to determine the density of states, and show that it agrees with the formulas in [J. Math. Phys. 42, 2961 (2001)] and [JHEP 04, 014 (2002)], hence providing a direct `physical measurement' of the associated reflection amplitude.

preprint2012arXiv

Conformal boundary state for the rectangular geometry

We discuss conformal field theories (CFTs) in rectangular geometries, and develop a formalism that involves a conformal boundary state for the 1+1d open system. We focus on the case of homogeneous boundary conditions (no insertion of a boundary condition changing operator), for which we derive an explicit expression of the associated boundary state, valid for any arbitrary CFT. We check the validity of our solution, comparing it with known results for partition functions, numerical simulations of lattice discretizations, and coherent state expressions for free theories.

preprint2012arXiv

Critical manifold of the kagome-lattice Potts model

Any two-dimensional infinite regular lattice G can be produced by tiling the plane with a finite subgraph B of G; we call B a basis of G. We introduce a two-parameter graph polynomial P_B(q,v) that depends on B and its embedding in G. The algebraic curve P_B(q,v) = 0 is shown to provide an approximation to the critical manifold of the q-state Potts model, with coupling v = exp(K)-1, defined on G. This curve predicts the phase diagram both in the ferromagnetic (v>0) and antiferromagnetic (v<0) regions. For larger bases B the approximations become increasingly accurate, and we conjecture that P_B(q,v) = 0 provides the exact critical manifold in the limit of infinite B. Furthermore, for some lattices G, or for the Ising model (q=2) on any G, P_B(q,v) factorises for any choice of B: the zero set of the recurrent factor then provides the exact critical manifold. In this sense, the computation of P_B(q,v) can be used to detect exact solvability of the Potts model on G. We illustrate the method for the square lattice, where the Potts model has been exactly solved, and the kagome lattice, where it has not. For the square lattice we correctly reproduce the known phase diagram, including the antiferromagnetic transition and the singularities in the Berker-Kadanoff phase. For the kagome lattice, taking the smallest basis with six edges we recover a well-known (but now refuted) conjecture of F.Y. Wu. Larger bases provide successive improvements on this formula, giving a natural extension of Wu's approach. The polynomial predictions are in excellent agreement with numerical computations. For v>0 the accuracy of the predicted critical coupling v_c is of the order 10^{-4} or 10^{-5} for the 6-edge basis, and improves to 10^{-6} or 10^{-7} for the largest basis studied (with 36 edges).

preprint2012arXiv

Exact exponents for the spin quantum Hall transition in the presence of multiple edge channels

Critical properties of quantum Hall systems are affected by the presence of extra edge channels - present, in particular, at higher plateau transitions. We study this phenomenon for the case of the spin quantum Hall transition. Using supersymmetry we map the corresponding network model to a classical loop model, whose boundary critical behavior was recently determined exactly. We verify predictions of the exact solution by extensive numerical simulations.

preprint2012arXiv

Logarithmic observables in critical percolation

Although it has long been known that the proper quantum field theory description of critical percolation involves a logarithmic conformal field theory (LCFT), no direct consequence of this has been observed so far. Representing critical bond percolation as the Q = 1 limit of the Q-state Potts model, and analyzing the underlying S_Q symmetry of the Potts spins, we identify a class of simple observables whose two-point functions scale logarithmically for Q = 1. The logarithm originates from the mixing of the energy operator with a logarithmic partner that we identify as the field that creates two propagating clusters. In d=2 dimensions this agrees with general LCFT results, and in particular the universal prefactor of the logarithm can be computed exactly. We confirm its numerical value by extensive Monte-Carlo simulations.

preprint2012arXiv

Potts critical frontiers of inhomogeneous and asymmetric bow-tie lattices

We study the critical frontiers of the Potts model on two-dimensional bow-tie lattices with fully inhomogeneous coupling constants. Generally, for the Potts critical frontier to be found exactly, the underlying lattice must be a 3-uniform hypergraph. A more general class of lattices are the 4-uniform ones, with unit cells contained within four boundary vertices. We demonstrate that in some cases, such lattices can be decomposed into triangular cells, and solved using a modification of standard techniques. This leads to the exact inhomogeneous Potts critical frontiers on various lattices, such as the bow-tie lattice with five different couplings, and critical points for asymmetric bow-tie lattices.

preprint2012arXiv

Rectangular amplitudes, conformal blocks, and applications to loop models

In this paper we continue the investigation of partition functions of critical systems on a rectangle initiated in [R. Bondesan et al, Nucl.Phys.B862:553-575,2012]. Here we develop a general formalism of rectangle boundary states using conformal field theory, adapted to describe geometries supporting different boundary conditions. We discuss the computation of rectangular amplitudes and their modular properties, presenting explicit results for the case of free theories. In a second part of the paper we focus on applications to loop models, discussing in details lattice discretizations using both numerical and analytical calculations. These results allow to interpret geometrically conformal blocks, and as an application we derive new probability formulas for self-avoiding walks.

preprint2012arXiv

The Hintermann-Merlini-Baxter-Wu and the Infinite-Coupling-Limit Ashkin-Teller Models

We show how the Hintermann-Merlini-Baxter-Wu model (which is a generalization of the well-known Baxter-Wu model to a general Eulerian triangulation) can be mapped onto a particular infinite-coupling-limit of the Ashkin-Teller model. We work out some mappings among these models, also including the standard and mixed Ashkin-Teller models. Finally, we compute the phase diagram of the infinite-coupling-limit Ashkin-Teller model on the square, triangular, hexagonal, and kagome lattices.

preprint2012arXiv

The puzzle of bulk conformal field theories at central charge c=0

Non-trivial critical models in 2D with central charge c=0 are described by Logarithmic Conformal Field Theories (LCFTs), and exhibit in particular mixing of the stress-energy tensor with a "logarithmic" partner under a conformal transformation. This mixing is quantified by a parameter (usually denoted b), introduced in [V. Gurarie, Nucl. Phys. B 546, 765 (1999)], and which was first thought to play the role of an "effective" central charge. The value of b has been determined over the last few years for the boundary versions of these models: $b_{\rm perco}=-5/8$ for percolation and $b_{\rm poly} = 5/6$ for dilute polymers. Meanwhile, the existence and value of $b$ for the bulk theory has remained an open problem. Using lattice regularization techniques we provide here an "experimental study" of this question. We show that, while the chiral stress tensor has indeed a single logarithmic partner in the chiral sector of the theory, the value of b is not the expected one: instead, b=-5 for both theories. We suggest a theoretical explanation of this result using operator product expansions and Coulomb gas arguments, and discuss the physical consequences on correlation functions. Our results imply that the relation between bulk LCFTs of physical interest and their boundary counterparts is considerably more involved than in the non-logarithmic case.

preprint2012arXiv

Transfer matrix computation of critical polynomials for two-dimensional Potts models

In our previous work we have shown that critical manifolds of the q-state Potts model can be studied by means of a graph polynomial P_B(q,v), henceforth referred to as the critical polynomial. This polynomial may be defined on any periodic two-dimensional lattice. It depends on a finite subgraph B, called the basis, and the manner in which B is tiled to construct the lattice. The real roots v = e^K - 1 of P_B(q,v) either give the exact critical points for the lattice, or provide approximations that, in principle, can be made arbitrarily accurate by increasing the size of B in an appropriate way. In earlier work, P_B(q,v) was defined by a contraction-deletion identity, similar to that satisfied by the Tutte polynomial. Here, we give a probabilistic definition of P_B(q,v), which facilitates its computation, using the transfer matrix, on much larger B than was previously possible. We present results for the critical polynomial on the (4,8^2), kagome, and (3,12^2) lattices for bases of up to respectively 96, 162, and 243 edges, compared to the limit of 36 edges with contraction-deletion. We discuss in detail the role of the symmetries and the embedding of B. The critical temperatures v_c obtained for ferromagnetic (v>0) Potts models are at least as precise as the best available results from Monte Carlo simulations or series expansions. For instance, with q=3 we obtain v_c(4,8^2) = 3.742 489(4), v_c(kagome) = 1.876 459 7(2), and v_c(3,12^2) = 5.033 078 49(4), the precision being comparable or superior to the best simulation results. More generally, we trace the critical manifolds in the real (q,v) plane and discuss the intricate structure of the phase diagram in the antiferromagnetic (v<0) region.

preprint2012arXiv

Transfer matrix computation of generalised critical polynomials in percolation

Percolation thresholds have recently been studied by means of a graph polynomial $P_B(p)$, henceforth referred to as the critical polynomial, that may be defined on any periodic lattice. The polynomial depends on a finite subgraph $B$, called the basis, and the way in which the basis is tiled to form the lattice. The unique root of $P_B(p)$ in $[0,1]$ either gives the exact percolation threshold for the lattice, or provides an approximation that becomes more accurate with appropriately increasing size of $B$. Initially $P_B(p)$ was defined by a contraction-deletion identity, similar to that satisfied by the Tutte polynomial. Here, we give an alternative probabilistic definition of $P_B(p)$, which allows for much more efficient computations, by using the transfer matrix, than was previously possible with contraction-deletion. We present bond percolation polynomials for the $(4,8^2)$, kagome, and $(3,12^2)$ lattices for bases of up to respectively 96, 162, and 243 edges, much larger than the previous limit of 36 edges using contraction-deletion. We discuss in detail the role of the symmetries and the embedding of $B$. For the largest bases, we obtain the thresholds $p_c(4,8^2) = 0.676 803 329 ...$, $p_c(\mathrm{kagome}) = 0.524 404 998 ...$, $p_c(3,12^2) = 0.740 420 798 ...$, comparable to the best simulation results. We also show that the alternative definition of $P_B(p)$ can be applied to study site percolation problems.

preprint2011arXiv

Corner free energies and boundary effects for Ising, Potts and fully-packed loop models on the square and triangular lattices

We obtain long series expansions for the bulk, surface and corner free energies for several two-dimensional statistical models, by combining Enting's finite lattice method (FLM) with exact transfer matrix enumerations. The models encompass all integrable curves of the Q-state Potts model on the square and triangular lattices, including the antiferromagnetic transition curves and the Ising model (Q=2) at temperature T, as well as a fully-packed O(n) type loop model on the square lattice. The expansions are around the trivial fixed points at infinite Q, n or 1/T. By using a carefully chosen expansion parameter, q << 1, all expansions turn out to be of the form \prod_{k=1}^\infty (1-q^k)^{α_k + k β_k}, where the coefficients α_k and β_k are periodic functions of k. Thanks to this periodicity property we can conjecture the form of the expansions to all orders (except in a few cases where the periodicity is too large). These expressions are then valid for all 0 <= q < 1. We analyse in detail the q \to 1^- limit in which the models become critical. In this limit the divergence of the corner free energy defines a universal term which can be compared with the conformal field theory (CFT) predictions of Cardy and Peschel. This allows us to deduce the asymptotic expressions for the correlation length in several cases. Finally we work out the FLM formulae for the case where some of the system's boundaries are endowed with particular (non-free) boundary conditions. We apply this in particular to the square-lattice Potts model with Jacobsen-Saleur boundary conditions, conjecturing the expansions of the surface and corner free energies to arbitrary order for any integer value of the boundary interaction parameter r. These results are in turn compared with CFT predictions.

preprint2011arXiv

Critical properties of joint spin and Fortuin-Kasteleyn observables in the two-dimensional Potts model

The two-dimensional Potts model can be studied either in terms of the original Q-component spins, or in the geometrical reformulation via Fortuin-Kasteleyn (FK) clusters. While the FK representation makes sense for arbitrary real values of Q by construction, it was only shown very recently that the spin representation can be promoted to the same level of generality. In this paper we show how to define the Potts model in terms of observables that simultaneously keep track of the spin and FK degrees of freedom. This is first done algebraically in terms of a transfer matrix that couples three different representations of a partition algebra. Using this, one can study correlation functions involving any given number of propagating spin clusters with prescribed colours, each of which contains any given number of distinct FK clusters. For 0 <= Q <= 4 the corresponding critical exponents are all of the Kac form h_{r,s}, with integer indices r,s that we determine exactly both in the bulk and in the boundary versions of the problem. In particular, we find that the set of points where an FK cluster touches the hull of its surrounding spin cluster has fractal dimension d_{2,1} = 2 - 2 h_{2,1}. If one constrains this set to points where the neighbouring spin cluster extends to infinity, we show that the dimension becomes d_{1,3} = 2 - 2 h_{1,3}. Our results are supported by extensive transfer matrix and Monte Carlo computations.

preprint2011arXiv

Edge states and conformal boundary conditions in super spin chains and super sigma models

The sigma models on projective superspaces CP^{N+M-1|N} with topological angle theta=pi mod 2pi flow to non-unitary, logarithmic conformal field theories in the low-energy limit. In this paper, we determine the exact spectrum of these theories for all open boundary conditions preserving the full global symmetry of the model, generalizing recent work on the particular case M=0 [C. Candu et al, JHEP02(2010)015]. In the sigma model setting, these boundary conditions are associated with complex line bundles, and are labelled by an integer, related with the exact value of theta. Our approach relies on a spin chain regularization, where the boundary conditions now correspond to the introduction of additional edge states. The exact values of the exponents then follow from a lengthy algebraic analysis, a reformulation of the spin chain in terms of crossing and non-crossing loops (represented as a certain subalgebra of the Brauer algebra), and earlier results on the so-called one- and two-boundary Temperley Lieb algebras (also known as blob algebras). A remarkable result is that the exponents, in general, turn out to be irrational. The case M=1 has direct applications to the spin quantum Hall effect, which will be discussed in a sequel.

preprint2011arXiv

Finite-temperature phase transition in a class of 4-state Potts antiferromagnets

We argue that the 4-state Potts antiferromagnet has a finite-temperature phase transition on any Eulerian plane triangulation in which one sublattice consists of vertices of degree 4. We furthermore predict the universality class of this transition. We then present transfer-matrix and Monte Carlo data confirming these predictions for the cases of the union-jack and bisected hexagonal lattices.

preprint2011arXiv

Indecomposability parameters in chiral Logarithmic Conformal Field Theory

Work of the last few years has shown that the key algebraic features of Logarithmic Conformal Field Theories (LCFTs) are already present in some finite lattice systems (such as the XXZ spin-1/2 chain) before the continuum limit is taken. This has provided a very convenient way to analyze the structure of indecomposable Virasoro modules and to obtain fusion rules for a variety of models such as (boundary) percolation etc. LCFTs allow for additional quantum numbers describing the fine structure of the indecomposable modules, and generalizing the `b-number' introduced initially by Gurarie for the c=0 case. The determination of these indecomposability parameters has given rise to a lot of algebraic work, but their physical meaning has remained somewhat elusive. In a recent paper, a way to measure b for boundary percolation and polymers was proposed. We generalize this work here by devising a general strategy to compute matrix elements of Virasoro generators from the numerical analysis of lattice models and their continuum limit. The method is applied to XXZ spin-1/2 and spin-1 chains with open (free) boundary conditions. They are related to gl(n+m|m) and osp(n+2m|2m)-invariant superspin chains and to nonlinear sigma models with supercoset target spaces. These models can also be formulated in terms of dense and dilute loop gas. We check the method in many cases where the results were already known analytically. Furthermore, we also confront our findings with a construction generalizing Gurarie's, where logarithms emerge naturally in operator product expansions to compensate for apparently divergent terms. This argument actually allows us to compute indecomposability parameters in any logarithmic theory. A central result of our study is the construction of a Kac table for the indecomposability parameters of the logarithmic minimal models LM(1,p) and LM(p,p+1).

preprint2010arXiv

A tree-decomposed transfer matrix for computing exact Potts model partition functions for arbitrary graphs, with applications to planar graph colourings

Combining tree decomposition and transfer matrix techniques provides a very general algorithm for computing exact partition functions of statistical models defined on arbitrary graphs. The algorithm is particularly efficient in the case of planar graphs. We illustrate it by computing the Potts model partition functions and chromatic polynomials (the number of proper vertex colourings using Q colours) for large samples of random planar graphs with up to N=100 vertices. In the latter case, our algorithm yields a sub-exponential average running time of ~ exp(1.516 sqrt(N)), a substantial improvement over the exponential running time ~ exp(0.245 N) provided by the hitherto best known algorithm. We study the statistics of chromatic roots of random planar graphs in some detail, comparing the findings with results for finite pieces of a regular lattice.

preprint2010arXiv

Bulk, surface and corner free energy series for the chromatic polynomial on the square and triangular lattices

We present an efficient algorithm for computing the partition function of the q-colouring problem (chromatic polynomial) on regular two-dimensional lattice strips. Our construction involves writing the transfer matrix as a product of sparse matrices, each of dimension ~ 3^m, where m is the number of lattice spacings across the strip. As a specific application, we obtain the large-q series of the bulk, surface and corner free energies of the chromatic polynomial. This extends the existing series for the square lattice by 32 terms, to order q^{-79}. On the triangular lattice, we verify Baxter's analytical expression for the bulk free energy (to order q^{-40}), and we are able to conjecture exact product formulae for the surface and corner free energies.

preprint2010arXiv

Conformal field theory at central charge c=0: a measure of the indecomposability (b) parameters

A good understanding of conformal field theory (CFT) at c=0 is vital to the physics of disordered systems, as well as geometrical problems such as polymers and percolation. Steady progress has shown that these CFTs should be logarithmic, with indecomposable operator product expansions, and indecomposable representations of the Virasoro algebra. In one of the earliest papers on the subject, V. Gurarie introduced a single parameter b to quantify this indecomposability in terms of the logarithmic partner t of the stress energy tensor T. He and A. Ludwig conjectured further that b=-5/8 for polymers and b=5/6 for percolation. While a lot of physics may be hidden behind this parameter - which has also given rise to a lot of discussions - it had remained very elusive up to now, due to the lack of available methods to measure it experimentally or numerically, in contrast say with the central charge. We show in this paper how to overcome the many difficulties in trying to measure b. This requires control of a lattice scalar product, lattice Jordan cells, together with a precise construction of the state L_{-2}|0>. The final result is that b=5/6 for polymers. For percolation, we find that b=-5/8 within an XXZ or supersymmetric representation. In the geometrical representation, we do not find a Jordan cell for L_0 at level two (finite-size Hamiltonian and transfer matrices are fully diagonalizable), so there is no b in this case.

preprint2010arXiv

Demixing of compact polymers chains in three dimensions

We present a Monte Carlo algorithm that provides efficient and unbiased sampling of polymer melts consisting of two chains of equal length that jointly visit all the sites of a cubic lattice with rod geometry L x L x rL and non-periodic (hard wall) boundary conditions. Using this algorithm for chains of length up to 40 000 monomers and aspect ratios 1 <= r <= 10, we show that in the limit of a large lattice the two chains phase separate. This demixing phenomenon is present already for r=1, and becomes more pronounced, albeit not perfect, as r is increased.

preprint2009arXiv

A Temperley-Lieb quantum chain with two- and three-site interactions

We study a quantum chain of spin-1/2 particles whose world lines form a loop gas. In addition to the usual interaction that projects a pair of neighboring spins onto quantum-group spin 0, we introduce a three-site interaction that projects onto spin 3/2. Changing the relative strengths of these interactions leads to a rich phase diagram containing gapped and critical phases, and exactly solvable points. For loop weight $n=2 \cos [π/(p+1)]$, the model corresponds to N=p-3 species of interacting anyons, of which only the Fibonacci case N=1 was studied previously. Our study applies to real n, and provides further exact results which were unnoticed even for the special case N=1.

preprint2009arXiv

The $Z_2$ staggered vertex model and its applications

New solvable vertex models can be easily obtained by staggering the spectral parameter in already known ones. This simple construction reveals some surprises: for appropriate values of the staggering, highly non-trivial continuum limits can be obtained. The simplest case of staggering with period two (the $Z_2$ case) for the six-vertex model was shown to be related, in one regime of the spectral parameter, to the critical antiferromagnetic Potts model on the square lattice, and has a non-compact continuum limit. Here, we study the other regime: in the very anisotropic limit, it can be viewed as a zig-zag spin chain with spin anisotropy, or as an anyonic chain with a generic (non-integer) number of species. From the Bethe-Ansatz solution, we obtain the central charge $c=2$, the conformal spectrum, and the continuum partition function, corresponding to one free boson and two Majorana fermions. Finally, we obtain a massive integrable deformation of the model on the lattice. Interestingly, its scattering theory is a massive version of the one for the flow between minimal models. The corresponding field theory is argued to be a complex version of the $C_2^{(2)}$ Toda theory.

preprint2009arXiv

The semiflexible fully-packed loop model and interacting rhombus tilings

Motivated by a recent adsorption experiment [M.O. Blunt et al., Science 322, 1077 (2008)], we study tilings of the plane with three different types of rhombi. An interaction disfavors pairs of adjacent rhombi of the same type. This is shown to be a special case of a model of fully-packed loops with interactions between monomers at distance two along a loop. We solve the latter model using Coulomb gas techniques and show that its critical exponents vary continuously with the interaction strenght. At low temperature it undergoes a Kosterlitz-Thouless transition to an ordered phase, which is predicted from numerics to occur at a temperature T \sim 110K in the experiments.

preprint2008arXiv

Critical points in coupled Potts models and critical phases in coupled loop models

We show how to couple two critical Q-state Potts models to yield a new self-dual critical point. We also present strong evidence of a dense critical phase near this critical point when the Potts models are defined in their completely packed loop representations. In the continuum limit, the new critical point is described by an SU(2) coset conformal field theory, while in this limit of the the critical phase, the two loop models decouple. Using a combination of exact results and numerics, we also obtain the phase diagram in the presence of vacancies. We generalize these results to coupling two Potts models at different Q.

preprint2007arXiv

Exact valence bond entanglement entropy and probability distribution in the XXX spin chain and the Potts model

By relating the ground state of Temperley-Lieb hamiltonians to partition functions of 2D statistical mechanics systems on a half plane, and using a boundary Coulomb gas formalism, we obtain in closed form the valence bond entanglement entropy as well as the valence bond probability distribution in these ground states. We find in particular that for the XXX spin chain, the number N_c of valence bonds connecting a subsystem of size L to the outside goes, in the thermodynamic limit, as <N_c> = (4/pi^2) ln L, disproving a recent conjecture that this should be related with the von Neumann entropy, and thus equal to 1/(3 ln 2) ln L. Our results generalize to the Q-state Potts model.

preprint2007arXiv

Unbiased sampling of globular lattice proteins in three dimensions

We present a Monte Carlo method that allows efficient and unbiased sampling of Hamiltonian walks on a cubic lattice. Such walks are self-avoiding and visit each lattice site exactly once. They are often used as simple models of globular proteins, upon adding suitable local interactions. Our algorithm can easily be equipped with such interactions, but we study here mainly the flexible homopolymer case where each conformation is generated with uniform probability. We argue that the algorithm is ergodic and has dynamical exponent z=0. We then use it to study polymers of size up to 64^3 = 262144 monomers. Results are presented for the effective interaction between end points, and the interaction with the boundaries of the system.

preprint2006arXiv

Transfer Matrices and Partition-Function Zeros for Antiferromagnetic Potts Models. IV. Chromatic polynomial with cyclic boundary conditions

We study the chromatic polynomial P_G(q) for m \times n square- and triangular-lattice strips of widths 2\leq m \leq 8 with cyclic boundary conditions. This polynomial gives the zero-temperature limit of the partition function for the antiferromagnetic q-state Potts model defined on the lattice G. We show how to construct the transfer matrix in the Fortuin--Kasteleyn representation for such lattices and obtain the accumulation sets of chromatic zeros in the complex q-plane in the limit n\to\infty. We find that the different phases that appear in this model can be characterized by a topological parameter. We also compute the bulk and surface free energies and the central charge.

preprint2005arXiv

Spanning forests and the q-state Potts model in the limit q \to 0

We study the q-state Potts model with nearest-neighbor coupling v=e^{βJ}-1 in the limit q,v \to 0 with the ratio w = v/q held fixed. Combinatorially, this limit gives rise to the generating polynomial of spanning forests; physically, it provides information about the Potts-model phase diagram in the neighborhood of (q,v) = (0,0). We have studied this model on the square and triangular lattices, using a transfer-matrix approach at both real and complex values of w. For both lattices, we have computed the symbolic transfer matrices for cylindrical strips of widths 2 \le L \le 10, as well as the limiting curves of partition-function zeros in the complex w-plane. For real w, we find two distinct phases separated by a transition point w=w_0, where w_0 = -1/4 (resp. w_0 = -0.1753 \pm 0.0002) for the square (resp. triangular) lattice. For w > w_0 we find a non-critical disordered phase, while for w < w_0 our results are compatible with a massless Berker-Kadanoff phase with conformal charge c = -2 and leading thermal scaling dimension x_{T,1} = 2 (marginal operator). At w = w_0 we find a "first-order critical point": the first derivative of the free energy is discontinuous at w_0, while the correlation length diverges as w \downarrow w_0 (and is infinite at w = w_0). The critical behavior at w = w_0 seems to be the same for both lattices and it differs from that of the Berker-Kadanoff phase: our results suggest that the conformal charge is c = -1, the leading thermal scaling dimension is x_{T,1} = 0, and the critical exponents are ν= 1/d = 1/2 and α= 1.

preprint2005arXiv

The arboreal gas and the supersphere sigma model

We discuss the relationship between the phase diagram of the Q=0 state Potts model, the arboreal gas model, and the supersphere sigma model S^{0,2} = OSP(1/2) / OSP(0/2). We identify the Potts antiferromagnetic critical point with the critical point of the arboreal gas (at negative tree fugacity), and with a critical point of the sigma model. We show that the corresponding conformal theory on the square lattice has a non-linearly realized OSP(2/2) = SL(1/2) symmetry, and involves non-compact degrees of freedom, with a continuous spectrum of critical exponents. The role of global topological properties in the sigma model transition is discussed in terms of a generalized arboreal gas model.

preprint2004arXiv

Exact Potts Model Partition Functions for Strips of the Triangular Lattice

We present exact calculations of the Potts model partition function Z(G,q,v) for arbitrary q and temperature-like variable v on n-vertex strip graphs G of the triangular lattice for a variety of transverse widths equal to L vertices and for arbitrarily great length equal to m vertices, with free longitudinal boundary conditions and free and periodic transverse boundary conditions. These have the form Z(G,q,v)=\sum_{j=1}^{N_{Z,G,λ}} c_{Z,G,j}(λ_{Z,G,j})^{m-1}. We give general formulas for N_{Z,G,j} and its specialization to v=-1 for arbitrary L. The free energy is calculated exactly for the infinite-length limit of the graphs, and the thermodynamics is discussed. It is shown how the internal energy calculated for the case of cylindrical boundary conditions is connected with critical quantities for the Potts model on the infinite triangular lattice. Considering the full generalization to arbitrary complex q and v, we determine the singular locus {\cal B}, arising as the accumulation set of partition function zeros as m\to\infty, in the q plane for fixed v and in the v plane for fixed q.

preprint2003arXiv

Transfer Matrices and Partition-Function Zeros for Antiferromagnetic Potts Models III. Triangular-Lattice Chromatic Polynomial

We study the chromatic polynomial P_G(q) for m \times n triangular-lattice strips of widths m <= 12_P, 9_F (with periodic or free transverse boundary conditions, respectively) and arbitrary lengths n (with free longitudinal boundary conditions). The chromatic polynomial gives the zero-temperature limit of the partition function for the q-state Potts antiferromagnet. We compute the transfer matrix for such strips in the Fortuin--Kasteleyn representation and obtain the corresponding accumulation sets of chromatic zeros in the complex q-plane in the limit n\to\infty. We recompute the limiting curve obtained by Baxter in the thermodynamic limit m,n\to\infty and find new interesting features with possible physical consequences. Finally, we analyze the isolated limiting points and their relation with the Beraha numbers.

preprint2002arXiv

Conformal field theory of the Flory model of polymer melting

We study the scaling limit of a fully packed loop model in two dimensions, where the loops are endowed with a bending rigidity. The scaling limit is described by a three-parameter family of conformal field theories, which we characterize via its Coulomb-gas representation. One choice for two of the three parameters reproduces the critical line of the exactly solvable six-vertex model, while another corresponds to the Flory model of polymer melting. Exact central charge and critical exponents are calculated for polymer melting in two dimensions. Contrary to predictions from mean-field theory we show that polymer melting, as described by the Flory model, is continuous. We test our field theoretical results against numerical transfer matrix calculations.

preprint2002arXiv

Parafermionic theory with the symmetry Z_5

A parafermionic conformal theory with the symmetry Z_5 is constructed, based on the second solution of Fateev-Zamolodchikov for the corresponding parafermionic chiral algebra. The primary operators of the theory, which are the singlet, doublet 1, doublet 2, and disorder operators, are found to be accommodated by the weight lattice of the classical Lie algebra B_2. The finite Kac tables for unitary theories are defined and the formula for the conformal dimensions of primary operators is given.

preprint2001arXiv

Phase diagram and critical exponents of a Potts gauge glass

The two-dimensional q-state Potts model is subjected to a Z_q symmetric disorder that allows for the existence of a Nishimori line. At q=2, this model coincides with the +/- J random-bond Ising model. For q>2, apart from the usual pure and zero-temperature fixed points, the ferro/paramagnetic phase boundary is controlled by two critical fixed points: a weak disorder point, whose universality class is that of the ferromagnetic bond-disordered Potts model, and a strong disorder point which generalizes the usual Nishimori point. We numerically study the case q=3, tracing out the phase diagram and precisely determining the critical exponents. The universality class of the Nishimori point is inconsistent with percolation on Potts clusters.

preprint2001arXiv

Transfer Matrices and Partition-Function Zeros for Antiferromagnetic Potts Models II. Extended Results for Square-Lattice Chromatic Polynomial

We study the chromatic polynomials for m \times n square-lattice strips, of width 9 <= m <= 13 (with periodic boundary conditions) and arbitrary length n (with free boundary conditions). We have used a transfer matrix approach that allowed us also to extract the limiting curves when n \to \infty. In this limit we have also obtained the isolated limiting points for these square-lattice strips and checked some conjectures related to the Beraha numbers.

preprint2001arXiv

Unusual corrections to scaling in the 3-state Potts antiferromagnet on a square lattice

At zero temperature, the 3-state antiferromagnetic Potts model on a square lattice maps exactly onto a point of the 6-vertex model whose long-distance behavior is equivalent to that of a free scalar boson. We point out that at nonzero temperature there are two distinct types of excitation: vortices, which are relevant with renormalization-group eigenvalue 1/2; and non-vortex unsatisfied bonds, which are strictly marginal and serve only to renormalize the stiffness coefficient of the underlying free boson. Together these excitations lead to an unusual form for the corrections to scaling: for example, the correlation length diverges as β\equiv J/kT \to \infty according to ξ\sim A e^{2β} (1 + bβe^{-β} + ...), where b is a nonuniversal constant that may nevertheless be determined independently. A similar result holds for the staggered susceptibility. These results are shown to be consistent with the anomalous behavior found in the Monte Carlo simulations of Ferreira and Sokal.

preprint2000arXiv

Duality relations for M coupled Potts models

We establish explicit duality transformations for systems of M q-state Potts models coupled through their local energy density, generalising known results for M=1,2,3. The M-dimensional space of coupling constants contains a selfdual sub-manifold of dimension D_M = [M/2]. For the case M=4, the variation of the effective central charge along the selfdual surface is investigated by numerical transfer matrix techniques. Evidence is given for the existence of a family of critical points, corresponding to conformal field theories with an extended S_M symmetry algebra.

preprint1997arXiv

Critical Behaviour of Random Bond Potts Models: A Transfer Matrix Study

We study the two-dimensional Potts model on the square lattice in the presence of quenched random-bond impurities. For q>4 the first-order transitions of the pure model are softened due to the impurities, and we determine the resulting universality classes by combining transfer matrix data with conformal invariance. The magnetic exponent beta/nu varies continuously with q, assuming non-Ising values for q>4, whereas the correlation length exponent nu is numerically consistent with unity. We present evidence for the correctness of a formerly proposed phase diagram, unifying pure, percolative and non-trivial random behaviour.