The development of a scalable and cost-effective nanofabrication method is of key importance for future advances in nanoelectronics. Thermal scanning probe lithography (t-SPL) is a growing nanopatterning method with potential for parallelization, offering unique capabilities that make it an attractive candidate for industrial nanomanufacturing. Here, we demonstrate the possibility to apply t-SPL for the fabrication of graphene devices. In particular, we use t-SPL to produce high performing graphene-based field effect transistors (FETs). The here described t-SPL process includes the fabrication of high-quality metal contacts, as well as patterning and etching of graphene to define the active region of the device. The electrical measurements on the t-SPL fabricated FETs indicate a symmetric conductance at the Dirac point and a low specific contact resistance without the use of any contact engineering strategy. The entire t-SPL nanofabrication process is performed without the need for masks, and in ambient conditions. Furthermore, thanks to the t-SPL in situ simultaneous patterning and imaging capability, no markers are required. These features substantially decrease fabrication time and cost.

An important objective of the microelectronics industry is to fabricate high-performance miniaturized devices on a large scale at a low cost and with a high throughput. The exploration of new materials, fabrication methods, and device architectures underpins this objective. Of particular interest in recent years has been the study of electronic devices where the active material is monolayer graphene, due to its intriguing properties, including remarkable carrier transport,1 high thermal conductivity,2 and outstanding mechanical stability.3 These properties have made graphene a promising candidate for a wide range of applications from flexible electronics to high-speed electronics.

The fabrication of graphene functional devices generally requires multiple patterning steps for defining the device active region and patterning metal contacts on graphene. Currently, electron-beam lithography (EBL) is the dominant technique for prototyping high-performance nanoscale devices out of graphene,4–8 producing better quality devices than photolithography (PL).9–11 However, the use of a focused electron beam in EBL is a major limitation in terms of cost, the electron beam induced sample damage and potential for parallelization to increase the throughput.

In recent years, novel nano-patterning techniques have been explored for fabricating nanoscale devices.12–15 Among those, thermal scanning probe lithography (t-SPL) is an attractive choice for multiple reasons.12,16–20 First, t-SPL is a maskless technique and is capable of patterning nanoscale features with sub-10 nm resolution.17,21 Second, the entire t-SPL nano-patterning process can take place in atmospheric condition or N2, which is a considerable advantage for achieving a cost-effective nano-patterning process. Third, t-SPL has a throughput comparable to EBL when using only one probe (∼105µm2 h−1),22 but by multiplexing with thermal nanoprobe arrays,23 it could reach a much larger throughput. Finally, and more importantly, a recent study has demonstrated that t-SPL can pattern high-performing and low-resistance metal contacts on monolayer MoS2.24 However, the application of t-SPL for the fabrication of graphene devices remains largely unexplored. The objective of this study is, hence, to study the effectiveness of t-SPL for fabricating graphene devices.

Here, we report the application of t-SPL performed in ambient conditions for the fabrication of graphene field-effect transistors (GFETs). In particular, we study the electrical properties of monolayer graphene using a bottom-gated device structure with no encapsulation layer. In this device structure, monolayer graphene in both channel and contact regions is in direct contact with the t-SPL resist during all the lithographic steps. Despite this, and without the use of contact engineering (see the supplementary material for details), we find that the GFETs fabricated using t-SPL exhibit a relatively low specific contact resistance of 600 Ω·μm. This value is within the typical range of specific contact resistance (100 Ω·μm–800 Ω·μm) for monolayer graphene devices built using EBL, including those employing contact engineering.4,8,25–28

Our experiments began with the mechanical exfoliation of monolayer graphene from a bulk graphite crystal. Graphene monolayers were exfoliated onto a heavily doped silicon substrate covered with 285 nm of thermally grown SiO2, which serves as the global back-gate in the final device structure. Monolayer graphene flakes were identified under an optical microscope and then verified using Raman spectroscopy. We then annealed the sample in an Ar/H2 ambient at 500 °C for 1 h. This annealing step removes the residual tape contaminants on the substrate. The cleanliness of the oxide surface is important for spin-coating the substrate with a uniform resist layer for performing nano-patterning using t-SPL. As we describe below, we used a commercial t-SPL system for all nano-patterning steps in the device fabrication process.

The first step in our device fabrication process was to define and etch the active region within the monolayer graphene flake. Figure 1 shows the schematic illustration and optical images of this processing step. For t-SPL nano-patterning, we used a two-polymer stack resist, consisting of a 210 nm thick PMGI (polymethylglutarimide) layer and a 15 nm thick thermosensitive polymer film, namely, PPA (polyphthalaldehyde). These layers were deposited sequentially onto the substrate through spin coating [Figs. 1(a) and 1(f)]. We then performed t-SPL to define the active region of interest. During the nano-patterning process, a heated nano-probe (typically heated to 200 °C at the probe-resist contact) thermally decomposes and evaporates PPA. The precise movement of the probe transfers a computer-generated pattern into the PPA film [Fig. 1(b)]. The decomposed PPA quickly evaporates without being re-deposited onto the surface of the sample. The pattern was then chemically etched into the underlying PMGI using a diluted TMAH (tetramethylammonium hydroxide) solution, as shown in Figs. 1(c) and 1(g), which exposes the unwanted monolayer graphene regions for removal. We used a brief oxygen plasma treatment (see the experimental method) for removing the exposed graphene regions, resulting in a rectangular active device region [Figs. 1(e) and 1(h)]. The PPA/PMGI stack resist during the oxygen plasma etching process serves as a hard mask for protecting the underneath monolayer graphene in the active region. Hence, it must be sufficiently thick to survive the plasma etch process. However, patterning high-resolution features using t-SPL requires the use of a thin PPA film (see Fig. S1 of the supplementary material), suggesting that the PMGI film must be made thick enough to protect the active device region. In Fig. 1(i), we show how the thickness of PPA and PMGI changes due to exposure to the oxygen plasma. For this experiment, we produced two sample groups. Each sample group consisted of multiple substrates coated with either PPA or PMGI. The initial thickness of the film within each sample group was identical. However, each substrate within the group was subjected to a different etching time. Each data point in Fig. 1(i) represents one substrate. The data show that a thin PPA film (with ∼12 nm initial thickness) withstood less than 15 s of exposure to oxygen plasma. In contrast, only 12 nm of the PMGI film was consumed after 25 s of etching under the same conditions, which is adequately long for etching monolayer graphene. This experiment explains the rationale for choosing the above-mentioned thicknesses for PMGI and PPA when patterning the active region of the graphene device. The resolution of the current pattern-etch transfer process is discussed in the supplementary material, where nanoscale etched graphene nano-ribbons are presented (Figs. S1 and S2).

FIG. 1.

Optimization of the t-SPL nano-patterning process for defining graphene regions. (a)–(e) Schematic illustration of the t-SPL process for patterning graphene active regions. (f) Optical microscope image of the starting graphene flake after coating with the PPA/PMGI resist. The green dotted box shows the target active region. (g) The same graphene flake after t-SPL patterning and chemical etching of the PMGI layer. (h) Final rectangular graphene ribbon obtained after 25 s oxygen plasma etching. The graphene ribbon has a length and width of 11.5 µm and 2.5 µm, respectively. Scale bars are 10 µm. (i) Etching rates of PPA and PMGI with oxygen plasma. Each data point represents one new substrate.

FIG. 1.

Optimization of the t-SPL nano-patterning process for defining graphene regions. (a)–(e) Schematic illustration of the t-SPL process for patterning graphene active regions. (f) Optical microscope image of the starting graphene flake after coating with the PPA/PMGI resist. The green dotted box shows the target active region. (g) The same graphene flake after t-SPL patterning and chemical etching of the PMGI layer. (h) Final rectangular graphene ribbon obtained after 25 s oxygen plasma etching. The graphene ribbon has a length and width of 11.5 µm and 2.5 µm, respectively. Scale bars are 10 µm. (i) Etching rates of PPA and PMGI with oxygen plasma. Each data point represents one new substrate.

Close modal

The second step in device fabrication after defining the active region on monolayer graphene was to pattern the metal electrodes. Figures 2(a)2(e) show schematic illustrations of the fabrication steps. First, a two-polymer stack of PPA/PMGI (15 nm PPA/210 nm PMGI) was spin-coated on monolayer graphene. Then, we used a heated t-SPL nano-probe to pattern the metal electrode regions in the PPA film. Subsequently, the patterns were transferred into the underlying PMGI through chemical etching in diluted TMAH [see Figs. 2(b) and 2(c)], which also produces the required undercut. In this process, the top thin PPA ensures high-resolution patterning, while the underlying PMGI layer eases the metal lift-off process. Finally, we deposited a stack of Cr/Au (10 nm/20 nm) metals using electron-beam evaporation, followed by the metal lift-off. Figure 2(h) shows an example of the optical image of the final graphene device structure. To demonstrate the capability of t-SPL to fabricate graphene field-effect transistors (GFETs) at the nanoscale, we fabricated metal contacts on graphene with a minimum channel length of 60 nm [Fig. S2(c)]. It is noteworthy that such a small spacing between metal contacts results from the pattern amplification due to the development of PMGI and partial non-line-of-sight metal deposition. With proper optimization of pattern amplification, PPA/PMGI thickness, wet etching duration, RIE etch conditions, probe size, pattern depth, and eventually using a different process without PMGI, sub-10 nm graphene nanoribbons and channel lengths can be obtained using the process described in this work.

FIG. 2.

t-SPL metal electrode patterning on graphene. (a)–(e) Schematic illustrations of the t-SPL patterning process for the fabrication of metal electrodes for GFETs. (f) in situ t-SPL imaging of monolayer graphene (rectangular ribbon with a length and width of 80 µm and 6.7 µm, respectively) after spin-coating the PPA/PMGI resist. (g) in situ t-SPL imaging of the structure, showing the patterned electrode features in the PPA layer. (h) Example of the optical image of a back-gated graphene device after lift-off. The spacings between electrodes are 0.6 µm, 2.3 µm, 4.3 µm, 6.1 µm, and 8.2 µm. Scale bars are 10 µm.

FIG. 2.

t-SPL metal electrode patterning on graphene. (a)–(e) Schematic illustrations of the t-SPL patterning process for the fabrication of metal electrodes for GFETs. (f) in situ t-SPL imaging of monolayer graphene (rectangular ribbon with a length and width of 80 µm and 6.7 µm, respectively) after spin-coating the PPA/PMGI resist. (g) in situ t-SPL imaging of the structure, showing the patterned electrode features in the PPA layer. (h) Example of the optical image of a back-gated graphene device after lift-off. The spacings between electrodes are 0.6 µm, 2.3 µm, 4.3 µm, 6.1 µm, and 8.2 µm. Scale bars are 10 µm.

Close modal

The t-SPL technique provides in situ simultaneous patterning and imaging, a capability that distinguishes t-SPL from other fabrication methods such as EBL or PL and provides important benefits, including the fact that there is no need for alignment marks.29,30 Furthermore, the here-shown graphene devices have been fabricated in the ambient environment without the need for ultra-high vacuum (UHV). Indeed, before nano-patterning, a cold t-SPL probe produces a thermal image of the monolayer graphene active region underneath the resist to determine the target patterning location [Fig. 2(f)]. A second benefit of in situ imaging is the ability to inspect the quality of the patterned features simultaneously with t-SPL nano-patterning, allowing for a true closed feedback loop.12Figure 2(g) shows an in situ image of the patterned metal electrodes in the PPA layer.

Next, we studied the electronic properties of the t-SPL fabricated GFETs at room temperature. Figure 3(a) shows the total resistance, Rtot, of three graphene devices as a function of the applied back-gate bias, indicating the increase in the device resistance with the channel length. The total resistance in Fig. 3(a) was calculated by taking the ratio of the drain–source voltage (Vds) to the current (Id). These devices [marked as 1–3 in the inset optical image in Fig. 3(a)] have different channel lengths. Note that the shorter length devices were not connected, possibly due to cracks in their channel regions. The Dirac point voltage (VDirac) for channels 1, 2, and 3 in Fig. 3(a) are 1.2 V, −1.4 V, and −5.6 V, respectively (Fig. S3). These Dirac voltage shifts correspond to residual carrier densities of 9.1 × 1010 cm−2, 10.6 × 1010 cm−2, and 4.2 × 1011 cm−2, which indicate the preservation of the intrinsic state of graphene after the t-SPL GEFT fabrication process.

FIG. 3.

Electrical characteristics of t-SPL fabricated devices. (a) Transfer characteristics of three GFETs in a TLM structure with the same channel width (3.5 µm) but different channel lengths. The contact width is 2.4 µm. The graphene ribbon has a full length of ∼50 µm. The measurements were made at Vds = 50 mV. (b) Plot of Rtot against the channel length for three different carrier densities for the GFETs in panel (a). The solid lines are linear fit to the data. The y-intercepts give 2Rc. (c) Extracted specific contact resistance vs the carrier density. The green data points are the extracted ρc and the red bars represent the error. The extracted ρc values near CNP (the gray shading region) have considerable error. (d) Intrinsic channel resistivity plotted against the carrier density. The inset shows the double-logarithmic plot of conductivity vs carrier density, giving an estimated upper bound for the residual carrier density n*. The black solid lines represent the extrapolated fits.

FIG. 3.

Electrical characteristics of t-SPL fabricated devices. (a) Transfer characteristics of three GFETs in a TLM structure with the same channel width (3.5 µm) but different channel lengths. The contact width is 2.4 µm. The graphene ribbon has a full length of ∼50 µm. The measurements were made at Vds = 50 mV. (b) Plot of Rtot against the channel length for three different carrier densities for the GFETs in panel (a). The solid lines are linear fit to the data. The y-intercepts give 2Rc. (c) Extracted specific contact resistance vs the carrier density. The green data points are the extracted ρc and the red bars represent the error. The extracted ρc values near CNP (the gray shading region) have considerable error. (d) Intrinsic channel resistivity plotted against the carrier density. The inset shows the double-logarithmic plot of conductivity vs carrier density, giving an estimated upper bound for the residual carrier density n*. The black solid lines represent the extrapolated fits.

Close modal

To analyze the electronic properties of these devices, we first extracted the contact resistance using the transmission-line method (TLM). In particular, the total resistance is the sum of the intrinsic channel resistance and the contact resistance, given by

(1)

where Rc is the contact resistance due to one electrode, ρch is the channel resistivity of monolayer graphene, and L and W are the channel length and width, respectively. In this equation, Rtot scales linearly with L, whereas the y-intercept is 2Rc. Moreover, the slope gives information about the intrinsic resistivity of graphene, which depends on the carrier concentration (n) in the graphene channel. In Fig. 3(b), we plot the total resistance of devices 1, 2, and 3 at three different electron carrier densities, giving an estimated total contact resistance (2Rc) of 504 Ω, 403 Ω, and 355 Ω at n = 1.1 × 1012 cm−2, 1.9 × 1012 cm−2, and 2.6 × 1012 cm−2, respectively. Note that we calculated the carrier density n using the following31 equation:

(2)

where Cox is the oxide capacitance (1.2 × 10−8 F/cm2), is the reduced Planck constant, υF is Fermi velocity of graphene (1 × 106 m/s), and kF=nπ is the Fermi wave vector. Figure 3(c) shows the extracted specific contact resistance, ρc, plotted against the carrier density for the range of the applied gate bias. The gray shading in this plot marks the region where uncertainty in the intercept of the linear fits to Rtot data is considerably large (>26% of the extracted Rc). This region corresponds to low carrier densities near the charge neutrality point (CNP). At high carrier densities, however, the error in estimated Rc is small, and ρc is as low as 600 Ω·μm in the electron branch. These values indicate that the here fabricated contacts to monolayer graphene are of good quality, considering that no supplementary contact engineering strategies were used. Table I compares ρc of our graphene devices fabricated by t-SPL in ambient conditions with some of the best results reported in literature for graphene FET fabricated by EBL in UHV and by photolithography with and without different supplementary contact engineering methods. We remark that similar methods of contact engineering can also be implemented together with the t-SPL process to achieve similar reductions of ρc (Table I in the supplementary material).

TABLE I.

Summary of state-of-art specific contact resistance ρc obtained by different fabrication methods.

ContactContact Specific contact Carrier density Gate
GrapheneLithographymaterialengineeringresistance (Ω·μm)(1012 cm−2)structureReferences
CVD PL Ti/Au UV ozone cleaning 200 … Global back-gate 38  
CVD PL Cr/Au CO2 cleaning 270 … Global back-gate 39  
None 960 … 
CVD EBL Pd/Au None 654 Global back-gate 40  
MoO3 doping 200 70 
CVD EBL Pd/Au High-purity Pd and high-vacuum deposition 100 3.75 Global back-gate 25  
Exfoliated EBL Cr/Pd/Au Edge contact 150 Global back-gate 41  
CVD EBL Au Holey contact 45 20 Global back-gate 27  
None 519 20 
Exfoliated EBL Ni Ni-etched contact 100 … Global back-gate 28  
None 630 
CVD EBL Ti/Pd/Au Double contact 260 4.75 Global back-gate 42  
Exfoliated t-SPL Cr/Au None 600 Global back-gate This work 
ContactContact Specific contact Carrier density Gate
GrapheneLithographymaterialengineeringresistance (Ω·μm)(1012 cm−2)structureReferences
CVD PL Ti/Au UV ozone cleaning 200 … Global back-gate 38  
CVD PL Cr/Au CO2 cleaning 270 … Global back-gate 39  
None 960 … 
CVD EBL Pd/Au None 654 Global back-gate 40  
MoO3 doping 200 70 
CVD EBL Pd/Au High-purity Pd and high-vacuum deposition 100 3.75 Global back-gate 25  
Exfoliated EBL Cr/Pd/Au Edge contact 150 Global back-gate 41  
CVD EBL Au Holey contact 45 20 Global back-gate 27  
None 519 20 
Exfoliated EBL Ni Ni-etched contact 100 … Global back-gate 28  
None 630 
CVD EBL Ti/Pd/Au Double contact 260 4.75 Global back-gate 42  
Exfoliated t-SPL Cr/Au None 600 Global back-gate This work 

Finally, we calculated the carrier mobility in the as-fabricated monolayer graphene FET. To do so, we extracted the channel resistivity of device 1 from the slope of the linear regressions in Fig. 3(b) at different carrier densities. Figure 3(d) shows the plot of ρch against the carrier density in graphene. The carrier mobility can be estimated by fitting the channel resistivity using ρch=neμL+σ01+ρs, where μL represents the mobility due to long-range scattering, σ0 is the minimum conductivity at CNP, and ρs indicates the contribution from short-range scattering.32,33 Fitting the data in Fig. 3(d) (pink solid curve) yields μL-Fit of ∼4500 cm2/V s and ∼6100 cm2/V s for the electron branch and the hole branch, respectively. These carrier mobilities are within the range of previously reported values for graphene on SiO2, which is between 2000 cm2/V s and 20 000 cm2/V s.34–37 

We also employed the theoretical graphene transport study by Adam et al.43 for evaluating the expected carrier mobility in the diffusive limit from the impurity concentration (nimp) at the graphene–oxide interface. In particular, the transport mobility follows μ = 20e/(hnimp), where nimp can be estimated from the residual carrier density (n*) using n* ≅ 0.3nimp. The plot in the inset of Fig. 3(d) shows the double-logarithmic plot of the conductivity vs carrier density, which gives an upper bound estimate of n* = 2.6 × 1011 cm−2 from the σ plateau. Using this information, we calculate carrier mobility of ∼5500 cm2/V s. The obtained mobility from the analytical solution is comparable with the fitting results, providing further confidence in the extracted specific contact resistance in Fig. 3(c).

In summary, we demonstrated the application of t-SPL for fabricating GFETs. We showed that it is possible to combine t-SPL with plasma etching for producing graphene structures of desired shapes. More importantly, the t-SPL fabricated metal electrodes resulted in low-resistance contacts on monolayer graphene without contact engineering. Further improvements of the contact resistance require an in-depth study to identify the factors that limit the contact resistance in the t-SPL process. t-SPL is very attractive compared to EBL because it offers in situ simultaneous imaging and patterning capabilities and it operates in ambient conditions, which is a considerable advantage for achieving a cost-effective nano-patterning process. The results presented here establish the prospects of t-SPL for the fabrication of graphene devices.

See the supplementary material for the following analyses and/or descriptions: the high resolution t-SPL pattern in the PPA resist, nanoscale etch-pattern transfer and metal deposition using the t-SPL process, unshifted transfer characteristics of GFET in Fig. 3(a), and the literature review of contact engineering of metal-graphene contacts.

X.L. and Z.H. contributed equally to this work.

We acknowledge the support of the National Science Foundation (Grant Nos. NSF CBET – 1914539 and CMMI 1914540) and the U.S. Army Research Office. The experiments were performed with the NanoFrazor, acquired through Grant No. NSF CMMI MRI – 1929453. This work was partially supported by the MRSEC Program of the National Science Foundation under Award No. DMR-1420073. This work was partially performed at the ASRC NanoFabrication Facility of CUNY in New York.

The data that support the findings of this study are available within the article.

1. Graphene sample preparation

Graphene layers were identified under an optical microscope and confirmed by Raman spectroscopy. Graphene flakes were exfoliated on 285 nm SiO2 on Si substrates from bulk graphite crystals (NGS Naturgraphit) using a Scotch tape. Monolayer graphene flakes were identified under an optical microscope and confirmed by Raman spectroscopy.

2. Defining graphene active region

Graphene samples were spin-coated with PMGI (polymethylglutarimide, Microchem) (2000 rpm for 35 s) followed by baking at 200 °C for 1 min. This step was repeated three times resulting in a film of 215 nm. Then, a PPA (polyphthalaldehyde, Sigma Aldrich) solution (0.9 wt. % in anisole) is spin-coated on PMGI (2000 rpm for 4 s and then 3000 rpm for 35 s) followed by baking at 90 °C for 3 min. This gives rise to a PPA/PMGI (15 nm/215 nm) polymer stack. The graphene active region was patterned on the PPA resist using a commercial t-SPL system, Nanofrazor (Heidelberg Instruments). After the t-SPL patterning, samples were immersed in a solution of TMAH in deionized water (tetramethylammonium hydroxide AZ726 MIF, MicroChemicals) (0.17 mol/L) for 400 s to chemically etch the exposed PMGI in the patterned region, then rinsed with deionized water (30 s) and IPrOH (30 s), and finally dried with N2. Then, samples were etched with a mild O2 reactive ion etching for 25 s (20 W, 100 mTorr; Oxford Instruments PlasmaPro NPG80 RIE) and transferred into the Remover PG (MicroChem) for a few hours to strip off the polymer resist, followed by rinsing (IPrOH) and drying (N2).

3. Fabrication of GFETs by t-SPL

PMGI and PPA were spin-coated with the same conditions as described for defining the graphene active region. The patterns of electrical contacts were generated by t-SPL in the PPA resist over the previously defined graphene active region. The same chemical etching by TMAH was also performed to expose the graphene in the contact regions, which generates an undercut profile to facilitate the metal lift-off. Metal deposition is performed using an AJA Orion 8E e-beam evaporator (deposition rate: 1 Å s−1 and pressure ∼10−8 Torr). Finally, samples were dipped in the Remover PG (MicroChem) for a few hours to lift off the metal and resist, followed by rinsing (IPrOH) and drying (N2). The metal electrodes of back-gated t-SPL GFETs presented in this paper have been fabricated using Cr/Au (10 nm/20 nm).

4. Electrical measurements

The electrical measurements of the TLM-structured GFETs were made inside a Lakeshore probe station at a pressure of 10−5 Torr and using a Keithley 4200-SCS parameter analyzer.

5. AFM measurements

All the thicknesses of polymer resists and metals were measured by AFM (Bruker Multimode 8) operating in the tapping mode.

1.
K. S.
Novoselov
,
A. K.
Geim
,
S. V.
Morozov
,
D.
Jiang
,
M. I.
Katsnelson
,
I. V.
Grigorieva
,
S. V.
Dubonos
, and
A. A.
Firsov
,
Nature
438
,
197
(
2005
).
2.
A. A.
Balandin
,
S.
Ghosh
,
W.
Bao
,
I.
Calizo
,
D.
Teweldebrhan
,
F.
Miao
, and
C. N.
Lau
,
Nano Lett.
8
,
902
(
2008
).
3.
C.
Lee
,
X.
Wei
,
J. W.
Kysar
, and
J.
Hone
,
Science
321
,
385
(
2008
).
4.
L.
Anzi
,
A.
Mansouri
,
P.
Pedrinazzi
,
E.
Guerriero
,
M.
Fiocco
,
A.
Pesquera
,
A.
Centeno
,
A.
Zurutuza
,
A.
Behnam
,
E. A.
Carrion
,
E.
Pop
, and
R.
Sordan
,
2D Mater.
5
,
025014
(
2018
).
5.
L.
Banszerus
,
M.
Schmitz
,
S.
Engels
,
M.
Goldsche
,
K.
Watanabe
,
T.
Taniguchi
,
B.
Beschoten
, and
C.
Stampfer
,
Nano Lett.
16
,
1387
(
2016
).
6.
A. S.
Mayorov
,
R. V.
Gorbachev
,
S. V.
Morozov
,
L.
Britnell
,
R.
Jalil
,
L. A.
Ponomarenko
,
P.
Blake
,
K. S.
Novoselov
,
K.
Watanabe
,
T.
Taniguchi
, and
A. K.
Geim
,
Nano Lett.
11
,
2396
(
2011
).
7.
Y.-C.
Lin
,
C.-C.
Lu
,
C.-H.
Yeh
,
C.
Jin
,
K.
Suenaga
, and
P.-W.
Chiu
,
Nano Lett.
12
,
414
(
2012
).
8.
F.
Xia
,
V.
Perebeinos
,
Y.-m.
Lin
,
Y.
Wu
, and
P.
Avouris
,
Nat. Nanotechnol.
6
,
179
(
2011
).
9.
J.
Fan
,
J. M.
Michalik
,
L.
Casado
,
S.
Roddaro
,
M. R.
Ibarra
, and
J. M.
De Teresa
,
Solid State Commun.
151
,
1574
(
2011
).
10.
C. A.
Chavarin
,
A. A.
Sagade
,
D.
Neumaier
,
G.
Bacher
, and
W.
Mertin
,
Appl. Phys. A
122
,
58
(
2016
).
11.
T.
Cusati
,
G.
Fiori
,
A.
Gahoi
,
V.
Passi
,
M. C.
Lemme
,
A.
Fortunelli
, and
G.
Iannaccone
,
Sci. Rep.
7
,
5
109
(
2017
).
12.
R.
Garcia
,
A. W.
Knoll
, and
E.
Riedo
,
Nat. Nanotechnol.
9
,
577
(
2014
).
13.
J. G.
Son
,
M.
Son
,
K.-J.
Moon
,
B. H.
Lee
,
J.-M.
Myoung
,
M. S.
Strano
,
M.-H.
Ham
, and
C. A.
Ross
,
Adv. Mater.
25
,
4723
(
2013
).
14.
A. N.
Abbas
,
G.
Liu
,
B.
Liu
,
L.
Zhang
,
H.
Liu
,
D.
Ohlberg
,
W.
Wu
, and
C.
Zhou
,
ACS Nano
8
,
1538
(
2014
).
15.
M. F.
El-Kady
and
R. B.
Kaner
,
ACS Nano
8
,
8725
(
2014
).
16.
R.
Szoszkiewicz
,
T.
Okada
,
S. C.
Jones
,
T.-D.
Li
,
W. P.
King
,
S. R.
Marder
, and
E.
Riedo
,
Nano Lett.
7
,
1064
(
2007
).
17.
X.
Liu
,
M.
Kumar
,
A.
Calo
,
E.
Albisetti
,
X.
Zheng
,
K. B.
Manning
,
E.
Elacqua
,
M.
Weck
,
R. V.
Ulijn
, and
E.
Riedo
,
ACS Appl. Mater. Interfaces
11
,
41780
(
2019
).
18.
Z.
Wei
,
D.
Wang
,
S.
Kim
,
S.-Y.
Kim
,
Y.
Hu
,
M. K.
Yakes
,
A. R.
Laracuente
,
Z.
Dai
,
S. R.
Marder
,
C.
Berger
,
W. P.
King
,
W. A.
de Heer
,
P. E.
Sheehan
, and
E.
Riedo
,
Science
328
,
1373
(
2010
).
19.
S. T.
Howell
,
A.
Grushina
,
F.
Holzner
, and
J.
Brugger
,
Microsyst. Nanoeng.
6
,
21
(
2020
).
20.
X.
Zheng
,
A.
Calo
,
T.
Cao
,
X.
Liu
,
Z.
Huang
,
P. M.
Das
,
M.
Drndic
,
E.
Albisetti
,
F.
Lavini
,
T.-D.
Li
,
V.
Narang
,
W. P.
King
,
J. W.
Harrold
,
M.
Vittadello
,
C.
Aruta
,
D.
Shahrjerdi
, and
E.
Riedo
,
Nat. Commun.
11
,
3463
(
2020
).
21.
Y. K. R.
Cho
,
C. D.
Rawlings
,
H.
Wolf
,
M.
Spieser
,
S.
Bisig
,
S.
Reidt
,
M.
Sousa
,
S. R.
Khanal
,
T. D. B.
Jacobs
, and
A. W.
Knoll
,
ACS Nano
11
,
11890
(
2017
).
22.
P. C.
Paul
,
A. W.
Knoll
,
F.
Holzner
,
M.
Despont
, and
U.
Duerig
,
Nanotechnology
22
,
275306
(
2011
).
23.
K. M.
Carroll
,
X.
Lu
,
S.
Kim
,
Y.
Gao
,
H.-J.
Kim
,
S.
Somnath
,
L.
Polloni
,
R.
Sordan
,
W. P.
King
,
J. E.
Curtis
, and
E.
Riedo
,
Nanoscale
6
,
1299
(
2014
).
24.
X.
Zheng
,
A.
Calò
,
E.
Albisetti
,
X.
Liu
,
A. S. M.
Alharbi
,
G.
Arefe
,
X.
Liu
,
M.
Spieser
,
W. J.
Yoo
,
T.
Taniguchi
,
K.
Watanabe
,
C.
Aruta
,
A.
Ciarrocchi
,
A.
Kis
,
B. S.
Lee
,
M.
Lipson
,
J.
Hone
,
D.
Shahrjerdi
, and
E.
Riedo
,
Nat. Electron.
2
,
17
(
2019
).
25.
H.
Zhong
,
Z.
Zhang
,
B.
Chen
,
H.
Xu
,
D.
Yu
,
L.
Huang
, and
L.
Peng
,
Nano Res.
8
,
1669
(
2015
).
26.
F.
Giubileo
and
A.
Di Bartolomeo
,
Prog. Surf. Sci.
92
,
143
(
2017
).
27.
V.
Passi
,
A.
Gahoi
,
E. G.
Marin
,
T.
Cusati
,
A.
Fortunelli
,
G.
Iannaccone
,
G.
Fiori
, and
M. C.
Lemme
,
Adv. Mater. Interfaces
6
,
1801285
(
2019
).
28.
W. S.
Leong
,
H.
Gong
, and
J. T. L.
Thong
,
ACS Nano
8
,
994
(
2014
).
29.
C.
Rawlings
,
H.
Wolf
,
J. L.
Hedrick
,
D. J.
Coady
,
U.
Duerig
, and
A. W.
Knoll
,
ACS Nano
9
,
6188
(
2015
).
30.
C.
Rawlings
,
U.
Duerig
,
J.
Hedrick
,
D.
Coady
, and
A. W.
Knoll
,
IEEE Trans. Nanotechnol.
13
,
1204
(
2014
).
31.
S.
Kim
,
J.
Nah
,
I.
Jo
,
D.
Shahrjerdi
,
L.
Colombo
,
Z.
Yao
,
E.
Tutuc
, and
S. K.
Banerjee
,
Appl. Phys. Lett.
94
,
062107
(
2009
).
32.
C. R.
Dean
,
A. F.
Young
,
I.
Meric
,
C.
Lee
,
L.
Wang
,
S.
Sorgenfrei
,
K.
Watanabe
,
T.
Taniguchi
,
P.
Kim
,
K. L.
Shepard
, and
J.
Hone
,
Nat. Nanotechnol.
5
,
722
(
2010
).
33.
L.
Banszerus
,
M.
Schmitz
,
S.
Engels
,
J.
Dauber
,
M.
Oellers
,
F.
Haupt
,
K.
Watanabe
,
T.
Taniguchi
,
B.
Beschoten
, and
C.
Stampfer
,
Sci. Adv.
1
,
e1500222
(
2015
).
34.
K.
Nagashio
,
T.
Yamashita
,
T.
Nishimura
,
K.
Kita
, and
A.
Toriumi
,
J. Appl. Phys.
110
,
024513
(
2011
).
35.
F.
Schwierz
,
Nat. Nanotechnol.
5
,
487
(
2010
).
36.
Y. W.
Tan
,
Y.
Zhang
,
K.
Bolotin
,
Y.
Zhao
,
S.
Adam
,
E. H.
Hwang
,
S.
Das Sarma
,
H. L.
Stormer
, and
P.
Kim
,
Phys. Rev. Lett.
99
,
246803
(
2007
).
37.
C.
Jozsa
,
M.
Popinciuc
,
N.
Tombros
,
H. T.
Jonkman
, and
B. J.
van Wees
,
Phys. Rev. Lett.
100
,
236603
(
2008
).
38.
W.
Li
,
Y.
Liang
,
D.
Yu
,
L.
Peng
,
K. P.
Pernstich
,
T.
Shen
,
A. R.
Hight Walker
,
G.
Cheng
,
C. A.
Hacker
,
C. A.
Richter
,
Q.
Li
,
D. J.
Gundlach
, and
X.
Liang
,
Appl. Phys. Lett.
102
,
183110
(
2013
).
39.
S.
Gahng
,
C.
Ho Ra
,
Y.
Jin Cho
,
J.
Ah Kim
,
T.
Kim
, and
W.
Jong Yoo
,
Appl. Phys. Lett.
104
,
223110
(
2014
).
40.
S.
Vaziri
,
V.
Chen
,
L.
Cai
,
Y.
Jiang
,
M. E.
Chen
,
R. W.
Grady
,
X.
Zheng
, and
E.
Pop
,
IEEE Electron Device Lett.
41
,
1592
(
2020
).
41.
L.
Wang
,
I.
Meric
,
P. Y.
Huang
,
Q.
Gao
,
Y.
Gao
,
H.
Tran
,
T.
Taniguchi
,
K.
Watanabe
,
L. M.
Campos
,
D. A.
Muller
,
J.
Guo
,
P.
Kim
,
J.
Hone
,
K. L.
Shepard
, and
C. R.
Dean
,
Science
342
,
614
(
2013
).
42.
A. D.
Franklin
,
S.-J.
Han
,
A. A.
Bol
, and
V.
Perebeinos
,
IEEE Electron Device Lett.
33
,
17
(
2011
).
43.
S.
Adam
,
E. H.
Hwang
,
V. M.
Galitski
, and
S.
Das Sarma
,
Proc. Natl. Acad. Sci. U. S. A.
104
,
18392
(
2007
).

Supplementary Material