Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model
The effect of correlated hopping on the charge and heat transport is investigated for the Falicov-Kimball model. An exact solutions for the electrical and thermal conductivities and thermoelectric power are obtained within the dynamical mean field theory. The temperature dependences of the transport...
Збережено в:
Дата: | 2014 |
---|---|
Автор: | |
Формат: | Стаття |
Мова: | English |
Опубліковано: |
Інститут фізики конденсованих систем НАН України
2014
|
Назва видання: | Condensed Matter Physics |
Онлайн доступ: | http://dspace.nbuv.gov.ua/handle/123456789/152911 |
Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Цитувати: | Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model / A.M. Shvaika // Condensed Matter Physics. — 2014. — Т. 17, № 4. — С. 43704: 1–14. — Бібліогр.: 33 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-152911 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1529112019-06-14T01:26:39Z Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model Shvaika, A.M. The effect of correlated hopping on the charge and heat transport is investigated for the Falicov-Kimball model. An exact solutions for the electrical and thermal conductivities and thermoelectric power are obtained within the dynamical mean field theory. The temperature dependences of the transport coefficients are analysed for the particular values of correlated hopping which correspond to the significant reconstruction of the density of states and transport function. The cases with strong enhancement of thermoelectric properties are elucidated. Дослiджено вплив корельованого переносу на перенос заряду i тепла в моделi Фалiкова-Кiмбала. В рамках теорiї динамiчного середнього поля отримано точнi розв’язки для електро- та теплопровiдностi i термо-е.р.с. Дослiджено температурнi залежностi коефiцiєнтiв переносу для певних значень корельованого переносу, що вiдповiдають суттєвiй перебудовi густини станiв та транспортної функцiї. Встановлено випадки, якi вiдповiдають значному покращенню термоелектричних властивостей. 2014 Article Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model / A.M. Shvaika // Condensed Matter Physics. — 2014. — Т. 17, № 4. — С. 43704: 1–14. — Бібліогр.: 33 назв. — англ. 1607-324X PACS: 72.15.Jf, 72.20.Pa, 71.27.+a, 71.10.Fd DOI:10.5488/CMP.17.43704 arXiv:1409.8495 http://dspace.nbuv.gov.ua/handle/123456789/152911 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
English |
description |
The effect of correlated hopping on the charge and heat transport is investigated for the Falicov-Kimball model. An exact solutions for the electrical and thermal conductivities and thermoelectric power are obtained within the dynamical mean field theory. The temperature dependences of the transport coefficients are analysed for the particular values of correlated hopping which correspond to the significant reconstruction of the density of states and transport function. The cases with strong enhancement of thermoelectric properties are elucidated. |
format |
Article |
author |
Shvaika, A.M. |
spellingShingle |
Shvaika, A.M. Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model Condensed Matter Physics |
author_facet |
Shvaika, A.M. |
author_sort |
Shvaika, A.M. |
title |
Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model |
title_short |
Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model |
title_full |
Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model |
title_fullStr |
Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model |
title_full_unstemmed |
Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model |
title_sort |
effect of correlated hopping on thermoelectric properties: exact solutions for the falicov-kimball model |
publisher |
Інститут фізики конденсованих систем НАН України |
publishDate |
2014 |
url |
http://dspace.nbuv.gov.ua/handle/123456789/152911 |
citation_txt |
Effect of correlated hopping on thermoelectric properties: Exact solutions for the Falicov-Kimball model / A.M. Shvaika // Condensed Matter Physics. — 2014. — Т. 17, № 4. — С. 43704: 1–14. — Бібліогр.: 33 назв. — англ. |
series |
Condensed Matter Physics |
work_keys_str_mv |
AT shvaikaam effectofcorrelatedhoppingonthermoelectricpropertiesexactsolutionsforthefalicovkimballmodel |
first_indexed |
2025-07-14T04:22:32Z |
last_indexed |
2025-07-14T04:22:32Z |
_version_ |
1837594787604070400 |
fulltext |
Condensed Matter Physics, 2014, Vol. 17, No 4, 43704: 1–14
DOI: 10.5488/CMP.17.43704
http://www.icmp.lviv.ua/journal
Effect of correlated hopping on thermoelectric
properties: Exact solutions for the Falicov-Kimball
model
A.M. Shvaika
Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine,
1 Svientsitskii St., 79011 Lviv, Ukraine
Received September 30, 2014, in final form November 18, 2014
The effect of correlated hopping on the charge and heat transport is investigated for the Falicov-Kimball model.
Exact solutions for the electrical and thermal conductivities and thermoelectric power are obtained within the
dynamical mean field theory. The temperature dependences of the transport coefficients are analysed for par-
ticular values of correlated hopping which correspond to the significant reconstruction of the density of states
and transport function. The cases with strong enhancement of thermoelectric properties are elucidated.
Key words: thermoelectric power, correlated hopping, Falicov-Kimbal model
PACS: 72.15.Jf, 72.20.Pa, 71.27.+a, 71.10.Fd
1. Introduction
Direct transformation of the heat flow into the electric current and vice versa attracts much attention
in science and engineering [1], but its applications are limited by the low thermoelectric figure of merit
of the traditional bulk materials—metals and semiconductors— good electric conductors are also good
heat conductors. On the other hand, bad metals or compounds with strong electron correlations, possess-
ing large variability of their band spectrum and density of states depending on the chemical structure
and doping, are considered as possible new candidates for thermoelectric materials [1, 2]. In many cases,
the anomalous properties of strongly correlated materials are ascribed to the on-site Coulomb or spin
interactions, which are the leading contributions.
It was already noticed by Hubbard in his seminal article [3] that besides the local Coulomb-type in-
teraction U
∑
i n̂i↑n̂i↓ there should be other nonlocal contributions: the intersite Coulomb interaction∑
i j Vi j n̂i n̂ j and the so-called correlated hopping∑
i jσ
t (2)
i j (n̂i σ̄+ n̂ j σ̄)c†
iσc jσ and
∑
i jσ
t (3)
i j n̂i σ̄c†
iσc jσn̂ j σ̄, (1.1)
which reflects the fact that different many-body states can overlap in different extent and, as a result, the
value of intersite hopping depends on the occupation of these states. The origin of the correlated hopping
can be either a direct intersite interaction or an indirect effective one [4, 5]. The effect of local Coulomb
interaction is the subject of the famous Hubbard model which is widely investigated in the theory of
strongly correlated electron systems.
The correlated hopping is less popular. Even the term used is not well established. Apart from the
term “correlated hopping”, many other terms circulate, such as “assisted hopping”, “bond-charge interac-
tion (repulsion)”, “occupation-dependent hopping”, “correlated hybridization”, etc. Similar contributions
in the theory of disordered systems are also known as an “off-diagonal disorder” [6]. Correlated hopping
was considered in connection with the new mechanisms of high temperature superconductivity [7, 8],
© A.M. Shvaika, 2014 43704-1
http://dx.doi.org/10.5488/CMP.17.43704
http://www.icmp.lviv.ua/journal
A.M. Shvaika
electron-hole asymmetry [9], and enhancement of magnetic properties [10]. Recently, the correlated hop-
ping has been examined in relation to the quantum dots [11–13] and optical lattices [14–16].
In this article we shall consider the effect of correlated hopping on the charge and heat transport
for the Falicov-Kimball model [17], which possesses exact solutions in the dynamical mean field the-
ory (DMFT) [18–21]. In section 2, we recall the main arguments of the linear response Kubo theory for
the charge and heat transport which will be important in further considerations. Section 3 provides the
DMFT solutions for the Falicov-Kimball model with correlated hopping and derivations of the charge and
energy current operators and transport coefficients in the homogeneous phase. It should be noted that at
low temperatures, the phase transitions to the ordered phases could occur, e.g., the ground state phase
diagrams for the one-dimensional (D = 1) and two-dimensional (D = 2) models display a variety of the
modulated phases [22–25], but phase diagrams for theD =∞ Falicov-Kimball model with correlated hop-
ping are unknown so far and we assume that a homogeneous solution is valid down to T = 0. In section 4
we consider peculiarities of the charge and heat transport with the change of correlated hopping value
and doping and we summarize in section 5.
2. Macroscopic and microscopic levels in the description of thermoelec-
tric effect
In the linear response Kubo theory, the charge j c (r ) and energy (heat) j Q (r ) currents are caused
by the electro-chemical potential gradient, including electrical field and charge distribution inhomo-
geneities, and temperature gradient and can be obtained from the following equations [26]
j c (r ) =−e
∫
L11(r ,r ′)∇µ̃(r ′)dr ′−e
∫
L12(r ,r ′)
∇T (r ′)
T (r ′)
dr ′, (2.1)
j Q (r ) =−
∫
L21(r ,r ′)∇µ̃(r ′)dr ′−
∫
L22(r ,r ′)
∇T (r ′)
T (r ′)
dr ′, (2.2)
where µ̃(r ) =µ−eV (r ) is an electro-chemical potential.
For the uniform (steady) dc charge and energy currents, we can define the dc electric conductivity
σdc = e2L11 , (2.3)
Seebeck coefficient (thermoelectric power E = S∇T )
S = 1
eT
L−1
11 L12 , (2.4)
and electronic contribution in thermal conductivity
κe = 1
T
[
L22 −L21L−1
11 L12
]
(2.5)
in terms of the generalized transport integrals Llm . The efficiency of the thermoelectric material is char-
acterized by the dimensionless figure of merit
Z T = T
σdcS2
κe+κph
. (2.6)
The standard route to introduce macroscopic currents in the microscopic lattice models is as fol-
lows [26]. First of all, one can define the charge polarization by
P̂ =∑
iα
R i zαn̂iα, (2.7)
where R i is the lattice site vector and n̂iα is the particle number operator at site i for particles of the kind
α with charge zα. Now, the charge current can be defined using the continuity equation as follows:
ĵ = dP̂
dt
= 1
i
[
P̂ , Ĥ
]
. (2.8)
43704-2
Effect of correlated hopping on thermoelectric properties
In a similar way, one can define the energy polarization by
P̂Q =∑
i
R i Ĥi , (2.9)
where Ĥi are the single site (local) contributions in total energy Ĥ =∑
i Ĥi , and the energy current by
ĵ Q = dP̂Q
dt
= 1
i
[
P̂Q , Ĥ
]
. (2.10)
Let us recall how it works for the case of noninteracting electrons. TheHamiltonian for noninteracting
electrons on the lattice can be written as follows:
Ĥ =∑
i
Ĥi , Ĥi = 1
2
∑
jσ
(
ti j c†
iσc jσ+ t j i c†
jσciσ
)
−µ∑
σ
c†
iσciσ , (2.11)
where Ĥi is a “local” contribution for site i . Notice that the intersite hopping term contributes to different
sites with a half weight.
The electrical charge polarization and current operators are defined by
P̂ = e
∑
iσ
R i c†
iσciσ , ĵ =−ie
∑
i jσ
(
R i −R j
)
ti j c†
iσc jσ (2.12)
and the energy polarization and current operators now take the form
P̂Q = 1
2
∑
i jσ
(
R i +R j
)
ti j c†
iσc jσ , ĵ Q =−i
∑
i jσ
(
R i −R j
)(∑
l
ti l tl j
)
c†
iσc jσ . (2.13)
On the other hand, the energy current operator can be rewritten as follows:
ĵ Q = ∑
i jσ
(
R i −R j
)
ti j
(
c†
iσ
dc jσ
dt
− dc†
iσ
dt
c jσ
)
or ĵ Q ∼ lim
t ′→t
1
2
(
d
dt
− d
dt ′
)
ĵ (t , t ′), (2.14)
where we have introduced
ĵ (t , t ′) =−ie
∑
i jσ
(
R i −R j
)
ti j c†
iσ(t )c jσ(t ′). (2.15)
Relation (2.14) is very important. In the cases when it holds, one can immediately write down an ex-
pressions for the generalized transport integrals Llm [27], the so-called Boltzmann relations, also known
as the Johnson-Mahan theorem in the theory of metals and semiconductors [28, 29]. In the simplest form,
when the dynamical screening effects can be ignored, the Boltzmann relations state that the generalized
transport integrals for electrons L11, L12 = L21, and L22 can be written in terms of one transport function
(relaxation time) I(ω):
Ll m = σ0
e2
+∞∫
−∞
dω
[
−d f (ω)
dω
]
I(ω)ωl+m−2, (2.16)
where f (ω) = 1/(eβω+1) is the Fermi distribution function. Once the screening effects and inelastic scat-
tering being taken into account, one have to replace these relations by the generalized one [27].
For the case of noninteracting electrons, the transport function has the form
Iαβ(ω) = 1
N
∑
kσ
∂εk
∂kα
∂εk
∂kβ
δ(ω+µ−εk ), (2.17)
where εk is the band energy (Fourier transform of the hopping integral ti j ), and we get the known rela-
tions
Lαβlm = σ0
e2
1
N
∑
kσ
∂εk
∂kα
∂εk
∂kβ
[
−d f (εk −µ)
dεk
]
(εk −µ)l+m−2. (2.18)
43704-3
A.M. Shvaika
The main consequences from equations (2.16) are as follows. The values of the transport integrals
Llm are determined by the features of the transport function I(ω) only within the Fermi window (n = 0)
of width ∼ 4T and its “moments”
(ω−µ)n
[
−d f (ω−µ)
dω
]
= (ω−µ)n
4T cosh2 ω−µ
2T
= T n−1
4
xn
cosh2 x/2
, x = ω−µ
T
(2.19)
which spread over a larger energy interval (see figure 1). For the metals, the transport function I(ω) is
-10 -5 0 5 10
x = (ω − µ)/Τ
-0.2
0
0.2
0.4
x
n /(
4
c
o
sh
x
/2
)2
n = 0
n = 1
n = 2
Figure 1. (Color online) Dimensionless Fermi window (n = 0) and its “moments” with extrema at points
x0 = 0, x1 =±1.543404638, and x2 =±2.399357280, respectively.
almost constant within the Fermi window which results in a weak temperature dependence of electrical
dc conductivity, the absence of thermoelectric power, and the linear temperature dependence of ther-
mal conductivity in agreement with the Wiedemann-Franz law. On the other hand, it is obvious from
equation (2.6) that large values of the figure of merit could be obtained when the thermal conductiv-
ity is strongly reduced, and it was shown by Mahan and Sofo [30] that these could be achieved for the
δ-function like transport function I (ω) ∼ δ(ω−ε0), when the electronic contribution in thermal conduc-
tivity (2.5) vanishes.
In order to get high values of the figure of merit Z T (2.6) for thermoelectric properties in some tem-
perature interval, one has to look for the systems with strongly asymmetric, within the Fermi window,
transport function I(ω)with sharp peak at (ω−µ)/T = 1.543404638 or −1.543404638 to enhance the See-
beck coefficient andminima at (ω−µ)/T =±2.399357280 to reduce the thermal conductivity. That is why,
the materials with strong electron correlations as well as disordered systems attract much attention as
prospective thermoelectric applications. This is due to the strong variability of their band structure and
density of states in the Fermi window depending on the chemical structure and doping. In particular,
it was already shown that the doping of Mott insulator can strongly improve its thermoelectric proper-
ties [31, 32]. In this case, the doping shifts the chemical potential into the lower or upper Hubbard band
which produces a strongly asymmetric density of states (d.o.s.) and transport function. Besides the Hub-
bard type local Coulomb interaction, one can also consider the non-local contributionswhich can produce
additional enhancements.
3. Charge and energy current operators and transport coefficients for
the Falicov-Kimball model with correlated hopping
A simplest model in the theory of strongly correlated electron systems is the Falicov-Kimball one [17],
which considers the local interaction between the itinerant d electrons and localized f electrons. It is
a binary alloy type model and it has an exact solution in the dynamical mean field theory (DMFT) [18].
43704-4
Effect of correlated hopping on thermoelectric properties
Its enhancement by correlated hopping was also considered, and the DMFT solutions with a nonlocal
self-energy were obtained [19–21].
The Hamiltonian of the Falicov-Kimball model with correlated hopping contains two contributions:
the local correlations Hloc and the nearest-neighbour intersite hopping Ht on the D-dimensional hyper-
cubic lattice
H = Hloc+Ht , Hloc =
∑
i
[
Uni d ni f −µ f ni f −µd ni d
]
,
Ht = 1
2
p
D
∑
〈i j 〉
[
t1d †
i d j + t2d †
i d j
(
ni f +n j f
)+ t3d †
i d j ni f n j f
]
. (3.1)
It is convenient to project on the occupation of the local f -states by introducing operators P+
i = ni f ,
P−
i = 1−ni f and to rewrite the Hamiltonian in matrix notations
Ht = 1
2
p
D
∑
〈i j 〉
[
t++i j P+
i d †
i d j P+
j + t−−i j P−
i d †
i d j P−
j + t+−i j P+
i d †
i d j P−
j + t−+i j P−
i d †
i d j P+
j
]
= 1
2
p
D
∑
〈i j 〉
d †
i ti j d j , d i =
(
P+
i di
P−
i di
)
, ti j =
(
t++i j t+−i j
t−+i j t−−i j
)
. (3.2)
The connection between the elements of the hopping matrix ti j and the initial hopping amplitudes is as
follows
t−− = t1 , t1 = t−− ,
t+− = t−+ = t1 + t2 , t2 = t+−(−+) − t−− ,
t++ = t1 +2t2 + t3 , t3 = t+++ t−−− t+−− t−+ . (3.3)
In infinite dimensions D →∞, the self-energy is nonlocal but an irreducible over the hopping part
Ξ(ω) of matrix Green’s function Gk (ω), when one performs an expansion around the atomic limit, is
local which allows one to develop the DMFT approach for the systems with correlated hopping [20, 21].
In terms of this irreducible part, the lattice Green’s function can be written as follows:
Gk (ω) = [
Ξ−1(ω)− tk
]−1
, (3.4)
where for the hopping matrix, we neglect the possible chirality
tk =
(
t++ t+−
t−+ t−−
)
εk , (3.5)
and assume that each term is proportional to the unperturbed band energy
εk = 1p
D
D∑
α=1
coskα . (3.6)
Below we shall use the hopping amplitude over empty states as an energy unit: t−− = t1 = 1.
An irreducible part Ξ(ω) as well as the matrix of the λ-fields Λ(ω) = ∥∥λαβ(ω)
∥∥ are solutions of the
system of equations
1
N
∑
k
[
Ξ−1(ω)− tk
]−1 = [
Ξ−1(ω)−Λ(ω)
]−1 = Gimp(ω), (3.7)
where the expressions for the Green’s function Gimp(ω) of the local impurity problem for the Falicov-
Kimball model are known
G++
imp
(ω) = 〈P+〉
ω+µd −U −λ++(ω)
, (3.8)
G−−
imp
(ω) = 〈P−〉
ω+µd −λ−−(ω)
, (3.9)
G+−
imp
(ω) =G−+
imp
(ω) = 0. (3.10)
43704-5
A.M. Shvaika
-6 -4 -2 0 2 4 6
ω
0
0.5
1
1.5
2
Im
G
(ω
)
Im G
++
(ω)
Im G
−−(ω)
Im G
imp
(ω)
a) t
2
/t
1
=0
-6 -4 -2 0 2 4 6
ω
0
0.5
1
1.5
2
Im
G
(ω
)
Im G
++(ω)
Im G
−−(ω)
Im G
imp
(ω)
c) t
2
/t
1
=0.5
-6 -4 -2 0 2 4 6
ω
0
0.5
1
1.5
2
Im
G
(ω
)
Im G
++(ω)
Im G
−−(ω)
Im G
imp
(ω)
e) t
2
/t
1
=1
-6 -4 -2 0 2 4 6
ω
0
0.5
1
1.5
2
Im
G
(ω
)
Im G
++(ω)
Im G
−−(ω)
Im G
imp
(ω)
b) t
2
/t
1
=−0.5
-6 -4 -2 0 2 4 6
ω
0
0.5
1
1.5
2
Im
G
(ω
)
Im G
++(ω)
Im G
−−(ω)
Im G
imp
(ω)
d) t
2
/t
1
=−1
-6 -4 -2 0 2 4 6
ω
0
0.5
1
1.5
2
Im
G
(ω
)
Im G
++(ω)
Im G
−−(ω)
Im G
imp
(ω)
f) t
2
/t
1
=−2
Figure 2. (Color online) Evolution of the one-particle d.o.s. with a change of the correlated hopping am-
plitude t2 (t3 = 0, t1 = 1,U = 2) (from reference [21]).
Previous investigations of the Falicov-Kimball model with correlated hopping [20, 21] have elucidated
several special cases with a strong reconstruction of the one-particle d.o.s., which in the considered case
is equal to:
Ad (ω) =− 1
π
∑
α,β=±
ImGαβ
imp
(ω). (3.11)
These crossover points correspond to the special forms of the hopping matrix (3.5):
1. The regular Falicov-Kimball model without correlated hopping corresponds to the case of t2 = t3 =
0, when all components of the hopping matrix (3.5) are the same tαβ = t1 and its determinant is
equal to zero dettk = 0 [figure 2 (a)].
2. The case when one or both of the diagonal components of the hopping matrix (3.5) are equal to
zero, e.g. t++ = 0 for t2 =− 1
2 (t1+t3), with strong reduction of the band widths and an orthogonality
singularity at the edge of one or both bands [figure 2 (b)].
3. The case of diagonal hopping matrix (3.5), when t2 = −t1 and t+− = t−+ = 0 and d.o.s. consists of
independent bands [figure 2 (d)].
In other cases, the correlated hopping results in the break of the electron-hole symmetry with strong
asymmetry of the d.o.s. and in the spread of the bands due to an effective increase of the hopping am-
plitude. The similar behaviour could be also predicted for the transport function and will be considered
below.
The Falicov-Kimball model with correlated hopping is similar to the binary alloy model with off-
diagonal disorder, and the first calculations of the thermoelectric power for the last one were done by
Hoshino and Niizeki [33] based on the Mott’s relation and for the special case of the zero determinant
dettk = 0 of hopping matrix (3.5). Below we shall provide a more general approach.
First of all, one has to rederive the charge and energy current operators in the presence of correlated
hopping. Starting from the expression for charge polarization which contains contributions both from
the d and f particles P̂ =∑
i R i
(
zd n̂i d + z f n̂i f
)
, one can get an expression for the charge current in the
presence of correlated hopping
ĵ c =−izd
∑
〈i j 〉
αβ=±
(
R i −R j
)
tαβi j Pα
i d †
i d j Pβ
j . (3.12)
The dc electric conductivity
σdc =−z2
d lim
Ω→0
1
Ω
Imχ(Ω) (3.13)
43704-6
Effect of correlated hopping on thermoelectric properties
is connected with the current-current Green’s function, which in the considered case has the following
form:
χ(iν) = T
2
∑
m
1
N
∑
k
∑
αβα′β′=±
tαβGβα′
k (iωm)tα
′β′
Gβ′α
k (iωm + iν), (3.14)
where for the D → ∞ hypercubic lattice with unperturbed band energy (3.6), a rigorous replacement
(dεk /dk)2 → 1
2 was used. After substituting the components of the lattice Green’s function (3.4) in (3.14)
and performing summation over the momentum k , one can get an expression for the dc conductivity in
a standard form through the transport function
σdc = z2
d L11 , L11 =
+∞∫
−∞
dω
[
−d f (ω)
dω
]
I (ω) . (3.15)
Analytic expressions for the transport function I (ω) in the homogeneous phase in terms of the solutions
of the DMFT impurity problem (3.7) were derived and are presented in the appendix A.
On the other hand, the energy polarization operator for the Falicov-Kimball model with correlated
hopping is equal to
P̂Q =∑
i
R i
[
Uni d ni f −µ f ni f −µd ni d
]+ 1
2
∑
〈i j 〉
αβ=±
(
R i +R j
)
tαβi j Pα
i d †
i d j Pβ
j . (3.16)
Now, from the continuity equation (2.10) we get an expression for the energy current operator
ĵ Q =− i
2
∑
〈i j 〉
αβ=±
(
R i −R j
)
Pα
i d †
i d j Pβ
j
{
tαβi j
[
U (ni f +n j f )−2µd
]+∑
lγ
tαγi l Pγ
l tγβl j
}
(3.17)
and one can check that it is equal to
ĵ Q = 1
2
∑
〈i j 〉
αβ=±
(
R i −R j
)
tαβi j
Pα
i d †
i
d(d j Pβ
j )
dt
− d(Pα
i d †
i )
dt
d j Pβ
j
. (3.18)
This means that the connection (2.14) between the energy and charge current operators holds in the case
of correlated hopping, and the Boltzmann relations (2.16) can be also used in this case. The numerical
results for the charge and heat transport in the systems with correlated hopping are presented in the
next section.
4. Charge and heat transport in presence of correlated hopping
It was already mentioned above that for different values of correlated hopping, different shapes of
the one-particle d.o.s. can be realized (figure 2) and the crossover from one regime to another takes place
at special forms of the hopping matrix (3.5): either its determinant or some matrix elements are equal to
zero. Let us check the behaviour of the charge and heat transport in the vicinity of these crossover points.
First of all, we consider the small values of correlated hopping t2 (below we shall put t3 = 0). In the
absence of correlated hopping t2 = 0 and for small values of the Coulomb interaction U (figure 3), the
shape of the d.o.s. Ad (ω) deviates by a small amount from the unperturbed one (Gaussian for theD →∞
hypercubic lattice). In this case, the transport function I (ω) slightly varies in the vicinity of the chemical
potential value µd = U /2 and approaches the constant value at large frequencies. The last feature is
caused by the Gaussian d.o.s. [31], when for the large energy values we still have an exponentially small
density of states with a finite relaxation rate (see appendix A), which cannot be observed in real systems.
Thus, we consider only the small and moderate temperature values when the Fermi window and its
moments (see figure 1) are inside the features of the d.o.s.
43704-7
A.M. Shvaika
-3 -2 -1 0 1 2 3 4
ω
0
0.5
1
A
d
(ω
)
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.5
1
1.5
2
σ
d
c
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.2
0.4
0.6
0.8
1
κ
e
t
2
= 0.1
t
2
= 0
t
2
= -0.1
-3 -2 -1 0 1 2 3 4
ω
0
5
10
15
20
I(
ω
)
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
-3
-2
-1
0
1
2
S
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.5
1
1.5
Z
T
t
2
= 0.1
t
2
= 0
t
2
= -0.1
Figure 3. (Color online) Density of states Ad (ω), transport function I (ω), and temperature dependences
of the dc electric σdc and thermal κe conductivities, Seebeck coefficient S, and figure of merit Z T at half
filling n f = nd = 1 forU = 0.5, t1 = 1, and t2 = 0.1, 0, −0.1.
The correlated hopping t2 being switched on leads to slight changes of the d.o.s., which become asym-
metric, and to tremendous changes in the shape and values of the transport function I (ω) with a strong
enhancement at the band edges. At low temperatures, the effect of correlated hopping on the electric
σdc and thermal κe conductivities is minor and leads to the appearance of the Seeback effect due to
electron-hole asymmetry. At higher temperature values, we observe an enhancement of all transport co-
efficients because the moments of the Fermi window (figure 1) now start to cover the sharp features of
the transport function at the band edges and this behaviour could be non-physical as a consequence of
the Gaussian d.o.s.
In the Mott insulator phase forU = 2, the effect of small values of correlated hopping is much smaller
(figure 4). Both the d.o.s. Ad (ω) and the transport function I (ω) change only slightly with respect to the
case without correlated hopping t2 = 0. The temperature dependence of the electric σdc and thermal
κe conductivities is almost the same and we observe only an increase of the Seeback effect at low tem-
peratures similar to the one in the doped Mott insulators [31] (non-zero correlated hopping breaks the
-2 -1 0 1 2 3 4
ω
0
0.5
1
A
d
(ω
)
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.05
0.1
0.15
0.2
σ
d
c
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.1
0.2
κ
e
t
2
= 0.1
t
2
= 0
t
2
= -0.1
-2 -1 0 1 2 3 4
ω
0
1
2
I(
ω
)
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
-10
-5
0
5
S
t
2
= 0.1
t
2
= 0
t
2
= -0.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.5
1
1.5
2
Z
T
t
2
= 0.1
t
2
= 0
t
2
= -0.1
Figure 4. (Color online) Same as in figure 3 forU = 2.
43704-8
Effect of correlated hopping on thermoelectric properties
-2 -1 0 1 2 3 4
ω
0
0.5
1
A
d
(ω
)
t
2
= -0.9
t
2
= -1
t
2
= -1.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.05
0.1
0.15
0.2
σ
d
c
t
2
= -0.9
t
2
= -1
t
2
= -1.1
0 0.1 0.2 0.3 0.4 0.5
T
0
0.1
0.2
κ
e
t
2
= -0.9
t
2
= -1
t
2
= -1.1
-2 -1 0 1 2 3 4
ω
0
1
2
I(
ω
)
t
2
= -0.9
t
2
= -1
t
2
= -1.1
0 0.1 0.2 0.3 0.4 0.5
T
-0.5
-0.25
0
0.25
0.5
S
t
2
= -0.9
t
2
= -1
t
2
= -1.1
0 0.1 0.2 0.3 0.4 0.5
T
0
5×10
-3
1×10
-2
Z
T
t
2
= -0.9
t
2
= -1
t
2
= -1.1
Figure 5. (Color online) Density of states Ad (ω), transport function I (ω), and temperature dependences
of the dc electric σdc and thermal κe conductivities, Seebeck coefficient S, and figure of merit Z T at half
filling n f = nd = 1 forU = 2, t1 = 1, and t2 =−0.9, −1, −1.1.
electron-hole symmetry like doping in the case without correlated hopping). It is known that the tem-
perature behaviour of thermoelectric properties in Mott insulator to a great extent is determined by the
temperature dependence of the chemical potential [31, 32], and due to the numerical issues we were not
able to determine the chemical potential values in the gap at very low temperatures with a precision
sufficient to get smooth dependences of the Seebeck coefficient.
For the opposite case of almost independent bands, when t2 ≈ −1 and off-diagonal elements of the
hopping matrix (3.5) are small and change the sign (t+− = t−+ → 0), the d.o.s. and transport function are
very smooth in the Fermi window, which results in themetallic behaviour of the electric and thermal con-
ductivity and in weak thermoelectric properties (figure 5). The Coulomb interactionU is less important
in this case.
In the vicinity of another crossover point at t2 =−0.5 (no direct hopping between the sites occupied
-3 -2 -1 0 1 2 3 4
ω
0
0.5
1
1.5
2
A
d
(ω
)
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
0
0.5
1
σ
d
c
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
0
0.1
0.2
κ
e
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
-3 -2 -1 0 1 2 3 4
ω
0
5
10
15
I(
ω
)
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
-3
-2
-1
0
1
2
S
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
0
0.5
1
1.5
Z
T
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
Figure 6. (Color online) Density of states Ad (ω), transport function I (ω), and temperature dependences
of the dc electric σdc and thermal κe conductivities, Seebeck coefficient S, and figure of merit Z T at half
filling n f = nd = 1 forU = 0.5, t1 = 1, and t2 =−0.4, −0.5, −0.6.
43704-9
A.M. Shvaika
by f -particles: t++ = 0) the behaviour is altogether different (figure 6). Exactly at the t2 =−0.5 value, the
d.o.s. Ad (ω) contains two bands separated by a small gap, and the upper band possesses a singularity at
the gap edge. It could be imagined that similar features should be observed for the transport function
I (ω). However, the transport function displays only a finite abrupt change instead of a singularity above
the gap. For the value of correlated hopping t2 =−0.6, the gap closes and the transport function is much
smoother. On the other side of the crossover point, for the t2 = −0.4 value, the gap also closes and the
shape of the d.o.s. does not change significantly in comparison with the case of t2 =−0.6 value, but now
the transport function displays a strong enhancement with a narrow peak. It should be noted that in
many theoretical simulations, the calculation of two-particle quantities, including the transport function,
is problematic, and sometimes an approximation I (ω) = πΓAd (ω), replacing by one-particle quantity,
is used instead (see, e.g. [13]). Our results definitely show that such an approximation is not valid and
produces altogether different transport functions for many cases.
Different shapes of the transport function in the vicinity of the t2 = −0.5 value manifest themselves
in different transport properties. Exactly at the crossover point t2 =−0.5, the temperature dependences
of the electric σdc and thermal κe conductivities as well as of the Seebeck coefficient are similar to the
one in the doped small gap Mott insulator. For the t2 = −0.6 value, a metallic behaviour is observed,
i.e., weak temperature dependence of the electric conductivity and small thermoelectric power. On the
other side of the crossover point, an enhancement of all transport coefficients is observed which is more
prominent for the electric conductivity and Seebeck coefficient than for the thermal conductivity and
causes an increase of the thermoelectric figure of merit Z T .
For large values of the Coulomb interaction U = 2, the behaviour is quite different (figure 7). Now,
-2 -1 0 1 2 3 4
ω
0
0.5
1
1.5
2
A
d
(ω
)
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
0
0.05
0.1
0.15
0.2
σ
d
c
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
0
0.1
0.2
κ
e
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
-2 -1 0 1 2 3 4
ω
0
2
4
6
8
10
I(
ω
)
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
-10
-5
0
5
S
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
0 0.1 0.2 0.3 0.4 0.5
T
0
2
4
6
8
10
Z
T
t
2
= -0.4
t
2
= -0.5
t
2
= -0.6
Figure 7. (Color online) Same as in figure 6 forU = 2.
the d.o.s. always possesses the gap and the width of the upper band strongly reduces. Exactly at the
crossover point t2 = −0.5, there is a singularity on the d.o.s. and there is only an abrupt change on the
transport function similar to the small-U case. For the t2 = −0.6 value, sharp features on the d.o.s. and
transport function become smoother. On the other hand, for the t2 = −0.4 value, the d.o.s. is almost the
same and the transport function develops a sharp peak. But now this peak is outside the Fermi window
and less effects the transport properties. For the t2 = −0.5 value, the chemical potential is sticked at
the upper edge of the lower band which results in a weak metallic behaviour of electric and thermal
conductivities though the thermoelectric power behaves like in the dopedMott insulator. For other values
of the correlated hopping, t2 =−0.6 and−0.4, the chemical potential shifts inside the gap, approaching the
features of the transport function in the upper band, and the electric and thermal conductivities manifest
the temperature dependences typical of the Mott insulator phase. Now, the Seebeck coefficient shows
an enhancement and anomalous behaviour at low temperatures, when the extrema of the moments of
Fermi window (figure 1) overlap with the features of the transport function following the change of the
43704-10
Effect of correlated hopping on thermoelectric properties
-2 -1 0 1 2 3 4
ω
0
0.5
1
1.5
A
d
(ω
)
n
f
= 0.1
n
f
= 0.2
n
f
= 0.3
n
f
= 0.4
n
f
= 0.5
n
f
= 0.6
n
f
= 0.7
n
f
= 0.8
n
f
= 0.9
0 0.1 0.2 0.3
T
0
0.05
0.1
0.15
σ
d
c
n
f
= 0.1
n
f
= 0.2
n
f
= 0.3
n
f
= 0.4
n
f
= 0.5
n
f
= 0.6
n
f
= 0.7
n
f
= 0.8
n
f
= 0.9
0 0.1 0.2 0.3
T
0
0.25
0.5
0.75
1
κ
e
n
f
= 0.1
n
f
= 0.2
n
f
= 0.3
n
f
= 0.4
n
f
= 0.5
n
f
= 0.6
n
f
= 0.7
n
f
= 0.8
n
f
= 0.9
-2 -1 0 1 2 3 4
ω
0
10
20
30
40
50
I(
ω
)
n
f
= 0.1
n
f
= 0.2
n
f
= 0.3
n
f
= 0.4
n
f
= 0.5
n
f
= 0.6
n
f
= 0.7
n
f
= 0.8
n
f
= 0.9
0 0.1 0.2 0.3
T
-10
-5
0
5
10
S
n
f
= 0.1
n
f
= 0.2
n
f
= 0.3
n
f
= 0.4
n
f
= 0.5
n
f
= 0.6
n
f
= 0.7
n
f
= 0.8
n
f
= 0.9
0 0.1 0.2 0.3
T
0
5
10
15
Z
T
n
f
= 0.1
n
f
= 0.2
n
f
= 0.3
n
f
= 0.4
n
f
= 0.5
n
f
= 0.6
n
f
= 0.7
n
f
= 0.8
n
f
= 0.9
Figure 8. (Color online) Density of states Ad (ω), transport function I (ω), and temperature dependences of
the dc electric σdc and thermal κe conductivities, Seebeck coefficient S, and figure of merit Z T forU = 2,
t1 = 1, t2 =−0.3, and different doping of f -states n f (nd +n f = 1).
temperature.
Such an enhancement of thermoelectric properties by correlated hopping is more prominent with
doping. In figure 8 we present the results for the case of Mott insulator U = 2 and for the intermediate
value of correlated hopping t2 = −0.3. The doping by f -particles over the half-filling n f Ê 0.5 and a
simultaneous decrease of the d -states occupation (nd = 1−n f ) reduces the peak on a transport function
below the chemical potential, while the peak above it enhances and moves inside the Fermi window. Due
to the fancy profile of the transport function, different trajectories of the chemical potential and Fermi
window cover its features, with the change of temperature at different doping levels, at different energies
and produce anomalous temperature dependences of the Seebeck coefficient and an enhancement of
thermoelectric properties at low temperatures. On the other hand, a decrease of the doping of f -states
n f É 0.5makes the d.o.s. and transport function more symmetric within the Fermi window and reduces
the thermoelectric properties.
5. Conclusions
In this article we have investigated the effect of correlated hopping on the charge and heat transport
and thermoelectric power in a correlated material described by the Falicov-Kimball model. Exact expres-
sions for the transport function are derived for a homogeneous phase, which is assumed to be valid at
low temperatures. We do not consider the possible phase transitions into the modulated phases [22–25],
which should be the subject of a separate investigation.
Depending on the value of correlated hopping, the crossover points which separate the regions with
different shapes of the d.o.s. and transport function are clarified. The temperature dependences of the
electrical and thermal conductivities and thermoelectric power are strongly effected by the presence of
singularities and peaks on the transport function and by the temperature evolution of the chemical po-
tential. The largest enhancement of thermoelectric properties is observed for the values of correlated
hopping −t1/2 < t2 < 0, when a direct hopping between the same many-body states at different sites
reduces and an indirect hopping becomes important. It should be noted that the calculations of the corre-
lated hopping amplitudes for different compounds, see e.g. [11, 12], provide absolute values close to this
interval but disagree by its sign, which is the main factor according to our results.
Unfortunately, the use of the Gaussian d.o.s. produces non-physical results for the transport function
outside the bands, which does not allow us to get reasonable values for the chemical potential at very low
and very high temperatures and calls for additional investigations with another unperturbed density of
states.
43704-11
A.M. Shvaika
A. Transport function in presence of correlated hopping
After substitution of equation (3.14) in (3.13) we get an expression for the transport function in the
form
I (ω) = 1
2π
∑
αβα′β′
tαβtα
′β′ 1
N
∑
k
ImGβα′
k (ω) ImGβ′α
k (ω). (A.1)
The components of the lattice Green’s function (3.4) can be written as follows:
Gβα
k (ω) = Aβα(ω)−Bβαεk
C (ω)−D(ω)εk +Pε2
k
, (A.2)
where
A++(ω) = Z−−(ω), A+−(ω) =−Z+−(ω),
A−+(ω) =−Z−+(ω), A−−(ω) = Z++(ω), (A.3)
B++ = t−−, B+− =−t+−,
B−+ =−t−+, B−− = t++, (A.4)
and
C (ω) = Z++(ω)Z−−(ω)−Z+−(ω)Z−+(ω) = detZ(ω),
D(ω) = Z++(ω)t−−+Z−−(ω)t++−Z+−(ω)t−+−Z−+(ω)t+−,
P = t++t−−− t+−t−+ = det
(
t++ t+−
t−+ t−−
)
, (A.5)
and we introduce a notation for the inverse irreducible part Z(ω) =Ξ−1(ω). The final analytic expression
depends on the value of parameter P .
In the case of P = 0, we have dettk = 0, and the expression for the transport function is as follows:
I (ω) = 1
4π
∑
αβα′β′
tαβtα
′β′
{
Re
[
Aβα′ (ω)Aβ′α(ω)
D2(ω)
F ′
∞
(
C (ω)
D(ω)
)]
− Re
[
Aβα′ (ω)Aβ′α(ω)
]
|D(ω)|2 ImF∞
(
C (ω)
D(ω)
) /
Im
C (ω)
D(ω)
}
, (A.6)
where
F∞(ζ) =
+∞∫
−∞
dε
ρ(ε)
ζ−ε , F ′
∞(ζ) = dF∞(ζ)
dζ
(A.7)
are a Hilbert transform of the unperturbed density of states and its first derivative, respectively, and for
the Gaussian d.o.s.
ρ(ε) = 1
W
p
π
e−ε
2/W 2
(A.8)
we have
F ′
∞(ζ) = 2
W 2 [1−ζF∞(ζ)] . (A.9)
43704-12
Effect of correlated hopping on thermoelectric properties
In the general case of P , 0, the transport function takes a more complicated form
I (ω) = 1
8πP 2
∑
αβα′β′
tαβtα
′β′
{
2Re
[[
Aβα′ (ω)−Bβα′E1(ω)
][
Aβ′α(ω)−Bβ′αE1(ω)
]
[E2(ω)−E1(ω)]2 F ′
∞
(
E1(ω)
)]
+2Re
[[
Aβα′ (ω)−Bβα′E2(ω)
][
Aβ′α(ω)−Bβ′αE2(ω)
]
[E1(ω)−E2(ω)]2 F ′
∞
(
E2(ω)
)]
− 1
ImE1(ω)
Im
[
F∞
(
E1(ω)
)
[E2(ω)−E1(ω)]
[
E∗
2 (ω)−E1(ω)
] ([
Aβα′ (ω)−Bβα′E1(ω)
][
A∗
β′α(ω)−B∗
β′αE1(ω)
]
+ [
A∗
βα′ (ω)−B∗
βα′E1(ω)
][
Aβ′α(ω)−Bβ′αE1(ω)
]) ]
− 1
ImE2(ω)
Im
[
F∞
(
E2(ω)
)
[E1(ω)−E2(ω)]
[
E∗
1 (ω)−E2(ω)
] ([
Aβα′ (ω)−Bβα′E2(ω)
][
A∗
β′α(ω)−B∗
β′αE2(ω)
]
+ [
A∗
βα′ (ω)−B∗
βα′E2(ω)
][
Aβ′α(ω)−Bβ′αE2(ω)
]) ]}
, (A.10)
where E1,2(ω) are roots of the denominator in equation (A.2)
E1,2(ω) = D(ω)
2P
±
√[
D(ω)
2P
]2
− C (ω)
P
. (A.11)
For large energy values and for unperturbed Gaussian d.o.s., the imaginary parts in the last terms of
equations (A.6) and (A.10) are small but their ratios are finite and produce a strong enhancement of the
transport function and make its shape altogether different from the one-particle d.o.s.
References
1. Mahan G., Sales B., Sharp J., Phys. Today, 1997, 50, No. 3, 42–47; doi:10.1063/1.881752.
2. Costi T.A., Hewson A.C., Zlatić V., J. Phys.: Condens. Matter, 1994, 6, No. 13, 2519; doi:10.1088/0953-8984/6/13/013.
3. Hubbard J., Proc. R. Soc. London, Ser. A, 1963, 276, No. 1365, 238–257; doi:10.1098/rspa.1963.0204.
4. Foglio M.E., Falicov L.M., Phys. Rev. B, 1979, 20, No. 11, 4554–4559; doi:10.1103/PhysRevB.20.4554.
5. Simón M.E., Aligia A.A., Phys. Rev. B, 1993, 48, No. 10, 7471–7477; doi:10.1103/PhysRevB.48.7471.
6. Blackman J.A., Esterling D.M., Berk N.F., Phys. Rev. B, 1971, 4, No. 8, 2412–2428; doi:10.1103/PhysRevB.4.2412.
7. Hirsch J.E., Physica C, 1989, 158, No. 3, 326–336; doi:10.1016/0921-4534(89)90225-6.
8. Bułka B.R., Phys. Rev. B, 1998, 57, No. 17, 10303–10306; doi:10.1103/PhysRevB.57.10303.
9. Didukh L., Skorenkyy Yu., Dovhopyaty Yu., Hankevych V., Phys. Rev. B, 2000, 61, No. 12, 7893–7908;
doi:10.1103/PhysRevB.61.7893.
10. Kollar M., Vollhardt D., Phys. Rev. B, 2001, 63, No. 4, 045107; doi:10.1103/PhysRevB.63.045107.
11. Meir Y., Hirose K., Wingreen N.S., Phys. Rev. Lett., 2002, 89, No. 19, 196802; doi:10.1103/PhysRevLett.89.196802.
12. Hübsch A., Lin J.C., Pan J., Cox D.L., Phys. Rev. Lett., 2006, 96, No. 19, 196401; doi:10.1103/PhysRevLett.96.196401.
13. Tooski S.B., Ramšak A., Bułka B.R., Žitko R., New J. Phys., 2014, 16, No. 5, 055001;
doi:10.1088/1367-2630/16/5/055001.
14. Eckholt M., García-Ripoll J.J., New J. Phys., 2009, 11, No. 9, 093028; doi:10.1088/1367-2630/11/9/093028.
15. Jürgensen O., Sengstock K., LühmannD.-S., Phys. Rev. A, 2012, 86, No. 4, 043623; doi:10.1103/PhysRevA.86.043623.
16. Liberto M.D., Creffield C.E., Japaridze G.I., Smith C.M., Phys. Rev. A, 2014, 89, 013624;
doi:10.1103/PhysRevA.89.013624.
17. Falicov L.M., Kimball J.C., Phys. Rev. Lett., 1969, 22, No. 19, 997–999; doi:10.1103/PhysRevLett.22.997.
18. Freericks J.K., Zlatić V., Rev. Mod. Phys., 2003, 75, No. 4, 1333–1382; doi:10.1103/RevModPhys.75.1333.
19. Schiller A., Phys. Rev. B, 1999, 60, No. 23, 15660–15663; doi:10.1103/PhysRevB.60.15660.
20. Shvaika A.M., Phys. Rev. B, 2003, 67, No. 7, 075101; doi:10.1103/PhysRevB.67.075101.
21. Shvaika A.M., Phys. Rev. B, 2010, 82, No. 11, 119901; doi:10.1103/PhysRevB.82.119901.
22. Gajek Z., Lemański R., Acta Phys. Pol. B, 2001, 32, No. 10, 3473–3476.
23. Wojtkiewicz J., Lemański R., Phys. Rev. B, 2001, 64, No. 23, 233103; doi:10.1103/PhysRevB.64.233103.
43704-13
http://dx.doi.org/10.1063/1.881752
http://dx.doi.org/10.1088/0953-8984/6/13/013
http://dx.doi.org/10.1098/rspa.1963.0204
http://dx.doi.org/10.1103/PhysRevB.20.4554
http://dx.doi.org/10.1103/PhysRevB.48.7471
http://dx.doi.org/10.1103/PhysRevB.4.2412
http://dx.doi.org/10.1016/0921-4534(89)90225-6
http://dx.doi.org/10.1103/PhysRevB.57.10303
http://dx.doi.org/10.1103/PhysRevB.61.7893
http://dx.doi.org/10.1103/PhysRevB.63.045107
http://dx.doi.org/10.1103/PhysRevLett.89.196802
http://dx.doi.org/10.1103/PhysRevLett.96.196401
http://dx.doi.org/10.1088/1367-2630/16/5/055001
http://dx.doi.org/10.1088/1367-2630/11/9/093028
http://dx.doi.org/10.1103/PhysRevA.86.043623
http://dx.doi.org/10.1103/PhysRevA.89.013624
http://dx.doi.org/10.1103/PhysRevLett.22.997
http://dx.doi.org/10.1103/RevModPhys.75.1333
http://dx.doi.org/10.1103/PhysRevB.60.15660
http://dx.doi.org/10.1103/PhysRevB.67.075101
http://dx.doi.org/10.1103/PhysRevB.82.119901
http://dx.doi.org/10.1103/PhysRevB.64.233103
A.M. Shvaika
24. Wojtkiewicz J., Lemański R., Acta Phys. Pol. B, 2001, 32, No. 10, 3467–3472.
25. Čenčariková H., Farkašovský P., Condens. Matter Phys., 2011, 14, No. 4, 42701; doi:10.5488/CMP.14.42701.
26. Mahan G.D., Many-Particle Physics, 3rd Edn., Kluwer Academic/Plenum Publishers, New York, 2000.
27. Freericks J.K., Zlatić V., Shvaika A.M., Phys. Rev. B, 2007, 75, No. 3, 035133; doi:10.1103/PhysRevB.75.035133.
28. Jonson M., Mahan G.D., Phys. Rev. B, 1980, 21, No. 10, 4223–4229; doi:10.1103/PhysRevB.21.4223.
29. Jonson M., Mahan G.D., Phys. Rev. B, 1990, 42, No. 15, 9350–9356; doi:10.1103/PhysRevB.42.9350.
30. Mahan G.D., Sofo J.O., Proc. Natl. Acad. Sci. USA, 1996, 93, No. 15, 7436–7439; doi:10.1073/pnas.93.15.7436.
31. Freericks J.K., Demchenko D.O., Joura A.V., Zlatić V., Phys. Rev. B, 2003, 68, No. 19, 195120;
doi:10.1103/PhysRevB.68.195120.
32. Zlatić V., Boyd G.R., Freericks J.K., Phys. Rev. B, 2014, 89, No. 15, 155101; doi:10.1103/PhysRevB.89.155101.
33. Hoshino K., Niizeki K., J. Phys. Soc. Jpn., 1975, 38, No. 5, 1320–1327; doi:10.1143/JPSJ.38.1320.
Вплив корельованого переносу на термоелектричнi
властивостi: точний розв’язок для моделi Фалiкова-Кiмбала
А.М.Швайка
Iнститут фiзики конденсованих систем НАН України, вул. I. Свєнцiцького, 1, 79011 Львiв, Україна
Дослiджено вплив корельованого переносу на перенос заряду i тепла в моделi Фалiкова-Кiмбала. В рам-
ках теорiї динамiчного середнього поля отримано точнi розв’язки для електро- та теплопровiдностi i
термо-е.р.с. Дослiджено температурнi залежностi коефiцiєнтiв переносу для певних значень корельова-
ного переносу,що вiдповiдають суттєвiй перебудовi густини станiв та транспортної функцiї. Встановлено
випадки, якi вiдповiдають значному покращенню термоелектричних властивостей.
Ключовi слова: термо-е.р.с., корельований перенос, модель Фалiкова-Кiмбала
43704-14
http://dx.doi.org/10.5488/CMP.14.42701
http://dx.doi.org/10.1103/PhysRevB.75.035133
http://dx.doi.org/10.1103/PhysRevB.21.4223
http://dx.doi.org/10.1103/PhysRevB.42.9350
http://dx.doi.org/10.1073/pnas.93.15.7436
http://dx.doi.org/10.1103/PhysRevB.68.195120
http://dx.doi.org/10.1103/PhysRevB.89.155101
http://dx.doi.org/10.1143/JPSJ.38.1320
Introduction
Macroscopic and microscopic levels in the description of thermoelectric effect
Charge and energy current operators and transport coefficients for the Falicov-Kimball model with correlated hopping
Charge and heat transport in presence of correlated hopping
Conclusions
Transport function in presence of correlated hopping
|