# Thermodynamics of Two-Dimensional Electron Gas in a Magnetic Field

Authors: V. I. Nizhankovskii
Journal: Physics Research International
 (2011)

---

## Abstract

Change of the chemical potential of electrons in aGaAs-AlxGa1−xAs heterojunction was measured in magnetic fields up to 6.5 T at several temperatures from 2.17 to 12.3 K. A thermodynamic equation of state of two-dimensional electron gas well describes the experimental results.


---

## Body

## 1. Introduction

One of the basic problems in understanding physical properties of a two-dimensional electron gas (2DEG) is to know whether the electron densityN is dependent on a magnetic field or not.Useful information about the electron densityN may be obtained from measurements of the thermodynamic properties of a 2DEG, the magnetization, the specific heat, and the chemical potential. Early results for the specific heat [1], the magnetization [2, 3], and the chemical potential [4, 5] were interpreted in N=const assumption. In order to do this, it was necessary to introduce a background density of states between the Landau levels. According to [6] an effective background density of states between the Landau levels may result from a spatial inhomogeneity of 2DEG. Interpretation of more recent measurements of the magnetization [7] allows change of N in response to the changing magnetic field.The first measurement of the chemical potentialμ performed on GaAs-AlxGa1-xAs heterojunction at temperature T=4.2 K has shown that the experimental dependence of μ on magnetic field was greatly different from the theoretical one calculated for N=const. Moreover, the Hall resistance calculated from experimental data for μ had plateaus typical for the quantum Hall effect (QHE) and was in good agreement with a direct measurement [8].In the present work measurements of the chemical potential are made at several temperatures from 2.17 to 12.3 K and a thermodynamic description of the results is proposed. The experiment was performed many years ago but not published, because an explanation was absent.


## 2. Experiment and Results

The modulation-doped sample was grown by molecular beam epitaxy method on a single-crystal semi-insulatingGaAs substrate in the following sequence: undoped GaAs; 0.8 μm; undoped GaAs-AlxGa1-xAs; 10 nm; GaAs-AlxGa1-xAs doped with Si to a concentration 5·1017-1018 cm−3; 50 nm; undoped GaAs; 20 nm.The method to measure the change of the chemical potentialΔμ=μ(H)-μ(0) is based on the determination of the change in charge on the measuring capacitor consisting of the sample under investigation and reference electrode. The idea is closely related to Lord Kelvin's investigations of the contact potential difference. Let us make a capacitor from two different metals and short its plates with an electrometer (see insert in Figure 1). The measured charge will be equal to Q=CΔφ, where C is the capacitance and Δφ is the the contact potential difference between the plates (equal to the difference between their work functions). In general, the contact potential difference differs from the difference of the chemical potentials by the magnitude of the potential difference of the double charged layers present at the surface of bulk metals. If the influence of a magnetic field on the capacitance C and on the work function of the reference electrode is negligibly small, then ΔQ=CΔμ/e, where e is the elementary electronic charge.Figure 1
Change of the chemical potential measured at different temperatures. Insert: principle of the measurement (S: sample, R: reference electrode, and Q: electrometer).The reference electrode was made from bronze. The sample with area 1 cm2 was insulated from the reference electrode with thin capacitor paper; resulting capacitance was about 100 pF. The measuring capacitor was mounted at the end of the coaxial line connected to the electrometer. The coaxial line was built with sapphire insulators, and the input resistance of the electrometer was greater than 1016 Ohm. Calibration of the measuring capacitor was performed in a way correspondent to the classical definition of the electrostatic capacitance, namely, the charge on the capacitor was plotted as a function of slowly increasing and decreasing voltage V(t) applied to the sample (dV/dt<1 mV/s). The capacitance was determined from a slope of the obtained graph. This calibration was done at different values of applied magnetic field H including, of course, the quantum Hall plateaus regions. No change in the capacitance was observed within 1% accuracy.The temperature was measured with a thermocouple Cu-Cu(0.01%Fe) and stabilized with a precision of 0.05 K.The experimental results forΔμ are depicted in Figure 1. The chemical potential oscillates and amplitude of oscillations decreases with increasing temperature.The behavior of the chemical potential calculated forN=const is shown in Figure 2. These dependencies are numerical solutions of an equation for the density of 2DEG (spin-splitting is neglected):(1)N=∑n∫0∞D(ε,εn)[1+exp⁡(ε-μkBT)]-1dε
with a Gaussian broadened Landau levels(2)D(ε,εn)=eHπchΓexp⁡[-(ε-εn)22Γ2],εn=(n+12)eHℏm*c.Figure 2
Change of the chemical potential calculated forN=const. Insert: results for T=2.17 K and Γ=0.56 meV (a), Γ=1.2 meV (b).The incoming parameters were taken as follows: the effective massm*=0.068m, the electron density N=4.1·1011 cm−2, and the broadening Γ=0.56 meV. The electron density was determined from the period of oscillations in reciprocal magnetic field. The effective mass and broadening were obtained from the cyclotron resonance line position and width measured on this sample at T=2 K and λ=119 μm.From a comparison between Figures1 and 2 we find that the assumption N=const is broken at all temperatures. By doubling the Landau level broadening it is possible to adjust the calculated amplitude of oscillations of the chemical potential to the measured one, but the shape of calculated curve remains strongly different from the experiment (see insert in Figure 2). Thus results calculated for N=const differ from experiment not only quantitatively but qualitatively.Similar disagreement between calculated forN=const dependence and measured one was observed in de Haas-van Alphen effect experiments [2, 3, 7] and in studies of the chemical potential [5, 8].We also tried to approximate calculated atN=const dependencies of the chemical potential on a magnetic field to those obtained experimentally by introducing an energy-independent background density of states D=xD0 where D0=m*/πℏ2=2.84·1010 meV-1cm−2 is the density of states at H=0 [6]:(3)N=xD0μ+(1-x)∑n∫0∞D(ε,εn)[1+exp⁡(ε-μkBT)]-1dε.Visual agreement with the experiment was obtained forD=0.5D0≃1.4·1010 meV−1cm−2 (Figure 3). This value of the background density of states is much larger than theoretical one D≈0.3·1010 meV−1cm−2 estimated in the μ-Gaussian statistical model [6] and experimental upper limit for D=0.036D0 is deduced directly from the photoluminescence spectra [9].Figure 3
Change of the chemical potential calculated forN=const with the background density of states D=0.5D0.Thus assumption ofN=const is of little use for explanation of experimental dependencies of the chemical potential.


## 3. Discussion

Recently we have proposed a thermodynamic approach to the magnetostriction and successfully applied it for explanation of the experiments performed with nickel [10], terbium molybdate [11], and manganese fluoride [12]. These results for 3-dimensional system gave an idea to apply thermodynamics for explanation of the behavior of the field and temperature dependencies of the chemical potential μ measured on GaAs-AlxGa1-xAs heterojunction.The thermodynamic description is based on the fact that magnetic field does not perform any work. The Lorentz force acting on an electric charge in magnetic field is directed perpendicular to its velocity. Free spinss precess around the direction of magnetic field H keeping an angle between s and H independent of the field strength. Net magnetic moment of the sample arises only due to the “friction” of spins with the environment (relaxation effects).Thus for 2DEG subjected to a magnetic field the first law of thermodynamics is the same as without the field:(4)ΔU=TΔS+αΔA,
where ΔU=U(H)-U(0) is the change of the energy, ΔS=S(H)-S(0) is the change of the entropy, α is the surface tension, and A is the area. Typical for two-dimensional systems term αΔA corresponds to -pΔV in three dimensions [13]. Recalling definition of the Gibbs energy G(=Nμ)=U-TS-αA, we obtain, for T=const and the unit area, ΔG=-Δα or(5)Δ(Nμ)=-Δα.
This is the thermodynamic equation of state of a 2DEG.Appearance of the Gibbs energy is compatible with the experimental fact that both thermodynamic parameters—the chemical potentialμ and the electron density N—are functions of a magnetic field. For example, the change of the free energy ΔF is equal to the work done at fixed T, A, and N and the change of the thermodynamic potential ΔΩ to the one done at fixed T, A, and μ [13]. The change of the Gibbs energy ΔG is equal to the work done at constant T, α and the term Δα adds the work produced by a change of the surface tension α.BecauseU(0)=N(0)μ(0)/2, S(0)=(π2/3)N(0)T/μ(0), α is negative, which is a natural consequence of the electron-electron repulsion. We know nothing about its dependence on a magnetic field and temperature—this problem should be elaborated by a microscopic theory. Assuming that for i>2 change of the surface tension Δα≃0 (it is difficult to imagine that the electron's repulsion is dependent on relatively low magnetic field when electron trajectories are well overlapped, because for a given i the cyclotron orbit area Si=i/2N) we calculated the dependencies Δμ(H) shown in Figure 4. As shown in Figure 5 the calculated Hall resistance has the plateaus that are typical of the QHE.Figure 4
Change of the chemical potential calculated forNμ=const.Figure 5
The Hall resistance calculated atT=2.17 K for Nμ=const. The horizontal strokes indicate the positions of the plateaus in units of h/e2. Insert: i=4 plateau in enlarged scale.Figures1 and 4 are very similar. Now we understand that experimental curves contain artifacts caused by too fast field sweep (about 0.5 T/minute). Due to the large time constant of the electrometer (an old model with a vibrating capacitor was used), oscillations of Δμ in field below 1.5 T were suppressed. Else, fast field sweep results in hysteresis of Δμ(H) dependencies in the regions of QHE plateaus, that was established 2 years after these measurements, were performed [14]. Hysteresis is aroused by a Hall effect involving weakly damped eddy currents induced in 2DEG upon a change in the magnetic field. This hysteresis is seen in Figure 1 on falling sections of Δμ: crossings of zero at H~4.3 T do not coincide. Besides, amplitude of oscillations at T=2.17 K is smaller than it was expected. Having in mind these artifacts, the agreement between the experiment and theory should be considered as almost perfect.Thus usual thermodynamics explains dependencies of the chemical potential on the magnetic field and temperature without any additional assumptions or adjustable parameters.Similarly the smoothed oscillations of the magnetization   and the specific heat could be explained.Besides the thermodynamic properties, the proposed description gives the natural explanation of the Hall plateaus observed at microwave frequencies [15, 16] (comparable with the cyclotron frequency ωc=eH/m*c)—the plateaus result from oscillations of the total number of electrons. (Theory of the QHE for N=const predicts a critical frequency that is much smaller than ωc   [17, 18].)Of course, the equilibrium thermodynamics could not pretend to explain all transport properties typical for the QHE.



---

*Source: 742158-2011-02-07.xml*