# Fermi Liquid Model of Radiation Induced Magnetoresistance Oscillations in GaAs/AlGaAs Heterostructure Two-Dimensional Electron System: Theoretical Evidence of an Electron Reservoir

###### Abstract

The magnetoresistance oscillations in GaAs/AlGaAs heterostructures induced by millimeterwave radiation recently observed by Zudov et al. [Phys. Rev. B 64, 201311 (2001)] is theoretically reproduced by introducing a model based on the finite temperature Fermi liquid theory of dc conductivity and the electron reservoir hypothesis. The obtained oscillation patterns show excellent agreement with the experimental result by Zudov et al. and are independent of the polarization of the radiation field in accordance with the experimental observation by Smet et al. [Phys. Rev. Lett. 95, 116804 (2005)].

###### pacs:

73.21.-b ; 73.21.Fg ; 73.43.Qt^{†}

^{†}preprint:

^{†}

^{†}thanks: email: toyoda@keyaki.cc.u-tokai.ac.jp

The new class of magnetoresistance oscillations in a high-mobility two-dimensional electron system (2DES) in a GaAs/AlGaAs heterostructure subjected to weak magnetic fields and millimeterwave radiation found by Zudov et al. Zudov2001 ()Zudov2003 ()Zudov2004 () and by Mani et al. Mani2002 ()Mani2004 () has revealed again the fascinating nature of the 2DES Studenikin2005 (); Dorozhkin2005 (); Simovic2005 (); Studenikin2007 (); Wirthmann2007 (); Hatke2008 (), evoking the discoveries of the quantum Hall effects (QHE) by von Klitzing vonKlitzing1980 () and by Tsui et al. Tsui1982 (). The period of these new oscillations is governed by the ratio of the millimeterwave to cyclotron frequencies, and the minima of the oscillations are characterized by an exponentially vanishing diagonal resistance. When the millimeterwave radiation is turned off, the magnetoresistance shows the well-established SdHvA oscillation whose period is governed by the ratio of the chemical potential to the cyclotron frequency. Hence this magnetoresistance oscillation is apparently induced by the illumination with millimeterwaves. Although various different theoretical models have been proposed Ryzhy1969 (); Ryzhii1970 (); Ryzhii2003a (); Ryzhii2003b (); Ryzhii2003c (); Ryzhii2004 (); Durst2003 (); Sedrakyan2008 (); Lei2004 (); Mikhailov2004 (); Andreev2003 (); Dietel2005 (); Torres2005 (); Robinson2004 (); Inarrea (); Inarrea2007 (); Vavilov2004a (); Vavilov2004b (); Dmitriev (); Dmitriev2005 (); Dmitriev2008 (), one crucial question remains to be solved. The experiment by Smet et al. Smet2005 () shows that the resistance oscillations are notably immune to the polarization of the radiation field. This observation is discrepant with these theories and seems to cast doubt on the validity of the theoretical models so far proposed Zhang (). In these circumstances it seems necessary to explore other theoretical possibilities to explain the phenomena.

Meanwhile, an extremely interesting observation was made by Holland et al., who measured the long-wavelength magnetoplasmon dispersion in a high mobility 2DES realized in a GaAs/AlGaAs heterostructure quantum well Holland2004 (). Using the coupling between the plasmon and THz radiation, they explored a wide range of filling factors for a fixed value of the wave-number vector to obtain an explicit filling-factor dependence of the dispersion. The observed dispersion shows a clear filling-factor dependent plateau-type dispersion, violently deviating from the well-established semiclassical formula at first glance. Their discovery suggests a previously unknown relation between the magnetoplasmon and the integer QHE (IQHE). The full theoretical explanation of the phenomenon was given very recently Toyoda2008 () by adopting the electron reservoir hypothesis (ERH) that was introduced more than two decades ago to explain the IQHE Baraff1981 (); Toyoda1984 (); Toyoda1985 (). That the experimental discovery by Holland et al. can be explained by adopting the ERH strongly indicates the significance of taking into account the electron reservoir in order to investigate the behavior of the electrons in the GaAs/AlGaAs heterostructure 2DES. In this paper a Fermi liquid model based on the ERH is proposed to explain the new class of the magnetoresistance oscillations within the framework of the standard quantum theory of electrical conductivity. The derived formula gives the magnetoresistance oscillation patterns which show excellent agreement with the experimental data observed by Zudov et al., including the SdHvA oscillation part.

The Fermi liquid theory of electrical conductivity was originally formulated by Eliashberg. The core of the theory is the analytic continuation of the finite temperature current correlation function with respect to the Matsubara frequency to obtain the retarded real-time current response function, which directly yields the conductivity. Here the formulation given in Ref.Toyoda1989 () is applied to an electron gas confined in the -plane and subjected to a magnetic field and a millimeter wave whose angular frequency is . The polarization of the radiation field is arbitrary. The electron charge and effective mass are denoted by and , respectively. The interactions between the electrons and impurities or defects are treated perturbatively by the Bethe-Salpeter (BS) equation for the vertex function. The quasi-electron spectrum is assumed to be , where is the principal quantum number for the Landau levels, is the cyclotron frequency, is the effective -factor, is the Bohr magneton with the electron rest mass , and the function sgn gives the sign of the spin variable . The Landau levels are degenerate and there is an additional quantum number , which corresponds to the -component of the electron momentum. The allowed range of this momentum is .

Although a microscopic model for the scattering mechanism responsible for the conductivity may be introduced via the proper vertex function in the BS equation, here we do not use a specific model but simply assume that the scattering mechanism is independent of and . Then applying the theory formulated in Ref.Toyoda1989 () to the two-dimensional quasi-electrons, the general expression of conductivity can be found as

(1) |

where , () is the retarded (advanced) Green function, and is the vertex function. The aim of this paper is to examine the general consequence of this conductivity formula when the ERH is adopted.

The fundamental assumption of the Fermi liquid model, i.e., the Fermi liquid hypothesis (FLH), can be mathematically formulated in terms of the retarded and advanced Green functions, and . In Ref.Toyoda1989 () it has been shown that by virtue of the FLH the product in the above conductivity formula can be written as

(2) |

where is the wave function renormalization constant, and is the damping parameter, which should satisfy the condition , and is the chemical potential.

In the theories so far proposed Ryzhy1969 (); Ryzhii1970 (); Ryzhii2003a (); Ryzhii2003b (); Ryzhii2003c (); Ryzhii2004 (); Durst2003 (); Sedrakyan2008 (); Lei2004 (); Mikhailov2004 (); Andreev2003 (); Dietel2005 (); Torres2005 (); Robinson2004 (); Inarrea (); Inarrea2007 (); Vavilov2004a (); Vavilov2004b (); Dmitriev (); Dmitriev2005 (); Dmitriev2008 (), the effects of the millimeterwave radiation only on the 2DES were considered. However, if there is an electron reservoir, the effects of the radiation on the electrons in the reservoir should also be taken into account. At present the mechanism of the electron reservoir is not known. Here electrons in the reservoir are assumed to be in bound states with binding energy . Since the amount of the energy that an excited electron can receive from millimeterwave radiation is , the condition that the electron can join the 2DES should be , which can also be written as . As the chemical potential is the minimum free energy to add an electron to the system, the emergence of such a process may be described by introducing another singularity in the retarded Green function at . This singularity of the Green function may be expressed as an effective chemical potential in the expression such as (2). These considerations lead to the following expression for :

(3) |

when the radiation is on, and

(4) |

when the radiation is off. Here we have defined , , and regularized step functions and . The limit corresponds to the sharp cut-off. By virtue of the delta functions the integration over in (1) for the two terms in (3) becomes trivial. Since the damping takes place for non-zero value of the momentum, the main contribution to the -integration comes from the vicinity of . That is, we can safely approximate . Then the conductivity formula (1) yields

(5) |

when the radiation is on, and

(6) |

when the radiation is off. The quantities for =1 and 2 depend on the microscopic mechanism of the scattering process responsible for the conductivity. The contribution from the non-singular term is denoted by . For simplicity, we assume . Using the explicit form of , the conductivity can be expressed as

(7) |

when the radiation is on, and

(8) |

when the radiation is off. Here is assumed to be a constant, and ’s are given as

(9) |

This conductivity formula is derived within the theoretical framework of the linear response approximation, in which the electric field is an externally controlled small perturbation. In the measurement by Zudov et al. Zudov2001 () the current is measured by controlling the electric field , while the current as well as the external electric field are kept zero. Therefore, the resistivity observed in their measurement should correspond to in this theory. The significance of considering the boundary conditions on the current expectation values has been discussed in Ref.Toyoda1985 (), where the IQHE formula is derived from a microscopic many-body Hamiltonian without depending on the linear response approximation. Thus, within the theoretical framework of the present formulation the resistivity corresponding to in Ref.Zudov2001 () is given as

(10) |

The function depends on three adjustable parameters , , and , while depends only on . These parameters will be adjusted to optimize the fit of the theoretical curve with the experimental data.

The chemical potential of 2DES without millimeterwaves is given as Toyoda2008 ()

(11) |

where is the electron number density in the limit of very weak magnetic field. This gives a good approximation for the range of the parameters in the experiments of Ref.Zudov2001 (). It should be noted here that the above approximation gives the classical Hall effect instead of the IQHE. This explains the experimental observation reported in Ref.Zudov2004 ().

The -dependence of for kG is shown by a solid line in Fig. 1, using the values of the physical parameters given in Ref.Zudov2001 (), i.e., K, cm, , and GHz. The adjustable parameters are chosen such that kG, , and . These values indicate (meV). The experimentally measured taken from Fig. 4 in Ref.Zudov2001 () is also shown by a dotted line for the same values. The oscillation pattern of shows excellent agreement with the experimental result. In Fig. 2 the -dependence of is shown by a solid line for kG, using the same values for the parameters. The experimentally measured taken from Fig. 4 in Ref.Zudov2001 () is also shown by a dotted line for the same values. The theoretical pattern shows excellent agreement with the expermental curve. The phase of the oscillation observed in the experimental curve is shifted toward positive direction. The shift becomes prominent for kG. Theoretically, an additional term to the quasi-electron energy spectrum, i.e., the proper self-energy, may cause such a shift. The magnetoresistivity data given in Fig. 1 of Ref.Zudov2001 () show a peak due to the magnetoplasmon resonance around kG. The existence of magnetoplasmon resonance peak in the magnetoresistivity in Hall bar structures was first observed by Vasiliadou et al. Vasiliadou1993 () in the microwave photoconductivity measurement. The magnetoplasmon resonance gives considerable contribution to the two-particle Green function, which is directly related to the proper self-energy of the single particle excitation spectrum via the finite temperature Ward-Takahashi relation Toyoda1989 (). In Fig. 3, the function is plotted for GHz. It exhibits clear SdHvA oscillation observed experimentally.

In conclusion, the -dependence of the oscillatory patterns of the millimeterwave induced magnetoresistance oscillations observed by Zudov et al. is almost perfectly reproduced from our theoretical model based on the FLH and the ERH. Furthermore, it should be noted that the present model is independent of the polarization of the raditation field in perfect accordance with the experimental observation by Smet et al. Smet2005 (). Since the model does not depend on any specific mechanisms of the electron scatterings that cause electric resistivity, this new magnetoresistivity oscillation seems to be a universal consequence of the Fermi liquid nature of the 2DES and the existence of an electron reservoir in GaAs/AlGaAs heterostructures.

In view of the fact that the present theory can reproduce excellently the experimental data including the immunity to the polarization of the radiation field, the next task should be to reveal the microscopic mechanism that realizes the electron reservoir. At present there seem to be no direct experimental clues to that question. Nevertheless, it should be noted that there has been no direct experimental evidence against the existence of an electron reservoir either. In other words there has been no direct experimental evidence that the electron number is fixed.

As the essential theoretical conjecture of this work, i.e., Eqs. (3) and (4), can be used to calculate various physically measurable quantities, the validity of the model presented in this work may be further tested.

###### Acknowledgements.

I thank Chao Zhang for pointing out the significance of this theory on the immunity of the radiation induced magnetoresistance oscillation to the polarization of the radiation field and bringing Ref. Smet2005 () to my attention. I also thank Masaki Yasuè and Wolfgang Bentz for valuable comments.References

- (1) M.A. Zudov, D.R. Du, J.A. Simmons, and J.L. Reno, Phys. Rev. B 64, 201311 (2001).
- (2) M. A. Zudov, R. R. Du, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 90, 046807 (2003).
- (3) M. A. Zudov, Phys. Rev. B 69, 041304(R) (2004).
- (4) R. G. Mani, J. H. Smet, K. von Klitzing, V. Narayanamurti, W. B. Johnson, and V. Umansky, Nature 420, 646 (2002).
- (5) R. G. Mani, J. H. Smet, K. von Klitzing, V. Narayanamurti, W. B. Johnson, and V. Umansky, Phys. Rev. Lett. 92, 146801 (2004).
- (6) S.A. Studenikin, M. Potemski, A. Sachrajda, M. Hilke, L.N. Pfeiffer, and K.W. West, Phys. Rev. B 71, 245313 82005).
- (7) S.I. Dorozhkin, J.H. smet, V. Umansky, and K. von Klitzing, Phys. Rev. B 71, 201306 (2005).
- (8) B. Simovic, C. Ellenberger, K. Ensslin, H.-P. Tranitz, and W. Wegscheider, Phys. Rev. B 71, 233303 (2005).
- (9) S.A. Studenikin, A.S. Sachrajda, J.A. Gupta, Z.R. Wasilewski, O.M. Fedorych, M. Byszewski, D.K. Maude, M. Potemski, M. Hilke, K.W. West, and L.N. Pfeiffer, Phys. Rev. B 76, 165321 (2007).
- (10) A. Wirthmann, B.D. McCombe, D. Heitmann, S. Holland, K.-J. Friedland, and C.-M. Hu, Phys. Rev. B 76, 195315 (2007).
- (11) A.T. Hatke, H.-S. Chiang, M.A. Zudov, L.N. Pfeiffer, and K.W. West, Phys. Rev. B 77, 201304 (2008).
- (12) K. von Klitzing, G. Dorda, and M. Pepper, Phys. rev. Lett. 45, 494 (1980).
- (13) D.C. Tsui, H.L. Stormer, and A.C. Gossard, Phys. Rev. Lett. 48, 1559 (1982).
- (14) V.I. Ryzhii, Sov. Phys. Solid State 10, 2286(1969) ; 11, 2078 (1970).
- (15) A.D. Gladun and V.L. Ryzhyi, Sov. Phys. JETP 30, 534 (1070).
- (16) V. Ryzhii, Phys. Rev. B 68, 193402 (2003).
- (17) V. Ryzhii and V. Vyurkov, Phys. Rev. B 68, 165406 (2003).
- (18) V. Ryzhii and R. Suris, J. Phys.: Condens. Matter 15, 6855 (2003).
- (19) V. Ryzhii, J. Phys. Soc. Jpn 73, 1539 (2004).
- (20) A.C. Durst, S. Sachdev, and S.M. Girvin, Phys. Rev. Lett. 91, 086803 (2003).
- (21) T.A. Sedrakyan and M.E. Raikh, Phys. Rev. Lett. 100, 086808 (2008).
- (22) X.L. Lei, J. Phys.: Condens. Matter 16, 4045 (2004).
- (23) S.A. Mikhailov, Phys. Rev. B 70, 165311 (2004).
- (24) A.V. Andreev, I.L. Aleiner, and A.J. Millis, Phys. Rev. Lett. 91, 056803 (2003).
- (25) J. Dietel, L.I. Glazman, F.W.J. Hekking, and F. von Oppen, Phys. Rev. B 71, 045329 (2005).
- (26) M. Torres and A. Kunold, Phys. Rev. B 71, 115313 (2005).
- (27) J.P. Robinson, M.P. Kennett, N.R. Cooper, and V.I. Fal’ko, Phys. Rev. Lett. 93, 036804 (2004).
- (28) J. Inarrea and G. Platero, Phys. Rev. Lett. 94, 016806 (2005); Phys. Rev. B 72, 193414 (2005); Phys. Rev. B 76, 073311 (2007).
- (29) J. Inarrea, Appl. Phys. Lett. 90, 172118 (2007).
- (30) M.G. Vavilov, I.A. Dmitriev, I.L. Aleiner, A.D. Mirlin, and D.G. Polyakov, Phys. Rev. B 70, 161306 (2004).
- (31) M.G. Vavilov and I.L. Aleiner, Phys. Rev. B 69, 035303 (2004).
- (32) I.A. Dmitriev, A.D. Mirlin, and D.G. Polyakov, Phys. Rev. Lett. 91, 226802 (2003); Phys. Rev. B 70, 165305 (2004); Phys. Rev. B 75, 245320 (2007); Phys. Rev. Lett. 99, 206805 (2007).
- (33) I.A. Dmitriev, M.G. Vavilov, I.L. Aleiner, A.D. Mirlin, and D.G. Polyakov, Phys. Rev. B 71, 115316 (2005).
- (34) I.A. Dmitriev, F. Evers, I.V. Gornyi, A.D. Mirlin, D.G. Polyakov, and P. Wolfe, phys. stat. sol. (b) 245, 239 (2008).
- (35) J. H. Smet, B. Gorshunov, C. Jiang, L. Pfeiffer, K. West, V. Umansky, M. Dressel, R. Meisels, F. Kuchar, and K. von Klitzing, Phys. Rev. Lett. 95, 116804 (2005).
- (36) C. Zhang, private communication.
- (37) S. Holland, Ch. Heyn, D. Heitmann, E. Batke, R. Hey, K. J. Friedland, and C.-M. Hu, Phys. Rev. Lett. 93, 186804 (2004).
- (38) T. Toyoda, N. Hiraiwa, T. Fukuda, and H. Koizumi, Phys. Rev. Lett. 100, 036802 (2008).
- (39) G. A. Baraff and D. C. Tsui, Phys. Rev. B 24, 2274 (1981).
- (40) T. Toyoda, V. Gudmundsson, and Y. Takahashi, Phys. Lett. 102A, 130 (1984)
- (41) T. Toyoda, V. Gudmundsson, and Y. Takahashi, Physica 132A, 164 (1985).
- (42) T. Toyoda, Phys. Rev. A 39, 2659 (1989).
- (43) E. Vaisliadou, G. Müller, D. Heitmann, D. Weiss, K. v. Klitzing, H. Nickel, W. Schlapp, and R. Lösch, Phys. Rev. 48, 17145 (1993).

Fgiure captions

Fig. 1

The -dependence of the function and the experimentally measured for when the system is illuminated with radiation. The parameters are: 100 GHz, K, cm . (a) Theoretical defined by Eq. (7) is plotted by the solid line. (b) Experimentally measured taken from Fig. 4 in Ref. Zudov2001 () is plotted by the dotted line.

Fig. 2

The -dependence of the function and the experimentally measured for when the system is illuminated with radiation. The parameters are same as Fig. 1. (a) Theoretical defined by Eq. (7) is plotted by the solid line. (b) Experimentally measured taken from Fig. 4 in Ref. Zudov2001 () is plotted by the dotted line.

Fig. 3

The -dependence of the function for when the system is not illuminated with radiation. It shows a typical SdHvA oscillation.