Raman fingerprints of atomically precise graphene nanoribbons

Abstract: Bottom-up approaches allow the production of ultra-narrow and atomically precise graphene nanoribbons (GNRs), with electronic and optical properties controlled by the specific atomic structure. Combining Raman spectroscopy and ab-initio simulations, we show that GNR width, edge geometry and functional groups all influence their Raman spectra. The low-energy spectral region below 1000 cm-1 is particularly sensitive to edge morphology and functionalization, while the D peak dispersion can be used to uniquely fingerprint the presence of GNRs, and differentiates them from other sp2 carbon nanostructures.

• The paper presents a combined experimental and theoretical Raman analysis that reveals unique edge-dependent vibrational fingerprints in graphene nanoribbons.
• It employs multi-energy Raman spectroscopy and DFPT simulations to quantify the influence of edge geometry and alkyl functionalization on low-energy modes.
• Key findings include distinct D peak dispersions and RLBM red-shifts, enabling differentiation of GNRs from defective graphene and carbon nanotubes.

Raman Signatures of Atomically Precise Graphene Nanoribbons

Introduction

This work presents a comprehensive experimental and theoretical investigation of the Raman spectra of atomically precise, ultra-narrow graphene nanoribbons (GNRs) synthesized via solution-based bottom-up approaches. Emphasis is placed on the relationships among GNR width, edge morphology, and edge functionalization, and their manifestation in distinct Raman features, extending beyond the established characteristics observed in pristine graphene and other sp² carbon nanostructures such as carbon nanotubes (CNTs). The study systematically differentiates the Raman response of GNRs from that of defective graphene and other related systems, with a focus on the low-energy vibrational modes and the behavior of the D peak under varying excitation energy.

Experimental and Theoretical Framework

Four different GNR systems are analyzed, covering variations in width, edge geometry, and alkyl functionalization patterns. These include cove-shaped ribbons labeled 4CNR, 6CNR, and 8CNR, as well as a chiral-edged m-ANR. Raman spectroscopy is performed using multiple laser excitation energies and cross-validated across independent spectrometers and laboratories. Ab initio simulations are carried out within the framework of density functional perturbation theory (DFPT), leveraging an effective width model and detailed consideration of edge passivation and alkyl chain attachment.

Low-Energy Spectral Features: Sensitivity to Structural Parameters

Unlike CNTs and ideal zigzag GNRs, which possess a well-defined, sharp, low-energy radially breathing mode (RBM or RLBM), the GNRs studied here display a broader, sometimes multi-component peak, with substantial frequency red-shifts and broadening depending on the edge structure and functionalization. Experimental measurements reveal RLBM frequencies ranging from approximately 230 cm⁻¹ in 4CNR to 130–150 cm⁻¹ in 6CNR, 8CNR, and m-ANR. Ab initio calculations confirm that the RLBM is not solely a function of ribbon width; rather, it is sensitively modulated by the periodic edge patterning (e.g., cove or chiral edge structures) and by the specific placement of alkyl side chains.

The simulations demonstrate that side-chain functionalization induces a significant red-shift of the RLBM, an effect not attributable to simple increases in width but rather to the dynamic involvement of the chains in the breathing vibration. This gives rise to broadening and, for longer chains or certain widths, the appearance of multiple sub-peaks due to chain–ribbon vibrational coupling. The RLBM in these GNRs therefore acts as a multidimensional probe, encoding not only the lateral confinement but also edge morphology and functionalization.

High-Energy Raman Modes: D and G Characterization

The high-energy Raman response is dominated by features analogously labeled as the G and D peaks. The G peak, corresponding to the zone-center E₂g phonon, is observed at elevated frequencies (~1605–1618 cm⁻¹) relative to bulk graphene and exhibits notable broadening. This is ascribed to relaxation of the momentum conservation typically imposed in 2D graphene, a consequence of finite lateral size and edge disorder. No splitting into longitudinal and transverse components is detected, in contrast with CNTs, as predicted by the computation for these specific edge geometries.

The D peak exhibits a structured line-shape with dominant sub-bands arising from breathing modes of six-membered rings. Importantly, both experiment and simulation reveal that, although the D/G intensity ratio and D peak profile superficially resemble those in defective graphene, a distinct Raman fingerprint emerges upon analysis of D peak dispersion with excitation energy. Unlike graphene, where the D peak exhibits a strong linear dispersion (~50 cm⁻¹/eV) due to the underlying Kohn anomaly and double-resonance process, the GNRs show much weaker and nontrivial dispersions (ranging from ~7 to 35 cm⁻¹/eV depending on energy regime and GNR type). These unique dispersions are linked to the electronic structure of GNRs, which are semiconducting with substantial bandgaps and lack the Dirac dispersion at K.

Additionally, first-order activation of the D mode is permitted in these systems due to Brillouin zone folding caused by periodic edge patterning, and the G+D combination mode—absent in defective graphene—emerges as a further distinguishing feature.

Distinction from Other sp² Carbon Nanostructures

This combined experimental-theoretical approach enables a robust set of Raman fingerprints for solution-synthesized, functionalized, narrow GNRs:

• The RLBM is broadened, red-shifted, and compositionally split due to strong coupling with edge geometry and functionalization, in contrast to the sharp, width-only dependent RBM of CNTs and ideal GNRs.
• The G peak is upshifted, broadened, and not longitudinal-transverse split.
• The D peak displays unique, structure- and energy-dependent dispersion and a multi-component line-shape, not solely explainable by sp² ring breathing nor typical disorder activation.
• The observation of modes forbidden or absent in graphene (such as G+D combinations) via first-order processes directly traces to the distinct vibrational and electronic properties induced by edge engineering.
• D peak dispersion and the structure of low-frequency Raman modes together enable unambiguous discrimination of GNRs from both defective graphene and CNTs.

Implications and Outlook

This study provides a rigorous and multi-parametric methodology for the fingerprinting of GNRs, essential for quality control, design, and deployment in nanoelectronic, spintronic, and optoelectronic applications where the role of edge and functional configuration is paramount. The findings establish that neither the zone-folding approximation nor width-only scaling laws are valid for functionalized, atomically controlled GNRs; full vibrational calculations including functional groups and edge topology are required for accurate spectral interpretation.

Practically, the demonstrated Raman fingerprints set the groundwork for the high-throughput, non-destructive characterization of GNRs with atomic-scale precision, integrating seamlessly with large-area or solution-processed production techniques. Theoretically, these results highlight the profound interplay between finite size, edge periodicity, chemical functionalization, and vibrational properties in low-dimensional carbon systems, opening questions on the engineering of phonon–electron interactions via edge chemistry for tailored quantum transport and optical response. Further extensions may explore more complex edge functionalizations, the influence of substrate interactions, and correlation with other spectroscopic or transport properties.

Conclusion

The analysis underlines the necessity of advanced vibrational and electronic modeling beyond standard approximations for GNR characterization. Raman spectroscopy, in conjunction with first-principles calculations, emerges as a robust diagnostic for disentangling structural parameters in ultra-narrow, functionalized GNRs, setting benchmarks for their identification and application in emergent device architectures. The work delineates the structural intricacies encoding the vibrational landscape of GNRs, providing foundational insight for future studies targeting the controlled exploitation of their edge and size-dependent properties.