Modulating the Electronic Properties of Graphene by Self-Organized

Sep 25, 2018 - State Key Laboratory of Low-Dimensional Quantum Physics and ... self-organization of S superlattices at the interface of graphene and S...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIVERSITY OF TOLEDO LIBRARIES

Article

Modulating the Electronic Properties of Graphene by Self-Organized S Identical Nanoclusters and Atomic Superlattices Confined at Interface Donglin Ma, Zhongqiu Fu, Xuelei Sui, Keke Bai, Jiabin Qiao, Chao Yan, Yu Zhang, Jingyi Hu, Qian Xiao, Xinrui Mao, Wenhui Duan, and Lin He ACS Nano, Just Accepted Manuscript • DOI: 10.1021/acsnano.8b04874 • Publication Date (Web): 25 Sep 2018 Downloaded from http://pubs.acs.org on September 26, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Modulating the Electronic Properties of Graphene by Self-Organized S Identical Nanoclusters and Atomic Superlattices Confined at Interface ⊥ Donglin Ma†,‡,§, Zhongqiu Fu†,§, Xuelei Sui , Keke Bai†, Jiabin Qiao†, Chao Yan†, Yu ⊥ † † † Zhang , Jingyi Hu , Qian Xiao , Xinrui Mao†, Wenhui Duan , and Lin He†,*



Center for Advanced Quantum Studies, Department of Physics, Beijing Normal

University, Beijing, 100875, People’s Republic of China ‡

Department of Physics, Capital Normal University, Beijing, 100048, People’s

Republic of China ⊥

State Key Laboratory of Low-Dimensional Quantum Physics and Collaborative

Innovation Center of Quantum Matter, Department of Physics, Tsinghua University, Beijing 100084, People’s Republic of China

ABSTRACT Ordered atomic-scale superlattices on surface hold great interest both for basic science and for potential applications in advanced technology. However, controlled fabrication of superlattices down to atomic scale has proven exceptionally challenging. Here we develop a segregation method to realize self-organization of S superlattices at the interface of graphene and S-rich Cu substrates. Via scanning tunneling microscope measurements, we directly image well-ordered identical nanoclusters superlattices and atomic superlattices under the cover of graphene. Scanning tunneling spectra show that the superlattices in turn could modulate the electronic structure of top-layer graphene. Importantly, a special-ordered S monoatomic superlattice commensurate with graphene lattice is found to drive semi-metal graphene into a symmetry-broken phase— —the electronic Kekulé distortion (KD) phase, which opens a bandgap of ~245 meV. KEYWORDS: graphene, interface, nanocluster superlattices, atomic superlattices, Kekulé distortion, scanning tunneling microscope

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Creation and control of periodic atomic-scale structures to build nanodevices and nanosystems meet the emerging need of nanoelectronics and data storage.1–3 Thanks to the invention of scanning tunneling microscope (STM),4–6 it becomes possible to fabricate atomic superlattices on metal surfaces by manipulation of STM tip.7–12 However, the drawback of tip manipulation is obvious: the atomic-scale precision is almost impossible to realize on a large scale. Another promising route is self-organized growth of periodic adatoms on designer surfaces,13–18 which has advantages in fabrication of large-area atomic superlattices with tunable size and periodicity.19 However, there are still two severe difficulties that need to overcome for the application of the atomic superlattices. One is that the atomic superlattices are usually only stable at low temperature and become unstable at high temperature due to thermal disturbance. The other is that the atomic superlattices on surface are easy to be destroyed because of absorption of other different atoms/molecules. Therefore, most of these delicately-built periodic atomic structures were realized in cryogenic environment and ultrahigh vacuum (UHV), which in principle hinder their potential application. The idea of self-organization of atomic superlattice at the interface between graphene and the supporting substrate naturally overcome the above-said two shortcomings of the atomic superlattices at the surface. Interface could be seen as a Z-confined two-dimensional space compared to the open space of surface. The confined space provides extra stability against thermal disturbance and protects the "trapped" atoms or molecules at the interface from contaminations.20 Intercalation of different atoms into the graphene-substrate interface have been reported in several graphene supported systems.21–41 However, the intercalation process has to be accomplished below the etching temperature of graphene (600~1000 K in UHV), which hinder the further evolution of intercalated atoms into large-area atomic superlattices.

ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

In this letter, we adopt an alternative strategy to introduce adatoms into graphene-substrate interface. Adatoms are designed to segregate from the bulk of substrate under the cover of graphene at graphene growth temperature (~1300 K), which process occurs above the temperature of the intercalation process reported previously. With the assist of high temperature, we realize self-organized S nanoclusters superlattices and atomic superlattices at interface between graphene and S-rich Cu substrate. Scanning tunneling spectra show that the superlattices in turn could modulate the electronic structure of top-layer graphene. Most intriguingly, our experiment demonstrates that the S monoatomic superlattice in special order can tune the band structure of graphene, driving graphene into the electronic Kekulé distortion (KD) phase with sizeable bandgap.

RESULTS AND DISCUSSION Figure 1 schematically shows chemical vapor deposition (CVD)-based growth process of graphene monolayer on Cu substrates. To introduce S adatoms confined at the interface between graphene and the substrates, we use S-rich Cu foils as the supporting substrates.16,42–44 In our experiment, the S-rich Cu substrates were annealed at high temperature as the first step to activate the S atoms both on the metal surface and in the bulk. In the second step, carbon sources were introduced into the system for graphene growth, then the sample was slowly cooled to room temperature (~20 ℃/min). The S atoms segregate from the Cu foils during the growth process, which is confirmed by our X-ray photoelectron spectroscopy (XPS) measurement (Fig. S1), and form two types of ordered S superlattices at the interface, as demonstrated subsequently in Fig. 2 and Fig. 4. The as-grown sample was delivered into UHV chamber from atmosphere and then annealed at ~650 K for degassing. All STM measurements were performed at 4.2 K unless otherwise noted. Figure 2 presents STM characterizations of the first type of S superlattice at the interface. As shown in Fig. 2a, large-area square-ordered superlattice of the intercalated S atoms over several hundreds of nanometers can be observed in our STM

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 22

measurement. The unit cell of the superlattice is rectangle with anisotropic periodicities of ~1.2 nm and ~1.5 nm respectively. In the zoom-in image (white dashed box in Fig. 2a), typical honeycomb lattices of the topmost graphene sheet were observed with blurred undercover features (Fig. 2b). The Fourier transform (FT) of the image (inset in Fig. 2b) clearly reflects the detailed periodicity of graphene lattice (the outer black-circled points) and the square-like S nanocluster superlattice (the inner red-circled points), as well as their spatial relationship that the superlattice is in line (or 90°) with graphene lattice in one crystalline direction. At low tunneling bias, atomic-resolved STM image of the intercalated S nanometer was obtained as shown in Fig. 2c. The nanometer cluster looks like the ordered assembly of S adatoms in a special “dice-five” shape (red circles in white box of Fig. 2c) with a size of ~0.5 nm, which is similar as the structure of S clusters on Cu(100).45 Such a result, along with the result obtained by our XPS measurement, demonstrates explicitly that the studied sample is the graphene monolayer on Cu substrate and the interface is intercalated with S atoms. It is also noted that the apparent height of S clusters is only ~30 pm, consistent with the apparent height of S atoms on Cu(100) (Height profile seen in Fig. S2).46 Based on the above analysis, detailed DFT calculations including the Cu(100) substrate were carried out to unravel the atomic-scale structures of nanoclusters. After an extensive structural search and comparison with high-resolution STM images, an energetically favorable model is found as shown in Fig. 2d. From this DFT relaxed model, we distinguish that the nanocluster could be a c(2×2)-like S5 assembly. The simulated STM image of this proposed structure (Fig. 2e) show identical “dice-five” topography, which fits well with the experiment. It is worth to mention that the observed S nanocluster superlattice on Cu(100) are topologically alike to p(2×2)-S reconstruction or clover-like √17 × √1714° -S reconstruction in previous literatures,45–49 but quite different in symmetry, periodicity and fine structures. It indicates that graphene templating effect take place in the formation of this new S nanocluster superlattice. All above evidences show that segregated S adatoms spontaneously assemble themselves into identical nanometer

ACS Paragon Plus Environment

clusters, and

Page 5 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

self-organize into large-area hierarchical superlattice at the interface between graphene and Cu substrate. Importantly, the obtained nanoclusters superlattice is quite robust and still stable after several cycles of air exposure and degas annealing (~650 K), as shown in Fig. S3. To further explore mechanism of the superlattice formation at interfaces, we carefully studied structures around boundaries of the superlattice. Figure 2f exhibits typical STM morphology around the edges of the nanocluster superlattices. It is surprisingly found that the isolated S nanocluster “vacancies” (red squares in Fig. 2f) or even individual S nanoclusters (Fig. S4) at the edges still strongly follow the rectangle symmetry of the superlattice. This observation indicates that the Cu(100) surface exerts a predominant templating effect on the periodic position of S nanoclusters. It should be noted that interaction between nanoclusters play role in the formation of the superlattice. We additionally conduct an analysis on nearest-neighboring average pair distances ̅ compared with the noninteracting adatom random separation distribution  ( ) for nanocluster superlattice.16 As shown in Fig. S5, the ̅ (~1.3 nm) obviously shift towards larger value compared with the peak position (~0.5 nm) in  ( ). The results show that long-range repulsive interaction exists in-between S nanoclusters, which also stabilizes the superlattice. Thus, it concludes that the substrate templating effect together with repulsive interaction in-between nanoclusters determine the periodicity of the nanocluster superlattice to a large extent. However, these two factors could only induce “square”-like nanocluster superlattice following the symmetry of Cu(100) substrate, which cannot explain the observed anisotropic periodicity of the superlattice (rectangle symmetry). Figure 2g shows a typical STM image around the boundary of the superlattice (indicated by blue lines). On the left region of the boundary, we do not observe the S nanoclusters superlattice. Instead, a typical one-dimensional moiré pattern generated by the lattices of graphene and Cu(100) facet is clearly observed. Interestingly, the one-dimensional moiré pattern on the left region is in line with the S superlattice on the right region (indicated by white dashed lines). This strongly indicates that the periodicity of the nanocluster superlattice is

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

possibly modulated by the moiré pattern of graphene and the underlying Cu facets. In Fig. 2h, we schematically model the three-layer structure of graphene, S superlattice and Cu(100). The graphene and Cu(100) is stacked in alignment (or 90°) and the moiré pattern could be visualized through the model simulation. The short-axis periodicity of the moiré supercell (~1.2 nm) is in coincidence with that of the S nanoclusters superlattice (white dashed lines). This observation implies that the one-dimensional graphene moiré pattern plays an extra templating role in stabilizing more condensed S nanoclusters arrays in parallel direction while no (or weak) templating effect on the vertical direction. In all, the substrate templating effect, repulsive interaction in-between nanoclusters together with graphene moiré pattern templating effect jointly determine the formation of S nanocluster superlattice. In turn, the underlying S nanocluster superlattices could modulate the electronic properties of on-top graphene. Figure 3a shows the tunneling spectra on graphene/nanocluster superlattice regions which is quite different from pristine regions (seen in Fig. S6). The Dirac point shifts upward to -160 meV, and a new dip arises at 350 meV. The dip is assigned as the emergence of new superlattice Dirac points at the conduction band, which are generated by the period potential of S nanocluster superlattice. The absence or faint feature of superlattice Dirac points at the valence band (theoretically at -670 meV) is possibly due to electron-hole asymmetry. The Line profile of tunneling spectra along the nanocluster superlattices (Fig. 3b) demonstrates that the energies of Dirac points and superlattice Dirac points remain constant whatever on and off the nanoclusters, proving the fact the superlattice potential generated by nanoclusters is a global effect. It is worth noting that the spectrum intensities along the nanocluster superlattice show periodic oscillation as indicated by the white dashed curve. Through tunneling spectra mappings (Fig. 3c, d), it could be clearly seen that the regions of nanoclusters always show a lower local density of states (LDOS) compared to around regions. This phenomenon indicates nanoclusters also have local effect on graphene which suppress electron density on the spot. But not like the global effect induced by superlattices, the local effect didn’t change the electronic band structure of

ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

graphene. Figure 4 presents STM characterizations of the second type of S superlattices at the interface. Figure 4a shows the main feature of the sample that a large amount of dispersed atomic protrusions with a density of (5.5±0.2) × 1014 cm-2 could be visualized on graphene lattices. These protrusions were attributed to individual S adatoms trapped at the interface. The S adatoms locally exhibit square or hexagonal order (inset of Fig. 4a). However, there is no long-range order in large area. Such a result is further confirmed by FT of the image (Fig. 4d), which shows scattering-ring feature (red dashed rings) besides the set of clear reciprocal lattice of graphene (black circles). We define the average radius of the rings as 1/ , where  corresponds to the average pair distance of the S adatoms and is estimated to be 0.5±0.15 nm. A relatively high kinetic energy or the short relaxation time of the S adatoms may be reason that limit the atoms forming long-range order in large area. This is especially true for areas with high coverage of S adatoms. For areas with low coverage of S adatoms, it is relatively easy to observe ordered S adatoms superlattices confined at the interface. As shown in Fig. 4b and 4c, ordered S adatoms superlattices with a S density of (2.6±0.2) × 1014 cm-2 and (3.5±0.2) × 1014 cm-2 were clearly observed at the area. The FT of the images (Fig. 4e, f) exhibit clear patterns of both the graphene lattice (black circles) and the S adatoms superlattice (first-order points marked by red circles). The FT patterns of superlattices show the symmetries of stripped hexagonal lattice and square lattice respectively, indicating of the diversity of atomic superlattices in the interface. Furthermore, the positions of the S adatoms have strong correlation with the high-symmetry points of graphene lattice, i.e., most of which are right underneath carbon atoms (red cycles in Fig. 4a inset). It indicates that the graphene may play a notable role of atomic templating in the formation of the S adatoms superlattices. Analysis on nearest-neighboring average pair distances ̅ compared with the noninteracting adatom random separation distribution  ( ) for atomic superlattices were also conducted in these regions (Fig. S7). The data also strongly indicates that long-range repulsive interaction exists in-between atoms in the superlattices.

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

It was proposed that specific-ordered exotic adatoms could tune the electronic properties of graphene into new quantum phases beyond the semimetal state of pristine graphene.50–58 One possible candidate of novel quantum phases in graphene is the Kekulé distortion phase.59–62 Exotic adatoms on graphene could scatter electrons between two non-equivalent Dirac cones (valleys) at the corners of the Brillouin zone (BZ) (Fig. 5a). The so-called intervalley scattering generates the Friedel oscillations (FO)—the electronic ripples in total charge density—in the vicinities of the adatoms. When the arrangement of the adatoms on graphene follows an (√3 × √3)R30° (R3) order (indicated with three different colors in inset of Fig. 5a), the induced FO of individual adatoms could be in a coherent phase and behave as the RKKY-type interaction between the adatoms. The coherence of the FO can significantly strengthen oscillations of LDOS and lead to the emergence of the KD phase in graphene. Fig. 4b shows a special area of herringbone-like morphology. The one-dimensional strips (red dashed lines marked in Fig. 5b) are likely to be the corrugation of underlying substrate. Figure 5c (white box in Fig. 5b) shows clear honeycomb lattices of graphene emerging over the whole area (white hexagons in Fig. 5c). Another important feature is that hexagonal contrasts of graphene lattice show a typical R3-ordered atomic feature (red hexagons in Fig. 5c), which is exactly alike to the predicted KD phase. The simultaneously-acquired dI/dV mapping gives more clearly evidence of the R3-ordered nature of the CDW state (Fig. 5d). Above results indicate that a typical KD phase was observed and the stripped underlay may be the origin of the special symmetry-broken phase of graphene. To elucidate the spatial configuration of graphene and substrate, FT of Fig. 5c is shown in Fig. 5e. Three sets of ordered patterns are marked in the image: the outer one corresponds to the reciprocal lattice of the graphene lattice (Qg, black circles); the intermediate pattern corresponds to the reciprocal lattice of the typical R3 superlattice, which is assigned to the KD phase (Qk, blue circles); the rest inner scattering peaks (red circles) can be seen as a largely stretched hexagonal pattern of underlying “substrate”, which is attributed to the reciprocal pattern of the ordered S adatoms superlattice, as

ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

demonstrated in Fig. 5d. Figure 5f schematically shows the atoms of the “substrate” (red balls) in real space superimposed on graphene lattice by extracting the details of the pattern from the FT. The atoms of the “substrate”, which is attributed to the S adatoms, are commensurate with graphene lattice and can be described by two unequal integer sums of graphene vectors  =   − 5" (~1.13 nm) and # =   + " (~0.43 nm), respectively. By further coloring the graphene lattice into three-color Kekulé texture (Fig. 5f upper right), it is surprisingly to find that the ordered adatoms sit in the same-color grid, which follow a hidden Kekulé order. We employed a phenomenological model of the KD phase to further confirm the proposed configuration (left panel of Fig. 5g). The simulated STM image (right panel of Fig. 5g) reproduces well the experimental topography shown in Fig. 5b and Fig. 5c, demonstrating the validity of the proposed configuration of graphene and adatoms superlattice. This well explains the origin of such strong KD phase in this area. The symmetry of graphene may play an important role in the formation of ordered adatoms superlattices at interface since that the adatoms superlattice to a large extent follows the order of graphene lattice. The KD phase is expected to possess a gapped band structure which depends on strength of intervalley scattering induced by adatoms. As characterized by our STS measurements (Fig. 5h), the spectra with different tip-sample distances (tuned by setpoints) show similar asymmetric bandgap of 245±5 meV ranging from -100 meV (top of the valence band) to 145 meV (bottom of the conduction band), which is totally different from the gapless feature of pristine graphene. Bandgap opening of graphene is of central importance for electronics applications of graphene. There are three general scenarios for the bandgap opening: (i) bandgap induced by chemical absorption; (ii) bandgap resulted from local strain;63,64 (iii) bandgap generated by strong interaction between graphene and substrate. We can easily rule out scenarios (i) and (ii) because there is neither chemical absorption nor local strain according to our STM measurements (strain analysis seen in Fig. S8). In scenario (iii), both experiment and theory demonstrated explicitly that there is no

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

observable bandgap for graphene on Cu substrate due to their weak interactions.65–68 Considering the observation of the unique R3 CDW morphology and special-ordered S superlattice in this region, it is reasonable to conclude that the bandgap opening originates from the KD of graphene. The origin of the gap in the KD phase was attributed to inequivalent electrons hopping t1 and t2 between the nearest sublattices (inset in Fig. 5h). According to the simulation, a similar gap would open at Γ point considering ~5% difference between t1 and t2, creating massive Dirac fermions with effective mass of mD= 0.28±0.02 me (see supplementary method and Fig. S9 for details of calculation). Compared with previous work, the observation of bandgap opening here may be attributed to two reasons: one is that a more ordered atomic array was formed under graphene compared to a hidden-Kekulé order, which profoundly enhance the coherence of FO; the other one is that the S atoms may have a stronger coupling with graphene lattices compared to Cu vacancies, which may strengthen the FO in the vicinity of atoms.

CONCLUSION In conclusion, our experiments describe a method to create ordered S nanoclusters superlattice and adatoms superlattice through self-organization at the interface of graphene and Cu substrate. Both the symmetries of graphene and Cu substrate affect the self-organization process, leading to various self-organized nanostructures. The superlattice with specific periodicity in turn modulate the electronic properties of graphene into the KD phase. The reported self-organization of atomic superlattices at interfaces could be extended to other systems, which may provide a route to realize exotic electronic states in graphene and other two-dimensional materials.

METHODS CVD preparation of graphene. The 25-µm Cu foil was purchased from Alfa Aesar. Before growth, Cu foil was first electropolished at 1.5 V DC voltage for 60 min, using

ACS Paragon Plus Environment

Page 10 of 22

Page 11 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

a mixture of phosphoric acid and ethylene glycol (volume ratio = 3:1) as the electrolyte. The pre-treated Cu foil was loaded into a 2-inch quartz tube of low-pressure chemical vapor deposition furnace for sample growth. The Cu foil was first heated from room temperature to 1,000 °C in 30 min and kept for another 30 min, with 50 sccm (standard cubic centimeter per minute) H2 and 50 sccm Ar as carrier gas. In the second step, the furnace temperature was set to 1000°C and then 5 sccm CH4 was introduced into the furnace to grow graphene for 10 min. After all growth, the sample was cooled down to room temperature slowly (~20 ℃/min). STM and STS measurements. STM/STS characterizations were performed in ultrahigh vacuum scanning probe microscopes (USM-1500 and USM-1400) from UNISOKU. The STM tips were obtained by chemical etching from a wire of Pt/Ir (80/20%) alloys. Lateral dimensions observed in the STM images were calibrated using a standard graphene lattice and a Si (111)-(7×7) lattice and Ag (111) surface. The dI/dV measurements were taken with a standard lock-in technique by turning off the feedback circuit and using a 793-Hz 5mV A.C. modulation of the sample voltage.

SUPPORTING INFORMATION The supporting information involves DFT calculations, tight-binding calculations, XPS data, supplementary STM images and strain analysis et al. The supporting information is available free of charge on the ACS Publications website at DOI:

AUTHOR INFORMATION * Corresponding Author E-mail: [email protected] § Donglin Ma and Zhongqiu Fu contributed equally to this work. The authors declare no competing financial interest.

ACKNOWLEDGEMENTS This work was supported by the National Natural Science Foundation of China (Grant ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nos. 11674029, 11422430, 11374035), the National Basic Research Program of China (Grants Nos. 2014CB920903, 2013CBA01603), the program for New Century Excellent Talents in University of the Ministry of Education of China (Grant No. NCET-13-0054), the China Postdoctoral Science Foundation (No. 212400207). L.H. also acknowledges support from the National Program for Support of Top-notch Young Professionals, support from “the Fundamental Research Funds for the Central Universities”, and support from “Chang Jiang Scholars Program”.

REFERENCES (1) (2) (3) (4) (5)

(6)

(7) (8) (9) (10) (11)

(12) (13)

Barth, J. V.; Costantini, G.; Kern, K. Engineering Atomic and Molecular Nanostructures at Surfaces. Nature 2005, 437, 671–679. Bohr, M. T. Nanotechnology Goals and Challenges for Electronic Applications. IEEE Trans. Nanotechnol. 2002, 1, 56–62. Thompson, D. A.; Best, J. S. The Future of Magnetic Data Storage Techology. IBM J. Res. Dev. 2000, 44, 311–322. Binnig, G.; Rohrer, H.; Gerber, C.; Weibel, E. Tunneling through a Controllable Vacuum Gap. Appl. Phys. Lett. 1982, 40, 178–180. Binnig, G.; Rohrer, H.; Gerber, C.; Weibel, E. Surface Studies by Scanning Tunneling Microscopy. Phys. Rev. Lett. 1982, 49, 57–61. Binnig, G.; Rohrer, H. Scanning Tunneling Microscopy—from Birth to Adolescence. Rev. Mod. Phys. 1987, 59, 615–625. Becker, R. S.; Golovchenko, J. A.; Swartzentruber, B. S. Atomic-Scale Surface Modifications Using a Tunnelling Microscope. Nature 1987, 325, 419–421. Foster, J. S.; Frommer, J. E.; Arnett, P. C. Molecular Manipulation Using a Tunnelling Microscope. Nature 1988, 331, 324–326. Crommie, M. F.; Lutz, C. P.; Eigler, D. M. Confinement of Electrons to Quantum Corrals on a Metal Surface. Science 1993, 262, 218–220. Gomes, K. K.; Mar, W.; Ko, W.; Guinea, F.; Manoharan, H. C. Designer Dirac Fermions and Topological Phases in Molecular Graphene. Nature 2012, 483, 306–310. Slot, M. R.; Gardenier, T. S.; Jacobse, P. H.; van Miert, G. C. P.; Kempkes, S. N.; Zevenhuizen, S. J. M.; Smith, C. M.; Vanmaekelbergh, D.; Swart, I. Experimental Realization and Characterization of an Electronic Lieb Lattice. Nat. Phys. 2017, 13, 672– 676. Eigler, D. M.; Schweizer, E. K. Positioning Single Atoms with a Scanning Tunnelling Microscope. Nature 1990, 344, 524–526. Song, C.-L.; Sun, B.; Wang, Y.-L.; Jiang, Y.-P.; Wang, L.; He, K.; Chen, X.; Zhang, P.; Ma,

ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

(14)

(15)

(16)

(17)

(18) (19) (20)

(21) (22) (23)

(24)

(25) (26) (27) (28)

(29)

(30)

X.-C.; Xue, Q.-K. Charge-Transfer-Induced Cesium Superlattices on Graphene. Phys. Rev. Lett. 2012, 108, 156803. Silly, F.; Pivetta, M.; Ternes, M.; Patthey, F.; Pelz, J. P.; Schneider, W.-D. Coverage-Dependent Self-Organization: From Individual Adatoms to Adatom Superlattices. New J. Phys. 2004, 6, 16. Yokoyama, T.; Takahashi, T.; Shinozaki, K.; Okamoto, M. Quantitative Analysis of Long-Range Interactions between Adsorbed Dipolar Molecules on Cu(111). Phys. Rev. Lett. 2007, 98, 206102. Knorr, N.; Brune, H.; Epple, M.; Hirstein, A.; Schneider, M. A.; Kern, K. Long-Range Adsorbate Interactions Mediated by a Two-Dimensional Electron Gas. Phys. Rev. B 2002, 65, 115420. Silly, F.; Pivetta, M.; Ternes, M.; Patthey, F.; Pelz, J. P.; Schneider, W.-D. Creation of an Atomic Superlattice by Immersing Metallic Adatoms in a Two-Dimensional Electron Sea. Phys. Rev. Lett. 2004, 92, 016101. Brune, H.; Giovannini, M.; Bromann, K.; Kern, K. Self-Organized Growth of Nanostructure Arrays on Strain-Relief Patterns. Nature 1998, 394, 451–453. Lu, W.; Sastry, A. M. Self-Assembly for Semiconductor Industry. IEEE Trans. Semicond. Manuf. 2007, 20, 421–431. Narayanan Nair, M.; Cranney, M.; Jiang, T.; Hajjar-Garreau, S.; Aubel, D.; Vonau, F.; Florentin, A.; Denys, E.; Bocquet, M.-L.; Simon, L. Noble-Metal Intercalation Process Leading to a Protected Adatom in a Graphene Hollow Site. Phys. Rev. B 2016, 94, 075427. Dresselhaus, M. S.; Dresselhaus, G. Intercalation Compounds of Graphite. Adv. Phys. 1981, 30, 139–326. Mu, R.; Fu, Q.; Jin, L.; Yu, L.; Fang, G.; Tan, D.; Bao, X. Visualizing Chemical Reactions Confined under Graphene. Angew. Chem. Int. Ed. 2012, 51, 4856–4859. Watcharinyanon, S.; Johansson, L. I.; Xia, C.; Virojanadara, C. Changes in Structural and Electronic Properties of Graphene Grown on 6H-SiC(0001) Induced by Na Deposition. J. Appl. Phys. 2012, 111, 083711. Sandin, A.; Jayasekera, T.; Rowe, J. E.; Kim, K. W.; Buongiorno Nardelli, M.; Dougherty, D. B. Multiple Coexisting Intercalation Structures of Sodium in Epitaxial Graphene-SiC Interfaces. Phys. Rev. B 2012, 85, 125410. Virojanadara, C.; Watcharinyanon, S.; Zakharov, A. A.; Johansson, L. I. Epitaxial Graphene on 6H-SiC and Li Intercalation. Phys. Rev. B 2010, 82, 205402. de Lima, L. H.; Landers, R.; de Siervo, A. Patterning Quasi-Periodic Co 2D-Clusters underneath Graphene on SiC(0001). Chem. Mater. 2014, 26, 4172–4177. Sutter, P.; Sadowski, J. T.; Sutter, E. A. Chemistry under Cover: Tuning Metal−Graphene Interaction by Reactive Intercalation. J. Am. Chem. Soc. 2010, 132, 8175–8179. Dong, A.; Fu, Q.; Wei, M.; Liu, Y.; Ning, Y.; Yang, F.; Bluhm, H.; Bao, X. Facile Oxygen Intercalation between Full Layer Graphene and Ru(0001) under Ambient Conditions. Surf. Sci. 2015, 634, 37–43. Larciprete, R.; Ulstrup, S.; Lacovig, P.; Dalmiglio, M.; Bianchi, M.; Mazzola, F.; Hornekær, L.; Orlando, F.; Baraldi, A.; Hofmann, P.; Lizzit, S. Oxygen Switching of the Epitaxial Graphene–Metal Interaction. ACS Nano 2012, 6, 9551–9558. Grånäs, E.; Knudsen, J.; Schröder, U. A.; Gerber, T.; Busse, C.; Arman, M. A.; Schulte, K.;

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(31)

(32)

(33)

(34) (35)

(36) (37)

(38)

(39)

(40)

(41) (42) (43) (44) (45) (46)

Andersen, J. N.; Michely, T. Oxygen Intercalation under Graphene on Ir(111): Energetics, Kinetics, and the Role of Graphene Edges. ACS Nano 2012, 6, 9951–9963. Gierz, I.; Suzuki, T.; Weitz, R. T.; Lee, D. S.; Krauss, B.; Riedl, C.; Starke, U.; Höchst, H.; Smet, J. H.; Ast, C. R.; Kern, K. Electronic Decoupling of an Epitaxial Graphene Monolayer by Gold Intercalation. Phys. Rev. B 2010, 81, 235408. Romero-Muñiz, C.; Martín-Recio, A.; Pou, P.; Gómez-Rodríguez, J. M.; Pérez, R. Strong Dependence of Flattening and Decoupling of Graphene on Metals on the Local Distribution of Intercalated Oxygen Atoms. Carbon 2016, 101, 129–134. 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 2011, 84, 125423. 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, 246804. Schumacher, S.; Wehling, T. O.; Lazić, P.; Runte, S.; Förster, D. F.; Busse, C.; Petrović, M.; Kralj, M.; Blügel, S.; Atodiresei, N.; Caciuc, V.; Michely, T. The Backside of Graphene: Manipulating Adsorption by Intercalation. Nano Lett. 2013, 13, 5013–5019. 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, 086111. Schumacher, S.; Huttmann, F.; Petrović, M.; Witt, C.; Förster, D. F.; Vo-Van, C.; Coraux, J.; Martínez-Galera, A. J.; Sessi, V.; Vergara, I.; Rückamp, R; Grüninger, M.; Schleheck, N.; Meyer zu Heringdorf, F.; Ohresser, R.; Kralj, M.; Wehling, T. O.; Michely, T. Europium underneath Graphene on Ir(111): Intercalation Mechanism, Magnetism, and Band Structure. Phys. Rev. B 2014, 90, 235437. Petrović, M.; Šrut Rakić, I.; Runte, S.; Busse, C.; Sadowski, J. T.; Lazić, P.; Pletikosić, I.; Pan, Z.-H.; Milun, M.; Pervan, P.; Atodiresei, N.; Brako. R.; Šokčević, D.; Valla, T.; Michely, T.; Kralj, M. The Mechanism of Caesium Intercalation of Graphene. Nat. Commun. 2013, 4, 2772. Oida, S.; McFeely, F. R.; Hannon, J. B.; Tromp, R. M.; Copel, M.; Chen, Z.; Sun, Y.; Farmer, D. B.; Yurkas, J. Decoupling Graphene from SiC(0001) via Oxidation. Phys. Rev. B 2010, 82, 041411. Robinson, J. A.; Hollander, M.; LaBella, M.; Trumbull, K. A.; Cavalero, R.; Snyder, D. W. Epitaxial Graphene Transistors: Enhancing Performance via Hydrogen Intercalation. Nano Lett. 2011, 11, 3875–3880. Fu, Q.; Bao, X. Surface Chemistry and Catalysis Confined under Two-Dimensional Materials. Chem. Soc. Rev. 2017, 46, 1842–1874. Colaianni, M. L.; Chorkendorff, I. Scanning-Tunneling-Microscopy Studies of the S-Induced Reconstruction of Cu(100). Phys. Rev. B 1994, 50, 8798–8806. Crommie, M. F.; Lutz, C. P.; Eigler, D. M. Imaging Standing Waves in a Two-Dimensional Electron Gas. Nature 1993, 363, 524–527. Hinch, B. J.; Frenken, J. W. M.; Zhang, G.; Toennies, J. P. Sulfur Adatom Diffusion on the Cu(111) Surface. Surf. Sci. 1991, 259, 288–300. Colaianni, M. L.; Chorkendorff, I. Scanning-Tunneling-Microscopy Studies of the S-Induced Reconstruction of Cu(100). Phys. Rev. B 1994, 50, 8798–8806. Walen, H.; Liu, D.-J.; Oh, J.; Yang, H. J.; Spurgeon, P. M.; Kim, Y.; Thiel, P. A. Sulfur

ACS Paragon Plus Environment

Page 14 of 22

Page 15 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

(47)

(48)

(49) (50)

(51) (52)

(53)

(54) (55) (56) (57) (58)

(59) (60) (61) (62)

(63) (64)

Atoms Adsorbed on Cu(100) at Low Coverage: Characterization and Stability against Complexation. J. Phys. Chem. B 2018, 122, 963–971. Zeng, H. C.; McFarlane, R. A.; Mitchell, K. A. R. Investigation with Low-Energy Electron Diffraction of the Adsorbate-Induced Metal Relaxations in the Cu(100)-(2×2)-S Surface Structure. Phys. Rev. B 1989, 39, 8000–8002. Barton, J. J.; Bahr, C. C.; Hussain, Z.; Robey, S. W.; Tobin, J. G.; Klebanoff, L. E.; Shirley, D. A. Direct Surface Structure Determination with Photoelectron Diffraction. Phys. Rev. Lett. 1983, 51, 272–275. Zeng, H. C.; Sodhi, R. N. S.; Mitchell, K. A. R. A Leed Crystallographic Analysis for the Cu(100)-(2×2)-S Surface Structure. Surf. Sci. 1986, 177, 329–337. Grushin, A. G.; Castro, E. V.; Cortijo, A.; de Juan, F.; Vozmediano, M. A. H.; Valenzuela, B. Charge Instabilities and Topological Phases in the Extended Hubbard Model on the Honeycomb Lattice with Enlarged Unit Cell. Phys. Rev. B 2013, 87, 085136. Weeks, C.; Franz, M. Interaction-Driven Instabilities of a Dirac Semimetal. Phys. Rev. B 2010, 81, 085105. García-Martínez, N. A.; Grushin, A. G.; Neupert, T.; Valenzuela, B.; Castro, E. V. Interaction-Driven Phases in the Half-Filled Spinless Honeycomb Lattice from Exact Diagonalization. Phys. Rev. B 2013, 88, 245123. 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. Hou, C.-Y.; Chamon, C.; Mudry, C. Electron Fractionalization in Two-Dimensional Graphenelike Structures. Phys. Rev. Lett. 2007, 98, 186809. Capponi, S. Phase Diagram of Interacting Spinless Fermions on the Honeycomb Lattice. J. Phys. Condens. Matter 2017, 29, 043002. Min, H.; Hill, J. E.; Sinitsyn, N. A.; Sahu, B. R.; Kleinman, L.; MacDonald, A. H. Intrinsic and Rashba Spin-Orbit Interactions in Graphene Sheets. Phys. Rev. B 2006, 74. Nandkishore, R.; Levitov, L. S.; Chubukov, A. V. Chiral Superconductivity from Repulsive Interactions in Doped Graphene. Nat. Phys. 2012, 8, 158–163. Ding, J.; Qiao, Z.; Feng, W.; Yao, Y.; Niu, Q. Engineering Quantum Anomalous/Valley Hall States in Graphene via Metal-Atom Adsorption: An Ab-Initio Study. Phys. Rev. B 2011, 84, 195444. Cheianov, V. V.; Fal’ko, V. I.; Syljuåsen, O.; Altshuler, B. L. Hidden Kekulé Ordering of Adatoms on Graphene. Solid State Commun. 2009, 149, 1499–1501. Cheianov, V. V.; Syljuåsen, O.; Altshuler, B. L.; Fal’ko, V. Ordered States of Adatoms on Graphene. Phys. Rev. B 2009, 80, 233409. Castro Neto, A. H.; Kotov, V. N.; Nilsson, J.; Pereira, V. M.; Peres, N. M. R.; Uchoa, B. Adatoms in Graphene. Solid State Commun. 2009, 149, 1094–1100. Gutiérrez, C.; Kim, C.-J.; Brown, L.; Schiros, T.; Nordlund, D.; Lochocki, E. B.; Shen, K. M.; Park, J.; Pasupathy, A. N. Imaging Chiral Symmetry Breaking from Kekulé Bond Order in Graphene. Nat. Phys. 2016, 12, 950–958. Guinea, F.; Katsnelson, M. I.; Geim, A. K. Energy Gaps and a Zero-Field Quantum Hall Effect in Graphene by Strain Engineering. Nat. Phys. 2010, 6, 30–33. Chen, C.; Avila, J.; Arezki, H.; Nguyen, V. L.; Shen, J.; Mucha-Kruczyński, M.; Yao, F.;

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(65) (66)

(67)

(68)

Boutchich, M.; Chen, Y.; Lee, Y. H.; Asensio, M. C. Large Local Lattice Expansion in Graphene Adlayers Grown on Copper. Nat. Mater. 2018, 17, 450–455. Giovannetti, G.; Khomyakov, P. A.; Brocks, G.; Karpan, V. M.; van den Brink, J.; Kelly, P. J. Doping Graphene with Metal Contacts. Phys. Rev. Lett. 2008, 101, 026803. Khomyakov, P. A.; Giovannetti, G.; Rusu, P. C.; Brocks, G.; van den Brink, J.; Kelly, P. J. First-Principles Study of the Interaction and Charge Transfer between Graphene and Metals. Phys. Rev. B 2009, 79, 195425. Gutiérrez, C.; Brown, L.; Kim, C.-J.; Park, J.; Pasupathy, A. N. Klein Tunnelling and Electron Trapping in Nanometre-Scale Graphene Quantum Dots. Nat. Phys. 2016, 12, 1069– 1075. Bai, K.-K.; Zhou, J.-J.; Wei, Y.-C.; Qiao, J.-B.; Liu, Y.-W.; Liu, H.-W.; Jiang, H.; He, L. Generating Atomically Sharp p-n Junctions in Graphene and Testing Quantum Electron Optics on the Nanoscale. Phys. Rev. B 2018, 97, 045413.

ACS Paragon Plus Environment

Page 16 of 22

Page 17 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Figure 1. A general synthetic method for self-organized growth of sulfur superlattices at interface between graphene and substrates. In the first step, the growth substrates were annealed at high temperature to dissolve and activate S atoms in the bulk. In the following, carbon sources were introduced into the system for graphene growth, then the samples were slowly cooled to room temperature. During the slow cooling process, vast of S atoms begun to segregate from the bulk onto the metal surfaces due to the decrease of solubility. At this stage, the S adatoms had sufficient energy and time to rearrange beneath the as-grown graphene sheet to form kinds of ordered assemblies.

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Self-organization of S nanoclusters superlattice at interface. (a) STM topography of large-area square-ordered superlattices (Vb = 800 mV, I = 400 pA). (b) Zoom-in image clearly exhibiting honeycomb lattices of graphene (Vb = -450 mV, I = 200 pA) Inset: Fourier transform (FT) of the large-version STM image of (b).. (c) Atomic-resolved STM image of the superlattice under graphene (Vb = 50 mV, I = 150 pA). The nanoclusters are the ordered assembly of adatoms in a special “dice-five” shape, as illustrated by red circles in white box. (d) S5 nanocluster on Cu(100) as determined by DFT calculations. Sulfur atoms are indicated by yellow balls, while copper atoms are the orange balls. (e) Simulated STM image of the S5 nanocluster on Cu(100). (f) Typical STM morphology (Vb = 150 mV, I = 310 pA) around the edges of the nanocluster superlattices shows S nanocluster “vacancies” (red squares) at the edges. (g) A boundary of the nanoclusters superlattice (Vb = 50 mV, I = 200 pA). One-dimensional moiré patterns are indicated by red dashed lines. The missing-row

ACS Paragon Plus Environment

Page 18 of 22

Page 19 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

defects and adding-row defects of the superlattice around the boundary were marked in white dashed oval and rectangle, respectively. (h) Schematics of the three layer configuration of graphene, S superlattice and Cu(100) substrate. The moiré pattern generated by the misalignment of graphene and Cu(100) is illustrated by white dashed rectangle. The unit cell of S superlattice is illustrated by red dashed rectangle. The scale bars are 6 nm in (a), 1 nm in (b), 4 nm-1 in inset of (b), 3 nm in (c, f, g).

Figure 3. The electronic properties of graphene modulated by the S nanoclusters superlattice. (a) Tunneling spectra on nanocluster region (black) and off that (red) respectively. (b) Line profile of tunneling spectra along the nanocluster superlattices. (c, d) Representative STS mapping of nanocluster superlattices at 91 meV and -111 meV. The scale bars are 2 nm in (c, d).

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Self-organization of S adatoms superlattice at the interface. (a) STM topography of large-area graphene lattices with a large amount of dispersed adatoms underneath graphene (Vb = 200 mV, I = 250 pA). Inset shows a small area where the S adatoms locally exhibit a weak square order. (b, c) Large-area ordered S adatoms superlattice with clear graphene lattice (Vb = 650 mV, I = 250 pA; Vb = -500 mV, I = 200 pA). (d, e, f) Corresponding FT images of (a, b, c). The scale bars are 4 nm in (a), 3 nm in (b), 2 nm in (c), 4 nm-1 in (d, e, f).

Figure 5. Realization of Kekulé distortion (KD) phase on graphene modulated by the S adatoms superlattice. (a) Brillouin zone and real-space schematics (upper left) of graphene in the KD phase. (b) STM topographic image of strong herringbone-like charge-density wave (CDW) morphologies (Vb = 300 mV, I = 200 pA). (c) Low-bias

ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

image (white box area in (a)) showing honeycomb lattice of graphene and R3-ordered atomic features, as marked by white and red hexagons (Vb = 30 mV, I = 200 pA). (d) The simultaneously-acquired dI/dV mapping (Vb = 218 mV, I = 200 pA). (e) Fourier transform of (c). (f) Schematic illustrations of underneath adatoms superimposed on color-coded graphene lattice. (g) Phenomenological model of the KD phase on the configuration of (f). (h) Scanning tunnelling spectra on the KD phase with different tunneling current setpoints (from 200 to 1000 pA). Edges of valence band and conduction band are indicated by arrows. Inset is the honeycomb model for calculation, with inequivalent electrons hoppings t1and t2 between nearest sublattices. The scale bars are 2 nm in (b), 1 nm in (c, d), 4 nm-1 in (e).

ACS Paragon Plus Environment

ACS Nano 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

82x44mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 22