Strain-Induced Pseudo–Magnetic Fields Greater Than 300 Tesla in Graphene Nanobubbles

See allHide authors and affiliations

Science  30 Jul 2010:
Vol. 329, Issue 5991, pp. 544-547
DOI: 10.1126/science.1191700


Recent theoretical proposals suggest that strain can be used to engineer graphene electronic states through the creation of a pseudo–magnetic field. This effect is unique to graphene because of its massless Dirac fermion-like band structure and particular lattice symmetry (C3v). Here, we present experimental spectroscopic measurements by scanning tunneling microscopy of highly strained nanobubbles that form when graphene is grown on a platinum (111) surface. The nanobubbles exhibit Landau levels that form in the presence of strain-induced pseudo–magnetic fields greater than 300 tesla. This demonstration of enormous pseudo–magnetic fields opens the door to both the study of charge carriers in previously inaccessible high magnetic field regimes and deliberate mechanical control over electronic structure in graphene or so-called “strain engineering.”

Graphene, a single atomic layer of carbon, displays remarkable electronic and mechanical properties (1, 2). Many of graphene’s distinctive properties arise from a linear band dispersion at low carrier energies (3) that leads to Dirac-like behavior within the two-dimensional (2D) honeycomb lattice—charge carriers travel as if their effective mass is zero (1). An intriguing recent prediction is that a distortion of the graphene lattice should create large, nearly uniform pseudo–magnetic fields and give rise to a pseudo–quantum Hall effect (4). Whereas an elastic strain can be expected to induce a shift in the Dirac point energy from local changes in electron density, it is also predicted to induce an effective vector potential that arises from changes in the electron-hopping amplitude between carbon atoms (5). This strain-induced gauge field can give rise to large pseudo–magnetic fields (Bs) for appropriately selected geometries of the applied strain (1, 6). In such situations, the charge carriers in graphene are expected to circulate as if under the influence of an applied out-of-plane magnetic field (710). It has recently been proposed that a modest strain field with triangular symmetry will give approximately uniform, quantizing Bs upward of tens of tesla (4).

Here, we report the measurement of Landau levels (LLs) arising from giant strain-induced pseudo–magnetic fields in highly strained graphene nanobubbles grown on the Pt(111) surface. Landau quantization of the electronic spectrum was observed by scanning tunneling microscopy (STM), which revealed pseudo–magnetic fields in excess of 300 T. Such enormous strain-induced pseudo–magnetic fields may allow the electronic properties of graphene to be controlled through various schemes for applying strain (11), as well as the exploration of new high-field physical regimes.

Strained graphene nanobubbles were created by in situ growth of sub-monolayer graphene films in ultrahigh vacuum on a clean Pt(111) surface (12) in order to avoid external contamination and trapped gases. The epitaxial graphene was grown by exposure of Pt(111) to ethylene followed by annealing (13, 14). Graphene grown on Pt is expected to be minimally coupled to the substrate, compared to graphene grown on other catalytic metals (15, 16). A Dirac-like band structure is preserved for graphene on Pt(111), as verified by a recent photoemission study (17). An STM topograph of the graphene/Pt(111) surface prepared in this manner (Fig. 1A) reveals a flat graphene patch (partially surrounded by Pt) that encompasses five graphene nanobubbles. Graphene nanobubbles frequently appear near the edges of a graphene patch, but are also sometimes observed in the center of flat patches or near the boundaries between patches and are presumably pinned near these locations (Fig. 1A). These nanobubbles are likely related to the larger-scale “wrinkle” structures observed by low energy electron microscopy that form upon cooling as a result of the differing thermal expansion coefficients of graphene and the platinum surface (17).

Fig. 1

STM images and STS spectra taken at 7.5 K. (A) Graphene monolayer patch on Pt(111) with four nanobubbles at the graphene-Pt border and one in the patch interior. Unreacted ethylene molecules and a small hexagonal graphene patch can be seen in the lower right (Itunneling = 50 pA, Vsample = 350 mV, 3D z-scale enhanced 4.6×). (Inset) High-resolution image of a graphene nanobubble showing distorted honeycomb lattice resulting from strain in the bubble (Itunneling = 50 pA, Vsample = 200 mV, max z = 1.6 nm, 3D z-scale enhanced 2×). (B) STS spectra of bare Pt(111), flat graphene on Pt(111) (shifted upward by 3 × 10−11 ohm−1), and the center of a graphene bubble (shifted upward by 9 × 10−11 ohm−1). Vmod = 20 mV.

Individual nanobubbles often have a triangular shape (Fig. 1A, inset), reflecting the lattice symmetry of the graphene and the underlying Pt surface, and are typically 4 to 10 nm across and 0.3 to 2.0 nm tall. Atomic-resolution imaging of the nanobubbles confirms the honeycomb structure of graphene here (Fig. 1A, inset), although the lattice is distorted because of the large strain occurring in these structures. The expected strain-induced pseudo–magnetic field in a graphene nanobubble can be estimated by using the relation Φ=(βh2/la)Φ0 for the flux per ripple in a distorted graphene sheet (6), where h is the height, l is the width, a is on the order of the C-C bond length, and Φ0 is the quantum of flux. The parameter β=log(t)/log(a) relates the change in the hopping amplitude between nearest neighbor carbon atoms (t) to bond length and has a typical magnitude of 2<β<3 for graphene. For a nanobubble of l = 4 nm and h = 0.5 nm, this yields a Bs of order 100 T. The large curvature and correspondingly high strain incorporated into the triangular nanoscale bubbles observed here make them ideal candidates for the observation of pseudo-LL because of large strain-induced pseudo–magnetic fields.

The local electronic structure of strained graphene nanobubbles and surrounding graphene films was characterized by scanning tunneling spectroscopy (STS) performed at ~ 7.5 K by using standard lock-in techniques to obtain differential conductance (dI/dV). The measurement of dI/dV reflects features in the local density of states (LDOS) of the surface at the position of the STM tip (18). STS measurement of the bare Pt surface was used to calibrate the LDOS of the tip upon approach and between sequences of spectra. STS spectra measured over the bare Pt regions (Fig. 1B) are relatively featureless and show the expected Pt(111) surface state (19). Spectra recorded over the flat graphene patches show a subtly modified structure compared with the clean Pt(111) surface, and no clear signatures of the graphene Dirac point were observed in these regions (Fig. 1B). Spectra measured at the boundary between the flat graphene and the nanobubbles (fig. S2) exhibit features consistent with a Dirac point located ~300 mV above the Fermi energy, as recently observed by photoemission (17), as well as a gaplike feature with a full width at half maximum (FWHM) of 127 ± 9 mV centered at the Fermi energy (Vsample = 0) recently associated with phonon-assisted inelastic tunneling (20).

STS measurements made directly over the nanobubble regions exhibit a succession of relatively strong peaks, spaced by more than 100 meV, which do not appear in the spectra on other regions of the sample (Fig. 1B). The peaks are typically weaker at negative bias, which may be attributed to the expected shorter vertical extension of wavefunctions at lower energies. These peaks in the LDOS of the graphene nanobubble are unlike features seen previously in STS performed on graphene in the absence of a magnetic field on SiC (21, 22) and SiO2 (20, 23) substrates and often overwhelm the usual graphene features. Figure 2A shows a series of spectra taken at different positions across a single nanobubble (topography shown in Fig. 2B). These spectra display the typical nanobubble peak structure as well as the inelastic feature at the Fermi energy described above [a broader bias range is shown in (fig. S3)]. Similar peak structure was observed on 10 different nanobubbles with four different STM tips, with some variation in peak spacing and amplitude that we presume was caused by variations in strain-induced electronic structure arising from different nanobubble geometries.

Fig. 2

(A) Sequence of eight dI/dV spectra (T ~ 7.5 K, Vmod = 20 mV) taken in a line across a graphene nanobubble shown in the image in (B). Red lines are data with quartic background subtracted; black dotted lines are Lorentzian peak fits (center of peaks indicated by dots, with blue dots indicating n = 0). Vertical dash-dot lines follow the energy progression of each peak order. (C) Normalized peak energy versus Embedded Image for peaks observed on five different nanobubbles follow expected scaling behavior from Eq. 1 (dashed line).

These peaks observed in the nanobubble LDOS can be attributed to LLs originating from a strain-induced pseudo–magnetic field. Other possible origins that we rule out as unlikely include confinement effects and defect creation. Confinement of Dirac fermions in graphene is difficult because of Klein tunnelling and the suppression of backscattering (1). If confinement somehow occurred at the nanobubble edges, the peaks observed would be expected to follow a different progression, most notably missing the n = 0 peak observed in spectra taken at the center of most nanobubbles (24). Confinement would also result in strong nodal patterns in LDOS (24), which were not observed. The production of defect states by strain is also unlikely because of the high energetic barrier toward defect creation. It is possible that some nanobubbles were formed at defect sites in the as-grown graphene, but the nanobubbles that we measured did not show signatures of defect physics in spectroscopy, such as tip height dependence of peak positions associated with the charging of defect states.

The most likely explanation for the nanobubble peaks is that they arise from a large, relatively uniform strain-induced pseudo–magnetic field. This pseudo–magnetic field is expected to mimic the influence of a real magnetic field applied perpendicular to the graphene sheet and give rise to LLs (4). These appear as a series of peaks in STS as they do for the case of a real magnetic field (25, 26). Specifically, the 2D massless Dirac nature of charge carriers near the Fermi energy causes the progression of peaks in the LDOS to follow the expression (7, 8):En=sgn(n)ωc|n|+EDirac,ωc=2evF2Bs (1)where En is the position in energy of the nth LL with respect to the Fermi level, ωc is the cyclotron resonance frequency arising from Bs, and vF is the Fermi velocity. Both positive and negative n will appear symmetrically about the Dirac point, corresponding to electron and hole states, respectively, as well as an n = 0 state coincident with the Dirac point (1, 7, 8).

The energy progression expected for LLs in graphene can be compared to STS data taken on a nanobubble by fitting the spectra with a sequence of Lorentzian peaks following Eq. 1 with a simple polynomial background (Fig. 2A). The observed peak structure follows the expected progression well, and a value of Bs for position “1” in Fig. 2A was determined to be 350 ± 40 T (energies of the n = 0, 1, and 2 states from the fit are shown below the spectrum for this position). Additional peaks that largely follow the expected sgn(n)|n| progression were also observed over a wider bias range. Small deviations from this progression occur for higher energy peaks, because the graphene dispersion is not strictly linear in this range. A plot of normalized energy (EnEDirac)/ωc versus sgn(n)|n|, shown in Fig. 2C, compiled from spectra on five different nanobubbles, demonstrates the expected scaling behavior for LLs in graphene.

The offset and spacing of the peak progression changes over different positions on individual nanobubbles, indicating a spatial variation of Bs and EDirac. The Bs profile across a strained nanobubble was extracted from the spectral peak spacing at different positions across a nanobubble, shown in Fig. 3B. As seen in Fig. 3A, the Bs profile is reasonably flat across the center of the bubble, indicating a relatively uniform pseudo–magnetic field of 300 to 400 T for this particular geometry. EDirac, coincident with the n = 0 LL, ranges from 0.2 to 0.3 eV across this region of the nanobubble. This variation in EDirac indicates scalar potential variations across the nanobubble, as expected for an elastic deformation that does not constitute a pure shear strain (1, 5).

Fig. 3

(A) Experimental topographic line scan and experimentally determined Bs profile over the tip trajectory shown by black line in (B). (B) STM topography of graphene nanobubble. (C) Topography of theoretically simulated graphene nanobubble with calculated Bs color map. (D) Simulated topographic line scan and Bs profiles extracted from line shown in (C).

To compare the experimental spatial dependence of the strain-induced pseudo–magnetic field to theoretical predictions, we simulated a triangular nanobubble with similar geometry to that shown in Figs. 2A and 3B by using continuum elasticity theory (27) [see supporting information for details (12)]. In this calculation, the edges of a triangular graphene patch were brought in toward the center to simulate the strain arising from the different coefficients of thermal expansion of graphene and the underlying Pt(111) substrate during experimental sample preparation. The predicted pseudo–magnetic field arising from this strain field was also calculated following (4). A 3D plot of the simulated nanobubble shape with the pseudo–magnetic field strength shown as a color map is displayed in Fig. 3C (the corresponding experimental nanobubble is shown in Fig. 3B). The simulated and experimental topographic profiles agree well, as seen by the comparison between Fig. 3, A and D. The simulated Bs profile (Fig. 3D) also agrees reasonably well with the experimentally extracted Bs (Fig. 3A) given that there are some differences in shape and uncertainty in the exact boundary conditions. The main discrepancy is that the experimental Bs profile appears much more uniform than the simulated one, which may be attributed to spatial averaging of the LDOS over the magnetic length scale of ~1.5 nm. The simulated and experimental nanobubbles both exhibit a relatively uniform strain-induced pseudo–magnetic field of 200 to 400 T across the central region with increasing field at the edges.

The exceptional flexibility and strength of graphene membranes (28, 29) coupled with the large strain-induced fields observed suggest that strain engineering of nanoscale energy levels (11, 30) may be a viable means of controlling the electronic structure of graphene, even at room temperature. The experimental demonstration of these enormous pseudo–magnetic fields also provides a new basis for the study of extreme high magnetic field regimes in a condensed-matter environment.

Supporting Online Material

Materials and Methods

SOM Text

Figs. S1 to S7


  • * These authors contributed equally to this work.

  • Present address: Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA.

  • Present address: Department of Physics and Astronomy and Department of Chemistry, University of British Columbia, Vancouver, BC V6T 121, Canada.

References and Notes

  1. Materials and methods are available as supporting material on Science Online.
  2. Research supported by the Director, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, of the U.S. Department of Energy (DOE) under contract no. DE-AC02-05CH11231 (instrumentation development and materials synthesis) and DOE contract no. DE-FG02-08ER46512 (experimental data analysis), by the Office of Naval Research Multidisciplinary University Research Initiative (MURI) award no. N00014-09-1-1066 (experimental STM measurements) and MURI award no. N00014-09-1-1063 (electronic structure calculations), and by Ministerio de Ciencia e Innovación (MICINN) (Spain), award nos. FIS2008-00124/CSD2007-00010 (continuum elasticity simulations). A.H.C.N. acknowledges the support of the Miller Institute for Basic Research. S.A.B. and K.L.M. acknowledge fellowship support from Natural Sciences and Engineering Research Council (Canada) and NSF, respectively. The authors declare no competing financial interests. Requests for materials should be addressed to M.F.C.

Stay Connected to Science

Navigate This Article