1D Cluster–Ising Model

• The one-dimensional Cluster–Ising model is defined as an exactly solvable spin chain that blends three-spin cluster interactions with antiferromagnetic couplings, enabling the study of symmetry-protected topological phases.
• It employs a Jordan–Wigner transformation to map the system to free fermions, yielding an exact solution that reveals unique entanglement properties and unconventional critical exponents.
• The model exhibits a quantum phase transition from a nonlocal string-ordered cluster phase to an antiferromagnetic phase with staggered magnetization, serving as a prototypical system for exploring topological quantum phase transitions.

The one-dimensional Cluster–Ising model is an exactly solvable spin-$\tfrac12$ chain in which a three-spin cluster interaction competes with an antiferromagnetic nearest-neighbor Ising coupling. In its standard form,

$$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$

with periodic boundary conditions unless otherwise stated, it realizes a zero-temperature transition between a cluster phase with hidden nonlocal order and an antiferromagnetic Ising phase with staggered $y$-magnetization (Smacchia et al., 2011). A commonly used equivalent parameterization writes

$$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$

where the critical point is reported at $\alpha_c=\tfrac13$ (Giampaolo et al., 2014). The model is notable for combining exact free-fermion solvability, symmetry-protected topological structure in the cluster regime, vanishing two-spin entanglement across the entire phase diagram, and critical behavior beyond the ordinary Ising universality class.

1. Hamiltonian, parameterizations, and physical content

The defining interaction of the model is the three-body cluster term

$-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$

which competes with the antiferromagnetic Ising-like term

$+\lambda\,\sigma_j^y \sigma_{j+1}^y.$

At $\lambda=0$, the Hamiltonian reduces to commuting stabilizers

$K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$

and the ground state is the cluster state $|C\rangle$, defined by $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$0 for all $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$1 (Smacchia et al., 2011). The $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$2-term favors anti-alignment of neighboring $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$3-components, so increasing $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$4 drives the chain toward a Néel-like antiferromagnetic phase.

Two conventions are used in the literature summarized here. One employs the coupling $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$5 and places the transition at $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$6 (Smacchia et al., 2011); the other uses $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$7 and places the transition at $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$8 (Giampaolo et al., 2014). This suggests that the two notations differ by an overall rescaling of the Hamiltonian and a reparameterization of the coupling.

Boundary conditions are physically consequential. Periodic boundary conditions are standard in the thermodynamic treatment of the model (Smacchia et al., 2011), whereas open boundary conditions are used in finite-size entanglement calculations, with $$H(\lambda)=-\sum_{j=1}^N \sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x + \lambda \sum_{j=1}^N \sigma_j^y \sigma_{j+1}^y,$$9 and $y$0 taken as an integer multiple of $y$1 (Giampaolo et al., 2014). In the open-chain setting, the cluster regime exhibits the symmetry-protected structure associated with a $y$2 symmetry (Smacchia et al., 2011).

2. Exact free-fermion solution and dual descriptions

Despite the three-spin interaction, the model maps to free fermions under the Jordan–Wigner transformation. In the $y$3-parameterization, the fermionized Hamiltonian is quadratic,

$y$4

and is diagonalized by Fourier and Bogoliubov transformations into

$y$5

with dispersion

$y$6

The spectral gap closes at $y$7, giving the quantum critical point (Smacchia et al., 2011).

This exact solution yields the partition function and free-energy density in closed form: $y$8 The model is also self-dual. Introducing dual spins $y$9, one finds

$$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$0

apart from boundary terms in the thermodynamic limit, so the self-dual point $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$1 coincides with criticality (Smacchia et al., 2011).

A Majorana formulation clarifies the open-chain boundary structure. With

$$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$2

the open-chain Hamiltonian leaves four uncoupled Majorana modes at $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$3,

$$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$4

which were identified as the origin of the nontrivial ground-state degeneracy and the $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$5-protected topological character of the cluster limit (Smacchia et al., 2011).

3. Phase structure, order parameters, and symmetry

The zero-temperature phase diagram contains two phases in the Hermitian model. For $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$6, the system is in the cluster phase; for $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$7, it is in the antiferromagnetic Ising phase; and $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$8 is the continuous quantum critical point (Smacchia et al., 2011). In the $$H=(-1+\alpha) \sum_{i=-l}^{l} S_{i-1}^x S_i^z S_{i+1}^x+\alpha \sum_{i=-l-1}^{l} S_i^y S_{i+1}^y,$$9-parameterization, the same structure is reported as a cluster phase for $\alpha_c=\tfrac13$0 and an antiferromagnetic phase for $\alpha_c=\tfrac13$1 (Giampaolo et al., 2014).

The antiferromagnetic phase is characterized by the staggered local order parameter

$\alpha_c=\tfrac13$2

with the exact result

$\alpha_c=\tfrac13$3

Thus the $\alpha_c=\tfrac13$4 phase is a conventional symmetry-breaking phase ordered along $\alpha_c=\tfrac13$5 (Smacchia et al., 2011).

The cluster phase has no local Landau order parameter. Its defining diagnostic is a nonlocal string order parameter,

$\alpha_c=\tfrac13$6

for which

$\alpha_c=\tfrac13$7

This nonlocal order, together with the open-chain degeneracy and $\alpha_c=\tfrac13$8 symmetry protection, identifies the cluster regime as a symmetry-protected topological phase (Smacchia et al., 2011).

The finite-size analysis in the $\alpha_c=\tfrac13$9-parameterization adds a symmetry-based distinction between ground states. In the antiferromagnetic phase, symmetry can be spontaneously broken, while in the cluster phase no ground state breaks the Hamiltonian symmetry. The cluster phase is described there as having a fourfold degenerate ground-state space, whereas the antiferromagnetic phase has a twofold degenerate ground state (Giampaolo et al., 2014).

4. Correlation pattern and entanglement structure

The model has an unusual correlation structure inherited from the three-spin interaction. Defining the basic contraction $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$0, one finds the selection rule

$-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$1

As consequences, the longitudinal magnetization and $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$2-correlators vanish identically,

$-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$3

the $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$4-correlator is nonzero only when the separation is a multiple of $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$5, and the $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$6-correlator carries the antiferromagnetic sign pattern $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$7 for odd $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$8 and $-\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$9 for even $+\lambda\,\sigma_j^y \sigma_{j+1}^y.$0 (Smacchia et al., 2011).

This correlation pattern underlies one of the model’s most distinctive quantum-information properties: pairwise entanglement vanishes for every pair of spins, at every separation, in both phases and at finite temperature. In the formulation based on concurrence,

$+\lambda\,\sigma_j^y \sigma_{j+1}^y.$1

and the same conclusion is reported for any pair of spins in the $+\lambda\,\sigma_j^y \sigma_{j+1}^y.$2-parameterization [(Smacchia et al., 2011); (Giampaolo et al., 2014)].

The absence of bipartite entanglement does not imply weak quantum correlations. In the thermally symmetric ground-state mixture, the residual multipartite entanglement satisfies

$+\lambda\,\sigma_j^y \sigma_{j+1}^y.$3

For a symmetry-broken Néel ground state, concurrence still vanishes, but the residual tangle becomes

$+\lambda\,\sigma_j^y \sigma_{j+1}^y.$4

Hence the cluster phase is maximally multipartite-entangled in this measure, while the Ising phase shows a reduction tied to the growth of staggered magnetization (Smacchia et al., 2011).

A more local refinement was obtained from the exact reduced density matrix of three adjacent spins $+\lambda\,\sigma_j^y \sigma_{j+1}^y.$5. After a local unitary transformation to $+\lambda\,\sigma_j^y \sigma_{j+1}^y.$6-form, the genuine multipartite concurrence is

$+\lambda\,\sigma_j^y \sigma_{j+1}^y.$7

with

$+\lambda\,\sigma_j^y \sigma_{j+1}^y.$8

This quantity is nonzero throughout the antiferromagnetic phase except at $+\lambda\,\sigma_j^y \sigma_{j+1}^y.$9, and it becomes nonzero already within the cluster phase, turning on at $\lambda=0$0. The same value of $\lambda=0$1 is obtained for symmetry-preserving and symmetry-breaking ground states (Giampaolo et al., 2014).

5. Critical behavior and universality class

The Hermitian Cluster–Ising transition is continuous but not in the ordinary Ising universality class. From the excitation spectrum and thermodynamic singularities, the exact critical exponents are

$\lambda=0$2

and the second derivative of the zero-temperature free energy is logarithmically singular at $\lambda=0$3 (Smacchia et al., 2011).

The non-Ising character is also visible in entanglement scaling. At criticality, the block entropy of a contiguous interval of length $\lambda=0$4 behaves as

$\lambda=0$5

which implies central charge

$\lambda=0$6

This value was traced to an exact decomposition into three Ising-like sectors, so the critical theory is effectively the sum of three Ising critical modes rather than a single Ising conformal field theory (Smacchia et al., 2011).

The three-spin genuine entanglement reproduces the same critical data. Its derivative with respect to $\lambda=0$7 has a finite-size maximum that grows as $\lambda=0$8, and in the thermodynamic limit

$\lambda=0$9

with $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$0. Finite-size scaling of this singularity yields $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$1, in agreement with the correlation-function analysis (Giampaolo et al., 2014). Thus local genuine tripartite entanglement is not merely present near the transition; it reproduces the critical scaling of the topological-to-antiferromagnetic quantum phase transition.

6. Non-Hermitian and algebraic extensions

A non-Hermitian extension adds an onsite dissipative term,

$K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$2

with $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$3 as the energy unit. The model remains exactly solvable after Jordan–Wigner fermionization because the Hamiltonian is still quadratic in fermions. The Hermitian critical point at $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$4 splits into three critical lines as non-Hermiticity increases, producing four phases: cluster, gapless, paramagnetic, and antiferromagnetic (Guo et al., 2022).

These four phases are distinguished by the string order parameter and staggered magnetization together with correlation asymptotics. The cluster phase has nonzero string order and zero staggered magnetization; the antiferromagnetic phase has zero string order and nonzero staggered magnetization; the gapless and paramagnetic phases both have vanishing bulk order parameters, but the former has algebraically decaying string correlations while the latter has exponential decay. The resulting sequence with increasing $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$5 at sufficiently large $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$6 is

$K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$7

The first two transitions are described as “KT-like” with $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$8, whereas the PM–AF transition is Ising with $K_j=\sigma_{j-1}^x \sigma_j^z \sigma_{j+1}^x,$9 (Guo et al., 2022).

A different exactly solvable generalization introduces next-nearest-neighbor and composite five-spin interactions,

$|C\rangle$0

This model is diagonalized by an algebraic generalization of the Jordan–Wigner transformation, with exact dispersion

$|C\rangle$1

Criticality is encoded in the cubic equation

$|C\rangle$2

and the chain is gapless iff at least one root lies on the unit circle. The distribution of roots inside and outside the unit circle determines which long-range or string order parameter is nonzero, and the corresponding critical lines have central charges $|C\rangle$3 or $|C\rangle$4 depending on the root configuration and mode structure (Yanagihara et al., 2020).

7. Terminological clarification and related but distinct usages

The expression “one-dimensional Cluster–Ising model” refers in the condensed-matter literature summarized above to a quantum spin chain with Pauli operators and a three-body cluster term, not to a classical Ising chain with a cluster-based observable. This distinction matters because a separate line of work studies a generalized classical one-dimensional Ising model with pairwise couplings and introduces clustered order parameters

$|C\rangle$5

to diagnose ordering in ferromagnetic, antiferromagnetic, alternating-bond, and frustrated cases (Andriushchenko et al., 2016).

That classical model has Hamiltonian

$|C\rangle$6

uses Ising variables $|C\rangle$7, and contains neither Pauli matrices nor three-spin cluster stabilizers. Its “cluster” terminology refers to connected groups of spins in locally favorable energy states, not to the cluster interaction of the quantum Cluster–Ising chain (Andriushchenko et al., 2016). A common misconception is therefore to treat cluster-based classical ordering diagnostics as direct variants of the quantum Cluster–Ising model. The two topics are conceptually adjacent only in a broad sense: both study nonstandard signatures of order in one-dimensional Ising-like systems, but they address different Hamiltonians, different observables, and different physical questions.