pith. sign in

arxiv: 2510.21269 · v2 · submitted 2025-10-24 · ❄️ cond-mat.stat-mech

Tensor Renormalization-Group study of the surface critical behavior of a frustrated two-layer Ising model

Christophe Chatelain (LPCT) This is my paper

Pith reviewed 2026-05-18 05:08 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords frustrated Ising modelsurface critical behaviortensor renormalization groupAshkin-Teller universalityscaling dimensionsduality relationtwo-layer Ising model
0
0 comments X

The pith

In the frustrated two-layer Ising model the surface magnetic scaling dimensions split and satisfy the duality x1^s equals one quarter x2^s.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies the surface critical behavior of two replicas of the 2D Ising model coupled by frustrated interactions. It applies tensor renormalization group calculations to show that the two-fold degeneracy of the surface magnetic scaling dimension present in the Ashkin-Teller class is lifted. The splitting occurs because frustration breaks the symmetry that reverses spins in only one replica. The resulting distinct exponents obey the simple relation that one is exactly one quarter the other.

Core claim

In the frustrated two-layer Ising model the two-fold degeneracy of the surface magnetic scaling dimension of the Ashkin-Teller universality class is lifted due to the breaking of the Z2 symmetry under spin reversal of a single Ising replica. The two distinct surface magnetic scaling dimensions x1^s and x2^s satisfy the duality relation x1^s equals 1/4 x2^s.

What carries the argument

The Bond-Weight Tensor Renormalization Group algorithm extended to boundaries, which numerically extracts the surface magnetic scaling dimensions.

If this is right

  • Surface criticality in the F2LIM differs from the standard Ashkin-Teller case because of the broken replica symmetry.
  • The duality supplies a precise numerical prediction for the ratio of the two surface exponents.
  • The same symmetry-breaking mechanism can be used to predict surface behavior in related marginally coupled models.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The duality may reflect a hidden relation in the renormalization flow of surface operators that is not visible in the bulk.
  • Analogous splitting of surface exponents could appear in other multi-replica models with selective symmetry breaking.
  • The relation offers a concrete target for analytic calculations of boundary conformal field theories in the presence of frustration.

Load-bearing premise

The extended bond-weight tensor renormalization group algorithm accurately captures the surface magnetic scaling dimensions without significant finite-size or truncation artifacts that would mask or create the reported splitting.

What would settle it

An independent Monte Carlo simulation on large lattices that measures the two surface magnetic exponents and tests whether their ratio is exactly four.

Figures

Figures reproduced from arXiv: 2510.21269 by Christophe Chatelain (LPCT).

Figure 1
Figure 1. Figure 1: The different steps of the Tensor-Renormalization Group algorithm for a finite strip. The black dots correspond to the tensors. To each line is attached an index which is one of the indices of the tensors at the two edges of the line. A black dot at the crossing of four lines is therefore a rank-4 tensor. A summation over the index of each line is implicit. The pink dots correspond to the initial tensors. … view at source ↗
Figure 2
Figure 2. Figure 2: Estimates of the bulk (left) and surface (right) critical dimensions of the Ashkin-Teller model versus the inverse of the number of iterations of the BTRG algorithm at the point K = 0.04894351 (close to the Ising point) of the critical line. The dashed lines are the exact values. 0.1 0.2 0.3 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1/niter x n (PBC) 0.1 0.2 0.3 0.5 1 1.5 2 2.5 1/niter x n (OBC) [PITH_FULL_IMAGE:figure… view at source ↗
Figure 3
Figure 3. Figure 3: Estimates of the bulk (left) and surface (right) critical dimensions of the Ashkin-Teller model versus the inverse of the number of iterations of the BTRG algorithm at the point K = 0.2724481 (close to the 4-state Potts point) of the critical line. The dashed lines are the exact values [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Estimates of the bulk (left) and surface (right) critical dimensions of the Ashkin-Teller model versus the coupling J1. The dashed lines are the exact values. the exact value. Measuring the scaling dimensions after 6 or 7 iterations appears to be a reasonable choice at least for the four smallest scaling dimensions. This choice leads to stable surface critical dimensions too. However, at larger values of J… view at source ↗
Figure 5
Figure 5. Figure 5: Estimates of surface critical dimensions of the Ashkin-Teller model with Fixed Boundary Conditions (Identical on the left subplot and Mixed on the right one) versus the coupling J1. The dashed lines are the exact values. 20%, and results from a too small number of states kept at each Tensor Renormalization. The estimates of scaling dimensions with Open Boundary Conditions are plotted on the right subplot o… view at source ↗
Figure 6
Figure 6. Figure 6: Free energy difference ∆F of the two-layer Ising model with difference FBCs at J0 ≃ 0.3042 versus the variable x. The different curves corresponds to different numbers of BTRG iterations, up to 10 iterations for the steepest red curve. The main plot corresponds to χ = 32 states while in the inset the data for χ = 64 states are plotted. tension vanishes. In the ordered phase, the mixed BC favors the presenc… view at source ↗
Figure 7
Figure 7. Figure 7: Critical line in the region 0 < J1 < J2. The blue + symbols corresponds to χ = 32 states, the red × symbols to χ = 64 states and the green ◦ symbols were obtained in [8]. The critical line is presented on [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Estimates of the bulk (left) and surface (right) critical dimensions of the two-layer Ising model versus the coupling J1. The dashed lines are the conjectured values. not be employed in the following. A parabolic fit of the numerical estimate of xστ with J1 was performed (plotted as a dotted curve on [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Estimates of surface critical dimensions of the two-layer Ising model with the first set of Fixed Boundary Conditions (Identical on the left and Mixed on the right) versus the coupling J1. The dashed lines are the conjectured values. 0 0.1 0.2 0.3 0.4 0.5 0 0.5 1 1.5 2 2.5 3 J1 x n (Ident. FBC) 0 0.1 0.2 0.3 0.4 0.5 0 0.5 1 1.5 2 J1 x n (Mixed FBC) [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Estimates of surface critical dimensions of the two-layer Ising model with the second set of Fixed Boundary Conditions (Identical on the left and Mixed on the right) versus the coupling J1. The dashed lines are the conjectured values [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
read the original abstract

Two replicas of a 2D Ising model are coupled by frustrated spin-spin interactions. It is known that this inter-layer coupling is marginal and that the bulk critical behavior belongs to the Ashkin-Teller (AT) universality class, as the $J_1-J_2$ Ising model. In this work, the surface critical behavior is studied numerically by Tensor Renormalization-Group calculations. The Bond-Weight Tensor Renormalization Group algorithm is extended to tackle systems with boundaries. It is observed that the two-fold degeneracy of the surface magnetic scaling dimension of the AT model is lifted in the frustrated two-layer Ising model (F2LIM). The splitting is explained by the breaking of the ${\mathbb Z}_2$-symmetry under spin reversal of a single Ising replica in the F2LIM. The two distinct surface magnetic scaling dimensions $x_1^s$ and $x_2^s$ of the F2LIM satisfies a simple duality relation $x_1^s=1/4x_2^s$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript numerically investigates the surface critical behavior of the frustrated two-layer Ising model (F2LIM) using an extension of the Bond-Weight Tensor Renormalization Group (TRG) algorithm to systems with boundaries. Building on the known mapping of the bulk to the Ashkin-Teller (AT) universality class, the authors report that the two-fold degeneracy of the surface magnetic scaling dimension is lifted, yielding two distinct exponents x1^s and x2^s that obey the duality x1^s = (1/4) x2^s. This splitting is attributed to the explicit breaking of the single-replica Z2 symmetry under spin reversal.

Significance. If the numerical extraction is free of truncation or finite-size artifacts, the result supplies a concrete instance of how marginal inter-layer frustration and single-replica symmetry breaking modify surface operator content within an AT-like fixed point. The boundary-adapted TRG implementation itself constitutes a methodological advance that could be applied to other layered or surface-critical models.

major comments (3)
  1. [§3] §3 (Boundary-extended TRG): The description of the boundary-tensor contraction scheme and the schedule for increasing bond dimension D does not include systematic extrapolation of the extracted surface scaling dimensions to D→∞ or to infinite system size. Surface singular values are known to be more sensitive to cutoff than bulk ones; without such extrapolation the reported splitting and the precise ratio 1/4 cannot be distinguished from truncation-induced effective fields.
  2. [§4.2] §4.2 (Numerical results for surface exponents): The values of x1^s and x2^s are presented without quoted uncertainties, without a table of D-dependence, and without a direct comparison to the known degenerate AT surface exponent (x^s_AT ≈ 0.5). This omission makes it impossible to assess whether the observed duality is robust or an artifact of the chosen truncation.
  3. [§5] §5 (Symmetry-breaking argument): The claim that single-replica Z2 breaking produces exactly the ratio x1^s = (1/4) x2^s is stated as an empirical observation followed by a qualitative symmetry argument. No perturbative expansion around the AT fixed point or renormalization-group flow equation is supplied to show why the ratio must be 1/4 rather than some other number fixed by the residual symmetry.
minor comments (2)
  1. [Figure 3] Figure 3 caption: the legend does not specify the bond dimension D used for each data set; this information should be added for reproducibility.
  2. [Abstract] The abstract states that the duality is 'observed' but does not mention the range of bond dimensions or the estimated numerical precision; a brief quantitative statement would improve clarity.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which help improve the clarity and robustness of our numerical results on the surface critical behavior of the frustrated two-layer Ising model. We address each major comment point by point below.

read point-by-point responses

  1. Referee: §3 (Boundary-extended TRG): The description of the boundary-tensor contraction scheme and the schedule for increasing bond dimension D does not include systematic extrapolation of the extracted surface scaling dimensions to D→∞ or to infinite system size. Surface singular values are known to be more sensitive to cutoff than bulk ones; without such extrapolation the reported splitting and the precise ratio 1/4 cannot be distinguished from truncation-induced effective fields.

    Authors: We agree that systematic extrapolation to D→∞ is valuable for confirming the absence of truncation artifacts in surface exponents. In the revised manuscript we will add an appendix with extrapolations of x1^s and x2^s versus 1/D (or appropriate scaling variable), using the sequence of bond dimensions already computed. These extrapolations will be shown to converge to the same values and ratio reported in the main text, thereby strengthening the evidence that the observed splitting is not a cutoff effect. revision: yes

  2. Referee: §4.2 (Numerical results for surface exponents): The values of x1^s and x2^s are presented without quoted uncertainties, without a table of D-dependence, and without a direct comparison to the known degenerate AT surface exponent (x^s_AT ≈ 0.5). This omission makes it impossible to assess whether the observed duality is robust or an artifact of the chosen truncation.

    Authors: We will revise §4.2 to include a table listing x1^s and x2^s for successive bond dimensions D, together with uncertainties estimated from the linear fits to the scaling of the surface singular values. We will also add an explicit comparison to the AT surface exponent ≈0.5, highlighting that the two distinct values in the F2LIM bracket the degenerate AT value and thereby demonstrate the lifting of degeneracy. revision: yes

  3. Referee: §5 (Symmetry-breaking argument): The claim that single-replica Z2 breaking produces exactly the ratio x1^s = (1/4) x2^s is stated as an empirical observation followed by a qualitative symmetry argument. No perturbative expansion around the AT fixed point or renormalization-group flow equation is supplied to show why the ratio must be 1/4 rather than some other number fixed by the residual symmetry.

    Authors: The ratio 1/4 is observed consistently in our TRG data for multiple system sizes and is tied to the residual symmetry after single-replica Z2 breaking. The qualitative argument in the manuscript links this factor to the operator mapping inherited from the AT duality. While a perturbative RG analysis around the AT fixed point would provide a more rigorous derivation, such an expansion lies outside the numerical scope of the present work; we therefore retain the combination of high-precision numerics and symmetry considerations as the primary support for the reported relation. revision: partial

standing simulated objections not resolved
  • A full perturbative expansion or explicit renormalization-group flow equation that derives the exact ratio 1/4 from the broken symmetry without relying on numerical input.

Circularity Check

0 steps flagged

Numerical TRG extraction of surface exponents is self-contained with no circular derivation

full rationale

The paper extends the Bond-Weight TRG algorithm to boundaries and numerically extracts surface magnetic scaling dimensions x1^s and x2^s for the frustrated two-layer Ising model, observing that they satisfy the relation x1^s = (1/4) x2^s while comparing to the known Ashkin-Teller surface exponents. This is a direct simulation result from boundary tensor singular values and correlation functions; no parameters are fitted to the target duality, no self-definitional equations are used, and the central claim does not reduce to a self-citation chain or imported ansatz. The work remains self-contained against external benchmarks for the AT class.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The numerical result assumes that the extended TRG algorithm faithfully represents the surface operators of the model and that the inter-layer coupling remains marginal at the surface. No explicit free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The inter-layer coupling is marginal and the bulk critical behavior belongs to the Ashkin-Teller universality class.
    Stated in the abstract as background knowledge used to interpret the surface results.

pith-pipeline@v0.9.0 · 5712 in / 1226 out tokens · 27108 ms · 2026-05-18T05:08:15.077259+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages

  1. [1]

    Several Tensor Renormalization-Group calculations [1, 2, 3] questioned well-established results obtained earlier by means of Monte Carlo simulations [4, 5, 6, 7]

    Introduction The frustrated J1 − J2 Ising model has recently attracted a renewed interest. Several Tensor Renormalization-Group calculations [1, 2, 3] questioned well-established results obtained earlier by means of Monte Carlo simulations [4, 5, 6, 7]. In this context, we constructed a frustrated two-layer Ising model (F2LIM) model sharing the same scali...

    work page

  2. [2]

    Boundary BTRG algorithm The calculations were performed by means of a Tensor Renormalization-Group (TRG) algorithm [27, 28] adapted to a system with boundaries

    Numerical details 2.1. Boundary BTRG algorithm The calculations were performed by means of a Tensor Renormalization-Group (TRG) algorithm [27, 28] adapted to a system with boundaries. Such a modification was recently presented for the Higher-Order Tensor Renormalization-Group (HoTRG) algorithm [29] and the Tensor-Network Renormalization (TNR) [30]. Our al...

    work page

  3. [3]

    The Ashkin-Teller model The classical Ashkin-Teller model corresponds to two Ising replicas, denoted σA and σB in the following, coupled by their local energy

    Surface critical behavior of the Ashkin-Teller model 3.1. The Ashkin-Teller model The classical Ashkin-Teller model corresponds to two Ising replicas, denoted σA and σB in the following, coupled by their local energy. The Hamiltonian reads −βH = J2 X i,j σA i,j[σA i+1,j + σA i,j+1] + J2 X i,j σB i,j[σB i+1,j + σB i,j+1] + J1 X i,j σA i,jσB i,j[σA i+1,jσB ...

    work page

  4. [4]

    The intra-layer coupling J2 is ferromagnetic while the inter-layer one J1 is ferromagnetic in one direction of the lattice and anti-ferromagnetic in the other one

    Surface critical behavior of the J1 − J2 model We consider a frustrated two-layer Ising model with Hamiltonian −βH = J2 X i,j σA i,j[σA i+1,j + σA i,j+1] + J2 X i,j σB i,j[σB i+1,j + σB i,j+1] + J1 X i,j σA i,j[σB i+1,j − σB i,j+1] + J1 X i,j σB i,j[σA i+1,j − σA i,j+1] (13) introduced in [8]. The intra-layer coupling J2 is ferromagnetic while the inter-l...

    work page

  5. [5]

    17 is similar to the duality relation between electric scaling dimensionsxe(n) = n2/2g and magnetic ones xm(m) = m2g/2 in the Coulomb gas with n = m = 1

    Eq. 17 is similar to the duality relation between electric scaling dimensionsxe(n) = n2/2g and magnetic ones xm(m) = m2g/2 in the Coulomb gas with n = m = 1. However, the g coupling that can be computed from the numerical estimate of xs 1 is significantly smaller than the expected value 2 xε = xs σ for the AT model and therefore the F2LIM [38]. The third ...

    work page

  6. [6]

    Li and L.-P

    H. Li and L.-P. Yang, Tensor network simulation for the frustrated J1 − J2 Ising model on the square lattice, Phys. Rev. E 104, 024118 (2021)

    work page 2021

  7. [7]

    Yoshiyama and K

    K. Yoshiyama and K. Hukushima, Higher-Order Tensor Renormalization Group Study of the J1 − J2 Ising Model on a Square Lattice , Phys. Rev. E 108, 054124 (2023)

    work page 2023

  8. [8]

    A. A. Gangat, Weak First-Order Phase Transitions in the Frustrated Square Lattice J1 − J2 Classical Ising Model , Phys. Rev. B 109, 104419 (2024)

    work page 2024

  9. [9]

    A. Kalz, A. Honecker, and M. Moliner, Analysis of the Phase Transition for the Ising Model on the Frustrated Square Lattice, Phys. Rev. B 84, 174407 (2011)

    work page 2011

  10. [10]

    S. Jin, A. Sen, and A. W. Sandvik, Ashkin-Teller Criticality and Pseudo-First-Order Behavior in a Frustrated Ising Model on the Square Lattice , Phys. Rev. Lett. 108, 045702 (2012)

    work page 2012

  11. [11]

    Kalz and A

    A. Kalz and A. Honecker, Location of the Potts-Critical End Point in the Frustrated Ising Model on the Square Lattice , Phys. Rev. B 86, 134410 (2012)

    work page 2012

  12. [12]

    S. Jin, A. Sen, W. Guo, and A. W. Sandvik, Phase Transitions in the Frustrated Ising Model on the Square Lattice, Phys. Rev. B 87, 144406 (2013)

    work page 2013

  13. [13]

    Chatelain, First-order transition and marginal critical behavior in a two-dimensional frustrated Ising model with efficient tensor-network implementation , Phys

    C. Chatelain, First-order transition and marginal critical behavior in a two-dimensional frustrated Ising model with efficient tensor-network implementation , Phys. Rev. E 111, 024109 (2025)

    work page 2025

  14. [14]

    Adachi, T

    D. Adachi, T. Okubo, and S. Todo, Bond-Weighted Tensor Renormalization Group, Phys. Rev. B 105, L060402 (2022)

    work page 2022

  15. [15]

    L. P. Kadanoff, Connections between the Critical Behavior of the Planar Model and That of the Eight-Vertex Model, Phys. Rev. Lett. 39, 903 (1977)

    work page 1977

  16. [16]

    Nienhuis, Critical Behavior of Two-Dimensional Spin Models and Charge Asymmetry in the Coulomb Gas, J Stat Phys 34, 731 (1984)

    B. Nienhuis, Critical Behavior of Two-Dimensional Spin Models and Charge Asymmetry in the Coulomb Gas, J Stat Phys 34, 731 (1984)

    work page 1984

  17. [17]

    R. J. Baxter, Exactly solved models of statistical mechanics , Academic Press, London (1982)

    work page 1982

  18. [18]

    Binder and P

    K. Binder and P. C. Hohenberg, Surface effects on magnetic phase transitions , Phys. Rev. B 9, 2194 (1974)

    work page 1974

  19. [19]

    Cardy, Scaling and Renormalization in Statistical Physics , Cambridge University Press, Cambridge, (1996)

    J. Cardy, Scaling and Renormalization in Statistical Physics , Cambridge University Press, Cambridge, (1996)

    work page 1996

  20. [20]

    Pleimling, Critical phenomena at perfect and non-perfect surfaces , J

    M. Pleimling, Critical phenomena at perfect and non-perfect surfaces , J. Phys. A: Math. Gen. 37, R79 (2004)

    work page 2004

  21. [21]

    J. L. Cardy, Conformal invariance and surface critical behavior, Nuclear Physics B240, 514 (1984)

    work page 1984

  22. [22]

    J. L. Cardy, Operator content of two-dimensional conformally invariant theories , Nuclear Physics B 270, 186 (1986)

    work page 1986

  23. [23]

    J. L. Cardy, Effect of boundary conditions on the operator content of two-dimensional conformally invariant theories, Nuclear Physics B 275, 200 (1986)

    work page 1986

  24. [24]

    Henkel, Conformal Invariance and Critical Phenomena , Springer, New York, (1999)

    M. Henkel, Conformal Invariance and Critical Phenomena , Springer, New York, (1999)

    work page 1999

  25. [25]

    von Gehlen and V

    G. von Gehlen and V. Rittenberg, Superconformal invariance in the Ashkin-Teller quantum chain with free boundary conditions , J. Phys. A: Math. Gen. 19 L1039 (1986). TRG study of the surface critical behavior of a frustrated two-layer Ising model 18

    work page 1986

  26. [26]

    von Gehlen and V

    G. von Gehlen and V. Rittenberg, The Ashkin-Teller quantum chain and conformal invariance, J. Phys. A: Math. Gen. 20 227 (1987)

    work page 1987

  27. [27]

    Alcaraz, M.N

    F.C. Alcaraz, M.N. Barber, M.T. Batchelor, R.J. Baxter, and G.R.W. Quispel, Surface exponents of the quantum XXZ, Ashkin-Teller and Potts models , J. Phys. A: Math. Gen. 20, 6397 (1987)

    work page 1987

  28. [28]

    Ueda and M

    A. Ueda and M. Oshikawa, Finite-Size and Finite Bond Dimension Effects of Tensor Network Renormalization, Phys. Rev. B 108, 024413 (2023)

    work page 2023

  29. [29]

    Ueda, Renormalization Group Flow and Fixed-Point in Tensor Network Representations, Ph.D dissertation, univ

    A. Ueda, Renormalization Group Flow and Fixed-Point in Tensor Network Representations, Ph.D dissertation, univ. Tokyo, arXiv:2401.18068

    work page arXiv

  30. [30]

    Huang, S.-H

    C.-Y. Huang, S.-H. Chan, Y.-J. Kao et P. Chen, Tensor network based finite-size scaling for two- dimensional Ising model , Phys. Rev. B 107, 205123 (2023)

    work page 2023

  31. [31]

    Guo, and T.-C

    W. Guo, and T.-C. Wei, Tensor network methods for extracting conformal field theory data from fixed-point tensors and defect coarse graining , Phys. Rev. E 109, 034111 (2024)

    work page 2024

  32. [32]

    Levin and C

    M. Levin and C. P. Nave, Tensor Renormalization Group Approach to 2D Classical Lattice Models, Physical Review Letters 99, 120601 (2007)

    work page 2007

  33. [33]

    Xiang, From Density Matrix and Tensor Network Renormalization, Cambridge University Press, Cambridge, (2023)

    T. Xiang, From Density Matrix and Tensor Network Renormalization, Cambridge University Press, Cambridge, (2023)

    work page 2023

  34. [34]

    S. Iino, S. Morita, and N. Kawashima, Boundary Tensor Renormalization Group , Phys. Rev. B 100, 035449 (2019)

    work page 2019

  35. [35]

    S. Iino, S. Morita, and N. Kawashima, Boundary conformal spectrum and surface critical behavior of classical spin systems: A tensor network renormalization study , Phys. Rev. B 101, 155418 (2020)

    work page 2020

  36. [36]

    R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK Users’ Guide, Society for Industrial and Applied Mathematics (1998)

    work page 1998

  37. [37]

    J. L. Cardy, Conformal invariance and universality in finite-size scaling , J. Phys. A: Math. and Gen. 17, L385 (1984)

    work page 1984

  38. [38]

    J L Cardy, Logarithmic corrections to finite-size scaling in strips, J. Phys. A: Math. Gen. 19 L1093 (1986) ; J L Cardy 1987 J. Phys. A: Math. Gen. 20 5039 (1987)

    work page 1986

  39. [39]

    Yamanaka and M

    M. Yamanaka and M. Kohmoto, Line of continuously varying criticality in the Ashkin-Teller quantum chain, Phys. Rev. B 52, 1138 (1995)

    work page 1995

  40. [40]

    D. L. O’Brien, P. A. Pearce, and S. Ole Warnaar, Finitized conformal spectrum of the Ising model on the cylinder and torus , Physica A 228, 63 (1996)

    work page 1996

  41. [41]

    N. Sh. Izmailian, Universal amplitude ratios for scaling corrections on Ising strips , Phys. Rev. E 84, 051109 (2011)

    work page 2011

  42. [42]

    Zamolodchikov, Irreversibility of the flux of the renormalization group in a 2D field theory , JETP Lett

    A.B. Zamolodchikov, Irreversibility of the flux of the renormalization group in a 2D field theory , JETP Lett. 43, 730 (1986)

    work page 1986

  43. [43]

    Saleur, Correlation functions of the critical Ashkin-Teller model on a torus , J Stat Phys 50, 475 (1988)

    H. Saleur, Correlation functions of the critical Ashkin-Teller model on a torus , J Stat Phys 50, 475 (1988)

    work page 1988