From small fullerenes to the graphene limit: A harmonic force‐field method for fullerenes and a comparison to density functional calculations for Goldberg–Coxeter fullerenes up to C980

We introduce a simple but computationally very efficient harmonic force field, which works for all fullerene structures and includes bond stretching, bending, and torsional motions as implemented into our open‐source code Fullerene. This gives accurate geometries and reasonably accurate vibrational frequencies with root mean square deviations of up to 0.05 Å for bond distances and 45.5 cm−1 for vibrational frequencies compared with more elaborate density functional calculations. The structures obtained were used for density functional calculations of Goldberg–Coxeter fullerenes up to C980. This gives a rather large range of fullerenes making it possible to extrapolate to the graphene limit. Periodic boundary condition calculations using density functional theory (DFT) within the projector augmented wave method gave an energy difference between −8.6 and −8.8 kcal/mol at various levels of DFT for the reaction C60→graphene (per carbon atom) in excellent agreement with the linear extrapolation to the graphene limit (−8.6 kcal/mol at the Perdew–Burke–Ernzerhof level of theory). © 2015 Wiley Periodicals, Inc.

[1]  S. C. O'brien,et al.  C60: Buckminsterfullerene , 1985, Nature.

[2]  Michael Goldberg,et al.  A Class of Multi-Symmetric Polyhedra , 1937 .

[3]  Douglas J. Klein,et al.  Elemental carbon cages , 1988 .

[4]  Andreas Hirsch,et al.  Spherical Aromaticity in Ih Symmetrical Fullerenes: The 2(N+1)2 Rule. , 2000, Angewandte Chemie.

[5]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[6]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[7]  Damir Vukičević,et al.  Atlas of kekulé valence structures of buckminsterfullerene , 2005 .

[8]  Nicola Marzari,et al.  First-principles determination of the structural, vibrational and thermodynamic properties of diamond, graphite, and derivatives , 2004, cond-mat/0412643.

[9]  Dresselhaus,et al.  Force-constant model for the vibrational modes in C60. , 1992, Physical review. B, Condensed matter.

[10]  Broughton,et al.  Intramolecular-force-constant model for C60. , 1992, Physical review. B, Condensed matter.

[11]  Walter Thiel,et al.  MNDO Study of Large Carbon Clusters , 1991 .

[12]  Ante Graovac,et al.  The Mathematics and Topology of Fullerenes , 2011 .

[13]  Marcel Swart,et al.  The role of aromaticity in determining the molecular structure and reactivity of (endohedral metallo)fullerenes. , 2014, Chemical Society reviews.

[14]  Sven J. Cyvin,et al.  Molecular vibrations of footballene , 1988 .

[15]  D. Bowler,et al.  Chemical accuracy for the van der Waals density functional , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[16]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[17]  J. Boettger,et al.  Local-density-functional study of the fullerenes, graphene and graphite , 1996 .

[18]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[19]  Jerzy Cioslowski,et al.  Electronic Structure Calculations on Fullerenes and Their Derivatives , 1995 .

[20]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[21]  D. Bethune,et al.  Bond Lengths in Free Molecules of Buckminsterfullerene, C60, from Gas-Phase Electron Diffraction , 1991, Science.

[22]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[23]  Fernando Flores,et al.  Weak chemical interaction and van der waals forces between graphene layers : A combined density functional and intermolecular perturbation theory approach , 2006 .

[24]  Francesco Mauri,et al.  Impact of the electron-electron correlation on phonon dispersion:Failure of LDA and GGA DFT functionals in graphene and graphite. , 2008, 0808.2285.

[25]  Arnout Ceulemans,et al.  C60 vibrates as a hollow sphere , 1994 .

[26]  Kyuho Lee,et al.  Higher-accuracy van der Waals density functional , 2010, 1003.5255.

[27]  Peter Schwerdtfeger,et al.  The topology of fullerenes , 2014, Wiley interdisciplinary reviews. Computational molecular science.

[28]  Georg Kresse,et al.  Accurate density functional calculations for the phonon dispersion relations of graphite layer and carbon nanotubes , 2003 .

[29]  Brett I. Dunlap,et al.  Efficient quantum-chemical geometry optimization and the structure of large icosahedral fullerenes , 2006 .

[30]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[31]  Harold W. Kroto.  Symmetry, Space, Stars, and C60 (Nobel Lecture) , 1997 .

[32]  Stefan Grimme,et al.  Benchmarking of London Dispersion-Accounting Density Functional Theory Methods on Very Large Molecular Complexes. , 2013, Journal of chemical theory and computation.

[33]  Christoph Rüchardt,et al.  C70 Is More Stable than C60: Experimental Determination of the Heat of Formation of C70 , 1994 .

[34]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[35]  Douglas J. Klein,et al.  Resonance in C60, buckminsterfullerene , 1986 .

[36]  Milan Randic,et al.  Numerical Kekulé Structures of Fullerenes and Partitioning of π-Electrons to Pentagonal and Hexagonal Rings , 2007, J. Chem. Inf. Model..

[37]  Kosmas Prassides,et al.  Crystal structure and bonding of ordered C60 , 1991, Nature.

[38]  James Avery,et al.  Program Fullerene: A software package for constructing and analyzing structures of regular fullerenes , 2013, Journal of computational chemistry.

[39]  Thomas F. George,et al.  Vibrational motions of buckminsterfullerene , 1987 .

[40]  Bun Chan,et al.  Evaluation of the heats of formation of corannulene and C60 by means of high-level theoretical procedures. , 2013, The journal of physical chemistry. A.

[41]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[42]  J. L. Dunn,et al.  A complete nearest-neighbor force field model for C60. , 2004, The Journal of chemical physics.

[43]  P. Fowler,et al.  Molecular graphs, point groups, and fullerenes , 1992 .

[44]  James Avery,et al.  Structure and Properties of the Nonface-Spiral Fullerenes T-C380, D3-C384, D3-C440, and D3-C672 and Their Halma and Leapfrog Transforms , 2014, J. Chem. Inf. Model..

[45]  W. T. Tutte.  How to Draw a Graph , 1963 .

[46]  Patrick W. Fowler,et al.  Complete bond force fields for trivalent and deltahedral cages: Group theory and applications to cubane, closo-dodecaborane, and buckminsterfullerene , 2001 .

[47]  Wojciech Grochala,et al.  Diamond: electronic ground state of carbon at temperatures approaching 0 K. , 2014, Angewandte Chemie.

[48]  James C. Ellenbogen,et al.  Smooth Scaling of Valence Electronic Properties in Fullerenes: From One Carbon Atom, to C60, to Graphene , 2013 .

[49]  Ralf Tonner,et al.  A systematic study of rare gas atoms encapsulated in small fullerenes using dispersion corrected density functional theory , 2015, J. Comput. Chem..

[50]  Clémence Corminboeuf,et al.  Is C60 buckminsterfullerene aromatic? , 2012, Physical chemistry chemical physics : PCCP.

[51]  M. Dion,et al.  van der Waals density functional for general geometries. , 2004, Physical review letters.

[52]  David Moncrieff,et al.  Standard Enthalpies of Formation of Fullerenes and Their Dependence on Structural Motifs , 2000 .