Eur. Phys. J. B (2017) 90: 130 DOI: 10.1140/epjb/e2017-80130-8 THE EUROPEAN PHYSICAL JOURNAL B Regular Article
Density of states in the bilayer graphene with the excitonic pairing interaction
Vardan Apinyana and Tadeusz K. Kope´c Institute of Low Temperature and Structure Research, Polish Academy of Sciences, P.O. Box 1410, 50-950 Wroc law 2, Poland
Received 3 March 2017 / Received in final form 16 May 2017 Published online 5 July 2017 c The Author(s) 2017. This article is published with open access at Springerlink.com
Abstract. In the present paper, we consider the excitonic effects on the single particle normal density of states (DOS) in the bilayer graphene (BLG). The local interlayer Coulomb interaction is considered between the particles on the non-equivalent sublattice sites in different layers of the BLG. We show the presence of the excitonic shift of the neutrality point, even for the noninteracting layers. Furthermore, for the interacting layers, a very large asymmetry in the DOS structure is shown between the particle and hole channels. At the large values of the interlayer hopping amplitude, a large number of DOS at the Dirac’s point indicates the existence of the strong excitonic coherence effects between the layers in the BLG and the enhancement of the excitonic condensation. We have found different competing orders in the interacting BLG. Particularly, a phase transition from the hybridized excitonic insulator phase to the coherent condensate state is shown at the small values of the local interlayer Coulomb interaction.
1 Introduction ences [10–12]. A very large excitonic gap, found in a suffi- ciently broad interval of the repulsive interlayer Coulomb The problem of excitonic pair formation in the graphene interaction parameter, given in reference [12], suggests and bilayer graphene (BLG) systems is one of the long- that the level of intralayer density and chemical poten- standing and controversial problems in the modern solid tial fluctuations in the BLG system, discussed in refer- state physics [1–12]. The chiral invariance of the free quasi- ence [11], are sufficiently small and do not affect the ro- particle Hamiltonian, combined with electron Coulomb bust excitonic insulator state in the BLG. The excitonic interaction term brought the idea about the possibility condensation in the single BLG is principally possible in of spontaneous chiral symmetry breaking (CSB) in a sin- the case of the static interlayer screening regime [10]. gle layer and bilayer graphene, reflecting in the form of Moreover, it has been shown that there exists a crit- the gapped states in the fermionic quasiparticle spec- ical value of the interlayer hopping amplitude γ1 [15], trum [1,13–18]. The excitonic gap equation has been de- which provides an interesting energy cutoff, below which rived in references [5–11] at the Bardeen-Cooper-Schrieffer the electron-hole correlations do not drive the system to- (BCS) limit of the excitonic transition scenario. The gen- wards the CSB excitonic transition. More recently, it has eral idea used there is based on the supposition of the been shown [19,20] that the single-particle coherent den- weak electronic correlations in the BLG, due to its large sity of states (DOS) in the usual two-dimensional (2D) [19] number of the fermionic flavors. It has been suggested in and three-dimensional (3D) [20] semiconducting systems references [1,16,17] that even a partial account of the wave at the zero temperature limit is always finite, reflecting vector or energy dependence of the gap function rules out with the excitonic condensate regime in these systems. the constant gap solution in the BCS limit. Meanwhile, it For the 3D semiconducting systems, the coherent DOS has been shown that even undoped graphene can provide spectra survive also for the higher temperatures [20]. This a variety of electron-hole type pairing CSB orders espe- situation is also typical for the double layer electronic cially for the strong Coulomb coupling case [16,17], which structures at the half filling [21] when considering the renders the treatments in [1,13–15,18]tobeobscure. electron-hole pair formation and condensation. In addi- Concerning the excitonic condensation, the proper in- tion, it has been demonstrated that the excitonic insula- clusion of the chemical potential fluctuation effects on tor state and the excitonic condensation are two distinct the excitonic pairing gap and condensation in the BLG phase transitions in the solid state [22,23], and the conden- system, and also the role of the average chemical poten- sate states are due to the electronic phase stiffness [20,24], tial on the excitonic effects have been discussed in refer- mechanism. Meanwhile, in the electronic bilayer systems, those mentioned phase transitions are indistinguishable a e-mail: [email protected] as it was shown in references [12,21]. On the other hand, Page 2 of 12 Eur. Phys. J. B (2017) 90: 130 the inclusion of the dynamic screening and the full band electronic BLG structure with the equal chemical poten- structure favors the pairing and condensation [25–27]. tials and equal on-site quasienergies in each layer. When In this paper, we use the bilayer Hubbard model to switching the local Coulomb potential between the lay- calculate the single-particle DOS functions in the BLG ers, we keep the charge neutrality equilibrium across the and we examine the excitonic effects in the DOS. For the BLG, by imposing the half-filling condition in each layer. noninteracting layers, we show the principal modifications The interlayer Hubbard model with the intralayer U and to the usual tight-binding single layer graphene’s DOS be- local interlayer W Coulomb interaction terms is subjected havior (with differences between the A and B DOS struc- by the following Hamiltonian tures near the Dirac’s neutrality point) and we estimate the excitonic blue-shift values in the normal DOS for the H = −γ0 (¯aσ(r)bσ(r )+h.c.) reasonable values of the interlayer hopping amplitude. rr σ At finite values of the interlayer interaction parameter, the DOS shows always the remarkable four peaks struc- − γ0 a˜¯σ(r)˜bσ(r )+h.c. ture and a very large interband hybridization gap opens rr σ for intermediate values of the interaction parameter. We ¯ have found the critical value of the interlayer interaction − γ1 bσ(r)˜aσ(r)+h.c. − μnσ(r) Δ rσ rσ =1,2 parameter at which the hybridization gap Hybr opens in Δ the system. We estimate the values of Hybr for different + U [(nη↑ − 1/2) (nη↓ − 1/2) − 1/4] strengths of the interlayer interaction parameter and for r η=a,b different values of the interlayer hopping amplitude γ1. a,˜ ˜b We show the modifications of the van Hove singularity W n r − / n r − / − / . (vHs.) peaks in the DOS when augmenting the interlayer + [( 1bσ( ) 1 2) ( 2˜aσ ( ) 1 2) 1 4] interaction parameter. An interesting interlayer hopping rσσ mediated crossover, from the insulating pairing state to (1) the excitonic condensate state follows from our considera- tions. Moreover, we show that at any reasonable value of The first two terms describe the intralayer electron hop- the interlayer interaction parameter W , there exist a crit- ping with the hopping parameter γ0. The summation rr ical value of the interlayer hopping parameter γ1 at which , in these terms denotes the sum over the nearest the hybridization gap vanishes and the BLG pass to the neighbours lattice sites in the separated honeycomb lay- coherent excitonic condensate state. We calculate the hy- ers in the BLG structure. We kept here the small letters bridization gap in the system, and we show that it becomes a, b anda, ˜ ˜b for the electron operators on the lattice sites larger when increasing the interaction parameter W . A, B and A,˜ B˜ respectively. The third term corresponds Furthermore, we estimate the excitonic shift energies, to the hopping between two different monolayers in the the intraband vHs. peaks separations and the values of the BLG and the parameter γ1 is the interlayer hopping am- DOS at the Dirac’s neutrality points for the zero tempera- plitude. We neglect the nonlocal interlayer hopping terms ture case. The full interaction bandwidth, considered here, because their role is essentially unimportant in the consid- mimics different limits of correlations in the BLG, which ered problem about the excitonic effects. When working have been discussed only partially in the literature, where with the grand canonical ensemble, we have added also the the discussions have been restricted only to the strong in- chemical potential terms for each layer and for each sub- terlayer screening regime. lattice. We suppose here that the chemical potentials μ The paper is organized as follows: in Section 2,we corresponding to different layers =1, 2 and different sub- introduce the bilayer Hubbard model for our BLG system. lattices A, B and A,˜ B˜ are the same (this is, of course, true, In Section 3 we give the form of the fermionic action in the if we consider purely electronic layers and we neglect the Feynman’s path integral formalism. Next, in Section 4,we effect of the disorder, caused by the additionally charged discuss the single-particle DOS in the BLG, for different impurities). It is important to mention here that the chem- interlayer interaction regimes. In Section 5 we present the ical potentials of electrons on the nonequivalent sublattice numerical results on the excitonic effects in the normal sites, in the same layer, get different shifts in different lay- DOS and, in Section 6,wegiveashortconclusiontoour ers due to the stacking order of the BLG structure. Next, paper. nσ(r) is the electron density operator for the fermions in the layer and with the spin σ. The four U-terms in the Hamiltonian in equation (1) describe the on-site intralayer 2 The interlayer Hubbard model Coulomb interactions in the BLG. The summation index η in the U-terms, in equation (1), refers to different sublat- The Bernal Stacked (BS) BLG system is composed of two tice fermions in a given layer, i.e., η = a, b, for the bottom coupled honeycomb layers with sublattice sites A, B and layer with =1,andη =˜a,˜b, for the top layer with A˜, B˜, in the bottom and top layers respectively, arranged = 2. We will study the excitonic effects in the BLG at in the z-direction in such a way that the atoms on the sites the half-filling condition in each layer, i.e., n =1,for A˜ in the top layer lie just above the atoms on the sites B in =1, 2. The last term in equation (1), describes the lo- the bottom layer graphene, and each layer is composed of cal interlayer Coulomb repulsion, and the parameter W is two interpenetrating triangular lattices. We consider the the corresponding interlayer Coulomb interaction. We put Eur. Phys. J. B (2017) 90: 130 Page 3 of 12
γ0 = 1, as the unit of energy, and we set kB =1, =1 Here, Ωn = π (2n +1)/β with n =0, ±1, ±2,...,arethe through the paper. fermionic Matsubara frequencies [29]. The four component Dirac spinors ψσk(Ωn), in equation (6), are introduced at each discrete state k in the reciprocal space and for 3 The fermionic action and Dirac each spin direction σ =↑, ↓. Being the generalized Weyl representation spinors [30,31], they are given as T ˜ Here, we write the partition function of our electronic BLG ψσk(Ωn)= aσk(Ωn),bσk(Ωn), a˜σk(Ωn), bσk(Ωn) . system. For simplifying the notations we will introduce the (7) −1 two component fermionic fields, corresponding to different The matrix Gσk (Ωn), in equation (6), is the inverse type of fermions in the two sublattices of the graphene Green’s function matrix, of size 4 × 4. It is defined as monolayers: fc =(¯c, c), where c = a, b, for the layer with ⎛ ⎞ E1(Ωn) −γ˜1k 00 =1andc =˜a,˜b for the layer with = 2. Next, the ⎜ ⎟ ⎜ ∗ ¯ ⎟ partition function will be written as ⎜ −γ˜1k E2(Ωn) −γ1−Δσ 0 ⎟ G−1 Ω ⎜ ⎟ . σk ( n)=⎜ ⎟ (8) 0 −γ1−Δσ E2(Ωn) −γ˜2k −S f ,f ,f ,f ⎝ ⎠ Z Df Df Df Df e [ a b a˜ ˜b]. = [ aσ bσ] aσ˜ ˜bσ (2) −γ∗ E Ω σ 00 ˜2k 1( n) Indeed, the structure of the matrix does not changes when The fermionic action in the exponential, in equation (2), −1 −1 inverting the spin direction, i.e., G (Ωn) ≡ G (Ωn). is given by ↓k ↑k The diagonal elements of the matrix, in equation (8), are β the quasienergies S fa,fb,fa˜,f˜ = SB [fc]+ dτH (τ) . (3) b E Ω −iΩ − μeff , c=a,b 0 ( n)= n (9) a,˜ ˜b where the effective chemical potentials μ with =1, 2 eff eff are defined as μ1 = μ + U/4andμ2 = μ + U/4+W . We have introduced here the imaginary-time variables γ τ [28,29], at each lattice site r in both layers of the BLG. The parameters ˜k,inequation(8), are the renormal- τ ,β β /T ized (nearest neighbors) intralayer hopping amplitudes: The variables vary in the interval (0 ), where =1 γ zγ γ k γ with T being the temperature. The first term in equa- ˜k = k 0,wherethe -dependent parameters k are tion (3) describes the fermionic Berry-terms, correspond- the energy dispersions in the BLG layers. Namely, we have ing to all fermionic flavors in the BLG system. They are 1 γ e−ikδ . given as k = z (10) δ β ∂ S f dτc rτ c rτ . The parameter z, is the number of the nearest neigh- B [ c]= ¯σ( )∂τ σ( ) (4) rσ 0 bors lattice sites in the honeycomb lattice. The compo- δ nents of the nearest-neighbors vectors √, for the bot- Furthermore, in order to examine the excitonic effects in δ(1) d/ ,d / δ(2) tom layer√ 1, are given by 1 = 2 3 2 , 1 = the BLG, we perform the real-space linearization of the (3) d/2, −d 3/2 and δ1 =(−d, 0). For the layer 2, we four-fermionic terms, in equation (1). We do not present (1) (2) √ here the details of such a procedure and we refer to our have obviously δ2 =(d, 0), δ2 = −d/2, −d 3/2 ,and (3) √ recent work, in reference [12], where this procedure is pre- δ2 = −d/2,d 3/2 . Then, for the function γ1k,weget sented in more details. For the next, we will consider the √ homogeneous BLG structure, when the pairing occurs be- kxd 1 −ikxd i 3 γ1k = e +2e 2 cos kyd , (11) tween the particles with the same spin orientations, i.e., 3 2 Δσσ = Δσσδσσ and we can assume that the pairing gap is real Δσ = Δ¯σ ≡ Δ. Then the excitonic pairing gap where d is the carbon-carbon interatomic distance. It is ∗ ∗ parameter is given by not difficult to realize that γ2k = γ1k ≡ γk, and we have γ γ∗ ≡ γ∗ Δ W ¯b rτ a rτ . ˜2k =˜1k ˜k, where we have omitted the layer index . σσ = σ( )˜σ( ) (5) The form of the Green’s function matrix, given in equation (8), has been used recently in the work of ref- Next, we pass to the Fourier space represen- erence [12] in order to derive the self-consistent equa- tation, given by the transformation cσ(r,τ)= 1 i(kr−Ω τ) tions, which determine the excitonic gap parameter Δ cσk(Ωn)e n ,whereN is the total βN kΩn and the effective bare chemical potentialμ ¯ in the inter- η number of sites on the -type sublattice, in the layer , acting BLG system. Particularly, this last one plays an and we write the partition function of the system in the important role in the BLG theory and redefines the charge form neutrality point (CNP) in the context of the exciton for- mation in the interacting BLG [12]. Quite interesting ex- S ψ,¯ ψ 1 ψ¯ Ω Gˆ−1 Ω ψ Ω . = βN σk( n) σk ( n) σk( n) (6) perimental results on that subject are given recently in kΩn σ references [32–34]. Page 4 of 12 Eur. Phys. J. B (2017) 90: 130
4 The single-particle DOS partition function in equation (2), in Section 3, will be transformed as 4.1 The sublattice spectral functions − 1 ψ¯ (Ω )Gˆ−1(Ω )ψ (Ω ) βN kΩn σk n σk n σk n Z = Dψ¯Dψ e σ In this section, we calculate the single-particle normal 1 1 ¯ 1 ¯ J (Ωn)ψ (Ωn)+ ψ (Ωn)J (Ωn) DOS for the BLG system, given by the Hamiltonian, in βN kΩn [ 2 kσ kσ 2 kσ kσ ] × e σ equation (1) and basing on the form of the fermionic ac- 1 J¯ (Ω )Gˆ (Ω )J (Ω ) tion, given in equation (6). The normal single-particle 2 kΩn kσ n σk n kσ n ≈ e σ , DOS is straightforwardly defined as (18) J Ω 1 where the auxiliary source field vectors kσ( n)arethe ρ ω S k,ω , ψ c( )=N c( ) (12) subjects of the Dirac’s spinors defined similar to the - k fields, in Section 3 c T where the index corresponds to the sublattice type in the Jσk(Ωn)= jaσk(Ωn),jbσk(Ωn),jaσ˜ k(Ωn),j˜ (Ωn) . c A, B, A,˜ B˜ bσk BLG layers, i.e., = , and the spectral function (19) Sc(k,ω), in the right hand side in equation (12), is defined The matrix Gˆσk(Ωn) is defined as the inverse of the with the help of the retarded real-time Green’s function GR k,ω matrixgiveninequation(8) and we have Gˆσk(Ωn)= c ( )[28,29] −1 2 ˆ−1 βN Gσk (Ωn) . Furthermore, the calculation of the S k,ω − 1 GR k,ω . c( )= π c ( ) (13) thermal Matsubara Green’s function, in equation (17), is straightforward. For the sublattice-A, in the bottom layer Here, we have supposed that the spin variable σ is fixed of the BLG, we get in the direction spin-↑, being completely unimportant for δ2Z 1 the considered here problem. In turn, the retarded Green = − a¯k(Ωn)ak(Ωn) R δ¯jaσk Ωn δjaσk Ωn function Gc (k,ω), could be obtained after the analytical ( ) ( ) 4 continuation into the real frequency axis, in the expression βN 11 = − Gk (Ωn). (20) of the respective Matsubara Green’s function Gc(k,Ωn) 4 (see the similar procedures in Refs. [28,29]) Then it follows that R G (k,ω)=Gc(k,Ωn)| + . (14) 1 c iΩn →ω+i0 G11 Ω − a Ω a Ω . k ( n)= βN k( n)¯k( n) (21) The explicit calculation of the thermal normal Green’s functions Gc(k,Ωn) follows from the definition of the nor- After the inversion of the matrix given in equation (8), 11 mal Matsubara Green’s function [29]. As usually, in the the Green’s function matrix component Gk (Ωn)takesthe real space, they are defined as the statistical average of following form the product of an annihilation c and a creation typec ¯ 4 operators, i.e., for our fermions, we have αik G11 Ω , k ( n)= iΩ − ε (22) i=1 n ik Gc (rτ,r τ )=− c(rτ)¯c(r τ ) . (15) α(1) After transforming into the Fourier space the fermionic where the dimensionless coefficients ik are ⎧ operators, entering in equation (15) the local expression (3) ⎨ P (εik) 1 , i , , (ε −ε ) j=3,4 if =1 2 of the Green’s function (for the symmetry reasons of the i+1 1k 2k (εik−εjk) αik =(−1) (3) action, in Eq. (6), we consider the time- and space-local ⎩ P (εik) 1 , i , , (ε −ε ) j=1,2 if =3 4 expression of the single particle Green’s function), in equa- 3k 4k (εik−εjk) tion (15), is (23) P(3) ε ε G rτ,rτ − 1 G k,Ω , and ( ik) is a polynomial of third order in ik.Namely, c ( )= βN c ( n) (16) we have kΩn (3) 3 2 P (εik)=εik + ω1kεik + ω2kεik + ω3k (24) where the Fourier transforms Gc (k,Ωn)aregivenby with the coefficients ωik, i =1,...,3, given as 1 Gc (k,Ωn)= ck(Ωn)¯ck(Ωn) . (17) eff eff βN ω1k = −2μ2 − μ1 , (25) eff eff eff 2 2 ω2k = μ μ +2μ − (Δ + γ1) −|γ˜k| , (26) In order to calculate the statistical average on the right- 2 2 1 hand side in equation (17), we will perform the Hubbard- and Stratanovich transformation in the expression of the eff eff 2 eff 2 eff 2 partition function. In the Dirac’s spinor notations, the ω3k = −μ1 μ2 + μ1 (Δ + γ1) + μ2 |γ˜k| . (27) Eur. Phys. J. B (2017) 90: 130 Page 5 of 12
The quasiparticle energy parameters εik,inequation(23) are defined as follows 1 2 2 ε1,2k = Δ + γ1 ± (W − Δ − γ1) +4|γ˜k| − μ,¯ 2 (28) 1 2 2 ε3,4k = −Δ − γ1 ± (W + Δ + γ1) +4|γ˜k| − μ,¯ 2 (29) where we have introduced a new bare chemical potential
μeff + μeff μ¯ = 1 2 . (30) 2 As we will show later on, the bare chemical potential has a fundamental impact on the DOS behavior in the BLG, and provide the excitonic shift on the frequency axis. The energy parameters in equations (28)and(29) define the electronic band structure in the BLG system with the ex- Fig. 1. The exact solution of the effective chemical potential, citonic pairing interaction (see also in Ref. [12]). It is not normalized to the intralayer hopping amplitude γ0 as a func- difficult to verify that for the noninteracting BLG, i.e., tion of the interlayer interaction parameter. The inset shows when U =0,W =0andΔ = 0, the expressions in equa- the dependence ofη ¯ on the intralayer on-site interaction U,at tions (28)and(29) are reducing to the usual tight binding T = 0. The linear solution of the exact chemical potential μ is ε γ1 2 2 used in the calculations. The temperature is set at T =0. dispersion relations i = ± 2 ± (kγ0) +(γ1/2) − μ, 2 (with k = |γk| ), discussed in reference [35] in the context of the real-space Green’s function study of the noninter- BLG (see in Refs. [32–34]). Contrary, the intensity of the acting BLG. It is important to mention that the normal μ/γ spectral functions in different layers of the BLG, coincide function ¯ 0 is much higher in our case, due to the sin- with each other when interchanging the sublattice nota- gle BLG considered here, while in the gated double BLG tions in the monolayers. This follows from the symmetry measurements, discussed in reference [32], the interlayer of the action, given in equation (6). Particularly, we obtain interaction is much weaker as compared with the double monolayer graphene (see about in Ref. [10]), and the rea- S=2,c=˜a(k,ω)=S=1,c=b(k,ω), (31) son for this is the effect of the finite amount of carrier density nT induced in the top BLG reference [32]. S ˜(k,ω)=S=1,c=a(k,ω). (32) =2,c=b In Figure 2, the solution for the chemical potential is In Figure 1,wehavepresentedthevariationoftheeffec- shown as a function of the intralayer Coulomb interaction parameter U and for different values of the interlayer in- tive chemical potentialμ/γ ¯ 0 as a function of the interlayer W/γ teraction parameter W . The linear slope of the chemical Coulomb interaction parameter 0. The temperature κ . dependence has been also shown. In the inset, in Figure 1, potential corresponds to the coefficient =025, given in equation (33). In Figure 3, the exact solution for μ is we have shown the dependence ofμ/γ ¯ 0 on the local in- tralayer Coulomb interaction parameter U,givenby presented for different values of the intralayer interaction parameter U. The zero temperature case is considered in W the picture and the interlayer hopping amplitude is fixed μ¯(W, U, γ0,T)=μ(W, U, γ0,T)+κU + , (33) γ . γ 2 at the value 1 =0128 0. where μ(W, U, γ0,T) is the exact solution of the chem- ical potential in the BLG. Different fixed values of the 4.2 The sublattice density of states interlayer interaction parameters W are considered in the picture. The more detailed discussion on the role of the With the help of the analytical continuation, given in chemical potential is given in reference [12]. We see that equation (14) and by the use of the formula 1/(x + iδ)= the behavior of the effective chemical potential as a func- P(1/x)−iπδ(x), we get for the single particle normal DOS tion of W/γ0 shows a finite, very large jump at T/γ0 =0 for the sublattice A and the shifted CNP corresponds well with the recent ex- perimental observations of the behavior of bilayer’s av- ρ ω 1 S k,ω erage chemical potential, given in reference [32], where a A( )=N A( ) direct measurement of the chemical potential of BLG has k been done, as a function of its carrier density n. Here, it 4 1 α δ ω − ε . is important to mention that the excitonic effects in the = N ik ( ik) (34) BLG lead to the significant shift of the double CNP in the i=1 k Page 6 of 12 Eur. Phys. J. B (2017) 90: 130
Here, we can transform the summation over the wave vec- tors, into the integration over the continuous variables, by introducing the 2D density of states corresponding to the noninteracting graphene layers ρ x 1 δ x − γ . 2D( )=N ( k) (35) k
The noninteracting DOS in the monolayer graphene be- yond the Dirac’s approximation [36,37], can be analyti- cally expressed in terms of the elliptic integral of the first kind K(x)[38]. Namely, we have ⎧ ⎨√ 1 4|x/γ0| , < |x| <γ , |x| K ϕ(|x/γ |) 0 0 2 ϕ(|x/γ0|) 0 ρ2D(x)= π2|γ |2 ⎩ √ 1 ϕ(|x/γ0|) ,γ< |x| < ∞, 0 K 4|x/γ | 0 4|x/γ0| 0 (36)
where, formally, we have enlarged the domain of variation Fig. 2. The chemical potential solution of the BLG normal- of the argument, into ±∞. The function ϕ(x), in equa- ized to the intralayer hopping amplitude γ0 as a function if tion (36), is given by [36] the intralayer interaction parameter U. Different values of the W 2 2 interlayer interaction parameter are considered and the in- 2 x − 1 γ . γ ϕ(x)=(1+x) − . (37) terlayer hopping amplitude is set at 1 =0128 0.Thezero 4 temperature limit is considered. After the definition in equation (37), we get for the normal A DOS function Γ ρ2D [i (ω)] αj [i (ω)] T 0 0 ρ ω A( )= |Λ ω | i,j=1,2 − [ 1 ( )] 0 ρ ω α ω 2D [ i ( )] j [ i ( )] , + |Λ ω | (38) i,j=3,4 + [ 3 ( )] Γ Γ U 0 0 U 0 4 where the frequency dependent dimensionless parameters UΓ 1 UΓ 5 0 0 i(ω)withi =1,...,4, in equation (38), are given by the 2 UΓ0 2 UΓ0 3 following expressions 0 ⎧ Γ ⎪ ω μeff ω μeff − γ − Δ , Μ ⎪ + 1 + 2 1 5 ⎪ i+1 ⎨ if i =1, 2, 4 − ( 1) ω 3 i ( )= 4 0 |zγ | ⎪ 0 ⎪ Γ 2 ⎪ eff eff ⎪
ω μ ω μ γ1 Δ , ⎪ + 1 + 2 + + ⎩ Μ 1 0 if i =3, 4 1 (39) 2 Next, for the functions Λ∓ (x), in the denominators, in 6 0 1 2 3 4 5 6 equation (38), we have
WΓ 2 0 2|γ0| x Λ∓ (x)=∓ (40) 2 2 2 0 1 2 3 4 5 6 (W ∓ γ1 ∓ Δ) +4|γ0| x
WΓ0 and it is clear that |Λ− (ε1) | = |Λ− (ε2) | and |Λ+ (ε3) | = Fig. 3. The exact numerical solution for the chemical poten- |Λ+ (ε4) |. For the considered assumption of the half-filling tial, normalized to the intralayer hopping amplitude γ0,asa in each layer, and as the theoretical and numerical calcula- function of the interlayer Coulomb interaction parameter and tions show, the normal single-particle DOS is not the same for different values of the intralayer interaction parameter U. for different sublattices in the given layer with =1, 2and The inset: the effective chemical potentialμ ¯, calculated for the near the shifted neutrality points. For the next, we will same values of the intralayer on-site interaction U and given omit the layer indexes near the DOS functions notations by equation (33). (due to the relations in Eqs. (31)and(32), in Sect. 3). Eur. Phys. J. B (2017) 90: 130 Page 7 of 12
Similarly, for the sublattice B, we have the following ex- pression for the B DOS ρ ω 1 S k,ω B( )=N B( ) k 4 1 β δ ω − ε . = N ik ( ik) (41) i=1 k
The coefficients βik,inequation(41)withi =1,...,4are given by the following relations ⎧ P (3)(ε ) ⎨ ik 1 , if i =1, 2, i+1 (ε1k−ε2k) j=3,4 (εik−εjk) βik =(−1) ⎩ P (3)(ε ) ik 1 , if i =3, 4, (ε3k−ε4k) j=1,2 (εik−εjk) (42)