Kitaev-Dimer Phase: Model-Dependent Insights

• Kitaev-Dimer Phase is a model-dependent dimerized regime in Kitaev systems characterized by distinct local dimer formations, including interlayer singlets, isolated dimers, and SSH-like topological modes.
• It appears in varied settings such as bilayer spin-½ models, generalized honeycomb configurations, and spin-1 chains, with its properties revealed through techniques like high-order series expansions and exact diagonalization.
• The phase exhibits diverse transition behaviors—from first-order and second-order Ising transitions to triplon condensation—highlighting its sensitivity to geometry, anisotropy, and interaction strengths.

The Kitaev-Dimer phase is a model-dependent dimerized regime that appears in several classes of Kitaev and Kitaev-type systems. In the most specific usage, it denotes the dimer singlet phase of the bilayer spin-$\frac12$ Kitaev model, where the ground state is a product of interlayer singlets and competes with the Kitaev quantum spin liquid (Koga et al., 2018). In other works, the same label or closely related usage refers to an Abelian gapped phase in the isolated-dimer limit of generalized honeycomb Kitaev models (Kamfor et al., 2010), a spontaneously dimerized phase in spin-1 Kitaev chains (Zhang et al., 2023), or an SSH-like topological dimer phase in interacting dimerized Kitaev superconductors (Wang et al., 2017). The term therefore does not denote a single universal phase across all models; rather, it designates dimer-dominated phases that arise from bond-directional Kitaev interactions in distinct microscopic settings.

1. Terminology and scope

In the cited literature, the phrase “Kitaev-dimer” is used in multiple, non-equivalent senses. The common element is the dominance of a dimerized local structure over the undimerized Kitaev regime, but the microscopic content of the dimer varies: interlayer singlets in bilayers, effective isolated dimers in perturbative limits, bond-alternating order in spin-1 chains, or edge-state-carrying dimerization in 1D superconducting chains.

Context	Defining object	Characterization
Bilayer spin-$\frac12$ Kitaev model	Interlayer singlet dimers	Product of interlayer singlets
Generalized honeycomb / hyperbolic Kitaev models	Isolated dimers on strong bonds	Abelian gapped phase governed by plaquette operators
Spin-1 Kitaev chains	Alternating strong and weak bonds	Spontaneous dimerization and broken lattice symmetry
Dimerized Kitaev superconductors	Alternating hopping/pairing/interactions	SSH-like or topological dimer phase with edge zero modes

This model dependence is essential. In some systems the Kitaev-Dimer phase is a featureless, trivial quantum paramagnet, while in others it is an Abelian topological phase or a symmetry-broken dimer crystal. A plausible implication is that the term is best interpreted as a structural descriptor tied to the dominant low-energy degrees of freedom, not as a universal phase label.

2. Bilayer spin-$\frac12$ Kitaev model: the canonical dimer-singlet usage

The bilayer Kitaev model is defined by

$$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$

where $J_K>0$ is the ferromagnetic intralayer Kitaev coupling and $J_H>0$ the antiferromagnetic interlayer Heisenberg coupling. At $J_H=0$, the system consists of two independent Kitaev layers, while at $J_K=0$ it is a product of interlayer dimer singlets (Koga et al., 2018).

The dimer expansion is performed from the strong-interlayer-coupling limit $J_K/J_H\ll 1$, treating the intralayer Kitaev interaction as a perturbation. A decisive technical ingredient is the existence of global parity symmetries in the singlet/triplet sectors. These imply that only clusters with even numbers of triplet excitations contribute, so the ground-state energy and interlayer correlation admit expansions

$$\frac{E_g}{N} = \sum_{i=0}^{30} a_i \left( \frac{J_K}{J_H} \right)^i,\qquad \left\langle {\bf S}_1\cdot{\bf S}_2 \right\rangle = \sum_{i=0}^{30} b_i \left( \frac{J_K}{J_H} \right)^i,$$

with only even powers present. This permits calculations up to 30th order in $\frac12$0 (Koga et al., 2018).

The resulting series analysis, supplemented by Padé and first-order inhomogeneous differential extrapolations, shows that the dimer singlet state is realized over a wide parameter region, specifically

$\frac12$1

In this regime, the dimer expansion agrees with exact diagonalization for both the ground-state energy and the interlayer spin-spin correlation. For $\frac12$2, the ground state cannot be captured by the dimer expansion, and exact diagonalization indicates a first-order transition to the quantum spin liquid phase (Koga et al., 2018).

Earlier exact-diagonalization, bond-operator mean-field, and cluster-expansion work reached a consistent conclusion: the bilayer model exhibits a first-order quantum phase transition between the Kitaev QSL and singlet-dimer states at $\frac12$3, and one-triplet excitations in the singlet-dimer regime are localized because of a local conserved quantity $\frac12$4 on each bi-hexagon (Tomishige et al., 2017). Exact diagonalization also found that for antiferromagnetic interlayer coupling the transition is signaled by a peak in $\frac12$5, a rapid drop in $\frac12$6, and a sharp change in the low-energy dynamical spin structure factor near $\frac12$7, whereas for ferromagnetic interlayer coupling no singularity is seen and the $\frac12$8 QSL connects smoothly to an $\frac12$9 Kitaev QSL (Tomishige et al., 2019).

3. Bilayer generalizations and interlayer valence-bond phases

The broader bilayer literature shows that the dimer phase is not unique to the pure bilayer Kitaev Hamiltonian and that its realization depends sensitively on stacking geometry, anisotropy, and additional exchange terms.

In stacked bilayer Kitaev models with different registries, increasing $\frac12$0 destroys the Kitaev spin liquid in favor of a paramagnetic dimer phase. Majorana-fermion mean-field theory, expansion techniques, and effective low-energy mappings show that the phase diagrams depend strongly on stacking and anisotropy. In AA stacking at strong anisotropy, the KSL-to-dimer transition is captured by a dual pseudo-spin Ising model and is second order in the $\frac12$1D Ising universality class, with critical scaling $\frac12$2, whereas isotropic cases can display a direct transition at $\frac12$3 and may involve weakly first-order behavior or more intricate excitation condensation scenarios (Seifert et al., 2018).

An explicit bilayer Kitaev-Heisenberg calculation with large-scale iPEPS and high-order series expansions found a valence bond solid state in a relatively narrow parameter region between the AFM and stripy phases. This state is adiabatically connected to isolated Heisenberg dimers, has vanishing local magnetization, preserves translational symmetry, and is characterized by strong interlayer rung correlations and weak, uniform intralayer correlations. It is reported as a phase that appears only in the bilayer model and is absent in the monolayer Kitaev-Heisenberg system (Samimi et al., 2024).

These results establish two important points. First, in bilayer settings the Kitaev-Dimer phase is often an interlayer-rung-singlet state rather than an in-plane dimer crystal. Second, the route into and out of the dimer phase is not universal: it may be first order, continuous, or preempted by other phases such as macro-spin phases or a flux phase with spontaneous interlayer coherence (Seifert et al., 2018).

4. Isolated-dimer limits, Abelian topological phases, and effective plaquette theories

A different usage of “Kitaev-Dimer phase” arises in Kitaev models analyzed around the isolated-dimer limit. In generalized honeycomb models built from arbitrary dimer coverings satisfying the trivalent matching rule, the Abelian gapped phase at $\frac12$4 is described perturbatively by

$\frac12$5

where the $\frac12$6 are conserved $\frac12$7 plaquette operators. In this context the “Kitaev-Dimer phase” is a gapped quantum spin liquid with Abelian anyons, not a trivial product state. Its detailed vortex properties depend strongly on the dimer covering; in covering III, for example, one- and two-vortex gaps depend on whether the effective plaquette is triangular or hexagonal, and vortex-vortex interactions can be either attractive or repulsive (Kamfor et al., 2010).

On regular hyperbolic trivalent tilings, the isolated-dimer limit again yields effective Hamiltonians built from plaquette variables $\frac12$8, now valid for arbitrary polygon length $\frac12$9. For the Kitaev and Kekulé colorings, the resulting dimer phase is described as adiabatically connected to the toric code phase with Chern number $$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$0, i.e. an Abelian topological phase supporting Abelian anyons. In the $$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$1 Bethe-lattice limit, by contrast, the gapped phase is topologically trivial (Vidal et al., 22 Jun 2025).

The parent-Hamiltonian construction that unifies kagome dimer models, ruby-lattice spin liquids, and the Kitaev honeycomb model gives a further perspective. Its weak-field limit reproduces kagome dimer physics, while strong fields project the model into either the spin-$$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$2 Kitaev honeycomb Hamiltonian or a spin-1 quadrupolar Kitaev model. The paper states that the “Kitaev-dimer phase” emerges via anyon fluctuations, and that the phase remains a $$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$3 spin liquid adiabatically connected to the dimer liquid under a nonlocal mapping to the kagome transverse-field Ising model (Verresen et al., 2022).

A useful contrast is provided by the decorated honeycomb Kitaev-type model in the isolated dimer limit. There, the effective Hamiltonian

$$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$4

has a ground state that is exactly a chiral spin liquid with spontaneous breaking of time-reversal symmetry, rather than a simple dimer paramagnet (Nasu et al., 2015). This shows that an isolated-dimer expansion in a Kitaev-type model need not imply a trivial dimer phase.

5. One-dimensional Kitaev-dimer phases in spin-1 chains

In spin-1 Kitaev chains, the Kitaev-dimer phase typically denotes a spontaneously dimerized phase with broken lattice symmetry. For the spin-1 Kitaev chain with uniaxial single-ion anisotropy,

$$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$5

the ground state in the flux-free sector maps exactly to a detuned PXP model. Infinite time-evolving block decimation finds a quantum phase transition from the Kitaev spin liquid to a dimer phase at

$$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$6

for $$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$7. The dimer phase is identified by the order parameter

$$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$8

and by spontaneous breaking of translational symmetry. The transition is reported as second order (Zhang et al., 2023).

With general single-ion anisotropies, iTEBD reveals a phase diagram containing the Kitaev spin liquid, gapless dimer phases, and ferroquadrupole phases. The KSL-to-dimer transition driven by uniaxial SIA is described as analogous to the confinement–deconfinement transition in the lattice Schwinger model with topological $$\mathcal{H}= - J_K\sum_{\langle ij\rangle_\alpha,n} S_{i,n}^\alpha S_{j,n}^\alpha +J_H \sum_{i} {\bf S}_{i,1}\cdot{\bf S}_{i,2},$$9-angle $J_K>0$0, while rhombic SIA shifts the effective $J_K>0$1-angle away from $J_K>0$2 and can replace the critical line by a crossover or a different transition structure (Zhang et al., 2024).

A distinct spin-1 realization occurs in the bilinear-biquadratic-Kitaev chain,

$J_K>0$3

where DMRG finds a Kitaev-dimer phase that is gapped, twofold degenerate, and distinguished by spontaneous breaking of the screw symmetry $J_K>0$4. It selects either $J_K>0$5- or $J_K>0$6-spin bonding and coexists with a crystalline order of alternating $J_K>0$7 fluxes, encoded by

$J_K>0$8

The string order parameter vanishes, indicating no SPT character, and the transition into the adjacent gapless quadrupolar phase is marked by gap closing and a change in central charge from $J_K>0$9 to $J_H>0$0 (Wei et al., 7 Aug 2025).

6. Dimerized Kitaev superconductors, diagnostics, and conceptual distinctions

In interacting dimerized Kitaev topological superconductors, the Kitaev-dimer label refers to a topological dimer phase rather than a magnetic dimer state. At the exactly solvable point $J_H>0$1, $J_H>0$2, dimerization is introduced through

$J_H>0$3

and the phase diagram contains seven distinct phases separated by gap closings at

$J_H>0$4

In this setting, the system is topological when a fermionic many-body Majorana zero-energy edge state emerges, and the “Kitaev-Dimer phase” denotes a topological dimer region, particularly for negative dimerization $J_H>0$5, where the phase is adiabatically connected to the Kitaev chain and carries Majorana edge modes (Ezawa, 2017).

A closely related exact solution uses two edge correlation functions,

$J_H>0$6

to distinguish the trivial phase, a topological superconductor, and an SSH-like topological phase. In the thermodynamic limit, the SSH-like phase has both correlators nonzero and hosts Dirac edge zero modes; the paper explicitly identifies it as a novel interacting analog of the SSH topological insulator (Wang et al., 2017).

Across the full body of work, the diagnostics of Kitaev-dimer phases are correspondingly diverse. In bilayer spin models they are identified by high-order dimer expansions, exact diagonalization, interlayer correlations, and local conserved quantities (Koga et al., 2018). In stacked and Kitaev-Heisenberg bilayers they are tracked by series expansions, triplon gaps, and tensor-network order parameters (Seifert et al., 2018). In spin-1 chains they are diagnosed by dimer order parameters, iTEBD or DMRG, entanglement structure, and flux or bond-parity operators (Zhang et al., 2023). In superconducting chains they are detected by exact edge correlators and many-body Majorana or Dirac zero modes (Wang et al., 2017).

Several common misconceptions are corrected by this comparison. The first is that a Kitaev-Dimer phase is always trivial; generalized honeycomb, hyperbolic, and parent-Hamiltonian constructions show that it can be Abelian topological (Kamfor et al., 2010). The second is that it always breaks translation symmetry; bilayer rung-singlet phases preserve translational symmetry, whereas spin-1 chain realizations can spontaneously dimerize (Samimi et al., 2024). The third is that it has a fixed critical behavior; the literature instead reports first-order transitions, second-order Ising transitions, triplon-condensation transitions, and model-specific crossovers depending on geometry, anisotropy, and the microscopic definition of the dimer (Tomishige et al., 2017).