Transport evidence for a sliding two-dimensional quantum electron solid.

Ignited by the discovery of the metal-insulator transition, the behaviour of low-disorder two-dimensional (2D) electron systems is currently the focus of a great deal of attention. In the strongly interacting limit, electrons are expected to crystallize into a quantum Wigner crystal, but no definitive evidence for this effect has been obtained despite much experimental effort over the years. Here, studying the insulating state of a 2D electron system in silicon, we have found two-threshold voltage-current characteristics with a dramatic increase in noise between the two threshold voltages. This behaviour cannot be described within existing traditional models. On the other hand, it is strikingly similar to that observed for the collective depinning of the vortex lattice in type-II superconductors. Adapting the model used for vortexes to the case of an electron solid yields good agreement with our experimental results, favouring the quantum electron solid as the origin of the low-density state.

Introduction

Recently, there has been progress in understanding the metallic side of the metal-insulator transition in low-disorder strongly-interacting 2D electron systems[1-10], while the origin of the low-density state on the insulating side has remained a mystery. Experimental investigations of the transport and thermodynamic properties of 2D electrons in semiconductors have suggested that these systems approach a phase transition at low electron densities to a new, unknown state that could be a quantum Wigner crystal or a precursor[11-16]; the term quantum means that the corresponding kinetic energy of 2D electrons is determined by the Fermi energy in contrast to the well-known classical Wigner crystal[17] in which the kinetic energy of electrons is determined by temperature. The phase transition point in the least-disordered 2D electron systems in semiconductors was found to be close to the critical electron density for the metal-insulator transition, below which the 2D electrons become localized (immobile). Although the insulating side of the metal-insulator transition has been extensively studied[18-20], no definitive conclusion has been reached concerning the origin of the low-density state. Nonlinear current-voltage (I − V) curves were observed and interpreted as either manifestation of the depinning of an electron solid[18,21,22] or the breakdown of the insulating phase within traditional scenarios such as strong electric field Efros-Shklovskii variable range hopping[23] or percolation (see, e.g., refs.[3,24] and references therein). The observation of broad-band voltage noise at the threshold I − V curves as well as the attempt to probe the low-density state in perpendicular magnetic fields also have not provided information that allows a choice between the depinning of the electron solid or traditional mechanisms[3,24]. It is worth noting that much confusion was introduced by the fact that many authors chose to interpret their data in terms of Wigner crystal, ignoring mundane interpretations.

Here we report a significant breakthrough in our understanding of the origin of the low-density state in a strongly interacting 2D electron system in silicon metal-oxide-semiconductor field-effect transistors (MOSFETs). We have observed two-threshold V − I characteristics with a dramatic increase in noise between the two threshold voltages at the breakdown of the insulating state. In the form of fluctuations with time, the noise in current increases dramatically above Vth1 and essentially disappears above Vth2. It is the sharp noise peak on the V − I curves that makes the two-threshold behaviour evident. The double threshold behaviour is very similar to that observed for the collective depinning of the vortex lattice in type-II superconductors (see, e.g., ref.[25]) provided the voltage and current axes are interchanged. This strongly favours the sliding 2D quantum electron solid, whereas the double threshold behaviour cannot be described within alternative scenarios such as percolation or overheating. We emphasize that rather than being an ideal Wigner crystal, the 2D electron system under study is likely to be closer to an amorphous solid, which is similar to the case of the vortex lattice in type-II superconductors, where the collective pinning was observed. It is important to note that in the earlier studies, where current was passed through the sample and the voltage between potential probes was measured, no distinct features were found on almost flat V(I) curves in the breakdown regime. The main problem with those studies is that the noise in the voltage signal is small, which prevents one from identifying a second threshold voltage. By contrast, we applied dc voltage between source and drain and measured the induced current. We found that in this measurement configuration, the noise in current is clear and dramatic. It is this experimental protocol that enabled the discovery of features that had not been observed before.

Results

Voltage-current characteristics

Figure 1 shows a set of low-temperature voltage-current curves at different electron densities in the insulating regime ns < nc (here nc ≈ 8 × 1010 cm−2 is the critical density for the metal-insulator transition in this electron system); the corresponding interaction parameter given by the ratio of the Coulomb and Fermi energies, rs = gv/(πns)1/2aB, exceeds rs ~ 20 (where gv = 2 is the valley degeneracy, ns is the areal density of electrons and aB is the effective Bohr radius in semiconductor). Two threshold voltages are observed at electron densities below ≈6 × 1010 cm−2: with increasing applied voltage, the current is near zero up to a voltage threshold Vth1, then increases sharply until a second threshold voltage Vth2 is reached, above which the slope of the V − I curve is significantly reduced, and the behaviour is linear although not ohmic (see also the top inset to Fig. 1). As the electron density is increased, the value of Vth1 decreases while the second threshold becomes less pronounced and eventually disappears. No hysteresis was observed for the range of electron densities studied. We point out that the observed behaviour (see also Fig. 2) is quite distinct from that reported in the insulating state in amorphous InO films, where the current was found to jump at the threshold voltage by as much as five orders of magnitude, and the V − I curves exhibit hysteresis consistent with bistability and electron overheating[26,27]. It is important to note also that in our experiment the power dissipated near the onset Vth1 is less than 10−16 W, which is unlikely to cause significant electron overheating; for comparison, the power dissipated near the threshold voltage in ref.[26] is more than three orders of magnitude higher. Note also that the occurrence of the double threshold cannot be explained within the percolation picture in which case a single threshold is expected[3]. Thus, the double threshold behaviour cannot be described within existing traditional models.

Fig. 1

Voltage-current characteristics. V − I curves are shown for different electron densities in the insulating state at a temperature of 60 mK. The dashed lines are fits to the data using Eq. (5). The top inset shows the V − I curve for ns = 5.20 × 1010 cm−2 on an expanded scale; also shown are the threshold voltages Vth1 and Vth2, the static threshold Vs = Vth2, and the dynamic threshold Vd that is obtained by the extrapolation of the linear region of the V − I curve to zero current. Bottom inset: activation energy Uc vs. electron density. Vertical error bars represent standard deviations in the determination of Uc from the fits to the data using Eq. (5). The dashed line is a linear fit

Fig. 2

The current as a function of time. Current is plotted as a function of time for ns = 5.2 × 1010 cm−2 and T = 60 mK at voltages V = 4.90 mV (which lies between Vth1 and Vth2; lower curve) and V = 5.44 mV (above Vth2). Inset: colour map of the broad-band noise at ns = 5.36 × 1010 cm−2 on a (V, T) plane

In the regime where both thresholds are present, the current measured at voltages between Vth1 and Vth2 exhibits strong fluctuations with time that are comparable to its value. Above the second threshold, Vth2, these anomalously large fluctuations disappear and the noise is barely perceptible. It is the sharp increase in noise on the V − I curves that makes the two-threshold behaviour evident. The noise is shown explicitly in Fig. 2, where the current is plotted as a function of time for density ns = 5.2 × 1010 cm−2.

Figure 3a shows the V − I characteristics for ns = 5.36 × 1010 cm−2 at different temperatures. As the temperature, T, is increased, the second threshold Vth2 becomes less pronounced and the threshold behaviour of the V − I curves eventually smears out due to the shrinkage of the zero-current interval.

Fig. 3

Voltage-current characteristics and noise. a V − I characteristics at ns = 5.36 × 1010 cm−2 for different temperatures. The dashed lines are fits to the data using Eq. (5). b The broad-band noise as a function of voltage for the same electron density and temperatures. The three upper curves are shifted vertically for clarity

Noise

The measured broad-band noise is shown as a function of voltage in Fig. 3b for different temperatures at electron density ns = 5.36 × 1010 cm−2. The inset to Fig. 2 is a colour map of the broad-band noise on a (V, T) plane. A large increase in the noise is observed between the thresholds Vth1 and Vth2 at the lowest temperature. This large noise decreases rapidly with increasing temperature in agreement with the two-threshold behaviour of the V − I curves of Fig. 3a.

The spectrum of the generated noise, measured at its largest value, is displayed in Fig. 4. The noise increases with decreasing frequency, f, according to the 1/f law with α = 0.6 ± 0.1, which is close to unity. This finding is consistent with the fact that the noise of the form 1/f with α close to unity was found in the linear regime of response in similar samples near the metal-insulator transition[28,29]. This indicates a universal behaviour of the noise spectrum.

Fig. 4

Noise spectrum. Noise is plotted as a function of frequency at ns = 5.36 × 1010 cm−2, T = 60 mK and V = 4.26 mV with resolution/bandwidth 2.5 Hz. The broad maxima at f ~ 10 and 60 Hz on the order of 1 pA Hz−1/2 are within the experimental uncertainty. The dashed line shows the 1/f0.6 dependence

Discussion

We will now analyze our results in light of a phenomenological theory based on pinned elastic structures. There is a striking similarity between the double-threshold V − I dependences in the low-density state of Si MOSFETs and those (with the voltage and current axes interchanged) known for the collective depinning of the vortex lattice in type-II superconductors (see, e.g., ref.[25]). The physics of the vortex lattice in type-II superconductors, in which the existence of two thresholds is well known, can be adapted for the case of an electron solid. Current flows for zero voltage in a superconductor; the depinning of the vortex lattice occurs when a non-zero voltage appears. In our case, the situation is reciprocal: a voltage is applied but at first the current is zero in the limit of zero temperature; the depinning of the electron solid is signalled by the appearance of a non-zero current. The transient region between the dynamic (Vd) and static (Vs) thresholds corresponds to the collective pinning of the solid. In this region the pinning occurs at the centres with different energies and the current is thermally activated:where U(V) is the activation energy. The static threshold Vs = Vth2 signals the onset of the regime of solid motion with friction. This corresponds to the conditionwhere E is the electric field and L is the characteristic distance between the pinning centres with maximal activation energy Uc. From the balance of the electric, pinning, and friction forces in the regime of solid motion with friction, one expects a linear V − I characteristic that is offset by the threshold Vd corresponding to the pinning forcewhere σ0 is a coefficient; note that near Vd, one can in general expect a power-law behaviour of (V − Vd). Assuming that the activation energy for the electron solid is equal towe obtain the expression for the currentThe fits to the data using Eq. (5) are shown by dashed lines in Figs. 1 and 3a. As seen from the figures, the experimental two-threshold V − I characteristics are described well by Eq. (5). The value of Uc decreases approximately linearly with electron density and tends to zero at ns ≈ 6 × 1010 cm−2 (the bottom inset to Fig. 1). This is in contrast to the vanishing activation energy of electron-hole pairs at nc obtained by measurements of the resistance in the limit of zero voltages/currents[30]. Presumably, the vanishing Uc is related to the minimum number of the strong pinning centres for which the collective pinning is still possible. The fact that the coefficient σ0 is approximately constant (σ0 ≈ 1.6 × 10−7 Ohm−1) indicates that the solid motion with friction is controlled by weak pinning centres[25]. We argue that the large noise in the regime of the collective pinning of the solid between Vd and Vs should be suppressed in the regime of solid motion with friction at V > Vs. Indeed, in the regime of the collective pinning of the solid between Vd and Vs, the solid deforms locally when the depinning occurs at some centre and then this process repeats at another centre etc., which leads to the generation of a large noise. In contrast, in the regime of solid motion with friction above the second threshold Vs, the solid slides as a whole due to the overbarrier motion, resulting in the suppression of noise. Thus, the physics of pinned periodic/elastic objects is relevant for the low-density state in a 2D electron system in silicon MOSFETs.

One can estimate the characteristic range of frequencies of the generated noise f ~ vd/L, where vd is the drift velocity. Using the parameters corresponding to the data shown in Fig. 3a (Vs ≃ 4.46 mV, Uc/kB ≃ 1.7 K, I ~ 20 pA, ns = 5.36 × 1010 cm−2, and L ~ 10−3 cm as determined from Eq. (2)), one obtains f ~ 500 Hz, which is in reasonable agreement with the experiment. Note that in this experiment we did not observe the narrow band noise that is expected at frequencies close to the washboard frequency related to the motion of the solid (either the noise signal might be too small or the frequency of the noise might be too high). Although the model used describes the experiment successfully, further in-depth theoretical consideration is needed.

Methods

Samples

Measurements were made in an Oxford dilution refrigerator with a base temperature of ≈50 mK on (100)-silicon MOSFETs with a peak electron mobility close to 3 m2 V−1 s−1 at T < 0.1 K similar to those described in detail in ref.[31]. The electron density was controlled by applying a positive dc voltage to the gate relative to the contacts; the oxide thickness was 150 nm. Samples had a Hall bar geometry of width 50 μm. To overcome the main experimental obstacle in the low-density low-temperature limit—high contact resistance—thin gaps were introduced in the gate metallization that allow a high electron density to be maintained near the contacts regardless of the density in the main part of the sample. As a result, contact resistances did not exceed ~10 kOhm and could be disregarded in the insulating state.

V − I characteristics and noise

Measurements of the V − I characteristics and noise were obtained in the main part of the sample with length 180 μm. The dc current and noise were measured using an ultra-low-noise current-voltage converter FEMTO DHPCA-100 (with 10 kOhm input resistance) connected to a digital voltmeter HP34401A Agilent Multimeter or lock-in SR830.