Li-Intercalated Graphene on Ir(111) - American Chemical Society

Feb 16, 2016 - Institut für Physik, Technische Universität Ilmenau, D-98693 Ilmenau, Germany. ABSTRACT: Graphene on Ir(111) was intercalated by Li. Th...
1 downloads 0 Views 9MB Size
Article pubs.acs.org/JPCC

Filling the Gap: Li-Intercalated Graphene on Ir(111) Johannes Halle,* Nicolas Néel, and Jörg Kröger Institut für Physik, Technische Universität Ilmenau, D-98693 Ilmenau, Germany ABSTRACT: Graphene on Ir(111) was intercalated by Li. The intercalant formed a variety of superstructures with increasing coverage. For all superstructures, the moiré pattern of pristine graphene on Ir(111) was discernible. At low coverage, the intercalant exhibited a disperse pattern of elementary assemblies that served as building blocks for Li-induced structures at higher coverage. In contrast to previous reports on other intercalants, Li intercalation gave rise to isotropic assemblies in the lowcoverage range. Pattern formation was qualitatively rationalized in terms of different energy costs and gains and a corresponding ranking of preferred intercalation regions. While for intermediate Li coverage an ordered intercalant superstructure was observed, higher coverage led to the decoration of step edges with adsorbed Li clusters.



fied, which served as building blocks for intercalation structures at higher coverage. An ordered phase of the intercalant was reflected by the atomically resolved p(√3 × √3)R 30° superstructure at monolayer coverage. The influence of the graphene moiré lattice was clearly discernible for all superstructures. Based on the observed intercalant structures, a ranking of preferred intercalation regions was proposed. At elevated exposures, Li formed clusters adsorbed atop graphene in the vicinity of step edges.

INTRODUCTION Pristine graphene exhibits a wealth of intriguing properties that render the sp2-bonded two-dimensional lattice of C atoms interesting for fundamental research and applications alike. The exceptionally high mechanical strength,1 ballistic electron transport at room temperature,2,3 and presence of massless Dirac Fermions4 represent only a few of the appealing characteristics. Excellent review articles thoroughly summarize the hallmarks of pristine and adsorbed graphene.5−8 Chemical modifications of graphene open the alley for tailoring additional properties and for novel applications. Epitaxially grown graphene on surfaces led to site-specific chemical reactivity.9−12 Doping graphene with donors or acceptors is one possible route to controlled modification of the electronic structure.13,14 Local pn junctions were realized in that way.15 Intercalation, or the transfer of material between epitaxial graphene and a substrate, is an emerging procedure to engineer the bonding to the substrate16−18 or to open a band gap.19 Intercalation is likewise relevant to the use of graphene in spin-electronic devices. Ferromagnetic materials,20−22 alkali metals, alkaline earth metals, and lanthanides23 were employed to tune the magnetic properties of graphene. In this context, graphite intercalation compounds are noteworthy since addition of alkali metals or rare earths was reported to induce superconductivity.24 Lithium adsorbed to pristine, or freestanding, graphene has been predicted to become superconducting at ∼8 K.25 Here a low-temperature scanning tunneling microscope (STM) was used to unravel structural aspects of Li-intercalated graphene on Ir(111). The main impetus to our investigations was the small size of Li compared to the size of other intercalants reported so far. Deviations from already observed intercalation structures were therefore expected and awaited their experimental verification. At the lowest Li coverage, small and isotropically shaped intercalation assemblies were identi© XXXX American Chemical Society



EXPERIMENTAL PROCEDURES Experiments were performed with an STM operated at 6 K and in ultrahigh vacuum (10−9 Pa). Clean Ir(111) was prepared by Ar+ bombardment and annealing. Each cycle was succeeded by O2 exposure at 10−5 Pa partial pressure and a sample temperature of 1270 K. A single layer of graphene was fabricated by exposing Ir(111) to C2H2 (purity 99.99%) with a partial pressure of 10−4 Pa and a sample temperature of 1470 K.26−29 While Li intercalation was observed already at room temperature30−32 the graphene-covered Ir(111) sample was exposed to Li at 630 K in the present experiments in order to ensure that the intercalant explores the entire surface for finding the energetically favored adsorption site. Higher sample temperatures were not used in order to avoid desorption of Li. For intercalated Li, desorption was reported for a temperature of 870 K.31 The coverage was increased by successively depositing Li for ∼4 min. Coverage was calibrated by assigning a coverage of 100% to the occupation of each graphene C6 ring with a single Li atom. Accordingly, the closed Li layer exhibiting the p(√3 × √3)R 30° superstructure reflects a coverage of 33%. All other coverage levels were then determined by Received: January 22, 2016 Revised: February 11, 2016

A

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

moiré lattice is 25.3 ± 0.2 Å, which indicates the alignment of ⟨112̅0⟩ graphene directions with crystallographic ⟨11̅0⟩ Ir(111) directions.35 On the basis of previous findings,34 we assign the highest protrusions [yellow dot in Figure 1a, yellow hexagon in Figure 1b] to C6 rings of the graphene lattice residing at facecentered cubic (fcc) Ir(111) sites. Slightly lower protrusions (red) correspond to C6 rings at hexagonal close-packed (hcp) Ir(111) sites. C6 rings at on-top Ir(111) sites give rise to depressions (blue). The sketch depicted in Figure 1b will hereafter be used to characterize the Li intercalation structures. Exposing graphene-covered Ir(111) to Li at 630 K led to STM images presented in Figure 2a,b. At low Li coverage of 3.7% (Figure 2a) and 6.1% (Figure 2b), a regular array of protrusions appeared in STM images. These protrusions exhibited an apparent height of 42 ± 3 pm and a diameter of 9.2 ± 1.0 Å, as extracted as the full width at half-maximum (fwhm) from cross-sectional profiles. The regular array of protrusions showed the same periodicity and orientation as the moiré lattice of pristine graphene. We assign these protrusions to intercalated Li at one specific moiré site and refer to them as dots in the following. It is reasonable to assume that this intercalation site is preferred over other regions of the moiré lattice. However, from STM images alone, the identification as fcc, hcp, or on-top region is difficult with present knowledge. Recent density functional calculations of Li-intercalated graphene on Ir(111) unraveled Ir(111) hcp hollow sites as preferred Li adsorption sites.30 In addition, these calculations showed that C6 rings of the graphene lattice are centered atop Li atoms,30 in agreement with predictions for freestanding graphene25 and graphite intercalation compounds.36,37 This arrangement corresponds to the hcp stacking region in the moiré unit cell. Therefore, we identify the moiré hcp regions as the preferred intercalation site. Figure 2c illustrates this assignment by enframing the occupied hcp sites (red) with thick lines. The number of atoms in an intercalant dot was estimated by multiplying the total number of carbon rings in the moiré unit cell by the determined coverage. As a result, dots at Li coverage of 3.7% contain assemblies with ∼4 Li atoms. The inset to Figure 2a shows the proposed arrangement of Li atoms (solid circles) beneath C6 rings (circles). Based on the atomic Li superstructure at monolayer coverage, 12 Li atoms represent the maximum occupation of a dot. Increasing the Li coverage led to an increase in lateral dot size (Figure 2b). The observed compact assemblies at this coverage exhibited a diameter of 11.5 ± 0.2 Å (fwhm of crosssectional profiles) and an apparent height of 44 ± 2 pm. The estimated average number of Li atoms occupying the dots at this coverage is ∼7. In addition, intercalant structures coalesced at step edges and gave rise to elongated Li accumulation regions at the lower sides of the steps. In contrast, pairs of dots and cloverlike assemblies formed at the upper side of step edges. The different Li intercalation behavior observed at lower and upper step edges may be due to the different sign of the graphene curvature in both regions.38 In particular, upward bending of graphene at lower edges facilitates accumulation of Li due to the interplay of different energy contributions as explained below. An additional asymmetry between these two regions is caused by the Smoluchowski effect,39 which leads to a charge accumulation (depletion) at lower (upper) step edges. Alkali metals are known to exhibit ionic bonding to the substrate surface; that is, the weakly bonded valence electron is readily transferred to the substrate.40 Indeed, recent calculations showed that Li intercalation is accompanied by

calculating the fraction of surface areas occupied by Li over the total area multiplied by 33%. For Li intercalation at room temperature, a (1 × 1) superstructure with respect to Ir was reported.31 This superstructure reflects Li coverage up to 100%. All STM images were recorded at constant current and with the bias voltage applied to the sample. STM data were processed by use of WSxM.33



RESULTS AND DISCUSSION Pristine Graphene and Low Li Coverage. Figure 1a shows an STM image of clean graphene on Ir(111). Due to the mismatch of lattice constants [graphene 2.452 Å, Ir(111) 2.715 Å],34 a moiré superstructure is formed whose lattice vectors are indicated as m1 and m2 in Figure 1a. The periodicity of the

Figure 1. (a) STM image of graphene on Ir(111), showing the moiré lattice (−0.2 V, 0.1 nA, 20 × 20 nm2) with lattice vectors m1 and m2. Yellow, red, and blue dots indicate fcc, hcp, and on-top regions of the moiré pattern, respectively. (b) Sketch of graphene on Ir(111). C6 rings are indicated as orange circles. Ir atoms appear as light gray (first layer) and dark gray (second layer) dots. Colored hexagons indicate high-symmetry regions of the moiré lattice. Crystallographic directions of Ir and graphene are shown in the bottom left corner. B

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 2. (a, b) STM images of Li-intercalated graphene at Li coverage of (a) 3.7% and (b) 6.1% (0.1 V, 0.1 nA, 40 × 40 nm2). Regular arrays of dots have the same periodicity and orientation as the moiré pattern of pristine graphene on Ir(111) and are assigned to Li assemblies residing at hcp sites of the moiré lattice. White arrows indicate the moiré ⟨112̅0⟩ direction. (a, inset) Suggested arrangement of Li atoms (solid circles) beneath C6 rings (open circles) in dots. The suggestion is based on the Li superstructure determined at higher coverage. (c) Sketch of moiré lattice (colored hexagons); hcp regions occupied by Li assemblies are indicated by thick lines.

substantial charge transfer from the intercalant to the Ir substrate as well as to graphene.30 The resulting partial positive charge of the intercalated Li atoms may therefore explain their preferred diffusion to lower step edges. Similar findings were reported for molecule adsorption on noble metal surfaces.41,42 Intermediate Li Coverage. In the coverage range presented next, the formation of cloverlike intercalation assemblies and the pairing of dots to elongated stripes was observed on terraces (Figure 3a). At lower coverage, they occurred solely at the upper side of step edges (Figure 2b). The occupation of the moiré lattice at intermediate Li coverage is illustrated in Figure 3b. The cloverlike assemblies exhibited an apparent height of their central part of 77 ± 2 pm and a distance between dots at their periphery of 24.8 ± 0.7 Å. Figure 4a demonstrates the preference for a uniform orientation. Only in rare cases was the opposite orientation observed, as shown by the close-up view in Figure 4b. Here, two cloverlike assemblies that are rotated by 180° with respect to the prevailing orientation are presented. Li atoms at the center of these assemblies are proposed to occupy a second moiré lattice region. As illustrated in the inset to Figure 4b, the orientation of such an assembly may be rotated by 180° when the central moiré site changes from fcc to on-top. According to the previously determined orientation of the moiré lattice (Figure 1) and the prevailing orientation of the cloverlike structures (Figure 4), the second occupied moiré region corresponds to the moiré fcc site. The observed stripes displayed an apparent height of 54 ± 4 pm, and the fwhm of cross-sectional profiles was 9.2 ± 0.2 Å. Stripe orientations correspond to all three high-symmetry directions of the moiré lattice. Occasionally coalesced stripes enclosed an angle of 120°, while angles of 60° were not observed. In the intermediate coverage range, stripes preponderantly formed in the vicinity of defects, which likely induced an anisotropic strain distribution in graphene. In Figure 4a, the top left corner of the STM image shows a graphene blister, which led to intercalation stripes oriented along the moiré lattice directions and pointing toward the blister. This indicates a directional influence of anisotropic strain on stripe orientation. At distances from the defect exceeding ∼30 nm, the cloverlike assemblies prevailed again, indicating isotropically strained areas. Similar conclusions were drawn for intercalation of Eu under graphene on Ir(111), where the density of intercalation stripes proved to be an indicator of preexisting strain in the

Figure 3. (a) STM image of Li-intercalated graphene at 14.8% Li coverage (0.1 V, 0.1 nA, 80 × 80 nm2). The central region of cloverlike assemblies is the moiré fcc site, which represents the second favorable Li intercalation site.

graphene sheet.43 However, in that case isotropically shaped intercalant assemblies were not observed, even in isotropically C

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 4. (a) STM image of Li-intercalated graphene at Li coverage of 16.7% (0.1 V, 0.1 nA, 80 × 80 nm2). The top left corner shows a graphene blister, which represents a defect of the regular moiré lattice. In its vicinity, that is, within distances of ∼30 nm, stripe patterns along moiré lattice directions formed, which indicate the locally anisotropic strain in graphene. The strain becomes isotropic again in regions where cloverlike intercalation assemblies dominate. (b) STM image of Li-intercalated graphene at Li coverage of 14.8% (0.1 V, 0.1 nA, 25 × 25 nm2). Two cloverlike intercalation assemblies that are rotated by 180° with respect to the prevailing orientation are marked by dashed white circles. (Inset) Sketch of two cloverlike Li intercalation assemblies with opposite orientation. By exchanging the central region of the assembly from fcc (yellow, right) to on-top (blue, left) the orientation of the assembly rotates by 180°. White arrows indicate the moiré ⟨112̅0⟩ direction.

Figure 5. (a) STM image of Li-intercalated graphene at Li coverage of 25.1% (0.2 V, 10 pA, 100 × 100 nm2). An Ir(111) step edge separates the upper (left) from the lower (right) terrace. Dark lines are due to regions beneath graphene that are not occupied by Li intercalants at this coverage. The regular array of circular protrusions is due to the graphene moiré lattice. (b) Close-up view of an intercalated region (−0.1 V, 1 μA, 5.2 × 5.2 nm2). The ⟨112̅0⟩ direction is indicated by the white arrow. Ordered patterns with different periodicities are visible. The smallest periodicity is due to the graphene atomic lattice. The hexagonal arrangement of depressions indicated by white circles is a Li-induced superstructure explained in the text. The large circular protrusions are due to the graphene moiré pattern. (c) Two-dimensional Fourier transform of the STM image in panel b. Fourier maxima marked with orange circles are due to the graphene lattice. Additional spots (dashed red circles) are due to the Li superstructure. Both lattices enclose an angle of 30°.

local strain relief, Er, and an energy cost due to graphene bending, Eb:

strained areas. Here, the small size of Li leads to a decrease in the energy gained due to strain relief at edges of intercalated assemblies. Therefore, formation of isotropic assemblies with cloverlike shape is enabled. The observed formation of dots, pairs of dots, stripes, and cloverlike intercalation structures may be explained in terms of different energy contributions.43 To this end, we define a site energy per intercalated Li atom:

Ee = Er − E b

With the help of the orientation preference of cloverlike Li assemblies, a ranking of energetically favored intercalation regions is suggested. The peripheral dots correspond to the initially occupied moiré sites at low Li coverage and thus represent the most favorable intercalation site. On the basis of recent calculations,30 this intercalation region is identified as the moiré hcp site. The center of the cloverlike arrangement is then the second favorable intercalation region and corresponds to the moiré fcc site, leaving the remaining on-top site as the least favorable intercalation region. In terms of Es, this ranking on‑top > Efcc . corresponds to Ehcp s s > Es

Es = Ea − Ed which includes the adsorption energy, Ea, as an energy gain and the delamination energy, Ed, as an energy cost. The delamination energy takes into account the lifting of the graphene sheet due to intercalation. In addition, intercalation creates graphene edges surrounding the Li assemblies. The edge energy per edge length, Ee, includes an energy gain due to D

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 6. STM image of adsorbed Li clusters at the lower edge of a step at (a) 15% coverage (0.3 V, 0.4 nA, 26 × 26 nm2) and (b) 32% coverage (0.1 V, 0.1 nA, 26 × 26 nm2). (b) The clusters grow with the same periodicty and orientation as the moiré lattice. Compared to the moiré lattice on the flat terrace, whose unit cells are indicated by the white lines, the cluster lattice is shifted by two moiré sites (≈29 Å) in ⟨011̅0⟩ direction, indicated by the dashed arrow. (Inset) Illustration of the shift from preferred intercalation (hcp) to preferred adsorption (fcc) site. Solid arrows indicate the moiré ⟨112̅0⟩ direction.

with Li atoms. An almost closed intercalation film was formed at a coverage of ∼33%. The intercalation film was imaged with atomic resolution. Figure 5b shows a close-up view of a representative region. The pattern with the smallest spatial periodicity (circular protrusions) is due to the atomic graphene lattice, while the structure with the largest periodicity (bright circular regions) is due to the moiré lattice. An additional superstructure is visible as depressions, which are indicated by the white circles in Figure 5b. With respect to the graphene lattice, this superstructure is described by p(√3 × √3)R 30°. The lattice vectors spanning the unit cell exhibit a length of 4.3 ± 0.1 Å, which deviates by 1% from the expected p(√3 × √3)R 30° lattice constant of 4.24 Å. The rotation angle was best determined from the Fourier transform of the STM image (Figure 5c). Graphene and superstructure Fourier spots are indicated by solid and dashed circles, respectively. The p(√3 × √3)R 30° superstructure was retained for the intercalated Li monolayer at ∼33% coverage. The same superstructure was reported for Li intercalation of graphene on SiC32 and graphite intercalation compounds with alkali metal interlayers.37 In addition, the corrugation of graphene in these intercalated regions is 9 ± 1 pm only, which is about half the value obtained for pristine graphene on Ir(111), 19 ± 1 pm.47 Previously, the moiré corrugation was identified as an appropriate measure of graphene−substrate hybridization.17,48 Therefore, the observed low corrugation hints at a low hybridization of graphene with the Li intercalation layer. This is in contrast to many other metal intercalants and substrates,17,20−22 which showed large moiré corrugations of 60−180 pm. In these reports the intercalants adjusted their lattice constants to the substrate lattice constant, giving rise to pseudomorphic growth accompanied by a (1 × 1) superstructure with respect to the substrate lattice, as opposed to the nonpseudomorphic p(√3 × √3)R 30° superstructure with respect to graphene reported here. Starting from Li coverage of ∼15%, regularly spaced protrusions of 100−150 pm apparent height and average

Our observations suggest that intercalation sites and orientation of Li assemblies in the absence of anisotropic strain are determined by Es, while Ee defines size and shape of the intercalated assemblies. In areas with large anisotropic strain, Er forces the intercalated assemblies into stripe shape,43 with an orientation given by the strain distribution. The findings presented here deviate from observations for Eu43 and CO44 intercalation under graphene on Ir(111), where moiré on-top sites were identified as preferred intercalation regions. This conclusion was based on the derived ranking of hcp 43 < Efcc Consistently, delamination energies Eon‑top d d < Ed . density functional calculations of graphene on Ir(111) showed that fcc and hcp regions are characterized by weak covalent bonds, while moiré on-top regions exhibit repulsive interaction.45 In the case of Li intercalation, however, it was suggested that the energy costs for delamination, Ed, are significantly reduced,31 owing to the small size of Li compared to most intercalants examined so far.17,21,22,31,43,44,46 Thus, a larger influence of the adsorption energy Ea on Es may be expected and, as shown by our intercalation experiments, leads to a different ranking of preferred moiré lattice sites for intercalation. For the same reason, Li requires less energy Eb for bending graphene, which leads to the appearance of disperse intercalation patterns consisting of many elementary Li assemblies at low to intermediate coverage, rather than to compact islands.43,44 High Li Coverage. Increasing Li coverage led to the coalescence of dots, stripes, and cloverlike assemblies. The resulting intercalation structure at a coverage of 25.1% is shown in the STM image of Figure 5a. It appears rather flat with embedded elongated depressions, which exhibit a width of 14.7 ± 0.1 Å, corresponding approximately to the average width of a moiré site (14.6 Å). We assign the elongated depressions to nonintercalated regions. Some direction preference was observed, which is probably due to the influence of anisotropic strain again. Upon further increase of Li coverage, the number of elongated depressions was reduced due to their occupation E

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C diameter of 15.6 ± 0.4 Å began to decorate the lower side of step edges (Figure 6a). Their mutual distance was virtually identical with the graphene moiré periodicity. With increasing coverage, the protrusions formed a hexagonal lattice with the moiré periodicity and orientation (Figure 6b). Compared to the moiré superstructure on the flat (lower) terrace, the superimposed new lattice was laterally shifted by two moiré sites (≈29 Å) in ⟨011̅0⟩ direction (dashed arrow in Figure 6b). Due to their variation in size and shape, we assign the protrusions to nonintercalated, i.e., adsorbed Li clusters. An adsorption phase was likewise reported for Eu23,43 and Cs23,31 on graphene-covered Ir(111). On the basis of calculations and photoelectron spectroscopy data for Li-intercalated graphene on Ir(111), however, adsorbed Li was not indicated.31 Since the calculations did not consider step edges, and due to the laterally averaged photoemission data, we think that this previous report31 is not necessarily at odds with our findings, which concern adsorbed Li clusters in the vicinity of step edges. The observed lateral shift of the array of clusters with respect to the moiré lattice (Figure 6) implies that adsorbed Li clusters preferably occupy fcc regions of the moiré superstructure. Ordered adsorption phases were likewise reported for Ir49 and H50 on nonintercalated graphene on Ir(111). For Ir clusters, hcp moiré regions were identified as preferred adsorption sites. In summary, the main findings of our paper reveal that, at low coverage, the small size of the Li intercalant leads to a disperse distribution of intercalation assemblies. In addition, the energy gain due to strain relief upon edge formation is low in isotropically strained areas, which enables isotropic cloverlike shapes. This observation is in contrast to previously studied larger intercalants, such as Eu23,43 and Cs.23,31 In accordance with recent calculations30 our results are compatible with Li preferably occupying hcp stacking regions of the moiré lattice. Based on the prevailing orientation of intercalated assemblies, fcc and on-top regions represent the second and third favorite intercalation sites, respectively. Larger intercalants were reported to adsorb at on-top regions first and then at fcc and hcp regions.43,44 For Li the graphene delamination energy is reduced,31 owing to the small size of the intercalant, which leads to adsorption energy rather than delamination energy being decisive for the different ranking of preferred adsorption regions. Similar to results reported for Li intercalation of graphene-covered SiC32 and for graphite intercalation compounds with alkali metal interlayers,37 a p(√3 × √3)R 30° superstructure was found with respect to graphene. However, the graphene corrugation in this elevated-coverage range is nearly an order of magnitude lower than for intercalants displaying pseudomorphic growth.17,21,22,46 Low hybridization of graphene with the Li layer and/or the absence of pseudomorphic Li growth may be at the origin of this observation.17 At high coverage, Li additionally exhibits an adsorbed phase, which is reflected by the decoration of step edges with regularly spaced Li clusters residing at fcc stacking regions of the moiré lattice.

leads to a low edge energy of intercalated assemblies reflected by their isotropic shape at low coverage. Moreover, our observations evidence that an adsorbed phase may form preferably at step edges. The presented findings may serve as a structural basis for forthcoming experimental and theoretical investigations of electronic and dynamic properties of Li-intercalated graphene as well as intercalation of other two-dimensional materials. Disperse distributions of intercalated assemblies could be of interest for nanopatterning of an adsorbed phase. The emerging low hybridization of graphene with the Li intercalant layer may render Li-intercalated graphene on Ir(111) particularly appealing for the efficient decoupling of adsorbed molecules.



AUTHOR INFORMATION

Corresponding Author

*E-mail [email protected]. Notes

The authors declare no competing financial interest.

■ ■

ACKNOWLEDGMENTS Financial support by the Deutsche Forschungsgemeinschaft through Grant KR 2912/10-1 is acknowledged. REFERENCES

(1) Lee, C.; Wei, X.; Kysar, J. W.; Hone, J. Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene. Science 2008, 321, 385−388. (2) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666−669. (3) Bolotin, K. I.; Sikes, K. J.; Hone, J.; Stormer, H. L.; Kim, P. Temperature-Dependent Transport in Suspended Graphene. Phys. Rev. Lett. 2008, 101, No. 096802. (4) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Katsnelson, M. I.; Grigorieva, I. V.; Dubonos, S. V.; Firsov, A. A. TwoDimensional Gas of Massless Dirac Fermions in Graphene. Nature 2005, 438, 197−200. (5) Rao, C. N. R.; Sood, A. K.; Subrahmanyam, K. S.; Govindaraj, A. Graphene: The New Two-Dimensional Nanomaterial. Angew. Chem., Int. Ed. 2009, 48, 7752−7777. (6) Neto, A. H. C.; Guinea, F.; Peres, N. M. R.; Novoselov, K. S.; Geim, A. K. The Electronic Properties of Graphene. Rev. Mod. Phys. 2009, 81, 109−162. (7) Wintterlin, J.; Bocquet, M.-L. Graphene on Metal Surfaces. Surf. Sci. 2009, 603, 1841−1852. (8) Kotov, V. N.; Uchoa, B.; Pereira, V. M.; Guinea, F.; Castro Neto, A. H. Electron-Electron Interactions in Graphene: Current Status and Perspectives. Rev. Mod. Phys. 2012, 84, 1067−1125. (9) Donner, K.; Jakob, P. Structural Properties and Site Specific Interactions of Pt with the Graphene/Ru(0001) Moiré Overlayer. J. Chem. Phys. 2009, 131, 164701. (10) Pan, Y.; Gao, M.; Huang, L.; Liu, F.; Gao, H.-J. Directed SelfAssembly of Monodispersed Platinum Nanoclusters on Graphene Moiré Template. Appl. Phys. Lett. 2009, 95, 093106. (11) Zhou, Z.; Gao, F.; Goodman, D. W. Deposition of Metal Clusters on Single-Layer Graphene/Ru(0001): Factors that Govern Cluster Growth. Surf. Sci. 2010, 604, L31−L38. (12) Altenburg, S. J.; Kröger, J.; Wang, B.; Bocquet, M.-L.; Lorente, N.; Berndt, R. Graphene on Ru(0001): Contact Formation and Chemical Reactivity on the Atomic Scale. Phys. Rev. Lett. 2010, 105, No. 236101. (13) McChesney, J. L.; Bostwick, A.; Ohta, T.; Seyller, T.; Horn, K.; González, J.; Rotenberg, E. Extended van Hove Singularity and Superconducting Instability in Doped Graphene. Phys. Rev. Lett. 2010, 104, No. 136803.



CONCLUSIONS Our results show that the size of an intercalant matters for pattern formation and selection of intercalation sites. A reduction of graphene delamination energy costs for small intercalants like Li enhances the influence of the adsorption energy gain. Thus, the adsorption geometry with respect to graphene and substrate becomes the main factor in determining preferred intercalation regions. In addition, the small size of Li F

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C (14) Liu, H.; Liu, D.; Zhu, Y. Chemical Doping of Graphene. J. Mater. Chem. 2011, 21, 3335−3345. (15) Emtsev, K. V.; Zakharov, A. A.; Coletti, C.; Forti, S.; Starke, U. Ambipolar Doping in Quasifree Epitaxial Graphene on SiC(0001) Controlled by Ge Intercalation. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 84, No. 125423. (16) Riedl, C.; Coletti, C.; Iwasaki, T.; Zakharov, A. A.; Starke, U. Quasi-Free-Standing Epitaxial Graphene on SiC Obtained by Hydrogen Intercalation. Phys. Rev. Lett. 2009, 103, No. 246804. (17) Pacilé, D.; Leicht, P.; Papagno, M.; Sheverdyaeva, P. M.; Moras, P.; Carbone, C.; Krausert, K.; Zielke, L.; Fonin, M.; Dedkov, Y. S.; et al. Artificially Lattice-Mismatched Graphene/Metal Interface: Graphene/Ni/Ir(111). Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, No. 035420. (18) Papagno, M.; Moras, P.; Sheverdyaeva, P. M.; Doppler, J.; Garhofer, A.; Mittendorfer, F.; Redinger, J.; Carbone, C. Hybridization of Graphene and a Ag Monolayer Supported on Re(0001). Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, No. 235430. (19) Enderlein, C.; Kim, Y. S.; Bostwick, A.; Rotenberg, E.; Horn, K. The Formation of an Energy Gap in Graphene on Ruthenium by Controlling the Interface. New J. Phys. 2010, 12, 033014. (20) Gyamfi, M.; Eelbo, T.; Waśniowska, M.; Wiesendanger, R. Impact of Intercalated Cobalt on the Electronic Properties of Graphene on Pt(111). Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 85, No. 205434. (21) Sicot, M.; Leicht, P.; Zusan, A.; Bouvron, S.; Zander, O.; Weser, M.; Dedkov, Y. S.; Horn, K.; Fonin, M. Size-Selected Epitaxial Nanoislands Underneath Graphene Moiré on Rh(111). ACS Nano 2012, 6, 151−158. (22) Decker, R.; Brede, J.; Atodiresei, N.; Caciuc, V.; Blügel, S.; Wiesendanger, R. Atomic-Scale Magnetism of Cobalt-Intercalated Graphene. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 87, No. 041403. (23) Schumacher, S.; Wehling, T. O.; Lazić, P.; Runte, S.; Förster, D. F.; Busse, C.; Petrović, M.; Kralj, M.; Blügel, S.; Atodiresei, N.; et al. The Backside of Graphene: Manipulating Adsorption by Intercalation. Nano Lett. 2013, 13, 5013−5019. (24) Csányi, G.; Littlewood, P. B.; Nevidomskyy, A. H.; Pickard, C. J.; Simons, B. D. The Role of the Interlayer State in the Electronic Structure of Superconducting Graphite Intercalated Compounds. Nat. Phys. 2005, 1, 42−45. (25) Profeta, G.; Calandra, M.; Mauri, F. Phonon-Mediated Superconductivity in Graphene by Lithium Deposition. Nat. Phys. 2012, 8, 131−134. (26) Hattab, H.; N'Diaye, A. T.; Wall, D.; Jnawali, G.; Coraux, J.; Busse, C.; van Gastel, R.; Poelsema, B.; Michely, T.; zu Heringdorf, F.J. M.; et al. Growth Temperature Dependent Graphene Alignment on Ir(111). Appl. Phys. Lett. 2011, 98, No. 141903. (27) Endlich, M.; Molina-Sánchez, A.; Wirtz, L.; Kröger, J. Screening of Electron-Phonon Coupling in Graphene on Ir(111). Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, No. 205403. (28) Endlich, M.; Miranda, H. P. C.; Molina-Sánchez, A.; Wirtz, L.; Kröger, J. Moiré-Induced Replica of Graphene Phonons on Ir(111). Ann. Phys. 2014, 526, 372−380. (29) Endlich, M.; Gozdzik, S.; Néel, N.; da Rosa, A. L.; Frauenheim, T.; Wehling, T. O.; Kröger, J. Phthalocyanine Adsorption to Graphene on Ir(111): Evidence for Decoupling from Vibrational Spectroscopy. J. Chem. Phys. 2014, 141, 184308. (30) Pervan, P.; Lazić, P.; Petrović, M.; Šrut Rakić, I.; Pletikosić, I.; Kralj, M.; Milun, M.; Valla, T. Li Adsorption versus Graphene Intercalation on Ir(111): From Quenching to Restoration of the Ir Surface State. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, No. 245415. (31) Petrović, M.; Rakić, I. Š; Runte, S.; Busse, C.; Sadowski, J. T.; Lazić, P.; Pletikosić, I.; Pan, Z.-H.; Milun, M.; Pervan, P.; et al. The Mechanism of Caesium Intercalation of Graphene. Nat. Commun. 2013, 4, No. 2772.

(32) Virojanadara, C.; Watcharinyanon, S.; Zakharov, A. A.; Johansson, L. I. Epitaxial Graphene on 6H-SiC and Li Intercalation. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 82, No. 205402. (33) Horcas, I.; Fernández, R.; Gómez-Rodríguez, J. M.; Colchero, J.; Gómez-Herrero, J.; Baro, A. M. WSXM: A Software for Scanning Probe Microscopy and a Tool for Nanotechnology. Rev. Sci. Instrum. 2007, 78, 013705. (34) N’Diaye, A. T.; Coraux, J.; Plasa, T. N.; Busse, C.; Michely, T. Structure of Epitaxial Graphene on Ir(111). New J. Phys. 2008, 10, 043033. (35) Slight misalignments of the moiré lattice orientation were observed between different regions on the sample with a standard deviation of 5°. This is due to small angles between graphene and Ir(111) unit cell vectors, which are magnified by a factor of ≈10.6 for the moiré lattice orientation. (36) DiVincenzo, D. P.; Mele, E. J. Cohesion and Structure in Stage1 Graphite Intercalation Compounds. Phys. Rev. B: Condens. Matter Mater. Phys. 1985, 32, 2538−2553. (37) Dresselhaus, M. S.; Dresselhaus, G. Intercalation Compounds of Graphite. Adv. Phys. 1981, 30, 139−326. (38) Coraux, J.; N’Diaye, A. T.; Busse, C.; Michely, T. Structural Coherency of Graphene on Ir (111). Nano Lett. 2008, 8, 565−570. (39) Smoluchowski, R. Anisotropy of the Electronic Work Function of Metals. Phys. Rev. 1941, 60, 661−674. (40) Ertl, G. In Physics and Chemistry of Alkali Metal Adsorption; Materials Science Monographs, Vol. 57; Bonzel, H. P., Bradshaw, A. M., Ertl, G., Eds.; Elsevier Science Publishers: Amsterdam, 1989; pp 1−10. (41) Stranick, S.; Kamna, M.; Weiss, P. Interactions and Dynamics of Benzene on Cu111 at Low Temperature. Surf. Sci. 1995, 338, 41−59. (42) Rogero, C.; Pascual, J. I.; Gómez-Herrero, J.; Baró, A. M. Resolution of Site-Specific Bonding Properties of C60 Adsorbed on Au(111). J. Chem. Phys. 2002, 116, 832−836. (43) Schumacher, S.; Förster, D. F.; Rösner, M.; Wehling, T. O.; Michely, T. Strain in Epitaxial Graphene Visualized by Intercalation. Phys. Rev. Lett. 2013, 110, No. 086111. (44) Grånäs, E.; Andersen, M.; Arman, M. A.; Gerber, T.; Hammer, B.; Schnadt, J.; Andersen, J. N.; Michely, T.; Knudsen, J. CO Intercalation of Graphene on Ir(111) in the Millibar Regime. J. Phys. Chem. C 2013, 117, 16438−16447. (45) Busse, C.; Lazić, P.; Djemour, R.; Coraux, J.; Gerber, T.; Atodiresei, N.; Caciuc, V.; Brako, R.; N’Diaye, A. T.; Blügel, S.; et al. Graphene on Ir(111): Physisorption with Chemical Modulation. Phys. Rev. Lett. 2011, 107, No. 036101. (46) Sicot, M.; Fagot-Revurat, Y.; Kierren, B.; Vasseur, G.; Malterre, D. Copper Intercalation at the Interface of Graphene and Ir(111) Studied by Scanning Tunneling Microscopy. Appl. Phys. Lett. 2014, 105, 191603. (47) This average value was obtained from several images recorded at currents ranging from 50 nA to 5 μA. Measured corrugations varied by 20% at most. (48) Preobrajenski, A. B.; Ng, M. L.; Vinogradov, A. S.; Mårtensson, N. Controlling Graphene Corrugation on Lattice-Mismatched Substrates. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78, No. 073401. (49) N’Diaye, A. T.; Bleikamp, S.; Feibelman, P. J.; Michely, T. TwoDimensional Ir Cluster Lattice on a Graphene Moiré on Ir(111). Phys. Rev. Lett. 2006, 97, No. 215501. (50) Balog, R.; Jørgensen, B.; Nilsson, L.; Andersen, M.; Rienks, E.; Bianchi, M.; Fanetti, M.; Lægsgaard, E.; Baraldi, A.; Lizzit, S.; et al. Bandgap Opening in Graphene Induced by Patterned Hydrogen Adsorption. Nat. Mater. 2010, 9, 315−319.

G

DOI: 10.1021/acs.jpcc.6b00729 J. Phys. Chem. C XXXX, XXX, XXX−XXX