Skip to main content
PMC home page
Nature Portfolio logoLink to Nature Portfolio
. 2015 Sep 11;6:8119. doi: 10.1038/ncomms9119

The positive piezoconductive effect in graphene

Kang Xu 1,*, Ke Wang 2,3,*, Wei Zhao 1,*, Wenzhong Bao 4,5, Erfu Liu 1, Yafei Ren 2,3, Miao Wang 1, Yajun Fu 1, Junwen Zeng 1, Zhaoguo Li 1, Wei Zhou 1, Fengqi Song 1, Xinran Wang 6, Yi Shi 6, Xiangang Wan 1, Michael S Fuhrer 4,7, Baigeng Wang 1,c, Zhenhua Qiao 2,3,b, Feng Miao 1,a, Dingyu Xing 1
PMCID: PMC4579395  PMID: 26360786

Abstract

As the thinnest conductive and elastic material, graphene is expected to play a crucial role in post-Moore era. Besides applications on electronic devices, graphene has shown great potential for nano-electromechanical systems. While interlayer interactions play a key role in modifying the electronic structures of layered materials, no attention has been given to their impact on electromechanical properties. Here we report the positive piezoconductive effect observed in suspended bi- and multi-layer graphene. The effect is highly layer number dependent and shows the most pronounced response for tri-layer graphene. The effect, and its dependence on the layer number, can be understood as resulting from the strain-induced competition between interlayer coupling and intralayer transport, as confirmed by the numerical calculations based on the non-equilibrium Green's function method. Our results enrich the understanding of graphene and point to layer number as a powerful tool for tuning the electromechanical properties of graphene for future applications.

Inline graphicThe thinness and elastic properties of graphene make it an ideal candidate for nanoelectromechanical systems. Here, the authors demonstrate the positive piezoconductive effect in suspended bi- and multi-layer graphene, showing the most pronounced response for three layer samples.

Graphene1,2,3,4 is an ideal material candidate for nano-electromechanical systems5 due to many advantageous features, including unparalleled breaking strength6, ultrahigh carrier mobility7 and excellent controllability of electronic structures via mechanical strain8,9. Many intriguing phenomena have been experimentally observed on strained graphene10,11,12, including the observation of pseudo-magnetic fields exceeding 300 T (ref. 13). More fascinating phenomena have been theoretically predicted for strained graphene, but yet to be realized experimentally, such as the zero-field quantum Hall effect14, strain-induced superconductivity of graphene15 and potential applications in valleytronics16,17. Thus, strain-controllable transport measurements are critical for in-depth understanding and further applications of graphene.

In situ piezoconductive measurements on graphene provide an effective approach to study the correlations between electrical properties and mechanical strains. So far, studies have been focused on mono-layer graphene, and the negative piezoconductive effect has been widely reported, independent of the differences in graphene film synthesis, transfer methods and sample substrates18,19,20,21,22,23. Systematic studies of the electromechanical properties of graphene with different number of layers and alterable interlayer interactions have been lacking.

Here we investigate the piezoconductive effect of suspended graphene membranes with various layer numbers by applying in situ stress with a scanning probe. We observe positive piezoconductance in bi- and multi-layer graphene, with tri-layer graphene showing the most pronounced response. This intriguing phenomenon can be explained by the model of strain-induced competition between electronic interlayer coupling and intralayer transport, and further confirmed by numerical calculations based on non-equilibrium Green's function method.

Results

Device fabrication and piezoconductive measurements

Suspended graphene membranes are mechanically exfoliated and deposited on Si/SiO2 wafers with pre-etched trenches. The number of the graphene layers is identified via colour interference, and confirmed by Raman spectroscopy. Metal electrodes (5 nm Ti/50 nm Au) are deposited through home-made shadow masks, which effectively avoid wet process-induced device performance degradation24. A typical device image is shown in Fig. 1a. To perform in situ piezoconductive measurements, we introduce pressure-modulated conductance microscopy25,26, which utilizes a non-conducting atomic force microscopy (AFM) tip to apply adjustable local pressure on the top of the suspended graphene, yielding a topography image of the strained graphene. The conductance/resistance of the device is monitored simultaneously, and comparison of the conductance/resistance image and topography offers piezoconductive information of the suspended graphene membranes. The detailed experimental setup is schematically shown in Fig. 1b.

Figure 1. Suspended graphene device and PCM setup.

Figure 1

Open in a new tab

(a) Optical microscope image of a four-terminal suspended bilayer graphene device, which is fabricated by home-made shadow mask method. Inset: SEM image of a suspended device. Scale bar, 3 μm. (b) Schematic setup of pressure-modulated conductance microscopy (PCM) that performs piezoconductive measurements on suspended graphene devices.

The measurements are carried out on mono- and multi-layer suspended graphene devices. Typical piezoconductive results from mono-, bi-, tri- and tetra-layer devices are shown in Fig. 2a–d, respectively. For each figure, the left panel shows a line trace of the topography image and the right panel shows the corresponding relative conductance change represented by g. Here g(x) is defined byInline graphic, where G(x) is the device conductance at AFM tip position x (where local pressure is applied) and G0 is the undisturbed conductance value. For mono-layer suspended graphene, the device conductance drops upon local pressure applied (Fig. 2a), indicating negative piezoconductive (that is, positive piezoresistive) effect. This is consistent with all previous studies19,20. However, for multi-layer suspended graphene, device conductance jumps upon local pressure applied (Fig. 2b–d), indicating positive piezoconductive effect, which has never been reported on graphene devices. We note that the measured tri-layer graphene device is stacked with the common Bernal (ABA) structure, as confirmed by Raman spectroscopy.

Figure 2. Layer number-dependent positive piezoconductive effect of graphene.

Figure 2

Open in a new tab

(a–d) Left panels are the line traces of the topography (tip position) image of graphene devices from mono-layer (a) to tetra-layer (d) with applied strain ɛ=0.54‰, 0.51‰, 0.33‰ and 0.20‰, respectively. Right panels are the line traces of the corresponding relative conductance change g (to the undisturbed conductance with no local pressure applied). Mono-layer device shows negative piezoconductive effect (conductance drops upon local pressure applied). Bi-, tri- and tetra-layer devices show positive piezoconductive effect (conductance jumps upon local pressure applied) with the most pronounced effect in tri-layer device. (e) The plot of the maximum relative conductance change (when AFM tip approaches the centre of the suspended membranes) gmax as a function of strain ɛ for various suspended graphene devices (layer number, n=1, 2, 3, 4 and 6 represented by different colours).

To further study the observed layer number-dependent piezoconductive effect, we perform the same measurement on various suspended graphene devices (layer number, n=1, 2, 3, 4 and 6) with different strains (up to 1‰). The detailed results are shown in Fig. 2e. Here we plot the maximum relative conductance change gmax (which usually appears when AFM tip approaches the centre of the suspended membranes). According to the geometry of our devices, the strain ɛ can be calculated by the equation

graphic file with name ncomms9119-m2.jpg

where h is the maximum strain-induced deflection and l is the length of the suspended graphene (same as the width of the trench, around 3 μm). Here h has been corrected by subtracting the height of the tip at which the force begins to rise (which can be extracted from the deflection-force curve, see Supplementary Fig. 1 and Supplementary Note 1 for details). As shown in the gmax versus ɛ plot in Fig. 2e, several data points from a mono-layer graphene device fall in the negative regime. The gauge factor (usually defined by relative resistance change divided by strain) is estimated to be about 0.6, similar to the values previously reported19,20. In sharp contrast, all multi-layer graphene devices show positive conductance response. More interestingly, for similar strain value, the n=3 (tri-layer) device shows much larger gmax than the other multi-layer devices (n=2, 4 and 6).

Theoretical interpretation of the underlying physical origin

The positive piezoconductive effect is difficult to explain in terms of strain-induced decrease of Fermi velocity, the model that has been applied to strained mono-layer graphene19,20. In the system of back-gated suspended mono-layer graphene, a model of inhomogeneous carrier density redistribution predicted a positive piezoconductive effect27 in contrast to our observation of negative piezoconductance for mono-layer graphene, and cannot explain our observation as well. The fact that positive piezoconductance is only observed in bi- and multi-layer graphene suggests that interlayer coupling is key; indeed, strain should modify the interaction between graphene layers, modifying the band structure and hence transport properties.

To understand the observed positive piezoconductive effect, we numerically calculate the transport properties of strained multi-layer graphene devices, which can be calculated by applying non-equilibrium Green's function technique and using two-terminal Landauer–Büttiker formula, with details described in the Methods section. To numerically study the piezoconductive effect in the multi-layer graphene systems in the presence of an external pressure, we consider a π-orbital tight-binding model Hamiltonian, which is written as4,28,29:

graphic file with name ncomms9119-m3.jpg
graphic file with name ncomms9119-m4.jpg
graphic file with name ncomms9119-m5.jpg
graphic file with name ncomms9119-m6.jpg

where Inline graphic and ci are the creation and annihilation operators on the site i. The first term H0 describes the pristine multi-layer graphene sheets, with μ being the chemical potential and tintra,inter denoting the intra- and interlayer nearest-neighbour hopping strength, respectively. In our consideration, we take tintra =2.60 eV and tinter =0.34 eV, respectively. The second term HD describes the influences of external disorders and the dephasing effect that is used to recover the macroscopic behaviour from a finite-sized quantum system and can be modelled by attaching individual virtual leads at each site i. Here ɛi is the on-site Anderson disorder strength that is uniformly distributed in the interval of (−W/2, W/2) with W characterizing the strength of disorder and Inline graphic and aik are the creation and annihilation operators for the virtual lead attached at the i-th site. In our consideration, we set the disorder strength to be 0.1 eV, because the prepared multi-layer graphene sheets are rather clean and it can better smear off the room temperature-induced thermal fluctuation. The last term HS represents the applied strain-induced effects including the induced site potentials and the variation of the intra- and interlayer nearest-neighbour hopping energies in the lattice-deformed region. The AFM tip-contacted lattices should experience noticeable lattice deformations, while the lattice deformation farther away from the AFM tip gradually decreases. Thus it is reasonable to assume that the variation of the intra- and interlayer hopping energy due to the lattice deformation takes the form of Inline graphic, where r0 is the centre of deformed region and L is the system length of the disordered scattering region, and the varying site potential takes the similar form of Inline graphic. In our calculations, the conductance is averaged over 1,000 ensembles at each point.

We first focus on the tri-layer graphene that shows the most pronounced positive piezoconductive effect. As illustrated in Fig. 3a, the tri-layer structure is constructed by mono-layer graphene in a Bernal ABA-stacking manner. Its conductance includes contributions from electron transport in both horizontal and vertical directions, and gmax is determined by the strain-induced competition between the intralayer hopping and interlayer interactions. In the presence of an applied pressure from the AFM tip, the direct consequences on the multi-layer graphene include: (1) the extension of the lattice constant in the lattice-deformed area that correspondingly decreases the intralayer nearest-neighbour hopping; and (2) the local compression of the interlayer separation in the lattice-deformed area that effectively modulates the interlayer interactions, that is, increasing the interlayer hopping and inducing the local site energies. By applying the non-equilibrium Green's function technique and Landauer–Büttiker formula (see details in the Methods section) on a tri-layer graphene ribbon with length of 100 nm and width of 13 nm, we theoretically study the relation between the maximum relative conductance change gmax and the applied strain ɛ for tri-layer graphene. It is noteworthy to mention that all our obtained data are used to qualitatively explain and understand the underlying physical origin of the positive piezoconductive effect, and for the better presentation we have chosen certain parameters to compare with the experimental observations. Our numerical results show that gmax is always positive and increases with increasing strain, as plotted in the blue cycles in Fig. 3b. The experimental data are plotted in the same graph, showing the similar tendency. The nonlinear feature indicates for higher strain gmax increases slower, and suggests the relative contribution from interlayer modulation becomes less dominant. This model also explains why the piezoconductive property of bi- and multi-layer graphene is distinct from mono-layer graphene. For strained mono-layer graphene, negative peizoconductance is observed because no interlayer modulation but only intralayer contribution is present.

Figure 3. Theoretical model and numerical calculations of tri-layer graphene.

Figure 3

Open in a new tab

(a) Schematic of the lattice structure change of Bernal (ABA) stacking tri-layer graphene due to a vertical load F applied. (b) The maximum relative conductance change gmax as a function of strain ɛ for tri-layer graphene. The red squares are the experimental data while the blue cycles are the numerical simulation results.

Discussion

To further understand the intriguing phenomenon that tri-layer graphene shows the most pronounced positive piezoconductive effect, we have performed first-principle calculations to study the lattice configuration of different multi-layer graphene devices under strain. Our first-principle calculations are performed using the projected-augmented-wave method as implemented in the Vienna ab initio simulation package30 with details described in the Methods section. Figure 4a shows how the local lattice configuration is altered by the same local strain in different multi-layer graphene systems. For the bilayer graphene, one can observe that both the top and bottom layers exhibit the structural deformations with different amplitudes. We further plot the strain-induced lattice variation Δd between the nearest two layers as a function of the layer number n in Fig. 4b. Along with the increase of the layer number, for example, from bilayer to up to hexa-layer, only the top three layers show noticeable structural deformations. Therefore, the positive piezoconductive effects for multi-layer graphene systems due to the interlayer modulations should be dominated by the top three layers.

Figure 4. Simulation results of the dependence on layer number.

Figure 4

Open in a new tab

(a) The structurally relaxed configuration for different multi-layer graphene (layer number, n=2–6 from left to right) in the presence of the same strain strength. (b) The strain-induced lattice variation Δd between the nearest two layers as a function of the layer number n. (c) The dependence of the piezoconductive factor γ on layer number n, with experimental data represented by red squares and simulation results represented by blue squares. The error bars of experimental data originate from the fitting process, while the error bars of simulating results originate from different disorder configurations in the numerical calculations.

We can now introduce a characteristic factor, the piezoconductive factor γ, to describe the piezoconductive properties of bi- and multi-layer graphene devices:

graphic file with name ncomms9119-m12.jpg

The piezoconductive factor γ has the same sign as piezoconductive effect. The square root dependence on strain is phenomenological but describes the experimental and theoretical data well and allows us to parameterize each gmax(ɛ) curve by a single-value γ. By applying the same non-equilibrium Green's function technique as used in studying the tri-layer graphene, we calculate the gmax–ɛ relation for graphene with various layer number n (n=2 to 6) (see detailed parameters in Supplementary Table 1), and extract the values of γ by using the above definition. The results are shown by the blue squares in Fig. 4c. The values of γ are positive, with the maximum value around 1.6 for tri-layer graphene. To compare the results with experimental observation, we plot the experimental data in the same graph (red squares) and find they are consistent with theoretical results. For bi- and tri-layer graphene, γ increases with n, suggesting the enhancement of the contribution from strain-induced interlayer modulation. But while n increases further, γ decreases, suggesting the interlayer modulation is dominated by the top three layers and additional parallel conduction paths from the extra layers suppress the positive piezoconductive effect. Here we note that, the roles of external effects (see Supplementary Figs 2–5 and Supplementary Note 2 for details), the stacking order (Supplementary Figs 6–8 and Supplementary Note 3) and back gate (Supplementary Fig. 9 and Supplementary Note 4), as well as finite size effect in the theoretical calculations (Supplementary Figs 10–11 and Supplementary Note 5), have been carefully explored in our work. Although the observation of negative piezoconductive effect in multi-layer graphene was reported previously31, it is likely related to the devices' unusually large contact resistance, rather than the intrinsic piezoconductive properties of multi-layer graphene as revealed in this manuscript.

In summary, we have systematically studied the piezoconductive effect of suspended graphene with different number of layers. In contrast to the negative piezoconductance observed in mono-layer graphene, positive piezoconductance is observed in bi- and multi-layer graphene, with tri-layer graphene showing the most pronounced response. It can be explained by the model of strain-induced competition between electronic interlayer coupling and intralayer transport. This model is further confirmed by numerical calculations based on non-equilibrium Green's function method. Our results enrich the understanding of electromechanical properties of graphene and underscore their potential applications on the field of nano-electromechanical systems and flexible electronics.

Methods

Suspended graphene device preparation and characterization

The suspended graphene membranes are obtained by using mechanical exfoliation method on the pre-defined trenches on 300-nm thick SiO2 wafers. The trenches are defined by standard photolithography method, followed by dry etching in an Inductively Coupled Plasma (ICP) system, where CH4 and CHF3 are used as etching gases. The typical width of the trenches is 3 μm and the depth is around 250 nm. The number of layers of graphene membranes is first identified by an optical microscopy and further confirmed by Raman spectroscopy. To avoid common wet process-induced device performance degradation and yield loss, the electronic devices are fabricated by using a home-made shadow mask method24. The electrodes are made of 5-nm Ti covered by 50-nm Au.

Setup of pressure-modulated conductance microscopy

A Bruker Multimode 8 AFM is used to build up the pressure-modulated conductance microscopy setup. A non-conducting AFM tip is used to apply adjustable local pressure on the top of the suspended graphene, yielding topography/height image of the strained graphene. The conductance/resistance of the device is measured via lock-in technique and monitored simultaneously as an external input of the AFM. The comparison of the conductance/resistance image and topography image offers detailed piezoconductive information of the suspended graphene membranes. The scanning procedure is done in contact mode with a Bruker NP-S10 tip, of which the Young's modulus is 0.24 N m−1. The detailed experimental setup is schematically shown in Fig. 1b.

Details of non-equilibrium Green's function simulation

Using multi-probe Landauer–Büttiker formula, the current in the lead p (either real or virtual lead) can be expressed as:

graphic file with name ncomms9119-m13.jpg

where Vp/q is the spin-independent bias in the lead p/q. The electronic transmission coefficient from lead q to lead p is calculated as Tpq = Tr[ΓpGrΓqGa], in which the line-width function Γp is defined as Inline graphic, and the retarded and advanced Green's function are given by Inline graphic, where I is the unit matrix with the same dimension as that of H. The retarded and advanced self-energy due to the coupling to all the real leads can be obtained numerically32. For the virtual leads, we assume Inline graphic and the dephasing strength Γd is fixed to 0.01 eV. In our simulations, a small external bias is applied between the left and right lead with VL=−VR=0.5 V. For dephasing effect, electrons lose the phase memory by entering and leaving the virtual leads. Thus, for each virtual lead i, the current has the constraint that Ji=0, which ensures the current conservation. Combining the above equation of Jp together with all boundary conditions for the real and virtual leads, the voltage Vp and current Jp in each real lead can be obtained.

Details of first-principle calculation

The generalized gradient approximation combined the Vander Waals correction with the DFT-D2 method of Grimme is used. The kinetic energy cutoff is set to be 500 eV. During the structure relaxation, the edge atoms and the probed atoms are not allowed to relax while the others are. All parameters are chosen to converge the forces to be <0.01 eV Å−1. The first Brillouin-zone integration is carried out by using the 3 × 3 × 1 Gamma-centred grids. A vacuum buffer space of 20 Å is set to prevent the interaction between adjacent slabs.

Additional information

How to cite this article: Xu, K. et al. The positive piezoconductive effect in graphene. Nat. Commun. 6:8119 doi: 10.1038/ncomms9119 (2015).

Supplementary Material

Supplementary Information

Supplementary Figures 1-11, Supplementary Table 1, Supplementary Notes 1-5 and Supplementary References

ncomms9119-s1.pdf (1.2MB, pdf)

Acknowledgments

We thank Chun Ning Lau and Pingheng Tan for stimulating discussions. This work was supported in part by the National Key Basic Research Program of China (2015CB921600, 2011CB922103, 2013CBA01603), the National Natural Science Foundation of China (11374142, 61574076), the Natural Science Foundation of Jiangsu Province (BK20130544, BK20140017, BK20150055), the Specialized Research Fund for the Doctoral Program of Higher Education (20130091120040), Fundamental Research Funds for the Central Universities and the Collaborative Innovation Center of Advanced Microstructures. K.W., Y.R. and Z.Q. are supported by 100 Talents Program of Chinese Academy of Sciences, NNSFC (11474265) and Anhui Provincial Natural Science Foundation, The Supercomputing Center of USTC is gratefully acknowledged for the high-performance computing assistance.

Footnotes

Author contributions F.M. conceived the project and designed the experiments; K.X., W.B., E.L., M.W., Y.F. and J.Z. carried out device fabrication and electrical measurements; K.W., W.Zhao, Y.R. and Z.Q. carried out numerical simulations; Z.L. and W.Zhou carried out Raman spectroscopy measurements and analysis; F.M., M.S.F., B.W., Z.Q., K.X., K.W. and W.Zhao carried out data analysis and interpretation; F.M., M.S.F., Z.Q., K.X. and K.W. co-wrote the paper with all authors contributing to the discussion and preparation of the manuscript.

References

  1. Novoselov K. S. et al. Electric field effect in atomically thin carbon films. Science 306, 666–669 (2004). [DOI] [PubMed] [Google Scholar]
  2. Novoselov K. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197–200 (2005). [DOI] [PubMed] [Google Scholar]
  3. Zhang Y., Tan Y., Stormer H. L. & Kim P. Experimental observation of the quantum Hall effect and Berry's phase in graphene. Nature 438, 201–204 (2005). [DOI] [PubMed] [Google Scholar]
  4. Neto A. C., Guinea F., Peres N., Novoselov K. S. & Geim A. K. The electronic properties of graphene. Rev. Mod. Phys. 81, 109–162 (2009). [Google Scholar]
  5. Arden W. et al. “More-than-Moore” International Technical Roadmap for Semiconductors (2010). [Google Scholar]
  6. Lee C., Wei X., Kysar J. W. & Hone J. Measurement of the elastic properties and intrinsic strength of monolayer graphene. Science 321, 385–388 (2008). [DOI] [PubMed] [Google Scholar]
  7. Bolotin K. I. et al. Ultrahigh electron mobility in suspended graphene. Solid State Commun. 146, 351–355 (2008). [Google Scholar]
  8. Ni Z. H. et al. Uniaxial strain on graphene: Raman spectroscopy study and band-gap opening. ACS Nano 2, 2301–2305 (2008). [DOI] [PubMed] [Google Scholar]
  9. Huang M., Yan H., Heinz T. F. & Hone J. Probing strain-induced electronic structure change in graphene by Raman spectroscopy. Nano Lett. 10, 4074–4079 (2010). [DOI] [PubMed] [Google Scholar]
  10. Teague M. L. et al. Evidence for strain-induced local conductance modulations in single-layer graphene on SiO2. Nano Lett. 9, 2542–2546 (2009). [DOI] [PubMed] [Google Scholar]
  11. Zhu S., Ghatkesar M. K., Zhang C. & Janssen G. Graphene based piezoresistive pressure sensor. Appl. Phys. Lett. 102, 161904 (2013). [Google Scholar]
  12. Bao W. et al. Controlled ripple texturing of suspended graphene and ultrathin graphite membranes. Nat. Nanotechnol. 4, 562–566 (2009). [DOI] [PubMed] [Google Scholar]
  13. Levy N. et al. Strain-induced pseudo–magnetic fields greater than 300 tesla in graphene nanobubbles. Science 329, 544–547 (2010). [DOI] [PubMed] [Google Scholar]
  14. Guinea F., Katsnelson M. I. & Geim A. K. Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering. Nat. Phys. 6, 30–33 (2009). [Google Scholar]
  15. Si C., Liu Z., Duan W. & Liu F. First-principles calculations on the effect of doping and biaxial tensile strain on electron-phonon coupling in graphene. Phys. Rev. Lett. 111, 196802 (2013). [DOI] [PubMed] [Google Scholar]
  16. Wu Z., Zhai F., Peeters F. M., Xu H. Q. & Chang K. Valley-dependent Brewster angles and Goos-Hänchen effect in strained graphene. Phys. Rev. Lett. 106, 176802 (2011). [DOI] [PubMed] [Google Scholar]
  17. Jiang Y., Low T., Chang K., Katsnelson M. I. & Guinea F. Generation of pure bulk valley current in graphene. Phys. Rev. Lett. 110, 046601 (2013). [DOI] [PubMed] [Google Scholar]
  18. Kim K. S. et al. Large-scale pattern growth of graphene films for stretchable transparent electrodes. Nature 457, 706–710 (2009). [DOI] [PubMed] [Google Scholar]
  19. Huang M., Pascal T. A., Kim H., Goddard W. A. III & Greer J. R. Electronic−mechanical coupling in graphene from in situ nanoindentation experiments and multiscale atomistic simulations. Nano Lett. 11, 1241–1246 (2011). [DOI] [PubMed] [Google Scholar]
  20. Smith A. D. et al. Electromechanical piezoresistive sensing in suspended graphene membranes. Nano Lett. 13, 3237–3242 (2013). [DOI] [PubMed] [Google Scholar]
  21. Wang Y. et al. Super-elastic graphene ripples for flexible strain sensors. ACS Nano 5, 3645–3650 (2011). [DOI] [PubMed] [Google Scholar]
  22. Fu X. et al. Strain dependent resistance in chemical vapor deposition grown graphene. Appl. Phys. Lett. 99, 213107 (2011). [Google Scholar]
  23. Lee Y. et al. Wafer-scale synthesis and transfer of graphene films. Nano Lett. 10, 490–493 (2010). [DOI] [PubMed] [Google Scholar]
  24. Bao W. et al. Lithography-free fabrication of high quality substrate-supported and freestanding graphene devices. Nano Res. 3, 98–102 (2010). [Google Scholar]
  25. Lau C. N., Stewart D. R., Williams R. S. & Bockrath M. Direct observation of nanoscale switching centers in metal/molecule/metal structures. Nano Lett. 4, 569–572 (2004). [Google Scholar]
  26. Miao F., Ohlberg D., Stewart D. R., Williams R. S. & Lau C. N. Quantum conductance oscillations in metal/molecule/metal switches at room temperature. Phys. Rev. Lett. 101, 016802 (2008). [DOI] [PubMed] [Google Scholar]
  27. Medvedyeva M. V. & Blanter Y. M. Piezoconductivity of gated suspended graphene. Phys. Rev. B 83, 045426 (2011). [Google Scholar]
  28. Xing Y., Sun Q. & Wang J. Influence of dephasing on the quantum Hall effect and the spin Hall effect. Phys. Rev. B 77, 115346 (2008). [Google Scholar]
  29. Jiang H., Cheng S., Sun Q. & Xie X. C. Topological insulator: a new quantized spin hall resistance robust to dephasing. Phys. Rev. Lett. 103, 036803 (2009). [DOI] [PubMed] [Google Scholar]
  30. Liu F., Ming P. & Li J. Ab initio calculation of ideal strength and phonon instability of graphene under tension. Phys. Rev. B 76, 064120 (2007). [Google Scholar]
  31. Zheng X., Chen X., Kim J., Lee D. & Li X. Measurement of the gauge factor of few-layer graphene. MEMS MOEMS 12, 013009 (2013). [Google Scholar]
  32. Lee D. H. & Joannopoulos J. D. Simple scheme for surface-band calculations. II. The Green's function. Phys. Rev. B 23, 4997 (1981). [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplementary Information

Supplementary Figures 1-11, Supplementary Table 1, Supplementary Notes 1-5 and Supplementary References

ncomms9119-s1.pdf (1.2MB, pdf)

Articles from Nature Communications are provided here courtesy of Nature Publishing Group

RESOURCES