Peripheral Nitrogen Topology as a Defect-Chemical Switch for Electronic and Magnetic States in Graphene: A First-Principles Study of Pyridinic, Pyridazinic, Pyrrolic, and Pyrazolic Configurations

Technical Summary: Peripheral Nitrogen Topology as a Defect-Chemical Switch in Graphene

Problem Statement
While defect and heteroatom engineering are established methods for tuning the electronic and magnetic properties of graphene, the specific role of peripheral nitrogen topologies surrounding graphene voids remains insufficiently understood. Previous computational studies have largely focused on graphitic, pyridinic, and pyrrolic nitrogen structures. However, there is a lack of systematic comparison regarding more specific peripheral nitrogen configurations that mimic organic heterocyclic structures (specifically pyridazinic and pyrazolic motifs) around a central vacancy. The precise link between these local heterocyclic-like nitrogen topologies and the resulting structural reconstruction, charge redistribution, band-structure modulation, and magnetic moment formation requires clarification to enable the design of metal-free spintronic and semiconducting carbon devices.

Methodology
The study employs spin-polarized first-principles calculations based on Density Functional Theory (DFT) using the Quantum ESPRESSO package.

• Model Construction: A 5 × 5 supercell of pristine graphene (50 carbon atoms) was used as the parent framework. A central void was created by removing specific carbon atoms, followed by the incorporation of nitrogen atoms at the defect boundary to form four distinct configurations: pyridinic, pyridazinic, pyrrolic, and pyrazolic.
• Computational Parameters: Calculations utilized the Generalized Gradient Approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) functional and Projector Augmented Wave (PAW) pseudopotentials. Geometry optimizations were performed using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm with variable cell relaxation.
• Analysis Techniques: The study evaluated structural stability via defect formation energy, analyzed local lattice reconstruction through bond-length variations, and quantified charge redistribution using charge-density difference maps and Bader charge analysis. Electronic and magnetic properties were assessed through spin-resolved band structures, spin-polarized projected density of states (SPDOS), and spin-density mapping.
• Validation: The computational framework was validated against a standard substitutional nitrogen-doped graphene model, reproducing known formation energies and confirming a non-magnetic ground state for substitutional nitrogen.

Key Contributions and Results

1. Structural Stability and Reconstruction:

   • Among the four motifs, the pyridinic configuration exhibits the most favorable thermodynamic stability, possessing the lowest defect formation energy (0.8295 Ry).
   • Pyridazinic graphene shows the highest structural distortion, with C-C bond lengths varying by 0.403 Å, attributed to asymmetric nitrogen positioning and local strain relief.
   • Pyrazolic graphene, despite having perfectly symmetric C-N bond lengths (1.398 Å), displays the highest formation energy (1.1912 Ry), indicating that structural symmetry does not guarantee energetic stability.
   • Pyrrolic graphene shows the least C-C bond distortion (0.211 Å) but possesses a moderate formation energy.

2. Electronic Structure and Band-Gap Engineering:

   • Metallic/Near-Metallic Behavior: The pyridinic, pyridazinic, and pyrrolic configurations retain metallic or near-metallic character, with defect-derived states crossing the Fermi level. These states arise from the hybridization of N 2p and neighboring C 2p orbitals.
   • Semiconducting Switch: In contrast, the pyrazolic configuration opens a narrow band gap of 0.377 eV (at the PBE level). This transition is attributed to the symmetric coordination of the nitrogen pair, which effectively passivates dangling vacancy-edge states and separates occupied and unoccupied defect-derived states.

3. Charge Redistribution:

   • Peripheral nitrogen atoms act as electron-accumulating centers.
   • Pyridinic and pyrazolic configurations show relatively balanced charge transfer between nitrogen sites.
   • Pyrrolic configurations exhibit significant charge asymmetry between the two nitrogen sites (-1.07 and -2.61), leading to non-equivalent local electronic environments.
   • Pyridazinic configurations show localized charge accumulation due to the proximity of the nitrogen pair.

4. Magnetic Properties:

   • Finite Magnetism: Pyridinic (1.35 μB), pyridazinic (1.49 μB), and pyrrolic (1.22 μB) configurations possess finite net magnetic moments. The magnetism originates from spin-polarized states involving both N 2p and neighboring C 2p orbitals, with the carbon network playing a major role in carrying the local spin density.
   • Spin Compensation: The pyrazolic configuration is uniquely spin-compensated, resulting in zero net magnetization. This non-magnetic state correlates directly with the opening of the band gap.

Significance and Claims
The paper establishes that the topology of peripheral nitrogen atoms around a graphene void acts as a "defect-chemical switch." The study demonstrates that by controlling the local heterocyclic-like nitrogen coordination (specifically choosing between pyridinic/pyridazinic/pyrrolic vs. pyrazolic motifs), one can deterministically tune defective graphene between:

• Magnetic, metallic/near-metallic states (pyridinic, pyridazinic, pyrrolic).
• Non-magnetic, narrow-gap semiconducting states (pyrazolic).

The authors conclude that these findings provide a structure-charge-spin relationship linking local nitrogen arrangement to macroscopic properties. This offers a metal-free route to pattern active spintronic and semiconducting domains directly into the graphene lattice through controlled vacancy-edge nitrogen coordination, addressing the limitations of pristine graphene regarding intrinsic band gaps and magnetism. The work suggests that future experimental and theoretical efforts could focus on exchange coupling and Curie temperature estimates to bridge these local defect states with long-range magnetic behavior.