Showing 1–50 of 294 results for author: Jacobsen, J

Cultural Transmission, Technology, and Treatment of the Elderly

Authors: Matthew J. Baker, Joyce P. Jacobsen

Abstract: We discuss the interrelationship between the treatment of the elderly, the nature of production, and the transmission of culture. Respect for the elderly is endogenous. Parents cultivate an interest in consuming culture in their children; when they are older, children compensate their elders proportional to the degree to which their interests were previously cultivated. We show that this model is… ▽ More We discuss the interrelationship between the treatment of the elderly, the nature of production, and the transmission of culture. Respect for the elderly is endogenous. Parents cultivate an interest in consuming culture in their children; when they are older, children compensate their elders proportional to the degree to which their interests were previously cultivated. We show that this model is functionally equivalent to one in which cultural goods are transferred across generations. We focus on the relative well-being of the elderly and use the model to explain patterns in their relative well-being across societies. An important theme is that the cultivation of culture and norms for the respect and support of the elderly bear a nonlinear relationship with many economic variables, such as capital and or land intensity in production. We also discuss the interaction of property rights with production, assets such as productive resources, and relative treatment of the elderly. Insecurity of some types of property rights, such as rights over output, may benefit the elderly, while secure rights over productive resources may also benefit the elderly. We discuss how the elderly could be affected by demographic, technological and policy changes in both developing and developed economies. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Addressing Misspecification in Simulation-based Inference through Data-driven Calibration

Authors: Antoine Wehenkel, Juan L. Gamella, Ozan Sener, Jens Behrmann, Guillermo Sapiro, Marco Cuturi, Jörn-Henrik Jacobsen

Abstract: Driven by steady progress in generative modeling, simulation-based inference (SBI) has enabled inference over stochastic simulators. However, recent work has demonstrated that model misspecification can harm SBI's reliability. This work introduces robust posterior estimation (ROPE), a framework that overcomes model misspecification with a small real-world calibration set of ground truth parameter… ▽ More Driven by steady progress in generative modeling, simulation-based inference (SBI) has enabled inference over stochastic simulators. However, recent work has demonstrated that model misspecification can harm SBI's reliability. This work introduces robust posterior estimation (ROPE), a framework that overcomes model misspecification with a small real-world calibration set of ground truth parameter measurements. We formalize the misspecification gap as the solution of an optimal transport problem between learned representations of real-world and simulated observations. Assuming the prior distribution over the parameters of interest is known and well-specified, our method offers a controllable balance between calibrated uncertainty and informative inference under all possible misspecifications of the simulator. Our empirical results on four synthetic tasks and two real-world problems demonstrate that ROPE outperforms baselines and consistently returns informative and calibrated credible intervals. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Critical spin chains and loop models with $U(n)$ symmetry

Authors: Paul Roux, Jesper Lykke Jacobsen, Sylvain Ribault, Hubert Saleur

Abstract: Starting with the Ising model, statistical models with global symmetries provide fruitful approaches to interesting physical systems, for example percolation or polymers. These include the $O(n)$ model (symmetry group $O(n)$) and the Potts model (symmetry group $S_Q$). Both models make sense for $n,Q\in \mathbb{C}$ and not just $n,Q\in \mathbb{N}$, and both give rise to a conformal field theory in… ▽ More Starting with the Ising model, statistical models with global symmetries provide fruitful approaches to interesting physical systems, for example percolation or polymers. These include the $O(n)$ model (symmetry group $O(n)$) and the Potts model (symmetry group $S_Q$). Both models make sense for $n,Q\in \mathbb{C}$ and not just $n,Q\in \mathbb{N}$, and both give rise to a conformal field theory in the critical limit. Here, we study similar models based on the unitary group $U(n)$. We focus on the two-dimensional case, where the models can be described either as gases of non-intersecting orientable loops, or as alternating spin chains. This allows us to determine their spectra either by computing a twisted torus partition function, or by studying representations of the walled Brauer algebra. In the critical limit, our models give rise to a CFT with global $U(n)$ symmetry, which exists for any $n\in\mathbb{C}$. Its spectrum is similar to those of the $O(n)$ and Potts CFTs, but a bit simpler. We conjecture that the $O(n)$ CFT is a $\mathbb{Z}_2$ orbifold of the $U(n)$ CFT, where $\mathbb{Z}_2$ acts as complex conjugation. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 36 pages

Emerging Jordan blocks in the two-dimensional Potts and loop models at generic $Q$

Authors: Lawrence Liu, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: It was recently suggested -- based on general self-consistency arguments as well as results from the bootstrap (arXiv:2005.07708, arXiv:2007.11539, arXiv:2007.04190) -- that the CFT describing the $Q$-state Potts model is logarithmic for generic values of $Q$, with rank-two Jordan blocks for $L_0$ and ${\mkern 1.5mu\overline{\mkern-1.5mu L\mkern-1.5mu}\mkern 1.5mu}_0$ in many sectors of the theory… ▽ More It was recently suggested -- based on general self-consistency arguments as well as results from the bootstrap (arXiv:2005.07708, arXiv:2007.11539, arXiv:2007.04190) -- that the CFT describing the $Q$-state Potts model is logarithmic for generic values of $Q$, with rank-two Jordan blocks for $L_0$ and ${\mkern 1.5mu\overline{\mkern-1.5mu L\mkern-1.5mu}\mkern 1.5mu}_0$ in many sectors of the theory. This is despite the well-known fact that the lattice transfer matrix (or Hamiltonian) is diagonalizable in (arbitrary) finite size. While the emergence of Jordan blocks only in the limit $L\to\infty$ is perfectly possible conceptually, diagonalizability in finite size makes the measurement of logarithmic couplings (whose values are analytically predicted in arXiv:2007.11539, arXiv:2007.04190) very challenging. This problem is solved in the present paper (which can be considered a companion to arXiv:2007.11539), and the conjectured logarithmic structure of the CFT confirmed in detail by the study of the lattice model and associated "emerging Jordan blocks." △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:2403.09881

Lattice realization of complex CFTs: Two-dimensional Potts model with $Q>4$ states

Authors: Jesper Lykke Jacobsen, Kay Joerg Wiese

Abstract: The two-dimensional $Q$-state Potts model with real couplings has a first-order transition for $Q>4$. We study a loop-model realization in which $Q$ is a continuous parameter. This model allows for the collision of a critical and a tricritical fixed point at $Q=4$, which then emerge as complex conformally invariant theories at $Q>4$, or even complex $Q$, for suitable complex coupling constants. Al… ▽ More The two-dimensional $Q$-state Potts model with real couplings has a first-order transition for $Q>4$. We study a loop-model realization in which $Q$ is a continuous parameter. This model allows for the collision of a critical and a tricritical fixed point at $Q=4$, which then emerge as complex conformally invariant theories at $Q>4$, or even complex $Q$, for suitable complex coupling constants. All critical exponents can be obtained as analytic continuation of known exact results for $Q \le 4$. We verify this scenario in detail for $Q=5$ using transfer-matrix computations. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 5 pages, 6 figures

Journal ref: Phys. Rev. Lett. 133 (2024) 077101

Heisenberg spin chain with random-sign couplings

Authors: Michele Fava, Jesper Lykke Jacobsen, Adam Nahum

Abstract: We study the 1D quantum Heisenberg chain with randomly ferromagnetic or antiferromagnetic couplings (a model previously studied by approximate strong-disorder RG). We find that, at least for sufficiently large spin $S$, the ground state has ``spin glass'' order. The spin waves on top of this state have the dynamical exponent ${z=3/2}$, intermediate between the values $z=1$ of the antiferromagnet a… ▽ More We study the 1D quantum Heisenberg chain with randomly ferromagnetic or antiferromagnetic couplings (a model previously studied by approximate strong-disorder RG). We find that, at least for sufficiently large spin $S$, the ground state has ``spin glass'' order. The spin waves on top of this state have the dynamical exponent ${z=3/2}$, intermediate between the values $z=1$ of the antiferromagnet and ${z=2}$ of the ferromagnet. DMRG simulations are in good agreement with the analytical results for spins ${S=1}$ and ${S=3/2}$. The case ${S=1/2}$ shows large finite size effects: we suggest that this case is also ordered, but with a small ordered moment. △ Less

Submitted 23 August, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 7 pages + 14 pages appendices, 14 figures, published version

Universality of closed nested paths in two-dimensional percolation

Authors: Yu-Feng Song, Jesper Lykke Jacobsen, Bernard Nienhuis, Andrea Sportiello, Youjin Deng

Abstract: Recent work on percolation in $d=2$ [J. Phys. A {\bf 55} 204002] introduced an operator that gives a weight $k^{\ell}$ to configurations with $\ell$ `nested paths' (NP), i.e. disjoint cycles surrounding the origin, if there exists a cluster that percolates to the boundary of a disc of radius $L$, and weight zero otherwise. It was found that ${\rm E}(k^{\ell}) \sim L^{-X_{\rm NP}(k)}$, and a formul… ▽ More Recent work on percolation in $d=2$ [J. Phys. A {\bf 55} 204002] introduced an operator that gives a weight $k^{\ell}$ to configurations with $\ell$ `nested paths' (NP), i.e. disjoint cycles surrounding the origin, if there exists a cluster that percolates to the boundary of a disc of radius $L$, and weight zero otherwise. It was found that ${\rm E}(k^{\ell}) \sim L^{-X_{\rm NP}(k)}$, and a formula for $X_{\rm NP}(k)$ was conjectured. Here we derive an exact result for $X_{\rm NP}(k)$, valid for $k \ge -1$, replacing the previous conjecture. We find that the probability distribution ${\rm P}_\ell (L)$ scales as $ L^{-1/4} (\ln L)^\ell [(1/\ell!) Λ^\ell]$ when $\ell \geq 0$ and $L \gg 1$, with $Λ= 1/\sqrt{3} π$. Extensive simulations for various critical percolation models confirm our theoretical predictions and support the universality of the NP observables. △ Less

Submitted 2 September, 2024; v1 submitted 30 November, 2023; originally announced November 2023.

Comments: 20 pages, 15 figures, Version 2 is formatted according to the SciPost style, and has some small corrections, including proper acknowledgments

Critical loop models are exactly solvable

Authors: Rongvoram Nivesvivat, Sylvain Ribault, Jesper Lykke Jacobsen

Abstract: In two-dimensional critical loop models, including the $O(n)$ and Potts models, the spectrum is exactly known, as are a few structure constants or ratios thereof. Using numerical conformal bootstrap methods, we study $235$ of the simplest 4-point structure constants. For each structure constant, we find an analytic expression as a product of two factors: 1) a universal function of conformal dimens… ▽ More In two-dimensional critical loop models, including the $O(n)$ and Potts models, the spectrum is exactly known, as are a few structure constants or ratios thereof. Using numerical conformal bootstrap methods, we study $235$ of the simplest 4-point structure constants. For each structure constant, we find an analytic expression as a product of two factors: 1) a universal function of conformal dimensions, built from Barnes' double Gamma function, and 2) a polynomial function of loop weights, whose degree obeys a simple upper bound. We conjecture that all structure constants are of this form. For a few 4-point functions, we build corresponding observables in a lattice loop model. From numerical lattice results, we extract amplitude ratios that depend neither on the lattice size nor on the lattice coupling. These ratios agree with the corresponding ratios of 4-point structure constants. △ Less

Submitted 25 September, 2024; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: 68 pages, v3: corrected Table (2.68), added comments on rank 3 Jordan cells

Journal ref: SciPost Phys. 17, 029 (2024)

Integrability of rank-two web models

Authors: Augustin Lafay, Azat M. Gainutdinov, Jesper Lykke Jacobsen

Abstract: We continue our work on lattice models of webs, which generalise the well-known loop models to allow for various kinds of bifurcations [arXiv:2101.00282, arXiv:2107.10106]. Here we define new web models corresponding to each of the rank-two spiders considered by Kuperberg [arXiv:q-alg/9712003]. These models are based on the $A_2$, $G_2$ and $B_2$ Lie algebras, and their local vertex configurations… ▽ More We continue our work on lattice models of webs, which generalise the well-known loop models to allow for various kinds of bifurcations [arXiv:2101.00282, arXiv:2107.10106]. Here we define new web models corresponding to each of the rank-two spiders considered by Kuperberg [arXiv:q-alg/9712003]. These models are based on the $A_2$, $G_2$ and $B_2$ Lie algebras, and their local vertex configurations are intertwiners of the corresponding $q$-deformed quantum algebras. In all three cases we define a corresponding model on the hexagonal lattice, and in the case of $B_2$ also on the square lattice. For specific root-of-unity choices of $q$, we show the equivalence to a number of three- and four-state spin models on the dual lattice. The main result of this paper is to exhibit integrable manifolds in the parameter spaces of each web model. For $q$ on the unit circle, these models are critical and we characterise the corresponding conformal field theories via numerical diagonalisation of the transfer matrix. In the $A_2$ case we find two integrable regimes. The first one contains a dense and a dilute phase, for which we have analytic control via a Coulomb gas construction, while the second one is more elusive and likely conceals non-compact physics. Three particular points correspond to a three-state spin model with plaquette interactions, of which the one in the second regime appears to present a new universality class. In the $G_2$ case we identify four regimes numerically. The $B_2$ case is too unwieldy to be studied numerically in the general case, but it found analytically to contain a simpler sub-model based on generators of the dilute Birman-Murakami-Wenzl algebra. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 69 pages

The two upper critical dimensions of the Ising and Potts models

Authors: Kay Joerg Wiese, Jesper Lykke Jacobsen

Abstract: We derive the exact actions of the $Q$-state Potts model valid on any graph, first for the spin degrees of freedom, and second for the Fortuin-Kasteleyn clusters. In both cases the field is a traceless $Q$-component scalar field $Φ^α$. For the Ising model ($Q=2$), the field theory for the spins has upper critical dimension $d_{\rm c}^{\rm spin}=4$, whereas for the clusters it has… ▽ More We derive the exact actions of the $Q$-state Potts model valid on any graph, first for the spin degrees of freedom, and second for the Fortuin-Kasteleyn clusters. In both cases the field is a traceless $Q$-component scalar field $Φ^α$. For the Ising model ($Q=2$), the field theory for the spins has upper critical dimension $d_{\rm c}^{\rm spin}=4$, whereas for the clusters it has $d_{\rm c}^{\rm cluster}=6$. As a consequence, the probability for three points to be in the same cluster is not given by mean-field theory for $d$ within $4<d<6$. We estimate the associated universal structure constant as $C=\sqrt{6-d}+ {\cal O}(6-d)^{3/2}$. This shows that some observables in the Ising model have an upper critical dimension of 4, while others have an upper critical dimension of $6$. Combining perturbative results from the $ε=6-d$ expansion with a non-perturbative treatment close to dimension $d=4$ allows us to locate the shape of the critical domain of the Potts model in the whole $(Q,d)$ plane. △ Less

Submitted 21 April, 2024; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 31 pages, 10 figures

Journal ref: JHEP 05 (2024) 092

On currents in the $O(n)$ loop model

Authors: Jesper Lykke Jacobsen, Rongvoram Nivesvivat, Hubert Saleur

Abstract: Using methods from the conformal bootstrap, we study the properties of Noether currents in the critical $O(n)$ loop model. We confirm that they do not give rise to a Kac-Moody algebra (for $n\neq 2$), a result expected from the underlying lack of unitarity. By studying four-point functions in detail, we fully determine the current-current OPEs, and thus obtain several structure constants with phys… ▽ More Using methods from the conformal bootstrap, we study the properties of Noether currents in the critical $O(n)$ loop model. We confirm that they do not give rise to a Kac-Moody algebra (for $n\neq 2$), a result expected from the underlying lack of unitarity. By studying four-point functions in detail, we fully determine the current-current OPEs, and thus obtain several structure constants with physical meaning. We find in particular that the terms $:\!J\bar{J}\!:$ in the identity and adjoint channels vanish exactly, invalidating the argument made in \cite{car93-1} that adding orientation-dependent interactions to the model should lead to continuously varying exponents in self-avoiding walks. We also determine the residue of the identity channel in the $JJ$ two-point function, finding that it coincides both with the result of a transfer-matrix computation for an orientation-dependent correlation function in the lattice model, and with an earlier Coulomb gas computation of Cardy \cite{car93}. This is, to our knowledge, one of the first instances where the Coulomb gas formalism and the bootstrap can be successfully compared. △ Less

Submitted 28 March, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 47 pages, v2: clarifications on a few points and improved presentation

Journal ref: SciPost Phys. 16, 111 (2024)

Simulation-based Inference for Cardiovascular Models

Authors: Antoine Wehenkel, Jens Behrmann, Andrew C. Miller, Guillermo Sapiro, Ozan Sener, Marco Cuturi, Jörn-Henrik Jacobsen

Abstract: Over the past decades, hemodynamics simulators have steadily evolved and have become tools of choice for studying cardiovascular systems in-silico. While such tools are routinely used to simulate whole-body hemodynamics from physiological parameters, solving the corresponding inverse problem of mapping waveforms back to plausible physiological parameters remains both promising and challenging. Mot… ▽ More Over the past decades, hemodynamics simulators have steadily evolved and have become tools of choice for studying cardiovascular systems in-silico. While such tools are routinely used to simulate whole-body hemodynamics from physiological parameters, solving the corresponding inverse problem of mapping waveforms back to plausible physiological parameters remains both promising and challenging. Motivated by advances in simulation-based inference (SBI), we cast this inverse problem as statistical inference. In contrast to alternative approaches, SBI provides \textit{posterior distributions} for the parameters of interest, providing a \textit{multi-dimensional} representation of uncertainty for \textit{individual} measurements. We showcase this ability by performing an in-silico uncertainty analysis of five biomarkers of clinical interest comparing several measurement modalities. Beyond the corroboration of known facts, such as the feasibility of estimating heart rate, our study highlights the potential of estimating new biomarkers from standard-of-care measurements. SBI reveals practically relevant findings that cannot be captured by standard sensitivity analyses, such as the existence of sub-populations for which parameter estimation exhibits distinct uncertainty regimes. Finally, we study the gap between in-vivo and in-silico with the MIMIC-III waveform database and critically discuss how cardiovascular simulations can inform real-world data analysis. △ Less

Submitted 29 July, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Non-invertible symmetries and RG flows in the two-dimensional $O(n)$ loop model

Authors: Jesper Lykke Jacobsen, Hubert Saleur

Abstract: In a recent paper, Gorbenko and Zan [arXiv:2005.07708] observed that $O(n)$ symmetry alone does not protect the well-known renormalization group flow from the dilute to the dense phase of the two-dimensional $O(n)$ model under thermal perturbations. We show in this paper that the required "extra protection" is topological in nature, and is related to the existence of certain non-invertible topolog… ▽ More In a recent paper, Gorbenko and Zan [arXiv:2005.07708] observed that $O(n)$ symmetry alone does not protect the well-known renormalization group flow from the dilute to the dense phase of the two-dimensional $O(n)$ model under thermal perturbations. We show in this paper that the required "extra protection" is topological in nature, and is related to the existence of certain non-invertible topological defect lines. We define these defect lines and discuss the ensuing topological protection, both in the context of the $O(n)$ lattice model and in its recently understood continuum limit, which takes the form of a conformal field theory governed by an interchiral algebra. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 27 pages, 5 figures

From combinatorial maps to correlation functions in loop models

Authors: Linnea Grans-Samuelsson, Jesper Lykke Jacobsen, Rongvoram Nivesvivat, Sylvain Ribault, Hubert Saleur

Abstract: In two-dimensional statistical physics, correlation functions of the O(N) and Potts models may be written as sums over configurations of non-intersecting loops. We define sums associated to a large class of combinatorial maps (also known as ribbon graphs). We allow disconnected maps, but not maps that include monogons. Given a map with n vertices, we obtain a function of the moduli of the corres… ▽ More In two-dimensional statistical physics, correlation functions of the O(N) and Potts models may be written as sums over configurations of non-intersecting loops. We define sums associated to a large class of combinatorial maps (also known as ribbon graphs). We allow disconnected maps, but not maps that include monogons. Given a map with n vertices, we obtain a function of the moduli of the corresponding punctured Riemann surface. Due to the map's combinatorial (rather than topological) nature, that function is single-valued, and we call it an n-point correlation function. We conjecture that in the critical limit, such functions form a basis of solutions of certain conformal bootstrap equations. They include all correlation functions of the O(N) and Potts models, and correlation functions that do not belong to any known model. We test the conjecture by counting solutions of crossing symmetry for four-point functions on the sphere. △ Less

Submitted 24 August, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 40 pages, v2: clarified a few subtle points

Journal ref: SciPost Phys. 15, 147 (2023)

Algebraic Bethe Ansatz for the Open XXZ Spin Chain with Non-Diagonal Boundary Terms via $U_{\mathfrak{q}}\mathfrak{sl}_2$ Symmetry

Authors: Dmitry Chernyak, Azat M. Gainutdinov, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: We derive by the traditional algebraic Bethe ansatz method the Bethe equations for the general open XXZ spin chain with non-diagonal boundary terms under the Nepomechie constraint [J. Phys. A 37 (2004), 433-440, arXiv:hep-th/0304092]. The technical difficulties due to the breaking of $\mathsf{U}(1)$ symmetry and the absence of a reference state are overcome by an algebraic construction where the t… ▽ More We derive by the traditional algebraic Bethe ansatz method the Bethe equations for the general open XXZ spin chain with non-diagonal boundary terms under the Nepomechie constraint [J. Phys. A 37 (2004), 433-440, arXiv:hep-th/0304092]. The technical difficulties due to the breaking of $\mathsf{U}(1)$ symmetry and the absence of a reference state are overcome by an algebraic construction where the two-boundary Temperley-Lieb Hamiltonian is realised in a new $U_{\mathfrak{q}}\mathfrak{sl}_2$-invariant spin chain involving infinite-dimensional Verma modules on the edges [J. High Energy Phys. 2022 (2022), no. 11, 016, 64 pages, arXiv:2207.12772]. The equivalence of the two Hamiltonians is established by proving Schur-Weyl duality between $U_{\mathfrak{q}}\mathfrak{sl}_2$ and the two-boundary Temperley-Lieb algebra. In this framework, the Nepomechie condition turns out to have a simple algebraic interpretation in terms of quantum group fusion rules. △ Less

Submitted 17 July, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Journal ref: SIGMA 19 (2023), 046, 47 pages

ExoClock Project III: 450 new exoplanet ephemerides from ground and space observations

Authors: A. Kokori, A. Tsiaras, B. Edwards, A. Jones, G. Pantelidou, G. Tinetti, L. Bewersdorff, A. Iliadou, Y. Jongen, G. Lekkas, A. Nastasi, E. Poultourtzidis, C. Sidiropoulos, F. Walter, A. Wünsche, R. Abraham, V. K. Agnihotri, R. Albanesi, E. Arce-Mansego, D. Arnot, M. Audejean, C. Aumasson, M. Bachschmidt, G. Baj, P. R. Barroy , et al. (192 additional authors not shown)

Abstract: The ExoClock project has been created with the aim of increasing the efficiency of the Ariel mission. It will achieve this by continuously monitoring and updating the ephemerides of Ariel candidates over an extended period, in order to produce a consistent catalogue of reliable and precise ephemerides. This work presents a homogenous catalogue of updated ephemerides for 450 planets, generated by t… ▽ More The ExoClock project has been created with the aim of increasing the efficiency of the Ariel mission. It will achieve this by continuously monitoring and updating the ephemerides of Ariel candidates over an extended period, in order to produce a consistent catalogue of reliable and precise ephemerides. This work presents a homogenous catalogue of updated ephemerides for 450 planets, generated by the integration of $\sim$18000 data points from multiple sources. These sources include observations from ground-based telescopes (ExoClock network and ETD), mid-time values from the literature and light-curves from space telescopes (Kepler/K2 and TESS). With all the above, we manage to collect observations for half of the post-discovery years (median), with data that have a median uncertainty less than one minute. In comparison with literature, the ephemerides generated by the project are more precise and less biased. More than 40\% of the initial literature ephemerides had to be updated to reach the goals of the project, as they were either of low precision or drifting. Moreover, the integrated approach of the project enables both the monitoring of the majority of the Ariel candidates (95\%), and also the identification of missing data. The dedicated ExoClock network effectively supports this task by contributing additional observations when a gap in the data is identified. These results highlight the need for continuous monitoring to increase the observing coverage of the candidate planets. Finally, the extended observing coverage of planets allows us to detect trends (TTVs - Transit Timing Variations) for a sample of 19 planets. All products, data, and codes used in this work are open and accessible to the wider scientific community. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: Recommended for publication to ApJS (reviewer's comments implemented). Main body: 13 pages, total: 77 pages, 7 figures, 7 tables. Data available at http://doi.org/10.17605/OSF.IO/P298N

Spaces of states of the two-dimensional O(n) and Potts models

Authors: Jesper Lykke Jacobsen, Sylvain Ribault, Hubert Saleur

Abstract: We determine the spaces of states of the two-dimensional $O(n)$ and $Q$-state Potts models with generic parameters $n,Q\in \mathbb{C}$ as representations of their known symmetry algebras. While the relevant representations of the conformal algebra were recently worked out, it remained to determine the action of the global symmetry groups: the orthogonal group for the $O(n)$ model, and the symmetri… ▽ More We determine the spaces of states of the two-dimensional $O(n)$ and $Q$-state Potts models with generic parameters $n,Q\in \mathbb{C}$ as representations of their known symmetry algebras. While the relevant representations of the conformal algebra were recently worked out, it remained to determine the action of the global symmetry groups: the orthogonal group for the $O(n)$ model, and the symmetric group $S_Q$ for the $Q$-state Potts model. We do this by two independent methods. First we compute the twisted torus partition functions of the models at criticality. The twist in question is the insertion of a group element along one cycle of the torus: this breaks modular invariance, but allows the partition function to have a unique decomposition into characters of irreducible representations of the global symmetry group. Our second method reduces the problem to determining branching rules of certain diagram algebras. For the $O(n)$ model, we decompose representations of the Brauer algebra into representations of its unoriented Jones--Temperley--Lieb subalgebra. For the $Q$-state Potts model, we decompose representations of the partition algebra into representations of the appropriate subalgebra. We find explicit expressions for these decompositions as sums over certain sets of diagrams, and over standard Young tableaux. We check that both methods agree in many cases. Moreover, our spaces of states are consistent with recent bootstrap results on four-point functions of the corresponding CFTs. △ Less

Submitted 18 January, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 60 pages, v2: minor clarifications

Journal ref: SciPost Phys. 14, 092 (2023)

Robust Hybrid Learning With Expert Augmentation

Authors: Antoine Wehenkel, Jens Behrmann, Hsiang Hsu, Guillermo Sapiro, Gilles Louppe, Jörn-Henrik Jacobsen

Abstract: Hybrid modelling reduces the misspecification of expert models by combining them with machine learning (ML) components learned from data. Similarly to many ML algorithms, hybrid model performance guarantees are limited to the training distribution. Leveraging the insight that the expert model is usually valid even outside the training domain, we overcome this limitation by introducing a hybrid dat… ▽ More Hybrid modelling reduces the misspecification of expert models by combining them with machine learning (ML) components learned from data. Similarly to many ML algorithms, hybrid model performance guarantees are limited to the training distribution. Leveraging the insight that the expert model is usually valid even outside the training domain, we overcome this limitation by introducing a hybrid data augmentation strategy termed \textit{expert augmentation}. Based on a probabilistic formalization of hybrid modelling, we demonstrate that expert augmentation, which can be incorporated into existing hybrid systems, improves generalization. We empirically validate the expert augmentation on three controlled experiments modelling dynamical systems with ordinary and partial differential equations. Finally, we assess the potential real-world applicability of expert augmentation on a dataset of a real double pendulum. △ Less

Submitted 11 April, 2023; v1 submitted 8 February, 2022; originally announced February 2022.

Journal ref: Transaction on Machine Learning Research, 2023

Geometric Algebra and Algebraic Geometry of Loop and Potts Models

Authors: Janko Böhm, Jesper Lykke Jacobsen, Yunfeng Jiang, Yang Zhang

Abstract: 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 repr… ▽ More 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. △ Less

Submitted 7 February, 2022; originally announced February 2022.

Comments: 66 pages

Report number: USTC-ICTS/PCFT-22-03

Learning Invariant Representations with Missing Data

Authors: Mark Goldstein, Jörn-Henrik Jacobsen, Olina Chau, Adriel Saporta, Aahlad Puli, Rajesh Ranganath, Andrew C. Miller

Abstract: Spurious correlations allow flexible models to predict well during training but poorly on related test distributions. Recent work has shown that models that satisfy particular independencies involving correlation-inducing \textit{nuisance} variables have guarantees on their test performance. Enforcing such independencies requires nuisances to be observed during training. However, nuisances, such a… ▽ More Spurious correlations allow flexible models to predict well during training but poorly on related test distributions. Recent work has shown that models that satisfy particular independencies involving correlation-inducing \textit{nuisance} variables have guarantees on their test performance. Enforcing such independencies requires nuisances to be observed during training. However, nuisances, such as demographics or image background labels, are often missing. Enforcing independence on just the observed data does not imply independence on the entire population. Here we derive \acrshort{mmd} estimators used for invariance objectives under missing nuisances. On simulations and clinical data, optimizing through these estimates achieves test performance similar to using estimators that make use of the full data. △ Less

Submitted 8 June, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: CLeaR (Causal Learning and Reasoning) 2022

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

Authors: Linnea Grans-Samuelsson, Rongvoram Nivesvivat, Jesper Lykke Jacobsen, Sylvain Ribault, Hubert Saleur

Abstract: 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 fou… ▽ More 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. △ Less

Submitted 22 March, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 49 pages, v3: improved explanations on a few points

Journal ref: SciPost Phys. 12, 147 (2022)

$U_{\mathfrak{q}}(\mathfrak{sl}_3)$ web models: Locality, phase diagram and geometrical defects

Authors: Augustin Lafay, Azat M. Gainutdinov, Jesper Lykke Jacobsen

Abstract: We continue investigating the generalisations of geometrical statistical models introduced in [13], in the form of models of webs on the hexagonal lattice H having a U_q(sl_n) quantum group symmetry. We focus here on the n=3 case of cubic webs, based on the Kuperberg A_2 spider, and illustrate its properties by comparisons with the well-known dilute loop model (the n=2 case) throughout. A local ve… ▽ More We continue investigating the generalisations of geometrical statistical models introduced in [13], in the form of models of webs on the hexagonal lattice H having a U_q(sl_n) quantum group symmetry. We focus here on the n=3 case of cubic webs, based on the Kuperberg A_2 spider, and illustrate its properties by comparisons with the well-known dilute loop model (the n=2 case) throughout. A local vertex-model reformulation is exhibited, analogous to the correspondence between the loop model and a three-state vertex model. The n=3 representation uses seven states per link of H, displays explicitly the geometrical content of the webs and their U_q(sl_3) symmetry, and permits us to study the model on a cylinder via a local transfer matrix. A numerical study of the central charge reveals that for each q $\in \mathbb{C}$ in the critical regime, |q|=1, the web model possesses a dense and a dilute critical point, just like its loop model counterpart. In the dense $q=-e^{i π/4}$ case, the n=3 webs can be identified with spin interfaces of the critical three-state Potts model defined on the triangular lattice dual to H. We also provide another mapping to a $\mathbb{Z}_3$ spin model on H itself, using a high-temperature expansion. We then discuss the sector structure of the transfer matrix, for generic q, and its relation to defect configurations in both the strip and the cylinder geometries. These defects define the finite-size precursors of electromagnetic operators. This discussion paves the road for a Coulomb gas description of the conformal properties of defect webs, which will form the object of a subsequent paper. Finally, we identify the fractal dimension of critical webs in the $q=-e^{i π/3}$ case, which is the n=3 analogue of the polymer limit in the loop model. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Journal ref: Nucl. Phys. B 979 (2022) 115789

Nested Closed Paths in Two-Dimensional Percolation

Authors: Yu-Feng Song, Xiao-Jun Tan, Xin-Hang Zhang, Jesper Lykke Jacobsen, Bernard Nienhuis, Youjin Deng

Abstract: 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 connec… ▽ More 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. △ Less

Submitted 14 February, 2021; originally announced February 2021.

Comments: 5 pages, 5 figures plus supplemental material, 1 page 2 figures

MSC Class: 82B43; 60K35

$U_q(\mathfrak{sl}_n)$ web models and $\mathbb{Z}_n$ spin interfaces

Authors: Augustin Lafay, Azat M. Gainutdinov, Jesper Lykke Jacobsen

Abstract: This is the first in a series of papers devoted to generalisations of statistical loop models. We define a lattice model of $U_q(\mathfrak{sl}_n)$ webs on the honeycomb lattice, for $n \ge 2$. It is a statistical model of closed, cubic graphs with certain non-local Boltzmann weights that can be computed from spider relations. For $n=2$, the model has no branchings and reduces to the well-known O(… ▽ More This is the first in a series of papers devoted to generalisations of statistical loop models. We define a lattice model of $U_q(\mathfrak{sl}_n)$ webs on the honeycomb lattice, for $n \ge 2$. It is a statistical model of closed, cubic graphs with certain non-local Boltzmann weights that can be computed from spider relations. For $n=2$, the model has no branchings and reduces to the well-known O($N$) loop model introduced by Nienhuis. In the general case, we show that the web model possesses a particular point, at $q=e^{i π/(n+1)}$, where the partition function is proportional to that of a $\mathbb{Z}_n$-symmetric chiral spin model on the dual lattice. Moreover, under this equivalence, the graphs given by the configurations of the web model are in bijection with the domain walls of the spin model. For $n=2$, this equivalence reduces to the well-known relation between the Ising and O($1$) models. We define as well an open $U_q(\mathfrak{sl}_n)$ web model on a simply connected domain with a boundary, and discuss in particular the role of defects on the boundary. △ Less

Submitted 1 January, 2021; originally announced January 2021.

Comments: 31 pages

Journal ref: J. Stat. Mech. (2021) 053104

Lattice regularisation of a non-compact boundary conformal field theory

Authors: Niall F. Robertson, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: 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 restric… ▽ More 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. △ Less

Submitted 14 December, 2020; originally announced December 2020.

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

Authors: Linnea Grans-Samuelsson, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: 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 indec… ▽ More 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. △ Less

Submitted 22 February, 2021; v1 submitted 24 October, 2020; originally announced October 2020.

Comments: 61 pages, 13 figures, 12 tables v2: Minor edits. Corrected typo in eq (188)

Environment Inference for Invariant Learning

Authors: Elliot Creager, Jörn-Henrik Jacobsen, Richard Zemel

Abstract: Learning models that gracefully handle distribution shifts is central to research on domain generalization, robust optimization, and fairness. A promising formulation is domain-invariant learning, which identifies the key issue of learning which features are domain-specific versus domain-invariant. An important assumption in this area is that the training examples are partitioned into "domains" or… ▽ More Learning models that gracefully handle distribution shifts is central to research on domain generalization, robust optimization, and fairness. A promising formulation is domain-invariant learning, which identifies the key issue of learning which features are domain-specific versus domain-invariant. An important assumption in this area is that the training examples are partitioned into "domains" or "environments". Our focus is on the more common setting where such partitions are not provided. We propose EIIL, a general framework for domain-invariant learning that incorporates Environment Inference to directly infer partitions that are maximally informative for downstream Invariant Learning. We show that EIIL outperforms invariant learning methods on the CMNIST benchmark without using environment labels, and significantly outperforms ERM on worst-group performance in the Waterbirds and CivilComments datasets. Finally, we establish connections between EIIL and algorithmic fairness, which enables EIIL to improve accuracy and calibration in a fair prediction problem. △ Less

Submitted 15 July, 2021; v1 submitted 14 October, 2020; originally announced October 2020.

Critical percolation on the kagome hypergraph

Authors: Christian R. Scullard, Jesper Lykke Jacobsen, Robert M. Ziff

Abstract: 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 hexa… ▽ More 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. △ Less

Submitted 3 September, 2020; originally announced September 2020.

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

Authors: Linnea Grans-Samuelsson, Lawrence Liu, Yifei He, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: 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… ▽ More 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. △ Less

Submitted 22 February, 2021; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: 37 pages, 3 figures, 15 tables. v2: Minor edits to beginning of Section 7.1 and Acknowledgments v3: 38 pages. Addition of Appendix C, in which we correct some of our reasoning about the lattice evidence for the structure of the continuum modules. This correction does not change the results. Minor other edits

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

Authors: Xiaojun Tan, Youjin Deng, Jesper Lykke Jacobsen

Abstract: 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.… ▽ More 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. △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: Percolation, Logarithmic conformal field theory, Monte Carlo simulation

Journal ref: Front. Phys. 15(4), 41501 (2020)

Understanding and Mitigating Exploding Inverses in Invertible Neural Networks

Authors: Jens Behrmann, Paul Vicol, Kuan-Chieh Wang, Roger Grosse, Jörn-Henrik Jacobsen

Abstract: Invertible neural networks (INNs) have been used to design generative models, implement memory-saving gradient computation, and solve inverse problems. In this work, we show that commonly-used INN architectures suffer from exploding inverses and are thus prone to becoming numerically non-invertible. Across a wide range of INN use-cases, we reveal failures including the non-applicability of the cha… ▽ More Invertible neural networks (INNs) have been used to design generative models, implement memory-saving gradient computation, and solve inverse problems. In this work, we show that commonly-used INN architectures suffer from exploding inverses and are thus prone to becoming numerically non-invertible. Across a wide range of INN use-cases, we reveal failures including the non-applicability of the change-of-variables formula on in- and out-of-distribution (OOD) data, incorrect gradients for memory-saving backprop, and the inability to sample from normalizing flow models. We further derive bi-Lipschitz properties of atomic building blocks of common architectures. These insights into the stability of INNs then provide ways forward to remedy these failures. For tasks where local invertibility is sufficient, like memory-saving backprop, we propose a flexible and efficient regularizer. For problems where global invertibility is necessary, such as applying normalizing flows on OOD data, we show the importance of designing stable INN building blocks. △ Less

Submitted 24 December, 2021; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: AISTATS 2021

Revisiting money and labor for valuing environmental goods and services in developing countries

Authors: Habtamu Tilahun Kassahun, Jette Bredahl Jacobsen, Charles F. Nicholson

Abstract: Many Stated Preference studies conducted in developing countries provide a low willingness to pay (WTP) for a wide range of goods and services. However, recent studies in these countries indicate that this may partly be a result of the choice of payment vehicle, not the preference for the good. Thus, low WTP may not indicate a low welfare effect for public projects in developing countries. We argu… ▽ More Many Stated Preference studies conducted in developing countries provide a low willingness to pay (WTP) for a wide range of goods and services. However, recent studies in these countries indicate that this may partly be a result of the choice of payment vehicle, not the preference for the good. Thus, low WTP may not indicate a low welfare effect for public projects in developing countries. We argue that in a setting where 1) there is imperfect substitutability between money and other measures of wealth (e.g. labor), and 2) institutions are perceived to be corrupt, including payment vehicles that are currently available to the individual and less pron to corruption may be needed to obtain valid welfare estimates. Otherwise, we risk underestimating the welfare benefit of projects. We demonstrate this through a rural household contingent valuation (CV) survey designed to elicit the value of access to reliable irrigation water in Ethiopia. Of the total average annual WTP for access to reliable irrigation service, cash contribution comprises only 24.41 %. The implication is that socially desirable projects might be rejected based on cost-benefit analysis as a result of welfare gain underestimation due to mismatch of payment vehicles choice in valuation study. △ Less

Submitted 26 August, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: There is no difference between the previous version and the current version. This version is to link with the published version DOI

Journal ref: Ecological Economics 177 (2020) 106771

Geometrical four-point functions in the two-dimensional critical $Q$-state Potts model: The interchiral conformal bootstrap

Authors: Yifei He, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: Based on the spectrum identified in our earlier work [arXiv:1809.02191], we numerically solve the bootstrap to determine four-point correlation functions of the geometrical connectivities in the $Q$-state Potts model. Crucial in our approach is the existence of "interchiral conformal blocks", which arise from the degeneracy of fields with conformal weight $h_{r,1}$, with $r\in\mathbb{N}^*$, and ar… ▽ More Based on the spectrum identified in our earlier work [arXiv:1809.02191], we numerically solve the bootstrap to determine four-point correlation functions of the geometrical connectivities in the $Q$-state Potts model. Crucial in our approach is the existence of "interchiral conformal blocks", which arise from the degeneracy of fields with conformal weight $h_{r,1}$, with $r\in\mathbb{N}^*$, and are related to the underlying presence of the "interchiral algebra" introduced in [arXiv:1207.6334]. We also find evidence for the existence of "renormalized" recursions, replacing those that follow from the degeneracy of the field $Φ_{12}^D$ in Liouville theory, and obtain the first few such recursions in closed form. This hints at the possibility of the full analytical determination of correlation functions in this model. △ Less

Submitted 14 May, 2020; originally announced May 2020.

Comments: 53 pages, 30 figures

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

Authors: Etienne Granet, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: 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… ▽ More 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. △ Less

Submitted 12 May, 2020; originally announced May 2020.

Comments: 50 pages

Journal ref: JHEP article number 69 (2020)

Shortcut Learning in Deep Neural Networks

Authors: Robert Geirhos, Jörn-Henrik Jacobsen, Claudio Michaelis, Richard Zemel, Wieland Brendel, Matthias Bethge, Felix A. Wichmann

Abstract: Deep learning has triggered the current rise of artificial intelligence and is the workhorse of today's machine intelligence. Numerous success stories have rapidly spread all over science, industry and society, but its limitations have only recently come into focus. In this perspective we seek to distill how many of deep learning's problems can be seen as different symptoms of the same underlying… ▽ More Deep learning has triggered the current rise of artificial intelligence and is the workhorse of today's machine intelligence. Numerous success stories have rapidly spread all over science, industry and society, but its limitations have only recently come into focus. In this perspective we seek to distill how many of deep learning's problems can be seen as different symptoms of the same underlying problem: shortcut learning. Shortcuts are decision rules that perform well on standard benchmarks but fail to transfer to more challenging testing conditions, such as real-world scenarios. Related issues are known in Comparative Psychology, Education and Linguistics, suggesting that shortcut learning may be a common characteristic of learning systems, biological and artificial alike. Based on these observations, we develop a set of recommendations for model interpretation and benchmarking, highlighting recent advances in machine learning to improve robustness and transferability from the lab to real-world applications. △ Less

Submitted 21 November, 2023; v1 submitted 16 April, 2020; originally announced April 2020.

Comments: perspective article published at Nature Machine Intelligence (https://doi.org/10.1038/s42256-020-00257-z)

Topological defects in periodic RSOS models and anyonic chains

Authors: J. Belletête, A. M. Gainutdinov, J. L. Jacobsen, H. Saleur, T. S. Tavares

Abstract: We provide a lattice regularization of all topological defects in minimal models CFTs using RSOS and anyonic spin chains. For defects of type $(1,s)$, we connect our result with the "topological symmetry" initially identified in Fibonacci anyons [Phys. Rev. Lett. 98, 160409 (2007)], and the center of the affine Temperley-Lieb algebra discussed in [1811.02551]. We show that the topological nature o… ▽ More We provide a lattice regularization of all topological defects in minimal models CFTs using RSOS and anyonic spin chains. For defects of type $(1,s)$, we connect our result with the "topological symmetry" initially identified in Fibonacci anyons [Phys. Rev. Lett. 98, 160409 (2007)], and the center of the affine Temperley-Lieb algebra discussed in [1811.02551]. We show that the topological nature of the defects is exact on the lattice as well. Our defects of type $(r,1)$, in contrast, are only topological in the continuum limit. Identifications are obtained by a mix of algebraic and Bethe-ansatz techniques. Most of our discussion is framed in a Hamiltonian (or transfer matrix) formalism, and direct and crossed channel are both discussed in detail. For defects of type $(1,s)$, we also show how to implement their fusion, which turns out to reproduce the tensor product of the underlying monoidal category used to build the anyonic chain. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: 42 pages, 4 figures

Integrable boundary conditions in the antiferromagnetic Potts model

Authors: Niall F. Robertson, Michal Pawelkiewicz, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: 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… ▽ More 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. △ Less

Submitted 29 March, 2020; v1 submitted 6 March, 2020; originally announced March 2020.

Comments: 32 pages; typos corrected

Out-of-Distribution Generalization via Risk Extrapolation (REx)

Authors: David Krueger, Ethan Caballero, Joern-Henrik Jacobsen, Amy Zhang, Jonathan Binas, Dinghuai Zhang, Remi Le Priol, Aaron Courville

Abstract: Distributional shift is one of the major obstacles when transferring machine learning prediction systems from the lab to the real world. To tackle this problem, we assume that variation across training domains is representative of the variation we might encounter at test time, but also that shifts at test time may be more extreme in magnitude. In particular, we show that reducing differences in ri… ▽ More Distributional shift is one of the major obstacles when transferring machine learning prediction systems from the lab to the real world. To tackle this problem, we assume that variation across training domains is representative of the variation we might encounter at test time, but also that shifts at test time may be more extreme in magnitude. In particular, we show that reducing differences in risk across training domains can reduce a model's sensitivity to a wide range of extreme distributional shifts, including the challenging setting where the input contains both causal and anti-causal elements. We motivate this approach, Risk Extrapolation (REx), as a form of robust optimization over a perturbation set of extrapolated domains (MM-REx), and propose a penalty on the variance of training risks (V-REx) as a simpler variant. We prove that variants of REx can recover the causal mechanisms of the targets, while also providing some robustness to changes in the input distribution ("covariate shift"). By appropriately trading-off robustness to causally induced distributional shifts and covariate shift, REx is able to outperform alternative methods such as Invariant Risk Minimization in situations where these types of shift co-occur. △ Less

Submitted 25 February, 2021; v1 submitted 2 March, 2020; originally announced March 2020.

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

Authors: Yifei He, Linnea Grans-Samuelsson, Jesper Lykke Jacobsen, Hubert Saleur

Abstract: 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.0722… ▽ More 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. △ Less

Submitted 20 February, 2020; originally announced February 2020.

Comments: 47 pages, 12 figures

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

Authors: Zoltan Bajnok, Jesper Lykke Jacobsen, Yunfeng Jiang, Rafael I. Nepomechie, Yang Zhang

Abstract: 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… ▽ More 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. △ Less

Submitted 22 July, 2020; v1 submitted 20 February, 2020; originally announced February 2020.

Comments: More explanations added, a new appendix added, published version

Report number: CERN-TH-2020-019, UMTG-303, USTC-ICTS/PCFT-20-04

Learning the Stein Discrepancy for Training and Evaluating Energy-Based Models without Sampling

Authors: Will Grathwohl, Kuan-Chieh Wang, Jorn-Henrik Jacobsen, David Duvenaud, Richard Zemel

Abstract: We present a new method for evaluating and training unnormalized density models. Our approach only requires access to the gradient of the unnormalized model's log-density. We estimate the Stein discrepancy between the data density $p(x)$ and the model density $q(x)$ defined by a vector function of the data. We parameterize this function with a neural network and fit its parameters to maximize the… ▽ More We present a new method for evaluating and training unnormalized density models. Our approach only requires access to the gradient of the unnormalized model's log-density. We estimate the Stein discrepancy between the data density $p(x)$ and the model density $q(x)$ defined by a vector function of the data. We parameterize this function with a neural network and fit its parameters to maximize the discrepancy. This yields a novel goodness-of-fit test which outperforms existing methods on high dimensional data. Furthermore, optimizing $q(x)$ to minimize this discrepancy produces a novel method for training unnormalized models which scales more gracefully than existing methods. The ability to both learn and compare models is a unique feature of the proposed method. △ Less

Submitted 14 August, 2020; v1 submitted 13 February, 2020; originally announced February 2020.

Comments: ICML 2020

Fundamental Tradeoffs between Invariance and Sensitivity to Adversarial Perturbations

Authors: Florian Tramèr, Jens Behrmann, Nicholas Carlini, Nicolas Papernot, Jörn-Henrik Jacobsen

Abstract: Adversarial examples are malicious inputs crafted to induce misclassification. Commonly studied sensitivity-based adversarial examples introduce semantically-small changes to an input that result in a different model prediction. This paper studies a complementary failure mode, invariance-based adversarial examples, that introduce minimal semantic changes that modify an input's true label yet prese… ▽ More Adversarial examples are malicious inputs crafted to induce misclassification. Commonly studied sensitivity-based adversarial examples introduce semantically-small changes to an input that result in a different model prediction. This paper studies a complementary failure mode, invariance-based adversarial examples, that introduce minimal semantic changes that modify an input's true label yet preserve the model's prediction. We demonstrate fundamental tradeoffs between these two types of adversarial examples. We show that defenses against sensitivity-based attacks actively harm a model's accuracy on invariance-based attacks, and that new approaches are needed to resist both attack types. In particular, we break state-of-the-art adversarially-trained and certifiably-robust models by generating small perturbations that the models are (provably) robust to, yet that change an input's class according to human labelers. Finally, we formally show that the existence of excessively invariant classifiers arises from the presence of overly-robust predictive features in standard datasets. △ Less

Submitted 4 August, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: ICML 2020 (Supersedes the workshop paper "Exploiting Excessive Invariance caused by Norm-Bounded Adversarial Robustness", arXiv:1903.10484)

How to train your neural ODE: the world of Jacobian and kinetic regularization

Authors: Chris Finlay, Jörn-Henrik Jacobsen, Levon Nurbekyan, Adam M Oberman

Abstract: Training neural ODEs on large datasets has not been tractable due to the necessity of allowing the adaptive numerical ODE solver to refine its step size to very small values. In practice this leads to dynamics equivalent to many hundreds or even thousands of layers. In this paper, we overcome this apparent difficulty by introducing a theoretically-grounded combination of both optimal transport and… ▽ More Training neural ODEs on large datasets has not been tractable due to the necessity of allowing the adaptive numerical ODE solver to refine its step size to very small values. In practice this leads to dynamics equivalent to many hundreds or even thousands of layers. In this paper, we overcome this apparent difficulty by introducing a theoretically-grounded combination of both optimal transport and stability regularizations which encourage neural ODEs to prefer simpler dynamics out of all the dynamics that solve a problem well. Simpler dynamics lead to faster convergence and to fewer discretizations of the solver, considerably decreasing wall-clock time without loss in performance. Our approach allows us to train neural ODE-based generative models to the same performance as the unregularized dynamics, with significant reductions in training time. This brings neural ODEs closer to practical relevance in large-scale applications. △ Less

Submitted 23 June, 2020; v1 submitted 7 February, 2020; originally announced February 2020.

Comments: Accepted to ICML 2020

Design and performance of a 35-ton liquid argon time projection chamber as a prototype for future very large detectors

Authors: D. L. Adams, M. Baird, G. Barr, N. Barros, A. Blake, E. Blaufuss, A. Booth, D. Brailsford, N. Buchanan, B. Carls, H. Chen, M. Convery, G. De Geronimo, T. Dealtry, R. Dharmapalan, Z. Djurcic, J. Fowler, S. Glavin, R. A. Gomes, M. C. Goodman, M. Graham, L. Greenler, A. Hahn, J. Hartnell, R. Herbst , et al. (49 additional authors not shown)

Abstract: Liquid argon time projection chamber technology is an attractive choice for large neutrino detectors, as it provides a high-resolution active target and it is expected to be scalable to very large masses. Consequently, it has been chosen as the technology for the first module of the DUNE far detector. However, the fiducial mass required for "far detectors" of the next generation of neutrino oscill… ▽ More Liquid argon time projection chamber technology is an attractive choice for large neutrino detectors, as it provides a high-resolution active target and it is expected to be scalable to very large masses. Consequently, it has been chosen as the technology for the first module of the DUNE far detector. However, the fiducial mass required for "far detectors" of the next generation of neutrino oscillation experiments far exceeds what has been demonstrated so far. Scaling to this larger mass, as well as the requirement for underground construction places a number of additional constraints on the design. A prototype 35-ton cryostat was built at Fermi National Acccelerator Laboratory to test the functionality of the components foreseen to be used in a very large far detector. The Phase I run, completed in early 2014, demonstrated that liquid argon could be maintained at sufficient purity in a membrane cryostat. A time projection chamber was installed for the Phase II run, which collected data in February and March of 2016. The Phase II run was a test of the modular anode plane assemblies with wrapped wires, cold readout electronics, and integrated photon detection systems. While the details of the design do not match exactly those chosen for the DUNE far detector, the 35-ton TPC prototype is a demonstration of the functionality of the basic components. Measurements are performed using the Phase II data to extract signal and noise characteristics and to align the detector components. A measurement of the electron lifetime is presented, and a novel technique for measuring a track's position based on pulse properties is described. △ Less

Submitted 2 March, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 28 pages, 12 figures, accepted by JINST

Your Classifier is Secretly an Energy Based Model and You Should Treat it Like One

Authors: Will Grathwohl, Kuan-Chieh Wang, Jörn-Henrik Jacobsen, David Duvenaud, Mohammad Norouzi, Kevin Swersky

Abstract: We propose to reinterpret a standard discriminative classifier of p(y|x) as an energy based model for the joint distribution p(x,y). In this setting, the standard class probabilities can be easily computed as well as unnormalized values of p(x) and p(x|y). Within this framework, standard discriminative architectures may beused and the model can also be trained on unlabeled data. We demonstrate tha… ▽ More We propose to reinterpret a standard discriminative classifier of p(y|x) as an energy based model for the joint distribution p(x,y). In this setting, the standard class probabilities can be easily computed as well as unnormalized values of p(x) and p(x|y). Within this framework, standard discriminative architectures may beused and the model can also be trained on unlabeled data. We demonstrate that energy based training of the joint distribution improves calibration, robustness, andout-of-distribution detection while also enabling our models to generate samplesrivaling the quality of recent GAN approaches. We improve upon recently proposed techniques for scaling up the training of energy based models and presentan approach which adds little overhead compared to standard classification training. Our approach is the first to achieve performance rivaling the state-of-the-artin both generative and discriminative learning within one hybrid model. △ Less

Submitted 15 September, 2020; v1 submitted 6 December, 2019; originally announced December 2019.

Preventing Gradient Attenuation in Lipschitz Constrained Convolutional Networks

Authors: Qiyang Li, Saminul Haque, Cem Anil, James Lucas, Roger Grosse, Jörn-Henrik Jacobsen

Abstract: Lipschitz constraints under L2 norm on deep neural networks are useful for provable adversarial robustness bounds, stable training, and Wasserstein distance estimation. While heuristic approaches such as the gradient penalty have seen much practical success, it is challenging to achieve similar practical performance while provably enforcing a Lipschitz constraint. In principle, one can design Lips… ▽ More Lipschitz constraints under L2 norm on deep neural networks are useful for provable adversarial robustness bounds, stable training, and Wasserstein distance estimation. While heuristic approaches such as the gradient penalty have seen much practical success, it is challenging to achieve similar practical performance while provably enforcing a Lipschitz constraint. In principle, one can design Lipschitz constrained architectures using the composition property of Lipschitz functions, but Anil et al. recently identified a key obstacle to this approach: gradient norm attenuation. They showed how to circumvent this problem in the case of fully connected networks by designing each layer to be gradient norm preserving. We extend their approach to train scalable, expressive, provably Lipschitz convolutional networks. In particular, we present the Block Convolution Orthogonal Parameterization (BCOP), an expressive parameterization of orthogonal convolution operations. We show that even though the space of orthogonal convolutions is disconnected, the largest connected component of BCOP with 2n channels can represent arbitrary BCOP convolutions over n channels. Our BCOP parameterization allows us to train large convolutional networks with provable Lipschitz bounds. Empirically, we find that it is competitive with existing approaches to provable adversarial robustness and Wasserstein distance estimation. △ Less

Submitted 9 November, 2019; v1 submitted 3 November, 2019; originally announced November 2019.

Comments: 9 main pages, 31 pages total, 3 figures. Accepted at 33rd Conference on Neural Information Processing Systems (NeurIPS 2019)

Bond percolation thresholds on Archimedean lattices from critical polynomial roots

Authors: Christian R. Scullard, Jesper Lykke Jacobsen

Abstract: 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 o… ▽ More 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. △ Less

Submitted 27 October, 2019; originally announced October 2019.

Journal ref: Phys. Rev. Research 2, 012050 (2020)

On Generalized $Q$-systems

Authors: Zoltán Bajnok, Etienne Granet, Jesper Lykke Jacobsen, Rafael I. Nepomechie

Abstract: We formulate $Q$-systems for the closed XXZ, open XXX and open quantum-group-invariant XXZ quantum spin chains. Polynomial solutions of these $Q$-systems can be found efficiently, which in turn lead directly to the admissible solutions of the corresponding Bethe ansatz equations. We formulate $Q$-systems for the closed XXZ, open XXX and open quantum-group-invariant XXZ quantum spin chains. Polynomial solutions of these $Q$-systems can be found efficiently, which in turn lead directly to the admissible solutions of the corresponding Bethe ansatz equations. △ Less

Submitted 13 March, 2020; v1 submitted 17 October, 2019; originally announced October 2019.

Comments: 27 pages. New version: included code, modified sec 3.2, and a few minor changes

Journal ref: JHEP, article number 177 (2020)

On zero-remainder conditions in the Bethe ansatz

Authors: Etienne Granet, Jesper Lykke Jacobsen

Abstract: We prove that physical solutions to the Heisenberg spin chain Bethe ansatz equations are exactly obtained by imposing two zero-remainder conditions. This bridges the gap between different criteria, yielding an alternative proof of a recently devised algorithm based on $QQ$ relations, and solving its minimality issue. We prove that physical solutions to the Heisenberg spin chain Bethe ansatz equations are exactly obtained by imposing two zero-remainder conditions. This bridges the gap between different criteria, yielding an alternative proof of a recently devised algorithm based on $QQ$ relations, and solving its minimality issue. △ Less

Submitted 13 March, 2020; v1 submitted 17 October, 2019; originally announced October 2019.

Comments: 14 pages. New version: corrected typos, modified abstract and added references

Journal ref: JHEP, article number 178 (2020)

Logarithmic correlation functions for critical dense polymers on the cylinder

Authors: Alexi Morin-Duchesne, Jesper Lykke Jacobsen

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 $α\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 s… ▽ More 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 $α\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 arcs are attached to the boundary of the cylinder. For $Z(x)$, the boundary is also decorated with simple arcs, but it also has two 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 parameter $τ=\frac{x-1}n\in(0,1)$. For small $τ$, the leading behaviours are proportional to $τ^{1/4}$, $τ^{1/4}\log τ$, $\logτ$ and $\log^2τ$. 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 $Δ= -\frac18$ or $Δ=0$. We obtain differential equations satisfied by the conformal correlators, solve these equations, and find a perfect agreement with the lattice results. We 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. △ Less

Submitted 11 July, 2019; originally announced July 2019.

Comments: 38 pages