Correlated electron states in high mobility two-dimensional electron systems (2DESs), including charge density waves and microemulsion phases intermediate between a Fermi liquid and Wigner crystal, are predicted to exhibit complex local charge order. Existing experimental studies, however, have mainly probed these systems at micron to millimeter scales rather than directly mapping spatial organization. Scanning probes should be well-suited to study the spatial structure of these states, but high mobility 2DESs are found at buried semiconductor interfaces, beyond the reach of conventional scanning tunneling microscopy. Scanning techniques based on electrostatic coupling to the 2DES deliver important insights, but generally with resolution limited by the depth of the 2DES. In this letter, we present our progress in developing a technique called “virtual scanning tunneling microscopy” that allows local tunneling into a high mobility 2DES. Using a specially designed bilayer GaAs/AlGaAs heterostructure where the tunnel coupling between two separate 2DESs is tunable via electrostatic gating, combined with a scanning gate, we show that the local tunneling can be controlled with sub-250 nm resolution.

A local tunnel probe can complement traditional transport measurements for understanding mesoscopic systems.1–10 In many situations, where electron-electron interactions are important, non-trivial spatial ordering of the electrons is a predicted consequence that cannot be fully probed through transport.11 For example, quantum Hall stripes are predicted to form at half-integer filling factors due to a characteristic anisotropy in transport along different crystalline axes in GaAs.12 Anisotropic transport has been measured in this regime,13 likely because of the charge density wave (CDW) associated with these stripes. Scanning tunneling microscopy (STM) has been used extensively to study CDW behavior,14–16 especially in the context of high-temperature superconductivity, but the buried nature of high mobility two-dimensional electron systems (2DESs) precludes this type of measurement.

The work presented in this letter is aimed at developing a technique that can overcome the obstacles to using traditional STM on mesoscopic devices, as depicted in Figure 1. The key idea is to change the nature of the tunnel barrier, from the surface energy barrier, which is fixed by material properties, to a tunable tunnel barrier inside a semiconductor heterostructure that can be modulated by an electric field. Taking advantage of the flexibility of the GaAs/AlGaAs materials system, one can grow heterostructures that support two separate 2DESs,17 and by controlling the aluminum doping define the height of the tunnel barrier between them. The tunneling is controlled locally by scanning a sharp metal tip above the heterostructure, which serves to electrostatically change the height of the tunnel barrier and enhance tunneling preferentially beneath the tip. Therefore, instead of tunneling from a 3D metal into a 2DES as would be the case in STM, the tunneling takes place between two interfacial 2D systems. This technique is called “virtual scanning tunneling microscopy” (VSTM) because a “virtual” STM tip is induced inside the heterostructure between the two 2DESs.

FIG. 1.

Comparison between (a) performing STM on a buried 2DES and (b) performing VSTM on a bilayer 2DES. In (a), the surface energy barrier coupled with the depth of the 2DES results in a small tunnel current from an STM tip. In (b), the presence of the tip locally gates the heterostructure and induces tunneling between the top and bottom 2DESs preferentially beneath the tip. (c) Diagram of the wavefunction extension effect, which is used to modulate tunnel couplings with an applied electric field. Reducing the electron density in one quantum well results in a larger wavefunction overlap between states localized in the two quantum wells, leading to increased tunneling. (d) VSTM sample heterostructure. Two 20 nm wide quantum wells are located 100 nm from each surface, with a 2.1% aluminum barrier between them.

FIG. 1.

Comparison between (a) performing STM on a buried 2DES and (b) performing VSTM on a bilayer 2DES. In (a), the surface energy barrier coupled with the depth of the 2DES results in a small tunnel current from an STM tip. In (b), the presence of the tip locally gates the heterostructure and induces tunneling between the top and bottom 2DESs preferentially beneath the tip. (c) Diagram of the wavefunction extension effect, which is used to modulate tunnel couplings with an applied electric field. Reducing the electron density in one quantum well results in a larger wavefunction overlap between states localized in the two quantum wells, leading to increased tunneling. (d) VSTM sample heterostructure. Two 20 nm wide quantum wells are located 100 nm from each surface, with a 2.1% aluminum barrier between them.

Close modal

A schematic depicting the tunneling process in a VSTM device is shown in Figure 1(c). As the electron density is lowered in one of the 2DESs with an electric field, the subband energy will become less negative relative to the Fermi energy, and therefore the difference between the subband energy and the conduction band edge in the barrier region will be reduced. This implies that as the density is reduced, the tunneling conductance should increase. This effect is called wavefunction extension because the tail of the wavefunction for electrons in one quantum well extends further into the barrier as the density is reduced, thereby leading to a stronger overlap with the wavefunction localized in the other quantum well, and increasing tunneling.

The layer design for the heterostructure investigated in this work is shown in Figure 1(d). Two 2DESs are formed in two 20 nm wide GaAs quantum wells approximately 100 nm from the top and bottom surfaces of the device. The barrier region has 2.1% aluminum doping, bringing the conduction band approximately 1 meV above the Fermi energy in this region. The 2DESs have densities of n = (2,3) × 1011 cm−2 with mobilities of μ = (1,3) × 106 cm2V−1s−1 at 3.3 K, respectively. Large-area tunneling and lateral transport in this heterostructure was reported previously.17,18 The tunneling process was shown to scatter in-plane momentum in these devices, which allows for tunneling between 2DESs with different electron densities. To perform scanning measurements, the active GaAs/AlGaAs layer is thinned to approximately 400 nm using an epoxy bond and stop etch (EBASE) procedure,19 which allows Ti/Au gates to be placed in close proximity to both quantum wells.

A schematic of the device geometry used for VSTM measurements is shown in Figure 2. Depletion gates (DGs) located on the top and bottom of the device are biased such that one ohmic contact connects to the bottom 2DES and one to the top 2DES.20 Therefore, when conductance is measured between the ohmic contacts, current must flow by tunneling in the region between the depletion gates. The “tunnel region” bounded by the depletion gates is gated from the bottom by the tunnel gate. There is no large gate on the top of the device over the tunnel region; instead, this region is gated by a scannable metallized AFM tip. The active mesa region between the top and bottom depletion gates is 5 μm wide and 15 μm long, of which the middle 5 μm is gated by the tunnel gate, as shown in the optical micrograph in Figure 3(a). Not pictured are marks deposited elsewhere on the chip with a 20 μm grid spacing to help locate the device via AFM topography. The measurements are carried out in a home-built scanning gate microscope housed in a cryogen-free dilution refrigerator.21 Although the system is capable of achieving electron temperatures as low as 45 mK, for the discussion that follows all measurements are taken at 3.3 K.

FIG. 2.

Device schematic for VSTM measurements. A metallized AFM tip locally gates the device to enhance the tunnel conductance between the top and bottom 2DES. The white arrows depict the direction of current flow though the device.

FIG. 2.

Device schematic for VSTM measurements. A metallized AFM tip locally gates the device to enhance the tunnel conductance between the top and bottom 2DES. The white arrows depict the direction of current flow though the device.

Close modal
FIG. 3.

Realization of VSTM. Scans over the red bounded region in the optical image measure induced tunneling versus tip position for two tip voltages, −6 V and −10 V. Bright areas in the images indicate strong tunneling induced locally by the AFM tip. The tip height is 300 nm during the measurement, with a 30 nm RMS height oscillation.

FIG. 3.

Realization of VSTM. Scans over the red bounded region in the optical image measure induced tunneling versus tip position for two tip voltages, −6 V and −10 V. Bright areas in the images indicate strong tunneling induced locally by the AFM tip. The tip height is 300 nm during the measurement, with a 30 nm RMS height oscillation.

Close modal

Returning to the schematic of the measurement shown in Figure 2, we would like to determine the tunnel current induced by the tip, but there is also a tunnel current from the areas of the device not gated by the tip. Though this current is smaller per unit area, the area bounded by the depletion gates is 5 μm × 15 μm, much larger than the area gated by the tip, so this background tunneling is significant. To overcome this, we oscillate the height of the tip and measure the tunnel conductance response at the oscillation frequency. This serves to reduce the background large-area gating due to metal on the cantilever. Oscillating the height of the tip rather than the voltage on the tip leads to crisper imaging results. If we oscillate the tip voltage, metal on the cantilever capacitively modulates tunneling over a large area, whereas if we oscillate the tip height the tunneling is strongly modulated only where the tip is closest to the device.

We can now start analyzing the response of the device to local gating. In Figure 3, the induced tunneling over the red bounded region in the optical micrograph is measured under two different gating conditions, −6 V on the tip and −10 V on the tip in Figures 3(b) and 3(c), respectively. In both cases, the tip height is 300 nm above the surface with a 30 nm RMS height oscillation, the tunnel gate is set at −490 mV, and an interlayer bias of 15 mV is applied to the device. For this device geometry, a signal is measured only when the tip has induced tunneling locally in the device. The response of the device changes significantly under the different gating conditions. In both scans, there is a region of strongest induced tunneling, where ΔI > 40 pA. This region represents the gating condition where the device is being operated at maximum gain. This region shifts upward in the image from −6 V to −10 V due to the effect of the bottom depletion gate, which is underneath the device over the bottom half of the images.

The data for Vtip = −10 V demonstrate the device response when it is “over-gated,” meaning, the top 2DES has been depleted past the mobility edge. Comparing the width of the mesa as seen in both images, it is narrower in the −10 V image and there is even a region near the edges where the response is negative, implying that the top 2DES is being totally depleted at the edges. In addition, an emergent fine structure is visible in the lower half of the image, a consequence of the non-uniform potential in the top 2DES in the presence of disorder. However, at this large tip voltage, the density in the bottom 2DES can also be affected by the tip, since the top 2DES is not acting to screen the tip potential from the bottom layer.

One important parameter that can be extracted from the data in Figure 3 is the lateral resolution of the measurement. We can estimate the lateral resolution by measuring how sharp the edges of the mesa appear in the tunneling data. This is accomplished by taking a line cut perpendicular to the mesa and taking a derivative to measure the width of the peak at the edge of the mesa. This analysis is shown in Figure 4 for two different tip heights, 300 nm and 150 nm above the surface.

FIG. 4.

The lateral resolution of the measurement can be estimated by measuring how sharp the edge of the device appears in the data at a tip height of (a) 300 nm (Vtip = −6 V) and (b) 150 nm ( V tip = 4 V) above the surface. The width of the peak in the derivative along a line cut gives a measure of the sharpness of the mesa, which improves from 570 nm to 245 nm as the tip is moved closer to the device.

FIG. 4.

The lateral resolution of the measurement can be estimated by measuring how sharp the edge of the device appears in the data at a tip height of (a) 300 nm (Vtip = −6 V) and (b) 150 nm ( V tip = 4 V) above the surface. The width of the peak in the derivative along a line cut gives a measure of the sharpness of the mesa, which improves from 570 nm to 245 nm as the tip is moved closer to the device.

Close modal

As expected, the resolution improves as the tip is brought closer to the device since the electric field from the tip becomes sharper as the tip is moved closer. The average width improves from 570 nm to 245 nm as the tip is brought closer, where the width of the peak is measured at the point where the signal falls to e 1 of its peak value. From simulation results published previously,17 the most significant factors that determine spatial resolution are the radius of the tip and the distance between the tip and the top 2DES. Although the tip bias and tip vibration amplitude can change the magnitude of the tunneling signal, they are not the main factors dictating the spatial resolution in the limit of a weak perturbation from the tip. The trace in Figure 4(b) indicates that there are new features in the data (around a distance of 4 μm), not seen in Figure 4(a). These features are due to a combination of surface contamination and charge traps in the dielectric on the surface, challenges which must be addressed in the device fabrication before pushing the measurements to a higher resolution. A higher resolution should be achievable, since electrostatic simulations17 predict a spatial resolution of 40 nm. The spatial resolution required for the technique depends on the specific system to be studied, but one can use the Fermi wavelength at a given electron density as a suitable metric. At an electron density of 5 × 10 10 cm 2, which can be achieved in current devices before hitting the mobility edge, the Fermi wavelength is λ F = 2 π / n = 112 nm. To observe Friedel oscillations associated with electron scattering off defects and edges in the 2DES, as have been studied with STM for surface 2DESs,22 one would require a spatial resolution smaller than half the Fermi wavelength, or 56 nm in this example.

The discussion presented in this letter demonstrates the viability of VSTM as a powerful scanning probe technique to study mesoscopic electron systems. We have shown that the tunneling between two buried 2DESs can be modulated locally with a scanning gate on a length scale below 250 nm, but need remains to optimize the device design and fabrication procedure to achieve the expected spatial resolution of 40 nm. The future directions for VSTM are very broad, and in the near term can include measurements of local tunneling into quantum Hall edge states, along with observing the spatially oscillatory electron density of states near impurities and edges. Also, by varying the interlayer bias applied to a VSTM device, the energy dependent density of states can be measured locally, in analogy to STM. More generally, the VSTM technique can be extended to any interfacial materials system where the tunnel barrier can be engineered through growth.

This work was supported by the Department of Energy, Office of Basic Energy Sciences, under Contract No. DE-AC02-76SF00515. The original VSTM concept was developed under the Center for Probing the Nanoscale (NSF NSEC Grant Nos. PHY-0425897 and PHY-0830228) and a Mel Schwartz Fellowship from the Stanford Physics Department. The work at Princeton was partially funded by the Gordon and Betty Moore Foundation through Grant No. GBMF2719, and by the National Science Foundation MRSEC-DMR-0819860 at the Princeton Center for Complex Materials. M.P. acknowledges support from the Hertz Foundation, NSF, Stanford, and Intel.

1.
M. A.
Eriksson
,
R. G.
Beck
,
M.
Topinka
,
J. A.
Katine
,
R. M.
Westervelt
,
K. L.
Campman
, and
A. C.
Gossard
,
Appl. Phys. Lett.
69
,
671
(
1996
).
2.
J.
Martin
,
N.
Akerman
,
G.
Ulbricht
,
T.
Lohmann
,
J. H.
Smet
,
K.
von Klitzing
, and
A.
Yacoby
,
Nat. Phys.
4
,
144
(
2008
).
3.
S. H.
Tessmer
,
P. I.
Glicofridis
,
R. C.
Ashoori
,
L. S.
Levitov
, and
M. R.
Melloch
,
Nature
392
,
51
(
1998
).
4.
M. A.
Topinka
,
B. J.
LeRoy
,
R. M.
Westervelt
,
S. E.
Shaw
,
R.
Fleischmann
,
E. J.
Heller
,
K. D.
Maranowski
, and
A. C.
Gossard
,
Nature
410
,
183
(
2001
).
5.
M. P.
Jura
,
M. A.
Topinka
,
L.
Urban
,
A.
Yazdani
,
H.
Shtrikman
,
L. N.
Pfeiffer
,
K. W.
West
, and
D.
Goldhaber-Gordon
,
Nat. Phys.
3
,
841
(
2007
).
6.
B.
Hackens
,
F.
Martins
,
S.
Faniel
,
C. A.
Dutu
,
H.
Sellier
,
S.
Huant
,
M.
Pala
,
L.
Desplanque
,
X.
Wallart
, and
V.
Bayot
,
Nat. Commun.
1
,
39
(
2010
).
7.
D. V.
Pelekhov
,
J. B.
Becker
, and
G.
Nunes
,
Appl. Phys. Lett.
72
,
993
(
1998
).
8.
R.
Crook
,
C.
Smith
,
W.
Tribe
,
S.
O'Shea
,
M.
Simmons
, and
D.
Ritchie
,
Phys. Rev. B
66
,
121301(R)
(
2002
).
9.
J. A.
Hedberg
,
A.
Lal
,
Y.
Miyahara
,
P.
Grutter
,
G.
Gervais
,
M.
Hilke
,
L.
Pfeiffer
, and
K. W.
West
,
Appl. Phys. Lett.
97
,
143107
(
2010
).
10.
A.
Vinante
,
G.
Wijts
,
O.
Usenko
,
L.
Schinkelshoek
, and
T. H.
Oosterkamp
,
Nat. Commun.
2
,
572
(
2011
).
11.
H.
Pothier
,
S.
Guéron
,
N.
Birge
,
D.
Esteve
, and
M.
Devoret
,
Phys. Rev. Lett.
79
,
3490
(
1997
).
12.
M. P.
Lilly
,
K. B.
Cooper
,
J. P.
Eisenstein
,
L. N.
Pfeiffer
, and
K. W.
West
,
Phys. Rev. Lett.
82
,
394
(
1999
).
13.
E.
Fradkin
and
S. A.
Kivelson
,
Phys. Rev. B
59
,
8065
(
1999
).
14.
W. D.
Wise
,
M. C.
Boyer
,
K.
Chatterjee
,
T.
Kondo
,
T.
Takeuchi
,
H.
Ikuta
,
Y.
Wang
, and
E. W.
Hudson
,
Nat. Phys.
4
,
696
(
2008
).
15.
C.
Brun
,
Z.-Z.
Wang
,
P.
Monceau
, and
S.
Brazovskii
,
Phys. Rev. Lett.
104
,
256403
(
2010
).
16.
R.
Coleman
,
B.
Giambattista
,
P.
Hansma
,
A.
Johnson
,
W.
McNairy
, and
C.
Slough
,
Adv. Phys.
37
,
559
(
1988
).
17.
A.
Sciambi
,
M.
Pelliccione
,
S. R.
Bank
,
A. C.
Gossard
, and
D.
Goldhaber-Gordon
,
Appl. Phys. Lett.
97
,
132103
(
2010
).
18.
A.
Sciambi
,
M.
Pelliccione
,
M.
Lilly
,
S.
Bank
,
A.
Gossard
,
L.
Pfeiffer
,
K.
West
, and
D.
Goldhaber-Gordon
,
Phys. Rev. B
84
,
085301
(
2011
).
19.
M.
Weckwerth
,
J.
Simmons
,
N.
Harff
,
M.
Sherwin
,
M.
Blount
,
W.
Baca
, and
H.
Chui
,
Superlattices Microstruct.
20
,
561
(
1996
).
20.
J. P.
Eisenstein
,
L. N.
Pfeiffer
, and
K. W.
West
,
Appl. Phys. Lett.
57
,
2324
(
1990
).
21.
M.
Pelliccione
,
A.
Sciambi
,
J.
Bartel
,
A. J.
Keller
, and
D.
Goldhaber-Gordon
,
Rev. Sci. Instrum.
84
,
033703
(
2013
).
22.
M. F.
Crommie
,
C. P.
Lutz
, and
D. M.
Eigler
,
Nature
363
,
524
(
1993
).