Oriented graphene nanoribbons embedded in hexagonal boron nitride trenches

Abstract

Graphene nanoribbons (GNRs) are ultra-narrow strips of graphene that have the potential to be used in high-performance graphene-based semiconductor electronics. However, controlled growth of GNRs on dielectric substrates remains a challenge. Here, we report the successful growth of GNRs directly on hexagonal boron nitride substrates with smooth edges and controllable widths using chemical vapour deposition. The approach is based on a type of template growth that allows for the in-plane epitaxy of mono-layered GNRs in nano-trenches on hexagonal boron nitride with edges following a zigzag direction. The embedded GNR channels show excellent electronic properties, even at room temperature. Such in-plane hetero-integration of GNRs, which is compatible with integrated circuit processing, creates a gapped channel with a width of a few benzene rings, enabling the development of digital integrated circuitry based on GNRs.

ARTICLE

Received 9 Dec 2016 | Accepted 24 Jan 2017 | Published 9 Mar 2017

Lingxiu Chen1,2,*, Li He1,3,*, Hui Shan Wang1,4,*, Haomin Wang1, Shujie Tang1,5, Chunxiao Cong6,7, Hong Xie1,

Lei Li1,4, Hui Xia8, Tianxin Li8, Tianru Wu1, Daoli Zhang3, Lianwen Deng4, Ting Yu6, Xiaoming Xie1,2 & Mianheng Jiang1,2

Graphene nanoribbons (GNRs) are ultra-narrow strips of graphene that have the potential to be used in high-performance graphene-based semiconductor electronics. However, controlled growth of GNRs on dielectric substrates remains a challenge. Here, we report the successful growth of GNRs directly on hexagonal boron nitride substrates with smooth edges and controllable widths using chemical vapour deposition. The approach is based on a type of template growth that allows for the in-plane epitaxy of mono-layered GNRs in nano-trenches on hexagonal boron nitride with edges following a zigzag direction. The embedded GNR channels show excellent electronic properties, even at room temperature. Such in-plane hetero-integration of GNRs, which is compatible with integrated circuit processing, creates a gapped channel with a width of a few benzene rings, enabling the development of digital integrated circuitry based on GNRs.

1 State Key Laboratory of Functional Materials for Informatics, Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences, 865 Changning Road, Shanghai 200050, China. 2 School of Physical Science and Technology, ShanghaiTech University, 319 Yueyang Road, Shanghai 200031, China. 3 School of Optical and Electronic Information, Huazhong University of Science and Technology, Wuhan 430074, China. 4 School of Physics and Electronics, Central South University, Changsha 410083, China. 5 Graduate University of the Chinese Academy of Sciences, Beijing 100049, China. 6 Division of Physics and Applied Physics, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371, Singapore. 7 State Key Laboratory of ASIC & System, School of Information Science and Technology, Fudan University, Shanghai 200433, China. 8 National Laboratory for Infrared Physics, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, 500 Yu Tian Road, Shanghai 200083, China.* These authors contributed equally to this work. Correspondence and requests for materials should be addressed to H.W. (email: mailto:[email protected]

Web End [email protected] ) or to X.X. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:14703 | DOI: 10.1038/ncomms14703 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms14703 OPEN

Oriented graphene nanoribbons embedded in hexagonal boron nitride trenches

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14703

Ideal graphene nanoribbons (GNRs) have been shown to exhibit extreme chirality dependence as metals or semiconductors1. Therefore, the capability to precisely produce GNRs

with dened chirality at the atomic level is required in order to engineer their band gap and electrical properties2,3. Conventional lithography46 always results in ragged edges along the GNRs. Other GNR synthesis methods, including sidewall growth on SiC7, advanced lithography8, sono-chemical methods9,10 and carbon nanotube unzipping11,12, still present difculties for either chirality control or width-scaling down to 10 nm and less. Recently, bottom-up synthesis methods using catalytic substrates1315 were demonstrated to form GNRs with well-dened edge structures and atomic precision. However, transferring techniques to produce devices without degrading the quality of the GNRs remain a formidable challenge. It is obvious that earlier approaches have fundamental limitations for further electronic investigation. Electronics always require scalable transfer-free approaches for growing GNRs and conducting band gap engineering. Controlled fabrication of oriented GNRs embedded on hexagonal boron nitride (h-BN) has the capability to overcome the above difculties. With proper control, the band gap and magnetic properties can be precisely engineered. Most desired features for GNRs can be automatically attained using this approach.

Here, we demonstrate the successful growth of GNRs directly on h-BN substrates with smooth edges and controllable widths via templated growth using chemical vapour deposition (CVD). Transistors made with sub-10-nm GNRs demonstrate large on off ratios of more than 104 at room temperature and carrier mobility values of B750 cm2 V 1 s 1. For the narrowest GNRs, the band gaps extracted from the electrical transport data are 40.4 eV.

ResultsTemplated growth of oriented GNRs in h-BN trenches. Figure 1 demonstrates our conceptual design and its experimental verication for the synthesis of GNRs on h-BN via templated growth (details of the experimental procedures are given in the Methods section). Single-crystal h-BN akes exhibited a very smooth surface after multiple cleaning steps for removing possible surface contaminants (see the atomic force microscopy (AFM) images in Fig. 1d and Supplementary Fig. 1a,b). The nano-trenches on h-BN were synthesized by nickel particle-assisted etching (Fig. 1b,e and Supplementary Fig. 1c,d). After CVD growth, the

trenches were lled with graphene (Fig. 1c,f and Supplementary Fig. 1e,f).

AFM measurement of oriented trenches in h-BN crystals. The nano-trenches exhibited predominant separation angles of 60 relative to each other with no obvious correlation to the direction of the gas ow, indicating anisotropic etching along the different crystallographic directions (Supplementary Figs 2 and 3). Separation angles of 30 and 90 were rarely observed. The lattice structures shown in the inset of Supplementary Figs 2a,b and 3a conrm that the crystallographic orientations of the h-BN trenches/edges follow a zigzag (ZZ) pattern. The crystallographically selective chemical reaction indicates lower activation energy along the ZZ patterns of the h-BN.

Figure 2ac shows some examples of nano-trenches obtained on h-BN substrates. The width of the trenches and their distribution were strongly dependent on the etching parameters and solutions (see Supplementary Table 1). Nanometre-sized-trenches, narrower than 10 nm in width, could be reproduced via an optimized process. Wider nano-trenches (Fig. 2c) could be synthesized by increasing the etching duration or temperature (see Supplementary Table 1). The trench in Fig. 2c had a depth of 0.334 nm with a bottom roughness comparable to that of pristine h-BN, indicating that the trench was mono-layered. The extremely anisotropic two-dimensional (2D) nature of h-BN predominantly conned the etching to a single-atom layer. Due to the size limitation of the AFM tips, the depth of the narrower trenches could not be exactly determined, but these trenches were reasonably judged to also be mono-layered based on the measurements of wider nano-trenches. The tip size limitation may also result in uncertainty in the width measurements of the nano-trenches.

AFM measurement of oriented GNRs on h-BN crystals. Figure 2df presents typical AFM friction images of GNRs embedded in the nano-trenches of the h-BN substrate, and the corresponding height images are given in Supplementary Fig. 4. A super-lattice structure, also known as a Moir pattern, is clearly visible in Fig. 2e,f for the synthesized GNRs with widths of B58 and B56 nm, respectively, which was better contrasted in the friction images than in the height images. For narrow GNRs, a Moir pattern was not observed (Fig. 2d and Supplementary Fig. 4a) because the GNR width was comparable to or smaller than the periodicity of the Moir pattern. Height variations in the range of 2040 pm could be seen near the boundary of the

a

b

Etching GNR growth

c

Figure 1 | Formation of GNRs in h-BN trenches. (a) Smooth surface of the h-BN; (b) Synthesis of nano-trenches on h-BN by Ni particle-assisted etching; (c) In-plane epitaxial template growth of GNRs via CVD; (df) AFM friction images corresponding to the schematics shown in ac, respectively. The friction images showed better contrast than the height images, especially for GNRs embedded in the h-BN nano-trenches. Scale bars, 200 nm.

2 NATURE COMMUNICATIONS | 8:14703 | DOI: 10.1038/ncomms14703 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14703 ARTICLE

a

6 nm

b

c

53 nm

0.334 nm

8 nm

0.136 nm

0.2 nm

0.2 nm 0.3 nm

0.3 nm

0.2 nm

0.2 nm

d

e

f

2.3 mV 4.1 mV

4.3 mV

58 nm

56 nm

6 nm

2.4 mV

4.5 mV

4.6 mV

Figure 2 | Images of nano-trenches and GNRs on the h-BN substrate. AFM height images of mono-layered nano-trenches on the h-BN crystal with the following widths: (a) B6 nm, scale bar: 10 nm, (b) B8 nm, scale bar: 20 nm and (c) B53 nm, scale bar: 40 nm. AFM friction images of graphene nanoribbons (GNRs) embedded in the nano-trenches on the h-BN crystal via template growth with the following widths: (d) B6 nm, scale bar: 10 nm, (e) B58 nm, scale bar: 40 nm and (f) B56 nm, scale bar: 40 nm. The giant Moir pattern can be observed in e,f, indicating that the GNRs are precisely aligned with the h-BN.

graphene/h-BN, which is most likely due to lattice distortions caused by mismatched lattice constants and coefcients of thermal expansion between the graphene and h-BN. The small outof-plane distortions also excluded the possibility of the formation of multilayered GNRs. More height images of the GNRs are shown in Supplementary Fig. 5.

The existence of the giant Moir pattern indicates that the graphene was highly crystalline and precisely aligned with the h-BN underneath. It was noticed that the Moir pattern appeared to be stretched along the GNR, while it appeared relaxed laterally. This trend differs from regular hexagons with a periodicity of B14 nm, which have always been observed with well-aligned graphene domains on h-BN1618. This observation gives a strong indication of the in-plane epitaxy between the graphene and the h-BN at the edges of the trench, where the graphene is stretched by tensile strain along the ribbon, due to a lattice mismatch between the graphene and h-BN. Supplementary Fig. 6 illustrates how the stretched Moir pattern formed, and Supplementary Fig. 7 shows the dependence of its wavelength along the stretched direction on the strain level. According to the relationship shown in Supplementary Fig. 7a, the strain in the GNR along the ribbon was estimated to be B0.750.3% from the wavelength of the stretched Moir pattern measured in the AFM images. In addition, atomic-resolution AFM images conrmed that the in-plane connections between the GNRs and h-BN were continuous (see Supplementary Fig. 8).

Supplementary Fig. 9 shows AFM images taken before, during and after the GNR growth, illustrating that the graphene grew via a step-ow mechanism from two step-edges of the h-BN top-layer trench and coalesced into a complete GNR. The growth of the GNRs was also observed to occur at a one-sided atomic step-edge, which was formed by mechanical cleavage. Supplementary Fig. 10 shows a case in which the nanoribbon grew along the atomic step-edge of the top layer on the h-BN and developed laterally. Atomic-resolution images show that there were no discernible rotational misalignments between the three lattices (the graphene epi-ribbon, the top h-BN layer and the underlying h-BN lattice), indicating that while the GNR developed, the chirality of its edge was unchanged. Of course, in this case, the width of the GNR was determined by the growth time, not by the nano-trench connement.

Raman characterization of GNRs embedded in the h-BN. Raman measurements were carried out to investigate the structural and electronic properties of the nanoribbons. Figure 3 shows the AFM image of a 15-nm-wide GNR and its Raman spectrum. In the spectra of Fig. 3b, a prominent sharp peak appeared at B1,365 cm 1, which was attributed to the Raman-active LO phonon of h-BN (ref. 19). For the spectrum of the GNR, the G-,

D-, D0- and 2D-bands were tted with a single-Lorentzian line shape (see Supplementary Fig. 11). The Raman spectra show a prominent G-band (B1,572.1 cm 1) as well as a single-Lorentzian-shaped 2D-band (B2,668.4 cm 1), which were expected for monolayer graphene. Compared to that of the pristine graphene domain18, the red shift of B9 cm 1 in the G-band position indicates the existence of a tensile strain of B0.6% in the GNR20,21. Furthermore, the GNR spectrum shows a D-band at 1,334.6 cm 1. At 1,617.5 cm 1 (D0-band), a tiny shoulder appears on the right side of the G band. Both the Raman

D- and D0-bands may have primarily originated from the lattice distortions and disorder at the GNR/h-BN boundary. The 2D-band exhibits a characteristic of I(2D)/I(G)41, also indicating the single-layer nature of the GNR, where I(2D) and I(G) represent the intensity values of the 2D and G-bands, respectively. In addition, the observed broadening of the full width at half-maximum for the Raman G- and 2D-bands may have been due to strain variations at the nanometre-scale22.

Electronic transport properties of GNR transistors. To investigate the electrical properties of the GNRs, eld effect transistor (FET) devices were produced using p-type doped Si/SiO2 (300 nm thick for the SiO2 layer; for details of the device fabrication, see Methods section and Supplementary Figs 1116). Figure 4 shows the eld effect characteristics of representative GNRs at different widths. For the 15-nm GNR, the plot of conductance G versus Vgate at different temperatures, shown in Fig. 4a, exhibits obvious modulations due to the external electrical eld. The eld effect mobility extracted from the 15-nm GNR was B916.1 cm2 V 1 s 1 at 300 K. According to a simple two-band model (see Supplementary Note 3), the value of the band gap could be extracted by tting the resistancetemperature curve (Fig. 4b; for details of the tting, please see Supplementary

NATURE COMMUNICATIONS | 8:14703 | DOI: 10.1038/ncomms14703 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14703

1,334.6 (55.1)

a b

0.2 nm

0 nm

Intensity (a.u.)

GNR

1,200 1,600 2,400 2,800

Wavenumber (cm1)

GNR

Figure 3 | Raman spectrum of a GNR on h-BN. (a) AFM image of a GNR with a width of B15 nm. Scale bar, 30 nm; (b) Raman scattering of the GNR. The Raman spectrum of h-BN is also shown beside that of the GNR for comparison. The spectrum traces were normalized and shifted on the intensity axis for clarity. The inset shows the Raman spectrum of the GNR after subtracting the h-BN background. Because the position of the D-band of the GNR is very close to that of a prominent Raman peak of h-BN, such a subtraction could be utilized to identify the existence of weak Raman peaks. The full width at half-maximum (FWHM) for each peak is given in parentheses with the peak position. The band names of the main peaks are also indicated. The parameters of all Raman peaks were extracted via Lorentzian tting, and the wavelength of the exciting laser was 488 nm.

0.01 0.02 0.03 0.30 0.45

a

b

100

Resistance (M)

0.1

0.01

10

G(e2 h1 )

G(e2 h1 )

1E-3

1

1E-4

2 K

0.1

1E-5

30 20 10 0 10 20 30

Vgate (V)

Vgate (V)

1/T (K1)

0 10 20 30 40 50 60

c

0.8

d

0.6

0.6

0.4

0.4

E g(eV)

0.2

0.2

0.0

0.0

40 20 0 20 40

w (nm)

Figure 4 | Electronic transport through GNR devices on h-BN. (a) Conductance (G) of GNRs with a width of B15 nm as a function of the back gate voltage (Vgate) at different temperatures; (b) Arrhenius plot of the resistance of the 15-nm GNR FET under different Vgate values at temperatures from 2 to 250 K. The solid curves are ts based on a simple two-band (STB) model (see Supplementary Note 3); (c) Electronic transport through a narrow ribbon with a width of B5 nm. Its conductance can be completely switched off, even at 300 K. The inset shows the conductance at Vgate 30 V versus inverse temperature from 200 to 300 K. The dashed line is t to the experimental data according to Gon exp

E 2k T

; (d) Band gap Eg extracted from experimental data for GNRs versus their ribbon width (w). The orange curve is a t of our experimental data with Eg (eV)Ba/(w b) (a is in units of eVnm, both b and

w are in units of nm). The light purple (a 0.8 eVnm) and light green curves (a 1.0 eVnm) are empirical curves adapted from ref. 9 and refs 24,25,

respectively. Error bars for the experimental data represent standard deviation of uncertainty in AFM measurement or gap extraction.

Information). The tting includes contributions from both the normal thermal activation and electrical contact resistance. The extracted band gap for the 15-nm GNR was 12023 meV. Next, the narrowest GNRs were measured. Some devices exhibited obvious transistor behaviour, even at room temperature. As shown in Fig. 4c, a B5-nm GNR device had a Gon/Goff ratio of

4104. Its conductance appeared to be completely insulating, with

Go10 4 e2 h 1 10 10 S from Vgate 0 to 10 V (off-

state); then, it gradually switched on, exhibiting a relatively high G 0.7 e2 h 1 at room temperature for Vgate 40 V. The

scattering mean free path was estimated to be B50 nm. The mobility of the GNR was estimated to be B765 cm2 V 1 s 1 at 300 K. Accurately measuring the off-state resistance of the narrowest GNR FET became challenging due to noise and set-up

4 NATURE COMMUNICATIONS | 8:14703 | DOI: 10.1038/ncomms14703 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14703 ARTICLE

limitations. It was found that the conductance near the on-state exhibited an exponential relationship with inverse temperature from 200 to 300 K. This result indicates that a Schottky barrier dominated the conductance of the GNR FET because of the high work functions of Pd and Ni (refs 9,23). We estimated the band gaps (Eg) of the narrow GNRs by tting the conductance near the on-state Gon / IonIoff exp

Eg

2kBT

, where kB is Boltzmanns constant and T is temperature. The extracted band gap value for the 5-nm GNR was 489.419.0 meV (see the Supplementary Information for details of the tting process). Using similar methods, the band gap values for other narrow GNRs were also estimated. As shown in Fig. 4d, the band gaps, Eg, extracted from all measured GNRs were plotted with respect to the width of the corresponding GNRs. It is obvious that the band gap scaled inversely with ribbon width. The width dependence of the band gap t well with the function Eg (eV)Ba/(w b) (w is in units of

nm), where parameter aE1.99 eV nm and parameter b 1.28 nm.

DiscussionIn the GNR samples, the band gap exhibited a strong dependence on the width of the ribbons. In particular, the sub-10-nm GNRs exhibited an electronic band gap of B0.5 eV. It is important to understand the possible origins of the band gap. Normally, the electronic structure of GNRs is strongly dependent on the edge of the GNRs. For this study, it was challenging to characterize the exact geometry of the narrow ribbons obtained (for example, the edge chirality and termination state of the bonds). However, the GNRs were generally parallel to the ZZ pattern of the h-BN and free of detectable edge roughness. Therefore, it is likely that most segments of the GNR edges followed a ZZ pattern. Previous studies in the literature predicted that GNRs with pristine ZZ edges would always be metallic because of their peculiar at-band edge states localized near the Fermi level1. Recent theoretical papers have predicted that the at-band would split due to electronelectron interactions, followed by the opening of a band gap caused by antiferromagnetic coupling between the spins along opposite edges of the ZZ ribbons24,25. More recently, band gaps in narrow ZZ nanoribbons have been experimentally observed in scanning tunnelling microscopy and scanning tunnelling spectroscopy studies26,27. A similar high density of states at-band near the Fermi level was also observed in ZZ-terminated atomically sharp graphene BN interfaces28. We

believe that the gap opening in the GNRs embedded in h-BN may also have been due to ee interactions. It is noted that our narrow GNRs with similar widths exhibited relatively larger band gaps than those reported in previous literature9,24. The uniaxial strain29 from in-plane grapheneBN bonding and the Bernal stacking30 on h-BN may have provided additional contributions to the gap opening in the GNRs embedded in h-BN.

The applicability of graphene for future digital devices is often questioned due to its intrinsic gapless nature. Nanoribbons offer a potential solution, but both the width and edges must be precisely controlled. In-plane grapheneBN hetero-structured lms have been reported, however, only on metal surfaces3135. In addition, control over the dimensions of graphene and grapheneBN boundaries has not been fully achieved. By employing the in-plane epitaxy of graphene in nano-trenches of h-BN, we have realized ZZ-oriented GNRs with a controlled width and smooth edges. The GNRs feature a tunable band gap, enabling sub-10-nm GNR FETs with onoff ratios 4104. Our results demonstrate that it is possible to resolve the fundamental gapless limitation of graphene, paving the way for the realization of graphene-based digital electronics that can operate at room temperature.

Methods

Etching process on h-BN surface. To prepare h-BN samples for nano-particle-assisted etching, multiple cleaning steps were adapted to ensure the reliable production of the cuts. Before depositing the h-BN ake, quartz substrates were cleaned with acetone and isopropyl alcohol and then annealed at 600 C for 30 min to remove organic contaminants. Subsequently, the h-BN was mechanically exfoliated using semiconductor-grade tape. Next, the h-BN samples were heated in a quartz tube at 500 C for 15 min under an Ar:H2 ow (850:150 sccm) to remove the tape residue. After the heat cleaning process, a solution of NiCl2:H2O at a concentration of 0.01 mg ml 1 was spun at 1,800 r.p.m. for 100 s onto the substrate surface and then baked for 10 min at 80 C on a hot plate to evaporate the solvent. This NiCl2-treated sample was then submitted to a two-step process under an

Ar:H2 ow (850:150 sccm): annealing at 500 C for 20 min, which resulted in Ni nano-particle formation, and etching at 1,200 C for 60180 min. The optimal etching pressure was found to be B150 Pa.

Growth of GNRs. Before graphene growth, etched h-BN akes were cleaned separately in HCl solution, DI water and acetone. The substrate with the h-BN akes was then loaded into a growth chamber. GNR growth was carried out in a low-pressure CVD furnace at 1,280 C under an Ar ow of 10 sccm, corresponding to 15 Pa; the samples were then annealed for 5 min. Next, the Ar ow was turned off, and a C2H2 ow and a mixture of silane and Ar (5% mole ratio of silane to Ar)

were introduced into the system for GNR growth. The ratio of C2H2 to silane was B1:1. The pressure was maintained at 5 Pa during the growth process, and the growth time was B5 min. After growth, both the C2H2 and silane/Ar ow were turned off, and the system was cooled to room temperature with owing Ar. The samples of GNRs on h-BN grown on a quartz surface were moved to a highly p-type doped silicon wafer with a 300-nm-thick SiO2 capping layer for electrical transport studies.

Atomic force microscopy. Cleaned samples were characterized using one AFM (Dimension Icon, Bruker), while atomic-resolution images were taken by another AFM (Multimode IV, Veeco) under ambient conditions. AFM measurements were acquired in contact mode to obtain height and friction images. SNL-10 AFM tips from Bruker, which possess a nominal tip radius of o10 nm, were used in all measurements. The use of friction contrast was necessary because this mode gives clear information about the super-structure and atomic lattice. For atomic-resolution scanning, the force constant k of the cantilever tips was in the range of0.050.5 N m 1. The scan rate was set to 1060 Hz to minimize any noise from thermal drift. The integral gain and set point were adjusted to be as low as possible during the measurement. Several hours of pre-scanning were carried out to warm up the scanner to obtain good stability during imaging. To ensure a highly accurate atomic-resolution image, scanners with a travel range of o10 mm along the x and y directions were used. Calibration at atomic resolution was performed with newly cleaved highly ordered pyrolytic graphite before measurement.

Characterization of GNRs using Raman spectroscopy. Raman spectra were obtained with a commercially available confocal Raman instrument: model Alpha 300R from WITec. The Raman data were recorded using a laser wavelength of 488 nm oL2:54 eV

. An objective lens with 100 magnication and a 0.95

numerical aperture was used, producing a laser spot that was B500 nm in diameter. The laser power was maintained at o1 mW on the sample surface to avoid laser-induced heating.

Device fabrication and transport measurements. GNR devices were created by a standard electron beam lithographic technique with Ni or Pd as the source and drain contacts on p-doped silicon wafers with 300-nm-thick SiO2. Next, the devices were annealed in a hydrogen atmosphere at 200 C for 3 h to remove the resist residues and to reduce the contact resistance between the GNRs and metal electrodes before electrical measurements because the thickness of the h-BN akes on the silicon wafer was B15 nm. Electrical transport measurements were carried out using a physical property measurement system (PPMS from Quantum Design, Inc.) via a Keithley 4200 semiconductor characterization system.

Data availability. The authors declare that the main data supporting the ndings of this study are available within the article and its Supplementary Information les. Additional data are available from the corresponding author upon request.

References

1. Nakada, K., Fujita, M., Dresselhaus, G. & Dresselhaus, M. S. Edge state in graphene ribbons: nanometer size effect and edge shape dependence. Phys. Rev. B Condens. Matter 54, 1795417961 (1996).

2. Son, Y. W., Cohen, M. L. & Louie, S. G. Half-metallic graphene nanoribbons. Nature 444, 347349 (2006).

3. Novoselov, K. S. et al. A roadmap for graphene. Nature 490, 192200 (2012).4. Han, M. Y., Ozyilmaz, B., Zhang, Y. & Kim, P. Energy band gap engineering of graphene nanoribbons. Phys. Rev. Lett. 98, 206805 (2007).

NATURE COMMUNICATIONS | 8:14703 | DOI: 10.1038/ncomms14703 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14703

5. Chen, Z. H., Lin, Y. M., Rooks, M. J. & Avouris, P. Graphene nano-ribbon electronics. Phys. E 40, 228232 (2007).

6. Bai, J., Zhong, X., Jiang, S., Huang, Y. & Duan, X. Graphene nanomesh. Nat. Nanotechnol. 5, 190194 (2010).

7. Sprinkle, M. et al. Scalable templated growth of graphene nanoribbons on SiC. Nat Nanotechnol. 5, 727731 (2010).

8. Tapaszto, L., Dobrik, G., Lambin, P. & Biro, L. P. Tailoring the atomic structure of graphene nanoribbons by scanning tunnelling microscope lithography. Nat. Nanotechnol. 3, 397401 (2008).

9. Li, X., Wang, X., Zhang, L., Lee, S. & Dai, H. Chemically derived, ultrasmooth graphene nanoribbon semiconductors. Science 319, 12291232 (2008).

10. Wu, Z. S. et al. Efcient synthesis of graphene nanoribbons sonochemically cut from graphene sheets. Nano Res. 3, 1622 (2010).

11. Jiao, L., Zhang, L., Wang, X., Diankov, G. & Dai, H. Narrow graphene nanoribbons from carbon nanotubes. Nature 458, 877880 (2009).

12. Kosynkin, D. V. et al. Longitudinal unzipping of carbon nanotubes to form graphene nanoribbons. Nature 458, 872876 (2009).

13. Cai, J. et al. Atomically precise bottom-up fabrication of graphene nanoribbons. Nature 466, 470473 (2010).

14. Jacobberger, R. M. et al. Direct oriented growth of armchair graphene nanoribbons on germanium. Nat. Commun. 6, 8006 (2015).

15. Ago, H. et al. Lattice-oriented catalytic growth of graphene nanoribbons on heteroepitaxial nickel lms. ACS Nano 7, 1082510833 (2013).

16. Tang, S. et al. Precisely aligned graphene grown on hexagonal boron nitride by catalyst free chemical vapor deposition. Sci. Rep. 3, 2666 (2013).

17. Xue, J. et al. Scanning tunnelling microscopy and spectroscopy of ultra-at graphene on hexagonal boron nitride. Nat. Mater. 10, 282285 (2011).

18. Tang, S. et al. Silane-catalysed fast growth of large single-crystalline graphene on hexagonal boron nitride. Nat. Commun. 6, 6499 (2015).

19. Geick, R., Perry, C. H. & Rupprech, G. Normal modes in hexagonal boron nitride. Phys. Rev. 146, 543547 (1966).

20. Mohiuddin, T. M. G. et al. Uniaxial strain in graphene by Raman spectroscopy: G peak splitting, Grneisen parameters, and sample orientation. Phys. Rev. B 79, 205433 (2009).

21. Huang, M. et al. Phonon softening and crystallographic orientation of strained graphene studied by Raman spectroscopy. Proc. Natl Acad. Sci. USA 106, 73047308 (2009).

22. Neumann, C. et al. Raman spectroscopy as probe of nanometre-scale strain variations in graphene. Nat. Commun. 6, 8429 (2015).

23. Javey, A., Guo, J., Wang, Q., Lundstrom, M. & Dai, H. Ballistic carbon nanotube eld-effect transistors. Nature 424, 654657 (2003).

24. Son, Y. W., Cohen, M. L. & Louie, S. G. Energy gaps in graphene nanoribbons. Phys. Rev. Lett. 97, 216803 (2006).

25. Yang, L., Park, C. H., Son, Y. W., Cohen, M. L. & Louie, S. G. Quasiparticle energies and band gaps in graphene nanoribbons. Phys. Rev. Lett. 99, 186801 (2007).

26. Magda, G. Z. et al. Room-temperature magnetic order on zigzag edges of narrow graphene nanoribbons. Nature 514, 608611 (2014).

27. Li, Y., Chen, M., Weinert, M. & Li, L. Direct experimental determination of onset of electronelectron interactions in gap opening of zigzag graphene nanoribbons. Nat. Commun. 5, 4311 (2014).

28. Drost, R. et al. Electronic states at the graphenehexagonal boron nitride zigzag interface. Nano. Lett. 14, 51285132 (2014).

29. Pereira, V. M., Neto, A. H. C. & Peres, N. M. R. Tight-binding approach to uniaxial strain in graphene. Phys. Rev. B 80, 045401 (2009).

30. Hunt, B. et al. Massive Dirac fermions and Hofstadter buttery in a van der Waals heterostructure. Science 340, 14271430 (2013).

31. Ci, L. et al. Atomic layers of hybridized boron nitride and graphene domains. Nat. Mater. 9, 430435 (2010).

32. Liu, L. et al. Heteroepitaxial growth of two-dimensional hexagonal boron nitride templated by graphene edges. Science 343, 163167 (2014).

33. Levendorf, M. P. et al. Graphene and boron nitride lateral heterostructures for atomically thin circuitry. Nature 488, 627632 (2012).

34. Liu, Z. et al. In-plane heterostructures of graphene and hexagonal boron nitride with controlled domain sizes. Nat. Nanotechnol. 8, 119124 (2013).

35. Gong, Y. J. et al. Direct chemical conversion of graphene to boron- and nitrogen- and carbon-containing atomic layers. Nat. Commun. 5, 3193 (2014).

Acknowledgements

The work performed at the Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences, was partially supported by the National Science and Technology Major Projects of China (grant no. 2011ZX02707), the Chinese Academy of Sciences (Grant Nos. KGZD-EW-303 and XDB04040300), a project from the Science and Technology Commission of Shanghai Municipality (Grant No. 16ZR1442700) and the National Science Foundation of China (Grant No. 61136005). The work conducted at the Shanghai Institute of Technical Physics, Chinese Academy of Sciences, was partially supported by the National Science Foundation of China (Grant No. 91321311). T.Y. is grateful for the support of MOE2012-T2-049. H.W.and X.X. thank T. Taniguchi, K. Watanabe (National Institute for Materials Science, Japan) and J.H. Edgar (Kansas State University, USA) for supplying the partial h-BN crystals.

Author contributions

H.W. and X.X. directed the research work. H.W. conceived and designed the experiments. L.H. performed the etching processes on h-BN. L.C. performed the growth experiments for the GNRs. S.T. carried out both the etching and GNR growth during the beginning stage of this project. L.H., L.C., H.Xia and T.L. performed the AFM experiments. H.S.W., H.W. and H.Xie fabricated the electronic devices. H.S.W. and L.L. performed the transport measurements. H.W. and L.C. performed the calculations concerning the stretched Moir pattern. L.C., C.C. and T.Y. performed the Raman measurements. H.W., X.X., D.Z., L.D., T.W., H.S.W., C.C. and T.Y. analysed the experimental data and designed the gures. H.W., X.X., C.C., T.Y., D.Z., L.D. andM.J. co-wrote the manuscript, and all authors contributed to critical discussions of the manuscript.

Additional information

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

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

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

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Chen, L. et al. Oriented graphene nanoribbons embedded in hexagonal boron nitride trenches. Nat. Commun. 8, 14703 doi: 10.1038/ncomms14703 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

6 NATURE COMMUNICATIONS | 8:14703 | DOI: 10.1038/ncomms14703 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications