Document Detail

The Multiradical Character of One- and Two-Dimensional Graphene Nanoribbons.
Jump to Full Text
MedLine Citation:
PMID:  23355507     Owner:  NLM     Status:  Publisher    
When is an acene stable? The pronounced multiradical character of graphene nanoribbons of different size and shape was investigated with high-level multireference methods. Quantitative information based on the number of effectively unpaired electrons leads to specific estimates of the chemical stability of graphene nanostructures.
Felix Plasser; Hasan Pašalić; Martin H Gerzabek; Florian Libisch; Rafael Reiter; Joachim Burgdörfer; Thomas Müller; Ron Shepard; Hans Lischka
Related Documents :
21480767 - The effect of progesterone on the electron emission and degradation of testosterone.
21405427 - Ferromagnetically coupled magnetic impurities in a quantum point contact.
16626627 - Lateral energy transfer model for adjacent light-harvesting antennae rods of c-phycocya...
21264437 - Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly...
11580547 - Zero-magnetization ferromagnet proven by helicity-switching compton scattering.
19791857 - Maximum locality in occupied and virtual orbital spaces using a least-change strategy.
Publication Detail:
Type:  JOURNAL ARTICLE     Date:  2013-1-28
Journal Detail:
Title:  Angewandte Chemie (International ed. in English)     Volume:  -     ISSN:  1521-3773     ISO Abbreviation:  Angew. Chem. Int. Ed. Engl.     Publication Date:  2013 Jan 
Date Detail:
Created Date:  2013-1-28     Completed Date:  -     Revised Date:  -    
Medline Journal Info:
Nlm Unique ID:  0370543     Medline TA:  Angew Chem Int Ed Engl     Country:  -    
Other Details:
Languages:  ENG     Pagination:  -     Citation Subset:  -    
Copyright Information:
Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
Institute for Theoretical Chemistry, University of Vienna, Waehringerstrasse 17, 1090 Vienna (Austria).
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Angew Chem Int Ed Engl
Journal ID (iso-abbrev): Angew. Chem. Int. Ed. Engl
Journal ID (publisher-id): anie
ISSN: 1433-7851
ISSN: 1521-3773
Publisher: WILEY-VCH Verlag, Weinheim
Article Information
Download PDF
Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim
Received Day: 23 Month: 9 Year: 2012
Revision Received Day: 14 Month: 11 Year: 2012
Print publication date: Day: 25 Month: 2 Year: 2013
Electronic publication date: Day: 28 Month: 1 Year: 2013
Volume: 52 Issue: 9
First Page: 2581 Last Page: 2584
PubMed Id: 23355507
ID: 3648980
DOI: 10.1002/anie.201207671

The Multiradical Character of One- and Two-Dimensional Graphene Nanoribbons**
Felix Plasserau1
Hasan Pašalićau1
Martin H Gerzabekau2
Florian Libischau3
Rafael Reiterau3
Joachim Burgdörferau3
Thomas Müllerau4
Ron Shepardau5
Hans Lischka*au1au6
Institute for Theoretical Chemistry, University of ViennaWaehringerstrasse 17, 1090 Vienna (Austria)
Institute of Soil Research, University of Natural Resources and Life SciencesPeter-Jordan-Strasse 82, 1190 Vienna (Austria)
Institute for Theoretical Physics, Vienna University of TechnologyWiedner Hauptstrasse 8–10, 1040 Vienna (Austria)
Institute of Advanced Simulation, Jülich Supercomputer Centre, Forschungszentrum Jülich52425 Jülich (Germany)
Chemical Sciences and Engineering Division, Argonne National LaboratoryArgonne, IL 60439 (USA)
Department of Chemistry and Biochemistry, Texas Tech UniversityLubbock, TX 79409-1061 (USA)
Correspondence: Institute for Theoretical Chemistry, University of Vienna, Waehringerstrasse 17, 1090 Vienna (Austria) E-mail:
**This work was supported by the Austrian Science Fund (SFB F41, ViCoM, and Project P20893-N19), the Deutsche Forschungsgemeinschaft (SPP 1315, Project GE1676/1-1), the Robert A. Welch Foundation (grant no. D-0005), and the Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, U.S. Department of Energy (contract number DE-AC02-06CH11357 to H.L. and R.S.). We are grateful to T. J. Mach and T. D. Crawford for implementation of the D2 diagnostic into PSI4 and for performing the respective calculations, as well as to the Vienna Scientific Cluster (VSC), Project 70151 for computer time.
Supporting information for this article is available on the WWW under
Re-use of this article is permitted in accordance with the Terms and Conditions set out at

Since the first experimental realizations of graphene nanoribbons,1 graphene nanodevices2 have attracted enormous attention in the quest for future nanoscale technologies. Because of their small band gaps and high charge-carrier mobilities, n-acenes (Figure 1) and functionalized acenes are being considered as highly interesting building units for organic electronic materials.3 Zigzag nanoribbons feature remarkable spin-polarization and half-metallic properties4 and chemically fascinating challenges have to be solved in the synthesis of such extended polyaromatic hydrocarbons (PAHs).5 Moreover, the chemical properties of graphene nanoribbons may be linked to the reactivity of black carbon surfaces and PAHs in soils, affecting the accumulation of persistent organic pollutants with environmental implications.6

The availability of longer acenes is, however, hampered by their increasing reactivity, with pentacene being the largest well-characterized acene.3b In recent years substantial progress has resulted in the synthesis of n-acenes up to n=9 by matrix isolation techniques (see Ref. 7 and references therein). Nevertheless, these higher acenes are very reactive; for example, heptacene was found to be stable only for 4 h in a poly(methyl methacrylate) matrix.8 To overcome the stability problems, larger acenes were functionalized7a by adding protecting groups which inhibit the native high reactivity of the acenes.

Quantum chemical investigations play a major role in clarifying the outstanding electronic properties of carbon nanoflakes displaying bi- or even multiradical character. Density functional theory (DFT) constitutes a natural choice because of its good performance in terms of accuracy and computational efficiency with numerous applications to acenes9 and graphene flakes.10 However, because of the pronounced radical character, unrestricted DFT methods must be used,9a which are afflicted with spin contamination problems and the challenge of choosing a proper functional. As alternatives, density matrix renormalization group (DMRG)11 and the active-space variational two-electron reduced-density-matrix (2-RDM) approaches have been applied to one-12 and some two-dimensional13 systems. Additionally, spin-flip configuration interaction14 coupled cluster with singles and doubles (CCSD), along with more extended calculations, have been performed15 to investigate the electronic structure of the n-acenes.

For an adequate treatment of multiradical systems multireference (MR) methods16 are especially well suited. Moreover, MR methods are not affected by spin contamination or energy instability problems. In the present work we employ the MR averaged quadratic coupled cluster method (MR-AQCC)17 which includes crucial size-extensivity corrections at the MR level. It has already been successfully applied to several smaller challenging biradical systems.18 We show here its applicability to significantly larger molecular systems in the range of 100 carbon atoms in the graphene nanoflake.

The analysis of the radical character of PAHs as model systems for graphene nanoflakes will be performed by considering the natural orbital (NO) occupations as computed from the AQCC density by 1) following the deviations of individual NO occupations ni from zero (unoccupied) and two (doubly occupied), respectively, and by 2) computing a density and a number of effectively unpaired electrons NU as originally introduced by Takatsuka et al.19 as the distribution of “odd” electrons, which provides a measure for the splitting of an electron pair into different spatial regions (see also Ref. 20). We use the formalism of Ref. 21 and the nonlinear formula ni2(2−ni)2 given in this reference to compute the density of unpaired electrons; it more cleanly separates the dynamical correlation contributions to the wave function from the truly open-shell contributions of the radical centers. Specifically, relative to the linear function also suggested in Ref. 21, it emphasizes contributions from orbitals with occupations near 1, while also suppressing contributions with occupations near 0 or 2.

Two classes of molecular systems (Figure 1) have been investigated: 1) quasi-one-dimensional n-acenes and n-phenacenes and 2) two-dimensional periacenes and circumacenes. The latter will be characterized by the notation (ma, nz) where m and n count the number of benzene rings along each direction and a and z denote armchair and zigzag boundaries, respectively (see Figure 1 c,d). The latter two structures differ only in the details of their armchair edges leading, however, to significant changes in the radical character.

In our approach we follow the strategy of the DMRG11 and 2-RDM12 calculations in focusing on the conjugated π system. This is achieved by freezing the σ orbitals at the self-consistent field (SCF) level. The molecular orbitals (MOs) are obtained from multiconfiguration (MC) SCF calculations for use in the subsequent MR-AQCC approach. In all cases the totally symmetric singlet state is computed. The active orbital space for the reference configurations was chosen such that polyradical character and several open-shell electrons could be properly represented. Two reference spaces including up to 16 active orbitals have been selected for the MR-AQCC calculations to account for the exceptional requirements of the different compounds. The basis set used is of double-zeta quality. More details on the computational methods can be found in the Supporting Information.

The acene series (Figure 2) exhibits a dramatic increase of the open-shell NO occupations with chain length n, resulting in an entire cascade of open-shell occupations. For n=11 two NOs are almost degenerate (with occupations of 1.1 e and 0.90 e), but even the neighboring orbital occupations (1.6 e and 0.4 e) easily qualify as open-shell orbitals. The increasing instability of the acenes with chain length correlates well with the increasing multiradical character displayed in Figure 2. On the other hand, the NO occupations for the phenacenes (Figure S1 in the Supporting Information) do not show any indication of biradical character up to n=10, a chain length at which the acenes already possess strong multiradical character.

The evolution of strong multiradical character in our calculations agrees well with the DMRG11 and 2-RDM12 results, but is in contrast to the recent CCSD calculations15 where from the analysis of the T1 diagnostics it was concluded that the acenes are of pure closed-shell character. However, it has been discussed intensively in the literature22 that the single excitations involved in this diagnostic are associated with orbital relaxation and that the double excitations must be analyzed to assess the multireference character of a CCSD calculation. CCSD calculations using the D2 diagnostic developed for that purpose22b have been performed23 using the PSI4 program system.20 The D2 diagnostic (see Table S2 in the Supporting Information) considerably exceeds the recommended threshold; therefore, by this measure, the CCSD method also displays significant multireference character, in good agreement with our MR-AQCC calculations.

Extending this analysis from quasi one-dimensional quantum wires to true two-dimensional nanoflakes, we investigate both (3a, nz) and (5a, nz) periacenes (Figure S2 and Figure 3) and the circumacene sequences (Figure S3 and Figure 5). Similar to the quasi-1D acenes, the two-dimensional periacene sequence develops a strong multiradical character with increasing zigzag chain length. This evolution is more rapid than in the acene case. Increasing the armchair length in the periacenes from m=3 to 5 (compare Figure S2 and Figure 3) enhances the open-shell character considerably.

For the 5a series (Figure 3) complete degeneracy in the NO occupations is found already for n=4. The radical character is predominantly located at the zigzag edges as shown for the (5a, 6z) case (Figure 4). In this case the number of effectively unpaired electrons NU is 5.0 e. Individual values attributed to the carbon atoms based on a Mulliken analysis24 of the density of unpaired electronsare given in Figure 4 as well. They illustrate the strong radical character and the related chemical reactivity at the zigzag edges. The four orbitals of strong open-shell character for (5a, 6z) periacene contain about three open-shell electrons. The difference to the NU value arises from the contribution of the remaining minor open-shell occupation included in the density of unpaired electrons. Spin-polarized DFT calculations on rectangular nanoribbons using the Heyd, Scuseria, and Ernzerhof (HSE06) hybrid density functional show10a that spin polarization sets in already at a relatively small size of (3a, 3z). The MR-AQCC NO occupations (Figure S2) lead to a similar finding concerning this onset but give, additionally, a quantitative picture of the extent of the radical generation along the periacene series.

A notably delayed appearance of biradical character is found for the (3a, nz) and (5a, nz) circumacenes (Figure S3 and Figure 5, respectively) in comparison to the periacene pattern (Figure 3). In the latter case an NO pair with occupations of roughly 0.75 e and 1.25 e appears already for n=3, whereas in the (5a, nz) circumacene series comparable occupation values are obtained only for n=5. For higher members of the circumacene series eventually a similar cascade of open-shell occupations is found as for the periacenes.

The density of unpaired electrons for the (5a, 6z) circumacene (Figure S4) closely resembles the one for the corresponding (5a, 6z) periacene. Consequently, we predict the circumacenes to approach with increasing length the radical character and reactivity of the periacenes. Jiang et al.10d have used the simple HOMO–LUMO gap criterion to estimate a critical value of nC=6 for (3a, nz) circumacenes to reach open-shell character. This result is in qualitative agreement with ours, even though the numbers of effectively unpaired electrons discussed below indicate the transition to significant open-shell character occurs already between n=4 and n=5.

According to the NU analysis (Figure 6) the graphene nanoflakes can be subdivided into three groups. The n-phenacenes show a practically horizontal line with NU values of approximately 0.3–0.5 e. They exemplify the closed-shell character within our analysis. Next, the n-acenes show the strong increase in unpaired character reaching an NU value of 4 e for n=11. This value is in good agreement with the four open-shell NOs shown in Figure 2. The third class of compounds is formed by the two-dimensional graphene flakes. They all show an increase in the unpaired electron character which is notably stronger than that of the acenes up to n=9. There are differences in the behavior of periacenes and circumacenes: the strong increase in NU with the chain length n starts immediately for the periacenes whereas in particular the (3a, nz) circumacenes start off more slowly.

It is now tempting to correlate the NU values of different classes of graphene nanoflakes with available experimental evidence for the reactivity. Accordingly, we can identify a region of stability below the critical line NU(1)≍1 e passing through pentacene, the largest well-accessible member in the acene series.3b This region includes the periacenes (3a, 2z) (perylene) and (5a, 2z), and the circumacenes (3a, 1z) (pyrene), (3a, 2z) (coronene, circumbenzene), (3a, 3z) (ovalene, circumnaphthalene), (5a, 1z) and (5a, 2z), all of which are known stable compounds.25 The synthesis of circumanthracene has also been reported.26 A somewhat more tentative assignment is the region delimited by the dashed line NU(2)=2.2 e passing through nonacene, the largest acene synthesized so far.7b It is interesting to note that the NU values for the acenes show an additional increase at n=9 indicating an even heightened difficulty in preparing acenes with n>9. This second region at or below NU(2) would suggest that two-dimensional flakes including the (3a, 3z) and (5a, 3z) periacenes and circumacenes up to (5a, 4z) should be synthesizable under similar precautions as those used for the higher n-acenes. For larger 2D nanoflakes currently available,1b we conjecture that the different scaling of circumference versus area results in increased stability in the bulk region, as the density of unpaired electrons is restricted to the edges (Figure 4). The comparison of the properties of the different PAHs investigated in this work illustrates the remarkable variety in open-shell character when different pathways are followed in extending the molecular size of PAHs. The analysis of different PAH topologies in terms of structural and electronic properties is of considerable chemical interest.27 we have shown that MR methods provide powerful tools for obtaining reliable information on these fascinating compounds which are so difficult to access experimentally.

1a. Novoselov KS,Jiang D,Schedin F,Booth TJ,Khotkevich VV,Morozov SV,Geim AK. Proc. Natl. Acad. Sci. USAYear: 2005102104511045316027370
1b. Han MY,Ozyilmaz B,Zhang YB,Kim P. Phys. Rev. LettYear: 20079820680517677729
2. Geim AK,Novoselov KS. Nat. MaterYear: 2007618319117330084
3a. Bendikov M,Wudl F,Perepichka DF. Chem. RevYear: 20041044891494515535637
3b. Anthony JE. Angew. ChemYear: 2008120460492
Angew. Chem. Int. EdYear: 200847452483
4. Son YW,Cohen ML,Louie SG. NatureYear: 200644434734917108960
5. Cai JM,Ruffieux P,Jaafar R,Bieri M,Braun T,Blankenburg S,Muoth M,Seitsonen AP,Saleh M,Feng XL,Müllen K,Fasel R. NatureYear: 201046647047320651687
6. Heymann K,Lehmann J,Solomon D,Schmidt MWI,Regier T. Org. GeochemYear: 20114210551064
7a. Zade SS,Bendikov M. Angew. ChemYear: 201012241044107
Angew. Chem. Int. EdYear: 20104940124015
7b. Tönshoff C,Bettinger HF. Angew. ChemYear: 201012242194222
Angew. Chem. Int. EdYear: 20104941254128
8. Mondal R,Shah BK,Neckers DC. J. Am. Chem. SocYear: 20061289612961316866498
9a. Bendikov M,Duong HM,Starkey K,Houk KN,Carter EA,Wudl F. J. Am. Chem. SocYear: 20041267416741715198569
9b. Jiang DE,Dai S. J. Phys. Chem. AYear: 200811233233518085758
10a. Hod O,Barone V,Scuseria GE. Phys. Rev. BYear: 200877035411
10b. Jiang DE,Sumpter BG,Dai S. J. Chem. PhysYear: 200712613470117430050
10c. Nagai H,Nakano M,Yoneda K,Kishi R,Takahashi H,Shimizu A,Kubo T,Kamada K,Ohta K,Botek E,Champagne B. Chem. Phys. LettYear: 2010489212218
10d. Jiang DE,Dai S. Chem. Phys. LettYear: 20084667275
11. Hachmann J,Dorando JJ,Aviles M,Chan GKL. J. Chem. PhysYear: 200712713430917919026
12. Gidofalvi G,Mazziotti DA. J. Chem. PhysYear: 200812913410819045079
13. Pelzer K,Greenman L,Gidofalvi G,Mazziotti DA. J. Phys. Chem. AYear: 20111155632564021563790
14. Casanova D,Head-Gordon M. Phys. Chem. Chem. PhysYear: 2009119779979019851557
15. Hajgató B,Szieberth D,Geerlings P,De Proft F,Deleuze MS. J. Chem. PhysYear: 200913122432120001050
16. Szalay PG,Müller T,Gidofalvi G,Lischka H,Shepard R. Chem. RevYear: 201211210818122204633
17. Szalay PG,Bartlett RJ. Chem. Phys. LettYear: 1993214481488
18a. Antol I,Eckert-Maksic M,Lischka H,Maksic ZB. Eur. J. Org. ChemYear: 200731733178
18b. Wang EB,Parish CA,Lischka H. J. Chem. PhysYear: 200812904430618681645
19. Takatsuka K,Fueno T,Yamaguchi K. Theor. Chim. ActaYear: 197848175183
20. Staroverov VN,Davidson ER. Chem. Phys. LettYear: 2000330161168
21a. Head-Gordon M. Chem. Phys. LettYear: 2003372508511
21b. Head-Gordon M. Chem. Phys. LettYear: 2003380488489
22a. Watts JD,Urban M,Bartlett RJ. Theor. Chim. ActaYear: 199590341355
22b. Nielsen IMB,Janssen CL. Chem. Phys. LettYear: 1999310568576
23. T. D. Crawford, D. C. Sherrill, F. Plasser, T. Mueller, R. Shepard, H. Lischka, unpublished results.
24. Mulliken RS. J. Chem. PhysYear: 19552318331840
25. Clar E. Polycyclic HydrocarbonsYear: 1964LondonAcademic Press
26. Broene RD,Diederich F. Tetrahedron LettYear: 19913252275230
27a. S. Fujii, T. Enoki, Acc. Chem. Res. 2013, DOI: 10.1021/ar300120y.
27b. Lambert C. Angew. ChemYear: 201112317941796
Angew. Chem. Int. EdYear: 20115017561758
Supplementary material

As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re-organized for online delivery, but are not copy-edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.


[Figure ID: fig01]
Figure 1 

Structures investigated: a) n-acene, b) n-phenacene, c) (ma, nz) periacene, and d) (ma, nz) circumacene.

[Figure ID: fig02]
Figure 2 

Evolution of MR-AQCC NO occupations with the length n of the acene chain.

[Figure ID: fig03]
Figure 3 

Dependence of MR-AQCC NO occupations on the zigzag length in (5a, nz) periacenes.

[Figure ID: fig04]
Figure 4 

Density of unpaired electrons for (5a, 6z) periacene (isovalue 0.005 e). NU=5.0 e, individual atomic values are given next to the respective carbon atoms.

[Figure ID: fig05]
Figure 5 

Dependence of NO occupations on zigzag length n in (5a, nz) circumacenes.

[Figure ID: fig06]
Figure 6 

Dependence of the number of effectively unpaired electrons NU with chain length n.

Article Categories:
  • Communications

Keywords: chemical stability, graphene, multiradical character, multireference methods, quantum chemistry.

Previous Document:  Isolation of Arthrobacter species from the phyllosphere and demonstration of their epiphytic fitness...
Next Document:  Chitosan Electrodeposition for Microrobotic Drug Delivery.