## Hexagonal Nanopits with the Zigzag Edge State on Graphite Surfaces Synthesized by Hydrogen-Plasma Etching

T. Matsui , 1 ∗ H. Sato , K. Kita , A. E. B. Amend , and Hiroshi Fukuyama 1 1 1 1 2 , † 1 Department of Physics, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan and 2 Cryogenic Research Center, The University of Tokyo, 2-11-16 Yayoi, Bunkyo-ku, Tokyo 113-0032, Japan (Dated: August 2, 2019)

We studied, by scanning tunneling microscopy, the morphology of nanopits of monolayer depth created at graphite surfaces by hydrogen plasma etching under various conditions such as H 2 pressure, temperature, etching time, and RF power of the plasma generation. In addition to the known pressure-induced transition of the nanopit morphology, we found a sharp temperature-induced transition from many small rather round nanopits of ∼ 150 nm size to few large hexagonal ones of 300-600 nm within a narrow temperature range. The remote and direct plasma modes switching mechanism, which was proposed to explain the pressure-induced transition, is not directly applicable to this newly found transition. Scanning tunneling spectroscopy (STS) measurements of edges of the hexagonal nanopits fabricated at graphite surfaces by this method show clear signatures of the peculiar electronic state localized at the zigzag edge (edge state), i.e. , a prominent peak near the Fermi energy accompanied by suppressions on either side in the local density of states. These observations indicate that the hexagonal nanopits consist of a high density of zigzag edges. The STS data also revealed a domain structure of the edge state in which the electronic state varies over a length scale of ∼ 3 nm along the edge. The present study will pave the way for microscopic understanding of the anisotropic etching mechanism and of spin polarization in zigzag nanoribbons which are promising key elements for future graphene nanoelectronics.

## I. INTRODUCTION

Reflecting the honeycomb lattice structure, graphene has two types of edges. One is the zigzag edge and the other the armchair one. Among various novel physical properties of graphene [1], the peculiar electronic state localized at the zigzag edge, the so called 'zigzag edge state' or simply 'edge state', is one of the most attractive ones. Such a localized state emerges due to symmetry breaking of the two equivalent sublattices in the honeycomb structure at the zigzag edge. Inspired by the prediction by Fujita and co-workers [2, 3], the edge state was experimentally found by scanning tunneling microscopy and spectroscopy (STM/S) measurements at naturally existing monatomic step-edges at highly oriented pyrolytic graphite (HOPG) surfaces [4, 5] and at similar step-edges in nanographites synthesized on HOPG [6]. The decay length of the edge state was measured as 1.2 ± 0.2 nm [5].

applied to carbon-based switching nanodevices because the system is expected to be semiconducting with a band gap in the AF configuration while it is metallic in the FM one [10, 11].

In principle, the zigzag edge state can be spin polarized even under infinitesimally small electron-electron interaction due to its flat band nature located near the Fermi energy ( E F ) [2]. In graphene nanoribbons with two parallel zigzag edges (zGNRs), the edge spins are predicted to order ferromagnetically (FM) along the edge and antiferromagnetically (AF) between the edges [2, 7-9]. By applying a magnetic field [10] or an in-plane electric field [11], the AF configuration will flip to the FM one. This can be

∗ matsui@phys.s.u-tokyo.ac.jp

† hiroshi@phys.s.u-tokyo.ac.jp

In order to experimentally test the theoretical predictions, it is crucial to realize nearly pure zigzag edges, because the spin polarized states become less stable against the non-magnetic state under edge disorders [1215]. After the first observations of the edge state [46], numerous trials have been made to shape graphene edges on an atomic scale in order to fabricate zGNRs with the spin-polarized edge state. They include reactive ion etching [16-18], anisotropic etching using metallic nanoparticles as catalysts [19-24], carbothermal etching [25-27], nanolithography with STM [28] and atomic force microscope (AFM) tips [29], unzipping of carbon nanotubes [30-34], and bottom-up fabrications [35-38]. However, non of them has succeeded in synthesizing zigzag edges which satisfy all the necessary conditions: (i) quality or purity (low concentration of the armchair edge), (ii) known edge termination, and (iii) scalability or reproducibility.

Among the previous attempts, the hydrogen (H) plasma etching has apparent advantages from the view point of the requirements (ii) and (iii). It is reported that this method can create hexagonal nanopits on graphite surfaces [39-41] with edges in parallel to the zigzag direction. In addition, the nanopits can be created in arbitrary spatial arrangements, with which one can synthesize zGNR at high densities [39, 42-44]. However, so far, no atomic-scale spectroscopy has been applied on the edges prepared by this method. Thus it is not clear if it

meets the requirement (i) so as to be used in 'edge-state devices'. In addition, there are conflicts among previous experimental results on monolayer graphene exfoliated on SiO . 2 Yang et al. [39] reported successful creation of hexagonal nanopits of uniform size, while Diankov et al. [40] and Hug et al. [41] reported circular ones. For a zGNR unzipped from a carbon nanotube, no nanopit creation by H-plasma etching has been found [32].

Recently, a new insight on H-plasma etching of graphite has been brought by Hug et al. [41]. They found that anisotropic etching, that grows hexagonal nanopits, proceeds when a graphite sample is located outside of a glow regime of the plasma; the 'remote' plasma mode. On the other hand, when a sample is located inside the glow, defect creation is promoted, resulting in multi-layer terraced nanopits with irregularly shaped edges; the 'direct' plasma mode. They demonstrated that the direct/remote mode can be switched by tuning the gas pressure. However, the role of etching temperature and other parameters are not known.

In this article, we report on comprehensive studies of the morphology of nanopits created at graphite surfaces by H-plasma etching. The shape, depth, and growth rate of the nanopits are examined on an atomic scale with STM for samples etched under various conditions. We found that the morphology suddenly changes from many small nanopits to few large hexagonal nanopits when the reaction temperature is increased from 450 to 500 ◦ C. Such a transition has not been reported before. It is also found that the morphology changes from many irregularly-shaped multi-layer terraced nanopits to few large hexagonal nanopits on increasing the hydrogen gas pressure from 13 to 60 Pa. This pressure variation is consistent with the previous report [41]. Most importantly, our STS measurements have verified that the edges of the hexagonal nanopits prepared by this method are zigzag edges of high purity with the most pronounced edge-state features ever reported. Spatial variations of the edge state studied by STS are also discussed.

## II. RESULTS AND DISCUSSION

## A. Pressure Dependence

Let us first show the H 2 gas pressure ( P ) dependence of H-plasma etching, where we obtained results that agree very well with those by Hug et al. [41]. Figures 1a-h are typical topographic STM images of graphite surfaces etched at several different P from 13 to 590 Pa, where etching temperature ( T ), etching time ( ), and RF power t to generate plasma ( W RF ) are fixed at 500 ◦ C, 40 min , . and 20 W, respectively. At the lowest pressure of 13 Pa, the surface is fully covered by irregularly-shaped nanopits with many steps down to the 10th layer from the top as shown in Figures 1a and b, where 1b is a magnified image of a region in 1a. The cross-sectional profile along the line indicated in Figure 1b is shown in Figure 1i, and each terrace height is replotted as a function of the step number in Figure 1j. From this, the averaged step height is 0 369 . ± 0 003 nm, which is 10 % larger than the interlayer . distance of graphite (0.335 nm). Possible intercalation of atomic H, one of the components of H plasma during the etching process, and its subsequent deintercalation in vacuum may expand the interlayer distance [45, 46]. Clearly, the defect formation is very active here.

As P increases to 60 Pa, the surface morphology suddenly changes, i.e. , the number of steps decreases rapidly, and the nanopit shape changes closer to hexagonal (Figure 1c). At a slightly higher P (= 110 Pa), all the nanopits become hexagonal and larger in size (Figure 1d) indicating that the defect formation is suppressed and the anisotropic lateral etching becomes dominant. Such tendencies steadily proceed up to P = 190 Pa (Figures 1ef) where only few very large hexagonal nanopits (500600 nm) grow. It is not easy to judge if the same tendencies still continue at P ≥ 150 Pa, partly because we cannot exclude the possibility that the top most layer has already been etched away, and partly because a new structure, which consists of deep linear-step-edges, appears as can be seen by the dark regions on the left and right sides of Figures 1g and h, respectively. The height of such step-edges varies randomly from monatomic to 40 layers, and their crystallographic direction is also randomly oriented with respect to the honeycomb lattice of graphite.

We analyzed these STM images quantitatively by evaluating relative terrace areas of different layers. The results are shown in Figures 1k and l where S n is an areal fraction of the n -th layer from the top and D max is the diameter, for irregularly-shaped nanopits, or the diagonal length, for hexagonal ones, of the largest nanopit. In cases where hexagonal nanopits substantially overlap, D max was estimated from the length of one of six sides of the hexagon, assuming its regular shape. Here, 0.22, 3.8, 12, 5.0, 17, 20 and 5 µ m 2 areas were surveyed for samples treated at P = 13, 60, 110, 150, 190, 350 and 590 Pa, respectively. In these figures, it is clear that the etching nature changes from irregular to anisotropic with increasing pressure within a narrow range between 60 and 110 Pa. This agrees very well with the result by Hug et al. [41] who observed a similar change of etching character between 40 and 70 Pa and attributed it to the direct/remote plasma mode change.

We have checked the applicability of the direct/remote plasma scenario to the present data. With our experimental setup, when the plasma glow edge is located inside the furnace, it is not possible to measure a distance ( l g ) between the glow edge and the end of the RF coil along the quartz tube. Thus, l g in this regime was estimated from the known relation, l g ∝ 1 / √ P [41], whose proportionality constant was determined by measuring P and l g when the glow edge is located outside the furnace on the downstream side. The estimation shows that l g reaches the sample position ( l s = 0 36 m) when . P ≥ 60 Pa, which is consistent with the direct/remote

FIG. 1. STM images of graphite surfaces etched by H-plasma at different P of (a)(b) 13, (c) 60, (d) 110, (e) 150, (f) 190, (g) 350, and (h) 590 Pa. The other parameters are fixed at T = 500 ◦ C, t = 40 min., and W RF = 20 W. (b) Magnified image of the region indicated by the square in (a). (i) Cross sectional profile along the solid line in (b), and (j) relation between the terrace height and the step number showing that the nanopit steps are of monolayer height. Note that the heights of the linear steps in (g) and (h) are much larger, i.e. , 10 and 4 layers height, respectively. (k) S n and (l) D max are plotted as a functions of P . It is clear that the etching character changes rapidly between 13 and 60 Pa from irregular to anisotropic as P increases, which is consistent with the observation in Ref. 41. However, in our case, the size distribution of the nanopits is much wider and D max is much larger.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-b1c62453-9901-404c-87e3-ee013e7cba47.jpg)

scenario.

The glowing region should roughly represent the spatial extension of H-ions such as H + 1 , H , and H + 2 + 3 . The life time of the H-ions shortens with increasing gas pressure because of more frequent collisions, resulting in a shorter l g . On the other hand, H-radicals created by ionization and recombination processes survive beyond the glowing region. Therefore, one can expect that the plasma species responsible for the defect creation and the isotropic etching (anisotropic etching) in the direct (remote) regime is H-ions (H-radicals) [41]. According to the measurement of H-plasma components as a function of pressure [47, 48], H + 3 seems to be the main contributor over the other H-ions.

It should be noted, however, that there is a quantitative difference between our result and the previous one by Hug et al. [41]. We found a large nanopit size distribution even in the remote plasma regime at P ≥ 110 Pa suggesting that the defect-creation is active over a wide P range, while they reported only hexagonal nanopits of uniform size suggesting that the nanopits grow only from preexisting defects. This difference may come from a slightly larger oxygen contamination in our gas flow system because oxygen is known to be strongly active in creating surface defects on graphite [49].

## B. Temperature Dependence

Next, we examined the temperature dependence of the H-plasma etching with fixed parameters of P = 110 Pa, t = 40 min., and W RF = 20 W. Surprisingly, we found an abrupt change of the etching character within a narrow temperature range between 450 and 500 ◦ C, as will be described below. Figures 2a-f show STM images of graphite surfaces etched at different T from 200 to 700 ◦ C. Here, the number of surface steps remains less than three layers at any T , since the etching conditions belong to the remote plasma regime in which heavy etching like Figure 1a never occurs.

4

FIG. 2. STM images of graphite surfaces etched by H-plasma at different T of (a) 200, (b) 300, (c) 400, (d) 500, (e) 600, and (f) 700 ◦ C. The other etching parameters are fixed at P = 110 Pa, t = 40 min , and . W RF = 20 W. (g) S n and (h) D max are plotted as a function of T . The dashed line in (h) is the result by Yang et al. [39] obtained with 30 Pa, 40 min., and 20 W. Etching characteristics drastically change between 450 and 500 ◦ C from the more isotropic lateral etching with active defect formation to the anisotropic lateral etching with reduced defect formation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-b2df0260-3a42-4743-b76c-41c2ada31247.jpg)

At T = 300 ◦ C, many small nanopits with rather round shape ( D max ≈ 110 nm) are created. At 400 and 450 ◦ C (not shown), the size distribution becomes wider and the shape becomes more hexagonal. D max is less than 150 nm here. On the other hand, at 500 ◦ C (Figure 2d), the situation changes drastically where fewer large hexagonal nanopits are found ( D max ≈ 300 or 600 nm, depending on the measurement series as described below). This can be seen more clearly in Figure 2h where D max jumps by a factor of two to three across the transition temperature, 450 ◦ C ≤ T c ≤ 500 ◦ C. In contrast, the relative ratio S 1 : S 2 : S 3 ( ≈ 30% : 55% : 15%) does not show such a jump across T c , as shown in Figure 2g. This suggests that the etching efficiency regardless of the anisotropy has a continuous dome-shaped T -dependence sketched by the dashed line in the figure. At T = 700 ◦ C (Figure 2f), only tiny hexagonal nanopits ( D max ≈ 50 nm) are found at a very low density at the surface layer.

ments (Series II) after the first series (I). As seen in Figure 2h, the two data sets agree with each other in terms of the abrupt jump in D max at the same T c . However, D max values at T &gt; T c are two times larger in Series I than II. Also, in Series II, S 3 and S 4 are larger than in Series I where no S 4 is seen. The defect creation rate in Series II seems to be slightly enhanced. One possible explanation for this is a change of oxygen contamination in the gas flow system.

Previous researchers [39, 40] reported similar domeshaped T -dependences of the etching rate but no characteristic change in the nanopit shape (see the dashed line in Figure 2h). In general, etching reactions are expected to be suppressed at sufficiently low T because of the lack of thermal energy. They will also be suppressed at sufficiently high T either due to the instability of CH 4 [50], which is expected to be the final product of the reaction [51], or due to thermal desorption of H before etching the surface [52]. These are probable reasons for such a T dependence.

We checked the reproducibility of the sharp transition at T c by carrying out an independent series of measure-

In order to test the applicability of the direct/remote mechanism to the T -driven mechanism found here, we have examined if l g changes by varying temperature. The glow edge location was again visually observed outside the furnace on the downstream side at P = 13 Pa and W RF = 20 W. Eventually, we found that l g changes but in the opposite direction to that expected from the model at first glance. Specifically, l g increases by 30 mm with increasing T from 300 to 600 ◦ C. Thus some unknown mechanism may control anisotropic etching independently from the direct/remote one. However, this cannot explain why Hug et al. [41] found the anisotropic etching even at T (= 400 ◦ C) ≤ T c . Alternatively, increasing T may affect the composition of the H-species in the sample region. For example, the fraction of Hradicals near the sample position may increase, which promotes anisotropic etching and elongates the glowing length.

In any case, to clarify the mechanism behind the the T -driven transitional change of the etching character we found, direct measurements of the H-plasma components in the sample region and theoretical investigations of

FIG. 3. STM images of graphite surfaces etched by H-plasma for different t of (a) 20, (b) 30, (c) 40, (d) 50, and (e) 60 min. The other parameters are fixed at T = 600 ◦ C, P = 110 Pa and W RF = 20 W. (f) S n and (g) D max are plotted as functions of t . There seems to exist a finite offset time of 10-15 min. before the anisotropic etching starts to proceed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-022f667a-e83f-40da-b708-f761528b5479.jpg)

roles of the temperature and H-radical are highly desired.

## C. Time Dependence

We studied the etching time dependence as well. Figures 3a-e show STM images of graphite surfaces etched for different t , at T = 600 ◦ C, P = 30 Pa, and W RF = 20 W. S 1 decreases linearly with increasing t as expected from a constant etching rate, while S 3 and S 4 remain very small, even when S 1 ≈ 50% at t = 60 min. (see Figure 3f), because the defect formation is rather suppressed at this temperature.

Although D max data are widely scattered, there seems to exist a finite time offset of the order of 10-15 min. before the nanopits start to grow (Figure 3g). Nanopits in the third layer also start to appear after a longer time offset of 30-40 min. The finite time offsets observed clearly in this inactive defect formation condition can give some insights into the kinematics of the surface defect formation. Recent molecular dynamics simulations [52-54] demonstrated that finite fluence is necessary for H ions to be adsorbed, break C-C bonds, and dissociate C atoms away on a graphene or graphite surface.

FIG. 4. STM images of graphite surfaces etched by H-plasma generated by different RF powers W RF of (a) 10, (b) 15, (c) 20, (d) 25, and (e)(f) 50 W. The etching time t is (a)(b)(e)(f) 60 min. and (c)(d) 40 min. The other etching parameters are kept at T = 600 ◦ C and P = 110 Pa. (g) S n and (h) Etching rate are plotted as functions of W RF , there seems to exist a finite offset power of 10-15 W until the anisotropic etching starts. The etching nature seems be changed between 25 and 50 W from the anisotropic lateral etching with reduced defect formation to the isotropic lateral etching with active defect formation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-80fc2830-c797-464f-99b8-55fa8b13c8e8.jpg)

## D. Plasma Power Dependence

Figures 4a-f show STM images of graphite surfaces etched at various W RF from 10 to 50 W. Other parameters are fixed at T = 600 ◦ C and P = 110 Pa except that t = 60 min. for (a)(b)(e)(f) and 40 min. for (c)(d). These parameters belong to the remote plasma mode. As expected, when W RF increases from 10 to 25 W, the change in the hexagonal nanopit formation is more or less the same as the change caused by increasing t in Figures 3a-e. In this range, S 1 decreases and S 2 and S 3 increase monotonically with a finite offset of 10-15 W (see Figures 4g,h). An interesting result was obtained at the highest RF power (= 50 W, Figures 4e,f). Here, the surface morphology completely changes to many multilayer terraced nanopits with irregularly shaped edges just like those in Figures 1a,b.

Again, we have checked a W RF dependence of l g when the glow edge is located outside the furnace, and found that when W RF increases from 25 to 50 W, l g increases by 40-50 mm at T = 600 ◦ C and P = 40 Pa and at T = 23 ◦ C and P = 14 Pa. Although we don't know exactly

FIG. 5. (a) STM image of an HOPG surface etched by H-plasma. The right dark side shows the corner of a hexagonal nanopit with zigzag edges of monolayer height. (b) Tunnel spectra obtained at three characteristically different positions that are near the edge (orange dots), on the upper terrace (black dots), and on the lower terrace (open circles) of the nanopit shown in (a). The specific positions are indicated in (a) by the orange dot ( x = 1 nm, 4 ≤ y ≤ 7 nm), the black dotted line ( x = 23 nm, -19 ≤ y ≤ 12 5 nm), and the white dotted line ( . x = -16 nm, -19 ≤ y ≤ 12 5 nm), respectively. . Only near the edge, a prominent peak at V ≈ 0 and adjacent d I /d V suppressions are found, which correspond to the zigzag edge state. (c)(d) Color plots of d I /d V as a function of y taken in the dotted rectangles in (a), which are at distances x = 10 (c) and 1 nm (d) from the edge, respectively. The edge state exists only in (d) with a spatial variation along the edge (domain structure; see main text). (e) Tunnel spectrum taken at a region of x = 1 nm and - ≤ 3 y ≤ 0 nm where V peak appears at a slightly positive value. All the data were taken at T = 4 7 K. The small dots are the spectrum near a zigzag step edge of monoatomic height at an exfoliated . graphite surface from Ref. 4. (f) Cartoon of hexagonal nanopits synthesized on a graphite surface by H-plasma etching and an STM tip for STS observation. The edge carbon atoms are expected to be terminated by single hydrogen atoms.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-9c8fda4a-fff4-4832-bf8d-bce412d3e571.jpg)

what is happening inside the furnace, it is reasonable to speculate that the observed change in surface morphology is caused when the plasma glow edge moves across the sample region with increasing W RF from 25 to 50 W.

## E. Scanning Tunneling Spectroscopy of Zigzag Edge State at Hexagonal Nanopit Edges

It is possible to judge if the linear edges of hexagonal nanopits are zigzag or armchair type by measuring relative angles ( θ ) between one of atomic rows of the B-site carbon atoms, which are selectively visible with STM, and the edges [4-6, 55]. The data indicate that all the edges of hexagonal nanopits synthesized by the present H-plasma etching technique are exclusively of zigzag type ( θ = 0 , 60 , or 120 ). ◦ ◦ ◦ It is, however, technically very difficult to identify exact arrangements of the endmost carbon atoms along the edges even with the STM technique. Instead, at present, STS observation of the edge state would be the most direct and sensitive detection of the edge quality or purity of zigzag edge structure. This is because the edge state is known to exist in the close vicinity of only the zigzag edge and not the armchair one [2, 4-6]. We have examined many different edges of monolayer depth with a low temperature STM/S apparatus operated at 4.7 K and 78 K in ultra high vacuum.

Figure 5a is an STM image of one of corners of a hexag- onal nanopit with a diagonal length of 160 nm created on an HOPG surface. Etching conditions applied here are the same as those for the sample shown in Figure 3d (600 ◦ C, 110 Pa, 50 min., 20 W). A typical tunneling spectrum taken in the vicinity of the long vertical edge in Figure 5a is shown by the orange circles in Figure 5b, while those taken far from the edge are shown by the black dots (on the upper terrace) and the open circles (on the lower terrace). The locations where the three spectra were taken are indicated by the orange dot (near edge), the black dotted line (upper terrace), and the white dotted line (lower terrace) in Figure 5a. Here, d I /d V , the differential tunnel conductance which is proportional to the local density of states (LDOS) of the graphite surface, is plotted as a function of the bias voltage ( V ) between the STM tip and the sample surface. V corresponds to an energy difference from the E F . In Figure 5b, only in the spectrum near the edge a prominent peak is found at V ≈ 0, which is characteristic of the edge state. The peak voltage ( V peak ) here is -11 mV. In addition to that, the peak is accompanied by clear LDOS suppressions on both sides of the peak. Such clear suppressions have never been reported before. As we reported recently [56], they should be indicative of a large sublattice imbalance [57, 58] indicating the high purity of the zigzag edge synthesized in this work. This is the first spectroscopic evidence for the graphene zigzag edge state created by H-plasma etching.

It is noted that the tunnel spectral results are consistent with Raman spectroscopy measurements of the D band peak which we made separately for multilayer graphene samples with and without the successive anisotropic etching (LP etching) after the initial isotropic etching (HP etching). See Supporting Information for more details.

Next we discuss the spatial variation of the edge state studied by the STS technique particularly along the edge. Let us define xy coordinate axes with the origin O as indicated in Figure 5a. Figures 5c and d are color plots of d /d I V as a function of y (parallel to the edge) taken at the two different distances x = 10 (c) and 1 nm (d) from the edge. The scanned areas for these STS data are indicated by the dotted rectangles in Figure 5a, in which the d /d I V spectra were averaged over 2 nm widths in the x direction. Since the edge state has a short decay length of only 1-2 nm [5, 56], one cannot detect an LDOS enhancement near V = 0 at all at x = 10 nm (see Figure 5c). On the other hand, in the close vicinity of the edge ( x = 1 nm, Figure 5d), the edge state exists almost all along the measurement range except at several positions near y = 2 5 5 , . , and 9 nm where the LDOS peak is less clear or almost disappears. In other words, the edge state is divided spatially into ∼ 3 nm long domains. While V peak is unchanged within each domain, it varies randomly between domains within -20 ≤ V peak ≤ 20 mV. For example, the spectrum averaged over the domain located at -3 ≤ y ≤ 0 nm has a peak at +14 mV (the green large circles in Figure 5e).

We observed similar domain structures for almost all edges we examined, and they may stem from their geometrical imperfections. In any case, the edges must be highly zigzag-rich and presumably of the best quality ever achieved judging from the observed prominent LDOS peak and clear adjacent suppressions. This can be seen by comparing it with the previously known edge state spectra [4, 5]. One of them is shown by the small dots in Figure 5e [4].

## III. CONCLUSION

In summary, by use of STM, we have investigated nanopits with monolayer depth created at graphite surfaces by H-plasma etching varying several key parameters such as temperature, H 2 pressure, etching time, and plasma generation RF power, independently. We found a new sharp T -induced transition of the surface morphology from many small nanopits of rather round shape to few large hexagonal nanopits within a narrow range between 450 and 500 ◦ C. In addition, we observed a P -induced transition from irregularly-shaped multilayer terraced nanopits to few large hexagonal nanopits, which is consistent with a recent report by other researchers. [41] Our own measurement of the pressure dependence of the glowing plasma edge location supports their proposal, i.e. , the lateral etching proceeds isotrop- ically (anisotropically) when the graphite sample is located inside (outside) of the glowing regime. However, the T -induced transition we found here cannot be explained straightforwardly by this model.

Edges of the hexagonal nanopits synthesized here are always aligned to the zigzag direction and not the armchair one. Moreover, characteristic LDOS features of the zigzag edge state, an electronic state localized at the zigzag edge, are clearly observed in STS measurements. This is spectroscopic evidence of the high purity zigzag edge. Thus the anisotropic H-plasma etching is a promising route for the fabrication of zGNR, one of the key device elements for future graphene nanoelectronics. The STS observations also show a spatial distribution of the edge state of the order of 2-3 nm along the edge. This demonstrates the high sensitivity of this tool to evaluate the edge quality on an atomic scale. The present work will hopefully stimulate theoretical works to elucidate the roles of the sample temperature and H-radicals in the anisotropic H-plasma etching of graphite surface and graphene.

## A. Methods

Hydrogen plasma etching was performed with a home made apparatus illustrated in Figure 6. Freshly cleaved graphite samples ( ∼ 5 × 2 × 0 3 mm ) are placed on a . 3 Grafoil plate (Grafoil, Graftech Inc.; 42 × 18 mm 2 with 3 mm high rims at three edges) and inserted in a quartz tube (42 mm i.d. , 1.2 m long), at the center of a muffle furnace. The radial location is 2 mm from the tube inner edge. We examined graphite supplied from three different companies (ZYA grade HOPG, Advanced Ceramic Corp.; PGC, Toshiba Mitsubishi-Electric Industrial Systems Corp.; Kish graphite type-B, Covalent Materials Corp.), and found no essential difference in the etching results. H 2 gas (7N purity) flows at a constant rate, while keeping P at the sample location constant by evacuating with a rotary vacuum pump with pumping speed 200 L/min through a needle valve. P was estimated from pressure values measured at both ends of the quartz tube by pirani gauges and the calculated conductance of the system. H-plasma is generated at 13.56 MHz with a copper coil wound around the quartz tube and an

FIG. 6. Schematic diagram of the experimental apparatus for H-plasma etching. Graphite samples are located in a quartz tube and a muffle furnace under a continuous downstream of the H-plasma.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-4e43fe17-c54f-428b-bde6-b7901b04b11d.jpg)

RF power source (model TX03-9001-00, ADTEC plasma Technology co., Ltd.).

After stabilizing the furnace temperature and the H 2 pressure, the H-plasma etching starts by turning the RF power on. After etching for a certain time duration, the power supplies both to the plasma generator and the furnace are turned off at the same time. It takes about 4 hours for the furnace to cool down to 200 ◦ C from 500 ◦ C while continuing to vacuum pump the quartz tube.

The morphology of etched sample surfaces is examined by two different STMs (UNISOK Co., Ltd.) operated by a commercial controller (model SPM100, RHK Technology, Inc.). All topographic images are obtained in the constant current mode ( I = 1 0 nA, . V = 500 mV) in

- [1] A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81 , 109 (2009).
- [2] M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, J. Phys. Soc. Jpn. 65 , 1920 (1996).
- [3] K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dresselhaus, Phys. Rev. B 54 , 17954 (1996).
- [4] Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. Tsukada, and H. Fukuyama, Appl. Surf. Sci. 241 , 43 (2005).
- [5] Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. Tsukada, and H. Fukuyama, Phys. Rev. B 73 , 085421 (2006).
- [6] Y. Kobayashi, K. Fukui, T. Enoki, K. Kusakabe, and Y. Kaburagi, Phys. Rev. B 71 , 193406 (2005).
- [7] H. Lee, Y.-W. Son, N. Park, S. Han, and J. Yu, Phys. Rev. B 72 , 174431 (2005).
- [8] Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 97 , 216803 (2006).
- [9] L. Pisani, J. A. Chan, B. Montanari, and N. M. Harrison, Phys. Rev. B 75 , 064418 (2007).
- [10] F. Mu˜ noz-Rojas, J. Fern´ndez-Rossier, a and J. J. Palacios, Phys. Rev. Lett. 102 , 136810 (2009).
- [11] Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature 444 , 347 (2006).
- [12] J. Fern´ndez-Rossier a and J. J. Palacios, Phys. Rev. Lett. 99 , 177204 (2007).
- [13] S. Bhowmick and V. B. Shenoy, J. Chem. Phys. 128 , 244717 (2008).
- [14] J. Jiang, W. Lu, and J. Bernholc, Phys. Rev. Lett. 101 , 246803 (2008).
- [15] O. V. Yazyev, R. B. Capaz, and S. G. Louie, Phys. Rev. B 84 , 115406 (2011).
- [16] M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. Rev. Lett. 98 , 206805 (2007).
- [17] C. Stampfer, J. Guettinger, F. Molitor, T. Ihn, and K. Ensslin, Appl. Phys. Lett. 92 , 012102 (2008).
- [18] J. B. Oostinga, B. Sacepe, M. F. Craciun, and A. F. Morpurgo, Phys. Rev. B 81 , 193408 (2010).
- [19] S. Konishi, W. Sugimoto, Y. Murakami, and Y. Takasu, Carbon 44 , 2338 (2006).
- [20] S. S. Datta, D. R. Strachan, S. M. Khamis, and A. T. C. Johnson, Nano Lett. 8 , 1912 (2008).

atmospheric conditions. STS data are taken at 4.7 K in ultra-high vacuum using a home made STM/S system which can work in multi-extreme environments [59].

We thank Y. Ohba and Y. L. Liu for their technical assistances in construction of the H-plasma etching apparatus and STS data analyses. This work was financially supported by Grant-in-Aid for Young Scientists (B) (Grant No. 25800191), Scientific Research (C) (Grant No. 15K05159), and Scientific Research (B) (Grant No. 18H01170) from JSPS. Authors acknowledge the free-touse software, WSxM [60]. We also appreciate the Cryogenic Research Center and the MERIT program of the University of Tokyo for their supply of liquid helium and for allowing us to use the laser Raman microscope, respectively.

- [21] L. Ci, Z. Xu, L. Wang, W. Gao, F. Ding, K. F. Kelly, B. I. Yakobson, and P. M. Ajayan, Nano Res. 1 , 116 (2008).
- [22] L. C. Campos, V. R. Manfrinato, J. D. SanchezYamagishi, J. Kong, and P. Jarillo-Herrero, Nano Lett. 9 , 2600 (2009).
- [23] Y. Sugiyama, O. Kubo, R. Omura, M. Shigehara, H. Tabata, N. Mori, and M. Katayama, Appl. Phys. Lett. 105 , 123116 (2014).
- [24] Y. Y. Li, M. X. Chen, M. Weinert, and L. Li, Nat. Commun. 5 , 4311 (2014).
- [25] P. Nemes-Incze, G. Magda, K. Kamaras, and L. P. Biro, Nano Res. 3 , 110 (2010).
- [26] B. Krauss, P. Nemes-Incze, V. Skakalova, L. P. Biro, K. von Klitzing, and J. H. Smet, Nano Lett. 10 , 4544 (2010).
- [27] F. Oberhuber, S. Blien, S. Heydrich, F. Yaghobian, T. Korn, C. Schueller, C. Strunk, D. Weiss, and J. Eroms, Appl. Phys. Lett. 103 , 143111 (2013).
- [28] G. Z. Magda, X. Jin, I. Hagym´si, P. Vancs´, Z. Osv´th, a o a P. Nemes-Incze, C. Hwang, L. P. Bir´, o and L. Tapaszt´, o Nature 514 , 608 (2014).
- [29] R. K. Puddy, P. H. Scard, D. Tyndall, M. R. Connolly, C. G. Smith, G. A. C. Jones, A. Lombardo, A. C. Ferrari, and M. R. Buitelaar, Appl. Phys. Lett. 98 , 133120 (2011).
- [30] L. Jiao, L. Zhang, Z. Wang, G. Diankov, and H. Dai, Nature 458 , 877 (2009).
- [31] L. Jiao, X. Wang, G. Diankov, H. Wang, and H. Dai, Nat. Nanotechnol. 5 , 321 (2010).
- [32] L. Xie, L. Jiao, and H. Dai, J. Am. Chem. Soc. 132 , 14751 (2010).
- [33] C. Tao, L. Jiao, O. V. Yazyev, Y.-C. Chen, J. Feng, X. Zhang, R. B. Capaz, J. M. Tour, A. Zettl, S. G. Louie, H. Dai, and M. F. Crommie, Nat. Phys. 7 , 616 (2011).
- [34] X. Zhang, O. V. Yazyev, J. Feng, L. Xie, C. Tao, Y.C. Chen, L. Jiao, Z. Pedramrazi, A. Zettl, S. G. Louie, H. Dai, and M. F. Crommie, ACS Nano 7 , 198 (2013).
- [35] X. Li, X. Wang, L. Zhang, S. Lee, and H. Dai, Science 319 , 1229 (2008).
- [36] M. Pan, E. C. Girao, X. Jia, S. Bhaviripudi, Q. Li, J. Kong, V. Meunier, and M. S. Dresselhaus, Nano Lett. 12 , 1928 (2012).
- [37] S. Wang, L. Talirz, C. A. Pignedoli, X. Feng, K. Muellen, R. Fasel, and P. Ruffieux, Nat. Commun. 7 , 11507

(2016).

- [38] P. Ruffieux, S. Wang, B. Yang, C. Sanchez-Sanchez, J. Liu, T. Dienel, L. Talirz, P. Shinde, C. A. Pignedoli, D. Passerone, T. Fumslaff, X. Feng, K. Muellen, and R. Fasel, Nature 531 , 489 (2016).
- [39] R. Yang, L. Zhang, Y. Wang, Z. Shi, D. Shi, H. Gao, E. Wang, and G. Zhang, Adv. Mater. 22 , 4014 (2010).
- [40] G. Diankov, M. Neumann, and D. Goldhaber-Gordon, ACS Nano 7 , 1324 (2013).
- [41] D. Hug, S. Zihlmann, M. K. Rehmann, Y. B. Kalyoncu, T. N. Camenzind, L. Marot, K. Watanabe, T. Taniguchi, and D. M. Zumb¨hl, npj 2D Mater. Appl. u 1 , 21 (2017).
- [42] Z. Shi, R. Yang, L. Zhang, Y. Wang, D. Liu, D. Shi, E. Wang, and G. Zhang, Adv. Mater. 23 , 3061 (2011).
- [43] R. Yang, Z. Shi, L. Zhang, D. Shi, and G. Zhang, Nano Lett. 11 , 4083 (2011).
- [44] G. Xie, Z. Shi, R. Yang, D. Liu, W. Yang, M. Cheng, D. Wang, D. Shi, and G. Zhang, Nano Lett. 12 , 4642 (2012).
- [45] Z. Waqar, J. Mater. Sci. 42 , 1169 (2007).
- [46] B. Eren, D. Hug, L. Marot, R. Pawlak, M. Kisiel, R. Steiner, D. M. Zumb¨ uhl, and E. Meyer, Beilstein J. Nanotechnol. 3 , 852 (2012).
- [47] S. Nunomura and M. Kondo, J. Appl. Phys. 102 , 093306 (2007).
- [48] A. Felten, D. McManus, C. Rice, L. Nittler, J.-J. Pireaux, and C. Casiraghi, Appl. Phys. Lett. 105 , 183104 (2014).
- [49] D. C. Kim, D.-Y. Jeon, H.-J. Chung, Y. Woo, J. K. Shin, and S. Seo, Nanotechnology 20 , 375703 (2009).
- [50] B. J. Wood and H. Wise, J. Phys. Chem. 73 , 1348 (1969).
- [51] Z. J. Pan and R. T. Yang, J. Catal. 123 , 206 (1990).
- [52] A. Davydova, E. Despiau-Pujo, G. Cunge, and D. B. Graves, J. Phys. D: Appl. Phys. 48 , 195202 (2015).
- [53] E. Despiau-Pujo, A. Davydova, G. Cunge, and D. B. Graves, Plasma Chem. Plasma Process. 36 , 213 (2016).
- [54] A. Harpale, M. Panesi, and H. B. Chew, Phys. Rev. B 93 , 035416 (2016).
- [55] Y. Kobayashi, K. Fukui, T. Enoki, and K. Kusakabe, Phys. Rev. B 73 , 125415 (2006).
- [56] A. E. B. Amend, T. Matsui, H. Sato, and H. Fukuyama, e-J. Surf. Sci. Nanotechnol. 16 , 72 (2018).
- [57] A. Cresti, F. Ortmann, T. Louvet, D. V. Tuan, and S. Roche, Phys. Rev. Lett. 110 , 196601 (2013).
- [58] V. M. Pereira, J. M. B. L. dos Santos, and A. H. C. Neto, Phys. Rev. B 77 , 115109 (2008).
- [59] H. Kambara, T. Matsui, Y. Niimi, and H. Fukuyama, Rev. Sci. Instrum. 78 , 073703 (2007).
- [60] I. Horcas, R. Fern´ndez, a J. M. G´ omez-Rodr´ ıguez, J. Colchero, J. G´ omez-Herrero, and A. M. Baro, Rev. Sci. Instrum. 78 , 013705 (2007).
- [61] Y. You, Z. Ni, T. Yu, and Z. Shen, Appl. Phys. Lett. 93 , 163112 (2008).
- [62] C. Casiraghi, A. Hartschuh, H. Qian, S. Piscanec, C. Georgi, A. Fasoli, K. S. Novoselov, D. M. Basko, and A. C. Ferrari, Nano Lett. 9 , 1433 (2009).
- [63] R. Yang, S. Wu, D. Wang, G. Xie, M. Cheng, G. Wang, W. Yang, P. Chen, D. Shi, and G. Zhang, Nano Res. 7 , 1449 (2014).

## SUPPORTING INFORMATION

## Raman Spectroscopy of Multilayer Graphene Etched by H-Plasma

Raman spectroscopy measurements were carried out to cross-check the conclusion drawn from the scanning tunneling spectroscopy (STS) measurements. We need to use sufficiently thin samples (less than 10 layers) otherwise surface information can easily be masked by bulk information since Raman spectroscopy detects molecular vibrations within a few µ m depth from the surface. Such multilayer graphene samples were prepared by exfoliating graphite on a silicon wafer which has a 285 nm thick SiO 2 toplayer.

Figure S1a is an optical microscope image of such a multilayer graphene flake (Sample A) obtained by Hplasma etching under a low H 2 pressure of 13 Pa at T = 600 ◦ C and W RF = 20 W for 15 min, hereafter the LP etching. LP etching creates multi-layer terraced nanopits with high densities. Figure S1b is an atomic force microscope (AFM) image of the region indicated by the black square and arrow in Figure S1a. The maximum thickness of this region is 7 layers or 2.5 nm. The heavily etched surface is very close to the graphite surface etched under similar conditions (see Figure 1a).

By applying successive H-plasma etching at a middle pressure of 190 Pa at T = 600 ◦ C and W RF = 20 W, hereafter the MP etching, following the LP etching, we were able to smoothen the rough surface to a flat one with few large hexagonal nanopits of monolayer depth. Figure S1c is an optical microscope image of a sample obtained by MP etching for 30 min. following LP one for 15 min. (Sample B). The maximum thickness here is 8 layers or 2.7 nm. Figure S1d is an AFM image of the region indicated in (c), where large and well shaped hexagonal nanopits are created as expected. It is interesting to note that the nanopit size is rather uniform here compared to that in graphite (Figure 1f).

Figure S1e shows Raman spectra obtained at the same regions where the AFM images were taken for Sample A and B. Clearly, a significant D band peak appears in Sample A indicating an irregular edge structure (the blue dots and line in the figure). The D band peak originates from intervalley scattering at irregularly shaped edges, including armchair edges [61, 62]. On the other hand, in Sample B, it almost disappears indicating high-purity zigzag edges around the hexagonal nanopits recovered by the successive MP etching (the red dots and line), even though there still remain many edges. This supports the conclusion drawn from the STS measurements. By tak- ing a more detailed look, it can be seen that the D band peak intensity does not return exactly to the background intensity, presumably due to remnant irregularities of the edges.

The present Raman results are essentially the same as those obtained by the layer-by-layer thinning technique developed by Yang et al. [63] where a mild oxygen plasma etching at room temperature was employed instead of the LP H-plasma etching used here.

FIG. S1. (a) Optical microscope image of a multilayer graphene flake (Sample A) etched by H-plasma at T = 600 ◦ C, P = 13 Pa, t = 15 min, and W RF = 20 W (LP etching). (b) AFM image of the region indicated by the square and arrow in (a). (c) Optical microscope image of a multilayer graphene flake (Sample B) etched at T = 600 ◦ C, P = 190 Pa, t = 30 min, and W RF = 20 W (MP etching) after LP etching. (d) AFM image of the region indicated by the square and arrow in (c). (e) Raman spectra for Sample A (blue: region (b)), Sample B (red: region (d)), and an as-exfoliated graphite flake (black). Inset: Magnification of the spectra around the D band peak. All the samples were exfoliated on SiO /Si substrates. 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00182v1-figure-90a2a8b0-97c7-4414-9a25-68e13562754c.jpg)
Tomohiro Matsui Hideki Sato Kazuma Kita André E. B. Amend Hiroshi Fukuyama
Hexagonal Nanopits with the Zigzag Edge State on Graphite Surfaces Synthesized by Hydrogen-Plasma Etching
We studied, by scanning tunneling microscopy, the morphology of nanopits of
monolayer depth created at graphite surfaces by hydrogen plasma etching under
various conditions such as H$_2$ pressure, temperature, etching time, and RF
power of the plasma generation. In addition to the known pressure-induced
transition of the nanopit morphology, we found a sharp temperature-induced
transition from many small rather round nanopits of ~150 nm size to few large
hexagonal ones of 300-600 nm within a narrow temperature range. The remote and
direct plasma modes switching mechanism, which was proposed to explain the
pressure-induced transition, is not directly applicable to this newly found
transition. Scanning tunneling spectroscopy (STS) measurements of edges of the
hexagonal nanopits fabricated at graphite surfaces by this method show clear
signatures of the peculiar electronic state localized at the zigzag edge (edge
state), i.e., a prominent peak near the Fermi energy accompanied by
suppressions on either side in the local density of states. These observations
indicate that the hexagonal nanopits consist of a high density of zigzag edges.
The STS data also revealed a domain structure of the edge state in which the
electronic state varies over a length scale of ~3 nm along the edge. The
present study will pave the way for microscopic understanding of the
anisotropic etching mechanism and of spin polarization in zigzag nanoribbons
which are promising key elements for future graphene nanoelectronics.