Zigzag-Chain Model Overview

• The zigzag-chain model is a one-dimensional framework characterized by a non-straight, alternating geometry that introduces competing nearest-neighbor and next-nearest-neighbor interactions.
• It enables the study of frustrated spin Hamiltonians, orbital orders, and topological superconductivity using techniques such as DMRG and exact analytical solutions.
• The model’s versatility in capturing magnetic, electronic, and structural transitions under symmetry-breaking conditions offers a robust platform for exploring novel quantum phases.

The zigzag-chain model is a family of one-dimensional and quasi-one-dimensional theoretical constructions in which the underlying geometry is not a straight chain but a zigzag arrangement of sites or bonds. In the literature, this geometry is used to encode competing nearest-neighbor and next-nearest-neighbor processes, sublattice staggering, locally noncentrosymmetric environments, and bond-dependent anisotropies. As a result, zigzag chains appear in quantum spin systems, multi-orbital Hubbard models, magnetoelectric and multipolar settings, topological superconductors, exciton-polariton lattices, trapped-ion structural transitions, and colloidal chains. The common theme is that the zigzag geometry reorganizes otherwise standard couplings into frustrated, staggered, or symmetry-lowered forms, thereby stabilizing phases such as dimer singlets, nematic-dimer states, block antiferromagnets, incommensurate helices, magnetization plateaus, flat-band compactons, and Majorana sectors (Saito et al., 2024, Hayami, 2023, Kumar et al., 5 Jul 2026, Johansson et al., 2018, Shimshoni et al., 2010).

1. Geometry, unit cells, and symmetry reduction

A zigzag chain typically introduces at least a two-site structure, because consecutive bonds are not equivalent. In the exciton-polariton realization, the links make angles of $\pm \pi/4$ with respect to the chain axis, and the factors $(-1)^n$ in the generalized DNLS equations encode the alternating sublattice structure directly (Johansson et al., 2018). In magnetic materials, the same geometry is often described as a chain with nearest-neighbor bonds along the zigzag rungs and next-nearest-neighbor bonds along the chain legs, so that frustration is intrinsic rather than added perturbatively.

A central symmetry feature is the distinction between global and local inversion. In the electric-axial-moment problem, the zigzag chain preserves global inversion and time-reversal symmetries while breaking local inversion symmetry at each sublattice. The site symmetry is reduced from globally centrosymmetric $D_{2h}$ to locally polar $C_{2v}$, which permits odd-parity crystalline electric fields, sublattice-dependent hybridization between orbitals of opposite parity, and site-dependent antisymmetric SOC (Hayami, 2023). In that setting, the minimal microscopic odd-parity term is

$$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$

with $p_i=+1$ on one sublattice and $-1$ on the other.

The zigzag geometry also modifies exchange pathways through bond angles and ligand environments. For YbCuS$_2$, the derived quantum spin Hamiltonian is obtained by incorporating octahedral crystal field effects, atomic spin-orbit coupling, and strong Coulomb interactions on Yb ions, followed by fourth-order perturbation theory. A small but finite anisotropic exchange coupling called the $\Gamma$-term appears, and two factors enhancing it are identified: the splitting of excited $f$-states with two holes, and the tilting of octahedra or distortion of the S-Yb-S bond angle, which break the symmetry of the Slater-Koster overlap (Saito et al., 2023). In BaCoTe$(-1)^n$0O$(-1)^n$1, the zigzag arrangement likewise makes next-nearest-neighbor processes unusually important because the next-nearest-neighbor bond distance is close to the nearest-neighbor distance, allowing realistic NNN hopping to compete strongly with NN hopping (Lin et al., 2024).

2. Frustrated spin Hamiltonians and magnetic phases

In quantum magnetism, the canonical zigzag-chain construction is the frustrated $(-1)^n$2-$(-1)^n$3 spin-$(-1)^n$4 chain, often supplemented by anisotropic terms. A representative Hamiltonian is the Heisenberg-$(-1)^n$5 model

$(-1)^n$6

whose DMRG phase diagram contains Tomonaga-Luttinger-liquid, vector-chiral, dimer-singlet, nematic singlet-dimer, and Ising-type ferromagnetic or antiferromagnetic phases, together with multicritical points and several exactly solvable points (Saito et al., 2024). In that framework, a finite $(-1)^n$7 term does not immediately destroy the dimer singlet; instead it produces a nematic singlet-dimer phase with a nonzero spin gap and a substantial zero-energy peak in the nematic dynamical structure factor. The paper interprets this as dilute but robust concentration of nematicity on top of singlets on dimers.

Rare-earth zigzag chains supply material realizations of this anisotropic physics. In YbCuS$(-1)^n$8, high-resolution powder neutron diffraction finds an elliptic helical incommensurate magnetic structure with propagation vector $(-1)^n$9 along the zigzag chain and a magnetic moment at least one-third smaller than expected for the Yb$D_{2h}$0 Kramers-doublet ground state. Under magnetic field, an up-up-down state appears at $D_{2h}$1 T with $D_{2h}$2. These observations agree with DMRG calculations for a zigzag spin-$D_{2h}$3 model containing isotropic Heisenberg interactions and off-diagonal symmetric $D_{2h}$4-type exchange interactions, and the incommensurate magnetism is interpreted as remnant off-diagonal spin correlations in a nematic dimer-singlet state (Onimaru et al., 8 Jan 2025). The same material class motivated the microscopic derivation of $D_{2h}$5-type superexchange noted above (Saito et al., 2023).

Other rare-earth zigzag chains emphasize different limits. PrTiNbO$D_{2h}$6 realizes a zigzag pseudospin-$D_{2h}$7 antiferromagnetic chain in which random crystal electric fields from Ti$D_{2h}$8/Nb$D_{2h}$9 site mixing generate a non-Kramers quasi-doublet with Ising moment, $C_{2v}$0 and $C_{2v}$1. Despite antiferromagnetic intersite coupling of about $C_{2v}$2 K, no magnetic freezing is detected down to $C_{2v}$3 K, while a sizable zero-field gap of about $C_{2v}$4 K is ascribed to off-diagonal anisotropy terms in the pseudospin Hamiltonian (Li et al., 2018). TbTi$C_{2v}$5Bi$C_{2v}$6 represents a quasi-one-dimensional Ising magnet in which crystalline-electric-field effects align the Tb$C_{2v}$7 moments along the zigzag chain direction, and magnetic field along that direction produces multiple metamagnetic transitions between a $C_{2v}$8 plateau and other plateaus (Guo et al., 2024).

Frustrated zigzag chains also interpolate between one-dimensional fractionalization and three-dimensional order. In $C_{2v}$9-TeVO$$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$0, the excitation spectrum in the ordered ground state contains both magnon dispersion and a spinon-like continuum that prevails above $$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$1 meV. The minimal exchange-network model has ferromagnetic $$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$2 K, antiferromagnetic $$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$3 K, weaker interchain couplings $$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$4 K and $$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$5 K, and substantial exchange anisotropies. The observed spectrum is reproduced by combining linear spin-wave theory with pre-calculated spinon-continuum results, making the coexistence of magnons and spinons a direct signature of dominant intrachain frustration plus weak interchain ordering (Pregelj et al., 2018).

An exactly constructed integrable zigzag ladder furnishes a complementary perspective. Its Hamiltonian,

$$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$6

shares some thermodynamical properties with the Heisenberg XXZ chain but exhibits different ordering and critical exponents. For positive nearest-neighbor Ising coupling it supports two gapped phases, dimerized antiferromagnetic order and usual antiferromagnetic (Néel) order, separated by an extended critical region that is a quantum spin liquid with broken parity symmetry and oscillatory long-distance $$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$7 correlations (Tavares et al., 2023).

3. Orbital, multipolar, and correlated-electron zigzag chains

Zigzag geometry is equally important in electronic and multi-orbital models, where it restructures hopping networks and local parity. In BaCoTe$$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$8O$$\mathcal{H}^{\rm odd} = -V \sum_{i \sigma} p_i (c^\dagger_{i s \sigma} c_{i p_y \sigma} + \text{H.c.}),$$9, first-principles calculations identify strongly anisotropic one-dimensional Co $p_i=+1$0 bands near the Fermi level and motivate a three-orbital Hubbard model on the zigzag chain. The dominant NNN hopping is the $p_i=+1$1-$p_i=+1$2 matrix element $p_i=+1$3 eV, much larger than the corresponding NN value. DMRG then shows that a model with only NN hopping favors a staggered antiferromagnetic region, whereas the realistic model including NN and NNN hopping stabilizes a block AFM order $p_i=+1$4 and a Mott insulating state with three half-filled orbitals (Lin et al., 2024).

In the $p_i=+1$5 zigzag spin-orbital chain relevant to CaV$p_i=+1$6O$p_i=+1$7, a $p_i=+1$8-orbital Hubbard model with crystal-field splittings $p_i=+1$9 and $-1$0 produces orbital-state transitions that change the effective spin problem itself. Depending on orbital occupation, the resulting low-energy magnet can be an $-1$1 antiferromagnetic chain, a zigzag $-1$2 chain with ferromagnetic NN and antiferromagnetic NNN couplings, an antiferromagnetic ladder, or a non-magnetic singlet regime with double occupancy in the $-1$3 orbital (Onishi, 2013). In this class of models, orbital order is not a secondary feature; it selects the topology and sign structure of the exchange network.

Multipolar zigzag chains highlight a different use of local asymmetry. In the locally noncentrosymmetric four-orbital model, odd-parity crystalline electric fields generate effective cross-product couplings between the electric dipole and the electric toroidal dipole, the latter corresponding to the electric axial moment. A linear-response calculation shows that an external electric field along the chain induces a staggered electric axial moment, while staggered electric-dipole ordering induces a uniform electric axial moment. Uniform electric quadrupole ordering also accompanies a uniform electric axial moment, and the resulting uniform axial moment permits transverse magnetization responses such as $-1$4 and $-1$5 (Hayami, 2023).

Zigzag-chain order can also describe lattice distortions rather than spin or charge Hamiltonians. In LiVSe$-1$6, the low-temperature structure is best modeled in monoclinic $-1$7 and exhibits zigzag chain-like local distortions with finite correlation length rather than vanadium trimerization. The zigzag motif is competitive with the trimer order in the LiV$-1$8 series, persists away from the trimer phase-transition boundary, and is more stable in LiVSe$-1$9 than in LiVS$_2$0 in terms of the temperature variation of atomic displacement and correlation length (Kojima et al., 2023). This provides a structural counterpart to the more familiar spin and orbital zigzag models.

4. Exact, integrable, and topological realizations

Several zigzag-chain models are analytically tractable because the geometry can be absorbed into staggered or dimerized couplings. A prominent example is the spin-$_2$1 $_2$2 magnetoelectric chain on a zigzag lattice, where the Katsura-Nagaosa-Balatsky mechanism produces a staggered Dzyaloshinskii-Moriya interaction. The Hamiltonian is reduced by non-uniform spin rotations to a dimerized $_2$3 chain and then solved exactly through the Jordan-Wigner transformation (Baran et al., 2018). In this model, the electric field can open a gapped zero-plateau phase absent in the straight-chain analogue, the polarization need not align with the external electric field, and the electric field may enhance the magnetocaloric effect.

Topological superconductivity admits an equally direct zigzag formulation. The zigzag Kitaev chain is constructed from two diagonally coupled one-dimensional Kitaev chains with diagonal inter-chain hoppings $_2$4 and $_2$5 and phase difference $_2$6. In Bogoliubov-de Gennes form, the model supports distinct topological sectors indexed by the winding number

$_2$7

For zero phase difference, $_2$8, $_2$9, and $\Gamma$0 correspond respectively to phases with four, two, and zero Majorana zero modes; at $\Gamma$1, the degeneracy of the Majorana modes is partially lifted and the topological boundaries are modified (Kumar et al., 5 Jul 2026). Because the two constituent chains contribute independently to the total winding number, the zigzag geometry provides a minimal setting in which MZM multiplicity is tunable rather than fixed.

The integrable IRF zigzag ladder noted above belongs in the same exact-solvability tradition, although its physics is not topological in the superconducting sense. There the quantum transfer matrix and non-linear integral equations determine the free energy, correlation length, and oscillation momentum. The most distinctive feature is that the critical phase is parity-broken and the leading long-distance spin correlations oscillate with momentum $\Gamma$2 at $\Gamma$3 and $\Gamma$4, rather than with the $\Gamma$5 oscillations characteristic of the XXZ chain (Tavares et al., 2023). This makes the zigzag chain not merely an alternative parameterization of a standard one-dimensional model, but a distinct universality setting.

5. Bosonic and structural zigzag transitions

In bosonic lattice realizations, zigzag geometry can act as a synthetic spin-orbit texture. For weakly coupled exciton-polariton condensates, a tight-binding reduction of Gross-Pitaevskii-type equations yields generalized coupled DNLS equations in which the TE-TM coupling acquires alternating signs,

$\Gamma$6

because the links are at $\Gamma$7 relative to the chain axis (Johansson et al., 2018). The linear dispersion

$\Gamma$8

opens a gap of width $\Gamma$9 at $f$0, and at $f$1 becomes exactly flat. At that flat-band point, the model supports exact compactons localized on two neighboring sites with spin-up and spin-down parts $f$2 out of phase; away from it, these continue into exponentially localized gap modes.

A different class of zigzag-chain models concerns structural instabilities in trapped ions. At finite quench rate, reducing the transverse confinement of an ion chain with periodic boundary conditions drives a linear-to-zigzag phase transition, producing zigzag domains oriented along different transverse planes. The complex order parameter is

$f$3

and under periodic boundary conditions twists between domains can relax into helical chains with nonzero winding number (Nigmatullin et al., 2011). The Kibble-Zurek analysis predicts $f$4 and $f$5, which molecular-dynamics simulations confirm quantitatively in the underdamped mean-field regime.

At zero temperature, the same linear-zigzag instability becomes a quantum phase transition. For trapped ions or repulsively interacting particles, the effective transverse displacement field maps to a one-dimensional quantum Ising model in a transverse field,

$f$6

with a quantum critical point shifted from the classical instability by tunneling effects (Shimshoni et al., 2010). Large-scale DMRG on the corresponding discretized $f$7 chain extracts critical exponents $f$8, $f$9, $(-1)^n$00, and central charge $(-1)^n$01, thereby confirming Ising universality with high precision (Silvi et al., 2013).

Colloidal chains furnish a classical soft-matter analogue. Paramagnetic colloids trapped in a linear optical array form a zigzag pattern when an external magnetic field induces repulsive dipolar interactions. For harmonic traps, the infinite-chain equilibrium problem has a pitchfork bifurcation at

$(-1)^n$02

and normal-mode analysis shows that the mode first becoming unstable is the transverse zigzag mode at $(-1)^n$03 (Straube et al., 2013). When the traps are switched off, the same symmetry controls nonequilibrium expanding patterns, even below the field needed for the equilibrium zigzag state.

6. Coupled-chain extensions and conceptual scope

Zigzag chains also function as building blocks for higher-dimensional theories. In the nearest-neighbor spin-$(-1)^n$04 Kitaev-Heisenberg-Gamma model on the honeycomb lattice, strong bond anisotropy permits a coupled-chain description in which the decoupled limit is a gapless Luttinger liquid and the inter-chain couplings are treated by self-consistent mean field theory (Yang et al., 2022). In the AFM-$(-1)^n$05 region this procedure finds three ordered states—$(-1)^n$06, commensurate counter-rotating spiral, and zigzag—and the two first-order phase boundaries merge at $(-1)^n$07 and $(-1)^n$08, which is predicted to be a quantum critical point. The same analysis connects the spiral order of $(-1)^n$09-Li$(-1)^n$10IrO$(-1)^n$11 and the zigzag order of Na$(-1)^n$12IrO$(-1)^n$13 within a single minimal model.

Taken together, these results show that the zigzag-chain model is best understood as a geometrical platform rather than a single Hamiltonian. Depending on the microscopic degrees of freedom and symmetry constraints, the geometry can manifest as competing NN and NNN exchanges, staggered spin-orbit terms, bond-dependent magnetoelectric couplings, locally odd-parity hybridization, or a soft transverse structural mode. Accordingly, zigzag-chain systems are not restricted to one canonical phase: the same geometry supports dimer singlets, nematic singlet-dimer behavior, block antiferromagnetism, up-up-down plateaus, oscillatory critical correlations, electric axial moments, Majorana multiplicity, compactons, and helical structural defects (Saito et al., 2024, Lin et al., 2024, Onimaru et al., 8 Jan 2025, Tavares et al., 2023, Hayami, 2023, Kumar et al., 5 Jul 2026). A plausible implication is that zigzag geometry should be treated on the same footing as interaction strength or dimensionality when classifying quasi-one-dimensional phases, because in the systems surveyed here it is the geometry itself that converts otherwise conventional local ingredients into qualitatively new effective theories.