Spin-1/2 Hexagonal-Plaquette Chain

• The paper demonstrates that the spin-1/2 hexagonal-plaquette chain exhibits exactly solvable Kitaev interactions with conserved Z₂ plaquette fluxes leading to Majorana fermion solutions.
• It reveals distinctive band structures with coexisting flat and dispersive Majorana bands and identifies a topological transition marked by zero-energy edge modes.
• The study contrasts the Kitaev model with a molecular realization featuring Ising-like anisotropy, highlighting effective pseudospin behavior, finite excitation gaps, and interchain-induced Néel order.

Searching arXiv for the cited papers and closely related work on hexagonal-plaquette and six-site spin chains. A spin-$\tfrac12$ hexagonal-plaquette chain is a quasi-one-dimensional quantum spin system whose repeating structural motif is a six-site hexagon, or an equivalent six-spin hexagonal plaquette, arranged in series along a single spatial direction. In the recent arXiv literature, the term refers most directly to two distinct but related settings: an exactly solvable Kitaev-like chain derived from a truncated honeycomb geometry, with bond-directional $x,y,z$ couplings, local $\mathbb Z_2$ plaquette integrals of motion, and an exact Majorana-fermion solution (Majumder, 18 Jun 2026); and a molecular realization of a spin-$\tfrac12$ hexagonal-plaquette chain with Ising-like anisotropy in $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$, whose low-energy physics maps onto an effective spin-$\tfrac12$ Ising-like chain with a finite excitation gap and weak interchain-coupling-driven Néel order (Yamaguchi et al., 10 Sep 2025). Across these realizations, the central theme is that hexagonal plaquette geometry produces nontrivial internal structure within each unit cell and supports low-energy phenomena not present in ordinary uniform Heisenberg chains.

1. Geometry and defining characteristics

The most direct theoretical realization is a quasi-one-dimensional Kitaev spin model built from a truncated honeycomb lattice, naturally described as an array of coupled hexagonal plaquettes connected in series (Majumder, 18 Jun 2026). The chain is obtained by retaining a single quasi-1D strip from the honeycomb Kitaev lattice, and the paper describes it as a “quasi-one-dimensional hexagonal chain layer embedded within the honeycomb geometry” and as an “array of coupled hexagonal plaquettes connected through anisotropic inter-plaquette links.” Each unit cell is one hexagonal plaquette with six spin-$\tfrac12$ sites labeled

$A,B,C,D,E,F.$

In that model, the microscopic variables are explicitly spin-$\tfrac12$,

$S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$

with $x,y,z$0 the Pauli matrices, and the analysis later sets $x,y,z$1 (Majumder, 18 Jun 2026). The geometry is quasi-1D because translational invariance survives only along one direction, while the transverse structure inside each repeat unit remains nontrivial and honeycomb-derived. For this reason, the model is described equivalently as a hexagonal plaquette chain, a truncated honeycomb strip, or a decorated 1D chain or ladder (Majumder, 18 Jun 2026).

An experimentally realized hexagonal-plaquette chain has also been reported in the molecular compound

$x,y,z$2

where local radical–Co trimers form linked hexagonal motifs and the chain is built from edge-sharing magnetic units (Yamaguchi et al., 10 Sep 2025). Each chain unit $x,y,z$3 contains one central effective Co spin $x,y,z$4 and two radical spins $x,y,z$5 and $x,y,z$6. The intramolecular bonds connect $x,y,z$7 to both radicals in the same unit, while the intermolecular couplings connect radicals on neighboring units along two parallel legs. This produces a quasi-1D chain of linked hexagons rather than a standard single spin chain or ladder (Yamaguchi et al., 10 Sep 2025).

The term should be distinguished from “hexamer” chains whose six-site periodicity is translational rather than geometric. The spin-$x,y,z$8 Heisenberg chain with hexamer modulation of exchange has a six-site repeating unit but is not a chain of actual geometric hexagonal plaquettes (Naseri et al., 2011). Likewise, a plaquette-centered spin chain derived from truncated SU(2) lattice gauge theory can be a chain of square plaquettes rather than hexagons (Yao, 2023). These distinctions are structurally important because connectivity determines the local integrals of motion, low-energy cluster states, and admissible topological or fractionalized excitations.

2. Kitaev hexagonal-plaquette chain

The exactly solvable Kitaev realization is defined by the Hamiltonian

$x,y,z$9

Here $\mathbb Z_2$0 labels plaquettes, $\mathbb Z_2$1 label the six sites within a plaquette, $\mathbb Z_2$2 are bond-dependent Kitaev couplings, and the final term $\mathbb Z_2$3 is the inter-plaquette link (Majumder, 18 Jun 2026). The bond sequence around each hexagon is $\mathbb Z_2$4, and neighboring hexagons are coupled by an extra $\mathbb Z_2$5-bond from $\mathbb Z_2$6 to $\mathbb Z_2$7. A common parameter choice is

$\mathbb Z_2$8

or simply $\mathbb Z_2$9 when $\tfrac12$0 (Majumder, 18 Jun 2026).

A defining property of the model is the existence of one conserved plaquette operator per hexagon,

$\tfrac12$1

with $\tfrac12$2 (Majumder, 18 Jun 2026). For a chain with $\tfrac12$3 plaquettes, there are $\tfrac12$4 independent conserved $\tfrac12$5 variables, splitting the Hilbert space into $\tfrac12$6 plaquette-flux sectors. The $\tfrac12$7 sector is the flux-free or positive-plaquette sector, and $\tfrac12$8 is the negative-plaquette sector. A notable difference from the full honeycomb Kitaev model is that, in this truncated chain, the plaquettes are independent and there are no larger composite loops, so flux is localized to each hexagon (Majumder, 18 Jun 2026).

The exact solution follows the standard Kitaev spin-to-Majorana representation,

$\tfrac12$9

with four Majorana fermions per spin. The physical spin subspace is enforced by

$$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$0

with the physical sector given by $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$1 (Majumder, 18 Jun 2026). The bond variables

$$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$2

act as static $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$3 gauge fields, and the interacting spin model reduces exactly, in each fixed gauge sector, to a quadratic Majorana Hamiltonian of free itinerant Majorana fermions hopping on the six-site unit-cell chain (Majumder, 18 Jun 2026).

A geometry-specific feature is that not all Majorana flavors become itinerant. Only the sites $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$4 and $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$5, associated with inter-plaquette connectivity, have all relevant flavors participating dynamically; sites $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$6 each retain one extra Majorana flavor that does not enter the Hamiltonian. Per plaquette there are four such non-dynamic Majoranas, equivalent to two complex fermions, giving an extra degeneracy

$$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$7

for $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$8 plaquettes (Majumder, 18 Jun 2026). The ground state is taken to lie in the flux-free sector $$(\mathrm{p\!-\!Py\!-\!V})_2[\mathrm{Co}(\mathrm{NO}_3)_2]$$9 for all $\tfrac12$0, which permits a gauge choice with intra-plaquette bond variables fixed to $\tfrac12$1 (Majumder, 18 Jun 2026).

3. Band structure, topology, and localized modes

Because the unit cell has six sites, the flux-free Majorana problem has six bands. In the uniform positive-plaquette sector, the spectrum consists of four dispersive bands and two flat bands (Majumder, 18 Jun 2026). For

$\tfrac12$2

the flat bands occur at

$\tfrac12$3

and the coexistence of flat and dispersive Majorana bands is a central feature of the hexagonal plaquette geometry (Majumder, 18 Jun 2026).

The topological analysis uses the bipartition

$\tfrac12$4

under which the Bloch Hamiltonian has chiral off-diagonal form

$\tfrac12$5

The system belongs to the BDI symmetry class, with chiral symmetry made manifest by this block structure, and inversion symmetry is also invoked to quantize the Zak phase (Majumder, 18 Jun 2026). The winding number is

$\tfrac12$6

with $\tfrac12$7 in the trivial phase and $\tfrac12$8 in the topological phase (Majumder, 18 Jun 2026).

For the symmetric case $\tfrac12$9, the physically important topological transition occurs at

$\tfrac12$0

or equivalently $\tfrac12$1 (Majumder, 18 Jun 2026). The phase structure is a trivial gapped phase for $\tfrac12$2 except for the isolated $\tfrac12$3 gapless limit, a critical point at $\tfrac12$4, and a topological or chiral phase for $\tfrac12$5 with protected zero-energy edge states (Majumder, 18 Jun 2026). Under open boundary conditions, for $\tfrac12$6, two isolated zero-energy states appear in the single-particle spectrum and their wavefunctions are exponentially localized at the two ends of the chain. These are Majorana-like edge modes protected by chiral symmetry (Majumder, 18 Jun 2026).

The model also supports a second class of low-energy states tied to plaquette-flux structure rather than global topology. Introducing a contiguous domain of negative plaquettes inside a positive background produces two zero modes localized at the domain boundaries (Majumder, 18 Jun 2026). A single flipped plaquette already produces isolated near-zero or zero modes. Two or more adjacent negative plaquettes still produce only two zero modes, localized at the left and right boundaries of the negative domain, and multiple disconnected negative domains generate multiple pairs of near-zero modes, one pair per domain (Majumder, 18 Jun 2026). These domain-wall plaquette modes are distinct from the topological chain-end edge modes: they exist in the trivial regime $\tfrac12$7 and disappear in the topological phase, where low-energy physics is instead dominated by extended bulk states together with the chain-end Majorana modes (Majumder, 18 Jun 2026).

4. Fractional excitations and dynamical signatures

In the Kitaev hexagonal-plaquette chain, “fractional excitations” have the standard Kitaev meaning: a spin flip fractionalizes into itinerant Majorana fermions and $\tfrac12$8 flux or plaquette excitations (Majumder, 18 Jun 2026). Consequently, the spin dynamics is not described by sharp magnon poles as in conventional ordered magnets or simple Heisenberg chains, but by broad multiparticle continua (Majumder, 18 Jun 2026).

The single-spin dynamical response studied in the model is

$\tfrac12$9

It was computed numerically by DMRG for a 24-site chain with open boundary conditions (Majumder, 18 Jun 2026). For representative couplings $A,B,C,D,E,F.$0, $A,B,C,D,E,F.$1, and $A,B,C,D,E,F.$2, the response exhibits broad continua in all regimes and no sharp magnon branch (Majumder, 18 Jun 2026). In the trivial phase $A,B,C,D,E,F.$3, most spectral weight lies near $A,B,C,D,E,F.$4 at low but finite frequency. At the critical point $A,B,C,D,E,F.$5, the bulk gap closing leads to very low-energy spectral weight near zero frequency, with a small residual gap attributed to finite size. In the topological phase $A,B,C,D,E,F.$6, the response is essentially gapless with the strongest weight around $A,B,C,D,E,F.$7, and momentum weight is spread over several $A,B,C,D,E,F.$8-values rather than concentrated only at $A,B,C,D,E,F.$9; this low-energy weight is associated with edge Majorana modes (Majumder, 18 Jun 2026).

These signatures distinguish the system from ordinary spin-$\tfrac12$0 Heisenberg chains. The experimentally relevant fingerprints proposed for a quasi-1D hexagonal-plaquette Kitaev chain are broad dynamical continua rather than conventional spin-wave modes, additional low-energy or zero-frequency spectral weight in the topological phase from edge Majoranas, and signatures of localized flux-domain or plaquette modes near zero energy (Majumder, 18 Jun 2026). A plausible implication is that spectroscopy in such systems is expected to probe both the internal plaquette-flux sector structure and the 1D topological boundary physics.

5. Ising-like hexagonal-plaquette chain in a molecular material

A distinct realization of a spin-$\tfrac12$1 hexagonal-plaquette chain with Ising-like anisotropy has been reported in

$\tfrac12$2

(Yamaguchi et al., 10 Sep 2025). The magnetic building block contains two verdazyl radical spins, each carrying spin $\tfrac12$3, and one Co$\tfrac12$4 ion whose low-temperature Kramers doublet is treated as an effective spin $\tfrac12$5. The dominant exchanges are an intraplaquette or intramolecular ferromagnetic coupling

$\tfrac12$6

from molecular-orbital calculation, with the caution that the actual $\tfrac12$7 is likely about half of that, and an intermolecular antiferromagnetic coupling

$\tfrac12$8

between radical spins along the chain (Yamaguchi et al., 10 Sep 2025). The crystal is orthorhombic, with space group $\tfrac12$9 and lattice constants

$S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$0

and the dominant chain direction is the $S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$1-axis (Yamaguchi et al., 10 Sep 2025).

The microscopic Hamiltonian is given as

$S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$2

with $S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$3 the Ising-like anisotropy, $S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$4 the radical $S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$5-value, and $S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$6 the Co $S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$7-tensor (Yamaguchi et al., 10 Sep 2025). The anisotropy originates from the spin-orbit-coupled low-energy Kramers doublet of high-spin Co$S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$8 in a distorted octahedral environment. From ESR and single-ion analysis, the reported principal values are

$S^\alpha=\frac{\hbar}{2}\sigma^\alpha,$9

and the projected exchange anisotropy is estimated as

$x,y,z$00

(Yamaguchi et al., 10 Sep 2025).

The low-energy theory is built in the regime

$x,y,z$01

where the system is viewed as weakly coupled trimers $x,y,z$02 (Yamaguchi et al., 10 Sep 2025). Projecting onto the lowest trimer doublet yields an effective pseudospin-$x,y,z$03 chain in variables $x,y,z$04. To second order in $x,y,z$05,

$x,y,z$06

while at third order one obtains an XXZ-type interaction with a small transverse component,

$x,y,z$07

with

$x,y,z$08

(Yamaguchi et al., 10 Sep 2025). The effective chain is therefore Ising-like, with dominant antiferromagnetic Ising coupling and much weaker transverse exchange. The gap in the effective model is described by

$x,y,z$09

(Yamaguchi et al., 10 Sep 2025).

Experimentally, the specific heat at zero field shows a sharp peak at

$x,y,z$10

indicating a transition to an antiferromagnetically ordered state, together with a broader rounded feature near $x,y,z$11 K attributed to 1D short-range correlations and a low-temperature Schottky-like shoulder below $x,y,z$12, suggesting a finite gap in the ordered phase (Yamaguchi et al., 10 Sep 2025). Fitting the low-temperature heat capacity with

$x,y,z$13

gives

$x,y,z$14

(Yamaguchi et al., 10 Sep 2025). The transition at $x,y,z$15 K is attributed to weak but finite interchain interactions, and a coupled-chain estimate gives

$x,y,z$16

with the explicit caveat that this should be regarded only as an indicative energy scale rather than a precise microscopic parameter (Yamaguchi et al., 10 Sep 2025).

6. Relation to adjacent one-dimensional models and recurrent misconceptions

The spin-$x,y,z$17 hexagonal-plaquette chain is often conflated with other six-site or plaquette-chain constructions, but the distinctions are substantive. The hexamer-modulated Heisenberg chain has a six-site unit cell produced by bond alternation and period-3 modulation on strong bonds, not by actual hexagonal connectivity (Naseri et al., 2011). Its principal results concern strong-coupling dimer mapping, an effective XXZ chain in uniform and spatially modulated fields, and magnetization plateaus at $x,y,z$18 and $x,y,z$19 of saturation within a seven-phase field-induced phase diagram (Naseri et al., 2011). These results are useful as background for six-site commensurability physics, but they do not define a geometric hexagonal-plaquette chain.

A second common confusion concerns plaquette-centered pseudospin chains from lattice gauge theory. The SU(2) plaquette-chain construction obeying ETH is built from a chain of square plaquettes, not hexagons, and its effective spin-$x,y,z$20 variables live on plaquette centers rather than on the vertices of a hexagonal graph (Yao, 2023). Its effective Hamiltonian contains a nearest-neighbor Ising term, a longitudinal field, and a three-site dressed transverse-flip term,

$x,y,z$21

which is structurally informative for constrained plaquette dynamics but square-specific in its coefficients and local recoupling structure (Yao, 2023).

The diamond chain is likewise a close analogue rather than a hexagonal case. It is a frustrated three-site-unit-cell plaquette chain built from repeated diamond units, and its Jordan–Wigner treatment shows that gauge or string factors are indispensable on nontrivial quasi-1D graphs (Verkholyak et al., 2010). That model exhibits a robust $x,y,z$22 plateau, strong-dimer regimes, and field-induced jumps, but these features are tied to a three-site unit cell and diamond connectivity rather than a six-site hexagonal plaquette (Verkholyak et al., 2010).

These comparisons clarify two recurrent misconceptions. First, “hexamer” and “hexagonal plaquette” are not interchangeable terms: the former denotes period-6 translational structure, whereas the latter implies a literal six-site loop geometry (Naseri et al., 2011). Second, not every plaquette chain with spin-$x,y,z$23 variables is a hexagonal-plaquette chain; square plaquette, diamond plaquette, and honeycomb-derived hexagonal plaquette systems realize different local combinatorics, different constrained subspaces, and different low-energy excitations [(Yao, 2023); (Verkholyak et al., 2010); (Majumder, 18 Jun 2026)].

Within current arXiv work, the most direct meanings of a spin-$x,y,z$24 hexagonal-plaquette chain are therefore a quasi-1D chain of coupled hexagonal plaquettes with bond-dependent Kitaev interactions and conserved plaquette fluxes (Majumder, 18 Jun 2026), or a quasi-1D chain of linked hexagonal magnetic units with strong Ising-like anisotropy and an effective low-energy pseudospin description (Yamaguchi et al., 10 Sep 2025). In the first case, the characteristic phenomena are exact solvability, Majorana fractionalization, flat and dispersive bands, a BDI-class topological transition, and zero-energy boundary or domain-wall modes. In the second, the defining features are trimer-based effective spins, easy-axis anisotropy, a finite excitation gap, and interchain-coupling-driven Néel order.