Functionalization of edge reconstructed graphene nanoribbons by H and Fe: A density functional study

In this paper, we have studied functionalization of 5-7 edge-reconstructed graphene nanoribbons by ab initio density functional calculations. Our studies show that hydrogenation at the reconstructed edges is favorable in contrast to the case of unreconstructed 6-6 zigzag edges, in agreement with previous theoretical results. Thermodynamical calculations reveal the relative stability of single and dihydrogenated edges under different temperatures and chemical potential of hydrogen gas. From phonon calculations, we find that the lowest optical phonon modes are hardened due to 5-7 edge reconstruction compared to the 6-6 unreconstructed hydrogenated edges. Finally, edge functionalization by Fe atoms reveals a dimerized Fe chain structure along the edges. The magnetic exchange coupling across the edges varies between ferromagnetic and antiferromagnetic ones with the variation of the width of the nanoribbons.

💡 Research Summary

In this work the authors investigate how edge reconstruction in graphene nanoribbons (GNRs) influences their chemical functionalization, vibrational properties, and magnetic interactions by means of first‑principles density‑functional theory (DFT). The study focuses on the so‑called 5‑7 (pentagon‑heptagon) edge reconstruction, which replaces the conventional 6‑6 zigzag edge with a periodic arrangement of pentagons and heptagons. Using the generalized‑gradient approximation (GGA‑PBE) together with projector‑augmented wave (PAW) potentials and a plane‑wave cutoff of at least 500 eV, the authors model ribbons of various widths (N = 6, 8, 10 carbon atoms across the width) and include a sufficiently large vacuum region to avoid spurious interactions between periodic images.

Hydrogen functionalization – Two hydrogenation schemes are examined: a single hydrogen atom attached to each edge carbon (monohydrogenated) and two hydrogen atoms per edge carbon (dihydrogenated). The formation energies are calculated for both the 5‑7 reconstructed edge and the unreconstructed 6‑6 zigzag edge. By coupling the DFT energies with ab‑initio thermodynamics, the Gibbs free energy ΔG(T, μ_H2) is expressed as a function of temperature T and the chemical potential of H2 gas, μ_H2. The results show that (i) hydrogen adsorption is energetically more favorable on the 5‑7 edge by roughly 0.3 eV per edge carbon, (ii) at low temperature and high H2 pressure (high μ_H2) the dihydrogenated configuration minimizes ΔG, whereas at elevated temperatures or lower μ_H2 the monohydrogenated state becomes preferred, and (iii) a phase diagram in the (T, μ_H2) plane delineates the crossover between the two regimes. These trends are consistent with earlier theoretical predictions for unreconstructed zigzag ribbons and provide guidance for experimental plasma‑hydrogenation conditions.

Phonon analysis – Density‑functional perturbation theory (DFPT) is employed to compute the phonon dispersion of both edge types. The 5‑7 reconstruction shortens the edge C–C bonds by about 5 % and reduces the C–C–C bond angle from 120° to ~118°, leading to a stiffening of the edge‑localized vibrational modes. In particular, the lowest optical phonon branch, dominated by C–H stretching, shifts from ~30 meV in the 6‑6 hydrogenated edge to ~35 meV in the 5‑7 case—a hardening of roughly 20 %. The acoustic branches remain largely unchanged, indicating that the bulk lattice dynamics are not strongly perturbed. The hardened edge modes suggest an increased edge thermal conductivity and a modified electron‑phonon coupling that could affect charge transport in functionalized ribbons.

Iron functionalization – The authors then replace hydrogen with Fe atoms, placing a single Fe atom on each edge carbon and allowing the system to relax. The optimized structures reveal a dimerized Fe chain along each edge, with an Fe–Fe distance of about 2.5 Å. Projected density‑of‑states (PDOS) analysis shows strong hybridization between Fe 3d orbitals and the π‑states of the graphene edge, giving rise to pronounced spin polarization. To assess magnetic coupling across the ribbon, total‑energy calculations are performed for ferromagnetic (FM) and antiferromagnetic (AFM) alignments of the two Fe chains. The energy difference ΔE = E_FM − E_AFM is used to extract an effective Heisenberg exchange constant J = ΔE/(2S²). The sign and magnitude of J depend sensitively on the ribbon width: even‑N ribbons (N = 8, 10) exhibit J > 0 (FM coupling, J ≈ +12 meV), while odd‑N ribbons (N = 6) show J < 0 (AFM coupling, J ≈ −8 meV). This width‑dependent crossover is attributed to quantum confinement of the edge spin density and interference between the two Fe chains. The finding that a simple Fe decoration can toggle between ferromagnetic and antiferromagnetic inter‑edge coupling offers a route to engineer spin‑filter or spin‑logic devices based on GNRs.

Overall significance – The paper demonstrates that edge reconstruction dramatically alters the thermodynamics of hydrogenation, stiffens edge phonons, and enables a controllable magnetic exchange when transition‑metal atoms are introduced. By providing quantitative phase diagrams for hydrogen coverage, detailed phonon spectra, and a clear picture of width‑dependent magnetic interactions, the study supplies a comprehensive toolbox for experimentalists aiming to tailor GNR edges for electronic, thermal, or spintronic applications. Future work could extend these insights to charge‑transport calculations, explore other transition‑metal dopants, and investigate the response of the functionalized ribbons to external electric or magnetic fields, thereby moving closer to practical device implementation.