Indication of band flattening at the Fermi level in a strongly correlated electron system.

Using ultra-high quality SiGe/Si/SiGe quantum wells at millikelvin temperatures, we experimentally compare the energy-averaged effective mass, m, with that at the Fermi level, m F , and verify that the behaviours of these measured values are qualitatively different. With decreasing electron density (or increasing interaction strength), the mass at the Fermi level monotonically increases in the entire range of electron densities, while the energy-averaged mass saturates at low densities. The qualitatively different behaviour reveals a precursor to the interaction-induced single-particle spectrum flattening at the Fermi level in this electron system.

Introduction

The creation and investigation of flat-band materials is a forefront area of modern physics[1-4]. The interest is ignited, in particular, by the fact that, due to the anomalous density of states, the flattening of the band may be important for the construction of room temperature superconductivity. The appearance of a flat band is theoretically predicted[5-7] in a number of systems including heavy fermions, high-temperature superconducting materials, 3He, and two-dimensional (2D) electron systems. As the strength of fermion-fermion interaction is increased, the single-particle spectrum becomes progressively flatter in the vicinity of the Fermi energy eventually forming a plateau. The flattening of the spectrum is related to the increase of the effective fermion mass m at the Fermi level and the corresponding peak in the density of states.

The role of electron-electron interactions in the behaviour of two-dimensional (2D) electron systems increases as the electron density is decreased. The interaction strength is characterized by the Wigner-Seitz radius, r  = 1/(πn )1/2 a (here n is the electron density and a is the effective Bohr radius in semiconductor), which in the single-valley case is equal to the ratio of the Coulomb and kinetic energies. At high electron densities, where r ~1, interactions are not prevalent; 2D electron systems exhibit metallic conductivity and their behaviour can be well described by the Fermi liquid theory. In the opposite limit of low electron densities (r  > 35) interactions are dominant, and the electrons are expected to form a Wigner crystal[8,9]. At intermediate values of r , where the energy of interactions exceeds all other energy scales but is still not high enough to cause the electrons to crystallize, different scenarios of the system behaviour were suggested theoretically, and numerous experiments on different electron systems were performed (see reviews[10-13] and references therein). Nevertheless, the behaviour of the strongly interacting 2D electron system is still not well understood. One of the important points is the residual scattering potential. Even in quite perfect electron systems the electron properties are influenced by the residual disorder[14-16]. At present, it is unclear whether the clean limit (i.e., the limit in which the influence of disorder on the electron properties is negligible) has been reached in experiments on the least disordered electron systems. Finally, there exist experimental facts that contradict both intuitive expectations and calculations, e.g., a decrease of the effective mass in fully spin-polarized single-valley electron systems[17,18]. In view of these discrepancies it becomes imperative to obtain new experimental data on the behaviour of strongly interacting 2D electron systems as close as possible to the ideal ones in the sense of both strength of electron-electron interactions and weakness of the scattering potential.

Raw experimental data obtained in strongly correlated 2D electron systems can be divided into two groups: (i) data describing the electron system as a whole, like the magnetic field required to fully polarize electron spins, thermodynamic density of states, or magnetization of the electron system, and (ii) data related solely to the electrons at the Fermi level, like the amplitude of the Shubnikov-de Haas oscillations yielding the effective mass m and Lande g-factor g at the Fermi level. As a rule, the data in the first group are interpreted using the quasiparticle language in which the energy-averaged values of effective mass, m, and Lande g-factor, g, are used. To determine the values, the formulas that hold for the case of non-interacting electrons are employed. Although this approach is ideologically incorrect, the results for m and g often turn out to be the same as the results for m and g . Particularly, simultaneous increase of the energy-averaged effective mass and that at the Fermi level was reported in earlier publications[10-12,19,20]. The most pronounced effects were observed in the 2D electron system in silicon metal-oxide-semiconductor field-effect transistors (MOSFETs) in which the effective mass is strongly enhanced at low densities while the g-factor stays close to its value in bulk silicon, the exchange effects being small[10,11]. The strongly enhanced effective mass in Si MOSFETs was previously interpreted in favor of the occurrence of Wigner crystal or an intermediate phase (like ferromagnetic liquid). In fact, the origin and presence of possible intermediate phases preceding the formation of Wigner crystal can depend on the degree of disorder in the electron system. Because the mobility in SiGe/Si/SiGe quantum wells is two orders of magnitude higher compared to Si MOSFETs, the origin of the low-density phases in these electron systems can be different. Note that the experimental results obtained in the least-disordered Si MOSFETs exclude the metal-insulator transition driven by localization. The disorder effects in higher mobility SiGe/Si/SiGe quantum wells should be yet smaller.

In this work we have investigated the region of complete polarization of the electron spins in the ultra-high-mobility 2D electron system in SiGe/Si/SiGe quantum wells at millikelvin temperatures, yielding the energy-averaged product g m, and Shubnikov-de Haas oscillations, yielding the product g m at the Fermi level (Fig. 1). We find that with decreasing electron density (or increasing interaction strength), the product g m at the Fermi level monotonically increases in the entire range of electron densities, while the energy-averaged product g m saturates at low densities. Taking into account the negligibility of the exchange effects in the 2D electron system in silicon, this difference can only be attributed to the different behaviours of the two effective masses. Their qualitatively different behaviour reveals a precursor to the interaction-induced single-particle spectrum flattening at the Fermi level in this electron system.

Figure 1

Product of the Lande factor and effective mass as a function of electron density determined by measurements of the field of full spin polarization, B , (squares) and Shubnikov-de Haas oscillations (circles) at T ≈ 30 mK. The empty and filled symbols correspond to two samples. The experimental uncertainty corresponds to the data dispersion and is about 2% for the squares and about 4% for the circles. (g 0 = 2 and m 0 = 0.19 m are the values for noninteracting electrons). The top inset shows schematically the single-particle spectrum of the electron system in a state preceding the band flattening at the Fermi level (solid black line). The dashed violet line corresponds to an ordinary parabolic spectrum. The occupied electron states at T = 0 are indicated by the shaded area. Bottom inset: the effective mass m versus electron density determined by analysis of the temperature dependence of the amplitude of Shubnikov-de Haas oscillations, similar to ref.[34]. The dashed line is a guide to the eye.

The parallel-field magnetoresistance (ı.e., magnetoresistance measured in the configuration where the magnetic field is parallel to the 2D plane) has been studied in the perpendicular orientation of the magnetic field and the current. In this case, the shape of the experimental dependences is the closest to that expected theoretically[21]; the detailed data taken in both perpendicular and parallel orientations of the current with respect to the in-plane magnetic field will be published elsewhere. The field of the full spin polarization, B , corresponds to a distinct “knee” of the experimental dependences followed by the saturation of the resistance (see the bottom inset to Fig. 2). The magnetic field where the spin polarization becomes complete is plotted as a function of electron density, B (n ), in Fig. 2 for two samples. Over the electron density range 0.7 × 1015 m−2 < n  < 2 × 1015 m−2, the data are described well by a linear dependence that extrapolates to zero at n  ≈ 0.14 × 1015 m−2 (dashed black line). However, at lower electron densities down to n  ≈ 0.2 × 1015 m−2 (up to r  ≈ 12), the experimental dependence B (n ) deviates from the straight line, and it linearly extrapolates to the origin. Note that the observed behaviour of the polarization field B is quite opposite to that in the presence of the orbital effects due to the finite thickness of a 2D electron system, which are known to lead to a downward deviation in the dependence B (n ) with increasing electron density[22,23]. Therefore, the observed behaviour of B is not related to the finite thickness effects.

Figure 2

Dependence of the field of complete spin polarization, B , on electron density at a temperature of 30 mK for two samples (dots and squares). The dashed black line is a linear fit to the high-density data which extrapolates to zero at a density n . The solid red line corresponds to the calculation[24] for the clean limit. Top inset: the low density region of the main figure on an expanded scale. Also shown by the dotted red line is the calculation[24] taking into account the electron scattering. Bottom inset: the parallel-field magnetoresistance at a temperature of 30 mK at different electron densities indicated in units of 1015 m−2. The polarization field B determined by the crossing of the tangents is marked by arrows.

The solid red line in Fig. 2 shows the polarization field B (n ) calculated using the quantum Monte Carlo method[24]. The experimental results are in good agreement with the theoretical calculations for the clean limit (here k is the Fermi wavevector and l is the mean free path), assuming that the Lande g-factor, renormalized by electron-electron interactions, is equal to 2.4. Although in ref.[24] Lande g-factor was equal to 2, the reason for the 20% discrepancy between the theory and experiment may be due to the finite size of the electron wave function in the direction perpendicular to the interface. Besides, the product k l decreases with decreasing electron density, which leads to a downward deviation in the theoretical dependence, as shown by the dotted red line in the upper inset to Fig. 2.

To check whether or not the residual disorder affects the results for the magnetic field of complete spin polarization, we compare our data with those previously obtained on samples with an order of magnitude lower mobility[25]. At high electron densities, the dependence B (n ) in ref.[25] is also linear and extrapolates to zero at a finite density; the slope of the dependence is equal to 6 × 10−15 T · m2 and is close to the slope 5.4 × 10−15 T · m2 observed in our experiment. However, the offset equal approximately to 0.3 × 1015 m−2 of ref.[25] is appreciably higher compared to our case. Therefore, the behaviour of the polarization field B is affected by the disorder potential in accordance with refs.[24,26]. Good agreement between our experimental data for B and the calculations for the clean limit[24] gives evidence that the electron properties studied in our samples are only weakly sensitive to the residual disorder and the clean limit has been reached in our samples.

The product g m that characterizes the whole 2D electron system can be determined in the clean limit from the equality of the Zeeman splitting and the Fermi energy of the spin-polarized electron systemwhere g  = 2 is the valley degeneracy and μ is the Bohr magneton.

On the other hand, the Lande g-factor g and effective mass m at the Fermi level can be determined by the analysis of the Shubnikov-de Haas oscillations in relatively weak magnetic fieldswhere T is the Dingle temperature, T is the temperature, m is the free electron mass, ħω is the cyclotron splitting, Δ is the Zeeman splitting, and Δ is the valley splitting. It is clear from the relation (2) that as long as one sets in the range of magnetic fields studied, the fitting parameters are T m , m , and g m [27]. The values T m and m are obtained in the temperature range where the spin splitting is insignificant; the detailed results for the effective mass will be published elsewhere. Being weakly sensitive to these two fitting parameters, the shape of the fits at the lowest temperatures turns out to be very sensitive to the product g m . The quality of the fits is demonstrated in Fig. 3. The magnetoresistance δρ  = ρ  − ρ 0 normalized by ρ 0 (where ρ 0 is the monotonic change of the dissipative resistivity with magnetic field) is described well using Eq. (2).

Figure 3

Fits of the normalized magnetoresistance δρ /ρ 0 at a temperature of ≈30 mK (dots) using Eq. (2) with (a) g m /m  = 0.905, T  = 0.12 K, m  = 0.25 m , and γ = 2.5% and (b) g m /m  = 1.11, T  = 0.15 K, m  = 0.33m , and γ = −0.5%. The filling factors ν = n hc/eB ⊥ at minima are indicated. Inset: the asymmetry coefficient γ versus electron density for two samples. The Dingle temperature for two spin subbands is found to be different in our samples, similar to the results for Si MOSFETs of ref.[27]. Although the effect is appreciably weaker in our case, we have to introduce another fitting parameter γ for . The difference between the Dingle temperatures for two spin subbands does not exceed 6%, the Dingle temperature for energetically favorable spin direction being smaller over the range of electron densities 0.6 × 1015 m−2 < n  < 2 × 1015 m−2, whereas at lower densities the quantity γ changing sign.

The main result of our study shown in Fig. 1 is that the products of the average g m and g m at the Fermi level behave similarly at high electron densities, where electron-electron interactions are relatively weak, but differ at low densities, where the interactions become especially strong. The product g m monotonically increases as the electron density is decreased in the entire range of electron densities, while the product g m saturates at low n . We emphasize that it is the qualitative difference in the behaviours of the two sets of data that matters, rather than comparison of the absolute values. Taking into account the negligibility of the exchange effects in the 2D electron system in silicon[10,11], this difference can only be attributed to the different behaviours of the two effective masses. Their qualitatively different behaviour indicates the interaction-induced band flattening at the Fermi level in this electron system. To add confidence in our results and conclusions, we show in bottom inset in Fig. 1 the data for the effective mass m determined by the analysis of the temperature dependence of the amplitude of Shubnikov-de Haas oscillations. The similar behaviour of m and g m with electron density allows one to exclude any possible influence of the g-factor on the behaviour of the product of the effective mass and g-factor, which is consistent with the previously obtained results for the 2D electron system in silicon.

The experimental results are naturally interpreted within the concept of the fermion condensation[28-30] that occurs at the Fermi level in a range of momenta, unlike the condensation of bosons. With increasing strength of electron-electron interactions, the single-particle spectrum flattens in a region Δp near the Fermi momentum p (top inset to Fig. 1). At relatively high electron densities n  > 0.7 × 1015 m−2, this effect is not important since the single-particle spectrum does not change noticeably in the interval Δp and the behaviours of the energy-averaged effective mass and that at the Fermi level are practically the same. Decreasing the electron density in the range n  < 0.7 × 1015 m−2 gives rise to the flattening of the spectrum so that the effective mass at the Fermi level, m  = p /V , continues to increase (here V is the Fermi velocity). In contrast, the energy-averaged effective mass does not, being not particularly sensitive to this flattening. Near the density n , the electron system is in a critical region in which the effective mass at the Fermi level is expected to be temperature dependent. A weak decrease of the value g m with temperature is indeed observed at the lowest-density point in Fig. 1. In the critical region, the increase of m is restricted by the limiting value determined by temperature: m  < p Δp/4k T. In our experiments, the increase of m reaches a factor of about two at n  = 0.3 × 1015 m−2 and T ≈ 30 mK, which allows one to estimate the ratio . It is the smallness of the interval Δp that provides good agreement between the calculation[24] and our experiment.

Methods

Samples were prepared based on a 15 nm SiGe/Si/SiGe quantum well grown in an ultrahigh-vacuum chemical-vapor-deposition (UHVCVD) apparatus[31,32]. The Hall-bar samples with width 50 μm and the distance between potential probes 150 μm were patterned using standard photo-lithography. Contacts consisted of AuSb alloy, deposited in a thermal evaporator in vacuum and annealed for 5 minutes in N2 atmosphere at 440 °C. Next, approximately 300 nm thick layer of SiO was deposited in a thermal evaporator and a >20 nm thick Al gate was deposited on top of SiO. No mesa etching was used, and the 2D electron gas was created in a way similar to silicon MOSFETs (for details, see refs.[33,34]). The maximum electron mobility in our samples reached 240 m2/Vs which is the highest mobility reported for this electron system[34].

Measurements were made in an Oxford TLM-400 dilution refrigerator in a temperature range 0.03–1.2 K. The resistance was measured by a standard four-terminal lock-in technique in a frequency range 1–11 Hz. The applied currents varied in the range 0.5–4 nA. We used a saturating infra-red illumination to improve the quality of the contacts and increase the electron mobility. This did not affect the electron density at a fixed gate voltage.