Kekulé-Patterned Lattice Overview

• Kekulé-patterned lattices are periodic modulations of a honeycomb structure that enlarge the unit cell and induce Brillouin-zone folding leading to valley mixing.
• They appear in varied textures—including Kekulé-O, Kekulé-Y, and Y-Kekulé distortions—that affect electronic band gaps, Dirac cone behavior, and topological properties.
• Experimental realizations in molecular graphene, graphene/TMD heterostructures, and Sn-intercalated systems validate these models and suggest novel avenues for quantum applications.

A Kekulé-patterned lattice is a lattice in which nearest-neighbor bonds, effective hoppings, densities, or analogous couplings are modulated periodically so that the original unit cell is enlarged—most often to a $\sqrt{3}\times\sqrt{3}$ supercell or another tripled-cell structure—and the original Brillouin zone is folded so that valleys such as $K$ and $K'$ are brought to the reduced-zone center. In honeycomb systems this modulation appears as alternating strong and weak bonds, as Kekulé-O or Kekulé-Y textures, or as Y-Kekulé distortions that preserve $C_{6v}$; in broader usage, related three-sublattice superstructures also occur as valence-bond order, exchange-field textures, and moiré reconstructions (Wu et al., 2015, Bergman, 2012, Jiang et al., 2024, Ye et al., 2023).

1. Structural definition and real-space motifs

The canonical Kekulé construction on a honeycomb lattice is a periodic modulation of nearest-neighbor couplings into two values. In one standard tight-binding form, the lattice is viewed as a triangular array of hexagons, with intra-hexagon hopping $t_0$ and inter-hexagon hopping $t_1$,

$$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$

so that six neighboring sites are grouped into “hexagonal artificial atoms” (Wu et al., 2015). In graphene-oriented formulations, the same logic is often described as a $\sqrt{3}\times\sqrt{3}$ bond modulation that enlarges the primitive cell from two sites to six (Bergman, 2012).

The real-space texture is not unique. In graphene/TMD heterostructures, four distinct textures are explicitly distinguished: Kekulé-Y-1, where two bonds out of six in each hexagon have high density; Kekulé-Y-2, where four bonds out of six are modulated in a different Y-like arrangement; Kekulé-O, where bonds form an O-shaped modulation; and Kekulé-M, where the bonds are equal but the on-site energies differ on the three sublattices (Zhang et al., 2022). This multiplicity is important because “Kekulé-patterned lattice” denotes a symmetry-related family of tripled-cell reconstructions rather than a single bond picture.

A particularly explicit honeycomb realization is the Y-Kekulé distortion. There, the stronger and weaker nearest-neighbor bonds are arranged in Y-shaped motifs, with half of the motifs pointing upward and the other half downward, so that the lattice remains globally sixfold symmetric. In the $p_z$-orbital tight-binding model the unit cell contains eight sites, and the bond modulation reorganizes the lattice into a network of hexamers and Y-shaped tetramers (Jiang et al., 2024). The corresponding Hamiltonian is

$$H = -t_0\sum_{\langle i,j\rangle}c^\dagger_i c_j - t_1\sum_{\langle i',j'\rangle}c^\dagger_{i'}c_{j'} .$$

2. Brillouin-zone folding, intervalley coupling, and order-parameter structure

The defining momentum-space consequence of a Kekulé pattern is Brillouin-zone folding. In graphene, a $K$0 modulation breaks the original translation symmetry, folds the two Dirac points into the same reduced-zone point, and introduces valley mixing (Bergman, 2012). In the language of ordered phases, the Kekulé bond-order parameter lives at wavevector $K$1, and its phase distinguishes the symmetry-related domains. A representative expression is

$K$2

with three inequivalent Kekulé branches related by $K$3 rotation, so the broken symmetry is effectively $K$4 (Xu et al., 2018).

The low-energy field theory is often written in terms of a complex Kekulé mass. In molecular graphene the three domains are parameterized as

$K$5

with relative phase shifts of $K$6 between domains (Bergman, 2012). In a lattice-gauge formulation of monolayer graphene, the external Kekulé modulation appears as a valley-mixing Dirac mass-like term $K$7, and for free fermions yields

$K$8

so the gap is proportional to $K$9 (Araki, 2011).

A common misconception is that every Kekulé texture necessarily opens a Dirac gap. The literature contains a clear counterexample. In Kekulé-Y graphene, the continuum Hamiltonian contains a valley-momentum-locking term, and the low-energy spectrum consists of two nested Dirac cones centered at $K'$0,

$K'$1

with distinct propagation velocities for the two isospin branches (Andrade et al., 2022). By contrast, Kekulé-O reconstructions on quasi-free-standing graphene or on CdS do open gaps through intervalley folding and mixing (Ngo et al., 23 Dec 2025, Betancur-Ocampo et al., 2024). The distinction is therefore texture-dependent, not terminological.

3. Topological phases, invariants, and boundary states

One major line of work treats the Kekulé-patterned honeycomb lattice as a topological medium. When six neighboring sites are grouped into hexagons and the inter-hexagon hopping exceeds the intra-hexagon hopping, the gap reopens at $K'$2, the $K'$3 and $K'$4 orbital sectors invert, and the system enters a $K'$5 topological phase characterized by an emergent pseudo time-reversal symmetry associated with $K'$6 symmetry (Wu et al., 2015). In that construction, the pseudospin is carried by the angular-momentum eigenstates of the hexagonal artificial atoms, and interface or edge spectra exhibit in-gap states with pseudo-Kramers structure.

The Y-Kekulé honeycomb lattice supports a distinct form of higher-order topology. The two coupling regimes are the expanded distortion $K'$7 and the shrunken distortion $K'$8. In the adiabatic limits, the expanded phase connects to isolated hexamers, while the shrunken phase connects to isolated Y-shaped tetramers (Jiang et al., 2024). The higher-order invariant is diagnosed by a Berry phase obtained from a local twist on one chosen hexamer; because of the $K'$9 symmetry of the hexamer, the phase is quantized as

$C_{6v}$0

which defines a $C_{6v}$1 index (Jiang et al., 2024).

In the expanded Y-Kekulé phase, the Berry phase remains quantized to $C_{6v}$2 for the first bulk gap below the Fermi level and to $C_{6v}$3 for the second as long as the bulk gap stays open. This nontrivial quantization reflects adiabatic connection to an isolated hexamer whose Wannier center lies at the center of the hexamer rather than at a trivial atomic position. The shrunken phase is topologically trivial because it is adiabatically connected to an isolated Y-shaped tetramer whose Wannier center sits at the central lattice site of the tetramer (Jiang et al., 2024).

Bulk-boundary correspondence in this setting is realized in finite flakes whose termination breaks the hexamers. For $C_{6v}$4, the expanded phase exhibits two triplets of in-gap states in the first gap and one triplet in the second. These corner modes are protected by the $C_{6v}$5 symmetry of the finite flake and are sharply localized at the corners; they are distinct from the edge states associated with weak first-order topology (Jiang et al., 2024). The literature therefore distinguishes at least three topological regimes within Kekulé-patterned lattices: $C_{6v}$6 pseudospin topology, higher-order $C_{6v}$7 topology, and weak first-order edge topology.

4. Correlation-driven Kekulé order and competing many-body phases

Kekulé patterning may be imposed externally, but it also arises spontaneously from interactions. In a cluster-charge extended Hubbard model on the honeycomb lattice, determinantal quantum Monte Carlo finds an intermediate Kekulé valence bond solid (KVBS) between the Dirac semimetal and the antiferromagnetic Mott insulator. The KVBS breaks translation symmetry, enlarges the unit cell to a $C_{6v}$8 supercell, and is described by an order parameter at $C_{6v}$9 whose histogram in the complex plane forms three clusters (Xu et al., 2018). The semimetal-to-KVBS transition is continuous, with a critical coupling $t_0$0 for $t_0$1, exponents $t_0$2 and $t_0$3, and numerically observed emergent $t_0$4 symmetry consistent with the chiral XY universality class (Xu et al., 2018). By contrast, the KVBS-to-AFMI transition appears first order.

A different interaction mechanism appears in the breathing kagome lattice near higher-order Van Hove singularities. There, anisotropic band flattening together with finite broadening yields approximate nesting at wavevector $t_0$5, producing a $t_0$6 bond-dominant density wave with some site-density modulation. Functional renormalization group identifies the dominant particle-hole instability at $t_0$7 and $t_0$8, and symmetry analysis classifies the order as the $t_0$9 irreducible representation of the extended point group $t_1$0 (Beck et al., 28 May 2025). In the line-graph reconstruction of kagome into an effective honeycomb structure, this ordered state becomes a Kekulé-Y pattern on the emergent honeycomb (Beck et al., 28 May 2025).

Kekulé ordering also competes with other symmetry-breaking channels. In strong-coupling lattice gauge theory for monolayer graphene, Coulomb interactions alone generate a nonzero chiral condensate $t_1$1, but an externally applied Kekulé distortion suppresses that order and restores chiral symmetry through a second-order transition at $t_1$2. The quasiparticle gap remains finite because the Kekulé term itself is gap-generating (Araki, 2011). In a different direction, nearest-neighbor attraction on the honeycomb lattice favors a Kekulé superconducting state whose bond order parameters oscillate with wavevector $t_1$3, producing a tripled unit cell; the semimetal first enters a p-Kekule superconductor and then undergoes a discontinuous transition into an s-Kekule state deeper in the superconducting phase (Roy et al., 2010). These examples establish that “Kekulé-patterned lattice” is not restricted to single-particle band engineering: it also names interaction-selected ordered states.

5. Experimental realizations and observational signatures

Several platforms now realize Kekulé-patterned lattices directly rather than as purely theoretical constructions.

Platform	Kekulé form	Reported consequence
Molecular graphene on Cu (Bergman, 2012)	Three Kekulé domains meeting at a Y-junction	Vortex, subgap zero mode, fractional charge $t_1$4 per spin
Graphene/TMD heterostructures (Zhang et al., 2022)	Sample-wide Kekulé-Y-1, Kekulé-Y-2, Kekulé-M	Long-range $t_1$5 charge-density wave
Sn-intercalated graphene on SiC (Ngo et al., 23 Dec 2025)	Kekulé-O	Gap $t_1$6
Graphene on CdS (Betancur-Ocampo et al., 2024)	Kekulé-O and Kekulé-Y/QBCP	Gap or quadratic band crossing, SOC $t_1$7

In molecular graphene, CO molecules on Cu define an artificial honeycomb lattice in which three inequivalent Kekulé domains can be positioned so as to meet at a Y-junction. Because the domains differ by phase shifts of $t_1$8, the total phase winding around the junction is $t_1$9, realizing a Kekulé vortex. The vortex core supports a zero-energy mode, and the bound charge is effectively fractionalized to $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$0 if spin is ignored (Bergman, 2012).

In graphene/TMD heterostructures, the TMD substrate supplies periodic scattering centers that couple $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$1 and $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$2, and STM/STS observes a sample-wide $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$3 charge-density wave with sharp Kekulé peaks in Fourier space. Three distinct textures—Kekulé-Y-1, Kekulé-Y-2, and Kekulé-M—are directly imaged at different energies, showing that multiple symmetry-allowed Kekulé bond textures can coexist within the same heterostructure family (Zhang et al., 2022).

In Sn-intercalated quasi-free-standing graphene on SiC(0001), LT-STM/STS, SPA-LEED, and DFT identify coexisting conventional and Kekulé-ordered graphene domains. The Kekulé phase is specifically Kekulé-O graphene with a $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$4 unit cell and a real-space period of about $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$5; STS shows a gap of about $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$6, present at both $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$7 K and $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$8 K (Ngo et al., 23 Dec 2025). On Cd-terminated CdS, DFT and an effective model identify two $$H= \varepsilon_0 \sum_i c_i^\dagger c_i + t_0\sum_{\langle i,j\rangle} c_i^\dagger c_j+ t_1\sum_{\langle i', j'\rangle} c_{i'}^\dagger c_{j'} ,$$9 graphene superlattices, a gapped Kekulé-O phase and a Kekulé-Y phase with a quadratic band crossing point, both with in-plane spin polarization induced by the substrate (Betancur-Ocampo et al., 2024).

Spectroscopy and transport provide additional signatures. In uniformly Y-shaped Kekulé-patterned graphene, optical conductivity at finite chemical potential shows a strong band nesting resonance that remains visible up to $\sqrt{3}\times\sqrt{3}$0, and this has been proposed as an optical fingerprint of the Y-shaped distortion (Mohammadi, 2022). In circularly gated Kekulé-Y regions, the split inner dispersion produces valley birefringence, separated caustic cusps, and a non-vanishing skew cross section associated with a valley Hall effect (Andrade et al., 2022).

6. Generalizations to phonons, magnons, solitons, and moiré systems

The Kekulé framework extends beyond electronic bond order. For phonons on a honeycomb lattice with a $\sqrt{3}\times\sqrt{3}$1 Kekulé modulation of force constants, the folded double-Dirac structure at $\sqrt{3}\times\sqrt{3}$2 supports a new phonon pseudospin with quantized Berry phases and pseudoangular momenta. The pseudospin Berry phase is

$\sqrt{3}\times\sqrt{3}$3

and the associated pseudospin Chern number is

$\sqrt{3}\times\sqrt{3}$4

leading to topologically protected pseudospin-polarized interface states, a phonon pseudospin Hall effect, and pseudospin-contrasting optical selection rules (Liu et al., 2017).

For magnons, a ferromagnetic honeycomb lattice with a Kekulé coupling texture and Dzyaloshinskii–Moriya interaction becomes a six-band bosonic problem with three critical lines and four topological phases. Chiral magnon edge states are determined by band Chern numbers, while Tamm-like edge states arise separately from intrinsic on-site interactions along the boundary sites (Pantaleon et al., 2018). In twisted bilayer easy-plane Néel antiferromagnets, the lattice elements themselves can be topological defects: meron and antimeron cores form hexagonal “meron hexads” with alternating short intracell distance $\sqrt{3}\times\sqrt{3}$5 and long intercell distance $\sqrt{3}\times\sqrt{3}$6, quantified by the Kekulé constant

$\sqrt{3}\times\sqrt{3}$7

so the soliton crystal realizes a genuine Kekulé-O type structure (Kim et al., 2024).

Moiré engineering supplies another generalization. In Kekulé moiré superlattices, a first layer is reconstructed into a $\sqrt{3}\times\sqrt{3}$8 supercell analogous to graphene’s Kekulé distortion so that it becomes nearly commensurate with a second hexagonal layer whose lattice constant is larger by roughly $\sqrt{3}\times\sqrt{3}$9. In MoTe$p_z$0/MnPSe$p_z$1, this folds remote valleys to the moiré $p_z$2 point, enables intervalley coupling through the substrate, produces valley pseudospin textures controlled by Néel-vector direction and external fields, and can yield a Chern insulator at one hole per moiré supercell (Ye et al., 2023). A plausible implication is that the Kekulé-patterned lattice has become a general symmetry-engineering principle for valley coupling, topology, and collective order across electronic, bosonic, magnetic, and moiré systems.