Triangulene: Graphene Nanoflake & π-Magnetism

• Triangulene is a triangular graphene nanoflake with zigzag edges and intrinsic open-shell character that underlies its unique π-magnetic properties.
• Studies reveal that its zero modes and sublattice polarization, explained by tight-binding and Hubbard models, define its spin configuration and magnetic behavior.
• On-surface synthesis has enabled the integration of triangulene into dimers, chains, and 2D networks, paving the way for designer quantum spin devices.

Triangulene is an equilateral triangular graphene nanoflake with zigzag edges and intrinsic open-shell character. In the $[n]$ notation, each edge contains $n$ fused benzene units; isolated $[n]$ triangulenes host $n-1$ zero modes and, at half filling, a ground-state spin $S=(n-1)/2$. The archetypal triangulene is the smallest triangular zigzag-edged graphene fragment with intrinsic open-shell character, composed of six fused benzene rings, and it has become a canonical system for carbon $\pi$-magnetism, on-surface synthesis, molecular Heisenberg physics, and triangulene-derived two-dimensional quantum materials (Ortiz et al., 2022, Mishra et al., 2020).

1. Molecular topology and nomenclature

Triangulene is variously described as a triangular zigzag-edged graphene fragment, a triangular graphene quantum dot, and a non-Kekulé triangular polycyclic aromatic hydrocarbon. These descriptions emphasize the same structural fact: the honeycomb $\pi$-network is cut into an equilateral triangle whose three sides are zigzag edges, so the two graphene sublattices are populated unequally. In the modern size notation, $[2]$, $[3]$, and $[4]$ triangulenes carry $n$0, $n$1, and $n$2, respectively, while the commonly discussed “triangulene” corresponds to $n$3triangulene, with two unpaired $n$4-electrons and intrinsic $n$5 (Catarina et al., 2023, Hieulle et al., 2021).

The non-Kekulé designation is central. For triangulene, no Kekulé valence structure can pair all $n$6 electrons into ordinary $n$7 bonds without leaving unpaired electrons, even though the neutral molecule has an even number of atoms and electrons. This impossibility is the origin of its open-shell character and of its status as the smallest triplet-ground-state polybenzoid or Clar’s hydrocarbon in the older literature (Sharma et al., 2017, Mishra et al., 2020).

A persistent misconception is that triangulene magnetism is merely a generic “edge effect.” The literature instead treats it as a consequence of a specific combination of zigzag boundary geometry, bipartite-lattice sublattice imbalance, and zero-mode degeneracy. The edge localization is real, but it is topology-driven rather than incidental (Rodrigues et al., 2024).

2. Zero-energy shell, sublattice polarization, and spin rules

In single-particle theory, isolated $n$8 triangulenes host $n$9 non-bonding zero modes. At charge neutrality those modes are half filled, and Coulomb repulsion favors maximal spin, giving

$[n]$0

in agreement with Lieb’s theorem and Ovchinnikov’s rule for bipartite lattices (Ortiz et al., 2022, Catarina et al., 2023). In the nearest-neighbor tight-binding description,

$[n]$1

the zero-energy shell is defined by

$[n]$2

For triangular graphene quantum dots, projecting the zero-energy equation onto a site yields the local constraint

$[n]$3

and the resulting shell states are sublattice-polarized: in the ideal model they reside only on the majority sublattice (Rodrigues et al., 2024).

This sublattice polarization is not a minor technicality. It governs where edge-state spectral weight appears, why minority-sublattice perturbations can be ineffective, and how degeneracies are lifted under local chemical modification. For $[n]$4triangulene, for example, the zero-energy shell contains four states, and the analytical construction shows explicitly that the minority-sublattice solution is trivial once boundary conditions are imposed (Rodrigues et al., 2024).

Correlation effects then convert the zero-mode shell into magnetism. For triangulene-like flakes on Au(111), the low-energy sector is routinely modeled with a $[n]$5-electron Hubbard Hamiltonian,

$[n]$6

whose mean-field solutions recover singly occupied zero modes, localized edge spin density, and exchange-stabilized high-spin states (Li et al., 2019). Collectively, these results establish that triangulene zero modes are flat-band-like, partially filled, and unusually susceptible to correlation-driven spin polarization.

3. On-surface synthesis and single-molecule spin spectroscopy

The decisive experimental advance in the field was the transition from predicted open-shell nanographenes to atomically precise triangulene-like molecules fabricated directly on metal surfaces. On Au(111), an extended triangulene (ETRI) with $[n]$7 and $[n]$8 was synthesized from brominated precursors deposited onto a gold surface held at $[n]$9. Low-temperature $n-1$0 spectroscopy revealed a narrow zero-bias resonance with full width at half maximum of about $n-1$1 meV, a Kondo temperature $n-1$2 K, magnetic-field splitting visible already at $n-1$3 T, and a fitted $n-1$4, consistent with an underscreened $n-1$5 Kondo state rather than a conventional $n-1$6 impurity. The same work showed atom-by-atom switching between nonmagnetic, $n-1$7, and $n-1$8 states by sequential hydrogen passivation and depassivation of radical sites (Li et al., 2019).

Nitrogen substitution introduced a different control axis: interfacial charge transfer. Aza-triangulene, synthesized on Au(111) and Ag(111) by reduction of ketone-substituted precursors using atomic hydrogen followed by thermal annealing and tip manipulation, behaves as a positively charged open-shell triplet on Au(111) but as a negatively charged closed-shell species on Ag(111). On Au(111), Kondo resonances and orbital maps match positively charged aza-triangulene; on Ag(111), the absence of Kondo features and the observed HOMO/LUMO pattern match the anion (Wang et al., 2021).

Larger heteroatom-doped triangulenes retain the same basic interplay between topology and charging, but at higher total spin. For aza-$n-1$9-triangulene on Au(111), bond-resolved imaging uncovered edge-distributed radical states, and weak zero-bias resonances with field-induced splitting up to $S=(n-1)/2$0 T were interpreted as partial Kondo screening of a cationic $S=(n-1)/2$1 state rather than a neutral $S=(n-1)/2$2 state. Orbital maps identified the relevant frontier channel and supported the conclusion that positive charging on Au(111) restores the four singly occupied frontier states characteristic of all-carbon $S=(n-1)/2$3triangulene (Vilas-Varela et al., 2023).

These studies collectively overturned the earlier concern that metallic substrates would inevitably quench triangulene magnetism. They instead showed that substrate coupling can preserve, partially screen, or even chemically select the realized spin state, depending on molecular topology and charge transfer (Li et al., 2019, Wang et al., 2021).

4. Dimers, chains, rings, nanostars, and trimers

Once individual triangulenes were established as magnetic building blocks, the next question was whether they could be coupled into designer spin Hamiltonians. Covalently bonded triangulene dimers provided the first direct answer. Two architectures were synthesized on Au(111): a directly linked dimer and a dimer separated by a 1,4-phenylene spacer. Both retain open-shell triangulene units but realize open-shell singlet ground states, with singlet-triplet inelastic thresholds of $S=(n-1)/2$4 mV and $S=(n-1)/2$5 mV, respectively. The low-energy sector is captured by an effective Heisenberg dimer of two spin-1 objects,

$S=(n-1)/2$6

with the threshold directly measuring $S=(n-1)/2$7 (Mishra et al., 2020).

Larger assemblies revealed genuinely collective spectra. A triangulene nanostar made from six $S=(n-1)/2$8triangulene units was synthesized on Au(111) by converting a cyclic hexamer precursor through on-surface cyclodehydrogenation at $S=(n-1)/2$9. The resulting macrocycle carries 12 unpaired $\pi$0-electrons and behaves as a ring of six antiferromagnetically coupled $\pi$1 sites. Inelastic tunneling spectroscopy resolved three spin excitations at

$\pi$2

well reproduced by a six-site Heisenberg ring with $\pi$3 (Hieulle et al., 2021). At the other end of the size spectrum, TTAT, built from pristine and N-doped triangulene motifs, realizes a symmetric ferromagnetic Heisenberg-like spin trimer with three singly occupied natural orbitals, a quartet ground state $\pi$4, and an effective

$\pi$5

with $\pi$6 meV in the fitted model (Vegliante et al., 14 May 2025).

The most developed one-dimensional platform uses $\pi$7triangulene as a delocalized spin-$\pi$8 building block. Finite covalent chains on Au(111) implement the antiferromagnetic Heisenberg Hamiltonian

$\pi$9

with $\pi$0 meV extracted from the dimer singlet-triplet step. Even chains have singlet ground states, odd chains have doublet ground states, the excitation gap decreases approximately as

$\pi$1

and Fourier analysis of site-resolved $\pi$2 maps reveals an m-shaped dispersion interpreted as confined spinons in a one-dimensional quantum box (Yuan et al., 2024).

Periodic boundary conditions add another layer of control. Using pristine, unsubstituted $\pi$3triangulene units on Au(111), cyclic five- and six-membered rings were constructed by stepwise on-surface synthesis and STM tip-induced dehydrogenation. The six-membered ring is planar and realizes a nearly uniform antiferromagnetic $\pi$4 ring, whereas the five-membered ring is buckled, has nonuniform exchange couplings, and shows distortion-lifted degeneracy with asymmetric Kondo and spin-flip maps (Li et al., 12 Feb 2026). Collectively, these systems establish triangulene as a modular route from two-spin exchange to finite-size spinons, antiferromagnetic spin rings, and ferromagnetic molecular trimers.

5. Two-dimensional triangulene crystals and competing correlated phases

Triangulene also functions as a “superatom” for extended lattices. A unified theory of triangulene two-dimensional crystals considered honeycomb networks $\pi$5 whose unit cell contains two triangulenes. Projecting onto the zero-mode manifold yields a generalized honeycomb Hamiltonian with internal $\pi$6-symmetric pseudospin. Depending on triangulene size, the resulting low-energy bands realize graphene-like Dirac physics for $\pi$7, spin-1 Dirac electrons for $\pi$8, $\pi$9-orbital honeycomb physics for $[2]$0, and a gapped flat-band system for $[2]$1 (Ortiz et al., 2022).

Once interactions are included, broken-symmetry magnetic phases become central. For centrosymmetric crystals such as $[2]$2, $[2]$3, and $[2]$4, mean-field Hubbard theory and spin-unrestricted DFT predict antiferromagnetic order, while non-centrosymmetric systems such as $[2]$5 become ferrimagnetic. In $[2]$6, DFT gives

$[2]$7

and the gap increases from

$[2]$8

to

$[2]$9

with a moment per triangulene of about $[3]$0, close to the isolated $[3]$1 value $[3]$2. Exchange constants extracted from DFT, mean field, CAS, and RPA lie in the $[3]$3–$[3]$4 meV range for $[3]$5 and $[3]$6–$[3]$7 meV for $[3]$8 (Catarina et al., 2023).

A distinct extended architecture, a $[3]$9triangulene Kagome covalent organic framework on Au(111), was interpreted differently. In that work, the flat valence and conduction bands of the triangulene-based Kagome lattice were argued to support an excitonic insulator ground state: STS showed a semiconducting gap of about $[4]$0 eV, GW-BSE yielded a triplet exciton eigenenergy of $[4]$1 eV, and the inferred excitonic order parameter was about $[4]$2 eV (Delgado et al., 2023). The literature therefore contains both broken-symmetry magnetic and flat-band excitonic interpretations for narrow-band triangulene materials, depending on lattice geometry and theoretical treatment. What is common to both is that triangulene zero modes generate unusually narrow, often flat, bands for which interaction effects are not perturbative.

6. Heteroatom engineering, superatomic graphene, and quantum platforms

Heteroatom substitution turns triangulene into a local spectroscopic probe as well as a magnetic building block. For aza-$[4]$3triangulene, a single substitutional nitrogen impurity acts as a rank-one perturbation within the zero-mode shell: $[4]$4 Only one level splits from the degenerate shell, the remaining $[4]$5 states stay at zero energy, and the splitting tracks the local zero-mode weight at the impurity site. In the ideal model, placing the impurity on the minority sublattice has no effect because the zero-mode amplitudes vanish there (Rodrigues et al., 2024). Fusion topology adds another control parameter. In aza-triangulene dimers on Au(111), all major fused products are oxidized by roughly one electron, but only asymmetric dimers show Kondo fingerprints. The proposed mechanism is that asymmetric, non-benzenoid fusion localizes frontier orbitals, increases Coulomb repulsion, and promotes single occupancy, whereas the symmetric dimer remains more delocalized and nonmagnetic; on Ag(111), the corresponding products are closed-shell (Calupitan et al., 2023).

A more radical extension replaces each graphene lattice site by a triangulene-derived molecular “superatom.” Superatomic graphene built from phosphorus-doped triangulene was reported by bottom-up on-surface synthesis, with STM revealing metallic band structures and magnon excitation under varying magnetic fields and DFT attributing the spin-polarized bands to $[4]$6-like orbitals of the triangulene-derived units (Ding et al., 2024). A later, more detailed study of phosphorus-doped triangulene on Ag(111) resolved a honeycomb superlattice with Dirac and flat bands; DFT predicted intrinsic half-metallic ferromagnetism in pristine PT-based superatomic graphene, with the spin-down channel metallic and the spin-up channel gapped by about $[4]$7 eV. Oxygen functionalization changed the isolated molecular spin from $[4]$8 to $[4]$9, shifted the flat bands away from $n$00, opened an idealized gap of about $n$01 eV, and converted the network into a semiconducting antiferromagnet (Zhai et al., 18 Jun 2026).

Triangulene has also entered the quantum-technology literature. Advanced first-principles calculations on triangulene, its aza-cation derivative, and a substituted triangulene crystal found triplet ground states separated from the first singlet by about $n$02 eV, spin-lattice relaxation as long as $n$03 ms at 300 K for the aza-cation prototype, and $n$04 ms at 10 K for a deuterated triangulene crystal in a nuclear-spin-free environment, together with highly spin-selective intersystem crossing relevant to optical readout and initialization (Sarkar et al., 9 Jul 2026). In parallel, a numerical renormalization group proposal argued that triangulene spin chains on superconducting substrates could host a superconducting spin-singlet qubit formed by two isolated low-lying singlet states undergoing an avoided crossing, with a mesoscopic triple-quantum-dot device proposed as an emulator and control architecture (Huang et al., 31 Mar 2025).

Across these developments, triangulene has shifted from a historically elusive non-Kekulé hydrocarbon to a chemically programmable platform spanning single-molecule $n$05-magnetism, local zero-mode engineering, Heisenberg spin networks, flat-band quantum materials, half-metallic superatomic graphene, and molecular quantum-device concepts.