Folded Ising Chain Constructions

• Folded Ising chain is a framework that transforms 1D Ising models via nonlocal changes into multi-leg, decorated, or topological geometries.
• It employs diverse constructions including mapping to cylindrical models, rarefied string dualities, chain-array renormalization, and boundary-engineered reductions.
• These methods yield exact solutions and accurate phase transition characterizations in both classical and quantum Ising systems.

“Folded Ising chain” denotes several distinct but closely related constructions in which a one-dimensional Ising system is reorganized into a multi-leg, higher-dimensional, decorated, or topological object without abandoning exact or controlled chain-based methods. In the literature, folding may mean a nonlocal change of variables that maps a multispin chain to a rectangular Ising model with helical boundary conditions, a duality that turns rarefied string operators into local spins on an $m$-leg Ising tube, an arraying of exactly solved chains to reconstruct two- and three-dimensional criticality, or a reduction of a transverse-field Ising chain to protected boundary or low-energy sectors described by Kitaev or SSH chains (Turban, 2016, Timonin et al., 2017, Kecoglu et al., 2020, Zhang et al., 2020, Ma et al., 10 May 2026).

1. Principal meanings of the folded-chain construction

The expression does not refer to a single canonical Hamiltonian. Instead, it appears in several technically precise settings. In one classical setting, a one-dimensional Ising model with $m$-spin interactions in a field is mapped to a zero-field rectangular Ising model of size $m\times N/m$ with the topology of a cylinder and helical boundary conditions (Turban, 2016). In another, a classical antiferromagnetic Ising chain in a field is dualized to an $m$-leg Ising tube, where rarefied string operators become ordinary tube spins (Timonin et al., 2017). In a chain-array construction, parallel Ising chains are coupled through a self-consistent field and iterated through exact one-dimensional renormalization-group flows to produce sharp two- and three-dimensional phase transitions (Kecoglu et al., 2020). In quantum settings, a transverse-field Ising chain is folded into a macroscopic topological qubit or into boundary sectors equivalent to one or two Kitaev chains (Zhang et al., 2020, Ma et al., 10 May 2026).

Construction	Starting object	Folded or equivalent object
Multispin chain in field	1D Ising chain with $m$-spin interactions	Rectangular Ising model on a cylinder with helical BC
Rarefied-string duality	Classical AF Ising chain in a field	$m$-leg Ising tube
Chain-array RG	Parallel exactly solved Ising chains	Effective 2D or 3D magnetization curves
Topological TFIC reduction	Open transverse-field Ising chain	Effective pseudospin or Kitaev/SSH sectors

A common structural feature is that local spin variables on the original chain are replaced by nonlocal variables, replicated sectors, or block degrees of freedom. This suggests that a folded Ising chain is best understood as a representation scheme rather than a uniquely fixed model.

2. Classical foldings into cylinders and tubes

A literal folded-chain construction is given by the one-dimensional Ising model with multispin interactions and longitudinal field

$$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$

with $\sigma_k=\pm 1$. For free boundary conditions and $N=mp$, the nonlocal change of variables

$\tau_k=\prod_{i=k}^N \sigma_i$

transforms the Hamiltonian into

$m$0

The bulk part is a rectangular Ising model with first-neighbour interactions $m$1 and $m$2, size $m$3, and the topology of a cylinder with helical boundary conditions. In the thermodynamic limit $m$4, $m$5, a two-dimensional critical singularity develops on the self-duality line $m$6 (Turban, 2016). At $m$7, the same model decomposes into $m$8 sublattices: for free boundary conditions,

$m$9

with correlation length

$m\times N/m$0

The folded geometry is therefore already encoded in the correlation structure.

A second classical realization starts from the antiferromagnetic Ising chain in a uniform field,

$m\times N/m$1

and introduces rarefied string operators

$m\times N/m$2

With shifted strings $m\times N/m$3, one defines tube spins

$m\times N/m$4

which map the chain exactly onto an $m\times N/m$5-leg Ising tube with nearest-neighbour couplings along the legs and four-spin plaquette interactions on neighbouring legs: $m\times N/m$6 The rarefied string correlation function becomes an ordinary two-point correlator along a tube leg,

$m\times N/m$7

For every odd $m\times N/m$8, there is a disorder line on which the asymptotics of $m\times N/m$9 changes from monotonic exponential decay to exponentially damped incommensurate oscillations, and these disorder lines are tied to Lee–Yang zeros in an $m$0-periodic complex magnetic field (Timonin et al., 2017). This is the most explicit realization of folding as an exact conversion of a nonlocal observable on a chain into a local observable on a multi-leg geometry.

3. Coupled and decorated classical chains

A different use of folding appears in chain-array reconstructions of higher-dimensional criticality. Starting from exactly solved one-dimensional Ising chains with Hamiltonian

$m$1

one restricts to chains along one spatial direction and writes the bond form

$m$2

Each isolated chain has no finite-temperature phase transition and no spontaneous magnetization for finite $m$3, but lateral couplings are represented at the start of each RG step by the self-consistent field

$m$4

with $m$5 in $m$6 and $m$7 in $m$8. Exact one-dimensional decimation yields the recursion relations for $m$9, and the chain magnetization is propagated through the density recursion

$m$0

Iterating the map $m$1 reconstructs sharp two- and three-dimensional magnetization curves and phase boundaries from coupled chains alone. The method produces approximate order-parameter exponents $m$2 in $m$3 with $m$4 and $m$5 in $m$6 with $m$7, while the phase boundaries compare well with the exact anisotropic $m$8 line $m$9 and with $m$0 Monte Carlo or RG estimates for the isotropic case (Kecoglu et al., 2020). Here folding means stacking or layering one-dimensional building blocks into a higher-dimensional cooperative magnet, with the exact chain RG determining whether inter-chain fields are amplified or suppressed.

Decoration produces another classical folded-chain geometry in the Ising model on the Toblerone lattice, a decorated two-leg ladder whose unit cell contains two leg spins $m$1 and one top spin $m$2,

$m$3

Tracing out the top spins yields the temperature-dependent effective rung coupling

$m$4

The transfer matrix has four nonzero eigenvalues $m$5, and a crossing of the sub-leading eigenvalues $m$6 and $m$7 produces a bifurcation of correlation lengths. Above the crossing, leg–leg correlations are governed by $m$8, whereas correlations involving a top spin are governed by $m$9, because symmetry forces the $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$0 prefactor to vanish in channels containing $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$1 (Chapman et al., 2024). This is a folded or decorated chain in the strict transfer-matrix sense: still one-dimensional, but with multiple internal degrees of freedom per cell and frustration-driven multiple length scales.

4. Quantum transverse-field foldings and topological sectors

The baseline quantum chain is the transverse-field Ising model

$$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$2

which is exactly solvable by the Jordan–Wigner transformation to noninteracting spinless fermions. Its excitation gap closes linearly at the quantum critical point $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$3, with $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$4. In the quasi-one-dimensional ferromagnet CoNb$$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$5O$$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$6, ultra-low-temperature thermal conductivity identifies the paramagnetic gap through the onset temperature $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$7, and a linear extrapolation gives the critical transverse field $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$8 T; at $$-\beta \mathcal{H}_N[\{\sigma\}] = K \sum_k \sigma_k \sigma_{k+1}\cdots \sigma_{k+m-1} + H \sum_k \sigma_k,$$9 T, phonon scattering persists down to $\sigma_k=\pm 1$0 mK, consistent with gapless magnetic excitations (Dai et al., 2011). This exactly solved chain provides the reference point for more elaborate folded constructions.

In the topological phase $\sigma_k=\pm 1$1, an open transverse-field Ising chain supports a nonlocal fermionic operator $\sigma_k=\pm 1$2 satisfying, in the thermodynamic limit,

$\sigma_k=\pm 1$3

The ground-state doublet $\sigma_k=\pm 1$4 is therefore a protected low-energy qubit. Real-space renormalization makes the folding explicit: coupling an $\sigma_k=\pm 1$5-site chain to an additional spin produces the effective Hamiltonian

$\sigma_k=\pm 1$6

where $\sigma_k=\pm 1$7 acts in the two-dimensional subspace spanned by the $\sigma_k=\pm 1$8-site ground doublet. The entire chain is thus reduced to a single effective spin-$\sigma_k=\pm 1$9, and the authors show that the $N=mp$0-site ground and first-excited states can be generated adiabatically from those of the $N=mp$1-site chain by turning on a boundary coupling. The time-dependent “dynamic crystallization” protocol uses

$N=mp$2

and numerically prepares the robust quasidegenerate doublet with high fidelity from a product state of noninteracting spins (Zhang et al., 2020). In this quantum usage, folding is not geometric replication but compression of an interacting many-body chain into a topological boundary qubit.

5. Boundary-engineered and isospectral foldings

Boundary engineering supplies another precise meaning. For the periodic transverse-field Ising chain

$N=mp$3

a type-II open boundary condition removes both the end bond and the transverse fields on the end sites: $N=mp$4 This choice implies edge spin conservation,

$N=mp$5

so the Hilbert space splits into sectors labeled by boundary spins. In each sector, a domain-wall basis maps the chain exactly onto a Kitaev-chain-like fermion Hamiltonian

$N=mp$6

for $N=mp$7. The same model also maps, under the standard Jordan–Wigner transformation, to an SSH-like Majorana chain with missing onsite couplings at the ends. The result is a switch in the topological degeneracy region: for type I open boundaries, topological degeneracy occurs for $N=mp$8, whereas for type II it occurs for $N=mp$9. For type II, all eigenstates are twofold degenerate for $\tau_k=\prod_{i=k}^N \sigma_i$0 and fourfold degenerate for $\tau_k=\prod_{i=k}^N \sigma_i$1, and the paper attributes the switch to the gauge dependence of the winding number in the SSH chain (Ma et al., 10 May 2026). Folding here is accomplished by boundary-sector decomposition and by reinterpreting boundary conditions as distinct cuts of an SSH ring.

A spectral version of folding appears in the construction of Ising analogues for multispin free-fermionic chains. The target model is a standard quantum Ising chain with open boundaries and inhomogeneous couplings,

$\tau_k=\prod_{i=k}^N \sigma_i$2

which becomes a nearest-neighbour Majorana chain

$\tau_k=\prod_{i=k}^N \sigma_i$3

For a family of multispin models with free-fermionic spectra, the quasienergies are encoded in polynomials $\tau_k=\prod_{i=k}^N \sigma_i$4. After the change of variables

$\tau_k=\prod_{i=k}^N \sigma_i$5

the modified Euclidean algorithm reconstructs an antisymmetric tridiagonal matrix whose off-diagonal entries determine the couplings $\tau_k=\prod_{i=k}^N \sigma_i$6. The resulting inhomogeneous Ising chain is isospectral to the original multispin model. For the $\tau_k=\prod_{i=k}^N \sigma_i$7 Fendley case, the numerical phase diagram contains gapped phases separated by critical lines, and the order–disorder transitions depend on the parity of the total number of energy density operators in the Hamiltonian (Alcaraz et al., 2023). This is a folded-chain construction in the sense that extended multispin interactions are compressed into a nearest-neighbour Ising chain with engineered inhomogeneities.

6. Material realizations and conceptual status

Bond-decorated Ising chains are realized in Dy-based molecular magnets, where Dy$\tau_k=\prod_{i=k}^N \sigma_i$8 supplies the Ising degree of freedom. The required microscopic condition is a ground Kramers doublet close to $\tau_k=\prod_{i=k}^N \sigma_i$9, which gives $m$00 and $m$01, with a large enough crystal-field gap that exchange and Zeeman couplings do not significantly admix excited doublets. In the $m$02 chain, ab initio calculations yield $m$03 with $m$04, and the fitted exchange constants are

$m$05

The model is an Ising–Heisenberg chain in which each Dy–Dy bond is decorated by a Cu–Mo–Cu trimer (Heuvel et al., 2010). In the Dy$m$06Cr$m$07 ring, the fitted parameters are

$m$08

with $m$09 and $m$10 for Dy$m$11, again producing an exactly solvable decorated Ising structure (Heuvel et al., 2010). These systems show that folded or decorated chain constructions are not purely formal; they can be embedded in real crystal-field environments and solved with transfer matrices after the decorating cluster is diagonalized.

An algebraic notion of folding also exists at the level of operator maps on integrable spin chains. For the Heisenberg XXX and Inozemtsev hyperbolic chains, a folding map identifies sites $m$12 and $m$13 and produces half-line Hamiltonians with integrable boundary terms and twisted Yangian symmetry (Gomez et al., 2016). That work is not an Ising construction, but it provides a broader meaning of “folding” as an operator-level projection from a line to a half-line.

Taken together, these constructions support a precise but plural usage. Folding may mean a nonlocal variable transformation that turns multispin couplings into nearest-neighbour couplings on a cylinder, a duality that converts string observables into local tube spins, a self-consistent coupling of exactly solved chains into higher-dimensional arrays, a reduction of a quantum chain to a protected boundary qubit, or a boundary condition that splits an Ising chain into Kitaev or SSH sectors. This suggests that the folded Ising chain is best regarded as a family of exact or controlled chain-based reorganizations of Ising degrees of freedom, unified by the replacement of a strictly linear description by a multi-leg, decorated, higher-dimensional, or topological one.