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., $Q / 2\pi=\tau_G=(1+\sqrt{5})/2$), two phases $\phi_1$, $\phi_2$, and amplitudes $h_J, h_g \ge 0$ with means $\bar J, \bar g$: $$J_i = \bar J + h_J \cos\bigl(Q(i+\tfrac12) + \phi_1\bigr),\qquad g_i = \bar g + h_g \cos\bigl(Qi + \phi_1 + \phi_2\bigr)$$ Modulation is called "strong" if $h_J, h_g > \bar J, \bar g$, 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: $$\sigma_j^x = \prod_{k<j} (1-2c_k^\dagger c_k)(c_j^\dagger + c_j), \qquad \sigma_j^z = 1-2c_j^\dagger c_j$$ The spin Hamiltonian becomes quadratic in $(c_i, c_i^\dagger)$ with pairing up to next-nearest neighbor: $$H = \frac{1}{2}\Psi^\dagger \mathcal{H}_{\rm BdG}\Psi + \text{const.}, \qquad \Psi = (c_1,\dots, c_N, c_1^\dagger,\dots,c_N^\dagger)^T$$ where the Bogoliubov–de Gennes (BdG) matrix $\mathcal{H}_{\rm BdG}$ is $2N\times 2N$ real symmetric-antisymmetric.

Diagonalization yields single-particle energies $\epsilon_k$ 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 $G_{ij}=\langle c_i c_j^\dagger\rangle$, $F_{ij}=\langle c_i c_j\rangle$.

The Majorana fermion representation, $c_i = (\gamma_{2i-1} + i\gamma_{2i})/2$, recasts the Hamiltonian as: $$H = i \sum_i J_i \gamma_{2i} \gamma_{2i+1} + i\sum_i g_i \gamma_{2i} \gamma_{2i+3}$$ Edge-localized zero-modes satisfy a simple recurrence; their localization-delocalization transition (phase boundary) is set by $\langle\ln|g_i/J_i|\rangle=0$.

3. Bulk Critical Properties

Bulk criticality shows universal features that interpolate between clean and random systems. For large size $N=2q$ (with $q$ a rational approximant to $Q/2\pi$):

• Entanglement entropy:

$$\overline{S_{\rm vN}(q)} \approx \frac{c_{\rm eff}}{3}\ln q + \text{const.}$$

• Clean-like regime: $c_{\rm eff} = 1/2$
• Strongly modulated QP-Ising: $c_{\rm eff} = 0.63(2)$
  • Energy gap (finite-size scaling):

$\overline{\delta_e}(q) \sim q^{-z}$

• Weak QP: $z = 1.00(1)$
• Strong QP: $z = 1.8(1)$
  • Bulk spin-spin correlator:

$$\overline{C_{\rm FM}(r)} \sim \frac{1}{(r/q)^{2\Delta_\sigma^{\rm bulk}}},\qquad \Delta_\sigma^{\rm bulk}=0.176(2)$$

• Wandering of reduced coupling:

$$S_\ell(j)=\sum_{i=j}^{j+\ell-1}\ln\left|\frac{J_i}{g_i}\right|, \qquad \mathrm{Var}[S_\ell] \sim w\ln \ell, \quad w \approx 1.3$$

This logarithmic wandering, with $w$ nonzero, places QP criticality intermediate between clean ($w=0$) and strong random ($\propto\ell$).

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 ($\mathbb{Z}_2$-neutral):

  $$O_{\rm PM}(r) = \left\langle\prod_{k=i}^{i+r-1}\sigma_k^z\right\rangle$$ - Cluster/SPT side ($\mathbb{Z}_2^T$-charged):

  $$O_{\rm SPT}(r) = \left\langle \sigma_i^x \sigma_{i+1}^y \prod_{k=i+2}^{i+r-1}\sigma_k^z\,\sigma_{i+r}^y\sigma_{i+r+1}^x \right\rangle$$

• Bulk topological invariant: The $\mathbb{Z}_2^T$ charge of $\mu$ (the disorder operator) at criticality.

• Boundary operator scaling (OBC):

$$\overline{C_{\rm FM}^{\rm bdy\!-\!bulk}(r)} \sim r^{-(\Delta_\sigma^{\rm bulk} + \Delta_\sigma^{\rm bdy})}$$

• Trivial QP-Ising: $\Delta_\sigma^{\rm bdy}=0.59(2)$, no entanglement degeneracy.

• Topological QP-Ising: $\Delta_\sigma^{\rm bdy}=1.66(3)$, robust twofold degeneracy in all low-lying entanglement levels.

  • Robustness: Small symmetry-preserving perturbations $-h'\sum \sigma_i^z$ ($h'\approx 10^{-3}$) 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 $\bar J = \bar g = 1/2$, the $(h_J, h_g)$ 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 $h_J = h_g = \bar J = \bar g$:

• Vertical ($h_J=h_g$): QP "Ising"-type critical line.
• Curved phase boundaries, exactly given by $\langle\ln|g_i/J_i|\rangle = 0$ (arising from the average over cosine modulations). For $\bar J = \bar g$, the transition is analytically:

$$\frac{\bar J}{h_g} = \frac{1+(h_J/h_g)^2}{2},\quad h_J<h_g;\qquad \frac{\bar J}{h_J} = \frac{1+(h_g/h_J)^2}{2},\quad h_g<h_J$$

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 $(c_{\rm eff}, \Delta_\sigma^{\rm bulk})$	Boundary Exponent $\Delta_\sigma^{\rm bdy}$	Entanglement Structure
Clean Ising CFT	$(1/2,\;1/8)$	$1/2$ (trivial), $2$ ($T$-enriched)	None
Infinite-randomness Ising (IRFP)	$(0.347,\;0.176)$	$0$	None
QP-Ising (Crowley et al.)	$(0.63,\;0.176)$	$0.59$	None
Topological QP-Ising (cluster-Ising)	$(0.63,\;0.176)$	$1.66$	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 ($\Delta_\sigma^{\rm bdy}=1.66$), 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 $\mathbb{Z}_2^T$ symmetry and cannot be removed without a phase transition or breaking said symmetry (Yang et al., 1 Feb 2026).