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.

g F 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 F and g F . Particularly, simultaneous increase of the energy-averaged effective mass and that at the Fermi level was reported in earlier publications [10][11][12]19,20 . The most  Dependence of the field of complete spin polarization, B c , 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 c . 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 10 15 m −2 . The polarization field B c determined by the crossing of the tangents is marked by arrows. 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 F m, and Shubnikov-de Haas oscillations, yielding the product g F m F at the Fermi level (Fig. 1). We find that with decreasing electron density (or increasing interaction strength), the product g F m F at the Fermi level monotonically increases in the entire range of electron densities, while the energy-averaged product g F 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.
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 c , 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 c (n s ), in Fig. 2 for two samples. Over the electron density range 0.7 × 10 15 m −2 < n s < 2 × 10 15 m −2 , the data are described well by a linear dependence that extrapolates to zero at n c ≈ 0.14 × 10 15 m −2 (dashed black line).
33m e , and γ = −0.5%. The filling factors ν = n s 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 × 10 15 m −2 < n s < 2 × 10 15 m −2 , whereas at lower densities the quantity γ changing sign.
However, at lower electron densities down to n s ≈ 0.2 × 10 15 m −2 (up to r s ≈ 12), the experimental dependence B c (n s ) deviates from the straight line, and it linearly extrapolates to the origin. Note that the observed behaviour of the polarization field B c 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 c (n s ) with increasing electron density 22,23 . Therefore, the observed behaviour of B c is not related to the finite thickness effects.
The solid red line in Fig. 2 shows the polarization field B c (n s ) calculated using the quantum Monte Carlo method 24 . The experimental results are in good agreement with the theoretical calculations for the clean limit  k l 1 F (here k F 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 F 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 c (n s ) 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 · m 2 and is close to the slope 5.4 × 10 −15 T · m 2 observed in our experiment. However, the offset equal approximately to 0.3 × 10 15 m −2 of ref. 25 is appreciably higher compared to our case. Therefore, the behaviour of the polarization field B c is affected by the disorder potential in accordance with refs. 24,26 . Good agreement between our experimental data for B c 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 F 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 system where g v = 2 is the valley degeneracy and μ B is the Bohr magneton.
On the other hand, the Lande g-factor g F and effective mass m F at the Fermi level can be determined by the analysis of the Shubnikov-de Haas oscillations in relatively weak magnetic fields where T D is the Dingle temperature, T is the temperature, m e is the free electron mass, ħω c is the cyclotron splitting, Δ Z is the Zeeman splitting, and Δ v is the valley splitting. It is clear from the relation (2) that as long as one sets = Z 1 i v in the range of magnetic fields studied, the fitting parameters are T D m F , m F , and g F m F 27 . The values T D m F and m F 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 F m F . The quality of the fits is demonstrated in Fig. 3. The magnetoresistance δρ xx = ρ xx − ρ 0 normalized by ρ 0 (where ρ 0 is the monotonic change of the dissipative resistivity with magnetic field) is described well using Eq. (2).
The main result of our study shown in Fig. 1 is that the products of the average g F m and g F m F 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 F m F monotonically increases as the electron density is decreased in the entire range of electron densities, while the product g F m saturates at low n s . 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 F determined by the analysis of the temperature dependence of the amplitude of Shubnikov-de Haas oscillations. The similar behaviour of m F and g F m F 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 F SciENtific REPORTS | 7: 14539 | DOI:10.1038/s41598-017-15117-x (top inset to Fig. 1). At relatively high electron densities n s > 0.7 × 10 15 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 s < 0.7 × 10 15 m −2 gives rise to the flattening of the spectrum so that the effective mass at the Fermi level, m F = p F /V F , continues to increase (here V F is the Fermi velocity). In contrast, the energy-averaged effective mass does not, being not particularly sensitive to this flattening. Near the density n c , 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 F m F with temperature is indeed observed at the lowest-density point in Fig. 1. In the critical region, the increase of m F is restricted by the limiting value determined by temperature: m F < p F Δp/4k B T. In our experiments, the increase of m F reaches a factor of about two at n s = 0.3 × 10 15 m −2 and T ≈ 30 mK, which allows one to estimate the ratio Δ ∼ . p p / 006 F . 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 N 2 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 m 2 /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.