Zigzag Graphene Nanoribbons (zGNRs)

• zGNRs are quasi-one-dimensional graphene strips with zigzag edges that exhibit robust spin-polarized states and tunable narrow band gaps.
• Edge symmetry distinguishes even-N (symmetric) and odd-N (asymmetric) zGNRs, affecting quantum transport and conduction gap formation under bias.
• Chemical functionalization and doping enable band gap engineering and emergence of quantum phases, paving the way for advanced spintronic devices.

A zigzag graphene nanoribbon (zGNR) is a quasi-one-dimensional strip of graphene obtained by cutting the two-dimensional honeycomb lattice along the so-called zigzag direction, yielding edges characterized by the alternation of carbon atoms from the two graphene sublattices. zGNRs are distinguished by their unique edge-localized electronic and magnetic properties, including robust spin-polarized states, narrow band gaps whose size and nature are tunable by structural symmetry, chemical functionalization, and electrostatic fields, as well as ballistic and spin-selective quantum transport phenomena. These features place zGNRs at the forefront of carbon-based nanoelectronics and spintronics.

1. Atomic Structure, Edge Symmetry, and Electronic Topology

A zGNR consists of a finite width of N zigzag carbon chains, with the boundaries formed by zigzag edge carbon atoms. Each edge carbon terminates either with hydrogen (passivated), a functional group, or remains covalently bonded to a substrate or adjacent lattice. The lattice is bipartite, comprising A and B sublattices along the edge.

The distinction between "symmetric" and "asymmetric" zGNRs is set by the presence or absence of a mirror plane through the ribbon center. Even-N zGNRs are symmetric (have a central mirror plane), while odd-N zGNRs are asymmetric. This symmetry categorization profoundly impacts the electronic structure and transport characteristics under finite bias, since the symmetry dictates the allowed mixing of transverse subbands and the opening or closing of conduction gaps (Li et al., 2010).

The electronic structure of pristine, hydrogen-terminated zGNRs is dominated by two nearly flat edge bands at the Fermi level. These originate from non-bonding π orbitals localized on the outermost carbon rows, with wavefunctions that decay exponentially into the ribbon, leading to a high local density of states at the edges (Blackwell et al., 2021, Ruffieux et al., 2015). The presence of interactions triggers a spontaneous spin-polarization instability—ferromagnetic ordering develops along each edge (all spins aligned), antiferromagnetically coupled between opposing edges, following mechanisms elucidated by mean-field Hubbard models and first-principles density functional theory (DFT).

2. Symmetry, Quantum Transport, and Band Gap Engineering

Symmetry governs quantum transport in zGNRs under an applied source-drain bias. In symmetric (even-N) ribbons, the π and π* subbands have opposite parity and cannot couple in the presence of a bias, opening a parity-protected conductance gap near the Fermi level. Current remains essentially zero until the bias exceeds a critical threshold $V_c \sim \Delta/e$, with $\Delta \sim 1/N$ the subband energy separation (Li et al., 2010). Odd-N (asymmetric) zGNRs lack this restriction, allowing subband mixing and supporting linear, ohmic conductance above the gap.

The transport formalism is built on the non-equilibrium Green’s function (NEGF) approach combined with DFT. The spin-dependent current is given by:

$$I_\sigma(V) = \frac{e}{h} \int_{\mu_R}^{\mu_L} T_\sigma(E,V) dE\,,$$

where $T_\sigma$ is the spin-resolved transmission, and $\Gamma_{L,R}$ are broadening matrices derived from the self-energies of the leads (Wu et al., 2012). This approach also incorporates the screening response, spatial variation of bond and local edge currents, and effects of mean-field electron-electron interactions (Cheraghchi, 2012).

Band gap engineering in zGNRs can be realized through edge patterning, chemical functionalization, or introduction of edge defects. When the number of type-A and type-B edge atoms is equal (no sublattice imbalance), all zero-energy edge states are quenched, and a significant bandgap opens—this general boundary-condition is key for designing semiconducting zGNRs as field-effect transistor (FET) channels (Zhang et al., 2011). The magnitude and character of the gap can be further tuned by edge passivation (e.g., sp³-hybridized O-functionalization), with nonmonotonic dependence on defect ratio and arrangement (Wadhwa et al., 2018).

3. Magnetism, π-Magnetism, and Edge State Manipulation

The flat band edge states at $E_F$ create a peak in the local density of states and facilitate π-magnetism, where ferromagnetic order develops along each zigzag edge, antiferromagnetically coupled across the ribbon width. DFT shows local moments of $\mu_C \sim 0.3$–$0.4\,\mu_B$ per edge carbon, with the exchange gap $\Delta_E$ between conduction and valence edge bands scaling as $\Delta_E \propto 1/W$ (ribbon width) (Blackwell et al., 2021, Ruffieux et al., 2015).

Magnetic interaction parameters extracted from total-energy DFT calculations show antiferromagnetic coupling energies in the range of $10$–$20$ meV per unit cell (Pizzochero et al., 2022). The presence of edge defects, chemical functionalization (e.g., hydrogen, OH, COOH), and substrate effects can disrupt or engineer the magnetic order, leading to the controlled emergence of paramagnetic centers (spin-½) and modifying spin transport (Pizzochero et al., 2022, Li et al., 2010).

In Janus GNRs, a designed sublattice imbalance $\Delta N=1$ per unit cell (realized by decorating one edge with topological benzene-derived defects) generates a robust, single ferromagnetic edge state predicted by Lieb's theorem. This facilitates the realization of a 1D Heisenberg spin-½ chain delocalized on a single edge, with the ground state verified experimentally by STS and DFT (Song et al., 2024).

4. Chemical Functionalization, Doping, and Substrate Effects

Substitutional doping and chemical termination of edges modulate the local potential and magnetic environment, enabling further tuning of spin and charge transport. Doping with group III elements (e.g., Be, B, Al) on one edge breaks the mirror symmetry, induces spin-dependent localization, and can result in strongly spin-polarized transport properties. In symmetric even-N zGNRs, Be doping leads to negative differential resistance (NDR) behavior only for the down-spin channel at $V \sim 1.0$ V, while up-spin current rises monotonically—a mechanism promising for spin-based oscillator devices (Wu et al., 2012, Jiang et al., 2014).

Asymmetric edge terminations, such as H–OH or H–COOH, lift the degeneracy of the edge bands, reducing the critical transverse electric field required to achieve half-metallicity (metallic for one spin, semiconducting for the other). This can result in nearly zero-bias, spin-polarized transport, and Dirac-cone-like dispersions for one spin channel at the Fermi level (Li et al., 2010).

The placement of zGNRs on substrates (SiO₂, h-BN) further alters their properties through hybridization, surface-induced distortion, or in-plane edge bonding. Embedding zGNRs in h-BN can stabilize edges via covalent C–B/C–N linkages, building in a transverse electric field that further controls band structure, edge state splitting, and offers robust, high-temperature magnetic conduction channels (Jiang et al., 17 Nov 2025, Kim et al., 2015).

5. Ballistic and Thermoelectric Transport Regimes

The quantum transport properties of zGNRs—particularly under finite bias—are governed by the subband structure, symmetry, and scattering phenomena. Ultra-narrow zGNRs show selection-rule-protected NDR regions, quantized transmission plateaus, and long coherence lengths of edge states. The on/off ratio of current can increase exponentially with ribbon length, reaching values above $10^5$ for device-scale ribbons (Cheraghchi, 2012, Bhalla et al., 2017).

Thermoelectric performance can be optimized by exploiting miniband formation in periodically patterned or defected zGNRs (graphene quantum dot arrays), contact geometry (line vs. surface), and valley degeneracy engineered via staggered sublattice potentials or transverse fields. Valley degeneracy, induced by substrate or field-driven band inversion, enhances conductance and power factor, yielding figures of merit $ZT > 3$ at room temperature under suitable contact configurations (Kuo, 2022, Kuo, 2024).

6. Quantum Phases: Ferromagnetism, Half-Metallicity, and Topological Superconductivity

Edge engineering enables the creation of new quantum phases in zGNRs. Terminating one edge with a π-conjugated biradical (e.g., trimethylenemethane) or fabricating Janus GNRs with an asymmetric edge defect array can induce a robust ferromagnetic semiconductor phase, characterized by a large spin-split gap ($\sim 1$ eV) and 100% spin polarization in a window around the Fermi level—conditions sufficient for room temperature spin transistors and quantum information applications (Hou et al., 2010, Song et al., 2024).

Under transverse electric fields or h-BN embedding, zGNRs enter a half-metallic state, with one spin species conducting and the other gapped. Coupling a half-metallic zGNR to an Ising superconductor induces odd-parity triplet superconductivity and a gate-tunable topologically nontrivial phase supporting Majorana zero modes, as described by a Kitaev-chain-like effective Hamiltonian (Ma et al., 17 Jun 2025).

7. Experimental Realization, Characterization, and Outlook

Bottom-up atomic precision synthesis via on-surface polymerization and cyclodehydrogenation has enabled the fabrication of zGNRs with atomically perfect edges of specified width and topology. Scanning probe microscopies (STM, nc-AFM) and tunneling spectroscopy (STS) reveal the spatial and energetic structure of edge states and confirm their magnetic nature (Ruffieux et al., 2015, Blackwell et al., 2021, Song et al., 2024). Magnetic order at the edges has now been directly detected in zGNRs embedded in h-BN using NV-center magnetometry, with room-temperature stability and ballistic transport established (Jiang et al., 17 Nov 2025).

Outlook for zGNRs is dictated by their ability to realize robust quantum spin chains, reconfigurable spintronic elements, high-performance FETs, and potentially topologically protected superconducting states with electrically controlled Majorana modes. The control and understanding of edge magnetism, chemical reactivity, quantum coherence, and device integration remain areas of active research.