Quasiperiodic Cluster-Ising Chain

• Quasiperiodic Cluster-Ising Chain is an exactly solvable quantum spin chain model exhibiting topological quasiperiodic fixed points at criticality, bridging clean and infinite-randomness universality classes.
• It employs the Jordan–Wigner transformation to map spins to free fermions, enabling exact diagonalization of the Bogoliubov–de Gennes Hamiltonian and precise studies of entanglement and correlators.
• Distinct boundary critical exponents and robust topological edge modes, protected by Z2 symmetries, delineate a unique phase in modulated quantum spin systems.

The quasiperiodic cluster-Ising chain is an exactly solvable quantum spin chain model exhibiting a novel class of topological quasiperiodic (QP) fixed points at criticality. These fixed points interpolate between clean and infinite-randomness critical behavior and are characterized by indistinguishable bulk universal properties but distinct, robust topological edge features. The model provides a comprehensive framework for studying topological classification in modulated quantum critical systems, particularly in aperiodic (quasiperiodic) environments where conventional clean or random universality paradigms no longer suffice (Yang et al., 1 Feb 2026).

1. Model Definition and Hamiltonian

The chain consists of $N$ sites, each with spin-$\frac12$ degrees of freedom represented by Pauli operators $\sigma_i^{x,y,z}$. The system is invariant under two global $\mathbb{Z}_2$ symmetries: the spin-flip operator $P = \prod_i \sigma_i^z$, and time-reversal symmetry $T$, implemented as complex conjugation.

The Hamiltonian is given by: $$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$ (for open chains, the last terms are omitted).

Coupling strengths $J_i$ (nearest-neighbor) and $g_i$ (next-nearest, cluster interaction) are modulated quasiperiodically. Define a Diophantine irrational $Q$ (e.g., $\frac12$0), two phases $\frac12$1, $\frac12$2, and amplitudes $\frac12$3 with means $\frac12$4: $\frac12$5 Modulation is called "strong" if $\frac12$6, and "weak" otherwise (irrelevant to RG at weak amplitude).

2. Exact Solution Methodology

The system admits an exact solution via the Jordan–Wigner transformation, mapping spins to free spinless fermions: $\frac12$7 The spin Hamiltonian becomes quadratic in $\frac12$8 with pairing up to next-nearest neighbor: $\frac12$9 where the Bogoliubov–de Gennes (BdG) matrix $\sigma_i^{x,y,z}$0 is $\sigma_i^{x,y,z}$1 real symmetric-antisymmetric.

Diagonalization yields single-particle energies $\sigma_i^{x,y,z}$2 and the many-body ground state is the Bogoliubov vacuum. All correlators (spin correlations, nonlocal string orders, entanglement) are computed via Wick’s theorem from the two-point functions $\sigma_i^{x,y,z}$3, $\sigma_i^{x,y,z}$4.

The Majorana fermion representation, $\sigma_i^{x,y,z}$5, recasts the Hamiltonian as: $\sigma_i^{x,y,z}$6 Edge-localized zero-modes satisfy a simple recurrence; their localization-delocalization transition (phase boundary) is set by $\sigma_i^{x,y,z}$7.

3. Bulk Critical Properties

Bulk criticality shows universal features that interpolate between clean and random systems. For large size $\sigma_i^{x,y,z}$8 (with $\sigma_i^{x,y,z}$9 a rational approximant to $\mathbb{Z}_2$0):

• Entanglement entropy:

$\mathbb{Z}_2$1

• Clean-like regime: $\mathbb{Z}_2$2
• Strongly modulated QP-Ising: $\mathbb{Z}_2$3
  • Energy gap (finite-size scaling):

$\mathbb{Z}_2$4

• Weak QP: $\mathbb{Z}_2$5
• Strong QP: $\mathbb{Z}_2$6
  • Bulk spin-spin correlator:

$\mathbb{Z}_2$7

• Wandering of reduced coupling:

$\mathbb{Z}_2$8

This logarithmic wandering, with $\mathbb{Z}_2$9 nonzero, places QP criticality intermediate between clean ($P = \prod_i \sigma_i^z$0) and strong random ($P = \prod_i \sigma_i^z$1).

4. Topological Distinction and Edge Structure

Topological features are manifest in both nonlocal string order parameters and boundary critical exponents:

• Nonlocal disorder/string operators:
  • Ising side ($P = \prod_i \sigma_i^z$2-neutral):

  $P = \prod_i \sigma_i^z$3 - Cluster/SPT side ($P = \prod_i \sigma_i^z$4-charged):

  $P = \prod_i \sigma_i^z$5

• Bulk topological invariant: The $P = \prod_i \sigma_i^z$6 charge of $P = \prod_i \sigma_i^z$7 (the disorder operator) at criticality.

• Boundary operator scaling (OBC):

$P = \prod_i \sigma_i^z$8

• Trivial QP-Ising: $P = \prod_i \sigma_i^z$9, no entanglement degeneracy.

• Topological QP-Ising: $T$0, robust twofold degeneracy in all low-lying entanglement levels.

  • Robustness: Small symmetry-preserving perturbations $T$1 ($T$2) do not affect edge degeneracy or relative decay rates of string order parameters at criticality.

5. Phase Diagram and Boundary Characterization

At fixed mean couplings $T$3, the $T$4 parameter space encompasses four phases:

1. FM: clean ferromagnetic state
2. SPT: clean cluster symmetry-protected topological phase
3. QP-FM: quasiperiodically modulated ferromagnet
4. QP-SPT: gapless but area-law entangled quasiperiodic SPT

Three boundary lines converge at $T$5:

• Vertical ($T$6): QP "Ising"-type critical line.
• Curved phase boundaries, exactly given by $T$7 (arising from the average over cosine modulations). For $T$8, the transition is analytically:

$T$9

These describe transitions between FM and SPT phases as the nature and strength of quasiperiodic modulation is tuned.

6. Comparison with Established Universality Classes

A summary of universality classes relevant to the QP cluster-Ising chain is presented below:

Universality Class	Bulk Exponents $$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$0	Boundary Exponent $$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$1	Entanglement Structure
Clean Ising CFT	$$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$2	$$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$3 (trivial), $$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$4 ($$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$5-enriched)	None
Infinite-randomness Ising (IRFP)	$$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$6	$$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$7	None
QP-Ising (Crowley et al.)	$$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$8	$$H = -\sum_{i=1}^{N-1} J_i\,\sigma_i^x\,\sigma_{i+1}^x -\sum_{i=1}^{N-2} g_i\,\sigma_i^x\,\sigma_{i+1}^z\,\sigma_{i+2}^x - \bigl(J_N \sigma_N^x\sigma_1^x + g_{N-1}\sigma_{N-1}^x\sigma_N^z\sigma_1^x + g_N\sigma_N^x\sigma_1^z\sigma_2^x\bigr)$$9	None
Topological QP-Ising (cluster-Ising)	$J_i$0	$J_i$1	Robust twofold degeneracy

The topological QP-Ising fixed point discovered in the cluster-Ising chain has identical bulk exponents to previously studied QP systems, but features distinct boundary scaling ($J_i$2), a robust entanglement spectrum degeneracy, and pronounced SPT string order at criticality. These features confirm that boundary phenomena differentiate QP-Ising universality classes even when bulk criticality appears indistinguishable. The topological distinction is protected by $J_i$3 symmetry and cannot be removed without a phase transition or breaking said symmetry (Yang et al., 1 Feb 2026).