Ab initio calculations using density functional theory (DFT) have been performed in order to explore structure and energy gap opening of graphene with bridged-trivacancy and single adsorbed with Fe atom. Compared to the previous reconstructed trivacancy adsorbed with Fe atom with the energy gap of 0.10 eV, one interesting structure for the Fe-doped bridged-trivacancy complex has been identified, with one Fe atom above the graphene plane, and possesses energy gap with the value of 0.32 eV in the bridged circumstance. The band gap can be explained by the decrease of the free electrons. These results provide insights to engineer graphene’s properties through defect addition and manipulation for industrial semiconductor applications such as the photocatalytic technology and graphene based electronics.
1. Introduction
Graphene is among the most potential utilized materials currently studied [1–3]. It exhibits extraordinary properties in electrology [4, 5], magnetism [6, 7], and dynamics [8], proved by various theoretical and experimental studies. However, graphene is a semimetal, which greatly limits the application in the graphene-based electronic field [9–11]. Consequently, in the recent decades, many approaches, such as by employing a staggered AB sublattice potential [12], strain effect [13], intrinsic spin-orbit coupling [14] in single-layer graphene, functionalization with chemical groups [15–17], or applying a perpendicular electric field in bilayer graphene [18], have been explored to open the gap of graphene for designing semiconductor device.
Chan et al. have studied the structures and magnetism about pure graphene adsorbed with many elements such as Li, Na, K, Ca, and Al [19]. However, precise and stable structures of the transition metal (TM) adsorbed on graphene are hardly obtained because TM adsorbents can move freely on the surface of pristine graphene, which limits the use of surface adatom as robust and reliable dopants for nanoelectronic and magnetic applications [20]. Recently, it is forecasted that attaching metal atoms to defected graphene is much firmer than to pristine graphene by theoretical studies [21–25]. Kim et al. have used X-ray magnetic circular dichroism (XMCD) to study the electronic and magnetic properties of transition metal (TM), for example, Fe, Co, and Ni clusters on monolayer graphene [26]. Also a very high magnetic anisotropy energy (MAE) has been found in single-vacancy graphene adsorbed with Ir dimer [27]. Dynamics of single Fe adatom on graphene vacancies have been studied using focused electron beam irradiation [28]. Vacancy-assisted doping method is one available technique for introducing dopants into graphene.
As a matter of fact, iron in graphene nanostructures are very common in experiment, such as Fe pairs doped in graphene [20], single Fe atoms in graphene vacancies [29, 30], and Fe3O4/graphene nanocomposites [31–34]. In previous work, band gap of graphene could be opening via Fe atom adsorption [11]. Recently, an interesting structure of graphene trivacancy, which was achieved by using an electron beam to bombard the graphene sheet, has been commonly observed by Robertson et al. [35]. The new structure has a bridging atom (thereafter called bridged-trivacancy), significantly different from the reconstructed trivacancy (r-trivacancy) observed by Wang et al. [36]. Consequently, a combination of Fe adsorption and b-trivacancy might contribute to extend the application of graphene by opening the band gap.
This work presents a detailed study of the structure and electronic properties of graphene with r-trivacancy and b-trivacancies trapping Fe atom, using the density functional theory (DFT). Our studies may verify the feasibility of improving graphene electrical properties by controlling TM atom doping and vacancies.
2. Calculation Details
First-principles calculations based on the spin-polarized density functional theory (DFT) were performed using the projector augmented wave (PAW) [37], as implemented in the Vienna Ab initio Simulation Package (VASP) program [38]. The generalized gradient approximation in the Perdew-Burke-Ernzerhof (PBE) [39] parameterization was employed for the exchange-correlation functional. The plane-wave cutoff energy was set to 500 eV. The geometry optimization and total energy calculations were performed until the residual force was within 0.01 eV/Å and the total energy converged to 10−5 eV. The k-point mesh for sampling the Brillouin zone was generated including the gamma point. A 4 × 4 × 1 k-point mesh was used for geometry optimization, and a 6 × 6 × 1 k-point mesh was used for total energy calculations.
We modeled the defect systems using supercells large enough (5 × 5 and 6 × 6 supercells in the lateral direction with a 20 Å thick vacuum in the vertical direction) to minimize boundary effects on the energetics of the systems. Graphene’s lattice constant was calculated to be 2.470 Å, which is similar to the calculated value (2.468 Å) by Dai et al. and the experimental value of 2.46 Å [21].
3. Structures for System
Three kinds of trivacancies in graphene have been discovered in previous research, the r-trivacancy observed by Wang in 2012 [35] and the b-trivacancy observed by Robertson in 2014 [36], and the optimized equilibrium configurations for the r-trivacancy and b-trivacancy in our simulations using DFT are shown in Figure 1. The b-trivacancy is not just stabilized by simple bond reconstructions between undercoordinated carbon atoms, as exhibited by r-trivacancy. Reconstructions consist of undercoordinated bridging carbon atoms spanning the vacancy to saturate edge atoms in b-trivacancy, which matches well with the previous work [36].
(a) Optimized 6 × 6 graphene supercell with r-trivacancy. (b) Optimized 6 × 6 graphene supercell with b-trivacancy.
All the systems with Fe adsorbate on perfect graphene and defective graphene (containing a r-trivacancy or b-trivacancy) are studied. Three possible adsorption sites of Fe on the perfect graphene are labeled as follows: hollow (H), top (T), and bridge (B). And the complexes of perfect graphene adsorbed with Fe which locates in the three possible sites are indicated as Fe@G(H), Fe@G(T), and Fe@G(B), respectively. For adsorption on defect-free graphene, Fe@G(H) (on the hollow site) is the most stable complex, which is in good agreement with the previous theoretical results [5, 19]. In the case of the adsorption of Fe on the defective graphene, Fe is positioned at different hollow, top, and bridge sites in the trivacancy region as the initial positions and then fully relaxed. The configurations of the adsorption of Fe on defective graphene are shown in Figure 2, in both top view and side view. For convenience, the complex of Fe adsorbed on the graphene with a r-trivacancy is labeled as Fe@RTV, and the one with b-trivacancy as Fe@BTV. From the side views, it is obvious that Fe atom is above the plane of graphene with b-trivacancy and the adatom-graphene distance d (Å) is 1.38 Å.
The optimized configurations of Fe adsorbed on the 6 × 6 graphene supercell with a trivacancy in both top view and side view. (a) Fe on the graphene with r-trivacancy (Fe@RTV) and (b) Fe on the graphene with b-trivacancy (Fe@BTV).
The formation energy Ef of a trivacancy in graphene, indicated by Ef, is defined as(1)Ef=Es-N-nNEg,where n=3 is the number of vacant C atoms in the system, and N is the total number of C atoms in supercells of the perfect graphene. Es and Eg are the energies of the graphene with r-trivacancy or b-trivacancy and the 5 × 5 or 6 × 6 graphene supercells, respectively. The adsorption energy of Fe on the perfect or defective graphene is defined as (2)Eads=Es-Eg-EFe,where Es, Eg, and EFe are the energies for the optimized systems of Fe-graphene with or without trivacancy, the 5 × 5 or 6 × 6 graphene supercells with or without trivacancy, and the isolated Fe atom, respectively.
Table 1 summarizes the energetic and structural properties of the 5 × 5 and 6 × 6 supercells of graphene with trivacancy adsorbed with Fe. For comparison, the formation energy of the trivacancy in graphene and the adsorption energy of Fe on the defect-free graphene are also included in Table 1. The adsorption energy for Fe@G(H) complex is the lowest, indicating that the stable position is the hollow site. The adsorption energy of Fe adsorbed on defective graphene with a trivacancy (Fe@TV) is much lower than that of adsorption on perfect graphene, indicating Fe@TV is more stable, which conforms to the previous studies [22]. It means that the introduction of a trivacancy on graphene will improve the adsorption of Fe atom. The formation energy of r-trivacancy is 11.68 eV, consistent with the result in [21], lower than that of b-trivacancy (14.86 eV), indicating that the b-trivacancy is a metastable structure. It is in good agreement with previous studies [35]. From Table 1, we also notice that graphene supercell with b-trivacancy adsorbed with Fe (Fe@BTV) has a lower adsorption energy in comparison with r-trivacancy adsorbed with Fe (Fe@RTV) in the same concentration, which relates with a better incorporation of the adsorbate into the graphene framework. It is worth noting that the band gaps of Fe@BTV are 0.32 eV and 0.25 eV for 5 × 5 and 6 × 6 graphene supercells, respectively. By contrast, the band gaps of Fe@RTV are 0.13 eV and 0.03 eV for 5 × 5 and 6 × 6 supercells, respectively. The system of Fe adsorbed on the perfect graphene has the maximum band gaps of 0.08 eV, which is similar to the result (100.8 meV) in [11]. Fe@BTV system may have a bearing on important applications such as graphene-based electronic device.
Energetic and structural properties of Fe on the defect-free graphene and defective graphene with r-trivacancy and b-trivacancy. The properties listed include the size of graphene supercells s, types of the Fe-graphene complex t, the Fe-graphene distance d (Å), formation energy Ef (eV) of trivacancy in graphene, adsorption energies Eads (eV), bulk band gap Δ (eV), and the magnetization of the adatom μag (μB). For comparison, the magnetization of the graphene with a trivacancy is included.
s
t
d (Å)
Ef (eV)
Eads (eV)
Δ (eV)
μag
(μB)
5 × 5
Fe@RTV
0.00
—
−6.35
0.13
2.00
Fe@BTV
1.24
—
−6.90
0.32
0.00
Fe@G(H)
1.65
—
−1.04
0.08
2.00
Fe@G(B)
2.23
—
−0.28
0.00
4.10
Fe@G(T)
2.06
—
−0.27
0.00
4.10
r-trivacancy
—
11.68
0.00
1.05
b-trivacancy
—
14.86
0.00
1.66
6 × 6
Fe@RTV
0.00
—
−6.61
0.03
2.00
Fe@BTV
1.38
—
−6.96
0.25
0.00
Fe@G(H)
1.50
—
−1.13
0.07
2.00
Fe@G(B)
2.20
—
−0.29
0.00
4.10
Fe@G(T)
2.12
—
−0.28
0.00
4.10
r-trivacancy
—
11.13
0.00
1.05
b-trivacancy
—
14.27
0.00
1.66
4. DOS for Fe@BTV
The graphene with b-trivacancy and r-trivacancy is metallic, as the existence of trivacancy in graphene will increase the number of free electrons in the system. The partial densities of states (PDOS) for spin-up and spin-down electrons of Fe atom and adjacent C atoms and the total density of states (TDOS) in 5 × 5 and 6 × 6 graphene are calculated as shown in Figure 3. The indices of C atoms are labeled in Figure 2(b) and C25 labels the bridge C atom. Obviously, for b-trivacancy, the spin-up and spin-down DOS of C25, C23, C31, Fe, and Fe@BTV is mirror symmetry, indicating no magnetic moment. The PDOS of adjacent C atoms match well with Fe 3d orbitals. In 5 × 5 graphene all the PDOS locate at −2~−0.75 eV, −6~0 eV, and 0.32~0.7 eV. Similar conclusions can be got in 6 × 6 graphene. In particular, the PDOS of C252p orbitals match so consistently with Fe 3d orbitals. It suggests that Fe atom has formed covalent bonds with adjacent C atoms. Stronger Fe-C interaction was observed between Fe-C25 than that between Fe-C23 and Fe-C31. The magnetization of Fe in b-trivacancy is 0μB while the Fe atom in r-trivacancy has magnetization of 2μB (Table 1). This is consistent with the fact that there are enough ten electrons in the 3d orbitals of Fe in b-trivacancy after two 4s electrons of Fe transferred into 3d orbitals and sharing two 2p electrons from neighbouring C atoms.
The partial densities of states (PDOS) for spin-up and spin-down electrons for Fe atom and C atoms vicinity of adsorbed Fe and the total density of states (TDOS) in 5 × 5 and 6 × 6 graphene: (a) 5 × 5 graphene, (b) 6 × 6 graphene.
Adsorbate opens the band gap of graphene with b-trivacancy. Fe adatom has formed covalent bonds with adjacent C atoms in the b-trivacancy region, and the number of free electrons is reduced, resulting in the opening of the energy gap.
The band structures for Fe@BTV are shown in Figure 4. It could be found that an energy gap is opened by trapping Fe to the b-trivacancy and the energy gap is larger in Figure 4(a) than that in Figure 4(b). The value of energy gap is 0.32 and 0.25 for (a) and (b), respectively.
Band structures of Fe@BTV complex. (a) Optimized 5 × 5 graphene supercell with b-trivacancy trapping Fe. (b) Optimized 6 × 6 graphene supercell with b-trivacancy trapping Fe.
5. Conclusions
In conclusion, density functional theory calculations have been performed to investigate the interaction of Fe adatom with bridged-trivacancy in graphene. The graphene with b-trivacancy trapping Fe atom has one stable structure, namely, Fe@BTV. The adsorption of Fe atom can open a finite energy gap in both r-trivacancy and b-trivacancy. The band gap of Fe@BTV is 0.32 eV, larger than the value of Fe@RTV (0.1 eV), which is close to the value for developing graphene-based electronics. Fe adatom has formed covalent bonds with adjacent C atoms in the b-trivacancy region, and the number of free electrons is reduced, resulting in the opening of the energy gap. The findings might provide a promising scheme for gap opening of graphene for industrial semiconductor applications and the photocatalytic technology.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors thank R. Yu for discussions. This work was supported by National Basic Research Program of China (2011CB606406) and the Fundamental Research Funds for the Central Universities (TP-A3: 06108170). This work used the resources of Shanghai Supercomputer Center.
DasS.KimM.LeeJ.ChoiW.Synthesis, properties, and applications of 2-D materials: a comprehensive review201439423125210.1080/10408436.2013.836075GeimA. K.NovoselovK. S.The rise of graphene20076318319110.1038/nmat18492-s2.0-33847690144YazyevO. V.HelmL.Defect-induced magnetism in graphene2007751212540810.1103/physrevb.75.1254082-s2.0-33847768178Castro NetoA. H.GuineaF.PeresN. M. R.NovoselovK. S.GeimA. K.The electronic properties of graphene200981, article 10910.1103/revmodphys.81.1092-s2.0-59949098337MaoY.YuanJ.ZhongJ.Density functional calculation of transition metal adatom adsorption on graphene2008201111520910.1088/0953-8984/20/11/1152092-s2.0-41849124219ZhangY.TalapatraS.KarS.VajtaiR.NayakS. K.AjayanP. M.First-principles study of defect-induced magnetism in carbon20079910720110.1103/physrevlett.99.1072012-s2.0-34548573442TopsakalM.AktürkE.SevinçliH.CiraciS.First-principles approach to monitoring the band gap and magnetic state of a graphene nanoribbon via its vacancies20087823543510.1103/physrevb.78.2354352-s2.0-59049096187SusiT.KotakoskiJ.KepaptsoglouD.ManglerC.LovejoyT. C.KrivanekO. L.ZanR.BangertU.AyalaP.MeyerJ. C.RamasseQ.Silicon–carbon bond inversions driven by 60-keV electrons in graphene20141131111550110.1103/physrevlett.113.115501NovoselovK. S.JiangZ.ZhangY.MorozovS. V.StormerH. L.ZeitlerU.MaanJ. C.BoebingerG. S.KimP.GeimA. K.Room-temperature quantum hall effect in graphene20073155817137910.1126/science.11372012-s2.0-33947176113ZhangY.TanY.-W.StormerH. L.KimP.Experimental observation of the quantum Hall effect and Berry's phase in graphene2005438706520120410.1038/nature042352-s2.0-27744475163DingJ.QiaoZ.FengW.YaoY.NiuQ.Engineering quantum anomalous/valley Hall states in graphene via metal-atom adsorption: an ab-initio study20118419544410.1103/physrevb.84.1954442-s2.0-82655179060XiaoD.YaoW.NiuQ.Valley-contrasting physics in graphene: magnetic moment and topological transport20079923680910.1103/physrevlett.99.236809GuineaF.KatsnelsonM. I.GeimA. K.Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering20106303310.1038/nphys1420KaneC. L.MeleE. J.Z2 topological order and the quantum spin hall effect2005951414680210.1103/PhysRevLett.95.146802BoukhvalovD. W.KatsnelsonM. I.Tuning the gap in bilayer graphene using chemical functionalization: density functional calculations200878808541310.1103/physrevb.78.0854132-s2.0-49649123498HammouriM.JhaS. K.VasilievI.First-principles study of graphene and carbon nanotubes functionalized with benzyne201511932187191872810.1021/acs.jpcc.5b04065YanA. J.ChouM. Y.Oxidation functional groups on graphene: structural and electronic properties20108212540310.1103/physrevb.82.125403OostingaJ. B.HeerscheH. B.LiuX.MorpurgoA. F.VandersypenL. M. K.Gate-induced insulating state in bilayer graphene devices2008715115710.1038/nmat2082ChanK. T.NeatonJ. B.CohenM. L.First-principles study of metal adatom adsorption on graphene2008772323543010.1103/physrevb.77.2354302-s2.0-45749127809HeZ.HeK.RobertsonA. W.KirklandA. I.KimD.IhmJ.YoonE.LeeG.-D.WarnerJ. H.Atomic structure and dynamics of metal dopant pairs in graphene20141473766377210.1021/nl500682j2-s2.0-84903976362DaiX. Q.ZhaoJ. H.XieM. H.TangY. N.LiY. H.ZhaoB.First-principle study of magnetism induced by vacancies in graphene201180334334910.1140/epjb/e2011-10955-xKrasheninnikovA. V.NieminenR. M.Attractive interaction between transition-metal atom impurities and vacancies in graphene: a first-principles study20111293–562563010.1007/s00214-011-0910-3CretuO.KrasheninnikovA. V.Rodríguez-ManzoJ. A.SunL.NieminenR. M.BanhartF.Migration and localization of metal atoms on strained graphene201010519610210.1103/physrevlett.105.196102HaldarS.PujariB. S.BhandaryS.CossuF.ErikssonO.KanhereD. G.SanyalB.Fen(n=1−6) clusters chemisorbed on vacancy defects in graphene: stability, spin-dipole moment, and magnetic anisotropy2014892020541110.1103/physrevb.89.205411KrasheninnikovA. V.LehtinenP. O.FosterA. S.PyykköP.NieminenR. M.Embedding transition-metal atoms in graphene: structure, bonding, and magnetism200910212412680710.1103/physrevlett.102.126807KimG.JhiS.LimS.ParkN.Effect of vacancy defects in graphene on metal anchoring and hydrogen adsorption20099417310210.1063/1.3126450EelboT.WaśniowskaM.ThakurP.GyamfiM.SachsB.WehlingT. O.FortiS.StarkeU.TiegC.LichtensteinA. I.WiesendangerR.Adatoms and clusters of 3d transition metals on graphene: electronic and magnetic configurations20131101313680410.1103/physrevlett.110.1368042-s2.0-84875691483HanY.GeG.-X.WanJ.-G.ZhaoJ.-J.SongF.-Q.WangG.-H.Predicted giant magnetic anisotropy energy of highly stable Ir dimer on single-vacancy graphene20138715540810.1103/physrevb.87.155408RobertsonA. W.MontanariB.HeK.KimJ.AllenC. S.WuY. A.OlivierJ.NeethlingJ.HarrisonN.KirklandA. I.WarnerJ. H.Dynamics of single Fe atoms in graphene vacancies20131341468147510.1021/nl304495v2-s2.0-84876038898KrasheninnikovA. V.LehtinenP. O.FosterA. S.PyykköP.NieminenR. M.Embedding transition-metal atoms in graphene: structure, bonding, and magnetism20091021212680710.1103/physrevlett.102.1268072-s2.0-64049105752Rodríguez-ManzoJ. A.CretuO.BanhartF.Trapping of metal atoms in vacancies of carbon nanotubes and graphene2010463422342810.1021/nn100356q2-s2.0-77955854646ZhouK.ZhuY.YangX.LiC.One-pot preparation of graphene/Fe3O4 composites by a solvothermal reaction201034122950295510.1039/c0nj00283f2-s2.0-78649377426WangG.LiuaT.XieaX.RenbZ.BaicJ.WangH.Structure and electrochemical performance of Fe3O4/graphene nanocomposite as anode material for lithium-ion batteries2011128333634010.1016/j.matchemphys.2011.03.049FanW.GaoW.ZhangC.TjiuW. W.PanJ.LiuT.Hybridization of graphene sheets and carbon-coated Fe3O4 nanoparticles as a synergistic adsorbent of organic dyes20122248251082511510.1039/c2jm35609k2-s2.0-84870024671RobertsonA. W.LeeG.-D.HeK.YoonE.KirklandA. I.WarnerJ. H.The role of the bridging atom in stabilizing odd numbered graphene vacancies20141473972398010.1021/nl501320a2-s2.0-84903957104WangH.WangQ.ChengY.LiK.YaoY.ZhangQ.DongC.WangP.SchwingenschlöglU.YangW.ZhangX. X.Doping monolayer graphene with single atom substitutions201212114114410.1021/nl20316292-s2.0-84855793049BlöchlP. E.Projector augmented-wave method1994501795310.1103/physrevb.50.17953KresseG.FurthmullerJ.Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set199661155010.1016/0927-0256(96)00008-0PerdewJ. P.BurkeK.ErnzerhofM.Generalized gradient approximation made simple199677, article 386510.1103/physrevlett.77.38652-s2.0-4243943295