Maple Leaf Lattice (MLL)

• Maple Leaf Lattice (MLL) is a quasi-two-dimensional Archimedean tiling formed by systematically depleting one-seventh of the triangular lattice sites, resulting in a unique structure of triangles and hexagons with coordination number 5.
• The MLL exhibits strong geometric frustration that fosters exotic magnetic phases, including canted 120° order, vector-chirality, and valence-bond states, as well as potential quantum spin-liquid behavior.
• Advanced numerical methods such as exact diagonalization, tensor networks, and coupled-cluster studies reveal detailed thermodynamic signatures and topological excitations, highlighting its relevance in studying complex antiferromagnetic systems.

The Maple Leaf Lattice (MLL) is a quasi-two-dimensional Archimedean tiling of the plane, defined by a regular one-seventh site depletion of the triangular lattice. Each site on the MLL is part of four triangles and one hexagon, conferring it a coordination number of $z=5$ and a unit cell containing six inequivalent sites. Strong geometric frustration arises from its network of edge-sharing triangles and hexagons, resulting in a rich landscape of classical and quantum phenomena in antiferromagnetic and related Hamiltonians.

1. Lattice Geometry, Symmetry, and Modeling Frameworks

The MLL is most commonly constructed by systematically removing one out of every seven lattice points from the triangular lattice, with each vacancy located at the center of a six-site hexagonal ring. The resulting structure preserves a triangular Bravais lattice (primitive vectors $a_1=(1,0)$, $a_2=(\frac{1}{2},\frac{\sqrt{3}}{2})$) with a six-site basis. Three symmetry-inequivalent types of nearest-neighbor bonds arise:

• "Triangle" bonds: forming corner-sharing triangles,
• "Hexagon" bonds: edges of the six-site rings,
• "Dimer" bonds: bridging remaining nearest neighbor pairs.

No reflection symmetry is present, but sixfold rotation, translation, and some inversion operations (depending on the three-dimensional embedding and specific compounds) are admitted (Aguilar-Maldonado et al., 2024, Schäfer et al., 26 Nov 2025). The generic spin Hamiltonian takes the form

$$H = \sum_{\langle ij\rangle_{d}} J_d\,\mathbf{S}_i\cdot\mathbf{S}_j + \sum_{\langle ij\rangle_{t}} J_t\,\mathbf{S}_i\cdot\mathbf{S}_j + \sum_{\langle ij\rangle_{h}} J_h\,\mathbf{S}_i\cdot\mathbf{S}_j + D\sum_i (S_i^z)^2$$

with antiferromagnetic $J_d,J_t,J_h>0$ and single-ion anisotropy $D$ (when relevant).

2. Classical and Quantum Magnetism: Phase Structure

Classically, the competition between the three types of AFM exchange causes nontrivial ground states:

• Canted 120° order: In the uniform case ($J_{d}=J_{t}=J_{h}$), spins on each triangle form coplanar 120° structures, but the lack of complete equivalence between triangles and hexagons leads to a nontrivial canting between neighboring triangles (Farnell et al., 2018, Nyckees et al., 23 Dec 2025).
• Vector-chirality phases: Chiral degrees of freedom emerge, and distinct ordering patterns arise:
  • Staggered vector chirality (SVC): adjacent triangles have opposing vector chiralities,
  • Uniform vector chirality (UVC): all triangles share the same chirality, stabilized only by further-neighbor AFM interactions (Ghosh, 16 Apr 2025, Aguilar-Maldonado et al., 2024).
• Valence bond order and dimerization: When one bond type dominates (e.g., $J_d \gg J_t, J_h$), the ground state becomes a product of singlets on dimer bonds, yielding an exact quantum dimer phase (Nyckees et al., 23 Dec 2025, Ghosh et al., 2023). This phase appears generically for large bond anisotropy and is separated by a first-order quantum phase transition from the canted 120°-ordered state.

Quantum treatments—via coupled-cluster, exact diagonalization, iPEPS, functional RG, and tensor networks—reveal an intricate phase diagram:

• At the isotropic point, quantum fluctuations substantially reduce magnetic order, leading to possible quantum spin-liquid (QSL) ground states or weakly ordered 120° structures depending on methodology and cluster size (Gresista et al., 2023, Sonnenschein et al., 2024, Schmoll et al., 2024, Ebert et al., 8 Jan 2026).
• Exact dimer and valence-bond crystal phases compete with chiral and collinear Néel orders in the vicinity of maximal frustration (Ebert et al., 8 Jan 2026, Nyckees et al., 23 Dec 2025).
• Gapped $\mathbb{Z}_2$ spin liquid and $U(1)$ Fermi surface spin liquid states are realized in slivers of parameter space, identified through projective symmetry group analysis, Gutzwiller-projected variational Monte Carlo, and direct comparison with pf-FRG-derived structure factors (Sonnenschein et al., 2024, Gresista et al., 2023, Ebert et al., 8 Jan 2026).

3. Thermodynamic and Dynamical Characterization

Thermodynamic investigations using high-temperature expansions, numerical linked-cluster expansions (NLCE), and exact diagonalization have established:

• The ground-state energy per site for the $S{=}1/2$ isotropic Heisenberg model converges to $E_0/N \approx -0.531$ (units of $J$), consistent across multiple advanced methods (Schäfer et al., 26 Nov 2025, Hutak, 16 May 2025, Nyckees et al., 23 Dec 2025).
• The specific heat $C(T)$ displays a double-peak structure in the isotropic case ($T_1 \approx 0.48 J$, $T_2 \approx 0.13 J$). This is interpreted as the release of entropy by gapped/Schottky-like triplon excitations and low-lying magnon modes, analogous to features found in the triangular lattice but strongly influenced by MLL's lack of reflection symmetry (Schäfer et al., 26 Nov 2025, Hutak, 16 May 2025).
• The uniform susceptibility $\chi(T)$ peaks at $T \approx 0.49 J$, and the zero-temperature value is $\chi_0 \approx 0.05-0.06$ for $S=1/2$ (Hutak, 16 May 2025).
• The equal-time spin structure factor $S(\mathbf{q})$ reveals broad maxima at $K$ points without sharp Bragg peaks for putative QSL or dimer phases, and evolves a characteristic "chiral twist" as temperature is lowered, directly tracking the absence of mirror symmetry (Schäfer et al., 26 Nov 2025, Schmoll et al., 2024).

Magnetization process studies reveal fractional magnetization plateaux due to the intricate interplay of localized singlet formation and correlated triplon kinetics; for example, in dimer-dominated limits, plateaux at $m = 1/6, 2/9, 2/7, 1/3$ emerge through classical and quantum density-wave crystallization mechanisms (Ghosh et al., 2023, Schmoll et al., 2024).

4. Realizations in Materials: Rare Earth, Transition Metal, and Mineral Platforms

Several families of materials closely approximate an ideal or distorted MLL:

• Ho$_3$ScO$_6$: Provides a rare-earth MLL realizing uniform vector-chirality ordering below $T_N = 4.1$ K, with experimental evidence for a 120° coplanar spin structure and positive vector chirality stabilized by a hierarchy $J_d \gg J_t > J_h$ (Aguilar-Maldonado et al., 2024, Ghosh, 16 Apr 2025).
• MgMn$_3$O$_7\cdot 3$H$_2$O: A $S{=}3/2$ realization displaying successive magnetic transitions at 15 K and 5 K, a complex chiral order landscape, and field-induced plateau phenomena up to 60 T, sensitive to Dzyaloshinskii–Moriya and further-neighbor interactions (Haraguchi et al., 2018).
• Bluebellite and Spangolite: Layered copper minerals with $S{=}1/2$ whose Cu$^{2+}$ ions realize the MLL or its distortions, though in practice ground states are often dimerized or cluster singlets as a consequence of bond anisotropy and weak ferromagnetic exchanges (Schmoll et al., 2024, Haraguchi et al., 2021).
• Na$_2$Mn$_3$O$_7$: Chemical/electrochemical control enables switching between 1D antiferromagnetic chains and glassy, frustrated magnetic states, demonstrating the tunability of competing orders by structural and disorder effects (Chao et al., 2023, Saha et al., 2023).

A subset of these compounds hosts exact dimer phases or valence-bond order and exhibits thermodynamic signatures that mimic one-dimensional quantum magnets despite the underlying two-dimensional topology (Haraguchi et al., 2021, Ghosh et al., 2023).

5. Topological Excitations, Thermal/Spin Transport, and Exotic Orders

The nontrivial geometry of the MLL supports a variety of unconventional excitations and transport:

• Triplon bands in the quantum dimer regime become topologically nontrivial, mapping onto Kagome magnon topological models with nonzero Chern numbers, $Z_2$ indices, and symmetry-protected edge modes (Esaki et al., 14 May 2025). These yield quantized spin Nernst and thermal Hall effects, with topological transitions and sign reversals as a function of Dzyaloshinskii–Moriya vectors and field.
• Altermagnetism: Non-collinear magnetic orders on the MLL—classified by selective breaking of parity ($\mathcal{P}$) and time reversal ($\mathcal{T}$)—give rise to prominent momentum-dependent, nonrelativistic magnon spin splitting. Both strong-coupling (Heisenberg) and weak-coupling (Hubbard) models generate a taxonomy of canted 120°, $q=0$ altermagnetic, and other exotic orders, whose experimental fingerprints manifest in magnon dispersion and even-parity spin gaps (Ghosh et al., 23 Jan 2026).

The rich competition of valence-bond, chiral, plaquette, and (gapped or gapless) quantum spin-liquid orders at and near the isotropic point produces a fertile setting for phenomena such as deconfined quantum criticality, multicriticality, and possible realization of Dirac or $Z_2$ quantum spin liquids (Gresista et al., 2023, Ebert et al., 8 Jan 2026, Sonnenschein et al., 2024).

6. Projective Symmetry Groups, Quantum Spin Liquids, and Theoretical Outlook

Analyses leveraging Projective Symmetry Group (PSG) constructions identify 17 distinct $U(1)$ and 12 distinct $\mathbb{Z}_2$ spin-liquid phases accessible to the MLL for various choices of nearest- and further-neighbor mean-field amplitudes (Sonnenschein et al., 2024). Among these, $U(1)$ Fermi surface and gapped $\mathbb{Z}_2$ spin-liquid states most closely reproduce the signature broad triangular structure factor maxima and near-absence of pinch points observed in functional RG and tensor network studies (Gresista et al., 2023, Sonnenschein et al., 2024, Schmoll et al., 2024).

Theoretical challenges remain in precisely characterizing the nature and stability of these QSL regimes—particularly their low-energy gauge structure ($U(1)$ vs $Z_2$), and potential for topologically ordered excitations—and their relationship to proximate ordered or valence-bond phases.

7. Comparison with Other Frustrated Lattices and Open Directions

The MLL stands at the unique intersection between the maximally frustrated kagome lattice ($z=4$, infinite ground-state degeneracy) and the less frustrated triangular lattice ($z=6$, robust 120° order). Quantum and classical studies consistently find the MLL is characterized by

• Intermediate levels of frustration and quantum reduction of order parameters (Farnell et al., 2018),
• Emergent one-dimensional behaviors in some mineral realizations,
• Absence or strong suppression of long-range order at the lowest accessible temperatures in certain compounds,
• Possibility of tuning between competing chiral, valence-bond, and spin-liquid orders by manipulating exchange anisotropy, further-neighbor coupling, lattice distortions, dimensionality (interlayer couplings), or disorder (Aguilar-Maldonado et al., 2024, Ebert et al., 8 Jan 2026, Chao et al., 2023).

Future experimental and theoretical directions include synthesis and inelastic spectroscopy on candidate rare-earth and transition-metal oxides, detailed studies of thermal/spin transport and thermal Hall effects, tensor network and large-scale DMRG computation for larger system sizes, and direct investigation of quantum criticality and topological orders in this frustrated, highly tunable quantum magnet platform.

---

Key References:

(Aguilar-Maldonado et al., 2024, Schäfer et al., 26 Nov 2025, Hutak, 16 May 2025, Farnell et al., 2018, Nyckees et al., 23 Dec 2025, Gresista et al., 2023, Ebert et al., 8 Jan 2026, Schmoll et al., 2024, Sonnenschein et al., 2024, Ghosh, 16 Apr 2025, Haraguchi et al., 2018, Chao et al., 2023, Saha et al., 2023, Schmoll et al., 2024, Haraguchi et al., 2021, Ghosh et al., 2023, Esaki et al., 14 May 2025, Ghosh et al., 23 Jan 2026)