arXiv:2510.21269 · Journal of Statistical Mechanics Theory and Experiment

Tensor renormalization-group study of the surface critical behavior of a frustrated two-layer Ising model

Christophe Chatelain

Published Feb 1, 2026• 0 citations

Abstract

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 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 / 4 x 2 s .

DOI: 10.1088/1742-5468/ae337f

Related papers

First-Order Transition and Marginal Critical Behavior in a 2D Frustrated Ising Model with Efficient Tensor-Network Implementation

Christophe Chatelain

The phase diagram of a novel two-dimensional frustrated Ising model with both anti-ferromagnetic and ferromagnetic couplings is studied using Tensor-Network Renormalization-Group techniques. This model can be seen as two anti-ferromagnetic Ising replicas coupled by non-local spin-spin interactions, designed in such a way that the continuum limit matches that of the still debated J1 -J2 model and induces a marginal critical behavior. Our model has the advantage of having more symmetries than the J1 -J2 model and of allowing a more straightforward implementation of Tensor-Network Renormalization-Group algorithms We demonstrate the existence of two transition lines, featuring both first and second-order regimes. In the latter, the central charge and the critical exponents are shown to be compatible with the Ashkin-Teller universality class. This picture is consistent with that given by Monte Carlo simulations of the J1 -J2 model but not with recent studies with Tensor-Network techniques.

Oct 2, 20240 citations

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

Christophe Chatelain

The phase diagram of a two-dimensional frustrated Ising model with both anti-ferromagnetic and ferromagnetic couplings is studied using tensor network renormalization group techniques. This model can be seen as two antiferromagnetic Ising replicas coupled by nonlocal spin-spin interactions, designed in such a way that the continuum limit matches that of the still-debated J_{1}-J_{2} model and induces a marginal critical behavior. Our model has the advantage of having more symmetries than the J_{1}-J_{2} model and of allowing a more straightforward implementation of tensor network renormalization group algorithms. We demonstrate the existence of two transition lines, featuring both first- and second-order regimes. In the latter, the central charge and the critical exponents are shown to be compatible with the Ashkin-Teller universality class. This picture is consistent with that given by Monte Carlo simulations of the J_{1}-J_{2} model but not with recent studies with tensor network techniques.

arXiv:2410.01460Physical review. EFeb 6, 20255 citations

Higher-order tensor renormalization group study of the J1−J2 Ising model on a square lattice

Kota Yoshiyama, Koji Hukushima

Phase transitions of the J_{1}-J_{2} Ising model on a square lattice are studied using the higher-order tensor renormalization group (HOTRG) method. This system involves a competition between the ferromagnetic interaction J_{1} and antiferromagnetic interaction J_{2}, and in previous studies, weak first-order and second-order transitions were observed near the ratio g=J_{2}/|J_{1}|=1/2. It has also been suggested that the universality class of the second-order phase transition connected to the first-order transition line for g>1/2 belongs to the Ashkin-Teller class, which is characterized by a continuously varying critical exponent with g, as predicted by field-theoretical and other studies. Our results, based on the HOTRG calculations for significantly larger sizes, indicate that the region of the first-order transition is marginally narrower than that in previous studies. Furthermore, it is suggested that the region where the critical exponent changes does not necessarily coincide with the Ashkin-Teller region.

arXiv:2303.07733Physical review. ENov 15, 202322 citations

Universality and non-universality in the Ashkin-Teller model

Alessandro Giuliani

The Ashkin-Teller model can be formulated as a pair of 2D Ising models, interacting via a four-spin interaction. I consider the case of weak anisotropy (slight a-symmetry between the two Ising layers) and weak coupling. I show that the system admits two critical temperatures whose difference varies continuously with the strength of the coupling, scaling with an anomalous exponent as one let the a-symmetry parameter go to zero. The specific heat diverges logarithmically at the critical points (as for Ising) but the constant in front of the logarithm is renormalized by an anomalous critical exponent. The logarithmic divergence of the specific heat dominates only in an exponentially small interval around the critical temperatures and outside it is modified into an anomalous power law behaviour. The proof is based on an exact mapping of Ashkin-Teller into a model of (1+1)D interacting fermions and on the implementation of constructive Renormalization group methods on the fermionic system. This PhD thesis includes: (1) a review of the exact solution of 2D Ising in terms of Grassmann functional integrals; (2) a review of the Grassmann representation for a class of interacting Ising models, including Ashkin-Teller and the 8 vertex model; (3) a detailed discussion of the multiscale analysis of the Ashkin-Teller model, based on fermionic Renormalization Group methods, including the study of the flow of the effective coupling constants (in particular a proof of vanishing of the beta function for this and similar Luttinger-like models, based on modified approximate Ward identities, is included). This thesis is based on joint work with V. Mastropietro.

arXiv:cond-mat/0412176Dec 7, 20040 citations

Renormalization-group approach to the Ising model with a free surface

Theodore W. Burkhardt, Erich Eisenriegler

We use simple position-space renormalization-group procedures involving only nearest-neighbor interactions to discuss the critical behavior of the ferromagnetically coupled semi-infinite Ising model with a free surface in two and three dimensions. In three dimensions for suitably enhanced surface couplings there are surface transitions and simultaneous surface-bulk transitions as well as the bulk transition. A phase diagram, critical exponents for the various transitions, and some representative magnetization profiles are calculated. Despite their simplicity, the renormalization transformations exhibit complete sets of fixed points for describing the various phase transitions and yield phase diagrams which are a significant improvement over mean-field theory. However, the calculated magnetization profiles show some unphysical features due to the finite cluster size on which the transformations are based.

Physical review. B, Solid stateOct 1, 1977124 citations

Boundary Criticality in the 2d Random Quantum Ising Model

Gaurav Tenkila, Romain Vasseur, Andrew C. Potter

The edge of a quantum critical system can exhibit multiple distinct types of boundary criticality. We use a numerical real-space renormalization group (RSRG) to study the boundary criticality of a 2d quantum Ising model with random exchange couplings and transverse fields, whose bulk exhibits an infinite randomness critical point. This approach enables an asymptotically numerically exact extraction of universal scaling data from very large systems with many thousands of spins that cannot be efficiently simulated directly. We identify three distinct classes of boundary criticality, and extract key scaling exponents governing boundary-boundary and boundary-bulk correlations and dynamics. We anticipate that this approach can be generalized to studying a broad class of (disordered) boundary criticality, including symmetry-enriched criticality and edge modes of gapless symmetry-protected topological states, in contexts were other numerical methods are restricted to one-dimensional chains.

arXiv:2410.19038Oct 24, 20240 citations