arxiv: v3 [cond-mat.mes-hall] 10 Apr 2013

Transcription

1 Efficiency of tree-terminal termoelectric transport under broken-time reversal symmetry Vinita Balacandran Center for Nonlinear and Comple Systems Università degli Studi dell Insubria Via Valleggio Como Italy arxiv: v3 cond-mat.mes-all] 10 Apr 013 Giuliano Benenti and Giulio Casati CNISM & Center for Nonlinear and Comple Systems Università degli Studi dell Insubria Via Valleggio Como Italy and Istituto Nazionale di Fisica Nucleare Sezione di Milano via Celoria Milano Italy (Dated: February 1 018) We investigate termoelectric efficiency of systems wit broken time reversal symmetry under a tree-terminal transport. Using a model of Aaronov-Bom interferometer formed wit tree noninteracting quantum dots we sow tat Carnot efficiency η C can be acieved wen te termopower is a symmetric function of te applied magnetic field. On te oter and te imal value of te efficiency at imum power is obtained for asymmetric termopower. Indeed we sow tat Curzon-Alborn limit is eceeded witin te linear response regime in our model. Moreover we investigate termoelectric efficiency for random Hamiltonians drawn from te Gaussian Unitary Ensemble and for a more abstract transmission model. In tis latter model we find tat te efficiency is improved using sarp energy-dependent transmission functions. PACS numbers: 7.0.Pa Ln I. INTRODUCTION Termoelectrics convert temperature gradients into electric voltages and vice versa. Strong demand for a cost-effective pollution-free form of energy conversion as resulted in a pletora of termoelectric-based applications. Unfortunately te practical usage is limited by etremely low performance of termoelectric materials. Increasing te efficiency of termoelectric materials is one of te main lines of current termoelectric researc 1 7. In te linear response regime te performance of termoelectric materials is caracterized by a single dimensionless parameter called figure of merit ZT wic is a combination of te main transport properties of a material i.e. te electric conductivity σ te termal conductivity κ and te termopower S as well as of te absolute temperature T: ZT = (σs /κ)t. Te imum efficiency is given by η = η C ZT +1 1 ZT (1) Carnot efficiency η C is reaced in te limit of ZT. Altoug termodynamics does not impose any upper bound on ZT te strong interdependence between te electric and termal transport properties makes it etremely ard to increase te value of ZT above 1. To compete wit te eisting mecanical based energy conversion applications development of termoelectric materialswitzt atleastabove3isrequired. Anattractive alternative to increase te performance is to use broken time reversal symmetry systems for wic te efficiency is determined by two parameters an asymmetry parameter anda figureofmerit y generalizingzt 8. Intesesystems imum efficiency η can be as ig as Carnot efficiency even wit very low values of figure of merit y provided te asymmetry parameter is very large. Note tat in time reversal symmetric systems = 1 and y = ZT. In a termodynamic system Carnot efficiency is obtained for fully reversible transformations. Tis requires quasistatic transformations and consequently te derived output power is zero. Hence te notion of efficiency at imum power was introduced and in te linear response regime is given by 9 η(ω ) = η C ZT ZT +. () Note tat in te limit of ZT η(ω ) takes te imum value of η C /. Tis upper bound is commonly referred to as te Curzon-Alborn limit In principle for broken time reversal symmetry systems tis limit can be eceeded wit values of te asymmetry parameter suc tat > 1 8. Moreover η(ω ) η C wen (always witin te linear response regime). Hence it is potentially of practical relevance to study te termoelectric transport in broken time reversal symmetry systems. Termopower for a broken time reversal symmetry system is in general asymmetric wit respect to te time reversibility breaking parameter B i.e. S(B) S( B). Te asymmetry parameter = S(B)/S( B). In a non-interacting system inelastic scattering can result in asymmetric termopower 01. Conveniently tis can be acieved wit te introduction of noise by means of a tird terminal (probe) wit its temperature and cemical potential adjusted suc tat tere is no net average flu of particles and eat between te terminal and te system. Various aspects of tree-terminal termoelectric transport were investigated in Refs. 6 Large asymmetry in te termopower were obtained in Ref. 0 using

2 an Aaronov-Bom interferometer model first treated in te contet of termoelectric transport in Ref. 5. However te obtained efficiency wasvery low 0. Also following te above lines a classical deterministic tree terminal transport model was studied sowing large asymmetry wit very low efficiency 7. Hence it remains to be seen weter and under wat conditions efficiency close to Carnot can be obtained using broken time reversal symmetry system under a tree terminal transport. In tis paper we investigate te termoelectric efficiency of non-interacting systems wit broken time reversal symmetry under tree terminal transport in te linear response regime. First we consider a model of Aaronov-Bom interferometer formed wit tree dots as analyzed in Ref. 0. In contrast to te earlier work we focus on te efficiency obtained over te global optimization of all system and reservoir parameters following simulated annealing metod 8. Our results sow tat Carnot efficiency η C can be obtained for imum efficiency η wen te termopower is symmetric i.e. = 1. By adding asymmetry η decreases. However te efficiency at imum power η(ω ) is imum wen tere is asymmetry in termopower ( 1). To te best of our knowledge for te first time we sow tat te Curzon-Alborn limit η(ω ) = η C / wic is a rigorous upper bound for systems wit time-reversal symmetry is eceeded in te linear response regime. However bot te efficiency η and η(ω ) decrease wit increase in asymmetry of termopower at large values of. Terefore we study broader classes of models in an attempt to improve efficiency at large asymmetry. We consider random Hamiltonians drawn from Gaussian Unitary Ensemble (GUE) and finally an abstract model of transmission probabilities. Wile in tis latter model we sow tat η(ω ) can take values as ig as 7η C still we could not find after global optimization of all parameters of te model large efficiency at large asymmetry. Our results obtained for very broad classes of non-interacting models implies tat it is practically very ard if not impossible to acieve in suc models large efficiency wit large asymmetry in termopower. Te paper is structured as follows: In Sec. II we review te calculations of efficiency of broken time reversal symmetry systems under tree terminal transport. Dependence of efficiency on asymmetry of termopower for various models wit broken time reversal symmetry is analyzed in Sec. III. Finally Sec. IV summarizes our results. II. MODEL AND METHOD A. General setup Te general set up consists of a system in contact wit two reservoirs left (L) and rigt (R) at temperatures T L = T + TT R = T and cemical potentials µ L = µ + µµ R = µ. Inelastic scattering effects are simulated by means of a tird (probe) reservoir (P) at temperature T P = T + T P and cemical potential µ P = µ+ µ P. Let J ρk and J Ek denote te electric and energy currents from te kt reservoir (k=lrp) into te system wit te steady-state constraints of carge and energy conservation k J ρk = 0 k J Ek = 0. Te sum of te entropy production rates at te reservoirs reads Ṡ = k (J Ek µ k J ρk )/T k. Witin linear response Ṡ = J X 4 i=1 J ix i were J and X are four dimensional vectors defined as J = (ej ρl J ql ej ρp J qp ) (3) X = ( µ et T T µ P et T P ). (4) T Here e is te electron carge and J qk J Ek µj ρk is te eat current. Te relation between te flues J i and te termodynamic forces X i witin linear irreversible termodynamics is J = LX (5) were L is a 4 4 Onsager matri J and X are written as column vectors. Eq. (5) can be written in te block matri form as ( ) ( )( ) Jα Lαα L = αβ Xα (6) J β L βα L ββ X β were α stands for (1) and β for (34). TeprobereservoirisadjustedinsucawaytatJ β = J 3 = J 4 = 0 tat is te net electric and eat flow from te probe into te system vanises. Tis implies tat X β = L 1 ββ L βαx α and J α = L X α L L αα L αβ L 1 ββ L βα. (7) Tus te problem as been reduced to two coupled flues as ( ) ( )( ) J1 L = 11 L 1 X1 J L 1 L (8) X were te reduced Onsager matri L satisfies te Onsager-Casimir relations L ij (B) = L ji ( B) (ij = 1). (9) Here B is a magnetic field breaking te time reversal symmetry. Seebeck and Peltier coefficients are given by S = L 1 /(etl 11 ) and Π = L 1 /(el 11 ). Termopower is asymmetric wen L 1 L 1 i.e. Π ST. Maimum efficiency is given by y +1 1 η = η C (10) y +1+1 and depends on two parameters: te asymmetry parameter and te figure of merit y were L 1 L (11) 1 y L 1 L 1 detl. (1)

3 3 Efficiency at imum power is η(ω ) = η C y +y. (13) Altoug te termodynamics does not impose any restriction on te attainable values of asymmetry parameter te positivity of entropy production rate implies tat () y 0 if < 0 0 y () if > 0 (14) were () = 4/( 1). Maimum values ofbot η and η(ω ) are obtained for a given wen y = (). We denote suc imum values as η and η(ω ) respectively. From Eqs. (10) and (13) it follows tat te teoretical upper bounds are given for imum efficiency by η = η C if < 1 η = η C if 1 (15) and for efficiency at imum power by η(ω ) = η C +1. (16) It is clear from te above equation tat Curzon-Alborn limit can in principle be overcome for broken-time reversal symmetry systems wit asymmetry parameter > 1 8. µ k )/k B T k ] + 1} 1 is te Fermi distribution function. Analogous epressions can be written for J 3 and J 4 provided te terminal L is substituted by P. Te Onsager coefficients L ij are obtained from te linear response epansion of te currents J i as L 11 = e T L 1 = et L = T de k Lτ Lk (E) f(e) ] de(e µ) k Lτ Lk (E) f(e) ] de(e µ) Lk (E) k Lτ f(e) L 1 = L 1. (19) Analogous formulas are obtained for L 33 L 34 = L 43 and L 44 wit tep terminalusedinstead ofl. Similarly te off diagonal block elements are obtained as L 13 = e T L 14 = et L 4 = T deτ LP (E) f(e) de(e µ)τ LP (E) ] f(e) de(e µ) τ LP (E) f(e) ] ] ] L 3 = L 14. (0) B. Noninteracting systems Termoelectric efficiency can be calculated eactly for noninteracting models by means of Landauer-Büttiker formalism 9. Considerten anoninteractingsystemwit Hamiltonian H S. We model te reservoirs as ideal Fermi gases wit Hamiltonian H R = qk E qk d qk d qk. (17) Here E qk is te energy of an electron in te state q in te kt reservoir and d qk and d qk are te corresponding creation and anniilation operators. Te electric and eat currents from te left reservoir are given by J 1 = e de τ kl (E)f L (E) τ Lk (E)f k (E)] k J = 1 de (E µ L )τ kl (E)f L (E) k τ Lk (E)f k (E)] (18) were τ kl (E) is te transmission probability from reservoir l to reservoir k at energy E and f k (E) = {ep(e Using τ PL (E) instead of τ LP (E) in Eq. (0) L 31 L 3 = L 41 and L 4 are obtained. Te transmissionprobabilitiesτ pq (E) arecalculated as τ pq (E) = TrΓ p (E)G(E)Γ q (E)G (E)] (1) were te broadening matrices Γ k are defined in terms of te self-energies Σ k as Γ k iσ k (E) Σ k (E)] and te (retarded) system Green function G(E) E H S k Σ k(e)] 1. Eplicit epression for te matri elements of te reduced Onsager matri L can be derived from Eq. (7). We obtain L 11 = L 11 e di(e E 1 )(E E 3 ) D R 3 L = L 1 di(e 3 µ)(e 1 µ)(e E 1 )(E E 3 ) D R 3 L 1 = L 1 e di(e 3 µ)(e E 1 )(E E 3 ) D R 3 L 1 = L 1 e di(e 1 µ)(e E 1 )(E E 3 ). () D R 3

4 4 Here di and D are 3 di = e T 3 ( 3 de i f ) ] i i=1 τ LP (E 1 )τ PL (E )+τ PR (E )]τ PL (E 3 ) ( ) { ( et D = det(l ββ ) = de i f ) i R i=1 τ PL (E i )+τ PR (E i )]}(E 1 µ)(e 1 E ). (3) t LP P T P P P t RP III. RESULTS A. Aaronov-Bom interferometer L T L L L t LR R R R T R Here we discuss te termoelectric efficiency of a tree quantum dot ring structure pierced by an Aaronov- Bom flu and coupled to tree reservoirs wit eac dot connected independently to one reservoir. Figure 1 displays a sketc of te model. Te system is described by te Hamiltonian H S = k ǫ k c k c k +(t LR c R c Le iφ/3 + t RP c P c Re iφ/3 +t PL c L c Pe iφ/3 +H.c.) (4) were ǫ k are te on-site energies t ij are te opping strengts (ij = LRP) φ is te flu. c k (c k) are te creation (anniilation) operators of te electron in te kt dot. Reservoirs are ideal fermi gases wit Hamiltonian given by Eq. (17). Te dot-reservoir coupling Hamiltonian is H SR = q (t ql c L d ql +t qr c R d qr + t qp c P d qp +H.c.). (5) Here t qj is te tunneling amplitude of an electron in te state q into te kt reservoir. We assume te wide-band limit and ence te broadening matrices are given by Γ k = γ k c k c k were γ k = π q t qk δ(e E qk ). Note tat γ k measures te tunneling rate of electrons between te reservoir k and te system. We follow Landauer-Büttiker formalism and calculate Onsager coefficients from Eqs.(19) to (1). Wen tere is anisotropy in te system (ǫ k ǫ j ) and te Aaronov- Bom flu φ is non-zero te off diagonal elements of te reduced Onsager matri L are asymmetric functions of te flu i.e. L 1 (φ) L 1 (φ) = L 1 ( φ). Te asymmetry parameter defined in Eq. (11) is te ratio of te off diagonal Onsager matri elements and is in general different from te time reversal symmetric case were = 1. For tis model te efficiency is a function of 1 independent parameters: ǫ 1 ǫ ǫ 3 t LR t RP t PL φ γ L γ R γ P T µ. We imize te efficiency over all tese parameters using te simulated annealing metod. For numerical convenience te parameters are restricted to FIG. 1. Scematic picture of Aaronov-Bom interferometer model FIG.. Optimized values of imum efficiency η as a function of te asymmetry parameter for te Aaronov- Bom interferometer discussed in te tet wit a linear panel (a)] and a logaritmic panel (b)] scale for η. Note tat η is scaled in units of η C in tis figure and tereafter. (a) (b)

5 5 te intervals ǫi 5 5] tij 5 5] T ] µ 5 5] φ 1 1] γi ] (in units suc tat e = = kb = 1). Note tat te results are uncanged by varying te range of parameter values by a few times. Fig. (a) sows te dependence of te numerically obtained imum efficiency η on te asymmetry parameter 30. In te absence of asymmetry in te termopower i.e. wen = 1 Carnot efficiency is reaced. As te asymmetry is introduced te optimized imum efficiency is always less tan Carnot efficiency. 0.6 Landauer-Bu ttiker formalism te obtained results are eact. Hence our results are bound only by macine accuracy and te calculated efficiency is accurate up to 15 decimal points. Even toug tere is only a small improvement wit respect to te Curzon-Alborn limit suc result is interesting in tat it provides te first evidence in a concrete model of te fact tat te CurzonAlborn limit wic is a universal upper bound for timereversal systems witin linear response can be eceeded wen time reversibility is broken. According to our numerical optimization large asymmetry results in low termoelectric efficiency. Also we did not observe a significant improvement in efficiency by increasing te number of levels in te system up to si (for systems wit a larger number of levels it becomes difficult to obtain convergence) B. Gaussian Unitary Ensemble model (a) We also plot te dependence of te optimized imum efficiency η on in a lin-log scale in Fig. (b). It is clear from te figure tat for < 1 η increases wit from its zero value at = 0. Also te increase is more drastic wit positive values of. Indeed only for > 0 Carnot efficiency is reaced. For negative values of η increases initially (up to = 10 were η = 03) and decreases tereafter. Similar results were obtained by imizing te efficiency at imum power η(ω ). Tis is illustrated in Fig. 3. For = 1 η(ω ) = ηc / and te CurzonAlborn limit is recovered. From our results we find tat at large asymmetry η(ω ) decreases wit te increase of. It is interesting to remark tat in contrast to η te imum value of η(ω ) is obtained wen 6= 1. In particular we find tat very near to te symmetric value = 1 te Curzon-Alborn limit of efficiency at imum power η(ω ) = ηc / can be eceeded as sown in te inset of Fig. 3. Note tat in FIG. 3. Optimized values of efficiency at imum power η(ω ) versus asymmetry parameter for te AaronovBom interferometer. Te inset sows te variation of η(ω ) near te symmetric value = 1. Note tat te Curzon-Alborn limit η(ω ) = ηc / is eceeded in tis regime. η(ω ) is scaled in units of ηc in tis figure and tereafter (b) FIG. 4. Relation between imum efficiency η and asymmetry parameter panel (a)] and between efficiency at imum power η(ω ) and asymmetry parameter panel (b)] for 0000 realizations of 3 3 Hamiltonian drawn from te Gaussian Unitary Ensemble. To see weter termoelectric efficiency can be improved in more general systems we turn to Random Matri Teory (RMT) models. Indeed Hamiltonians drawn from Gaussian Unitary Ensemble (GUE) give a good description of comple pysical systems wit broken time reversal symmetry. To tis end in tis section we study te efficiency (and asymmetry) distribution of GUE Hamiltonians. We consider a random Hamiltonian HS drawn from a Gaussian Unitary Ensemble31. Tat is matri elements Hlm of te Hamiltonian HS are cosen suc tat te diagonal entries are normally distributed wit zero mean 3

6 6 and unit variance (N(0 1)) and te off-diagonal entries are drawn independently and identically from te normal distribution (subject to being Hermitian) wit mean zero and variance 1 (N(0 1 )). In oter words H lm = U lm + iv lm were U lm V lm N(0 1 ) for 1 l < m n H ml = Hlm and H mm N(01) for 1 m n wit n te dimension of te matri. To compare wit te results of te Aaronov-Bom interferometerwesettedimensionoftematriasn = 3. Also dot-reservoircouplingis takento be sameastat of te Aaronov-Bom interferometer discussed in te previous subsection. Results of one suc calculation wit 0000 different realizations are plotted in Fig. 4. Here we take T = 0.1 µ = 0.1 and γ L = γ R = γ P = 0.1. Note tat results are not sensitive to variations of te parameters T µ γ L γ R γ P. Top panel represents te dependence of imum efficiency η on te asymmetry parameter and bottom panel te dependence of efficiency at imum power η(ω ) on. Results sown in te figure reveal tat as te asymmetry in te termopower is increased te efficiency decreases. Tis implies tat it is practically ard to acieve large efficiency wit large asymmetry in termopower and confirms our results obtained by optimizing te efficiency of Aaronov-Bom interferometer in te previous subsection. Actually we find tat obtaining large values of efficiency is igly improbable. Tis is illustrated in Fig. 5 wit te logaritmic distribution of η and η(ω ). Bot η and η(ω ) follows an eponentially decaying distribution indicating tat large efficiency is eponentially improbable. Also te distribution of asymmetry parameter follows eponential distribution implying tat large asymmetry is also a rare situation in tis model. Note tat te imum value of efficiency is obtained around te symmetric point = 1. Hence to better clarify te role of asymmetry we ave plotted te logaritmic distribution of η and η(ω ) restricted to data for wic 1 > in te inset of Fig. 5(a) and 5(b) respectively. From te figure it is clear tat distribution of η and η(ω ) is eponential even in te absence of symmetry. On te oter and tere is no numerical evidence of a cut-off or a sarp border forbidding η = η C at > 1. In order to study te dependence of efficiency on te number of levels of te system we studied GUE Hamiltonians by increasing te dimension n. Our results wit n = 4 and n = 5 (data not presented) sow tat tere is little improvement in efficiency wit te increase in te levels of te system. Tus aving studied bot te Aaronov-Bom interferometer model and te broader class of GUE Hamiltonians we can conclude tat for Hamiltonian models it is difficult to obtain large values of efficiency wit large asymmetry. Hence in te net subsection we turn to a more abstract transmission model and investigate weter te efficiency can be improved. P P P 10 1 (a) (b) ) (c) FIG. 5. Distribution of (a) imum efficiency η (b) efficiency at imum power η(ω ) and (c) asymmetry parameter. Inset of (a) and (b) sows te distribution of η and η(ω ) taking into account only data suc tat 1 >. Note tat obtaining ig values of bot efficiency and asymmetry is igly improbable. C. Transmission model Here we assume tat te information about te system is inaccessible and only te scattering matri tat defines te transmission probabilities τ ij (E) at eac energy E are known. Onsager matri elements are calculated directly from Eqs. (19) and (0) by substituting te values of transmission probabilities τ ij (E). Tese transmission probabilities are bound to follow i τ ij(e) = 1 tat implies te conservation of probability and te sum rule i τ ij(e) = j τ ij(e) tat ensures tat te currents vanis at equilibrium. To simplify te calculations te transmission probabilities are taken to be constant over a window of energy. More precisely eac energy window is caracterized by two parameters: Ē k te center of te window and k te widt of te window. Tus in our model τ ij (E) = (c ij ) k (0 (c ij ) k 1) wen energy E Ēk k /Ēk + k /] and τ ij (E) = 0 oterwise. Efficiency is optimized over te parameters temperature T cemical potential µ energies Ēk widt k and transmission probabilities (c ij ) k. We consider n non-overlapping transmission windows (k = 1... n).

7 7 First we consider a single energy window of transmission(n = 1). Here teintegrationdomainr 3 ineq. () is restricted to te cube Ē1 1 / E i Ē1 + 1 / (i = 13) and 3 ( L 1 L 1 = A de i f ) ] i were te constant R 3 i=1 (E 3 E 1 )(E E 1 )(E E 3 ) (6) A = e3 T 3 3 D (c LP) 1 (c PL ) 1 +(c PR ) 1 ](c PL ) 1. (7) Te integral in Eq. (6) vanises since it is an odd function of for instance E 1 E 3 and te integration domain is symmetric under ecange of E 1 and E 3. Tus L 1 = L 1 and a symmetric termopower is obtained. For n = energy windows we found numerically tat te asymmetry parameter is always limited to a finite interval. Infinitely large values of asymmetry are obtained for transmission wit at least tree energy windows. For n = 3 windows we need to optimize over 0 parameters: T µ Ē k k and (c LR ) k (c LP ) k (c PL ) k (c PR ) k (k = 13). (Te oter transmissions (c ij ) k are ten determined from te conditions i (c ij) k = j (c ij) k = 1.) f ' (E) PL (E) LP (E) RL (E) LR (E) RP (E) PR (E) (a) (b) (c) (d) (e) (f) (g) -4-0 E FIG. 7. Optimized transmission functions for = 1.4 using te transmission model. Panel (a) sows f wit f derivative of te Fermi distribution function (f (E) = f(e) ) as a function of energy E wit optimal temperature T = 1 and cemical potential µ = Te transmission probabilities τ LR τ RL τ LP τ PL τ RP and τ PR at eac energy E are plotted in panels (b) (c) (d) (e) (f) and (g) respectively. Note tat optimal transmission probabilities are eiter 0 or 1. FIG. 6. Maimum efficiency η versus asymmetry parameter under tree terminal transport wit transmission model described in te tet. Solid curve corresponds to te teoretically predicted upperboundη (given byeq. (15)) for η for a generic broken time reversal symmetry system. Te optimized efficiency saturates te teoretical upper bound only near te symmetric value = 1. Fig. 6 sows te variation of te optimized imum efficiency η wit te asymmetry parameter. As in te case of te Aaronov-Bom interferometer Carnot efficiency η C is obtained only in te symmetric case = 1. Introducing asymmetry results in te reduction of imum efficiency from η C. In particular for positive values of asymmetry η increases quadratically from zero value to η C at = 1 and decreases tereafter. Wen te asymmetry parameter is negative η increases wit from = 0 to = 50 (were η = 46) and ten decreases for larger values of.

8 8 ) FIG. 8. Efficiency at imum power η(ω ) versus asymmetry parameter under tree terminal transport using te transmission model described in te tet. Teoretically predicted upper bound η(ω ) (given by Eq. (16)) of η(ω ) is plotted as solid curve. Dotted dased line corresponds to te Curzon-Alborn limit η(ω ) = η C/. Note tat te Curzon-Alborn limit is eceeded in te interval 1] wit te transmission model. Compared to te Aaronov-Bom interferometer model te rate of variation of imum efficiency wit asymmetry is muc less for te transmission model. Also te values of te optimized efficiency wit large asymmetry are greater tan in te Aaronov-Bom interferometer model by more tan one order of magnitude. For instance wen te asymmetry parameter = 10 te optimized value of η = 045η C for te Aaronov- Bom interferometer model and 8η C for te transmission model. To understand te large difference of efficiency obtained between te two models we take a closer look on te optimized values of te parameters. Our analysis sows tat optimal values of η are obtained wen te transmission probabilities (c ij ) k wic in principle could take any value between 0 and 1 are eiter 0 or 1. In particular we find tat te efficiency is imal for particular combinations of te transmission probabilities as described below. Out of te tree combinations (one for eac window) two correspond to symmetric transmission ((c ij ) k = (c ji ) k for ij = LRP) and one to perfectly asymmetric transmission ((c ij ) k = 1 and (c ji ) k = 0 for some values of i and j). Also for te symmetric cases only two reservoirs are involved in te transmission. For instance te optimized values of parameters for = 1.4 are plotted in Fig. 7. Panel (a) corresponds to te derivative of Fermi distribution function f(e) wit te optimized temperature T = 1 and cemical potential µ = Tree energy windows of transmission are 1) ] ) ] and3)93101]. Forte firstcase τ LR = τ RL = 0 τ LP = τ PL = 0 and τ RP = τ PR = 1. Here te transmission is perfectly symmetric and is only between te rigt reservoir and te probe (i.e. R P). Transmission is perfectly asymmetric in te second energy window. Tis is clear from te figure as τ LR = 0τ RL = 1 τ LP = 1τ PL = 0 and τ RP = 0τ PR = 1. Transport ofelectrons between te reservoirsis from L R P L. For te tird energy window te transmission probabilities are τ LR = τ RL = 0 τ LP = τ PL = 1 and τ RP = τ PR = 0 and are completely symmetric between te left reservoir and te probe (i.e. L P). We also obtained similar results using models wit tree delta peaks for transmission instead of transmission windows. Note tat for te delta peaks model te widt of transmission windows tends to zero i.e. k 0. Curzon-Alborn limit is sligtly eceeded for small values of asymmetry in te Aaronov-Bom interferometer model. Can te efficiency at imum power go significantly beyond te Curzon-Alborn limit? To ceck tis we optimize efficiency at imum power η(ω ). Results of our optimization are sown in Fig. 8. As epected Curzon-Alborn limit η(ω ) = η C / is eceeded. Wen > 0 η(ω ) increases quadratically wit increase in te asymmetry parameter until = At = it takes te imum value of 69η C. Beyond η(ω ) decreases. Moreoverfor te asymmetryparameterin te interval1] η(ω ) > η C /. Our numerical analysis sows tat te optimized values of transmission probabilities corresponding to te above interval is completely symmetric for two energy windows and completely asymmetric for te tird window as discussed in te Fig. 7. Also outside te interval te optimized values of te probabilities are completely symmetric for two windows wereas for te tird window it is not completely asymmetric between te different reservoirs. It is clear from te definition of transmission probabilities in Eq. (1) tat Aaronov-Bom interferometer model wit tree quantum dots discussed in te previous section cannot realize te above mentioned transmission probabilities tat leads to η(ω ) largely above te Curzon-Alborn limit. A more complicated Hamiltonian wit large number of levels migt approac a one to one correspondence between te results in te transmission and Aaronov- Bom interferometer model. Similarly to imum efficiency at large values of te efficiency at imum power decreases wit. Also for < 0 η(ω ) increases initially before te decrease (wit te imum value η(ω ) = 461 obtainedat = 50). Notetatforlargevaluesofasymmetry bot te imum efficiency and te efficiency at te imum power are muc smaller tan te upper bounds set by termodynamics (as discussed in Ref. 8 ) and sown as solid curves in te Figs. 6 and 8. We ave also increased te number of transmission

9 9 windows up to n = 7 (data not sown) and found only a sligt increase of te imum value of η(ω ) up to 709η C at n = 7. Our numerical results suggest tat a tree terminal termoelectric transport yields large asymmetry wit very low efficiency. Does te termodynamics impose bound on imum value of efficiency in non-interacting systems wit broken time reversal symmetry under tree terminal transport? Is it possible to acieve Carnot efficiency wit asymmetry in te system under tree terminal transport? Tese questions are addressed in te net subsection by considering a general model wit random values for te Onsager coefficients. D. Random Onsager matri model Consider ten a 4 4 random Onsager matri L. Te elements L ij are cosen from te uniform distribution 11]andaredirectlysubstituted ineq. (7)tocalculate te reduced Onsager matri L. In general te elements of te Onsager matri are not independent and are bound to fulfill certain conditions. First of all it follows from termodynamics tat te entropy production rate as to bepositivei.e. Ṡ = J 1X 1 +J X +J 3 X 3 +J 4 X 4 0. Tis requires tat te Onsager matri L is positive-definite. Now we turn our attention to specific case of noninteracting tree-terminal system. In tis case tere are furter restrictions on te elements of Onsager matri L. For instance it is clear from te Eqs. (19) and (0) tat L 13 L 4 L 31 L 4 0 (and tat as epected in general all diagonal elements L ii 0 i = 1...4). Also te off diagonal elements of eac block matri are even functions of magnetic field and ence are equal i.e. L 1 = L 1 L 14 = L 3 L 3 = L 41 L 34 = L 43. Note tat tis symmetry of te off diagonal elements is broken for te reduced Onsager matri L were we can ave L 1 L 1. Positivity of te transmission probabilities τ ji (E) 0 and of f/ and te sum rule i τ ij(e) = j τ ij(e) imply L 11 L 13 L 31 ; L L 4 L 4 ; L 33 L 13 L 31 ; and L 44 L 4 L 4. Tere are furter constraints on te values of Onsager matri elements from te structure of reservoirs wic we assume to be ideal Fermi gases at temperature T and cemical potential µ. Upper bound on te energy integrals appearing in te Onsager coefficients are ten given by et T e T de(e µ)τ ij (E) de(e µ) τ ij (E) deτ ij (E) f(e) ] T f(e) ] T ln f(e) ] T3 π 4. (8) Te first and te tird bound are saturated by setting τ ij (E) = 1 for all values of energy te second one by setting τ ij = 0 for E < µ and τ ij = 1 oterwise or vice versa τ ij = 0 for E > µ and τ ij = 0 oterwise. Even wit all tese restrictions we find tat a random Onsager matri model can saturate te upper bounds η and η(ω ) from termodynamics for efficiencies (solid curves in Figs. 6 and 8). In particular η(ω ) η C wen. However it is clear from te definition of Onsager coefficients in Eqs. (19) and (0) tat for a non-interacting system tere are furter correlations between te different Onsager coefficients. For instance te different integrands are weiged by (E µ) n n = 01. For small systems in particular for a tree level system tese correlations are significant and put bounds on imum acievable efficiency as implied by our optimization results. IV. CONCLUSIONS We ave investigated te termoelectric efficiency of broken time reversal symmetry systems under tree terminal transport. In tese systems te efficiency is determined by two parameters te asymmetry parameter and te figure of merit y. First we ave studied te optimized efficiency of a realistic model of Aaronov- Bom interferometer formed wit tree non-interacting quantum dots using simulated annealing. Our results sow tat Carnot efficiency η C can be obtained wen te termopower is symmetric ( = 1). Introducing asymmetry in termopower te imum efficiency η decreases from η C. However te efficiency at imum powerη(ω ) is imal wit asymmetry in termopower ( 1). In particular our studies illustrate tat Curzon-Alborn limit can be eceeded in te linear response regime for a realistic model wit broken time reversal symmetry. We note tat our results could be of eperimental relevance in view of te recent progress in te pase-coerent manipulation of eat in solid-state nanocircuits see 6333 and references terein. We ave also studied te termoelectric efficiency of a generic model using random Hamiltonians drawn from GUE. Our analysis sows tat it is igly improbable to obtain large values of bot asymmetry in termopower and efficiency. Furtermore we ave sown tat te efficiency can be improved using an energy dependent transmission. In particular optimizing te transmission matri elements attreeenergywindowsweavefound moretan oneorder of magnitude increase in te efficiency at large asymmetry over te tree level Aaronov-Bom interferometer model. Te optimal values of transmission probabilities at eac energy window are eiter 0 or 1. Also two energy windows correspond to symmetric transmission and one to perfectly asymmetric transmission. Te Curzon-Alborn limit is eceeded in te interval 1] of te asymmetry parameter and η(ω ) as large as 7η C is obtained. Similarly to te Hamiltonian models considered in tis paper te efficiency decreases at large values of asymmetry.

10 10 Our etensive and accurate numerical results suggest tat a tree terminal termoelectric transport is viable only for large asymmetry wit very low efficiency. On te oter and using a model wit random values of Onsager coefficients one may obtain Carnot efficiency for imum efficiency η and efficiency at imum power η(ω ) for arbitrary large values of asymmetry. However we argue tat for non-interacting systems wit a small number of levels tere are correlations between te Onsager coefficients tat bound te efficiency. In order to obtain large efficiency for large asymmetry we need to turn our attention to systems wit more tan tree terminals or to interacting systems or to go beyond te linear response. Tis remains to be analyzed in future. Note. After completion of our work we became aware of a related work 34 sowing te eistence of upper bounds on termodynamic efficiencies for tree-terminal transport as a consequence of te unitarity of te scattering matri. Te results from our transmission model depicted in Figs. 6 and 8 in practice saturate tese new bounds. 1 G. Maan B. Sales and J. Sarp Pys. Today 50 4 (1997). A. Majumdar Science (004). 3 M. S. Dresselaus G. Cen M. Y. Tang R. G. Yang H. Lee D. Z. Wang Z. F. Ren J. -P. Fleurial and P. Gogna Adv. Mater (007). 4 G. J. Snyder and E. S. Toberer Nature Mater (008). 5 A. Sakouri Annu. Rev. Mater. Res (011). 6 Y. Dubi and M. Di Ventra Rev. Mod. Pys (011). 7 G. Benenti and G. Casati Pil. Trans. R. Soc. A (011). 8 G. Benenti K. Saito and G. Casati Pys. Rev. Lett (011). 9 C. Van den Broeck Pys. Rev. Lett (005). 10 J. Yvon Proceedings of te International Conference on Peaceful Uses of Atomic Energy (United Nations Geneva 1955) p P. Cambadal Les Centrales Nucléaires (Armand Colin Paris 1957). 1 I. I. Novikov J. Nucl. Energy 7 15 (1958). 13 F. Curzon and B. Alborn Am. J. Pys 43 (1975). 14 M. Esposito K. Lindenberg and C. Vanden Broeck Pys. Rev. Lett (009). 15 B. Gaveau M. Moreau and L.S. Sculman Pys. Rev. Lett (010). 16 M. Esposito R. Kawai K. Lindenberg and C. Van den Broeck Pys. Rev. Lett (010). 17 N. Nakpatomkun H. Q. Xu and H. Linke Pys. Rev. B (010). 18 U. Seifert Pys. Rev. Lett (011). 19 Y. Apertet H. Ouerdane C. Goupil and P. Lecoeur Pys. Rev. E ibid (01). 0 K. Saito G. Benenti G. Casati and T. Prosen Pys. Rev. B (R) (011). 1 D. Sáncez and L. Serra Pys. Rev. B (R) (011). O. Entin-Wolman Y. Imry and A. Aarony Pys. Rev. B (010). 3 R. Sáncez and M. Büttiker Pys. Rev. B (011). 4 J.-H. Jiang O. Entin-Wolman and Y. Imry Pys. Rev. B (01). 5 O. Entin-Wolman and A. Aarony Pys. Rev. B (01). 6 B. Sotmann R. Sáncez A. N. Jordan and M. Büttiker Pys. Rev. B (01). 7 M. Horvat T. Prosen G. Benenti and G. Casati Pys. Rev. E (01). 8 W. H. Press B. P. Flannery S. A. Teukolsky and W. T. Vetterling Numerical Recipes in Fortran 77: Te Art of Scientific Computing (Cambridge University Press Cambridge 199). 9 S. Datta Electronic Transport in Mesoscopic Systems (Cambridge University Press Cambridge 1995). 30 In te optimization η at a particular value of c of te asymmetry parameter is obtained by subtracting a cost function c from te optimizing function η. 31 See for instance F. Haake Quantum Signatures of Caos nd. ed. (Springer-Verlag Berlin 000). 3 F. Giazotto T. T. Heikkilä A. Luukanen A. M. Savin and J. P. Pekola Rev. Mod. Pys (006). 33 F. Giazotto and M. J. Martínez-Pérez Nature (01). 34 K. Brandner K. Saito and U. Seifert Pys. Rev. Lett (013).