SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Open Badges Nano Review

Effects of interdot hopping and Coulomb blockade on the thermoelectric properties of serially coupled quantum dots

David M T Kuo1 and Yia-Chung Chang2*

Author affiliations

1 Department of Electrical Engineering and Department of Physics, National Central University, Chungli, 320, Taiwan

2 Research Center for Applied Sciences, Academic Sinica, Taipei, 115, Taiwan

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2012, 7:257  doi:10.1186/1556-276X-7-257

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/257

Received:4 February 2012
Accepted:16 May 2012
Published:16 May 2012

© 2012 Kuo and Chang; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We have theoretically studied the thermoelectric properties of serially coupled quantum dots (SCQDs) embedded in an insulator connected to metallic electrodes. In the framework of Keldysh Green’s function technique, the Landauer formula of transmission factor is obtained using the equation of motion method. Based on such analytical expressions of charge and heat currents, we calculate the electrical conductance, Seebeck coefficient, electron thermal conductance, and figure of merit (ZT) of SCQDs in the linear response regime. The effects of interdot hopping and electron Coulomb interactions on ZT are analyzed. We demonstrate that ZT is not a monotonic increasing function of interdot electron hopping strength (tc). We also show that in the absence of phonon thermal conductance, SCQD can reach the Carnot efficiency as tcapproaches zero.



Recently, many considerable studies have been devoted to seeking efficient thermoelectric materials with the figure of merit (ZT) larger than 3 because there exist potential applications of solid-state thermal devices such as coolers and power generators [1-6]. Some theoretical efforts have pointed out that a single quantum dot (QD) junction system can have a very impressive ZT in the absence of phonon conductance [7-9]. However, in practice, it is difficult to maintain a large temperature gradient which is needed to produce sufficient temperature difference across the nanoscale junction. To reduce the temperature gradient across the QD junction, it is essential to consider many serially coupled quantum dots (SCQDs) [1,5]. The transport property of a junction involving N serially coupled QDs with strong electron Coulomb interactions is one of the most challenging topics of condensed matter physics. To gain some insight, we investigate in the present paper the thermoelectric effect of serially coupled quantum dots (SCQDs) as shown in the inset of Figure 1a.

thumbnailFigure 1. Functions of temperature. (a) (ZT)0 and (b) ZT as functions of temperature for various interdot hopping strengths (tc = 0.1,0.5,1,3, and 4Γ0). E= EF + 30Γ0, U= 30 Γ0, U,j =10 Γ0, and ΓL = Γ = Γ0. STM, scanning tunneling microscope.

It has been shown that the transport properties of the SCQD system exhibit several interesting features, including current rectification (due to the Pauli spin blockade), negative differential conductance, nonthermal broadening of tunneling current, and coherent tunneling in the Coulomb blockade regime [10]. Although many theoretical investigations of the above phenomena have been reported, most of them did not investigate the thermoelectric properties of SCQDs [11-13]. This study investigates the ZT of a SCQD embedded in a semiconductor nanowire with small phonon thermal conductance [4]. It is expected that the SCQD system has a potential to enhance the ZT of nanowires. Here, we consider nanoscale semiconductor QDs, in which the energy level separations are much larger than their on-site Coulomb interactions and thermal energies. Thus, only one energy level for each quantum dot needs to be considered. A two-level Anderson model [13] is employed to simulate the SCQD junction system.

Theoretical model

Using Keldysh-Green’s function technique [13], the charge and heat currents of SCQD connected to metallic electrodes are given by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M1">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M2">View MathML</a>


where <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M3">View MathML</a> is the transmission factor. <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M4">View MathML</a> denotes the Fermi distribution function for the left (right) electrode. The left (right) chemical potential is given by μL(μR). TL(R) denotes the equilibrium temperature of the left (right) electrode. e and h denote the electron charge and Planck’s constant, respectively. <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M5">View MathML</a> denotes the transmission function, which can be calculated by evaluating the on-site retarded Green’s function (GF) and lesser GF [13]. The indices and j denote the th QD and the jth QD, respectively. Based on the equation of motion method, we can obtain analytical expressions of all GFs in the Coulomb blockade regime. Details are provided in [13]. The transmission function in the weak interdot limit (tc/U≪ 1, where tc and U denote the electron interdot hopping strength and on-site Coulomb interaction, respectively) can be recast into the following form:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M6">View MathML</a>


where Im means taking the imaginary part of the function that follows and

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M7">View MathML</a>


Γ=L(1),R(2)(ε) denotes the tunnel rate from the left electrode to dot A (E1) and the right electrode to dot B (E2), which is assumed to be energy- and bias-independent for simplicity. μ= εE + iΓ/2. We can assign the following physical meaning to Equation 3. The sum in Equation 3 is over eight possible configurations labeled by m. We consider an electron (of spin σ) entering level , which can be either occupied (with probability <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M8">View MathML</a>) or empty (with probability <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M9">View MathML</a>). For each case, the electron can hop to level j, which can be empty (with probability <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M10">View MathML</a>), singly occupied in a spin <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M11">View MathML</a> state (with probability <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M12">View MathML</a>) or spin σ state (with probabilitybj,σ = Nj,σ cj), or a double-occupied state (with probabilitycj). Thus, the probability factors associated with the eight configurations appearing in Equation 4 become <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M13">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M14">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M15">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M16">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M17">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M18">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M19">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M20">View MathML</a>. Σm in the denominator of Equation 4 denotes the self-energy correction due to Coulomb interactions and coupling with level j (which couples with the other electrode) in configuration m. We have <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M21">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M22">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M23">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M24">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M25">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M26">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M27">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M28">View MathML</a>. E, U, and U,j denote, respectively, the energy levels of dots, intradot Coulomb interactions, and interdot Coulomb interactions. Here, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M29">View MathML</a>Im Σj denotes the effective tunneling rate from level l to the other electrode through level j in configuration m. For example, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M30">View MathML</a>Im <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M31">View MathML</a>. It is noted that <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M32">View MathML</a> has a numeratorΓj for all configurations. Furthermore, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M33">View MathML</a> is just the on-site single-particle retarded GF for level as given in Equation (A16) of [13], and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M34">View MathML</a> corresponds to its partial GF in configuration m. The transmission function written this way has the same form as Landauer’s formula for a single QD with multiple energy levels including intralevel and interlevel electron Coulomb interactions [14,15].

The probability factors of Equation 3 are determined by the thermally averaged one-particle occupation number and two-particle correlation functions, which can be obtained by solving the on-site lesser Green’s functions [13]:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M35">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M36">View MathML</a>


Note that j in Equations 3, 5, and 6. In the linear response regime, Equations 1 and 2 can be rewritten as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M37">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M38">View MathML</a>


where ΔV = μL μR and ΔT = TL TR are the voltage and temperature differences across the junction, respectively. Thermoelectric response functions in Equations 7 and 8 are given by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M39">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M40">View MathML</a>


<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M41">View MathML</a>



<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M42">View MathML</a>


Here, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M43">View MathML</a> and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M44">View MathML</a> are evaluated in the equilibrium condition. It can be shown that the Onsager relation <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M45">View MathML</a> is preserved. These thermoelectric response functions can also be found in [7], where authors investigated the thermoelectric properties of a single QD.

If the system is in an open circuit, the electrochemical potential will form in response to a temperature gradient; this electrochemical potential is known as the Seebeck voltage (Seebeck effect). The Seebeck coefficient (amount of voltage generated per unit temperature gradient) is defined as <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M46">View MathML</a>. To judge whether the system is able to generate power or refrigerate efficiently, we need to consider the figure of merit, which is given by

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M47">View MathML</a>


Here, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M48">View MathML</a> is the electrical conductance, and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M49">View MathML</a> is the electron thermal conductance. (ZT)0 represents the ZT value in the absence of phonon thermal conductance, κph. For simplicity, we assume κph = κph,0Fs[16-18]. <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M50">View MathML</a> is the universal phonon thermal conductance arising from acoustic phonon confinement in a nanowire [16-18], which was confirmed in the phonon wave guide [19]. The expression of κph = κph,0Fs with Fs = 0.1 can explain well the phonon thermal conductance of silicon nanowire with surface states calculated by the first-principles method [16]. The dimensionless scattering factor Fs arises from phonon scattering with surface impurities or surface defects of quantum dots [1,16]. Here, we adopt Fs = 0.02, which is smaller than Fs = 0.1 because QDs can enhance the phonon scattering rates and reduce phonon thermal conduction as pointed out in [1].

Results and discussion

Here, we consider the case of identical QDs in the optimization of ZT, although it is understood that the size fluctuation of QDs can suppress ZT [13]. In Figure 1a,b we plot (ZT)0 and ZT as a function of temperature for various electron hopping strengths. We adopt the following physical parameters: E=EF + 30Γ0, U= 30Γ0, U,j = 10Γ0, and ΓL = ΓR = Γ = 1Γ0. All energy scales are in the units of the characteristic energy, Γ0. In Figure 1a, we see that (ZT)0 increases with decreasing tc and diverges as tc → 0. This behavior can be proved rigorously as we shall illustrate below. It implies that SCQD can reach the Carnot efficiency in the limit of extremely weak interdot coupling, if one can fully suppress κph, for example, by inserting a nanoscale vacuum layer to block the phonon heat current. Although it would be a challenging task to implement a vacuum layer between one of the electrodes and SCQD, it may be possible to test this idea out via a scanning tunneling microscopic experiment using a setup as shown in the inset of Figure 1a. In Figure 1b, we see that ZT is enhanced with increasing tc until tc reaches 3Γ0, and it becomes reduced for higher tc.

The diverging behavior of (ZT)0 with respect to tc is further illustrated in Figure 2d. The maximum ZT is suppressed in the presence of κph, which is much larger than κe for small tc. The behaviors of ZT shown in Figure 1b are mostly determined by the power factor (S2Ge). Once tc is larger than 3Γ0, the reduction of S2 is faster than the increase of Ge. This explains why the maximum ZT at tc = 4Γ0 is smaller than that at tc = 3Γ0. The location of ZTmax is nearly independent of tc, and it occurs near kB T = 8.8Γ0. For comparison, we also show the results (curves with triangle marks) for the case without electron Coulomb interactions in Figure 1b. It is seen that the maximum ZT is enhanced when we turn off the electron Coulomb interactions. Such a behavior is similar to that of a single QD with multiple energy levels [7,8]. The effect of electron Coulomb interactions is significant only for temperature between 6Γ0 and 50Γ0. Namely, the electron Coulomb interactions are negligible when U/(kBT)≫1or U/(kBT)≪1.

To further understand the behavior of ZT with respect to tc, we plot the electrical conductance (Ge), Seebeck coefficient (S), electrical conductance κe, and (ZT)0 as functions of tc in Figure 2 for various detuning energies, ΔEEF. When E is close to the Fermi energy, Ge and κe are enhanced, whereas S and (ZT)0 are suppressed. The behavior of (ZT)0 at Δ = 30Γ0 in the absence of Coulomb interactions is also shown by the curve with triangles, which has a similar trend as the solid line. Thus, it is instructive to analyze (ZT)0 in the absence of Coulomb interactions. Keeping the leading order of <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M51">View MathML</a>, we have <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M52">View MathML</a>, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M53">View MathML</a>, and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M54">View MathML</a>. Therefore, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M55">View MathML</a>, S = − Δ/eT is independent on tc, and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M56">View MathML</a> vanishes up to <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M57">View MathML</a>. Thus, the leading order of κe is <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M58">View MathML</a>. This indicates that (ZT) <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/257/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/257/mathml/M59">View MathML</a> in the limit of weak interdot hopping.

thumbnailFigure 2. Functions oftcat kBT = 5Γ0. (a) Electrical conductance (Ge), (b) Seebeck coefficient (S), (c) electrical thermal conductance (κe), and (d) (ZT)0 as functions of tcat kBT = 5Γ0for Δ = 10Γ0(dotted curves), 20Γ0(dashed curves), and 30Γ0(solid curves). Other parameters are the same as those of Figure 1.

Figure 3 shows ZT as a function of Δ = EEF for various electron hoping strengths at kBT = 10Γ0. Other physical parameters are kept the same as those for Figure 1. When tc = 0.1Γ0, the maximum ZT (ZTmax) occurs at near Δ = 27Γ0. The peak position only shifts slightly to higher Δwith increasing tc. We have ZTmax = 2.79 and 3.18 for tc = 1Γ0 and 3Γ0, respectively. However, at tc = 4Γ0, we have ZTmax = 3.07, which is smaller than ZTmax for tc = 3Γ0. Thus, it also illustrates that ZT is not a monotonically increasing function of tc. We further calculated ZT as a function of tc for Δ = 10,20,30Γ0 and kBT = 10Γ0 in the presence of κph and found that again, ZT is not a monotonically increasing function of tc (not shown here). We conclude that as long as κph dominates over κe, the tc dependence of ZT is mainly determined by the power factor S2Ge, where the behaviors of Ge and S are similar to the results shown in Figure 2a,b. When tc/Γ0≤1, Ge increases much faster than the reduction ofS2 for increasing tc, and the power factor slowly reaches the maximum when tc approaches 3Γ0. When tc > 3Γ0, the power factor decreases due to the fast reduction ofS2 which prevails over the increase of Ge. The curve with triangle marks is for tc = 3Γ0 in the absence of Coulomb interaction. We see that ZTmax is larger when U=U,j = 0. Based on the results of Figure 3, we conclude that it is important to control the detuning energy Δ for the optimization of ZT.

thumbnailFigure 3. ZT as a function ofΔfor different electron hopping strength at kBT = 10Γ0. Other parameters are the same as those of Figure 1.

In Figures 1, 2, and 3 we have considered the case with EF below QD energy levels. It would be interesting to investigate the case with EF above the energy levels of QDs. Figure 4 shows Ge, S, κe, and ZT of an SCQD with tc = 3Γ0 as functions of applied gate voltage for various temperatures. Once tc>(ΓL + ΓR) = 2Γ0, the eight peaks for Ge can be resolved at kBT = 1Γ0. These eight peaks correspond to the following resonant channels: Etc, E + tc, E + U,jtc, E + U,j + tc, E + U,j + U tc, E + U,j + U + tc, E + 2U,j + U tc, and E + 2U,j + U + tc, which are tuned by the gate voltage to be aligned with EF. These eight channels result from the four configurations of p1, p3, p6, and p8 in Equation 4. Such a result implies that SCQD with identical QDs acts as a QD with effective two levels of Etc and E + tc and satisfying Hund’s rule. These eight peaks are smeared out with increasing temperature. The sign changes of S with respect to the gate voltage result from the bipolar effect, i.e., the competition between electrons and holes, where holes are defined as the unoccupied states below EF[13]. The electronic thermal conductance (κe) also exhibits eight peaks, and we noticed that the local maxima of the κe curve nearly coincide with the local minima of the Ge curve. We see that ZT values are still larger than 3 even when Eis deeply below EF(say, at eVg = 70Γ0). This is attributed to the electron Coulomb interaction. To illustrate that, we also show the results with U=U,j = 0 at kBT = 3Γ0(see the curve with triangle marks). The oscillation of ZT in the case of U=U,j = 0 is attributed to the sign change of S at Vg = 10Γ0. Note that S goes to zero at Vg = 10Γ0, which results from the electron-hole symmetry (with E + tc and Etc straddling EF symmetrically). We see that ZT vanishes for eVg ≥ 40Γ0in the absence of electron Coulomb interactions. Unlike the case of EF <E, where the finite U causes reduction of ZT, here, the electron Coulomb interaction leads to enhancement of ZT when EF >E.

thumbnailFigure 4. Function of applied gate voltage. (a) Ge, (b) S, (c) κe, and (d) ZT as a function of applied gate voltage for kBT = 1Γ0(solid), 2Γ0(dashed), and 3Γ0(dotted). E= EF + 10Γ0 and tc = 3Γ0. Other parameters are the same as those of Figure 1. The curves with triangle marks are for the case without electron Coulomb interactions for kBT = 3Γ0.


In summary, the thermoelectric properties including Ge, S, κe, and ZT of the SCQD junction system are investigated theoretically. We demonstrate that the Carnot efficiency can be reached when tc approaches zero in the absence of phonon thermal conductance. When the phonon contribution dominates the thermal conductance of the SCQD junction, the optimization of ZT can be obtained by the thermal power defined as S2Ge. We also found that the presence of electron Coulomb interactions can lead to either reduction or enhancement of ZT, depending on whether the Fermi level is below or above the QD level.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DMTK and Y-CC established the theoretical formalism. DMTK carried out the numerical calculations and drafted the manuscript. Y-CC conceived this study and participated in its coordination. All authors read and approved the final manuscript.

Authors’ information

DMTK received his Ph.D. degree from National Taiwan University in 1996. He joined the Department of Electrical Engineering, National Central University in 2003 as an assistant professor, became an associate professor in 2005, and professor in 2008. His main research interests include nanodevices and quantum transport.

Y-CC received his Ph.D. degree from the California Institute of Technology in 1980. He joined the Physics Department, University of Illinois at Urbana-Champaign in 1980 as a visiting research assistant professor and became an assistant professor in 1982, associate professor in 1986, and professor in 1991. In 2005, he joined Academia Sinica, Taiwan as a Distinguished Research Fellow of the Research Center for Applied Sciences. His main research interests include condensed matter theory, semiconductor electronics, photonic materials, and optoelectronic devices.


This work was supported in part by the National Science Council, Taiwan under contract nos. NSC 99-2112-M-008-018-MY2 and NSC 98-2112-M-001-022-MY3.


  1. Minnich AJ, Dresselhaus MS, Ren ZF, Chen G: Bulk nanostructured thermoelectric materials: current research and future prospects.

    Energy Environ Sci 2009, 2:466-479. Publisher Full Text OpenURL

  2. Mahan G, Sales B, Sharp J: Thermoelectric materials: new approaches to an old problem.

    Physics Today 1997, 50(3):42-47. Publisher Full Text OpenURL

  3. Venkatasubramanian R, Siivola E, Colpitts T, O’Quinn B: Thin-film thermoelectric devices with high room-temperature figures of merit.

    Nature 2001, 413:597-602. PubMed Abstract | Publisher Full Text OpenURL

  4. Boukai AI, Bunimovich Y, Tahir-Kheli J, Yu JK, Goddard IIIWA, Heath JR: Silicon nanowires as efficient thermoelectric materials.

    Nature 2008, 451:168-171. PubMed Abstract | Publisher Full Text OpenURL

  5. Harman TC, Taylor PJ, Walsh MP, LaForge BE: Quantum dot superlattice thermoelectric materials and devices.

    Science 2002, 297:2229-2232. PubMed Abstract | Publisher Full Text OpenURL

  6. Hsu KF, Loo S, Guo F, Chen W, Dyck JS, Uher C, Hogan T, Polychroniadis EK, Kanatzidis MG: Cubic AgPbmSbTe2+m: bulk thermoelectric materials with high figure of merit.

    Science 2004, 303:818-821. PubMed Abstract | Publisher Full Text OpenURL

  7. Murphy P, Mukerjee S, Moore J: Optimal thermoelectric figure of merit of a molecular junction.

    Phys Rev B 2008, 78:161406-161410. OpenURL

  8. Kuo DMT, Chang YC: Thermoelectric and thermal rectification properties of quantum dot junctions.

    Phys Rev B 2010, 81:205321-205331. OpenURL

  9. Dubi Y, Di Ventra M: Heat flow and thermoelectricity in atomic and molecular junctions.

    Rev Modern Phys 2011, 83:131-155. Publisher Full Text OpenURL

  10. Ono K, Austing DG, Tokura Y, Tarucha S: Current rectification by Pauli exclusion in a weakly coupled double quantum dot system.

    Science 2002, 297:1313-1317. PubMed Abstract | Publisher Full Text OpenURL

  11. Fransson J, Rasander M: Pauli spin blockade in weakly coupled double quantum dots.

    Phys Rev B 2006, 73:205333-205342. OpenURL

  12. Sun QF, Xing Y, Shen SQ: Double quantum dot as detector of spin bias.

    Phys Rev B 2008, 77:195313. OpenURL

  13. Kuo DMT, Shiau SY, Chang YC: Theory of spin blockade, charge ratchet effect, and thermoelectrical behavior in serially coupled quantum dot system.

    Phys Rev B 2011, 84:245303-245314. OpenURL

  14. Kuo DMT, Chang YC: Tunneling current spectroscopy of a nanostructure junction involving multiple energy levels.

    Phys Rev Lett 2007, 99:086803-086807. PubMed Abstract | Publisher Full Text OpenURL

  15. Chang YC, Kuo DMT: Theory of charge transport in a quantum dot tunnel junction with multiple energy levels.

    Phys Rev B 2008, 77:245412-245428. OpenURL

  16. Markussen T, Jauho AP, Brandbyge M: Surface-decorated silicon nanowires: a route to high-ZT thermoelectrics.

    Phys Rev Lett 2009, 103:055502-055506. PubMed Abstract | Publisher Full Text OpenURL

  17. Santamore DH, Cross MC: Effect of phonon scattering by surface roughness on the universal thermal conductance.

    Phys Rev Lett 2001, 87:115502-115506. PubMed Abstract | Publisher Full Text OpenURL

  18. Rego LGC, Kirczenow G: Quantized thermal conductance of dielectric quantum wires.

    Phys Rev Lett 1998, 81:232-236. Publisher Full Text OpenURL

  19. Schwab K, Henriksen EA, Worlock JM, Roukes ML: Measurement of the quantum of thermal conductance.

    Nature 2000, 404:974-977. PubMed Abstract | Publisher Full Text OpenURL