# Spatially resolved one-dimensional boundary states in graphene–hexagonal boron nitride planar heterostructures

Jewook Park, Jaekwang Lee, Lei Liu, Kendal Clark, Corentin Durand, Changwon Park, Bobby Sumpter, Arthur Baddorf, Ali Mohsin, Mina Yoon, et al.

► To cite this version:

Jewook Park, Jaekwang Lee, Lei Liu, Kendal Clark, Corentin Durand, et al. Spatially resolved one-dimensional boundary states in graphene–hexagonal boron nitride planar heterostructures. Nature Communications, 2014, 5 (1), pp.5403. 10.1038/ncomms6403. hal-01923494

HAL Id: hal-01923494

https://hal.science/hal-01923494

Submitted on 24 Jun 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

Received 5 Jun 2014 | Accepted 29 Sep 2014 | Published 7 Nov 2014

DOI: 10.1038/ncomms6403

# Spatially resolved one-dimensional boundary states in graphene-hexagonal boron nitride planar heterostructures

Jewook Park<fnref n="1" />,*, Jaekwang Lee<fnref n="1" />,*, Lei Liu², Kendal W. Clark³, Corentin Durand<fnref n="1" />, Changwon Park<fnref n="1" />, Bobby G. Sumpter<fnref n="1" />, Arthur P. Baddorf<fnref n="1" />, Ali Mohsin², Mina Yoon<fnref n="1" />, Gong Gu² & An-Ping Li<fnref n="1" />

Two-dimensional interfaces between crystalline materials have been shown to generate unusual interfacial electronic states in complex oxides. Recently, a one-dimensional interface has been realized in hexagonal boron nitride and graphene planar heterostructures, where a polar-on-nonpolar one-dimensional boundary is expected to possess peculiar electronic states associated with edge states of graphene and the polarity of boron nitride. Here we present a combined scanning tunnelling microscopy and first-principles theory study of the graphene-boron nitride boundary to provide a first glimpse into the spatial and energetic distributions of the one-dimensional boundary states down to atomic resolution. The revealed boundary states are about 0.6 eV below or above the Fermi level depending on the termination of the boron nitride at the boundary, and are extended along but localized at the boundary. These results suggest that unconventional physical effects similar to those observed at two-dimensional interfaces can also exist in lower dimensions.

<footnote n="1">¹Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA. ²The University of Tennessee, Knoxville, Tennessee 37996, USA. ³Central Methodist University, Fayette, Missouri 65248, USA. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to M.Y. (email: myoon@ornl.gov) or to G.G. (email: ggu1@utk.edu) or to A.-P.L. (email: apli@ornl.gov).</footnote>
The ability to control the interfacial properties at the boundary between different materials has formed one of the foundations of modern condensed matter physics and device technology¹⁻⁵. With the advent of graphene and other two-dimensional (2D) crystals, one-dimensional (1D) equivalents of conventional interfaces are envisioned: linear boundaries separating dissimilar materials joint in a single 2D atomic sheet. Graphene and hexagonal boron nitride (hBN) offer a prototypical system to build such a 2D heterostructure⁶⁻¹¹. They are isostructural, nearly lattice-matched and isoelectronic, yet their different band structures and the polar-on-nonpolar boundary promise interesting properties arising from their 2D heterostructures¹²⁻¹⁴, just as for interfaces in the three-dimensional (3D) counterparts¹⁻⁵. Indeed, the graphene-hBN zigzag boundary is predicted to possess peculiar boundary states¹²⁻¹⁴. Half-metallicity has been reported at this boundary based on theoretical calculations¹⁵⁻¹⁹, although zigzag graphene nanoribbons have insulating edge states¹³,²⁰⁻²³. Despite the enormous interest in graphene-hBN heterostructures, the theoretically expected boundary states have not yet been validated experimentally. Some early attempts to address the boundary states have been focused on the connections between graphene and a metallic supporting substrate. For example, when a graphene island is bonded to the edges of a Pt(111) step, a boundary state is shown to be associated with specific sublattices of the graphene structure²⁴. On the other hand, graphene nanoislands grown on Ir(111) do not show boundary states as the graphene edges are hybridized with the substrate, and the edge states at Ir(111) are quenched by the strong interaction between the metal and the graphene²⁵. Moreover, in a 3D polar-on-nonpolar heterostructure of such as Ge-GaAs or LaAlO₃-SrTiO₃, the polar discontinuity introduces an electric field, which could create an unsustainable electrostatic potential at the interface (termed as the 'polar catastrophe'²⁶), leading to electronic or atomic reconstructions³,⁵,²⁶,²⁷ and even preventing an atomically sharp boundary in the latter case⁵,²⁷. It is unclear whether such a polar catastrophe would occur in a 2D polar-on-nonpolar heterostructure of graphene-hBN.

Here we combine scanning tunnelling microscopy (STM) and first-principles theory methods to study the spatial and energetic distributions of the electronic boundary states in the graphene-hBN in-plane heterostructures. The interfaces studied here are crystallographically coherent with sharp transitions from graphene zigzag edges to B (or N)-terminated monolayer hBN on a Cu foil substrate⁶. The boundary states are revealed down to the atomic level with energy levels about 0.6 eV below and above the Fermi level, respectively, for the B- and N-terminated boundaries. The boundary states are shown to be extended along the boundary, and exponentially decay into the bulk of graphene and hBN with decay lengths in a range of 0.24-0.51 nm. The boundary states come from the mixing of π orbitals of C, B and N at each edge, with the occupied bonding states (π<sub>CB</sub>) at B-terminated edges and the unoccupied anti-bonding states (π*<sub>CN</sub>) at N-terminated edges. The effect of the electric field induced by the polarity discontinuity at the interface is discussed, which can be effectively screened out by the charge carriers from graphene and Cu substrate.

## Results

**STM topography of graphene-hBN heterostructures.** An in-plane, single-atomic layer graphene-hBN heterostructure grown on a Cu foil⁶ is used to experimentally investigate the boundary states in this work owing to the weak substrate interaction²⁸ and a sharp boundary at the atomic scale⁶. Figure 1b shows an STM image of two adjacent graphene-hBN junctions formed by an hBN region between two graphene regions, with magnified images of the respective boundaries shown in Fig. 1a and Fig. 1c. The zigzag boundary orientations are indicated in Fig. 1a,c. The image contrast in Fig. 1a is dominated by a substrate step, which appears to be close to an armchair orientation. However, superimposing the schematic atomic lattice on the image and comparing with dI/dV map (details shown in Supplementary Fig. 1) reveal that the structures of the boundary are actually zigzag. The simultaneously acquired scanning tunnelling spectroscopy (STS) map (Fig. 1d) shows distinctively different

<figure><img src="image_3.png" /><figcaption>**Figure 1 |** STM topography of a planar graphene-hBN-graphene heterostructure. (a-c) STM images, with large-scale image in **b** showing two adjacent boundaries (50 × 25 nm², 0.5 V, 200 pA), and **a** and **c** showing atomically resolved images near each boundary (6 × 6 nm²). The edge terminations of hBN are marked with dotted lines. Insets: fast Fourier transform from the images. (d) Simultaneously obtained dI/dV map with image in **b** (AC modulation 500 Hz, 50 mV). (e) Schematic illustration of two boundaries on Cu(100) substrate projected along <011> direction. The height difference of the two junctions is due to an underlying Cu(100) step (Δh = 1.8 Å). 'B' and 'N' indicate the hBN region termination, determined by STS in conjunction with first-principles calculations. Gr and BN represent graphene and hBN regions for simplicity.</figcaption></figure>
contrasts between graphene (blue) and hBN (red) regions. Our assignments of graphene and hBN regions are based on STS measured across the boundaries (Supplementary Figs 2 and 3), where graphene and hBN show distinctively different characteristics with graphene showing a dip at Dirac point of -0.35 V and the hBN showing a large band gap of ~4 eV (refs 6,29–32). The boundary on the right side (Fig. 1c) shows a small height difference (~0.7 Å) originating from the difference in the electronic density of states between the graphene and hBN, and the left boundary (Fig. 1a) has a height difference of 2.5 Å that consists of a substrate step (1.8 Å) on Cu(100), the predominant surface orientation of the Cu foil substrate⁶,²⁹ (Supplementary Fig. 4). Figure 1e depicts a schematic of these two boundaries on the substrate.

The boundary structures are further confirmed by simultaneously acquired atomic-resolution STM images and STS maps. As shown in Supplementary Fig. 2, all the graphene-hBN boundaries are found to contain some segments oriented 120° away from the primary edge direction, which have different terminations than the primary edges, considering the lattice symmetry of hBN and that a zigzag interface is more energetically favourable (about 1.37 eV per edge carbon atom) than a Klein-type interface³³ for this graphene-hBN heterostructure on a Cu substrate. By comparing measured boundary states with theoretical analysis (discussed later), we identify boundary terminations as indicated in Fig. 1a,c.

**Energetic distributions of boundary states.** We now examine the electronic properties of both B- and N-terminated graphene-hBN boundaries by measuring the differential conductance dI/dV, which is proportional to the local electronic density of states (LDOS) distribution, especially at low bias voltages. Figure 2a shows the dI/dV map at a bias voltage of –0.6 V for the two boundaries displayed in Fig. 1b. The electronic states can be seen at both boundaries as intense red-coloured regions. Enhanced electronic states are observed only at the B-terminated boundaries (indicated by green dotted lines), not at the N-terminated boundary (indicated by black dotted lines) at this bias voltage. The measured dI/dV curves in Fig. 2b-d show LDOS variations along three different lines going from graphene to hBN as indicated in Fig. 2a. Remarkably, the dI/dV curves at the B-terminated boundaries (red curves) display enhanced LDOS at around –0.6 V that corresponds to the intense red-coloured region in the dI/dV map (Fig. 2a). On the other hand, the N-terminated boundary (Fig. 2d) does not show an increase in LDOS at –0.6 V; instead, it shows highly localized electronic states at positive bias as explained below.

**Spatial distributions of boundary states.** The spatial confinement of boundary states can be visualized clearly by plotting the measured dI/dV intensity contours at a constant bias voltage. Figure 3b shows the atomically resolved spatial distributions of measured dI/dV across a graphene-hBN boundary, acquired simultaneously with the STM image (Fig. 3a) at 0.6 V. At this bias, the LDOS enhancement appears only at N-terminated edges, and is confined within a few atomic lattices near the boundary. Moreover, the dI/dV intensity at constant bias decays exponentially from the boundary into the graphene and hBN regions (Fig. 3c). We find that the boundary state decays slightly faster in hBN than in graphene, with a larger decay length on the graphene side (0.47 ± 0.08 nm) than the hBN side (0.32 ± 0.05 nm). In contrast, the B-terminated boundary states become more pronounced at –0.6 V as shown in Fig. 3d-f with decay lengths of 0.51 ± 0.17 and 0.24 ± 0.05 nm on the graphene and the hBN side, respectively.

<figure><img src="image_4.png" /><figcaption>**Figure 2 | Electronic boundary states measured by STS.** (a) dI/dV conductance map (45 × 12.5 nm²) of the Gr-hBN-Gr regions measured at -0.6 V to the sample. The boundary states (between green dotted lines) are noticeable as intense red colour. Circles on arrowed lines mark positions of STS from graphene to the hBN, where green lines go across the B-terminated regions and black lines across the N-terminated region. (b-d) dI/dV curves acquired along three lines marked in **a**. The red curves correspond to those acquired between green dotted lines in **a**. Spectroscopy data are offset for clarity.</figcaption></figure>

These boundary states are reproducibly observed in an energy range of 0.45–0.78 eV below the Fermi level for the B-terminated boundaries and 0.60–0.90 eV above the Fermi level for the N-terminated boundaries (see Supplementary Figs 2,5,6 for more details), confirming that the appearance of these states are not experimental artefacts but an intrinsic boundary property. Moreover, all the boundary states are localized at and extended along the boundary, although the electronic states at B- and N-terminated boundaries have clearly different energy levels.

**First-principles calculations of boundary states.** We now compare the experimental observations with first-principles density functional theory (DFT) calculations (see Methods for the details). As illustrated in Fig. 4a, we consider a supercell in which an hBN ribbon of 7 zigzag rows joint by a 12-row graphene ribbon, on a 5-layer Cu slab, oriented (100) with a 0° rotation angle. In view of the experimentally identified superstructure (six graphene or hBN unit cells match to five Cu(100) surface lattice unit cells in the direction perpendicular to the boundary⁶), graphene-hBN supercell is not unique. We consider a flat graphene-hBN on Cu(100) while excluding local lattice distortions that vary with supercell configurations. We project the density of states on each hBN or graphene zigzag row in the unit cell (see Supplementary Fig. 7). To facilitate the comparison of the calculated LDOS with the measured dI/dV, we plot the LDOS contour within a small energy region (from –0.8 eV to 0.4 eV) near the Fermi level in Fig. 4b,c for the B- and N-terminated boundaries, respectively, with the contribution
<figure><img src="image_5.png" /><figcaption>**Figure 3** | Spatial distributions of boundary states. (a,d) Atomically resolved topography image of planar junction showing the seamless connections between graphene (right) and hBN (left) regions. (b,e) Simultaneously obtained dI/dV conductance map with the images (a,d) displaying distinct boundary states in red (with the same colour scheme as shown in Fig. 2a) at different bias voltages with the same set point (200 pA). Image sizes: **a** and **b**, 4.2 × 4.2 nm²; **d** and **e**, 4 × 4 nm². Hexagonal lattice of hBN is overlaid on topography. Green and black dotted lines indicate B- and N-terminated boundaries, respectively. (c,f) Averaged line profiles across the boundary regions (marked with boxes in **b** and **e**) showing exponential decay of the interfacial state into the both graphene and hBN regions. Red and blue lines are exponential fitting curves.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4** | Calculated boundary states of the planar graphene-hBN-graphene heterostructure on a Cu(100) substrate. (a) Schematic of the heterostructure (top view). Seven hBN zigzag rows are laterally sandwiched by 12 graphene rows. Five Cu layers are considered as the substrate. B atoms are specified in blue, N atoms in green and carbon atoms in black. (b) Contour map related to the B-terminated state made up of LDOS projected onto 11 zigzag rows (5 graphene and 6 hBN rows) across the B-terminated boundary. B-terminated states are located around 0.6 eV below the Fermi level. (c) Contour map related to the N-terminated state across the N-terminated boundary. The boundary states are located around 0.1 eV above the Fermi level.</figcaption></figure>

from Cu excluded for clarity. For the B-terminated boundary (Fig. 4a), the boundary states are located about 0.6 eV below the Fermi level and localized within only a few rows (row-to-row distance is 2.165 Å, given a 1.44-Å C–C or B–N length) of

the interface. The spatial distribution of LDOS extends into the graphene region more deeply than into the hBN region, with the decay length in graphene 3.1 Å and in hBN 2.0 Å, respectively. The highest intensity of LDOS appears at the C site for both
B- and N-terminated boundary states, consistent with the measured *dI/dV* (Fig. 3c,f). Both the calculated energy position of the B-terminated boundary state and its decay length are in good agreement with the STS measurement.

The calculations employing the DFT exchange and correlation (XC) functional of the local density approximation (LDA) also yield highly localized states at the N-terminated boundary located around 0.1 eV above the Fermi level (Fig. 4c), lower than that measured by STS mapping. As LDA calculations usually underestimate the band gaps and especially the energies of unoccupied states, we further performed Heyd–Scuseria–Ernzerhof (HSE) hybrid functional calculations³⁴ on a small supercell (see Methods for details) to account for the discrepancy in the N-terminated boundary state level. It is found that the energy level of B-terminated states remain unchanged from the LDA result on the same supercell, but the N-terminated states move up to ~0.6 eV above the Fermi level (see Supplementary Fig. 8); and the shift accounts for the difference in the N-terminated boundary states between the experiment and the LDA results in Fig. 4c.

Although the interactions with the Cu substrate are very weak²⁸, we find that there are still hybridizations between the *d*-orbitals of Cu and the π-orbitals of graphene and hBN near the B-terminated boundary (see Supplementary Fig. 9). However, no appreciable interaction is seen between N-terminated states and Cu states. These dissimilar interactions arise from the different electronegativities of B and N. As a comparison, we have compared the calculated LDOS with and without a Cu substrate (Supplementary Fig. 10) and found that the N-terminated boundary states remain largely unchanged, but the B-terminated states are modified by the Cu substrate due to the hybridization with Cu substrate.

In comparison, we carried out DFT calculations on an in-plane graphene–hBN junction on Cu(111). As shown in Supplementary Fig. 11, the calculated results for both B- and N-terminated boundaries are essentially the same as those on Cu(100), confirming that crystallographic orientation of Cu has little effect on the boundary states.

## Discussion

We now clarify the nature of the observed boundary states on the basis of the combined experimental and theoretical study. When two semi-infinite sheets of graphene and hBN form a zigzag heterojunction on a Cu substrate, the mixing of π-orbitals of C, B and N at each edge gives rise to four sets of bands that correspond to the bonding and anti-bonding states between C–N and C–B (Supplementary Fig. 12). Although the boundary states are spin polarized in the free-standing case¹³ as are the edge states of free-standing zigzag hBN nanoribbons³⁵, calculations by us and others³⁶ indicate that the boundary magnetism disappears upon adsorption onto a Cu substrate. The relevant bands, close to the Fermi level, are the occupied bonding states ($\pi_{CB}$) at B-terminated edges and the unoccupied anti-bonding states ($\pi^*_{CN}$) at N-terminated edges, both localized at the interface. The localized states lay ~0.6 eV below and above the Fermi level, respectively (Fig. 3).

In addition, different electronegativities of B and N make the hBN strip the 2D equivalent to a polar slab. We calculated that a free-standing eight-row hBN strip experiences a polar field of ~0.11 eV Å⁻¹, which is reduced to about 0.06 eV Å⁻¹ by screening when joint with graphene. The Cu substrate brings about an additional screening, as corroborated in our calculations. In a 3D polar-on-nonpolar heterostructure, 2D charge sheets create a constant overall electric field and the electrostatic potential would diverge with increasing thickness, leading to a

'polar catastrophe'²⁶. The polar catastrophe can be averted by charge reconstruction²⁶ or atom rearrangement⁵. For (100) or (111) Ge–GaAs interfaces, which are the higher-dimensional analogies of the zigzag graphene–hBN boundaries, atom rearrangement prevents an atomically sharp interface⁵,²⁷. In a 2D polar-on-nonpolar heterostructure, the electric field induced by the polar discontinuity at the boundary is inversely proportional to the distance, and the electrostatic potential would also diverge although not as fast as in 3D. Furthermore, the charge carriers from graphene and Cu substrate can effectively screen out the polar field. As such, the polarity discontinuity in 2D may not be as catastrophic. Nevertheless, planar graphene–hBN heterostructures provide an excellent platform for exploring novel physics in a 1D system—the interface in 2D space.

## Methods

**Sample preparation and STM/STS experiments.** The graphene–hBN planar junction on Cu foil sample was prepared using the 2D heteroepitaxial growth technique⁶. Moreover, (100) has been found to be the predominant Cu surface orientation as confirmed by low-energy electron microscopy and low-energy electron diffraction⁶ and X-ray diffraction and electron backscatter diffraction (Supplementary Fig. 4). STM/STS were performed using a variable-temperature STM (Omicron) with Nanonis (SPECS) controller and an electrochemically etched tungsten tip. All measurements were carried out in ultra-high vacuum (<3 × 10⁻¹⁰ Torr) at room temperature. We annealed the sample at 430 °C for 48 h for cleaning. We used constant current mode for topography images and a lock-in technique for differential conductance measurement with 500 Hz, 40–100 mV AC modulations. All STM images and spectroscopy data were analysed and processed with WSxM software.

**First-principles calculations.** First-principles calculations, based on DFT, were performed using the Vienna Ab initio Simulation Package. The projector-augmented wave method was used to mimic the ionic cores, while the LDA considering the Ceperly–Alder–Perdew and Zunger (CA-PZ) functional was employed for the XC functional. In addition, an HSE-screened hybrid XC functional³⁴ was used to more properly determine the energies of interface states. The overall fraction of Fock exchange and the length scale $\omega^{-1}$ for exchange screening are adopted as the HSE12s (ref. 37) form of $a = 0.425$ and $\omega = 0.408$ Å⁻¹. Atomic positions, as well as lattice parameters, were optimized using a conjugate gradient algorithm. The ionic and electronic relaxations were performed by applying a convergence criterion of $5 \times 10^{-2}$ eV Å⁻¹ per ion and $10^{-4}$ eV per electronic step, respectively. The rectangular graphene–BN hybrid structure, 43.3 × 2.5 Å, is used for the LDOS calculations, and a vacuum of 16 Å between the hybrids is considered. Also, $1 \times 10 \times 1$ Monkhorst–Pack meshes were used to perform the integration over the Brillouin zone.

## References

1. Bert, J. A. et al. Direct imaging of the coexistence of ferromagnetism and superconductivity at the LaAlO₃/SrTiO₃ interface. *Nat. Phys.* **7**, 767–771 (2011).

2. Li, L., Richter, C., Mannhart, J. & Ashoori, R. C. Coexistence of magnetic order and two-dimensional superconductivity at LaAlO₃/SrTiO₃ interfaces. *Nat. Phys.* **7**, 762–766 (2011).

3. Ohtomo, A. & Hwang, H. Y. A high-mobility electron gas at the LaAlO₃/SrTiO₃ heterointerface. *Nature* **427**, 423–426 (2004).

4. Reyren, N. et al. Superconducting interfaces between insulating oxides. *Science* **317**, 1196–1199 (2007).

5. Kroemer, H. Nobel Lecture: Quasielectric fields and band offsets: teaching electrons new tricks. *Rev. Mod. Phys.* **73**, 783–793 (2001).

6. Liu, L. et al. Heteroepitaxial growth of two-dimensional hexagonal boron nitride templated by graphene edges. *Science* **343**, 163–167 (2014).

7. Han, G. H. et al. Continuous growth of hexagonal graphene and boron nitride in-plane heterostructures by atmospheric pressure chemical vapor deposition. *ACS Nano* **7**, 10129–10138 (2013).

8. Sutter, P., Cortes, R., Lahiri, J. & Sutter, E. Interface formation in monolayer graphene-boron nitride heterostructures. *Nano Lett.* **12**, 4869–4874 (2012).

9. Ci, L. et al. Atomic layers of hybridized boron nitride and graphene domains. *Nat. Mater.* **9**, 430–435 (2010).

10. Gao, Y. et al. Toward single-layer uniform hexagonal boron nitride–graphene patchworks with zigzag linking edges. *Nano Lett.* **13**, 3439–3443 (2013).

11. Kim, S. M. et al. Synthesis of patched or stacked graphene and hBN flakes: a route to hybrid structure discovery. *Nano Lett.* **13**, 933–941 (2013).

12. Okada, S., Igami, M., Nakada, K. & Oshiyama, A. Border states in heteroheets with hexagonal symmetry. *Phys. Rev. B* **62**, 9896–9899 (2000).
13. Pruneda, J. M. Origin of half-semimetallicity induced at interfaces of C-BN heterostructures. *Phys. Rev. B* **81**, 161409 (2010).

14. Liu, Y., Bhowmick, S. & Yakobson, B. I. BN white graphene with 'colorful' edges: the energies and morphology. *Nano Lett.* **11**, 3113-3116 (2011).

15. Bhowmick, S., Singh, A. K. & Yakobson, B. I. Quantum dots and nanoroads of graphene embedded in hexagonal boron nitride. *J. Phys. Chem. C* **115**, 9889-9893 (2011).

16. Liu, Y., Wu, X., Zhao, Y., Zeng, X. C. & Yang, J. Half-metallicity in hybrid graphene/boron nitride nanoribbons with dihydrogenated edges. *J. Phys. Chem. C* **115**, 9442-9450 (2011).

17. Xiao, H. P. et al. Size effect of half-metallic properties of BN/C hybrid nanoribbons. *Phys. B: Cond. Matt.* **407**, 4770-4772 (2012).

18. Bristowe, N. C., Stengel, M., Littlewood, P. B., Artacho, E. & Pruneda, J. M. One-dimensional half-metallic interfaces of two-dimensional honeycomb insulators. *Phys. Rev. B* **88**, 161411 (2013).

19. Dutta, S., Manna, A. K. & Pati, S. K. Intrinsic half-metallicity in modified graphene nanoribbons. *Phys. Rev. Lett.* **102**, 096601 (2009).

20. Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S. & Geim, A. K. The electronic properties of graphene. *Rev. Mod. Phys.* **81**, 109-162 (2009).

21. Son, Y.-W., Cohen, M. L. & Louie, S. G. Half-metallic graphene nanoribbons. *Nature* **444**, 347-349 (2006).

22. Ritter, K. A. & Lyding, J. W. The influence of edge structure on the electronic properties of graphene quantum dots and nanoribbons. *Nat. Mater.* **8**, 235-242 (2009).

23. Tao, C. et al. Spatially resolving edge states of chiral graphene nanoribbons. *Nat. Phys.* **7**, 616-620 (2011).

24. Merino, P. et al. Sublattice localized electronic states in atomically resolved graphene-Pt(111) edge-boundaries. *ACS Nano* **8**, 3590-3596 (2014).

25. Li, Y. et al. Absence of edge states in covalently bonded zigzag edges of graphene on Ir(111). *Adv. Mater.* **25**, 1967-1972 (2013).

26. Nakagawa, N., Hwang, H. Y. & Muller, D. A. Why some interfaces cannot be sharp. *Nat. Mater.* **5**, 204-209 (2006).

27. Harrison, W. A., Kraut, E. A., Waldrop, J. R. & Grant, R. W. Polar heterojunction interfaces. *Phys. Rev. B* **18**, 4402-4410 (1978).

28. Khomyakov, P. A. et al. First-principles study of the interaction and charge transfer between graphene and metals. *Phys. Rev. B* **79**, 195425 (2009).

29. Tian, J. et al. Graphene induced surface reconstruction of Cu. *Nano Lett.* **12**, 3893-3899 (2012).

30. Gao, L., Guest, J. R. & Guisinger, N. P. Epitaxial graphene on Cu(111). *Nano Lett.* **10**, 3512-3516 (2010).

31. Hwang, B. et al. Electron-beam assisted growth of hexagonal boron-nitride layer. *Curr. Appl. Phys.* **13**, 1365-1369 (2013).

32. Joshi, S. et al. Boron nitride on Cu(111): an electronically corrugated monolayer. *Nano Lett.* **12**, 5821-5828 (2012).

33. Klein, D. J. Graphitic polymer strips with edge states. *Chem. Phys. Lett.* **217**, 261-265 (1994).

34. Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. *J. Chem. Phys.* **118**, 8207-8215 (2003).

35. Zheng, F. et al. Half metallicity along the edge of zigzag boron nitride nanoribbons. *Phys. Rev. B* **78**, 205415 (2008).

36. Li, Y. & Mazzarello, R. Structural and electronic properties of hybrid graphene and boron nitride nanostructures on Cu. *Phys. Rev. B* **88**, 045317 (2013).

37. Moussa, J. E., Schultz, P. A. & Chelikowsky, J. R. Analysis of the Heyd-Scuseria-Ernzerhof density functional parameter space. *J. Chem. Phys.* **136**, 204117 (2012).

## Acknowledgements

This research was conducted at the Center for Nanophase Materials Sciences, which is sponsored at Oak Ridge National Laboratory by the Scientific User Facilities Division, Office of Basic Energy Sciences, US Department of Energy (DOE). A portion of theory work (J.L.) is supported by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US DOE. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the US DOE under contract no. DE-AC02-05CH11231. The work at UTK was supported by NSF (ECCS-1231808) and DARPA (approved for public release; distribution is unlimited).

## Author contributions

J.P. and J.L. contributed equally to this work. A.-P.L. conceived the project and designed the STM experiments and M.Y. designed the theory tasks. Measurements were carried out by J.P., K.W.C, C.D. and A.-P.L.; samples were prepared by L.L. and G.G.; and the first-principles calculations were performed by J.L., C.P. and M.Y. All authors participated in the data interpretation and paper writing.

## Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Park, J. et al. Spatially resolved one-dimensional boundary states in graphene-hexagonal boron nitride planar heterostructures. *Nat. Commun.* 5:5403 doi: 10.1038/ncomms6403 (2014).