Nov 21, 2015 - 19 See, e.g., the review, Philip Phillips, Rev. Mod. Phys. 82,. 1719 (2010). 20 Shiping Feng, Jihong Qin, and Tianxing Ma, J. Phys.: Condens.
Charge dynamics in doped Mott insulators on a honeycomb lattice Xixiao Ma Department of Physics, Beijing Normal University, Beijing 100875, China
Yu Lan College of Physics and Electronic Engineering, Hengyang Normal University, Hengyang 421002, China
Ling Qin College of Physics and Engineering, Chengdu Normal University, Chengdu 611130, China
Shiping Feng
arXiv:1511.06822v1 [cond-mat.supr-con] 21 Nov 2015
Department of Physics, Beijing Normal University, Beijing 100875, China Within the framework of the fermion-spin theory, the charge transport in the doped Mott insulators on a honeycomb lattice is studied by taking into account the pseudogap effect. It is shown that the conductivity spectrum in the low-doped regime is separated by the pseudogap into a low-energy non-Drude peak followed by a broad midinfrared band. However, the decrease of the pseudogap with the increase of doping leads to a shift of the position of the midinfrared band towards to the low-energy non-Drude peak, and then the low-energy Drude behavior recovers in the high-doped regime. The combined results of both the doped honeycomb-lattice and square-lattice Mott insulators indicate that the two-component conductivity induced by the pseudogap is a universal feature in the doped Mott insulators. PACS numbers: 74.72.Kf, 74.25.Gz, 74.25.F-
I.
INTRODUCTION
In a Mott insulator, the spin configuration arrangement attempts to minimize the energy by aligning the spins on the lattice sites. However, in some cases, the spin configuration arrangement makes it impossible to minimize the energy on all spins on the lattice sties at the same time - a situation known as the geometrical spin frustration. The parent compounds of cuprate superconductors are a Mott insulator on a square lattice1 , where there is no geometrical spin frustration. Superconductivity is then realized when charge carriers are doped into this parent Mott insulator, where the phase diagram is characterized by two salient phenomena, the high temperature superconductivity and the occurrence of the pseudogap2 . In particular, the pseudogap is particularly obvious in the underdoped regime, where the doping concentration is too low for the optimal superconductivity. A number of experimental probes3,4 show that below a characteristic temperature T ∗ , which can be well above the superconducting (SC) transition temperature Tc in the underdoped regime, the physical response of cuprate superconductors can be interpreted in terms of the formation of a pseudogap by which it means a suppression of the spectral weight of the low-energy excitation spectrum. In this case, an important issue is whether the pseudogap phenomenon occurred in the doped squarelattice Mott insulators is universal or not. By universality we refer to the pseudogap phenomenon that does not depend on the details of the geometrical spin frustration. The doped Mott insulators on a honeycomb lattice in many ways provide an ideal arena in which to explore this question. This follows a fact that the honeycomb lattice is a loosely coupled lattice with only three nearestneighbor (NN) sites, where the strongly geometrical spin frustration and quantum fluctuation exist5–11 , and there-
fore allowing a test of the effect of the geometrical spin frustration on the pseudogap state. Although there are only a few examples of materials where the spins are located on a honeycomb lattice, many fascinating phenomena have been reported in experiments5–11 , such as the quantum spin-liquid state and superconductivity. On the other hand, the Mott insulators on a honeycomb lattice are also of interests in their own right, with many unanswered questions12–16 . In the doped square-lattice Mott insulator, the physical quantity which most evidently displays the signature for the pseudogap is the charge transport, which is manifested by the optical conductivity and resistivity, where the characteristic features are that the conductivity is divided into two-components by the pseudogap with a non-Drude-like narrow band centered around energy ω ∼ 0 followed by a broadband centered in the midinfrared region, while the resistivity grows nearly linearly with temperature3,4 . In particular, this two-component conductivity extends to the pseudogap boundary in the phase diagram at T ∗ . However, this anomalous structure of the conductivity spectrum in the pseudogap state can not be described by the standard Landau Fermi-liquid theory. In our earlier studies17 , the charge transport in the doped square-lattice Mott insulators has been discussed by considering the pseudogap effect, and then all main features of the optical measurements on the doped square-lattice Mott insulators are qualitatively reproduced3 . In this paper, we try to study the charge transport in the doped honeycomb-lattice Mott insulators along with this line. We show that the qualitative behavior of the conductivity in the doped honeycomblattice Mott insulators is similar with that obtained in the doped square-lattice Mott insulators. In particular, we show that the transfer of the part of the low-energy spectral weight in the conductivity in the low-doped regime
2 to the higher energy region to form the midinfrared band can be also attributed to the effect of the pseudogap on the infrared response in the doped honeycomb-lattice Mott insulators. The paper is organized as follows. The basic formalism is presented in Sec. II, where we generalize the theory of the charge transport from the previous case in the doped square-lattice Mott insulators to the present case in the doped honeycomb-lattice Mott insulators, and then evaluate explicitly the conductivity. Within this theoretical framework, the quantitative characteristics of the conductivity in the doped honeycomb-lattice Mott insulators are discussed in Sec. III. Finally, we give a summary in Sec. IV.
II.
FORMALISM
y
† A, B is the sublattice index, Cυlσ (Cυlσ ) is the electron y x z creation (annihilation) operator, Sυl = (Sυl , Sυl , Sυl ) are spin operators, t is the NN hopping integral, J is the NN spin-spin antiferromagnetic (AF) exchange coupling constant, and µ is the chemical potential. It has been shown that the local constraint of no double electron occupancy in the t-J model (3) can be treated properly within the fermion-spin theory20,21 , where the constrained electron operators Cυl↑ and Cυl↓ are decoupled as, − Cυl↑ = h†υl↑ Sυl ,
+ Cυl↓ = h†υl↓ Sυl ,
(4)
respectively, with the fermion operator hυlσ = e−iΦlσ hυl that keeps track of the charge degree of freedom together with some effects of spin configuration rearrangements due to the presence of the doped charge carrier itself, while the spin operator Sυl represents the spin degree of freedom, and then the local constraint of no double electron occupancy is always satisfied in actual calculations. In this fermion-spin representation (4), the original t-J model (3) can be rewritten as, H = t
X † + − − + (hBl+ˆη↑ hAl↑ SAl SBl+ˆη + h†Bl+ˆη↓ hAl↓ SAl SBl+ˆη lˆ η
+
+ − h†Al+ˆη↑ hBl↑ SBl SAl+ˆη
− µ
x FIG. 1: The honeycomb lattice with the primitive lattice vectors τ1 and τ2 . Each unit cell contains two sites (indicates by open and solid circles) belonging to the sublattices A and B, respectively.
The t-J model, or equivalently the large U Hubbard model, are prototypes to study the strong correlation effects in solids, especially in connection with the doped Mott insulators18,19 . On the other hand, the honeycomb lattice as shown in Fig. 1 is formed by the primitive lattice vectors12 , √ √ τ1 = [3/2, 3/2], τ2 = [3/2, − 3/2], (1) where the lattice constant has been set as a = 1. Since each unit cell contains two sites, belonging to different sublattices A and B, the corresponding reciprocal lattice vectors in momentum-space are given by, √ √ (2) k1 = [1/3, 1/ 3]2π, k2 = [1/3, −1/ 3]2π.
X υlσ
− + + h†Al+ˆη↓ hBl↓ SBl SAl+ˆη ) X h†υlσ hυlσ + Jeff (SAl · SBl+ˆη + SBl · SAl+ˆη ), (5) lˆ η
with Jeff = (1 − δ)2 J, and δ = hh†υlσ hυlσ i = hh†υl hυl i is the doping concentration. Within the framework of the fermion-spin theory (4), the charge transport in the doped square-lattice Mott insulators has been discussed17 , and the result shows that the striking behavior of the low-energy non-Drude peak and unusual midinfrared band is closely related to the emergence of the pseudogap. In the following discussions, we generalize these analytical calculations from the case in the doped square-lattice Mott insulators17 to the case in the doped honeycomb-lattice Mott insulators, and then evaluate explicitly the conductivity. However, as we have mentioned above, there are two sublattices A and B in the unit cell, and in this case, the chargecarrier and spin Green’s functions are two-dimensional matrices. Following our previous discussions for the case in the doped square-lattice Mott insulators17 , the corresponding intra-sublattice and inter-sublattice parts of the full charge-carrier Green’s functions can be evaluated explicitly as (see the Appendix A),
In this case, the t-J model on a honeycomb lattice can 1 1 X , (6a) gintra (k, ω) = (ν) be expressed as, 2 ν=1,2 ω − ξ − Σ(h) (k, ω) k (ν) X † X † † 1 iθk X 1 H = −t (CAlσ CBl+ˆησ + CBlσ CAl+ˆησ ) + µ Cυlσ Cυlσ g e (−1)ν . (6b) inter (k, ω) = (ν) (h) 2 lˆ ησ υlσ ω − ξk − Σ(ν) (k, ω) ν=1,2 X + J (SAl · SBl+ˆη + SBl · SAl+ˆη ), (3) where ν = 1, 2, and the mean-field (MF) charge-carrier lˆ η (ν) excitation spectrum ξk = (−1)ν Ztχ|γk | − µ, the spin P † + − + − correlation function χ = hSAl SBl+ˆη i = hSBl SAl+ˆη i, Z = supplemented by a local constraint σ Cυlσ Cυlσ ≤ 1 to 3 is the number of the NN sites on a honeycomb lattice, remove double electron occupancy, where the summation P γk = (1/Z) ηˆ e−ik·ˆη = eiθk |γk |, while the charge-carrier is over all sites l, and for each l, over its NN sites ηˆ, υ =
3 self-energies are obtained as, (h)
with the electron current-current correlation function,
(h)
(h)
Σ(1) (k, ω) = Σintra (k, ω) + e−iθk Σinter (k, ω),
(7a)
(h) Σ(2) (k, ω)
(7b)
=
(h) Σintra (k, ω)
−
(h) e−iθk Σinter (k, ω),
Π(τ − τ 0 ) = −hTτ j(τ ) · j(τ 0 )i,
with the corresponding intra-sublattice and inter(h) sublattice parts of the self-energy, Σintra (k, ω) and (h) Σinter (k, ω), repectively, that are evaluated explicitly in the Appendix A. In a doped Mott insulator, the strong correlation effect and the related charge-carrier quasiparticle coherence are closely related to the charge(h) carrier self-energy Σ(ν) (k, ω), and then the behavior of the charge-carrier quasiparticle is most fully described in terms of the charge-carrier Green’s function (6). Now we turn to calculate the conductivity of the doped honeycomb-lattice Mott insulators in terms of the full charge-carrier Green’s functions (6). The conductivity in the system is expressed as22 , σ(ω) = −
j =
ImΠ(ω) , ω
(9)
where j is the electron current operator. The external magnetic field can be coupled to the electrons, which are − + now represented by Cυl↑ = h†υl↑ Sυl and Cυl↓ = h†υl↓ Sυl in the fermion-spin representation (4). For the calculation of the electron current density, we need to obtain the electron polarization operator, which is a summation over all the particles and their positions, and can be evaluated in the fermion-spin representation as P P † P = −e Rl Cυlσ Cυlσ = e Rl h†υl hυl . The electron υlσ
υl
current operator on the other hand is obtained by evaluating the time-derivative of the polarization operator j = ∂P/∂t = (i/~)[H, P] as22 ,
(8)
iet X + − − + + − − + ηˆ(hAl↑ h†Bl+ˆη↑ SAl SBl+ˆη + hAl↓ h†Bl+ˆη↓ SAl SBl+ˆη − hBl↑ h†Al+ˆη↑ SBl SAl+ˆη − hBl↓ h†Al+ˆη↓ SBl SAl+ˆη ). (10) ~ lˆ η
In the fermion-spin approach, this electron current operator (10) can be decoupled as, j =
ieχt X ieφt X + − − + + − − + ηˆ(h†Al+ˆησ hBlσ − h†Bl+ˆησ hAlσ ) + ηˆ(SAl SBl+ˆη + SAl SBl+ˆη − SBl SAl+ˆη − SBl SAl+ˆη ), ~ ~ lˆ ησ
(11)
lˆ η
where φ = hh†Blσ hAlσ i = hh†Alσ hBlσ i is the charge-carrier’s particle-hole parameter. The second term in the right-hand side of Eq. (11) refers toPthe contribution to the electron current density from the electron spin, however, it has been + − − + + − − + shown explicitly17 that ηˆ(SAl SBl+ˆη + SAl SBl+ˆη − SBl SAl+ˆη − SBl SAl+ˆη ) ≡ 0, i.e., there is no direct contribution to lˆ η
the electron current density from the electron spin, and then the majority contribution to the electron current-current correlation comes from the charge carriers (then the electron charge). However, the strong interplay between the charge carriers and spins has been considered through the spin’s order parameters entering in the charge-carrier part of the contribution to the current-current correlation. In this case, the electron current-current correlation function can be evaluated in terms of the full charge carrier Green’s function (6) as, Π(ipm ) = 2(Ze)2
1 X (g) (g) [γintra (k)Πintra (k, ipm ) + γinter (k)Πinter (k, ipm )], N
(12)
k
where the intra-sublattice and inter-sublattice current vertices are given by, √ √ 4 (χt)2 [sin2 ( 3ky /2) − cos( 3ky /2) cos(3kx /2) + 1], (13a) 9 √ √ √ √ √ 2 γinter (k) = − (χt)2 [2 cos(2 3ky ) − 2 cos( 3ky ) + 4 cos(3kx ) cos( 3ky ) + 4 cos(3kx /2) cos( 3ky /2) cos( 3ky ) 9 √ √ √ − 6 cos(3kx /2) cos( 3ky /2) + cos(3kx ) − 7]/[2 cos( 3ky ) + 4 cos( 3ky /2) cos(3kx /2) + 3], (13b)
γintra (k) =
and the intra-sublattice charge-carrier bubble Πintra (k, ω) is obtained from the corresponding intra-sublattice part of the charge-carrier Green’s function (6a) as, (g)
Πintra (k, ipm ) =
1X gintra (k, iωn + ipm )gintra (k, iωn ), β iω
(14)
n
and is closely related to the quasiparticle scattering within the same sublattice, while the inter-sublattice charge-carrier bubble Πinter (k, ω) is obtained from the corresponding inter-sublattice part of the charge-carrier Green’s function (6b)
4 as, (g)
Πinter (k, ipm ) = e−2iθk
1X ginter (k, iωn + ipm )ginter (k, iωn ), β iω
(15)
n
and therefore is directly associated with the quasiparticle scattering between the sublattices A and B. Substituting the electron current-current correlation function (12) and the intra-sublattice and inter-sublattice charge-carrier bubbles (14) and (15) into Eq. (8), we thus obtain the conductivity as σ(ω) = σintra (ω) + σinter (ω), where the intra-sublattice conductivity σintra (ω) and inter-sublattice conductivity σinter (ω) can be evaluated explicitly as, Z (Ze)2 X ∞ dω 0 nF (ω 0 + ω) − nF (ω 0 ) (h) (h) γintra (k)Aintra (k, ω 0 + ω)Aintra (k, ω 0 ) , (16a) σintra (ω) = N ω −∞ 2π k Z (Ze)2 X −2iθk ∞ dω 0 nF (ω 0 + ω) − nF (ω 0 ) (h) (h) σinter (ω) = e γinter (k)Ainter (k, ω 0 + ω)Ainter (k, ω 0 ) , (16b) N ω −∞ 2π k
respectively, where the intra-sublattice and inter-sublattice components of the charge-carrier spectral functions are (h) (h) given by Aintra (k, ω) = −2Imgintra (k, ω) and Ainter (k, ω) = −2Imginter (k, ω), respectively, while nF (ω) is the fermion distribution function. III.
QUANTITATIVE CHARACTERISTICS
In Fig. 2, we plot the conductivity σ(ω) as a function of energy at doping levels (a) δ = 0.09 and (b) δ = 0.15 for parameter t/J = 2.5 with temperature T = 0.002J, where the charge e has been set as the unit. It is remarkable that as the case of the doped square-lattice Mott insulators17 , σ(ω) at the low-doped regime is separated into two bands, with the higher energy band, corresponding to the unusual midinfrared band, that shows a broad peak around ω ∼ 0.6t for δ = 0.09, while the low-energy band forms a sharp peak at ω ∼ 0, and deviates strongly from the Drude behavior. To show this deviation from the Drude behavior clearly, we have fitted our present result of σ(ω) for the doping concentration δ = 0.09, and the fitted result is also shown in Fig. 2a (dashed line). Obviously, the fitted result of the lower-energy peak decays as → 1/ω. This 1/ω decay of the conductivity at low energies is closely related with the linear temperature resistivity, since it reflects an anomalous frequency dependent scattering rate proportional to ω instead of ω 2 as would be expected in the standard Landau Fermiliquid. Furthermore, as seen from Fig. 2b, although the weight of the midinfrared band increases with the increase of doping, the midinfrared band moves towards to the low-energy non-Drude band. In particular, we find that the low-energy non-Drude peak incorporates with the midinfrared band in the high-doped regime, and then the low-energy Drude type behavior of the conductivity recovers. All these results are qualitatively similar to the case in the doped square-lattice Mott insulators17 , reflecting that the two-component conductivity is a universal feature for the doped Mott insulators. In Eq. (16), there are two parts of the quasiparticle contribution to the redistribution of the spectral weight in the conductivity σ(ω): the contribution from σintra (ω) comes from the spectral function obtained in terms of the intra-sublattice part of the charge-carrier Green’s function (6a), and therefore is closely associated with the quasiparticle scattering within the same sub-
lattice, while the additional contribution from σinter (ω) originates from the spectral function obtained in terms of the inter-sublattice part of the charge-carrier Green’s function (6b), and is closely related to the quasiparticle scattering between the sublattices A and B. Within the framework of the fermion-spin theory, the essential physics of the charge transport in the doped honeycomblattice Mott insulators is the same as that in the doped square-lattice case17 , and can be attributed to the emergence of the pseudogap. This follows a fact that there is an intrinsic link between the self-energy and pseudogap. (h) The charge-carrier self-energy Σ(ν) (k, ω) in Eq. (7) also can be rewritten as21,23 , (ν)
(h)
Σ(ν) (k, ω) ≈
¯ pg (k)]2 [∆ (ν)
,
(17)
ω + εk
(ν)
(h)
where εk is the energy spectrum of Σ(ν) (k, ω), while ¯ (ν) ∆ pg (k) is the pseudogap, and is therefore identified as being a region of the self-energy effect in which the pseudogap suppresses the spectral weight. This pseu¯ (ν) dogap ∆ pg (k) and the corresponding energy spectrum (ν) εk can be obtained explicitly in terms of the q self-energy (h)
(ν)
(ν)
¯ pg (k) = L (k)/ Σ(ν) (k, ω) in Eq. (7) as ∆ 2 (ν)
(ν)
(ν)
(ν)
L1 (k) and
εk = L2 (k)/L1 (k), respectively, with the functions (h) (ν) (h) (ν) L1 (k) = −Σ(ν)o (k, ω = 0) and L2 (k) = Σ(ν)e (k, ω = 0). In this case, the pseudogap parameter is obtained as P ¯ (ν) ¯ (ν) ∆ pg = (1/N ) k ∆pg (k). In Fig. 3, we plot the pseudo¯ (1) ¯ (2) gap parameters ∆ pg (solid line) and ∆pg (dashed line) as a function of doping for t/J = 2.5 with T = 0.002J. It is ¯ (1) ¯ (2) shown clearly that although both ∆ pg and ∆pg decrease monotonously with the increase of doping, the decrease ¯ (1) of the magnitude with increasing doping is faster in ∆ pg (2) (1) (2) ¯ pg . In particular, ∆ ¯ pg ≈ ∆ ¯ pg near the halfthan in ∆ ¯ (1) ¯ (2) filling. Since both the pseudogaps ∆ pg and ∆pg appear in the intra-sublattice and inter-sublattice Green’s functions in Eq. (6), they together induce the transfer of the
5
I\ 4
(a) 4 ..--圄、
""..圄、
4ω 4
呻ω d
c
c
2
2
.
DLE
m
.
DLE
2
2
。
、-
(b)
、-田.,. .,-.、
."......、
e
s
、』国,
、_.".
。
b
。。
0.2
0.4
0.6
0.8
1.0
00
0.2
0.4
0.6
1.0
0.8
FIG. 2: The conductivity as a function of energy at (a) δ = 0.09 and (b) δ = 0.15 with T = 0.002J for t/J = 2.5. The dashed line at (a) is obtained from a numerical fit σ(ω) = A/(ω + B), with A ∼ 0.058 and B ∼ −0.017.
part of the low-energy spectral weight in the conductivity in the low-doped regime to the higher energy region to form the midinfrared band. Moreover, the onset of the region to which the spectral weight is transferred is ¯ pg ∼ ∆ ¯ (1) ¯ (2) close to the sum of the pseudogaps ∆ pg + ∆pg , and then the decrease of the pseudogaps with increasing doping therefore leads to a shift of the position of the midinfrared band towards to the low-energy non-Drude peak.
、
、
、
、
、
、
、
「,\。巳叫
0.6 0.4 0.2 0
0
0.05
0.10
0.15
Acknowledgments
The authors would like to thank L¨ ulin Kuang and Huaisong Zhao for helpful discussions. This work was supported by the funds from the Ministry of Science and Technology of China under Grant No. 2012CB821403, and the National Natural Science Foundation of China under Grant Nos. 11274044 and 11574032, and YL was supported by the Science Foundation of Hengyang Normal University under Grant No. 13B44, and Hunan Provincial Natural Science Foundation of China under Grant No. 2015JJ3027.
、
、
、
、
、
0.8
、
、
、
、
、
、
、
、
、
1.0
type behavior of the conductivity recovers in the highdoped regime. Incorporating the present result with that obtained in the doped square-lattice Mott insulators17 , it is thus shown that the two-component conductivity induced by the pseudogap is a universal feature in the doped Mott insulators.
0.20 Appendix A: Charge-carrier and spin Green’s functions in the t-J model on a honeycomb lattice
¯ (1) FIG. 3: The pseudogap parameters (∆ pg ) (solid line) and (2) ¯ (∆pg ) (dashed line) as a function of doping with T = 0.002J for t/J = 2.5.
We first define the charge-carrier and spin Green’s functions as, g(k, ω) =
IV.
CONCLUSIONS
D(k, ω) =
Within the framework of the fermion-spin theory, we have studied the charge transport of the doped Mott insulators on a honeycomb lattice by taking into account the pseudogap effect. Our result shows that the conductivity spectrum in the low-doped regime is separated into a non-Drude peak centered around energy ω ∼ 0 followed by a broad midinfrared band. However, the decrease of the pseudogaps with increasing doping leads to a shift of the position of the midinfrared band towards to the lowenergy non-Drude peak, and then the low-energy Drude
D(z) (k, ω) =
gintra (k, ω) ginter (k, ω) ∗ ginter (k, ω) gintra (k, ω)
Dintra (k, ω) Dinter (k, ω) ∗ Dinter (k, ω) Dintra (k, ω)
,
(A1a) ,
D(z)intra (k, ω) D(z)inter (k, ω) ∗ D(z)inter (k, ω) D(z)intra (k, ω)
(A1b) . (A1c)
Since the spin system in the t-J model (5) is anisotropic away from half-filling, therefore two spin Green’s functions D(k, ω) and Dz (k, ω) have been defined to describe properly spin propagations21 . In the MF approximation, the t-J model (5) can be decoupled explicitly as HMFA = Ht + HJ − 2N Ztφχ, with
6
Ht = χt
X † X † (hBl+ˆησ hAlσ + h†Al+ˆησ hBlσ ) − µ hυlσ hυlσ , lˆ ησ
HJ =
(A2a)
υlσ
X 1 + − − + + − − + z z z z Jeff [(SAl SBl+ˆη + SAl SBl+ˆη + SBl SAl+ˆη + SBl SAl+ˆη ) + 2(SAl SBl+ˆ η + SBl SAl+ˆ η )], 2
(A2b)
lˆ η
where = 1 + 2tφ/Jeff . Based on the equation of motion method, it is easy to find the intra-sublattice and intersublattice parts of the MF charge-carrier Green’s function, 1 1 X (0) , (A3a) gintra (k, ω) = 2 ν=1,2 ω − ξ (ν) k X 1 1 (0) ginter (k, ω) = eiθk , (A3b) (−1)ν (ν) 2 ω − ξk ν=1,2 and the intra-sublattice and inter-sublattice parts of the MF spin Green’s functions, ! X B (ν) 1 1 (0) k Dintra (k, ω) = − , (ν) (ν) (ν) ω − ωk ω + ωk ν=1,2 2ωk (ν)
(0) Dinter (k, ω)
= e
iθk
X
(−1)
ν
X B (ν) zk
(0)
(0) D(z)inter (k, ω)
= e
(ν)
ω−
X B (ν) zk (ν)
ν=1,2
(ν)
(ν)
ω − ωk
1
(ν) ν=1,2 2ωzk iθk
1
2ωk
ν=1,2
D(z)intra (k, ω) =
Bk
2ωzk
(ν) ωzk
!
1
−
(ν)
,
(A4b)
ω + ωk !
1
−
(A4a)
(ν)
,
(A4c)
ω + ωzk
1 (ν)
ω − ωzk
−
!
1 (ν)
,
(A4d)
ω + ωzk
(ν)
with Bk = λ[(−1)ν A1 |γk | − A2 ], Bzk = −λχ[1 − (−1)ν |γk |]/2, λ = 2ZJeff , A1 = χz + χ/2, A2 = χz + χ/2, the z z spin correlation function χz = hSAl SBl+ˆ η i, and the MF spin excitation spectra, 1 1 1 A3 ) + 2 A3 + A4 ] − (−1)ν λ2 [αA2 (1 − ) + + A4 ]|γk |, (A5a) Z 2 Z 2 1 1 (ν) (ωzk )2 = λ2 [αχ(|γ|2k − ) + A3 ] − (−1)ν λ2 [αχ(1 − ) + A3 ]|γk |, (A5b) Z Z where P A3 = αC + (1 − α)/(2Z), A4P= αCz + (1 − α)/(4Z), and the spin correlation functions C = + − z 2 z z (1/Z 2 ) ηˆηˆ0 hSAl+ˆ η SAl+ˆ η 0 i. In order to satisfy the sum rule of the correlation ηˆηˆ0 hSAl+ˆ η SAl+ˆ η 0 i, C = (1/Z ) + − function hSυl Sυl i = 1/2 in the case without an AF long-range order, the important decoupling parameter α has been introduced in the above calculation17,21 , which can be regarded as the vertex correction. As in the case of the doped square-lattice Mott insulators17 , the self-consistent equation that is satisfied by the full charge-carrier Green’s function is obtained as22 , (ν)
(ωk )2 = λ2 [αA1 (|γ|2k −
g(k, ω) = g (0) (k, ω) + g (0) (k, ω)Σ(h) (k, ω)g(k, ω),
(A6)
where the intra-sublattice and inter-sublattice parts of the self-energy Σ(h) (k, ω) is obtained from the spin bubble as, 1 X 1X (h) Σintra (k, iωn ) = |Λp+q+k |2 gintra (k + p, iωn + ipm )Πintra (p, q, ipm ), (A7a) 2 N pq β ip m
(h) Σinter (k, iωn )
1 X 2 1 X ∗ = g (k + p, iωn + ipm )Π∗inter (p, q, ipm ), Λp+q+k 2 N pq β ip ,iq inter m
(A7b)
n
with Λk = Ztγk , and the corresponding intra-sublatticel and inter-sublattice parts of the spin bubble, 1 X (0) (0) Πintra (p, q, ipm ) = Dintra (q, iqm )Dintra (p + q, ipm + iqm ), β iq
(A8a)
m
Πinter (p, q, ipm ) =
1 X (0) (0) Dinter (q, iqm )Dinter (p + q, ipm + iqm ). β iq
(A8b)
m
Since the self-energy Σ(h) (k, ω) is not an even function of ω, we follow common practice to separate
Σ(h) (k, ω) into the symmetric and antisymmetric parts,
7 (h)
(h)
i.e., Σ(h) (k, ω) = Σe (k, ω) + ωΣo (k, ω), and then (h) (h) both Σe (k, ω) and Σo (k, ω) are an even function (h) of ω. Moreover, Σo (k, ω) is directly associated with the charge-carrier quasiparticle coherent weight22 as (h) −1 ZhF (k, ω) = 1 − ReΣo (k, ω). In the static limit, although ZhF (k) still is a function of momentum, however, the wave vector k in ZhF (k) can be chosen as k0 = [2π/3, 0], and then the corresponding intra-sublattice and inter-sublattice parts of the quasiparticle coherent weight can be obtained as, (ν)−1 ZhF
=1−
(h) ReΣ(ν)o (k, ω
= 0) |k0 .
(A6) can be evaluated explicitly as,
gintra (k, ω) =
1 = 8N 2
X
|Λp+q+k |
pqνν 0 ν 00
2
ginter (k, ω) =
(3)
Fpqk
+ ω+ (h) Σinter (k, ω)
(ν 00 ) ωq+p
1 = 8N 2
+
X
(ν 0 ) ωq
(−1)
(ν) − ξ¯p+k
+
(ν)
ω−
(1)
(ν 00 )
(ν 0 )
+ ωq
e
(ν) − ξ¯p+k
(ν 0 )
ω − ωq+p − ωq
−
(ν) ξ¯p+k
+
ω + ωq+p +
(ν 0 ) ωq
−
(ν) ξ¯p+k
0
0
(ν 00 )
(ν 0 )
ω − ωq+p − ωq
(ν) − ξ¯p+k
(A11a)
(ν 0 ) (ν 00 ) (ν) Bq Bq+p 2 Λp+q+k ZhF (ν 0 ) (ν 00 ) ωq ωq+p
Fpqk (ν 00 )
Fpqk
+
,
(ν) − ξ¯p+k
(3)
Fpqk
00
(ν 0 )
ω + ωq+p − ωq ! (4)
pqνν 0 ν 00
(ν 00 )
(2)
Fpqk
ν+ν 0 +ν 00 −i(θq+p +θq +θp+k )
(2)
+
(ν) (ν)
where ξ¯k = ZhF ξk . With the help of the Green’s functions in Eq. (A10) and the spin Green’s functions in Eq. (A4), the intra-sublattice and inter-sublattice parts of the self-energy in Eq. (A7) can be evaluated explicitly as,
Fpqk (ν 00 ) ωq+p
(ν) ZhF 1 iθk X , (A10b) e (−1)ν (ν) 2 ω − ξ¯ ν=1,2 k
(A9)
(ν 0 ) (ν 00 ) (ν) Bq Bq+p ZhF (ν 0 ) (ν 00 ) ωq ωq+p
(A10a)
k
In this case, the intra-sublattice and inter-sublattice parts of the full charge-carrier Green’s functions in Eq.
(h) Σintra (k, ω)
(ν) 1 X ZhF , 2 ν=1,2 ω − ξ¯(ν)
+
(1)
Fpqk (ν 00 )
Fpqk ω−
(ν 00 ) ωq+p
(ν 0 )
ω + ωq+p − ωq ! (4) (ν 0 )
+ ωq
00
(ν)
− ξ¯p+k
(ν) − ξ¯p+k
,
(A11b)
0
00
(ν) (ν ) (ν ) (2) (ν ) (ν ) (1) (ν ) (ν ) where Fpqk = nB (ωq+p )[1 + nB (ωq )] + nF (ξ¯p+k )[nB (ωq ) − nB (ωq+p )], Fpqk = [1 + nB (ωq )][1 + nB (ωq+p )] − (4) (ν 0 ) (ν) (ν 0 ) (ν 00 ) (3) (ν 00 ) (ν) (ν 0 ) (ν 00 ) )], F = )nB (ωq ) + nF (ξ¯ )[1 + nB (ωq ) + nB (ω )], F = nB (ω nF (ξ¯ )[1 + nB (ωq ) + nB (ω p+k (ν 0 ) nB (ωq )[1
q+p
+
(ν 00 ) nB (ωq+p )]
q+p
pqk
0
q+p
p+k
In this case, the self-consistent equations that are satisfied by the intra-sublattice and inter-sublattice parts of the quasiparticle coherent weight in Eq. (A9) must be solved simultaneously with following self-consistent equations,
C =
(ν) Bk 1 X 1 (ν) |γk |2 (ν) coth[ βωk ], 2N 2 ω ν,k
χ
z
(ν) Bzk 1 X 1 (ν) ν = (−1) |γk | (ν) coth[ βωzk ], 2N 2 ω
Cz =
1 4N
φ =
1 4N
ν,k
(ν) 1 1 X Bk 1 (ν) = coth[ βωk ], (ν) 2 2N 2 ν,k ωk
χ =
(ν) Bk 1 X 1 (ν) (−1)ν |γk | (ν) coth[ βωk ], 2N 2 ω ν,k
k
(A12c)
(A12g)
zk
then all the order parameters, the decoupling parameter α, and the chemical potential µ are determined by the self-consistent calculation without using any adjustable parameters. With the above self-consistent calculation of the order parameters, the intra-sublattice and intersublattice parts of the full charge-carrier Green’s functions can be evaluated explicitly as, gintra (k, ω) =
1X 1 , (ν) 2 ν ω − ξ − Σ(h) (k, ω) k
(A12d)
(A12f)
zk
(ν) Bzk 1 (ν) 1 X |γk |2 (ν) coth[ βωzk ], 2N 2 ω ν,k
X (ν) 1 (ν) ZhF 1 − tanh[ β ξ¯k ] , (A12a) 2 ν,k X 1 (ν) (ν) (−1)ν |γk |ZhF 1 − tanh[ β ξ¯k ] ,(A12b) 2
(A12e)
k
ν,k
δ =
pqk
00
(ν) (ν ) (ν ) − nF (ξ¯p+k )[nB (ωq ) − nB (ωq+p )], and nB (ω) is the boson distribution function.
ginter (k, ω) =
(A13a)
(ν)
1 iθk X 1 e (−1)ν , (A13b) (ν) (h) 2 ω − ξ − Σ (k, ω) ν k
(ν)
8 as quoted in Eq. (6).
1 2
3
4
5
6
7
8
9
10
11
J. G. Bednorz and K. A. M¨ uller, Z. Phys. B 64, 189 (1986). See, e.g., the review, S. H¨ ufner, M. A. Hossain, A. Damascelli, and G. A. Sawatzky, Rep. Prog. Phys. 71, 062501 (2008). See, e.g., the review, Tom Timusk and Bryan Statt, Rep. Prog. Phys. 62, 61 (1999). B. Batlogg, H. Y. Hwang, H. Takagi, R. J. Cava, H. L. Kao, and J. Kwo, Physica C 235-240, 130 (1994). S. Shamoto, T. Kato, Y. Ono, Y. Miyazaki, K. Ohoyama, M. Ohashi, Y. Yamaguchi, and T. Kajitani, Physica C 306, 7 (1998). R. Weht, A. Filippetti, and W. E. Pickett, Europhys. Lett. 48, 320 (1999). V. Kataev, A. M¨ oller, U. L¨ ow, W. Jung, N. Schittner, M. Kriener, and A. Freimuth, J. Magn. Magn. Mater. 290, 310 (2005). A. M¨ oller, U. L¨ ow, T. Taetz, M. Kriener, G. Andr´e, F. Damay, O. Heyer, M. Braden, and J. A. Mydosh, Phys. Rev. B 78, 024420 (2008). Soichiro Okumura, Hikaru Kawamura, Tsuyoshi Okubo, and Yukitoshi Motome, J. Phys. Soc. Jpn. 79, 114705 (2010). Y. J. Yan, Z. Y. Li, T. Zhang, X. G. Luo, G. J. Ye, Z. J. Xiang, P. Cheng, L. J. Zou, and X. H. Chen, Phys. Rev. B 85, 085102 (2012). D. E. McNally, J. W. Simonson, J. J. Kistner-Morris, G. J. Smith, J. E. Hassinger, L. DeBeer-Schmitt, A. I. Kolesnikov, I. A. Zaliznyak, and M. C. Aronson, Phys. Rev. B 91, 180407(R) (2015).
12
13
14
15
16
17
18 19
20
21
22
23
Sandro Sorella, Yuichi Otsuka, and Seiji Yunoki, Scientific Reports 2, 992 (2012). James M. Murray and Oskar Vafek, Phys. Rev. B 89, 205119 (2014). Annica M. Black-Schaffer and Carsten Honerkamp, J. Phys.: Condens. Matter 26, 423201 (2014). Zhenyue Zhu and Steven R. White, Mod. Phys. Lett. B 28, 1430016 (2014). Shenghan Jiang, Andrej Mesaros, and Ying Ran, Phys. Rev. X 4, 031040 (2014). Ling Qin, Jihong Qin, and Shiping Feng, Physica C 497, 77 (2014); Feng Yuan, Jihong Qin, Shiping Feng, and Wei Yeu Chen, Phys. Rev. B 67, 134504 (2003); Shiping Feng and Zhongbing Huang, Phys. Lett. A 232, 293 (1997). P. W. Anderson, Science 235, 1196 (1987). See, e.g., the review, Philip Phillips, Rev. Mod. Phys. 82, 1719 (2010). Shiping Feng, Jihong Qin, and Tianxing Ma, J. Phys.: Condens. Matter 16, 343 (2004); Shiping Feng, Z.B. Su, and L. Yu, Phys. Rev. B 49, 2368 (1994). See, e.g., the review, Shiping Feng, Yu Lan, Huaisong Zhao, L¨ ulin Kuang, Ling Qin, and Xixiao Ma, Int. J. Mod. Phys. B 29, 1530009 (2015). See, e.g., G. D. Mahan, Many-Particle Physics, (Plenum Press, New York, 1981). Shiping Feng, Huaisong Zhao, and Zheyu Huang, Phys. Rev. B. 85, 054509 (2012); Yu Lan, Jihong Qin, and Shiping Feng, Phys. Lett. A 377, 2210 (2013).