Electron enrichment of zigzag edges of armchair-oriented graphene nano-ribbons increases their stability and induces pinning of Fermi level

E. Louis; E. San-Fabian; G. Chiappe; J.A. Verges

arxiv: 1907.01467 · v2 · pith:DDF4XU7Hnew · submitted 2019-07-02 · ❄️ cond-mat.mes-hall · cond-mat.mtrl-sci

Electron enrichment of zigzag edges of armchair-oriented graphene nano-ribbons increases their stability and induces pinning of Fermi level

E. Louis , E. San-Fabian , G. Chiappe , J.A. Verges This is my paper

Pith reviewed 2026-05-25 10:48 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall cond-mat.mtrl-sci

keywords graphene nanoribbonszigzag edgeselectron enrichmentFermi level pinningspin polarizationdensity functional theoryPariser-Parr-Pople model

0 comments

The pith

Electron enrichment at zigzag edges stabilizes armchair graphene nanoribbons and pins the Fermi level.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Armchair-oriented graphene nano-ribbons have zigzag edges with localized states that act like free radicals and can capture electrons. Calculations show that adding electrons lowers the total energy, reaching a minimum that depends primarily on the size of the zigzag edges. This charging removes spin polarization from the ground state, places the fundamental gap between the polarized and non-polarized neutral cases, and causes the valence and conduction bands to approach mid-gap with slopes of 0.1 and 0.9, pinning the Fermi level in line with experimental observations. The compatibility of a model Hamiltonian with density functional theory enables study of larger systems.

Core claim

Zigzag edges of neutral armchair-oriented Graphene Nano-Ribbons show states strongly localized at those edges that behave as free radicals capturing electrons during processing. Total energy calculations using spin polarized DFT show charging is feasible, and Pariser-Parr-Pople model energies are compatible allowing larger systems. Total energy decreases with captured electrons to a minimum depending on zigzag edges size. Charging makes ground state non spin polarized, fundamental gap intermediate, and induces Fermi level pinning with valence and conduction band slopes of 0.1 and 0.9 versus the gap, unlike the symmetric 0.5 in neutral ribbons.

What carries the argument

Electron capture at localized zigzag edge states, studied via spin-polarized DFT and the Pariser-Parr-Pople model Hamiltonian.

If this is right

Total energy of the ribbons reaches a minimum at a number of captured electrons that mainly depends on zigzag edge size.
The ground state of charged ribbons is not spin polarized.
The fundamental gap of charged ribbons lies between the gaps of spin-polarized and non-polarized neutral solutions.
Charging changes the band onsets versus fundamental gap from symmetric linear approach with slope 0.5 to asymmetric with slopes 0.1 and 0.9, pinning the Fermi level.
Results agree with experimental data on the gap and Fermi pinning.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The stabilization by electron enrichment could be relevant for processing and device fabrication of graphene nanoribbons.
The size dependence suggests that optimal charging varies with ribbon dimensions, potentially allowing control of electronic properties.
Extension to even larger systems using the model could test scaling of the pinning effect.

Load-bearing premise

The energies obtained from the Pariser-Parr-Pople model Hamiltonian are compatible with density functional theory calculations, permitting the study of larger ribbon systems.

What would settle it

A measurement showing that the total energy does not minimize at an electron count dependent on zigzag edge size, or that the band slopes upon charging are not approximately 0.1 and 0.9.

Figures

Figures reproduced from arXiv: 1907.01467 by E. Louis, E. San-Fabian, G. Chiappe, J.A. Verges.

**Figure 2.** Figure 2: (Color online) Six spin up molecular orbitals arou [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: (Color online) Six spin up molecular orbitals arou [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

**Figure 4.** Figure 4: (Color online) Energies of charged {m, n}-GNRs referred to that of the respective neutral ribbons, calculated by means of HF-PPP versus those obtained within the DFT framework and the combination functional/wavefunction set B3LYP/6-31+G*-SD. The results correspond to the ribbons included in [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗

**Figure 5.** Figure 5: (Color online) Total energies (all in eV) of severa [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗

**Figure 6.** Figure 6: (Color online) HOMO (circles) and LUMO (squares) e [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗

**Figure 7.** Figure 7: (Color online) Fundamental gap for neutral and char [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗

read the original abstract

Zigzag edges of neutral armchair-oriented Graphene Nano-Ribbons show states strongly localized at those edges. They behave as free radicals that can capture electrons during processing, increasing ribbon's stability. Thus, charging and its consequences should be investigated.Total energy calculations of finite ribbons using spin polarized Density Functional Theory (DFT) show that ribbon's charging is feasible. Energies for Pariser-Parr-Pople (PPP) model Hamiltonian are compatible with DFT allowing the study of larger systems. Results for neutral ribbons indicate: i) the fundamental gap of spin polarized (non polarized) solutions is larger (smaller) than experimental data, ii) the ground state is spin polarized, a characteristic still not observed experimentally. Total energy of GNRs decreases with the number of captured electrons reaching a minimum for a number that mainly depends on zigzag edges size. The following changes with respect to neutral GNRs are noted: i) the ground state is not spin polarized, ii) fundamental gap is in-between that of spin polarized and non polarized solutions of neutral ribbons, iii) while in neutral ribbons valence and conduction band onsets vs. the fundamental gap, linearly and symmetrically approach mid-gap with slope 0.5, charging induces Fermi level pinning, i.e., the slopes of the valence and conduction bands being about 0.1 and 0.9, in agreement with experiment.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Charging lowers total energy in these GNRs to a minimum set by zigzag edge length and produces asymmetric Fermi pinning slopes of ~0.1/0.9, but the DFT-to-PPP step for charged cases lacks shown validation.

read the letter

The main point is that adding electrons to the zigzag edges of armchair GNRs lowers the total energy until it reaches a minimum that depends mainly on edge size, turns off spin polarization, and pins the Fermi level so the valence and conduction onsets move with slopes near 0.1 and 0.9 instead of the symmetric 0.5 seen in neutral ribbons. The authors also note that the charged gap sits between the spin-polarized and non-polarized neutral gaps, and they say this pinning matches experiment. They reach the larger systems by switching from spin-polarized DFT on small finite ribbons to the PPP model after checking that the energies line up on neutral cases. That combination is the concrete new observation here. The calculations are direct total-energy minimizations with no fitted parameters, which is a plus. The soft spots are straightforward: the abstract and methods give no convergence tests, basis-set information, or error estimates, and the compatibility statement between DFT and PPP is not shown to cover the charged total-energy curves or the linear band slopes. If that check was only done on neutral small ribbons, the minimum location and the 0.1/0.9 slopes cannot be taken as firmly established for the larger ribbons. This is a narrow but usable piece for the GNR edge-stability crowd. A reader who runs their own DFT or tight-binding calculations on similar ribbons will want to see the numbers and test the slopes themselves. It deserves a serious referee once the authors add the missing technical checks on the charged systems; the core claim is simple enough that review would not take long.

Referee Report

2 major / 0 minor

Summary. The paper claims that zigzag edges in armchair-oriented graphene nano-ribbons act as free radicals that capture electrons, lowering total energy to a minimum that depends primarily on zigzag edge size. Using spin-polarized DFT on small ribbons and the PPP model on larger ones (asserting compatibility), it reports that charging eliminates spin polarization, yields an intermediate fundamental gap, and induces Fermi-level pinning with valence and conduction band onsets vs. gap having slopes ~0.1 and ~0.9 (vs. symmetric 0.5 in neutral case), matching experiment.

Significance. If the PPP-DFT compatibility holds for charged systems and the reported band slopes are robust, the work would provide a direct total-energy explanation for GNR edge stabilization during processing and for the experimentally observed Fermi pinning, without introducing fitted parameters. The approach of minimizing total energy with respect to captured charge is a strength.

major comments (2)

[Abstract/Methods] Abstract/Methods: No convergence tests, basis-set details, functional choices, or error estimates are supplied for the spin-polarized DFT total-energy calculations that locate the energy minimum versus captured electrons; without these the location of the minimum and the subsequent band-slope claims rest on unverified numerical protocols.
[Results] Results: The assertion that 'Energies for Pariser-Parr-Pople (PPP) model Hamiltonian are compatible with DFT' is not accompanied by explicit comparisons for charged ribbons or for the linear band-onset slopes under charging; this validation is load-bearing for extrapolating the energy-minimum location and the 0.1/0.9 pinning slopes to larger systems.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important issues regarding methodological transparency and validation that we address below. We will revise the manuscript to incorporate the requested details and comparisons.

read point-by-point responses

Referee: [Abstract/Methods] Abstract/Methods: No convergence tests, basis-set details, functional choices, or error estimates are supplied for the spin-polarized DFT total-energy calculations that locate the energy minimum versus captured electrons; without these the location of the minimum and the subsequent band-slope claims rest on unverified numerical protocols.

Authors: We agree that the manuscript as submitted lacks explicit documentation of these computational parameters. In the revised version we will add a Methods subsection (or appendix) specifying the DFT functional, basis set, convergence criteria, and any error estimates associated with the total-energy minima versus electron count. This will allow readers to assess the numerical robustness of the reported energy minimum and derived band slopes. revision: yes
Referee: [Results] Results: The assertion that 'Energies for Pariser-Parr-Pople (PPP) model Hamiltonian are compatible with DFT' is not accompanied by explicit comparisons for charged ribbons or for the linear band-onset slopes under charging; this validation is load-bearing for extrapolating the energy-minimum location and the 0.1/0.9 pinning slopes to larger systems.

Authors: The referee is correct that the current text does not supply side-by-side PPP versus DFT data for charged ribbons or for the band-onset slopes. We will include new supplementary figures or a table in the revision that directly compare total energies and band onsets between the two methods for representative charged systems, thereby justifying the extrapolation to larger ribbons and the reported 0.1/0.9 slopes. revision: yes

Circularity Check

0 steps flagged

No circularity: results from direct total-energy computations

full rationale

The paper's central results (energy minimum vs. captured electrons, Fermi pinning with valence/conduction slopes ~0.1/0.9) are obtained from explicit total-energy evaluations on finite ribbons using spin-polarized DFT (small systems) and the PPP Hamiltonian (larger systems). The abstract states compatibility between PPP and DFT energies but does not define any output quantity in terms of a fitted parameter or reduce a prediction to an input by construction. No self-citations, uniqueness theorems, or ansatzes are invoked as load-bearing steps. The reported slopes and minimum electron count are computed outputs, not renamings or self-referential definitions. This matches the default expectation of a non-circular computational study.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the compatibility of the PPP model with DFT results and on the validity of spin-polarized DFT for neutral ribbons; no free parameters or new entities are introduced in the abstract.

axioms (2)

domain assumption Energies for the Pariser-Parr-Pople (PPP) model Hamiltonian are compatible with DFT allowing the study of larger systems
Explicitly invoked to justify extending DFT results to larger ribbons.
domain assumption Spin-polarized DFT solutions are the appropriate reference for neutral ribbons
Used to compare fundamental gaps and ground-state character with charged cases.

pith-pipeline@v0.9.0 · 5799 in / 1491 out tokens · 38782 ms · 2026-05-25T10:48:14.943409+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Total energy of GNRs decreases with the number of captured electrons reaching a minimum for a number that mainly depends on zigzag edges size... charging induces Fermi level pinning, i.e., the slopes of the valence and conduction bands being about 0.1 and 0.9
IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Energies for Pariser-Parr-Pople (PPP) model Hamiltonian are compatible with DFT allowing the study of larger systems

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

[1]

Eftekhari, H

A. Eftekhari, H. Garcia, The necessity of structural irr egularities for the chemical applications of graphene, Mat erials Today Chemistry 4 (2017) 1 – 16. doi:10.1016/j.mtchem.2017.02.003. URL http://www.sciencedirect.com/science/article/pii/S2468519417300034

work page doi:10.1016/j.mtchem.2017.02.003 2017
[2]

Solis-Fernandez, M

P . Solis-Fernandez, M. Bissett, H. Ago, Synthesis, stru cture and applications of graphene-based 2d heterostructu res, Chem. Soc. Rev. 46 (2017) 4572–4613. doi:10.1039/C7CS00160F. URL http://dx.doi.org/10.1039/C7CS00160F

work page doi:10.1039/c7cs00160f 2017
[3]

Y . Zhu, H. Ji, H.-M. Cheng, R. S. Ruo ﬀ, Mass production and industrial applications of graphene m aterials, National Science Review 5 (1) (2018) 90–101. doi:10.1093/nsr/nwx055. URL http://dx.doi.org/10.1093/nsr/nwx055

work page doi:10.1093/nsr/nwx055 2018
[4]

S. Wang, L. Talirz, C. A. Pignedoli, X. Feng, K. M¨ ullen, R . Fasel, P . Ru ﬃeux, Giant edge state splitting at atomically precise graph ene zigzag edges, Nature Communications 7 (2016) 11507. URL http://dx.doi.org/10.1038/ncomms11507http://10.0.4.14/ncomms11507https://www.nature.com/articles/ ncomms11507#supplementary-information

work page doi:10.1038/ncomms11507http://10.0.4.14/ncomms11507https://www.nature.com/articles/ 2016
[5]

Leopold, R

T. Leopold, R. Pascal, F. Roman, On-surface synthesis of atomically precise graphene nanoribbons, Advanced Materi als 28 (29) (2016) 6222–6231. doi:10.1002/adma.201505738. URL https://onlinelibrary.wiley.com/doi/abs/10.1002/adma.201505738

work page doi:10.1002/adma.201505738 2016
[6]

W.-X. Wang, M. Zhou, X. Li, S.-Y . Li, X. Wu, W. Duan, L. He, E nergy gaps of atomically precise armchair graphene sidewal l nanoribbons, Phys. Rev. B 93 (2016) 241403. doi:10.1103/PhysRevB.93.241403. URL https://link.aps.org/doi/10.1103/PhysRevB.93.241403

work page doi:10.1103/physrevb.93.241403 2016
[7]

Xu, T.-W

W. Xu, T.-W. Lee, Recent progress in fabrication techniq ues of graphene nanoribbons, Mater. Horiz. 3 (2016) 186–207 . doi:10.1039/C5MH00288E. URL http://dx.doi.org/10.1039/C5MH00288E

work page doi:10.1039/c5mh00288e 2016
[8]

Y eh, P .-Y

C.-N. Y eh, P .-Y . Lee, J.-D. Chai, Electronic and Optical Properties of the Narrowest Armchair Graphene Nanoribbons Studied by Density Functional Methods, Australian Journal of Chemistry 69 (9) (2016) 960–968. URL https://doi.org/10.1071/CH16187

work page doi:10.1071/ch16187 2016
[9]

Wu, J.-D

C.-S. Wu, J.-D. Chai, Electronic properties of zigzag gr aphene nanoribbons studied by tao-dft, Journal of Chemical Theory and Computation 11 (5) (2015) 2003–2011. doi:10.1021/ct500999m. URL https://doi.org/10.1021/ct500999m

work page doi:10.1021/ct500999m 2015
[10]

S. Li, C. K. Gan, Y .-W. Son, Y . P . Feng, S. Y . Quek, Anomalous length-independent frontier resonant transmission peaks in armchair graphene nanoribbon molecular wires, Carbon 76 (2014) 285 – 291. doi:10.1016/j.carbon.2014.04.079. URL http://www.sciencedirect.com/science/article/pii/S0008622314004114

work page doi:10.1016/j.carbon.2014.04.079 2014
[11]

Talirz, H

L. Talirz, H. S¨ ode, T. Dumsla ﬀ, S. Wang, J. R. Sanchez-V alencia, J. Liu, P . Shinde, C. A. Pig nedoli, L. Liang, V . Meunier, N. C. Plumb, M. Shi, X. Feng, A. Narita, K. M¨ ullen, R. Fasel, P . Ru ﬃeux, On-surface synthesis and characterization of 9-atom w ide armchair graphene nanoribbons, ACS Nano 11 (2) (2017) 1380–1388, pMID: 281295 07. doi:10.102...

work page doi:10.1021/acsnano.6b06405 2017
[12]

Y .-C. Chen, D. G. de Oteyza, Z. Pedramrazi, C. Chen, F. R. Fischer, M. F. Crommie, Tuning the band gap of graphene nanoribbons synthesized from molecular precursors, ACS Nano 7 (7) (2013) 6123–6128, pMID: 23746141. doi:10.1021/nn401948e. URL https://doi.org/10.1021/nn401948e

work page doi:10.1021/nn401948e 2013
[13]

Talirz, H

L. Talirz, H. Sde, J. Cai, P . Ru ﬃeux, S. Blankenburg, R. Jafaar, R. Berger, X. Feng, K. Mllen, D. Passerone, R. Fasel, C. A. Pignedoli, Termini of bottom-up fabricated graphene nanoribbons, Jou rnal of the American Chemical Society 135 (6) (2013) 2060–20 63, pMID: 23350872. doi:10.1021/ja311099k. URL https://doi.org/10.1021/ja311099k

work page doi:10.1021/ja311099k 2013
[14]

Ru ﬃeux, J

P . Ru ﬃeux, J. Cai, N. C. Plumb, L. Patthey, D. Prezzi, A. Ferretti, E . Molinari, X. Feng, K. Mllen, C. A. Pignedoli, R. Fasel, Electronic structure of atomically precise graph ene nanoribbons, ACS Nano 6 (8) (2012) 6930–6935, pMID: 2285 3456. doi:10.1021/nn3021376. URL https://doi.org/10.1021/nn3021376

work page doi:10.1021/nn3021376 2012
[15]

M. Koch, F. Ample, C. Joachim, L. Grill, V oltage-depend ent conductance of a single graphene nanoribbon, Nature Nan otechnology 7 (2012)

work page 2012
[16]

URL http://dx.doi.org/10.1038/nnano.2012.169https://www.nature.com/articles/nnano.2012.169

work page doi:10.1038/nnano.2012.169https://www.nature.com/articles/nnano.2012.169 2012
[17]

P . B. Bennett, Z. Pedramrazi, A. Madani, Y .-C. Chen, D. G . de Oteyza, C. Chen, F. R. Fischer, M. F. Crommie, J. Bokor, Bottom-up graphene nanoribbon ﬁeld-e ﬀect transistors, Applied Physics Letters 103 (25) (2013) 25 3114. doi:10.1063/1.4855116. URL https://doi.org/10.1063/1.4855116

work page doi:10.1063/1.4855116 2013
[18]

Hauptmann, The aromatic sextet: V on e

S. Hauptmann, The aromatic sextet: V on e. clar; john wil ey & sons ltd, london, new york, sydney, toronto 1972; 128 sei ten mit zahlreichen formelbildern; format 13 20 cm; broschiert £ 1, 50, Zeitschr ift fr Chemie 13 (5) (1973) 200–200. doi:10.1002/zfch.19730130531. URL https://onlinelibrary.wiley.com/doi/abs/10.1002/zfch.19730130531

work page doi:10.1002/zfch.19730130531 1972
[19]

Sol` a, Forty years of Clar’s aromatic π-sextet rule, Frontiers in chemistry 1 (2013) 22

M. Sol` a, Forty years of Clar’s aromatic π-sextet rule, Frontiers in chemistry 1 (2013) 22. doi:10.3389/fchem.2013.00022. URL https://www.ncbi.nlm.nih.gov/pubmed/24790950https://www.ncbi.nlm.nih.gov/pmc/PMC3982536/

work page doi:10.3389/fchem.2013.00022 2013
[20]

Y ang, C.-H

L. Y ang, C.-H. Park, Y .-W. Son, M. L. Cohen, S. G. Louie, Q uasiparticle energies and band gaps in graphene nanoribbon s, Phys. Rev. Lett. 99 (2007) 186801. doi:10.1103/PhysRevLett.99.186801. URL https://link.aps.org/doi/10.1103/PhysRevLett.99.186801

work page doi:10.1103/physrevlett.99.186801 2007
[21]

O. V . Y azyev, Emergence of magnetism in graphene materi als and nanostructures, Reports on Progress in Physics 73 (5 ) (2010) 056501. doi:10.1088/0034-4885/73/5/056501. URL https://doi.org/10.1088%2F0034-4885%2F73%2F5%2F056501

work page doi:10.1088/0034-4885/73/5/056501 2010
[22]

Golor, C

M. Golor, C. Koop, T. C. Lang, S. Wessel, M. J. Schmidt, Ma gnetic correlations in short and narrow graphene armchair n anoribbons, Phys. 14 Rev. Lett. 111 (2013) 085504. doi:10.1103/PhysRevLett.111.085504. URL https://link.aps.org/doi/10.1103/PhysRevLett.111.085504

work page doi:10.1103/physrevlett.111.085504 2013
[23]

A. D. Zdetsis, E. N. Economou, A pedestrian approach to t he aromaticity of graphene and nanographene: Signiﬁcance o f huckels (4n +2)π electron rule, The Journal of Physical Chemistry C 119 (29) ( 2015) 16991–17003. doi:10.1021/acs.jpcc.5b04311. URL https://doi.org/10.1021/acs.jpcc.5b04311

work page doi:10.1021/acs.jpcc.5b04311 2015
[24]

Merino-D´ ıez, A

N. Merino-D´ ıez, A. Garcia-Lekue, E. Carbonell-Sanro m` a, J. Li, M. Corso, L. Colazzo, F. Sedona, D. S´ anchez-Port al, J. I. Pascual, D. G. de Oteyza, Width-dependent band gap in armchair graphene na noribbons reveals fermi level pinning on au(111), ACS Nano 1 1 (11) (2017) 11661–11668, pMID: 29049879. doi:10.1021/acsnano.7b06765. URL https://doi.org/1...

work page doi:10.1021/acsnano.7b06765 2017
[25]

A. D. Becke, Density-functional exchange-energy appr oximation with correct asymptotic behavior, Phys. Rev. A 38 (6) (1988) 3098–3100

work page 1988
[26]

C. Lee, W. Y ang, R. G. Parr, Development of the colle-salvetti correlation-energy formula into a functional of the e lectron density, Phys. Rev. B 37 (1988) 785–789. doi:10.1103/PhysRevB.37.785. URL http://link.aps.org/doi/10.1103/PhysRevB.37.785

work page doi:10.1103/physrevb.37.785 1988
[27]

A. D. Becke, Density functional thermochemistry. iii. the role of exact exchange, J. Chem. Phys. 98 (7) (1993) 5648– 5652. doi:10.1063/1.464913. URL http://aip.scitation.org/doi/abs/10.1063/1.464913

work page doi:10.1063/1.464913 1993
[28]

P . J. Stephens, F. J. Devlin, C. F. Chabalowski, M. J. Fri sch, Ab initio calculation of vibrational absorption and ci rcular dichroism spectra using density functional force ﬁelds, J. Phys. Chem. 98 (45) (1994) 11623–11627. doi:10.1021/j100096a001. URL http://dx.doi.org/10.1021/j100096a001

work page doi:10.1021/j100096a001 1994
[29]

P . C. Hariharan, J. A. Pople, The inﬂuence of polarizati on functions on molecular orbital hidrogenation energies, Theor. Chim. Acta 28 (1973) 213

work page 1973
[30]

Clark, J

T. Clark, J. Chandrasekhar, G. W. Spitznagel, P . V . R. Schleyer, Eﬃcient diﬀuse function-augmented basis sets for anion calculations. iii. the 3-21+g basis set for ﬁrst-row elements, Li–F, J. Comput. Chem. 4 (3 ) (1983) 294 – 301. doi:10.1002/jcc.540040303. URL http://dx.doi.org/10.1002/jcc.540040303

work page doi:10.1002/jcc.540040303 1983
[31]

M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria , M. A. Robb, J. R. Cheeseman, G. Scalmani, V . Barone, B. Mennu cci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P . Hratchian , A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Had a, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y . Honda, O. Kitao, H. ...

work page
[32]

J. Stewart, Optimization of parameters for semiempiri cal methods v: Modiﬁcation of nddo approximations and appli cation to 70 elements, Journal of Molecular Modeling 13 (12) (2007) 1173–1213. doi:10.1007/s00894-007-0233-4 . URL http://dx.doi.org/10.1007/s00894-007-0233-4

work page doi:10.1007/s00894-007-0233-4 2007
[33]

M. J. S. Dewar, W. Thiel, Ground states of molecules. 38. the mndo method. approximations and parameters, J. Am. Chem . Soc. 99 (1977) 4899

work page 1977
[34]

J. A. V erg´ es, E. SanFabi´ an, G. Chiappe, E. Louis, Fit o f pariser-parr-pople and hubbard model hamiltonians to cha rge and spin states of polycyclic aromatic hydrocarbons, Physical Review B 81 (8) (2010) 085120, pT: J; TC: 6; UT: WOS:000275053300047

work page 2010
[35]

Chiappe, E

G. Chiappe, E. Louis, E. San-Fabi´ an and J. A. V erg´ es, Can model Hamiltonians describe the electron–electron inte raction in π -conjugated systems?: PAH and graphene, Journal of Physics: Condensed M atter 27 (46) (2015) 463001. URL http://stacks.iop.org/0953-8984/27/i=46/a=463001

work page 2015
[36]

A. D. Zdetsis, E. Economou, Rationalizing and reconcil ing energy gaps and quantum conﬁnement in narrow atomically precise armchair graphene nanoribbons, Carbon 116 (2017) 422 – 434. doi:10.1016/j.carbon.2017.02.006. URL http://www.sciencedirect.com/science/article/pii/S0008622317301240

work page doi:10.1016/j.carbon.2017.02.006 2017
[37]

J. A. V erg´ es, G. Chiappe, E. San-Fabi´ an, E. Louis, Conductance through the armchair graphene nanoribbons 9-agnr : Strong dependence on contact to leads, Phys. Rev. B 98 (2018) 155415. doi:10.1103/PhysRevB.98.155415. URL https://link.aps.org/doi/10.1103/PhysRevB.98.155415 15

work page doi:10.1103/physrevb.98.155415 2018

[1] [1]

Eftekhari, H

A. Eftekhari, H. Garcia, The necessity of structural irr egularities for the chemical applications of graphene, Mat erials Today Chemistry 4 (2017) 1 – 16. doi:10.1016/j.mtchem.2017.02.003. URL http://www.sciencedirect.com/science/article/pii/S2468519417300034

work page doi:10.1016/j.mtchem.2017.02.003 2017

[2] [2]

Solis-Fernandez, M

P . Solis-Fernandez, M. Bissett, H. Ago, Synthesis, stru cture and applications of graphene-based 2d heterostructu res, Chem. Soc. Rev. 46 (2017) 4572–4613. doi:10.1039/C7CS00160F. URL http://dx.doi.org/10.1039/C7CS00160F

work page doi:10.1039/c7cs00160f 2017

[3] [3]

Y . Zhu, H. Ji, H.-M. Cheng, R. S. Ruo ﬀ, Mass production and industrial applications of graphene m aterials, National Science Review 5 (1) (2018) 90–101. doi:10.1093/nsr/nwx055. URL http://dx.doi.org/10.1093/nsr/nwx055

work page doi:10.1093/nsr/nwx055 2018

[4] [4]

S. Wang, L. Talirz, C. A. Pignedoli, X. Feng, K. M¨ ullen, R . Fasel, P . Ru ﬃeux, Giant edge state splitting at atomically precise graph ene zigzag edges, Nature Communications 7 (2016) 11507. URL http://dx.doi.org/10.1038/ncomms11507http://10.0.4.14/ncomms11507https://www.nature.com/articles/ ncomms11507#supplementary-information

work page doi:10.1038/ncomms11507http://10.0.4.14/ncomms11507https://www.nature.com/articles/ 2016

[5] [5]

Leopold, R

T. Leopold, R. Pascal, F. Roman, On-surface synthesis of atomically precise graphene nanoribbons, Advanced Materi als 28 (29) (2016) 6222–6231. doi:10.1002/adma.201505738. URL https://onlinelibrary.wiley.com/doi/abs/10.1002/adma.201505738

work page doi:10.1002/adma.201505738 2016

[6] [6]

W.-X. Wang, M. Zhou, X. Li, S.-Y . Li, X. Wu, W. Duan, L. He, E nergy gaps of atomically precise armchair graphene sidewal l nanoribbons, Phys. Rev. B 93 (2016) 241403. doi:10.1103/PhysRevB.93.241403. URL https://link.aps.org/doi/10.1103/PhysRevB.93.241403

work page doi:10.1103/physrevb.93.241403 2016

[7] [7]

Xu, T.-W

W. Xu, T.-W. Lee, Recent progress in fabrication techniq ues of graphene nanoribbons, Mater. Horiz. 3 (2016) 186–207 . doi:10.1039/C5MH00288E. URL http://dx.doi.org/10.1039/C5MH00288E

work page doi:10.1039/c5mh00288e 2016

[8] [8]

Y eh, P .-Y

C.-N. Y eh, P .-Y . Lee, J.-D. Chai, Electronic and Optical Properties of the Narrowest Armchair Graphene Nanoribbons Studied by Density Functional Methods, Australian Journal of Chemistry 69 (9) (2016) 960–968. URL https://doi.org/10.1071/CH16187

work page doi:10.1071/ch16187 2016

[9] [9]

Wu, J.-D

C.-S. Wu, J.-D. Chai, Electronic properties of zigzag gr aphene nanoribbons studied by tao-dft, Journal of Chemical Theory and Computation 11 (5) (2015) 2003–2011. doi:10.1021/ct500999m. URL https://doi.org/10.1021/ct500999m

work page doi:10.1021/ct500999m 2015

[10] [10]

S. Li, C. K. Gan, Y .-W. Son, Y . P . Feng, S. Y . Quek, Anomalous length-independent frontier resonant transmission peaks in armchair graphene nanoribbon molecular wires, Carbon 76 (2014) 285 – 291. doi:10.1016/j.carbon.2014.04.079. URL http://www.sciencedirect.com/science/article/pii/S0008622314004114

work page doi:10.1016/j.carbon.2014.04.079 2014

[11] [11]

Talirz, H

L. Talirz, H. S¨ ode, T. Dumsla ﬀ, S. Wang, J. R. Sanchez-V alencia, J. Liu, P . Shinde, C. A. Pig nedoli, L. Liang, V . Meunier, N. C. Plumb, M. Shi, X. Feng, A. Narita, K. M¨ ullen, R. Fasel, P . Ru ﬃeux, On-surface synthesis and characterization of 9-atom w ide armchair graphene nanoribbons, ACS Nano 11 (2) (2017) 1380–1388, pMID: 281295 07. doi:10.102...

work page doi:10.1021/acsnano.6b06405 2017

[12] [12]

Y .-C. Chen, D. G. de Oteyza, Z. Pedramrazi, C. Chen, F. R. Fischer, M. F. Crommie, Tuning the band gap of graphene nanoribbons synthesized from molecular precursors, ACS Nano 7 (7) (2013) 6123–6128, pMID: 23746141. doi:10.1021/nn401948e. URL https://doi.org/10.1021/nn401948e

work page doi:10.1021/nn401948e 2013

[13] [13]

Talirz, H

L. Talirz, H. Sde, J. Cai, P . Ru ﬃeux, S. Blankenburg, R. Jafaar, R. Berger, X. Feng, K. Mllen, D. Passerone, R. Fasel, C. A. Pignedoli, Termini of bottom-up fabricated graphene nanoribbons, Jou rnal of the American Chemical Society 135 (6) (2013) 2060–20 63, pMID: 23350872. doi:10.1021/ja311099k. URL https://doi.org/10.1021/ja311099k

work page doi:10.1021/ja311099k 2013

[14] [14]

Ru ﬃeux, J

P . Ru ﬃeux, J. Cai, N. C. Plumb, L. Patthey, D. Prezzi, A. Ferretti, E . Molinari, X. Feng, K. Mllen, C. A. Pignedoli, R. Fasel, Electronic structure of atomically precise graph ene nanoribbons, ACS Nano 6 (8) (2012) 6930–6935, pMID: 2285 3456. doi:10.1021/nn3021376. URL https://doi.org/10.1021/nn3021376

work page doi:10.1021/nn3021376 2012

[15] [15]

M. Koch, F. Ample, C. Joachim, L. Grill, V oltage-depend ent conductance of a single graphene nanoribbon, Nature Nan otechnology 7 (2012)

work page 2012

[16] [16]

URL http://dx.doi.org/10.1038/nnano.2012.169https://www.nature.com/articles/nnano.2012.169

work page doi:10.1038/nnano.2012.169https://www.nature.com/articles/nnano.2012.169 2012

[17] [17]

P . B. Bennett, Z. Pedramrazi, A. Madani, Y .-C. Chen, D. G . de Oteyza, C. Chen, F. R. Fischer, M. F. Crommie, J. Bokor, Bottom-up graphene nanoribbon ﬁeld-e ﬀect transistors, Applied Physics Letters 103 (25) (2013) 25 3114. doi:10.1063/1.4855116. URL https://doi.org/10.1063/1.4855116

work page doi:10.1063/1.4855116 2013

[18] [18]

Hauptmann, The aromatic sextet: V on e

S. Hauptmann, The aromatic sextet: V on e. clar; john wil ey & sons ltd, london, new york, sydney, toronto 1972; 128 sei ten mit zahlreichen formelbildern; format 13 20 cm; broschiert £ 1, 50, Zeitschr ift fr Chemie 13 (5) (1973) 200–200. doi:10.1002/zfch.19730130531. URL https://onlinelibrary.wiley.com/doi/abs/10.1002/zfch.19730130531

work page doi:10.1002/zfch.19730130531 1972

[19] [19]

Sol` a, Forty years of Clar’s aromatic π-sextet rule, Frontiers in chemistry 1 (2013) 22

M. Sol` a, Forty years of Clar’s aromatic π-sextet rule, Frontiers in chemistry 1 (2013) 22. doi:10.3389/fchem.2013.00022. URL https://www.ncbi.nlm.nih.gov/pubmed/24790950https://www.ncbi.nlm.nih.gov/pmc/PMC3982536/

work page doi:10.3389/fchem.2013.00022 2013

[20] [20]

Y ang, C.-H

L. Y ang, C.-H. Park, Y .-W. Son, M. L. Cohen, S. G. Louie, Q uasiparticle energies and band gaps in graphene nanoribbon s, Phys. Rev. Lett. 99 (2007) 186801. doi:10.1103/PhysRevLett.99.186801. URL https://link.aps.org/doi/10.1103/PhysRevLett.99.186801

work page doi:10.1103/physrevlett.99.186801 2007

[21] [21]

O. V . Y azyev, Emergence of magnetism in graphene materi als and nanostructures, Reports on Progress in Physics 73 (5 ) (2010) 056501. doi:10.1088/0034-4885/73/5/056501. URL https://doi.org/10.1088%2F0034-4885%2F73%2F5%2F056501

work page doi:10.1088/0034-4885/73/5/056501 2010

[22] [22]

Golor, C

M. Golor, C. Koop, T. C. Lang, S. Wessel, M. J. Schmidt, Ma gnetic correlations in short and narrow graphene armchair n anoribbons, Phys. 14 Rev. Lett. 111 (2013) 085504. doi:10.1103/PhysRevLett.111.085504. URL https://link.aps.org/doi/10.1103/PhysRevLett.111.085504

work page doi:10.1103/physrevlett.111.085504 2013

[23] [23]

A. D. Zdetsis, E. N. Economou, A pedestrian approach to t he aromaticity of graphene and nanographene: Signiﬁcance o f huckels (4n +2)π electron rule, The Journal of Physical Chemistry C 119 (29) ( 2015) 16991–17003. doi:10.1021/acs.jpcc.5b04311. URL https://doi.org/10.1021/acs.jpcc.5b04311

work page doi:10.1021/acs.jpcc.5b04311 2015

[24] [24]

Merino-D´ ıez, A

N. Merino-D´ ıez, A. Garcia-Lekue, E. Carbonell-Sanro m` a, J. Li, M. Corso, L. Colazzo, F. Sedona, D. S´ anchez-Port al, J. I. Pascual, D. G. de Oteyza, Width-dependent band gap in armchair graphene na noribbons reveals fermi level pinning on au(111), ACS Nano 1 1 (11) (2017) 11661–11668, pMID: 29049879. doi:10.1021/acsnano.7b06765. URL https://doi.org/1...

work page doi:10.1021/acsnano.7b06765 2017

[25] [25]

A. D. Becke, Density-functional exchange-energy appr oximation with correct asymptotic behavior, Phys. Rev. A 38 (6) (1988) 3098–3100

work page 1988

[26] [26]

C. Lee, W. Y ang, R. G. Parr, Development of the colle-salvetti correlation-energy formula into a functional of the e lectron density, Phys. Rev. B 37 (1988) 785–789. doi:10.1103/PhysRevB.37.785. URL http://link.aps.org/doi/10.1103/PhysRevB.37.785

work page doi:10.1103/physrevb.37.785 1988

[27] [27]

A. D. Becke, Density functional thermochemistry. iii. the role of exact exchange, J. Chem. Phys. 98 (7) (1993) 5648– 5652. doi:10.1063/1.464913. URL http://aip.scitation.org/doi/abs/10.1063/1.464913

work page doi:10.1063/1.464913 1993

[28] [28]

P . J. Stephens, F. J. Devlin, C. F. Chabalowski, M. J. Fri sch, Ab initio calculation of vibrational absorption and ci rcular dichroism spectra using density functional force ﬁelds, J. Phys. Chem. 98 (45) (1994) 11623–11627. doi:10.1021/j100096a001. URL http://dx.doi.org/10.1021/j100096a001

work page doi:10.1021/j100096a001 1994

[29] [29]

P . C. Hariharan, J. A. Pople, The inﬂuence of polarizati on functions on molecular orbital hidrogenation energies, Theor. Chim. Acta 28 (1973) 213

work page 1973

[30] [30]

Clark, J

T. Clark, J. Chandrasekhar, G. W. Spitznagel, P . V . R. Schleyer, Eﬃcient diﬀuse function-augmented basis sets for anion calculations. iii. the 3-21+g basis set for ﬁrst-row elements, Li–F, J. Comput. Chem. 4 (3 ) (1983) 294 – 301. doi:10.1002/jcc.540040303. URL http://dx.doi.org/10.1002/jcc.540040303

work page doi:10.1002/jcc.540040303 1983

[31] [31]

M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria , M. A. Robb, J. R. Cheeseman, G. Scalmani, V . Barone, B. Mennu cci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P . Hratchian , A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Had a, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y . Honda, O. Kitao, H. ...

work page

[32] [32]

J. Stewart, Optimization of parameters for semiempiri cal methods v: Modiﬁcation of nddo approximations and appli cation to 70 elements, Journal of Molecular Modeling 13 (12) (2007) 1173–1213. doi:10.1007/s00894-007-0233-4 . URL http://dx.doi.org/10.1007/s00894-007-0233-4

work page doi:10.1007/s00894-007-0233-4 2007

[33] [33]

M. J. S. Dewar, W. Thiel, Ground states of molecules. 38. the mndo method. approximations and parameters, J. Am. Chem . Soc. 99 (1977) 4899

work page 1977

[34] [34]

J. A. V erg´ es, E. SanFabi´ an, G. Chiappe, E. Louis, Fit o f pariser-parr-pople and hubbard model hamiltonians to cha rge and spin states of polycyclic aromatic hydrocarbons, Physical Review B 81 (8) (2010) 085120, pT: J; TC: 6; UT: WOS:000275053300047

work page 2010

[35] [35]

Chiappe, E

G. Chiappe, E. Louis, E. San-Fabi´ an and J. A. V erg´ es, Can model Hamiltonians describe the electron–electron inte raction in π -conjugated systems?: PAH and graphene, Journal of Physics: Condensed M atter 27 (46) (2015) 463001. URL http://stacks.iop.org/0953-8984/27/i=46/a=463001

work page 2015

[36] [36]

A. D. Zdetsis, E. Economou, Rationalizing and reconcil ing energy gaps and quantum conﬁnement in narrow atomically precise armchair graphene nanoribbons, Carbon 116 (2017) 422 – 434. doi:10.1016/j.carbon.2017.02.006. URL http://www.sciencedirect.com/science/article/pii/S0008622317301240

work page doi:10.1016/j.carbon.2017.02.006 2017

[37] [37]

J. A. V erg´ es, G. Chiappe, E. San-Fabi´ an, E. Louis, Conductance through the armchair graphene nanoribbons 9-agnr : Strong dependence on contact to leads, Phys. Rev. B 98 (2018) 155415. doi:10.1103/PhysRevB.98.155415. URL https://link.aps.org/doi/10.1103/PhysRevB.98.155415 15

work page doi:10.1103/physrevb.98.155415 2018