## Electron Breakdown

An electron possesses charge and spin. In general, these properties are confined to the electron. However, in strongly interacting many-body electronic systems, such as one-dimensional wires, it has long been theorized that the charge and spin should separate. There have been tantalizing glimpses of this separation experimentally, but questions remain. **Jompol et al.** (p. 597) looked at the tunneling current between an array of one-dimensional wires and a two-dimensional electron gas and argue that the results reveal a clear signature of spin-charge separation.

## Abstract

In a one-dimensional (1D) system of interacting electrons, excitations of spin and charge travel at different speeds, according to the theory of a Tomonaga-Luttinger liquid (TLL) at low energies. However, the clear observation of this spin-charge separation is an ongoing challenge experimentally. We have fabricated an electrostatically gated 1D system in which we observe spin-charge separation and also the predicted power-law suppression of tunneling into the 1D system. The spin-charge separation persists even beyond the low-energy regime where the TLL approximation should hold. TLL effects should therefore also be important in similar, but shorter, electrostatically gated wires, where interaction effects are being studied extensively worldwide.

The effects of interactions are almost impossible to calculate in a general many-particle system, although they cannot be ignored. However, for a one-dimensional (1D) system, Luttinger, building on an approximation scheme of Tomonaga, constructed a soluble 1D model with infinite linear dispersion and a restricted set of interactions. The solution has the remarkable property that the excitations of spin and charge behave independently and move with different speeds. It has been argued (*1*) that all 1D metals are adiabatically continuous with the Tomonaga-Luttinger model at low energies, and hence spin-charge separation should be observable in real systems. Determining the extent of its applicability would provide a major test of more general methods of modeling interaction effects, with relevance to quantum devices and the theory of high-temperature superconductivity. Recent work (*2*) presents a more general theory of 1D systems with a nonlinear dispersion, but the effects of spin are not yet included.

Some properties of the Tomonaga-Luttinger liquid (TLL), such as power-law behavior, have been observed and studied in a variety of systems, such as carbon nanotubes (*3*) and edge states in the fractional quantum Hall regime (*4*), but these experiments have not directly resolved the dispersion of the excitations in a TLL. Only a few experiments have attempted to detect the spin-charge separation directly, for example, by photoemission spectroscopy (*5*–*7*) and tunneling spectroscopy between a pair of closely spaced cleaved-edge-overgrown quantum wires (*8*–*10*). The latter 1D-1D tunneling results are striking and provide some evidence of dispersing spinons and holons—the excitations of a TLL. However, in these experiments TLLs act as both probe and subject, so an independent study—in a different geometry—of the excitation spectrum is vital to be sure of the interpretation.

We use a 2D electron gas (2DEG) as the (well-understood) probe layer. Use of an array of highly regular wires averages out impurity, length-resonance, and charging effects, and measurements of power-law behavior and spin-charge separation can be made with just the lowest 1D subband, without electrons becoming localized. Interpretation of TLL results is thus much easier, and results obtained in the nonlinear regime, where the theory is much less well developed, can be directly interpreted as a modification of the 1D elementary excitations. Our 1D wires are formed using split gates, where the confinement is much weaker than in the overgrown wires, making the results relevant to a broad range of other devices. The tunneling current *I* between the 1D wires and an adjacent low-disorder layer containing a 2DEG depends on the overlap between the spectral functions of the two systems. This overlap is varied by using an in-plane magnetic field *B* perpendicular to the wires to offset the two spectral functions in *k* space by Δ*k* = *eBd*/*ħ* along the wires, where *d* is the center-to-center tunneling distance between the two systems (*11*). By applying a positive bias *V*_{dc} to the 2DEG relative to the wires, electrons tunnel into excited states of the 2DEG, from matching states below the Fermi energy in the 1D wires, allowing investigation of the energy dependence. Thus, the 2D system acts as a spectrometer, and the differential conductance *G* displays resonant-tunneling peaks corresponding to overlapping points in the offset dispersion relations, where energy and momentum are conserved.

Our devices contain an array of about 350 extremely regular quantum wires in the upper layer of a GaAs-AlGaAs double-quantum-well (DQW) structure (Fig. 1A). The two layers are separately contacted using a surface-gate depletion scheme (*12*–*14*). The wires, of length *L* = 17.5 μm and lithographic width 0.17 μm (device A) and 0.18 μm (device B), were squeezed by a negative gate voltage *V*_{wg}. There is an additional small ungated region p of width 0.9 μm, which provides a current path to the entrances to the 1D wires. The tunneling conductance *G* = d*I*/d*V*_{dc} between the 1D wires and the 2D layer was measured as a function of source-drain bias *V*_{dc} in an in-plane magnetic field *B* perpendicular to the wires at lattice temperatures *T* down to ~40 mK.

First, the effect of the wire-gate voltage *V*_{wg} was characterized (Fig. 1). Initially, at small *V*_{wg}, the upper 2DEG screens the lower 2DEG. The first drop in *G* occurs at −0.40 V when the top well depopulates, leaving narrow 1D wires between the gates. Then, the bottom well is gradually depleted until it pinches off at −0.80 V. *B* was then swept at each *V*_{wg} (with *V*_{dc} = −0.3 mV, to avoid the zero-bias anomaly described later), revealing a series of peaks (Fig. 1B), one pair of peaks at *B _{i}*

^{±}for the

*i*th 1D subband (

*i*≥ 1). For this device (A, 0.17 μm wire width), six peaks are observed at

*V*

_{wg}= −0.49 V, corresponding to three occupied 1D subbands (

*13*). With more negative

*V*

_{wg}, the 1D subband spacing increases and the 1D density decreases. Hence, the number of subbands can be reduced to just one. At

*V*

_{wg}~ −0.65 V, the wires become insulating. The widths of the wire gates and of the long narrow 2D p region were chosen carefully such that even with just a single 1D subband in the top 2DEG there was minimal modulation of the lower 2DEG and current could still reach the wires. However, the p region inevitably contributes a 2D-2D parasitic tunneling current that appears in the measurements, but this is small and independent of the tunnel current from the wires, even after wire pinch-off, so it can be measured and allowed for.

We now choose *V*_{wg} = −0.62 V, well into the region where there is just one 1D subband. Figure 2, A and B, shows *G* as a function of the dc interlayer bias (*V*_{dc}) and *B* as color-scale plots, at high (1 K) and low (~40 mK) lattice temperatures, respectively. In the absence of interactions, there should be peaks in *G* that track the 1D and 2D parabolic dispersion relations; these are indicated with dashed and solid black curves, respectively. The parasitic 2D-2D tunneling in the ungated p region is also shown, as green dash-dotted lines labeled p. The crossing points along *V*_{dc} = 0 occur at *B*^{−} = 1.08 T and *B*^{+} = 3.62 T (Fig. 2B). The 1D Fermi wave vector is *k*_{F1D} = *ed*(*B*^{+} − *B*^{−})/2*ħ*, giving the approximate electron density in the wires *n*_{1D} ≅ 40 μm^{−1} (from *k*_{F1D} = π*n*_{1D}/2). The upturned (dashed) parabola maps the energy of the 1D wires as a function of wave vector *k*.

The curves drawn in Fig. 2 are those expected for single-particle tunneling. However, there is an additional region of high conductance to the right of *B*^{+}. It drops off along a straight line moving diagonally up and to the left from *B*^{+} (indicated with an arrow). To show this more clearly, the data are differentiated with respect to *B* and plotted in Fig. 2C. We have taken detailed data in this region, sweeping *B*, in another device (B) in the same single-subband regime (Fig. 3B). We find a region of large negative d*G*/d*B* along a straight line to the right of the high-conductance region; this appears dark blue and is indicated by the red dashed line. We compare this with the theoretical prediction for noninteracting electrons, shown in Fig. 3A, where a similar dark blue feature also occurs. However, the latter tracks the 1D parabola along its length. In contrast, the feature in Fig. 3B disperses away from the 1D parabola. Because features in the conductance reflect singular features in the spectral function, we can conclude that the 1D parabola and the red dashed line track the dispersion of two independent excitations, which in the TLL framework correspond to the spinon and holon, for spin and charge excitations, respectively.

To confirm this interpretation, we have calculated the tunneling spectra for a TLL coupled to a 2D system of electrons. The framework for these calculations already exists in the literature (*15*, *16*), so we only describe the relevant details here. To compute the tunneling current, we require the spectral function of a TLL, which in general depends on four parameters: the spinon and holon velocities *v*_{σ} and *v*_{ρ}, respectively, plus two exponents γ_{σ} and γ_{ρ}. We assume repulsive, spin-rotation-invariant interactions, which implies γ_{σ} = 0 and *v*_{ρ} > *v*_{σ}; in the absence of back-scattering, *v*_{σ} ≈ *v*_{F1D}, the 1D Fermi velocity. For simplicity, we also neglect interbranch scattering processes, which sets γ_{ρ} = 0. The spectral function of the 2D system is taken to be a Lorentzian of width Γ, where Γ is the disorder-scattering rate in the 2D system (*17*).

This minimal model is convenient for analytical calculations and is expected to reproduce the main features of the tunneling spectra associated with spin-charge separation. It will not, however, capture the zero-bias anomaly, which is absent for γ_{ρ} = 0. Finite γ_{ρ} is also expected (*16*) to lead to a weakening of the singular features in the spectra, making our calculation the “optimum” case for observing spin-charge separation. Because the TLL model relies on the linearization of the spectrum at low energies, the model’s results are only formally applicable in the low-bias part of the spectrum, where this linear approximation holds.

Figure 3C shows the calculated conductance in the vicinity of the *B*^{+} point as a function of scaled magnetic field at fixed voltage for noninteracting electrons, a clean TLL, and a TLL with a little disorder broadening. The singular peak at high magnetic field is independent of the form of the excitations, and hence is common to both systems. The peak at lower magnetic field for the noninteracting system splits, in the TLL case, into a weakened singular peak (which for *v*_{σ} = *v*_{F1D} occurs at the same magnetic field as the noninteracting peak), plus a finite discontinuity away from the peak. This sudden drop in the conductance away from the peak is precisely what is observed experimentally in Fig. 2, A and B. Figure 3D shows d*G*/d*B* with a little disorder. The spinon is identified as a maximum/minimum pair, which disperse together in the *B*-*V*_{dc} plane, as shown in Fig. 3E, for more realistic disorder broadening. The holon is identified as a single minimum, which disperses away from the spinon. We note that an extra feature is also predicted to occur at the *B*^{−} point. However, because all three features from this point remain in proximity to each other, it is difficult to resolve them individually.

The experimental results in Fig. 3B are consistent with the predictions in Fig. 3E, at least in the low-bias regime where the linear approximation is reasonable. We have observed this extra feature in the conductance in three devices, at temperatures between ~60 mK and 1.8 K, and on several thermal cycles. We thus conclude that we are observing spin-charge separation, the hallmark of a TLL.

An additional feature that cannot be explained in the noninteracting model is the zero-bias anomaly (ZBA), the strong suppression of conductance along *V*_{dc} = 0, visible as a dark-blue line in Fig. 2, A and B. This is likely to be related to the energy cost for an electron to tunnel into or out of a 1D wire, as it disturbs the line of electrons on either side of it. It has previously been observed in 1D-1D tunneling (*8*), and, for a TLL, *G* is expected to have identical power-law dependences on *V*_{dc} and temperature *T*.

Figure 4A shows the ZBA in the tunneling conductance *G* as a function of *V*_{dc} at various *T*, at a field midway between *B*^{−} and *B*^{+}. Figure 4, B and C, shows *G*(*V*_{dc} = 0,*T*) and *G*(|*V*_{dc}|, *T* < 70 mK), respectively, on log-log plots. Both clearly vary as power laws, as labeled, over considerably more than one order of magnitude. The corresponding power-law exponents α_{T} ≈ 0.45 ± 0.04 and α_{V} ≈ 0.52 ± 0.04 are very similar, as expected. To illustrate that temperature and bias play a similar role in smearing the energy, *G*(|*V*_{dc}|) is plotted in Fig. 4D as a function of *V*′_{dc} = , the simplest way of adding forms of energy smearing together as noise in quadrature. (3*k*_{B}*T* is an estimate of the thermal energy spread, and the factor of 3 is chosen to match the offset between the curves in Fig. 4, B and C.)

The TLL is characterized by the “anomalous” exponent γ_{ρ} and the spinon and holon velocities. γ_{ρ} can be expressed in terms of an interaction parameter *g* [γ_{ρ} = (*g* + 1/*g* − 2)/8], which indicates the strength of the interactions (*g* = 1 for noninteracting particles). Accurate a priori calculation of *g* cannot be done in general. However, an estimate based on the 1D densities used suggests *g* ~ 0.7 (supporting online text). We now determine *g* from the various experimental results.

In TLL theory, the spin velocity *v*_{σ} is approximately equal to the Fermi velocity *v*_{F1D} for weak backscattering (*v*_{σ} < *v*_{F1D} for stronger interactions) (*18*, *19*). Thus, the spin mode should follow the low-bias part of the noninteracting 1D parabola (or even lie below it). The charge mode should propagate faster than the spin mode for repulsive interactions. Thus, we label the curve that follows the 1D parabola “spin” and the extra line “charge.” For device B, we deduce the velocities close to zero bias to be *v*_{σ} ≈ *v*_{F1D} = *ħk*_{F}/*m*_{1D} = 1.13 × 10^{5} ms^{−1} and *v*_{ρ} ≈ 1.4 *v*_{F1D}, respectively. In a Galilean invariant system, *g* = *v*_{F1D}/*v*_{ρ} (*19*), so this gives *g* ≈ 0.7, which is consistent with the estimate given above. For device A, *g* is very similar to that of device B, as is the 1D electron density. Another device, with a thinner (7 nm) barrier, yields *g* ≈ 0.65.

Within TLL theory, the power-law exponents extracted from the ZBA can also be directly related to *g*. The form of the exponent depends on whether the excitations are affected by the ends of the wire. This is determined by comparing the energy of the tunneling electron to Δ*E* = 2*ħv*_{F1D}/*gL*, which is related to the inverse time scale for the holon to travel to the end of a wire of length *L* (*8*, *20*). For *k*_{B}*T*, *eV*_{dc} >> Δ*E*, the ends are unimportant and the process is called “bulk” tunneling, with an exponent α_{bulk} = (*g* + 1/*g* − 2)/4. Taking *L* = 17μm and *g* = 0.7 gives Δ*E* ≈ 150 mK for our device, and therefore almost all of our ZBA data shown in Fig. 4 are expected to be in the bulk-tunneling regime. The measured values of α* _{T}* = 0.45 and α

*= 0.52 (device B) give*

_{V}*g*≈ 0.28 and 0.26, respectively.

The values of *g* found from the ZBA exponents are much smaller than that extracted from the holon branch in Fig. 3. We briefly offer possible explanations for this [with more detail in (*14*)]. One possibility is that impurities and imperfections may make the effective length of each wire shorter than the lithographic length. It is well known that an impurity in a TLL will effectively cut it into two at low energies. Although at finite energies, electrons are able to tunnel through such constrictions, one might expect that the constrictions will modify the form of the ZBA exponent in a similar way to “ends.” If we use the formula for ZBA exponent in the end-tunneling regime (*14*), the extracted values of *g* (0.53 from α* _{T}* and 0.49 from α

*) are more comparable to that extracted from the holon branch.*

_{V}Alternatively, it may not be valid to extrapolate the higher-energy properties of the branches to low energies where the ZBA is measured. The interactions between the wires and/or between the wires and the 2DEG may become important, and may alter the form of the ZBA exponent (*21*). The effective dielectric constant of the material may also be energy dependent, changing the strength of the interactions at low energies.

TLL theory is based on the assumption of a linear dispersion relation about the Fermi energy. We go beyond that regime at high dc bias. There is very little work on interactions in 1D wires at high energy where the TLL approximation breaks down. Haldane (*1*) argued that all 1D metals are adiabatically continuous with the TLL, that is, perturbations such as the band curvature are only expected to lead to a renormalization of the TLL parameters, so spin-charge separation should persist; this is backed up by recent renormalization-group calculations (*22*). Recent theoretical work on spinless fermions for a curved band has shown an intriguing interplay of Fermi-liquid and TLL behavior (*2*, *23*). Numerical calculations have been performed in this regime (*24*) using quantum Monte Carlo methods, with results that resemble our experimental results—a parabolic spin branch and a fairly straight charge branch.

## Supporting Online Material

## References and Notes

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵Materials and methods are available as supporting material on
*Science*Online. - ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- We acknowledge the UK Engineering and Physical Sciences Research Council for funding. Y.J. was supported by a scholarship from the Thai Ministry of Science. We thank F. Sfigakis for experimental assistance and B. E. Kardynal and C. H. W. Barnes for useful discussions.