The spectrum of allowed energy of electron in graphene superlattice in the quantizing magnetic field is investigated. Such spectrum consists of number of so-called magnetic minibands. The width of these minibands depends on the superlattice barriers power and on the magnetic field intensity. The explicit form of electron spectrum is derived in the case of weak magnetic field. The possibility of electron-electron and electron-phonon bound states is shown. The binding energies of these states are calculated. The binding energy is shown to be the function of magnetic field intensity.
1. Introduction
The interactions of electrons with elementary excitations in crystal lattice have fundamental implications on properties of materials and lead to such many-body phenomena as superconductivity and charge-density waves. Such interactions take an unusual form in graphene [1–3] which can be described by the Dirac-like form of the effective Hamiltonian. A new area for investigations of chiral massless fermions is induced by discovery of this 2D material. So, interactions of Dirac electrons with elementary excitations which lead to the emergence of the bound states (BS) are under the intensive theoretical [4–7] and experimental study [8–10] last time.
Renormalization of Dirac spectrum due to interactions of electrons with lattice vibrations was obtained in [11–15]. Formations of polarons, plasmon-phonon complexes, and electron-hole pares are investigated in [4, 16–20], where BSs were shown to appear if energy of quasiparticles interaction exceeds a threshold value. In [19], the modification of Dirac spectrum near the threshold of optical phonon emission was studied. The spectrum characteristics of the electron-phonon quasiparticles (in particular, the electron-phonon binding energy versus the coupling parameter gel-ph) were found in [19] within the theory, taking into account the singular vertex corrections beyond perturbation theory. Such corrections lead to the spectrum that remains below the optical phonon energy and corresponds to the electron-phonon BS in comparison with that obtained within the Wigner-Brillouin perturbation theory [11–13].
Influence of various factors (geometry parameters, spin-orbit interaction, charge impurities, etc.) on BS properties in graphene-based structures was studied in [4, 5, 21–27]. From an application technological point of view, the tunability of graphene electronic and optical properties by external fields is of particular importance [21, 28–32]. So, the effect of electric and magnetic fields on the properties of the BSs in graphene structures is of high interest among the researchers now [20, 21, 33–35]. The influence of Landau level mixing on the energy of magnetoexcitons and magnetoplasmons in graphene was studied in [36], where dispersions of the excitations were shown to be changed by virtual transitions between Landau levels, caused by Coulomb interaction (electron-hole channel). In [20], the effect of the magnetic field on the electron-phonon BSs was studied. The electron-phonon binding energy was shown to diverge at the resonant values of the magnetic field intensity. Such divergences corresponded to the electron-phonon hybrid states formed in the spectrum between the graphene Landau levels. These resonance states were obtained within the perturbation theory from the poles of the two-particle Green function and determine the structure of the magnetophonon resonance [37–39].
Presently, among the different graphene structures, the special attention is paid to graphene superlattices (GSLs) [40–44]. Intensive investigations of electric and optical properties of GSL are explained by their different possible experimental and technological applications [31, 45–49]. Below we obtain the binding energies of electron-electron and electron-phonon coupling in GSL in the quantizing magnetic field and determine the binding energy dependence on the magnetic field intensity and on the coupling parameters gel-el and gel-ph.
2. Electron Spectrum of GSL in Quantizing Magnetic Field
In this section, we obtain the effective electron spectrum of GSL in quantizing magnetic field. We consider a GSL which is in the plane xy (Figure 1). Near the Dirac point, the electron spectrum of GSL has the explicit form [40]:(1)εGSLp=uF2px2+Δ2sin2pyd2ħ,where Oy is the GSL axis (Figure 1), uF=υFsinη/η is the Fermi velocity along the axis Ox, η=V0d/4ħυF is the dimensionless power of SL barriers [40], V0 is the barrier height, d=2·10-6 cm is the GSL period, Δ=2ħυF/d, and υF=108 cm/s is the Fermi velocity in graphene. Linearized equation (equation analogous to the Dirac equation), which describes the motion of an electron in the GSL, takes the form [50, 51]:(2)uFσxπ^x+Δσysinπ^yd2ħψ=εψ.Here, σ=σx,σy,σz are the Pauli matrixes, ψ is spinor function describing the electron states in GSL, π^=p^+eA/c, p^=-iħ∂x,∂y, A=0,Hx is the vector potential, and H is the intensity of magnetic field which is supposed to be perpendicular to the GSL plane. Indeed, in the absence of magnetic field, the substitution of spinor for free charge carrier [52] into (2) gives the dispersion law (1).
Schematic of the studied system.
In nonzero quantizing magnetic field (ħυF/λH≫T, λH=cħ/eH is the Larmor radius, and T is the temperature) we represent ψ=eikyχ, where χ is spinor which is equal to(3)χ=χ1xχ2x.χ1 and χ2 describe electron states in the first and second graphene sublattice correspondingly. After acting on the spinor ψ twice with the operator(4)H^GSL=uFσxπ^x+Δσysinπ^yd2ħ,using the next correlations H^GSL2ψ=ε2ψ, σxσy=-σyσx=iσz, and making a shift of the argument kλH2+x→x, we derive the following:(5)-ħ2uF2∂x2χ+ΔUxχ=ε2-ΔΔHσzχ.Here, we define ΔH=2πħuF/ΛH, ΛH=4πλH2/d, and(6)Ux=Δsin22πxΛH-2ΔHsin2πxΛHσz.
To find the eigen values of ε in (5) for the considered direction of vector H with respect to GSL axis, we use the method developed in [53]. It is easy to see from (6) that Ux+ΛH=Ux and therefore χx satisfies the Bloch theorem. Thus, we can write(7)χx=∑neinqΛHχ0x-nΛH,where χ0 is the eigenspinor of the next operator H^02=-ħ2uF2∂x2+ΔUx and corresponds to the electron state near the nth minimum of the function Ux [53]. Generally, the spectrum of allowed energy ε consists of a number of minibands (magnetic minibands [53]). Having made some transformations, we derive(8)ε2-χ0H^02χ0=4ħ2uF2χ0+ΛH2∂xχ0ΛH2cosqΛH.
For the weak magnetic fields ΔH≪Δ, the minima of Ux are separated by the wide barriers. In this case, the components of eigenspinor χ0 for the operator H^02 are close to the functions of the harmonic oscillator localized near the nth minimum of the function Ux. Further, we consider that electron is in the lowest Landau level. Thus, we obtain instead (8)(9)εq=ΔΔH-DHcosqΛH,where(10)DH=8πΔexp-2π2ΔΔH.The width of allowed miniband is seen from (10) to depend on the GSL barriers power and on the magnetic field intensity.
3. Emergence of the BS of Two Electrons in GSL
Here, we show that in GSL with electron spectrum (9) in the presence of quantizing magnetic field the appearance of BS of two electrons is possible at arbitrarily weak repulsion between them. Equation that defines the dispersion law of the BS can be obtained by determining the poles of the two-particle Green function [54]:(11)gel-el2π∫-π/ΛH+π/ΛHdqE-ε1q-p/2-ε2q+p/2+i0=1,where gel-el is the effective coupling parameter of quasiparticles with dispersion lows ε1q and ε2q and p is the momentum of BS. We use (11) in the case of electron-electron interaction in GSL in the presence of quantizing magnetic field. Equation (11) has undamped solutions if gel-el>0 (repulsion) and E>2ΔΔH+2DH. After integration in (11), we obtain the dispersion law of BS:(12)Ep=2ΔΔH+gel-el2ΛH2+4DH2cos2pΛH2.It is seen from (12) that electrons can bind if energy E is close to the top of conduction band where effective masses are negative. This fact explains the possibility of electrons binding while gel-el>0. Binding energy is(13)E0=gel-el2ΛH2+4DH2-2DH.
For weak repulsion between electrons gel-el≪DHΛH, the binding energy is proportional to square of gel-el: E0=gel-el2/4DHΛH2. For intensive repulsion between electrons gel-el≫DHΛH, we have E0=-2DH+gel-el/ΛH. The binding energy dependence on the magnetic field intensity is shown in Figure 2 for different values of coupling parameter.
Binding energy E0 versus the magnetic field intensity H. (a) gel-el=2ħuF, (b) gel-el=1.5ħuF, and (c) gel-el=ħuF.
4. Emergence of the Electron-Phonon BS in GSL
Now, we obtain the possibility of electron-optical phonon BS in GSL in the quantizing magnetic field. The equation that determines the dispersion law of such BS is [54](14)gel-ph2π∫-π/ΛH+π/ΛHdqE-ħω0-εq+i0=1.Here, the coupling parameter gel-ph<0 and corresponds to the attraction between electron and phonon and ω0 is the frequency of optical phonon which is assumed to be dispersionless. The undamped solution of (14) is possible if E<ħω0+ΔΔH-DH. Dispersionless of optical phonon leads to the fact that electron-phonon BS appears at all values of quasiparticle momentum. In this case, the binding energy is(15)E0=gel-ph2ΛH2+DH2-DH.In comparison with that obtained for electron-phonon coupling in graphene in the absence of magnetic field [19], the binding energy (15) does not exhibit exponential dependence on the coupling constant gel-ph.
5. Conclusion
Above we have investigated the spectrum of allowed energy of electron in GSL in the presence of quantizing magnetic field ħυF/λH≫T. Generally, such spectrum consists of number of so-called magnetic minibands. The width of these minibands depends on the GSL barriers power and on the magnetic field intensity. For our calculations, we assume that electron is in the magnetic miniband which corresponds to the lowest harmonic oscillator state. The explicit form of electron spectrum in this case is (9). The formula (9) is adequate if magnetic field satisfies the next condition: ΔH≪Δ, which is true up to H=15 T (the power of GSL barriers is of η~3.136 [40]). For such magnetic fields, temperature should not exceed 100 K. Nitrogen boiling point is suitable to this end. Moreover, the upper limit of magnetic field can be increased by changing the power of GSL barriers η.
Thus, experimental measurements could confirm the suggested theory if one uses the GSL with the period d=2·10-6 cm (typical value for superlattices [41, 45]) and with the height of barrier V0~0.392 eV (in this case, η~3.136). To this end, the GSL obtained by a sheet of graphene deposited on a banded substrate formed by periodically alternating layers of SiC and hexagonal BN is suitable one [41]. Another suitable structure could be obtained by periodical insertion of hydrogenated graphene (graphane [55]) between gapless graphene strips. For such material barrier height consists of V0~2.7 eV [41, 55]. If V0~2.7 eV, then η~22 and the condition ΔH≪Δ will be also performed up to H=15 T.
Using (9), we have shown the possibility of electron-electron and electron-phonon BSs and have calculated the binding energy of these states, which is the function of magnetic field (Figure 2). In the case of electron-optical phonon coupling, the binding energy does not exhibit exponential dependence on the coupling constant gel-ph in comparison with that obtained for graphene in the absence of magnetic field [19]. Electron-phonon BS appears at all values of quasiparticle momentum. This feature is explained by the dispersionless of optical phonon.
In [20], the electron-phonon binding energy was shown to diverge in graphene at the resonant values of the magnetic field intensity. Such divergences correspond to the electron-phonon hybrid states formed in the spectrum between the graphene Landau levels. As it is shown above in GSL electron-phonon binding energy (15) has the finite maximum as well as electron-electron binding energy (Figure 2). It is obtained instead of divergence of the binding energy as it was for the case of single-layer graphene [20]. Finite maximum of the electron-electron and electron-phonon binding energy in GSL in quantizing magnetic field is due to the broadening of energy levels (Landau levels). Such broadening takes place if quantizing magnetic field is applied perpendicularly to the superlattice axis [53].
In the case of electron-electron coupling, BSs correspond to the Bose-Einstein statistics. At low temperatures (under the critical temperature), such electron system can form a superfluid gas [56]. The presence of additional potential in graphene in this case is essential for the following reason. Electrons can bind if energy E is close to the top of conduction band of GSL where effective masses are negative (12). This fact explains the possibility of electrons binding while gel-el>0. In the absence of additional potential (at the linear dispersion limit), there are no regions with negative effective masses. So, “repulsive” electron-electron interaction cannot form bound states in such limit.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The work was supported with the funding of the Ministry of Education and Science of the Russian Federation within the base part of the State task no. 2014/411 (Project code 522).
HwangE. H.Das SarmaS.Dielectric function, screening, and plasmons in two-dimensional graphene2007752020541810.1103/physrevb.75.2054182-s2.0-34347363039GangadharaiahS.FaridA. M.MishchenkoE. G.Charge response function and a novel plasmon mode in graphene20081001616680210.1103/PhysRevLett.100.1668022-s2.0-42549114664BursetP.HerreraW.YeyatiA. L.Proximity-induced interface bound states in superconductor-graphene junctions2009804041402(R)10.1103/physrevb.80.041402BreyL.FertigH. A.Elementary electronic excitations in graphene nanoribbons2007751212543410.1103/physrevb.75.1254342-s2.0-34047092825GuptaK. S.SenS.Bound states in graphene20092429910710.1142/S0217732309028886ZBL1165.813242-s2.0-65249161761BostwickA.SpeckF.SeyllerT.HornK.PoliniM.AsgariR.MacDonaldA. H.RotenbergE.Observation of plasmarons in quasi-freestanding doped graphene20103285981999100210.1126/science.11864892-s2.0-77953078917ApalkovV. M.ChakrabortyT.Gap structure of the Hofstadter system of interacting Dirac fermions in graphene20141121717640110.1103/PhysRevLett.112.1764012-s2.0-84899741388YanJ.ZhangY.KimP.PinczukA.Electric field effect tuning of electron-phonon coupling in graphene2007981616680210.1103/physrevlett.98.1668022-s2.0-34247189037PisanaS.LazzeriM.CasiraghiC.NovoselovK. S.GeimA. K.FerrariA. C.MauriF.Breakdown of the adiabatic Born-Oppenheimer approximation in graphene20076319820110.1038/nmat18462-s2.0-33847655495BostwickA.OhtaT.SeyllerT.HornK.RotenbergE.Quasiparticle dynamics in graphene200731364010.1038/nphys4772-s2.0-33845998961ParkC.-H.GiustinoF.CohenM. L.LouieS. G.Velocity renormalization and carrier lifetime in graphene from the electron-phonon interaction200799808680410.1103/physrevlett.99.0868042-s2.0-34548250859CalandraM.MauriF.Electron-phonon coupling and electron self-energy in electron-doped graphene: calculation of angular-resolved photoemission spectra2007762020541110.1103/physrevb.76.2054112-s2.0-36049027406TseW.-K.Das SarmaS.Phonon-induced many-body renormalization of the electronic properties of graphene2007992323680210.1103/PhysRevLett.99.2368022-s2.0-36849071047BaskoD. M.AleinerI. L.Interplay of Coulomb and electron-phonon interactions in graphene200877404140910.1103/physrevb.77.0414092-s2.0-38349148384CappellutiE.BenfattoL.Vertex renormalization in dc conductivity of doped chiral graphene200979303541910.1103/physrevb.79.0354192-s2.0-60949109436StauberT.PeresN. M. R.Effect of Holstein phonons on the electronic properties of graphene200820505500210.1088/0953-8984/20/5/0550022-s2.0-43049096822StojanovićV. M.VukmirovićN.BruderC.Polaronic signatures and spectral properties of graphene antidot lattices2010821616541010.1103/physrevb.82.1654102-s2.0-78149342806JablanM.SoljačićM.BuljanH.Unconventional plasmon-phonon coupling in graphene2011831616140910.1103/physrevb.83.1614092-s2.0-79961106869BadalyanS. M.PeetersF. M.Electron-phonon bound state in graphene2012852020545310.1103/physrevb.85.2054532-s2.0-84861823509ZhuJ.BadalyanS. M.PeetersF. M.Electron-phonon bound states in graphene in a perpendicular magnetic field20121092525660210.1103/PhysRevLett.109.2566022-s2.0-84871548390RecherP.NilssonJ.BurkardG.TrauzettelB.Bound states and magnetic field-induced valley splitting in gate-tunable graphene quantum dots200979808540710.1103/physrevb.79.0854072-s2.0-61549114667RyuS.MudryC.HouC.-Y.ChamonC.Masses in graphenelike two-dimensional electronic systems: topological defects in order parameters and their fractional exchange statistics2009802020531910.1103/physrevb.80.2053192-s2.0-77950853233XuJ. G.WangL.WengM. Q.Quasi-bound states and Fano effect in T-shaped graphene nanoribbons20131141515370110.1063/1.48241832-s2.0-84886550912KlöpferD.De MartinoA.EggerR.Bound states and supercriticality in graphene-based topological insulators201331142710.3390/cryst30100142-s2.0-84900465331AsmarM. M.UlloaS. E.Spin-orbit interaction and isotropic electronic transport in graphene20141121313660210.1103/PhysRevLett.112.1366022-s2.0-84898076320IadecolaT.NeupertT.ChamonC.Topological gaps without masses in driven graphene-like systems2014891111542510.1103/physrevb.89.1154252-s2.0-84897875960AsmarM. M.UlloaS. E.Symmetry-breaking effects on spin and electronic transport in graphene2015911616540710.1103/physrevb.91.165407AbergelD. S. L.ChakrabortyT.Long-range coulomb interaction in bilayer graphene2009102505680710.1103/PhysRevLett.102.0568072-s2.0-61349141303ApalkovV. M.ChakrabortyT.Electrically tunable charge and spin transitions in Landau levels of interacting Dirac fermions in trilayer graphene201286303540110.1103/physrevb.86.0354012-s2.0-84863647072ZareniaM.PartoensB.ChakrabortyT.PeetersF. M.Electron-electron interactions in bilayer graphene quantum dots2013882424543210.1103/physrevb.88.2454322-s2.0-84892386742PieperA.HeinischR. L.WelleinG.FehskeH.Dot-bound and dispersive states in graphene quantum dot superlattices2014891616512110.1103/physrevb.89.1651212-s2.0-84899751599MargapotiE.StrobelP.AsmarM. M.SeifertM.LiJ.SachsenhauserM.CeylanÖ.PalmaC.-A.BarthJ. V.GarridoJ. A.Cattani-ScholzA.UlloaS. E.FinleyJ. J.Emergence of photoswitchable states in a graphene-azobenzene-au platform201414126823682710.1021/nl503681z2-s2.0-84916608188PoundA.CarbotteJ. P.NicolE. J.Effects of electron-phonon coupling on Landau levels in graphene201184808512510.1103/physrevb.84.0851252-s2.0-80052497335SkinnerB.ShklovskiiB. I.VoloshinM. B.Bound state energy of a Coulomb impurity in gapped bilayer graphene201489404140510.1103/physrevb.89.0414052-s2.0-84893141437RakytaP.TóváriE.CsontosM.CsonkaS.CsordásA.CsertiJ.Emergence of bound states in ballistic magnetotransport of graphene antidots2014901212542810.1103/physrevb.90.1254282-s2.0-84907235922LozovikY. E.SokolikA. A.Influence of landau level mixing on the properties of elementary excitations in graphene in strong magnetic field20127, article 1341910.1186/1556-276x-7-1342-s2.0-84857731404AndoT.Magnetic oscillation of optical phonon in graphene200776202471210.1143/jpsj.76.0247122-s2.0-33847397507GoerbigM. O.FuchsJ.-N.KechedzhiK.Fal'koV. I.Filling-factor-dependent magnetophonon resonance in graphene20079908740210.1103/PhysRevLett.99.087402KimY.PoumirolJ. M.LombardoA.KaluginN. G.GeorgiouT.KimY. J.NovoselovK. S.FerrariA. C.KonoJ.KashubaO.Fal'KoV. I.SmirnovD.Measurement of filling-factor-dependent magnetophonon resonances in graphene using Raman spectroscopy20131102222740210.1103/physrevlett.110.2274022-s2.0-84878399706BarbierM.VasilopoulosP.PeetersF. M.Extra Dirac points in the energy spectrum for superlattices on single-layer graphene201081707543810.1103/physrevb.81.0754382-s2.0-77954877026RatnikovP. V.Superlattice based on graphene on a strip substrate200990646947410.1134/S00213640091801432-s2.0-76349105906Briones-TorresJ. A.Madrigal-MelchorJ.Martínez-OrozcoJ. C.Rodríguez-VargasI.Electrostatic and substrate-based monolayer graphene superlattices: energy minibands and its relation with the characteristics of the conductance curves2014739811210.1016/j.spmi.2014.05.0282-s2.0-84902438375SongJ. C. W.ShytovA. V.LevitovL. S.Electron interactions and gap opening in graphene superlattices20131112626680110.1103/physrevlett.111.2668012-s2.0-84891598598BolmatovD.MouC.-Y.Graphene-based modulation-doped superlattice structures2011112110210710.1134/S10637761110100432-s2.0-79952764148FalkovskyL. A.Optical properties of graphene and IV-VI semiconductors200851988789710.1070/pu2008v051n09abeh0066252-s2.0-58449105705PedersenT. G.PedersenJ. G.Transport in graphene antidot barriers and tunneling devices20121121111371510.1063/1.47688442-s2.0-84871246662PedersenJ. G.GunstT.MarkussenT.PedersenT. G.Graphene antidot lattice waveguides2012862424541010.1103/physrevb.86.2454102-s2.0-84871040286PedersenJ. G.PedersenT. G.Hofstadter butterflies and magnetically induced band-gap quenching in graphene antidot lattices2013872323540410.1103/physrevb.87.2354042-s2.0-84878992811ThomsenM. R.BrunS. J.PedersenT. G.Dirac model of electronic transport in graphene antidot barriers2014263333530110.1088/0953-8984/26/33/3353012-s2.0-84905454383KruecklV.RichterK.Bloch-Zener oscillations in graphene and topological insulators20128511543310.1103/physrevb.85.115433KryuchkovS. V.Kukhar'E. I.Solitary electromagnetic waves in a graphene superlattice under influence of high-frequency electric field201470708110.1016/j.spmi.2014.03.0082-s2.0-84897443729BursetP.YeyatiA. L.BreyL.FertigH. A.Transport in superlattices on single-layer graphene2011831919543410.1103/physrevb.83.1954342-s2.0-79961125478ZhylichA. G.The energy spectrum of electrons and optical properties of the superlattice in a magnetic field1994363792805PitaevskiiL. P.Weakly bound excitation states in a crystal1976432382388LebègueS.KlintenbergM.ErikssonO.KatsnelsonM. I.Accurate electronic band gap of pure and functionalized graphane from GW calculations2009792424511710.1103/physrevb.79.2451172-s2.0-67650113174UiminG. V.MaksimovL. A.BarabanovA. F.Bound states of electrons (holes) on a square lattice and their superfluidity1988485312316