Three-Terminal ZGNR Devices

• Three-terminal ZGNR devices are nanoscale graphene systems with three electrodes enabling quantum control of charge and spin.
• They utilize zigzag edge states, tunable magnetic fields, and gate voltages to achieve spin filtering, half-metallicity, and quantum interference.
• These devices exhibit nonlinear I–V characteristics and gate-tunable thermoelectric performance, promising for advanced spintronics and energy applications.

A three-terminal zigzag graphene nanoribbon (ZGNR) device is a nanoscale quantum transport system in which a finite-width, finite-length ZGNR is connected to three electrodes. These include two leads at the ribbon’s left and right zigzag edges and a third electrode, typically at the top edge or as a voltage probe, enabling advanced control and measurement of electronic and spintronic transport phenomena. Three-terminal ZGNR devices form a versatile platform for studying and exploiting ballistic charge and spin currents, quantum interference, spin filtering, half-metallicity, and thermoelectric performance, owing to the interplay of graphene’s unique edge states, tunable magnetic fields, and gate voltages.

1. Device Geometry and Quantum Hamiltonian

A canonical three-terminal ZGNR device features a central ZGNR of width $W$ (number of zigzag chains) and length $L$ (number of unit cells), contacted by three semi-infinite leads: left (L, source), right (R, drain), and top (T, voltage probe or gate). The tight-binding Hamiltonian governing the electronic structure includes nearest-neighbor hopping, an out-of-plane Zeeman term for spin splitting, a local gate voltage, and coupling to the leads:

$$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$

• $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$, with $t \approx 2.7$ eV.
• $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$ introduces spin splitting via Zeeman field ($M_z$).
• $$H_{\text{gate}} = \sum_{i,\sigma} V_i c_{i\sigma}^\dagger c_{i\sigma}$$ for gate control.
• Each lead is represented as a semi-infinite system ($H_{\text{lead}}^{(\alpha)}$) coupled to device edges by $H_{\text{coup}}^{(\alpha)}$.

The lead self-energies $L$0, incorporating contacts into the device Green’s function, are computed via surface Green’s function techniques. Rashba spin-orbit coupling (SOC), modeled by $L$1, is relevant when substrate-induced SOC is significant (Ganguly et al., 2018).

2. Ballistic Transport Formalism

Ballistic transport through three-terminal ZGNRs is analyzed with nonequilibrium Green function (NEGF) methods within the Landauer–Büttiker framework. The central retarded Green’s function is

$L$2

Spin-resolved transmission between leads $L$3 is

$L$4

Charge and spin currents in lead $L$5 under bias are

$L$6

with $L$7 the Fermi-Dirac distribution at the local chemical potential. The total charge and spin currents are $L$8 and $L$9. The Fano factor $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$0 quantifies shot noise and quantum interference:

$$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$1

This formalism allows for detailed analysis of zero-bias conductance, spin polarization, and current-voltage characteristics under varying magnetic and electrostatic conditions, as well as shot noise and quantum coherence (Tamuli et al., 29 Jul 2025, Farghadan et al., 2015, Ganguly et al., 2018).

3. Spectral Features: Transmission, DOS, and Fabry–Pérot Physics

In the absence of magnetization, the spectral and transport properties of three-terminal ZGNRs are dominated by subband quantization and edge-localized states.

• The transmission $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$2 exhibits narrow peaks at the Fermi level from quasi-flat edge-state bands, and quantized steps as new transverse modes are activated by increased energy.
• The density of states (DOS) shows van Hove singularities at subband onsets and a pronounced peak at $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$3 from flat-band edge modes.
• Finite ribbon length $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$4 leads to Fabry–Pérot–like interference, with conductance modulations set by the resonance condition $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$5.
• As bias is increased, longitudinal mode activation produces step-like features and oscillations in $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$6 and $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$7, reflecting the quantized opening of new channels.

The interplay of quantum confinement and resonant edge states establishes a transport gap at low bias, crossing over to nonlinear behavior as bias passes $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$8 V, with current saturation at large $$H = H_{\text{hop}} + H_{Z} + H_{\text{gate}} + \sum_\alpha \left(H_{\text{lead}}^{(\alpha)} + H_{\text{coup}}^{(\alpha)}\right)$$9 due to all available modes conducting (Tamuli et al., 29 Jul 2025).

4. Spin Filtering, Half-Metallicity, and Symmetry Effects

Zeeman field $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$0 breaks spin degeneracy, shifting spin subbands by $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$1 and enabling strong spin filtering.

• For moderate $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$2 (e.g., 0.5 eV), edge state peaks for spin-up and spin-down separate by $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$3.
• At critical $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$4 in narrow ribbons, one spin species’ channel becomes fully gapped at $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$5, while the opposite spin remains gapless—realizing half-metallicity.
• Spin-resolved conductance $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$6 shows widely separated resonances at large $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$7; total conductance near the Fermi energy is then carried by a single spin channel, with spin polarization $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$8.
• For wide ribbons, spin filtering persists but with reduced selectivity due to bulk mode contributions.

With SOC (e.g., Rashba), spin-polarized currents can be generated even in the absence of a Zeeman field. The position and symmetry of the output leads strongly affect spin-polarized transmission, such that in mirror-symmetric configurations the $$H_{\text{hop}} = -t \sum_{\langle i,j \rangle, \sigma} (c_{i\sigma}^\dagger c_{j\sigma} + \text{H.c.})$$9 and $t \approx 2.7$0 polarizations in the two drains are equal in magnitude and opposite in sign, whereas $t \approx 2.7$1 polarization is identical. Asymmetric placement of the terminals enables further tunability (Tamuli et al., 29 Jul 2025, Ganguly et al., 2018).

5. Nonlinear Transport, Thermoelectricity, and Noise

The I–V characteristics of three-terminal ZGNR devices reveal several distinct regimes:

• At $t \approx 2.7$2, current is suppressed up to $t \approx 2.7$3 V, then rises nonlinearly as additional modes contribute.
• Introduction of $t \approx 2.7$4 enhances current by activating spin-polarized bands, especially in narrow ribbons, producing non-monotonic dependence on $t \approx 2.7$5 due to competition between spin-split states.
• Spin current $t \approx 2.7$6 is negligible at small $t \approx 2.7$7 but grows with bias and saturates at large $t \approx 2.7$8 in narrow ribbons.
• Oscillations in $t \approx 2.7$9 at low bias reflect Fabry–Pérot resonances.

The Fano factor, derived from shot noise calculations, exhibits oscillatory behavior near subband edges and is suppressed ($$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$0) within the half-metallic window where transmission is nearly perfect ($$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$1).

Hybrid ZGNR-molecular devices with a third ZGNR top gate enable gate-voltage manipulation of thermoelectric coefficients, including the Seebeck coefficient $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$2 and the thermoelectric figure of merit $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$3. The three-terminal configuration enables tuning of $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$4 via gate voltage, achieving $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$5 at 300 K (with greater enhancement at low $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$6). Phonon transport is suppressed due to the molecule's disruption of graphene’s vibrational modes, further optimizing $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$7 (Saha et al., 2011).

6. Gate, Geometry, and Magnetic Field Tunability

Key device parameters—ribbon width $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$8, length $$H_Z = M_z \sum_i (c_{i\uparrow}^\dagger c_{i\uparrow} - c_{i\downarrow}^\dagger c_{i\downarrow})$$9, gate voltage $M_z$0, and Zeeman field $M_z$1—provide orthogonal knobs for tuning electronic and spintronic performance:

• $M_z$2 and $M_z$3 set subband spacing and edge state overlap, enabling design of energy windows for efficient spin filtering or logic.
• $M_z$4 shifts the chemical potential and can enhance spin-current by positioning $M_z$5 in regions of maximal spin asymmetry.
• $M_z$6 independently controls the onset and width of the half-metallic window, with moderate fields ($M_z$7–1.0 eV) accessible via proximity coupling to a ferromagnet.
• The third terminal enables nonlocal spin and charge readout, separation, and routing of spin-polarized currents, and may function as a top gate for thermoelectric or logic applications.

Temperature stability is strong: spin polarization remains $M_z$8 at $M_z$9 V up to $$H_{\text{gate}} = \sum_{i,\sigma} V_i c_{i\sigma}^\dagger c_{i\sigma}$$0 K; spin current remains above 50% of its value at 300 K, demonstrating robust room-temperature operation (Tamuli et al., 29 Jul 2025).

7. Applications and Outlook

Three-terminal ZGNR devices provide a highly tunable route to spintronic and quantum transport functionalities:

• Gate- and magnetic-field–controllable spin-filters, spin valves, and spin logic elements with quantized edge-state transport and half-metallic operation (Tamuli et al., 29 Jul 2025, Zeng et al., 2010).
• Nonlocal spin detection and manipulation, as well as spin-charge current separation in multiterminal settings (Ganguly et al., 2018).
• Gate-tunable nanoscale thermoelectric devices with suppressed phonon conduction and optimized $$H_{\text{gate}} = \sum_{i,\sigma} V_i c_{i\sigma}^\dagger c_{i\sigma}$$1 (Saha et al., 2011).
• All-electrical generation, control, and amplification of fully spin-polarized currents without the need for ferromagnetic contacts, harnessing edge-state magnetism and Hubbard correlations in engineered nanostructures (Farghadan et al., 2015).

Design rules emphasize the use of well-defined zigzag edges (to maximize edge magnetism), optimal lead and gate architectures (to select target spin-polarization or logic), and platform scalability (via chemical synthesis or advanced lithography). Room-temperature stability and scalability strongly support the promise of three-terminal ZGNR architectures for practical graphene-based spintronics and energy conversion devices (Tamuli et al., 29 Jul 2025).