arXiv:0804.0404v1 [cond-mat.mtrl-sci] 2 Apr 2008 h aaytcnb ihrcytlieo iud[1 during [11] liquid 13]. or [12, crystalline (ii) either growing rate- be process the can nucleation a of catalyst the The chirality be upon the to decided while is shown CNT 11], been [10, has step (i) limiting nanotube. fur- process growing the the diffusion into the by The carbon followed (ii) of fragment and incorporation graphitic the ther interior, a on its of carbon across two nucleation of the or by diffusion via surface followed the nanoparticle CVD controlled is (i) step the be processes: initial easily of important This can step and first thermochemistry. CNTs the of is growth surface nanoparticle the achieving the for can optimized extent CNTs? be what of catalyst growth To metal deterministic the process? of growth play choice catalyst CVD metal the the of in composition role chemical What the arise. does questions constants, important Therefore, lattice of [9]. fects variation the of result prop- d a chemical However, as and act [7] 6]. physical to erties [5, different able process show growth also metals CVD are these the Au), in Ag, catalysts and as (Cu, Pt) (Pd, metals transition coinage late nanoparticles inert alloys. heavy metallic metals, their that other and many shown of Ni) been Co, has (Fe, it metals Recently, iron-group to ited 4]. limited [3, yield chiralities to nanotube and of [2] for distributions temperatures allow low to at demonstrated growth been goal. CNT this already reaching has for method believed approach This promising is most nanoparticles the growth metallic be to (CVD) by catalyzed deposition CNTs electronic vapor into of chemical CNTs ob- The incorporating major to devices. the way is the on properties, stacle electronic of thus inability chirality and The predetermined indices, with CNTs [1]. of electronics growth generation performing particu- future in of technologies, lar many of constituent important bn nrisadocptos[] n eaiitcef- relativistic and [8], occupations and energies -band h eopsto fgspaecro rcroson precursors carbon gas-phase of decomposition The lim- been has composition catalyst the time, long a For an become to expected are (CNTs) nanotubes Carbon ntttRmn eRceceNmeiu nPyiu e Mat des Physique Num´erique en Recherche de Romand Institut ASnmes 84.k 10.e 11.h 82.33.Ya 81.15.Gh, 81.07.De, 68.43.Jk, distri chirality numbers: narrow PACS with and particul carbon temperatures in the low metals, at and Coinage CNTs catalyst zigza of metals. through the the than mon favorable of among of more significantly choice is binding the edges the armchair CN by growing of through controlled the terms of be stability in can the analyzed and pathways, is diffusion process their transition late The of nanoparticles metals. on (CNT) nanotubes carbon ffc fMtlEeeti aayi rwho abnNanotub Carbon of Growth Catalytic in Element Metal of Effect sn rtpicpe acltos emdlteceia v chemical the model we calculations, principles first Using nttt fTertclPyis H11 asne Switz Lausanne, CH-1015 Physics, Theoretical of Institute cl oyehiu ´deaed asne(EPFL), F´ed´erale Lausanne Polytechnique de Ecole lgV Yazyev V. Oleg Dtd oebr8 2018) 8, November (Dated: ∗ n lrd Pasquarello Alfredo and etitoso h rcsigtmeaue [2]. temperatures impose processing also the CNTs on restrictions of crys- applications stage. nucleation the technological the Certain maintaining during nanoparticles for of state necessary tem- of talline low growth thus Sufficiently the are for [14]. peratures chiralities templates specific as with act nanoparticles CNTs can metallic shape solid definite only of but process growth the 2] h ufcswr oee ytrelyrsaswith slabs three-layer by nanoparticles p modeled fcc-metal were surfaces of The surfaces [20], [21]. facets We the (111) dominate [17]. energy which nanoparticles lowest the metallic considered in smallest even primarily the facets, of for models case good the be to shown Infinite were surfaces respectively. slab models, two-dimensional bulk periodic three-dimensional and of means by modeled were low distributions. at chirality narrow CNTs with Cu, growing and particular for temperatures in opportunities metals, chiralities. attractive coinage definite grow- offer that for the show preference of results the edges Our infer armchair we and CNT, zigzag ing From of pathways. stability diffusion bulk the and in con- surface as energies both activation well sidering these calculate For as then we compared. facets species, are the carbon nanoparticles energies on the of dimers binding interior and the the atoms metal, pro- carbon considered growth of each nucleation CVD and For diffusion The the and steps. modeling by Pt) metals. analyzed Pd, Au) is cess In (Ni, Ag, transition scheme. (Cu, the late total-energy coinage consider in a we using nanoparticles study, CNTs our catalytic of the growth CVD of composition ical rcino h ufc rao c-ea uffclusters. Wulff fcc-metal of area surface largest (100) the second extended of the This for fraction to approximation thickness. corresponding good same facets, a nanoparticle the also of is surface model stepped (311) considered modeled a also we by which We step-edges, eV. surface 0.1 microfaceted (100) within energies binding tive (3 h ufc n h neiro ealcnanoparticles metallic of interior the and surface The chem- the of role the investigate we work, this In × e iu IRA,C-05Luan,Switzerland Lausanne, CH-1015 ´eriaux (IRRMA), )ui el,wihgaatecnegneo rela- of convergence guarantee which cells, unit 3) rC,aefudt ao V growth CVD favor to found are Cu, ar N,P,P)adciae(u g Au) Ag, (Cu, coinage and Pt) Pd, (Ni, .W n httedffso pathways diffusion the that find We T. ns u h eaiesaiiyvaries stability relative the but ones, g butions. prdpsto CD rwhof growth (CVD) deposition apor rcro.Bnigo h CNT the of Binding precursor. -addaoi abnspecies, carbon diatomic and o- radand erland es 2

where E and EM are the total energies of the systems with and without carbon, respectively, and nC corre- sponds to the number of carbon atoms. For reference, the chemical potential of the free carbon atom in the spin-triplet ground state is 7.7 eV. We considered the binding and the diffusion of both carbon atoms and dimers, which are initially produced upon the dehydrogenation of common carbon sources like CH4, C2H4, C2H2, and C2H5OH. For all considered po- sitions [Fig. 1(a)] and metals, the chemical potentials of monoatomic and diatomic carbon are positive, i.e. indi- cating lower stability than for graphene, and are generally higher for the inert coinage metals [Fig. 1(b)]. On coinage metals, the surface dimers are considerably more stable than all forms of monoatomic carbon. This implies that the diffusion process only involves surface dimers, when the carbon source corresponds to a diatomic carbon pre- FIG. 1: (Color online) (a) Schematic representation of studied cursor such as C2H2. At variance, atomic carbon has carbon positions at metal surfaces. (b) Chemical potentials the lowest chemical potential for late transition metals, of monoatomic and diatomic carbon in different positions ex- as a subsurface interstitial in Ni and Pd and as a surface pressed with respect to the chemical potential of carbon in graphene. (c) Surface diffusion mechanisms for the carbon adatom at the Pt surface. adatom (1) and the carbon dimer (2 – concerted mode; 3 Adatoms on a (111) surface are threefold coordinated – atom-by-atom mode) on a (111) fcc metal surface and (d) and occupy two different positions, ccp and hcp [17], respective diffusion barriers for different metals. with chemical potentials differing by less than 0.1 eV. Only the lowest adatom chemical potentials are shown in Fig. 1(b). The highest adatom chemical potential was In the two-dimensional models, the periodically repeated obtained for Ag, in agreement with the calculations of ˚ slabs were separated by at least 10 A to avoid spurious Abild-Pedersen et al. [27] This low stability results from × × interactions. The bulk metals were modeled with 3 3 3 the energy position of the d-band, which is lowest for this supercells. metal [8]. At surface steps, the carbon adatom assumes a The total energy calculations were performed within a semi-octahedral five-fold coordination showing the same selfconsistent density functional theory framework, as im- local neighborhood as at an adsorption site on the (100) plemented in the PWSCF plane-wave pseudopotential code surface. For Ni and Cu the adsorption site at the step is of the Quantum-ESPRESSO distribution [22]. The Perdew- preferred over that on a (111) terrace, as found for Ni in Burke-Ernzerhof exchange-correlation density functional a previous study [18]. However, the converse is true for [23] was used. The ultrasoft pseudopotentials [24] em- all heavier metals showing larger lattice constants. ployed in the present study treat the d-electrons of both Carbon atoms in subsurface and bulk positions occupy transition and coinage metals as valence electrons. The either tetrahedral or octahedral voids, in which they show one-electron valence wave functions and the electron den- coordination numbers of 4 and 6, respectively. Subsur- sity were described by plane-wave basis sets with kinetic face carbon atoms located in voids closest to the sur- energy cutoffs of 25 Ry and 300 Ry, respectively [25]. face are more stable than surface carbon adatoms for all The Fermi-Dirac function (kBT =0.01 Ry) was used for metals except Pt. The preferred occupation sites are the the population of the one-electron states in metallic sys- tetrahedral voids in Au and Pt characterized by large lat- × × k tems. For the slab models, a 4 4 1 point sampling tice constants, and octahedral voids in the lighter metals. was performed. The bulk supercells were sampled with Only the lowest-energy values are given in Fig. 1(b). Bulk × × k 4 4 4 point meshes. The atomic positions of the bot- interstitials, i.e. carbon atoms located deep inside the tom layer were kept fixed in ideal bulk positions during crystalline nanoparticle, have considerably higher chem- the relaxation in the slab models. The nudged elastic ical potentials than subsurface interstitials due to an en- band method with the climbing-image scheme [26] was hancement of the elastic response energy [19, 28]. Bulk used to find the transition states of the diffusion paths. interstitials preferentially reside in octahedral voids for Binding of monoatomic and diatomic carbon. The all metals studied. Subsurface and bulk dimers are found chemical potential of different forms of carbon expressed to be highly unstable and are not further discussed here. with respect to the total-energy per atom of graphene µG Activation barriers of diffusion. The relevant diffusion is defined as channels of monoatomic and diatomic carbon [Fig. 1(c)] E − EM on metallic nanoparticles are addressed through the cal- µ = − µG, (1) nC culation of their activation barriers Eact [Fig. 1(d)]. For 3 monoatomic carbon, we considered surface, subsurface and bulk diffusion modes. In case of diatomic carbon, only the surface diffusion mode was investigated. The re- ported barrier is the total-energy difference between the transition state and the lowest of the connected minima. The adatom diffusion barrier is lower for the coinage metals with respect to the late-transition metal belong- ing to the same row. Moving down the periodic ta- ble results in higher diffusion barriers. For the adatom diffusion on the Ni(111) surface, we found a barrier of 0.39 eV, in agreement with previous theoretical estimates [10, 15, 17, 19]. This value is also in fair agreement with the activation energy of the experimental CNT growth FIG. 2: (Color online) (a) Schematic representation of the process, suggesting that surface diffusion is the rate- binding of ZZ (left) and AC (right) carbon nanotube to a flat limiting step [10, 11]. The adatom diffusion on Cu and (111) surface of the catalyst nanoparticle. The edge fragments Ag surfaces is characterized by surprisingly low barri- used to model the interaction between the nanotube and the nanoparticle surface are shown in red. (b) Chemical poten- ers of 0.07 eV and 0.20 eV, respectively. This implies tials of edge atoms of AC and ZZ fragments bound to the that a significant acceleration of the diffusion step can metal surfaces and steps. The chemical potential of carbon be achieved for these metals when a monoatomic carbon atoms at free AC and ZZ edges are 2.13 eV and 2.83 eV, re- precursor (e.g. CH4) is used. spectively. (c) Thermodynamic chirality preference ratios at ◦ Subsurface and bulk diffusion occur via hopping of in- a synthesis temperature of 900 C for various metal catalysts. terstitial carbon atoms between adjacent tetrahedral and octahedral voids. The diffusion barriers for this channel such an investigation is computationally prohibitive be- is generally higher than for adatom diffusion. However, cause of the large configuration space involved [12, 13]. the activation energy for subsurface diffusion decreases To address this issue, we here focused on the binding with increasing lattice constant, following an opposite of minimal structural units of armchair (AC) and zigzag trend with respect to the adatom diffusion and leading (ZZ) edges to either flat or stepped surfaces. We modeled to competition between the two diffusion modes for the AC and ZZ edges by 1,2-didehydrobenzene (C H ) and heaviest metals. Activation barriers of subsurface dif- 6 4 dehydrobenzene (C H ) molecules having 2 and 1 edge fusion are always lower than those of the corresponding 6 5 atoms, respectively, as illustrated in Fig. 2(a). bulk diffusion due to the reduced elastic response in prox- The binding energies are quantified by the chemical imity of the nanoparticle surface. However, in the partic- potential per edge atom µ defined as ular case of Au, activation energies for surface, subsur- edge face, and bulk diffusion fall in a narrow range between − − free E Efrag EM 0.52 and 0.64 eV, implying that monoatomic carbon in µedge = µedge + , (2) nC gold can diffuse uniformly across the nanoparticle already free at low temperatures. where µedge is the chemical potential per edge atom for an When the carbon precursor gives rise to diatomic unbound edge, Efrag the total energy of an isolated model carbon, the dimer diffusion channel becomes relevant. fragment, and nC the number of edge atoms. We esti- free We considered two different mechanisms for the surface mated µedge through the calculations of ideal graphene free dimer diffusion, the concerted and the atom-by-atom and suitable graphene nanoribbons. For µedge of AC and modes, as illustrated pictorially in Fig. 1(c). For the ZZ edges, we found values of 2.13 eV and 2.83 eV, in metals studied, the concerted mode shows activation bar- good agreement with previously reported results [16]. riers exhibiting small variations. Calculated values range The calculated chemical potentials µedge correspond- between 0.43 eV for Ag and 0.70 eV for Pt, and are com- ing to optimally placed fragments bound to either facets parable to the diffusion barrier of the carbon adatom on or step-edges are given in Fig. 2(b) for the various metals. Ni(111). The atom-by-atom mode yields higher diffusion These results show that all studied metals share common barriers with energy differences ranging from almost zero qualitative features: (i) The chemical potententials µedge free for Pt to 0.62 eV for Ni. In Fig. 1(d), the activation ener- are always lower than µedge, implying that the binding gies reported for dimer diffusion correspond to the lowest of the graphitic fragment to the metallic nanoparticles is energy mode, i.e. to the concerted mode. energetically favored. Binding to step-edges is system- Binding of CNT to metal nanoparticle. To understand atically preferred over binding to flat terraces, in accord the relation between the chirality of the growing CNT with experimental indications [14, 29]. (ii) The chemical and the chemical composition of the catalyst, a compre- potentials µedge in Fig. 2(b) and those of mono-atomic hensive study of the interaction between graphitic frag- and di-atomic species in Fig. 1(b) exhibit similar behav- ments and metallic nanoparticles is required. However, ior as the metal is varied, consistent with the properties 4 of the metal d-band [8]. However, the former chemical catalytic nanoparticles of coinage metals. For these met- potentials are always lower for corresponding metals, in- als, the stability and the diffusion barriers of diatomic dicating that the aggregation of bound C atoms in the carbon allow one to restrict the diffusion pathways to form of graphitic fragments is energetically favored. (iii) the nanoparticle surface by choosing an appropriate gas- The chemical potentials µedge are always positive, i.e. phase carbon source. Furthermore, particularly low diffu- the C atoms are less stable than in the reference ideal sion barriers are found for carbon adatoms in the case of graphene, suggesting that there exists a driving force to- Cu and Ag, suggesting the possibility of realizing growth wards extending the size of graphitic aggregates. The at low growth temperatures. The Cu catalyst is also combination of these three features points to the exis- found to be the most promising candidate for growing tence of a negative chemical potential gradient for the carbon nanotubes of definite chirality. nucleation and further growth of CNTs from a carbon feedstock. This supports the catalytic activity of a wide We acknowledge useful discussions with L. Forr´o, ˇ ˇ range of metals, including the inert coinage metals. K. Lee, A. Magrez, and Z. Sljivanˇcanin. We used com- putational resources at CSCS. The chemical potentials µedge in Fig. 2(b) also pro- vide insight into the ability of different metal catalysts of growing CNTs with definite chiralities. For all considered metals, AC edges are found to be more stable than ZZ ∗ edges. In order to provide a more quantitative measure oleg.yazyev@epfl.ch of this preference at a growth temperature T , we define [1] R. H. Baughman, A. A. Zakhidov, and W. A. de Heer, the ratio Science 297, 787 (2002). AC − ZZ [2] M. Cantoro et al., Nano Lett. 6, 1107 (2006). T µ µ R = exp − edge edge , (3) [3] S. M. Bachilo et al., J. Am. Chem. Soc. 125, 11186 " kT # (2003). [4] Y. Miyauchi et al., Chem. Phys. Lett. 387, 198 (2004). assuming that the CNT nucleation occurs under thermo- [5] D. Takagi et al., Nano Lett. 6, 2642 (2006). dynamic control. For a graphitic CNT nucleus, the chem- [6] W. Zhou et al., Nano Lett. 6, 2987 (2006). ical potential of the edge carbon atoms is higher than [7] H. Okamoto, Desk Handbook. Phase Diagrams for Binary that of the other carbon atoms, leading to the growth of Alloys (ASM International, Materials Park, 2000). CNT with minimal edge perimeters. Under this assump- [8] B. Hammer, Y. Morikawa, and J. K. Nørskov, Phys. Rev. 76 tion, the ratio RT characterizes the distribution of CNT Lett. , 2141 (1996). [9] H. H¨akkinen, M. Moseler, and U. Landman, Phys. Rev. chiralities indicating the preference of growing AC-like Lett. 89, 033401 (2002). or ZZ-like CNTs. Ultimate selectivity is achieved for ZZ 95 T T [10] S. Hofmann et al., Phys. Rev. Lett. , 036101 (2005). and AC CNTs for R =0 and R →∞, respectively. [11] K. Bartsch et al., J. Appl. Phys. 97, 114301 (2005). T Figure 2(c) shows R for the different metal catalysts [12] S. Reich, L. Li, and J. Robertson, Phys. Rev. B 72, at a typical CVD synthesis temperature of 900◦C. All 165423 (2005). metals show a strong preference for growing AC rather [13] S. Reich, L. Li, and J. Robertson, Chem. Phys. Lett. 421, than ZZ nanotubes, in agreement with experimental re- 469 (2006). [14] H. Zhu et al., Small 1, 1180 (2005). ports for iron-group transition metals [3]. In Fig. 2(c), 427 900◦C [15] S. Helveg et al., Nature , 426 (2004). R ranges between 5 for Ni and 130 for Cu. The [16] X. Fan et al., Phys. Rev. Lett. 90, 145501 (2003). case of Cu is particularly noteworthy for several reasons. [17] Q.-M. Zhang et al., Phys. Rev. B 69, 205413 (2004). T First, its high value of R implies that the CNTs grown [18] H. S. Bengaard et al., J. Catal. 209, 365 (2002). on Cu nanoparticles are highly enriched by metallic AC [19] F. Abild-Pedersen et al., Phys. Rev. B 73, 115419 (2006). species, in accordance with preliminary experimental re- [20] L. Vitos et al., Surf. Sci. 411, 186 (1998). 77 ports [6]. Second, Cu features the lowest carbon diffu- [21] F. Baletto and R. Ferrando, Rev. Mod. Phys. , 371 (2005). sion barriers suggesting that CVD synthesis could take [22] S. Baroni et al., http://www.quantum-espresso.org/. place at much lower temperatures. In these conditions, T [23] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. the chirality preference ratio R would be further en- Lett. 77, 3865 (1996). hanced. Third, for Cu, the chemical potentials in Fig. [24] D. Vanderbilt, Phys. Rev. B 41, R7892 (1990). 2(b) show that graphitic fragments indifferently bind to [25] A. Pasquarello et al., Phys. Rev. Lett. 69, 1982 (1992); flat or stepped surfaces, suggesting that step-edge forma- K. Laasonen et al., Phys. Rev. B 47, 10142 (1993). tion on Cu nanoparticles is not likely to occur during the [26] G. Mills, H. J´onsson, and G. K. Schenter, Surf. Sci. 324, 305 (1995); G. Henkelman, B. P. Uberuaga, and nucleation of carbon nanotubes [14, 29]. Preserving the H. J´onsson, J. Chem. Phys. 113, 9901 (2000). shape of the nanoparticle is a prerequisite for the realiza- [27] F. Abild-Pedersen et al., Phys. Rev. Lett. 99, 016105 tion of chirality-specific CNT growth on predetermined (2007). catalytic templates. [28] M. Todorova et al., Phys. Rev. Lett. 89, 096103 (2002). In conclusion, our study of the diffusion and nucleation [29] J. A. Rodriguez-Manzo et al., Nat. Nanotechnol. 2, 307 stages of CVD growth of CNTs highlights the potential of (2007).