Numerical and analytical results of the investigation of the thermodynamic properties of a quasi-2D electron gas are presented. The density of states, the temperature derivative of the chemical potential, and the heat capacity of the gas at the resonance points and away from it are analyzed. It is shown that, in the dependence of the heat capacity on the chemical potential, there are additional steps at the resonance points. The width of the additional steps increases with temperature. With the increase in temperature, when the main steps are practically blurred, there are still visible marks from the additional steps.
State Research Foundation of UzbekistanВА-ФА-Ф-5-020(ВА-ФА-Ф-2-005)1. Introduction
A quasi-two-dimensional electron gas in heterostructural quantum wells causes both fundamental and applied interest, since the spatial quantization of the carrier energy leads to the manifestation of interesting low-dimensional effects [1–8]. Of particular interest is the behavior of a two-dimensional electron gas in a magnetic field [9–15], since the latter can significantly change the properties of the system.
Under certain conditions, the effects of size quantization appear in the observed thermodynamic quantities even in the absence of a magnetic field [16–23].
It is known that the main reason for the manifestation of size-quantum effects in a gas is a stepwise change in the density of states at critical points, i.e., at the points of resonance, where the chemical potential μ is compared with the energy levels En of spatial quantization [1]. Such a change in the density of states manifests itself in one form or another in all the observed thermodynamic quantities. The change in the chemical potential relative to the energy levels of spatial quantization can be observed depending on the concentration ns or width of the quantum well L.
In [18], the entropy, which is incident on a single particle s=(∂S/∂ns)T=-(∂μ/∂T)ns in a quasi-two-dimensional electron gas, is studied as a function of the chemical potential μ. It turned out that at low temperatures it has narrow peaks at the resonance points μ = En (n = 1, 2, 3...), and the height and width of these peaks depend on the temperature and the scattering event (degree of disorder).
In this paper, the heat capacity of a quasi-two-dimensional ideal gas is studied as a function of the chemical potential at different temperatures. Based on known thermodynamic relationships, the density of states, the temperature derivative of the chemical potential, and the heat capacity of a two-dimensional electron gas are analyzed. For these quantities, low-temperature asymptotic formulas are obtained.
Numerical simulations show that, in the dependence of the heat capacity on the chemical potential, at the resonance points μ = En (n = 1, 2, 3...) there are additional steps. The width of the additional steps is proportional to the temperature and they are connected by sharp changes in the value dμ/dT at the resonance point. In conclusion, the results of the work are briefly summarized.
2. Basic Relationships
Let us consider an ideal electron gas with a simple parabolic spectrum. Energy is measured from the bottom of the zone of a massive semiconductor. In a heterostructural quantum well, the energy of the transverse motion is quantized, and the electron gas becomes quasi-two-dimensional. One-particle spectrum can be represented in the form [1](1)E=ħ2k22m+En,where k2=kx2+ky2-wave vector, En are the energy levels of the transverse motion, and is the effective mass of the electron. The total two-dimensional concentration is(2)ns=g0∑n=1∫En∞fEdE=g0∑nTln1+eμ-En/T,where g0=m/πħ2 and f(E) is the Fermi-Dirac distribution function. The thermodynamic density of a state [9] is defined as(3)gs=∂ns∂μ=g0∑n1eEn-μ/T+1
Expressions for the heat capacity of a gas can be obtained by differentiating the total energy(4)Us=g0∑n∫En∞EfEdE
by temperature, i.e., Cs=∂Us/∂T. We represent it in the form(5)Cs=C1+dμdTC2
here x=(E-μ)/2T. A temperature derivative of the chemical potential dμ/dT=-B/A can be found from the conditions of constancy of the number of particles (2), where(8)B=∑nln1+eμ-En/T+En-μTeEn-μ/T+1(9)A=∑n1eEn-μ/T+1
First, let us analyze the low-temperature asymptotics of these formulas.
(i) Let the chemical potential μ be near the level EN, and the temperature is sufficiently low μ ~ EN, T << EN − EN-1. Then in the expressions (6) and (7) all the summands with μ ~ EN < En exponentially small and N−1 the term with μ ~ EN > En give a finite contribution. By retaining also the N - term (since μ ~ EN), from (6) and (7) we obtain (N > 1)(10)C1=g0N-1∫-∞∞xμ+2xTcoshx2dx+g0∫EN-μ/2T∞xμ+2xTcoshx2dxC2=g02N-1∫-∞∞μ+2xTcoshx2dx+g02∫EN-μ/2T∞μ+2xTcoshx2dx
The first integrals in these expressions are calculated analytically. Then(11)C1=g0N-1Tπ23+∫EN-μ/2T∞xμ+2xTcoshx2dx(12)C2=g0N-1μ+12∫EN-μ/2T∞μ+2xTcoshx2dx
The remaining integrals in (11) and (12) can also be calculated analytically, but the results are cumbersome. These integrals can be calculated numerically with good accuracy.
Similarly, near the resonance point μ~EN under the condition T<<EN-EN-1 of (8) and (9) we have(13)B=ln1+eμ-EN/T+EN-μTeEN-μ/T+1A=N-1+1eEN-μ/T+1Thus, near the resonance point μ~EN under the condition T<<EN-EN-1 the quantity is equal d μ / dT to(14)dμdT=-eEN-μ/T+1ln1+eμ-EN/T+EN-μ/TeEN-μ/T+1N-1+1In the approximation considered, for the density of states (3) we can obtain the following formula:(15)gs=g0N-1+1eEN-μ/T+1The quantity dμ/dT also determines the entropy per electron s=(∂S/∂ns)T=-(∂μ/∂T)ns, and the low-temperature dependence s(μ) is studied in detail in [18]. Using (11), (12), and (14), according to formula (5), we can find the specific heat value near the resonance point μ=EN, and also the line shape Cs(μ) near μ~EN.
(ii) Suppose that the chemical potential μ is between the levels EN-1 and EN, i.e., EN-1 < μ < EN and the temperature is low enough T << EN − EN-1. Then in the expressions (11) and (12) the values of the integrals are exponentially small. Neglecting them, we obtain(16)C1=g0N-1Tπ23(17)C2=g0N-1μ
Under the conditions considered, the quantity dμ/dT in formula (14) is exponentially small. Thus, according to (5), the specific heat in this limit has the form(18)Cs=C1=g0N-1Tπ23
and the density of state (15) reduces to the form(19)gs=g0N-1
It can be seen from the comparisons (18) and (19) that when the condition EN-1<μ<EN, T<<EN-EN-1 and T<<EN-μ, both the density of the state and the heat capacity of a quasi-two-dimensional ideal gas gs, vary abruptly with increasing chemical potential μ(20)3Csπ2T=gs=N-1(3) At low temperatures at the resonance points T<<EN-EN-1 and μ=EN from (11) and (12) we have(21)C1=g0N-1Tπ23+∫0∞xEN+2xTcoshx2dx(22)C2=g0N-1μ+12∫0∞EN+2xTcoshx2dx
The integrals in these expressions are calculated analytically. Then(23)C1=g0Tπ23N-12+ENln2(24)C2=g0ENN-12+Tln2
Similarly, it follows from (14) that(25)dμdT=-ln2N-1/2Thus, at the resonance point the heat capacity (5) is(26)Cs=C1+dμdTC2=g0N-12Tπ231-3ln22π2N-1/22≈g0N-12Tπ23For N=2, the second term in the curly bracket is approximately equal 0,065 but decreases rapidly with growth. It is approximately equal to half the sum of the successive steps of the heat capacity (18) obtained under the condition EN-1<μ<EN. In a similar way, from (15) we can find an expression for the density of states at the resonance point(27)gs=g0N-12And in this limit the jumps in the specific heat (26) are due to a change in the density of states (27)(28)3Csπ2T≈gs=N-12
3. Numerical Results and Discussion
At not too low temperatures, the dependence of the heat capacity on the chemical potential can be obtained only numerically, on the basis of relations (5) - (9). The explicit form of the electron spectrum in a heterostructural quantum well can be determined by solving the quantum-mechanical problem. In the future, in numerical calculations, we use its simple approximation En=E1n2, n=1,2,3…
First, consider the dependence of the density of states gs(μ) and the derivative dμ/dT of the chemical potential μ. The results of numerical simulations of these quantities using formulas (3) and (8) and (9) are E1=0.05eV shown in Figure 1.
Dependences (a) of the density of states gs(μ)/g0 and (b) a derivative of the chemical potential |dμ/dT| at E1=0.05eV different temperatures: T=0.001eV (11.6 K), solid lines, and T=0.01eV (116 K), dashed lines.
The nature of the dependence gs(μ) is well known [2]. According to (3), the total density of the state is equal to the sum of the density of states of sequentially located overlapping minibands. Sharp jumps are observed at sufficiently low temperatures, when the chemical potential crosses the bottom of the miniband.
With increasing temperature, sharp step changes gs(μ) turn into a smoother dependence (dashed line in Figure 1(a)). As can be seen from the relations (3), this is due to the smoothing of the Fermi-Dirac step near E ~ μ the temperature increase.
The nature of the dependence of dμ / dT on the chemical potential was studied in detail in [18]. As can be seen from the graph (solid line in Figure 1(b)), at sufficiently low temperatures at resonance point μ = En the value dμ / dT has sharp peaks. According to formula (25), the height of these peaks depends on the quantum level number N and is independent of temperature and effective mass. The height of the peaks decreases with increasing quantum number N. As the temperature increases, the widths of these peaks increase (the dashed line in Figure 1(b)). Inside the miniband EN-1<μ<EN, where the density of state is constant (at T→0) and the derivative of d μ / dT is exponentially small.
As can be seen from the figure, at low temperatures the steps of gs(μ) and the peaks in dμ/dT are well described by the asymptotic formulas (19), (25), and (27). Near the resonance point, the shape of the lines can also be described with the help of relations (14) and (15).
Let us now consider the dependence of the specific heat on the chemical potential Cs(μ). The results of numerical calculations based on the relations (5)-(9), at E1 = 0 .05 eV are shown in Figure 2.
Dependences of the heat capacity of a quasi-two-dimensional electron gas on the chemical potential at E1=0.05eV and different temperatures: (a) T=0.001eV (11.6 K), (b) T=0.005eV (58 K), (c) T=0.01eV (116 K), and (d) T=0.025eV (290 K).
T=0.001 eV (11.6 K)
T=0.005 eV (58 K)
T=0.01 eV (116 K)
T=0.025 eV (290 K)
According to Figure 2(a), at sufficiently low temperatures, the specific heat varies abruptly, i.e., practically repeating the dependence of the density of states (see Figure 1(a)). This is also seen from the analytical result (20). The peak width of derivative d μ / dT at the resonance point is also small ~ T and is not noticeable in the graph.
With increasing temperature, additional weak steps appear at resonance points μ=En, Figure 2(b). With further temperature increase, the width of these steps increases, which is more noticeable in the graph of Figure 2(c). We note that as the temperature increases, the steps blur simultaneously. According to Figure 2(d), at temperature T=290 K, the main steps are practically blurred and the traces of additional steps are still noticeable.
As can be seen from the figure, at low temperatures, far from the resonance point EN-1 < μ < EN, the steps in dependence Cs(μ) are well described by formula (18) and at the resonance point by formula (28).
The low-temperature value of the heat capacity at the resonance point 3Cs(μ=EN)/(π2T)≈N-1/2 in Figures 2(b)–2(d) is represented by a horizontal dot line, which are obtained for N=2,3. A small deviation from the numerical result is due to the approximation admitted in (26).
4. Conclusions
In this paper, based on a simple model, the results of an analysis of the heat capacity of a quasi-two-dimensional ideal gas are presented. The calculations do not take into account many factors, such as the scattering process, electron-electron interaction, and so on. which can significantly affect the results.
In spite of this, it was found that in a simple model of a quasi-two-dimensional ideal gas, depending on the specific heat from the chemical potential other than the basic ones-ladder-like steps that are connected by the peculiarities of the density-of-states variations, there are still additional steps-the values of dμ / dT due to negative peaks at the resonance points. The width of the additional steps increases with temperature. With an increase in temperature, when the main steps are practically blurred, there are still visible marks from the additional steps (see Figure 2(d)).
Asymptotic analyzes of the density of state gs(μ), derivative dμ/dT, and heat capacity Cs(μ) at low temperatures, both at resonance points and far from it, are carried out. Simple formulas are obtained for them. Such singularities of heat capacity, obtained on the basis of a simple model, serve as a test when considering more realistic models.
Data Availability
The results of this study are included in the article and can be made freely available. Also, the authors embedded original figure in the text manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the State Research Foundation of Uzbekistan (Grant no. ВА-ФА-Ф-5-020(ВА-ФА-Ф-2-005)). The authors would like to thanks Professor Gafur Gulyamov for useful discussion results of the manuscript.
AndoT.FowlerA. B.SternF.Electronic properties of two-dimensional systems198254243767210.1103/RevModPhys.54.437BastardG.Wave Mechanics Applied to Semiconductor Heterostructures1988BastardG.BruntJ. A.FerreiraR.1991Acad PressBeenakkerC.van HoutenH.Quantum Transport in Semiconductor Nanostructures199144Elsevier1228Solid State Physics10.1016/S0081-1947(08)60091-0ButcherP.MarchN. H.TosiM. P.1993New York, NY, USADattaS.1995Cambridge University PressImryY.1997Oxford Univ. PressBakuevaL.MusikhinS.SargentE. H.ShikA.Multi-layer contacts for organic light-emitting diodes with enhanced injection efficiency200170810.1557/PROC-708-BB3.29vo KlitzingK.DordaG.PepperM.New method for high-accuracy determination of the fine-structure constant based on quantized hall resistance198045, article 49410.1103/physrevlett.45.4942-s2.0-33745141474TsuiD. C.StormerH. L.GossardA. C.Two-dimensional magnetotransport in the extreme quantum limit19824822155915622-s2.0-3594901906510.1103/PhysRevLett.48.1559WharamD. A.ThorntonT. J.NewburyR.PepperM.AhmedH.FrostJ. E.HaskoD. G.PeacockD. C.RitchieD. A.JonesG. A.One-dimensional transport and the quantisation of the ballistic resistance1988218L209L21410.1088/0022-3719/21/8/002Van WeesB. J.Van HoutenH.BeenakkerC. W. J.WilliamsonJ. G.KouwenhovenL. P.Van Der MarelD.FoxonC. T.Quantized conductance of point contacts in a two-dimensional electron gas19886098488502-s2.0-000152985410.1103/PhysRevLett.60.848ZawadzkiW.53SpringerSeries in Solid-State SciencesVagnerD.Thermodynamics of Two-Dimensional Electron on Landau Levels20063102152RamosA. C. A.FariasG. A.AlmeidaN. S.Thermodynamics of a quasi-two-dimensional electron gas: Effects of magnetic fields, temperature and finite width20114310187818812-s2.0-7996116877510.1016/j.physe.2011.06.030GokhfeldV. M.On the thermodynamics of a quasi-2D electron gas200531769773GlazerN. V.ErmolaevA. M.RashbaG. I.TheM. A.saltman.The saltman. Heat capacity two-dimensional electron gas with a superlattice2011963151520VarlamovA. A.KavokinA. V.GalperinY. M.Quantization of entropy in a quasi-two-dimensional electron gas2016931510.1103/PhysRevB.93.155404KukushkinI. V.MeshkovS. V.TimofeevV. F.Density of states two-dimensional electrons in a transverse magnetic field1988155219264KuntsevichA. Y.TupikovY. V.PudalovV. M.BurmistrovI. S.Strongly correlated two-dimensional plasma explored from entropy measurements20156, article 72982-s2.0-84934951485TsaranV. Y.KavokinA. V.SharapovS. G.VarlamovA. A.GusyninV. P.Entropy spikes as a signature of Lifshitz transitions in the Dirac materials201771102712-s2.0-85028686951GrassanoD.PulciO.ShubnyiV. O.SharapovS. G.GusyninV. P.KavokinA. V.VarlamovA. A.Detection of topological phase transitions through entropy measurements: The case of germanene2018972010.1103/PhysRevB.97.205442ShubnyiV. O.GusyninV. P.SharapovS. G.VarlamovA. A.Entropy per particle spikes in the transition metal dichalcogenides201844656156610.1063/1.5037559