## Nonlinear Seebeck effect of SU( N ) Kondo impurity

D. B. Karki 1, 2 and Mikhail N. Kiselev 2

1 International School for Advanced Studies (SISSA), Via Bonomea 265, 34136 Trieste, Italy 2 The Abdus Salam International Centre for Theoretical Physics (ICTP) Strada Costiera 11, I-34151 Trieste, Italy

We develop a theoretical framework to study the influences of coupling asymmetry on the thermoelectrics of a strongly coupled SU( N ) Kondo impurity based on a local Fermi liquid theory. Applying non-equilibrium Keldysh formalism, we investigate charge current driven by the voltage bias and temperature gradient in the strong coupling regime of an asymmetrically coupled SU( N ) quantum impurity. The thermoelectric characterizations are made via non-linear Seebeck effects. We demonstrate that the beyond particle-hole (PH) symmetric SU( N ) Kondo variants are highly desirable with respect to the corresponding PH symmetric setups in order to have significantly improved thermoelectric performance. The greatly enhanced Seebeck coefficients by tailoring the coupling asymmetry of beyond PH symmetric SU( N ) Kondo effects are explored. Apart from presenting the analytical expressions of asymmetry dependent transport coefficients for general SU( N ) Kondo effects, we make a close connection of our findings with the experimentally studied SU(2) and SU(4) Kondo effects in quantum dot nano structures. Seebeck effects associated with the theoretically proposed SU(3) Kondo effects are discussed in detail.

## I. INTRODUCTION

The increasing demand of quantum technologies for energy harvesting has attracted growing attention towards the necessity for the nano material based energy converters [1]. The presence of quantization effects in nano scale systems allows the controllable comprehension and subsequent control of underlying transport processes [2]. In addition, nano scale systems can offer greatly enhanced thermoelectric response with respect to conventional bulk counterparts [1-3]. These properties of nano scale systems has rekindled the field of thermoelectricity [4]. Over the past years, several experiments has resulted in exciting thermoelectric measurement for nano scale systems, such as quantum dots (QDs), carbon nano-tubes (CNTs), quantum point contacts (QPCs), etc [4-6]. The rapid progress of nano technology has allowed the fine tuning of nano scale transport process, nonetheless, complete understanding of electron interactions in such small scale remains the most challenging task [7].

the SU(2) symmetry group. However, the conventional SU(2) Kondo effects, being protected by particle-hole (PH) symmetry, offer vanishingly small thermoelectric conversion [7]. To achieve appreciable thermopower, the occupation factor of the quantum impurity should be integer, while the PH symmetry should be lifted [16]. The SU( N ) Kondo model with integer occupancy m offers the possibility of avoiding half-filled regime so as to achieve the enhanced thermoelectric production over the conventional SU(2) Kondo correlated systems [1618].

A generic nano device consists of a quantum impurity with intrinsic spin S which is tunnel coupled to two electron reservoirs, the source and the drain. The low energy transport processes are then controlled by the strong interaction between localized spin S and itinerant electrons in the reservoirs. The spin S=1/2 impurity interacting with a single orbital channel of conduction electrons forms a fully screened ground state resulting in quasiparticle resonances at the Fermi level. This paradigmatic screening phenomenon is termed as Kondo effect [8] which is characterized by a low energy scale T K , the Kondo temperature. The many-body Kondo resonance at the Fermi level opens an effective path towards the enhancement of thermoelectric production at the nano scale level [9]. Recent experiments [10-15] have further expanded the scope of transport measurements in Kondo correlated nano scale systems. Most of these studies have been focused on the transport measurement for the spin S=1/2 Kondo impurity described by

The orbital degeneracy of the quantum impurity combines with the true spin symmetry to form the Kondo effect described by higher symmetry group SU( N ). Here the occupancy factor m takes all possible values starting from 1 to N -1. The paradigmatic SU(4) Kondo physics has been experimentally studied in CNTs [12, 19-23], double QDs [24] and single-atom transistor [25]. Various theoretical works [26-31] have contributed towards better understanding of SU(4) Kondo physics over past years. In addition, exciting proposals have been put forth for the experimental realization of different variants of SU( N ) Kondo systems. Possible realization of SU(3) Kondo effects using triple QDs with three and four edge states of the quantum Hall effects was suggested in Ref. [32], which been verified recently using numerical renormalization group study [33]. The proposals for the solid-state realization of SU(6) [34] and SU(12) [35] Kondo effects have likewise attracted considerable attentions both theoretically and experimentally. Beside obtaining the solid-state realization of these exotic SU( N ) Kondo effects, an increasing effort has been put in their cold atomic realization [36-39].

Most of the previous studies on SU( N ) Kondo effects have been focused solely on charge current measurements. However, thermoelectric characterization in a generic nano device usually involves the Seebeck effects. To the best of our knowledge very few studies have tried to uncover the thermoelectric measurements of Kondo effects described by higher symmetry group. The See-

beck effects with a SU(4) Kondo effects have been studied in Ref. [17] and a general theoretical framework for thermoelectric transport of a SU( N ) Kondo model has been developed recently in Ref. [16]. These studies are limited to the setup with perfectly symmetrical tunnel coupling, which is very rarely the case of an experiment. In fact, junction asymmetry could provide important informations about the underlying many body effects [40], thus, has to be taken into account in any calculation to compare its result with the experimental data [41]. Therefore unveiling the effects associated with the coupling asymmetry towards the thermoelectric characterization, the Seebeck effects, of a SU( N ) Kondo effects has remained a challenging problem for many years. In this contribution we develop a theoretical framework based on a local Fermi-liquid theory in combination with the out of equilibrium Keldysh approach to study the influences of coupling asymmetry on the thermoelectric transport of a strongly coupled SU( N ) Kondo impurity.

The paper is organized as follows. In Sec. II we discuss in detail on the formulation of the model, for thermoelectric transport calculations at the strong-coupling regime of SU( N ) Kondo effects, capturing the effects of coupling asymmetry and arbitrary temperature and chemical potentials of the electron reservoirs. We outline the charge current calculations for a SU( N ) Kondo impurity which account for both elastic and inelastic effects using the nonequilibrium Keldysh formalism in Sec. III. The Sec. IV is devoted to the summary of our results for the thermoelectric transport coefficients of SU( N ) Kondo correlated systems characterizing the non-linear Seebeck effects. In this section, apart from presenting the analytical expressions of coupling asymmetry dependent transport coefficients for general SU(N ) Kondo effects, we made a separate analysis of thermoelectrics with i) experimentally studied SU(2) and SU(4) Kondo effects and ii) theoretically proposed SU(3) Kondo effects. The last section V contains the conclusion of our work together with the possible future research plans based on the present work.

## II. MODEL DESCRIPTION

We consider a quantum impurity tunnel coupled to two conducting reservoirs as shown in Fig. 1. The impurity possess N -fold degeneracy by combining the spin and other degrees of freedom, such as the orbital degeneracy. In addition, there are N -species (orbitals) of electrons in both left (L) and right (R) reservoirs. The rotation of the reservoir's electrons is then described by the SU( N ) transformation. Therefore, to describe our system we start form SU( N ) impurity Anderson model [42, 43],

$$H = \sum _ { k, r } \varepsilon _ { k } \left [ c ^ { \dagger } _ { \mathbb { L }, k r } c _ { \mathbb { L }, k r } + c ^ { \dagger } _ { \mathbb { R }, k r } c _ { \mathbb { R }, k r } \right ] + H _ { \text{imp} } + H _ { \text{tun} }. \ ( 1 ) \quad \sin _ { \substack { t _ { \mathbb { L } } = t \complement \\..... } } \text{we e}$$

Here we introduce the notation 'r' to represent the orbital index that takes all possible values starting from 1 to N . The operator c † γ,k r creates an electron with

FIG. 1. Upper panel: Schematic representation of experimental setup for investigating Seebeck effect in nanostructures, where a SU( N ) quantum impurity is sandwiched between two conducting reservoirs. The left (red) and right (blue) reservoirs are in thermal equilibrium, separately, at temperature T L and T R respectively. The tunneling-matrix elements from the impurity to the left/right reservoirs are characterized by t L = cos t θ and t R = sin t θ with θ ∈ (0 , π/ 2). Lower panel: The asymmetry of the tunneling junction is accounted for by introducing a parameter C ≡ (Γ L -Γ R ) / (Γ L + Γ R )=cos2 θ with Γ L R / = πρ res | t L R / | 2 , ρ res being the density of states of the reservoirs. The magenta line represents the variation of asymmetry parameter C with respect to the asymmetry angle θ . We choose the Fermi level in such a way that the chemical potentials of the left and right reservoirs take some specific values µ L R / = ± e ∆ V 2 (1 ∓C ). This choice of chemical potentials amounts to greatly simplify the calculation of charge and heat current through a strongly coupled Kondo impurity (see text for detail). These chemical potentials are represented by the red and blue curve respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00415v1-figure-8b73b885-5d9e-4da9-8dfd-312ad6c8165e.jpg)

momentum k in r-th orbital of the γ (=L R) reservoir. , The energy of conduction electrons ε k is measured with respect to the chemical potential µ . The second term of Eq. (1) represents the Hamiltonian of the impurity possessing N degenerate flavors with single energy level ε d . Then we write the impurity Hamiltonian as

$$H _ { \text{imp} } = \varepsilon _ { d } \sum _ { r } d _ { r } ^ { \dagger } d _ { r } + U \sum _ { r < r ^ { \prime } } d _ { r } ^ { \dagger } d _ { r } d _ { r ^ { \prime } } ^ { \dagger } d _ { r ^ { \prime } }, \quad ( 2 )$$

where d † r is the electron creation operator of the impurity and U represents the charging energy which is assumed to be the largest energy scale of the model. The tunneling processes from the impurity to the reservoirs are accounted for by the very last term of the Eq. (1),

$$H _ { \text{tun} } = \sum _ { k, r } \left ( t _ { L } c _ { L, k r } ^ { \dagger } + t _ { R } c _ { R, k r } ^ { \dagger } \right ) d _ { r } + \text{H.c.}, \quad ( 3 )$$

We explicitly assume the tunneling asymmetry by assigning the tunneling-matrix elements t γ such that t L = cos t θ and t R = sin t θ with θ ∈ (0 , π/ 2). Then the intrinsic total local level width associated with the tunneling is given by Γ γ = πρ res | t γ | 2 with ρ res being the density of states of the reservoirs. For the sake of clarity, we introduce the parameter C ≡ (Γ L -Γ R ) / (Γ L +Γ )=cos2 R θ

to characterize the asymmetry of the tunneling junction. This asymmetry further appears in the Glazman-Raikh rotation [44] of Eq. (1) in the basis of reservoirs electrons

$$\left ( \begin{array} { c } b _ { k r } \\ a _ { k r } \end{array} \right ) = \left ( \begin{array} { c c } \cos \theta & \sin \theta \\ \sin \theta & - \cos \theta \end{array} \right ) \left ( \begin{array} { c } c _ { L, k r } \\ c _ { R, k r } \end{array} \right ). \quad ( 4 ) \quad \text{$\stackrel{ \lambda \circ \\ \text{cefl} } \\ \text{reset}. \end{array}$$

Note that the transformation Eq. (16) effectively decouples the operators a k r from the impurity degrees of freedom. Here we consider the general case of having arbitrary number of electrons m =1 2 , , · · · , N -1 in the impurity. Therefore, the specific choice of impurity level ε d = U (1 -m m/N -) provides the fundamental representation with ∑ r d d † r r ≡ ∑ r n r = m . We then perform the Schrieffer-Wolff transformation [45] followed by the rotation Eq. (16) of the Hamiltonian Eq. (1) to project out the charge states, which results in

$$\mathcal { H } = \sum _ { k, r } \varepsilon _ { k } \left ( a _ { k r } ^ { \dagger } a _ { k r } + b _ { k r } ^ { \dagger } b _ { k r } \right ) + H _ { K o n d o }. \quad ( 5 ) \quad \text{$\quad \text{the} \quad$} \\ \text{spec}.$$

The Kondo Hamiltonian is expressed in terms of antiferromagnetic coupling J K between the impurity spin glyph[vector] S and the spin operator of reservoir electrons placed at the origin glyph[vector] T as [40, 46, 47]

$$H _ { K o n d o } = J _ { K } \vec { \mathcal { S } } \cdot \vec { T }, \ J _ { K } = \frac { t ^ { 2 } } { U } \frac { N ^ { 2 } } { m ( N - m ) }. \quad ( 6 ) \quad \text{where} \ \text{i}$$

glyph[negationslash]

The N 2 -1 traceless components of impurity spin S i ( =1 2 i , · · · N 2 -1) are given by S r r , ′ = d d † r r ′ -m/Nδ r r , ′ with the constraint r , r = ′ N,N ′ . Likewise, the spin operator of reservoir electrons placed at the origin is expressed as glyph[vector] T = ∑ kk , ′ rr ′ b † k r Λ i rr ′ b k ′ r ′ , Λ i being the N × N generators of SU( N ) group. Note that S i are N ! m N !( -m )! × N ! m N !( -m )! matrices acting on states with m electrons.

The ground state of spin S =1 2 SU( / N ) impurity considered in this work is characterized by the complete screening of impurity spin which results in the formation of Kondo singlet. The low-energy regime of fullyscreened Kondo effect is consistently describe by FL theory [48-50]. Applying the standard point-splitting procedure [16, 40, 47, 49] to the Hamiltonian Eq. (6) imparts the low energy FL Hamiltonian of SU( N ) Kondo impurity,

$$\mathcal { H } _ { 0 } & = \nu \sum _ { r } \int _ { \varepsilon } \varepsilon \left [ a _ { \varepsilon r } ^ { \dagger } a _ { \varepsilon r } + b _ { \varepsilon r } ^ { \dagger } b _ { \varepsilon r } \right ], & ( 7 ) &$$

$$\mathcal { H } _ { \text{int} } = \sum _ { r < r ^ { \prime } } \overset { \cdot } { \int } _ { \varepsilon _ { 1 - 4 } } \left [ \frac { \phi _ { 1 } } { \pi \nu } + \frac { \phi _ { 2 } } { 4 \pi \nu } \sum _ { j = 1 } ^ { 4 } \varepsilon _ { j } \right ] \colon b _ { \varepsilon _ { 1 } r } ^ { \dagger } b _ { \varepsilon _ { 2 } r } b _ { \varepsilon _ { 3 } r ^ { \prime } } ^ { \dagger } b _ { \varepsilon _ { 4 } r ^ { \prime } } \colon. \quad \text{He}.$$

$$\mathcal { H } _ { 0 } & = \nu \sum _ { r } \int _ { \varepsilon } \varepsilon \left [ a _ { \varepsilon r } ^ { \cdot } a _ { \varepsilon r } + b _ { \varepsilon r } ^ { \cdot } b _ { \varepsilon r } \right ], \quad ( i ) \\ \mathcal { H } _ { \text{el} } & = - \sum _ { r } \int _ { \varepsilon _ { 1 - 2 } } \left [ \frac { \alpha _ { 1 } } { 2 \pi } ( \varepsilon _ { 1 } + \varepsilon _ { 2 } ) + \frac { \alpha _ { 2 } } { 4 \pi } ( \varepsilon _ { 1 } + \varepsilon _ { 2 } ) ^ { 2 } \right ] b _ { \varepsilon _ { 1 } r } ^ { \dagger } b _ { \varepsilon _ { 2 } r }, \\ \mathcal { H } _ { \text{int} } & = \sum _ { r < r ^ { \prime } } \int _ { \varepsilon _ { 1 - 4 } } \left [ \frac { \phi _ { 1 } } { \pi \nu } + \frac { \phi _ { 2 } } { 4 \pi \nu } \sum _ { i = 1 } ^ { 4 } \varepsilon _ { j } \right ] \colon b _ { \varepsilon _ { 1 } r } ^ { \dagger } b _ { \varepsilon _ { 2 } r } b _ { \varepsilon _ { 3 } r ^ { \prime } } ^ { \dagger } b _ { \varepsilon _ { 4 } r ^ { \prime } } \colon \quad \text{Here are} \quad \text{pri} i$$

The PH symmetric version of Eq. (7) is originally proposed by Nozieres [48] which is commonly known as Nozieres FL theory. In Eq. (7) the density of states per species for a one dimensional channel is represented by the symbol ν . The scattering (elastic) effects in the FL are accounted for by the Hamiltonian H el , where α 1 and α 2 are the first and second generations of Nozieres FL coefficients respectively. The four fermions term represents the interaction part of the Hamiltonian H int which is expressed in terms of FL parameters φ 1 and φ 2 . These FL parameters are related to the associated Kondo temperature of the corresponding SU( N ) impurity. The FL parameters characterizing the scattering effects are connected to those of interaction effects by the relation, α 1 =( N -1) φ 1 and α 2 =( N -1) φ / 2 4. In addition the Bethe ansatz provides further link between α 1 and α 2 [40, 47],

$$\det _ { \mathcal { A } } \frac { \alpha _ { 2 } } { \alpha _ { 1 } ^ { 2 } } = \frac { N - 2 } { N - 1 } \frac { \Gamma ( 1 / N ) \tan ( \pi / N ) } { \sqrt { \pi } \Gamma \left ( \frac { 1 } { 2 } + \frac { 1 } { N } \right ) } \cot \left [ \frac { m \pi } { N } \right ]. \ \ ( 8 )$$

Where Γ( x ) is the Euler's gamma-function. Therefore the low energy FL Hamiltonian Eq. (7) is completely specified by only one FL parameter, say α 1 . We make a connection of α 1 with the corresponding Kondo temperature such that T SU( N ) K =1 /α 1 , the N -dependence in FL parameters is implicit. Note that we have retained upto the four fermions term in Eq. (7), the higher-order terms produce the current correction beyond cubic order in applied bias and temperature gradient which is beyond the scope of present work.

It is then a straightforward procedure to proceed with the calculation of physical observables by treating the scattering Hamiltonian H el and interaction part H int perturbatively. However, in the spirit of Nozieres phenomenology, the scattering effects are fully described by an energy-dependent phase shift δ el r ( ε ). The Kondo singlet acts as the scatterer for the incoming electrons from the leads. Outgoing and incoming electrons are then differ from each other by the elastic phase shift δ el r ( ε ). The Nozieres FL parameters α 1 and α 2 are the first and second order coefficients in the Taylor-series expansion of elastic phase shift. While the scattering effects are easily accounted for by the elastic phase-shift, the perturbative treatment of H int produces complicated self-energy diagrams. This complication can be simplified a bit by including the Hartree contribution of self-energy in the elastic phase shift [40, 47]. Then the Taylor expansion of phase shift reads glyph[negationslash]

$$\dots \quad & \circ \Pr \infty \, \text{min} \, \text{max} \\ ( 7 ) & \delta _ { r } ( \varepsilon ) = \delta _ { 0 } + \alpha _ { 1 } \varepsilon + \alpha _ { 2 } \varepsilon ^ { 2 } - \sum _ { r ^ { \prime } \neq r } \left [ \phi _ { 1 } \int _ { - \infty } ^ { \infty } d \varepsilon \delta n _ { r ^ { \prime } } ( \varepsilon ) \\ \iota _ { \varepsilon _ { 2 } r }, & \quad & + \frac { \phi _ { 2 } } { 2 } \left ( \varepsilon \int _ { - \infty } ^ { \infty } d \varepsilon \delta n _ { r ^ { \prime } } ( \varepsilon ) + \int _ { - \infty } ^ { \infty } d \varepsilon \varepsilon \delta n _ { r ^ { \prime } } ( \varepsilon ) \right ) \right ]. \ ( 9 ) \\ & \quad - \quad. \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad.$$

Here the zero-energy phase shift of SU( N ) Kondo impurity with m electrons is

$$\delta _ { 0 } = \frac { m \pi } { N }. \quad \quad \quad ( 1 0 )$$

In Eq. (9) we used the definition of the actual FL quasiparticle distribution relative to the Fermi-energy ε F as

δn r ( ε ) ≡ n r ( ε ) -Θ( ε F -ε ) = 〈 b † k r b k r 〉 -Θ( ε F -ε ), Θ being the step function. Using Eq. (16) we expressed the average 〈 b † k r b k r 〉 in terms of the equilibrium Fermidistribution functions f γ ( ε )= [ 1 + exp ( ε -µ γ T γ )] -1 of the left and right reservoirs; 〈 b † k r b k r 〉 = cos 2 θf L + sin 2 θf R . In addition, we have implemented the specific choice of Fermi-level such that

$$\int _ { - \infty } ^ { \infty } d \varepsilon \delta n _ { \mathfrak { r } } ( \varepsilon ) = 0. \quad \quad ( 1 1 ) \quad \begin{matrix} T \mathfrak { h } \\ \text{cl} \mathfrak { h } \\ \text{ha} \end{matrix}$$

This equation is always satisfied as far as the condition µ L cos 2 θ + µ R sin 2 θ = ε F is full-filled. We then made the following specification for the chemical potentials of the reservoirs,

$$\mu _ { L } = \ e \Delta V \sin ^ { 2 } \theta \equiv \ \frac { \ e \Delta V } { 2 } \left ( 1 - \mathcal { C } \right ), \quad \left ( 1 2 \right ) \quad \text{The} \quad \text{plet} \quad \text{$\tau_{\tau_{\epsilon}}$}$$

$$\mu _ { \text{R} } = - e \Delta V \cos ^ { 2 } \theta \equiv - \frac { e \Delta V } { 2 } \left ( 1 + \mathcal { C } \right ), \quad \left ( 1 3 \right ) \quad \text{trac} \\ \text{erGV}$$

to make ε F =0 It is also noted that the details related to . the choice of the temperatures in the reservoirs do not affects the necessary condition to satisfy the Eq. (11). To be more general, we do not yet impose any restriction on the choice of T L and T R . Using these specification of chemical potentials and temperatures of the reservoirs, the straightforward integration of phase shift expression Eq. (9) leads

$$\delta _ { r } ( \varepsilon ) = \delta _ { 0 } + \alpha _ { 1 } \varepsilon + \alpha _ { 2 } \left ( \varepsilon ^ { 2 } - \mathcal { A } \right ). \quad \ \left ( 1 4 \right ) \quad \text{of $f$}$$

To obtained Eq. (14) we have made the use of FL identity α 2 = ( N -1) φ / 2 4 and the new definition,

$$\mathcal { A } { = } \frac { 1 } { 6 } \left [ ( \pi T _ { \mathrm L } ) ^ { 2 } ( 1 { + } \mathcal { C } ) { + } ( \pi T _ { \mathrm R } ) ^ { 2 } ( 1 { - } \mathcal { C } ) { + } \frac { 3 } { 2 } ( 1 { - } \mathcal { C } ^ { 2 } ) ( e \Delta V ) ^ { 2 } \right ]. \quad \text{consist} \quad \text{(s)}$$

In the following section the scattering effects in addition to the Hartree contribution to the self energy correction will be accounted for by the Eq. (14). To obtain the self energy correction beyond Hartree contribution we will be treating the interaction Hamiltonian H int perturbatively with the small parameters ( e ∆ V, T L , T R ) /T SU( N ) K .

## III. CURRENT CALCULATION

Using the basis of scattering states that includes the elastic effects and Hartree term, we cast the charge current expression into the form [40, 51]

$$\hat { I } { = } \frac { e } { 2 h \nu } \sum _ { r } \sin 2 \theta \left [ a _ { r } ^ { \dagger } ( x ) b _ { r } ( x ) { - } a _ { r } ^ { \dagger } ( { - } x ) \mathcal { E } b _ { r } ( { - } x ) { + } \text{H.c.} \right ], \\ / \tau \varepsilon \, \cdot$$

(15)

for b r ( x )= ∑ k b k r e ikx and E b r ( x )= ∑ k E k b k r e ikx . To write Eq. (15) we have also omitted the terms of the form ∑ r p= , ± p a † r (p x a ) r (p x ) since they do not produce finite contribution to the mean current. In addition, we expressed N × N scattering matrix E k in terms of phase shift expression Eq. (14) such that E k =exp[2 iδ r ( ε k )]. To compute the various observables from Eq. (15) we need the following averages directly obtained from the Glazman-Raikh rotation

$$\text{of} \quad \begin{pmatrix} \lambda \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \infty \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inneutral \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \innegative \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \innormal \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inyes \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \innumeric \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \innot \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \ininno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \, \inno \,
\text{(1.6)} \\ \text{The average of $E$} \, (15)$ provides the elastic current, in-$$

The average of Eq. (15) provides the elastic current (including the corresponding Hartree contribution) which have the compact form analogous to the LandauerB¨ uttiker expression

$$I _ { \text{el} } = \frac { e } { h } \sum _ { r } ^ { N } \int _ { - \infty } ^ { \infty } d \varepsilon \, \mathcal { T } _ { r } ( \varepsilon ) \, \left [ f _ { \text{L} } ( \varepsilon ) - f _ { \text{R} } ( \varepsilon ) \right ]. \quad ( 1 7 )$$

The effective transmission coefficient T r ( ε ) is completely specified by the phase shift expression Eq. (14); T r ( ε ) ≡ ( 1 -C 2 ) sin 2 [ δ r ( ε )]. To write T r ( ε ) into more tractable form, we perform its Taylor expansion in energy and retained upto the second order terms,

$$\intertext { l to } \mathcal { T } _ { r } ( \varepsilon ) = \left ( 1 - \mathcal { C } ^ { 2 } \right ) & \left [ \mathcal { T } _ { 0 } - \alpha _ { 2 } \mathcal { A } \sin 2 \delta _ { 0 } + \alpha _ { 1 } \sin 2 \delta _ { 0 } \, \varepsilon \\ \intertext { l 1 } \text{ion } \mathcal { T } _ { r } & \quad + \left ( \alpha _ { 1 } ^ { 2 } \cos 2 \delta _ { 0 } + \alpha _ { 2 } \sin 2 \delta _ { 0 } \right ) \varepsilon ^ { 2 } \right ].$$

Here T 0 =sin 2 δ 0 is the zero energy transmission coefficient. Then it is trivial procedure to compute the elastic current by plugging in Eq. (18) into Eq. (17). The exact computation of Eq. (17) follows from the consideration of following integrals [51, 52],

$$\mathcal { e } _ { n } = \int _ { - \infty } ^ { \infty } \varepsilon ^ { n } \left [ f _ { L } ( \varepsilon ) - f _ { R } ( \varepsilon ) \right ] d \varepsilon, \, n = 0, 1 \, \text{and} \, 2. \, \ ( 1 9 )$$

Conventional way of calculating the integrals in Eq. (19) consists of Sommerfeld expansion of ∆ f ( ε ) ≡ f L ( ε ) -f R ( ε ) in the small parameters ∆ T ≡ T L -T R and ∆ V . However, the Fourier-transform technique allows us to compute Eq. (19) exactly. Fourier transforming the function ∆ f ( ε ) into real time reads,

$$\Delta f ( t ) = \frac { 1 } { 2 \pi } \int _ { - \infty } ^ { \infty } d \varepsilon \, e ^ { - i \varepsilon t } \Delta f ( \varepsilon ). \quad \ \ ( 2 0 )$$

Performing the n-times partial differentiation of Eq. (20) and taking the limit t → 0 we get,

$$\frac { 2 \pi } { ( - i ) ^ { n } } \ \frac { \partial ^ { n } \Delta f ( t ) } { \partial t ^ { n } } \Big | _ { t = 0 } = \int _ { - \infty } ^ { \infty } d \varepsilon \varepsilon ^ { n } \Delta f ( \varepsilon ). \quad ( 2 1 )$$

Fourier transformation of the Fermi distributions of left and right reservoirs allows us to write

$$\Delta f ( t ) = \frac { i } { 2 } \left [ \frac { T _ { \text{L} } e ^ { - i \mu _ { \text{L} } t } } { \sinh ( \pi T _ { \text{L} } t ) } - \frac { T _ { \text{R} } e ^ { - i \mu _ { \text{R} } t } } { \sinh ( \pi T _ { \text{R} } t ) } \right ]. \quad ( 2 2 )$$

Plugging in Eq. (22) into Eq. (21) with the chemical potentials as specified in Eqs. (12) and (13) we obtain K 0 = ∆ e V , K 1 = ( [ πT L ) 2 -( πT R ) 2 - C 3 ( e ∆ ) V 2 ] / 6 and

$$\underset { \text{phase} } { \text{on, we} } \quad \mathcal { K } _ { 2 } = & \frac { e \Delta V } { 3 } \left [ \frac { ( e \Delta V ) ^ { 2 } } { 4 } \left ( 1 + 3 \mathcal { C } ^ { 2 } \right ) + \frac { 1 - \mathcal { C } } { 2 } ( \pi T _ { L } ) ^ { 2 } + \frac { 1 + \mathcal { C } } { 2 } ( \pi T _ { R } ) ^ { 2 } \right ].$$

For completeness we re-express the elastic current in terms of the integrals in Eq. (19) as

$$I _ { \text{el} } = \frac { N e \left ( 1 - \mathcal { C } ^ { 2 } \right ) } { h } & \left [ \left ( \mathcal { T } _ { 0 } - \alpha _ { 2 } \mathcal { A } \sin 2 \delta _ { 0 } \right ) \mathcal { K } _ { 0 } + \alpha _ { 1 } \sin 2 \delta _ { 0 } \mathcal { K } _ { 1 } \quad \begin{array} { c } \text{via} ; \, 2 \\ \text{eq.} \, \text{(2)} \end{array} \right ) \\ & + \left ( \alpha _ { 1 } ^ { 2 } \cos 2 \delta _ { 0 } + \alpha _ { 2 } \sin 2 \delta _ { 0 } \right ) \mathcal { K } _ { 2 } \right ].$$

Now we turn to the discussion of inelastic effects leaving aside the Hartree contributions, which has been already accounted for by the phase shift expressed in Eq. (14). As we anticipated earlier, the perturbative treatment of H int imparts the interaction corrections to the charge current. This approach requires the expressions of noninteraction Green's functions (GFs) described by H 0 . The matrices of the non-interacting GFs in Keldysh space [53] are given by

$$& \mathcal { G } _ { b b / a a } ( k, \varepsilon ) = & \frac { 1 } { \varepsilon - \varepsilon _ { k } } \tau _ { z } + i \pi \begin{pmatrix} \mathbb { F } _ { b / a } & \mathbb { F } _ { b / a } + 1 \\ \mathbb { F } _ { b / a } - 1 & \mathbb { F } _ { b / a } \end{pmatrix} \delta ( \varepsilon - \varepsilon _ { k } ), \quad \left ( \Sigma ^ { - + } \\ & \mathcal { G } _ { b a / a b } ( k, \varepsilon ) = & i \pi \begin{pmatrix} 1 & 1 \\ 1 & 1 \end{pmatrix} \mathbb { F } _ { a b } \ \delta ( \varepsilon - \varepsilon _ { k } ).$$

Here the parameters F b/a ( ε ) and F ab ( ε k ) are expressed in terms of different populations; F ( b ε k )=2 〈 b † k b k 〉-1, F ( a ε k )=2 〈 a a † k k 〉-1 and F ab =2 〈 b † k a k 〉 . The z-component of Pauli-matrix is represented by τ z . However, in the flat-band limit only the off-diagonal parts of G bb ( k, ε ), namely G + -bb ( k, ε ) and G -+ bb ( k, ε ) produce the finite contribution to the charge current. The straightforward mathematical steps provide the following Fouriertransformed real-time GFs

$$\text{anonymous recursing } \text{can use } \text{$\alpha$} \text{$\alpha$} \\ \mathcal { G } _ { b b } ^ { + - } ( t ) = & - \frac { \pi \nu } { 2 } \left [ \frac { T _ { L } ( 1 + \mathcal { C } ) e ^ { - i \mu _ { L } t } } { \sinh ( \pi T _ { L } t ) } + \frac { T _ { R } ( 1 - \mathcal { C } ) e ^ { - i \mu _ { R } t } } { \sinh ( \pi T _ { R } t ) } \right ], \quad \text{the re} \\ \mathcal { G } _ { a b / b a } ( t ) = & - \frac { \pi \nu } { 2 } \sqrt { 1 - \mathcal { C } ^ { 2 } } \left [ \frac { T _ { L } e ^ { - i \mu _ { L } t } } { \sinh ( \pi T _ { L } t ) } - \frac { T _ { R } e ^ { - i \mu _ { R } t } } { \sinh ( \pi T _ { R } t ) } \right ].$$

Here G + -bb ( ) and t G -+ bb ( ) are connected by causality ret lations. In practice, the GFs expressed in Eqs. (24) and (25) are sufficient for the calculation of charge current. To calculate the inelastic correction to the charge current we then apply the perturbation theory using Keldysh formalism [53],

$$\delta I _ { \text{in} } = \langle T _ { C } \hat { I } ( t ) e ^ { - i \int d t ^ { \prime } \mathcal { H } _ { \text{int} } ( t ^ { \prime } ) } \rangle, \quad \quad ( 2 6 )$$

where C denotes the double-side η = ± Keldysh contour and T C is corresponding time-ordering operator. We used the expression of charge current operator Eq. (15) and interaction Hamiltonian H int into Eq. (26) to obtain the interaction correction to the charge current

$$\delta I _ { \text{in} } = \mathcal { L } \int _ { - \infty } ^ { \infty } \frac { d \varepsilon } { 2 \pi } \left ( \Sigma ^ { - + } - \Sigma ^ { + - } \right ) ( \varepsilon ) i \pi \nu \Delta f ( \varepsilon ). \ \ ( 2 7 ) \quad \stackrel { \chi \cos } { \mathcal { K } _ { 0, 2 } } } { \text{charg} }$$

To arrived from Eq. (26) to Eq. (27) we have already subtracted the diverging terms, which amounts the renormalization of FL coefficients (see Ref. [40] for detail). In addition, we introduced the new notation via; Z = N N ( -1) h eπ ( 1 -C 2 ) cos 2 δ 0 . The self-energies in Eq. (27) are expressed in real-time as

$$\Sigma ^ { \eta _ { 1 } \eta _ { 2 } } ( t ) & = \left ( \frac { \phi _ { 1 } } { \pi \nu ^ { 2 } } \right ) ^ { 2 } \sum _ { k _ { 1 - 3 } } \left [ \mathcal { G } ^ { \eta _ { 1 } \eta _ { 2 } } _ { b b } ( k _ { 1 }, t ) \\ & \quad \mathcal { G } ^ { \eta _ { 2 } \eta _ { 1 } } _ { b b } ( k _ { 2 }, - t ) \mathcal { G } ^ { \eta _ { 1 } \eta _ { 2 } } _ { b b } ( k _ { 3 }, t ) \right ].$$

To compute the self-energies, now we specify the temperatures of the left and right reservoirs T R = T and T L = +∆ T T with ∆ T&gt; 0. In practice one can numerically solve for the self-energy using the GFs of Eq. (25). However, it is manageable to find the analytical expression of the self energy difference to the first order in ∆ T and second order in e ∆ V which reads

$$- \varepsilon _ { k } \right ), \quad \left ( \Sigma ^ { - + } - \Sigma ^ { + - } \right ) ( \varepsilon ) = \frac { \phi _ { 1 } ^ { 2 } } { i \pi \nu } \left [ \frac { 3 } { 4 } ( e \Delta V ) ^ { 2 } ( 1 - \mathcal { C } ^ { 2 } ) + \varepsilon ^ { 2 } + ( \pi T ) ^ { 2 } \\ ( 2 4 ) & & + \frac { \Delta T } { T } ( \pi T ) ^ { 2 } \left ( 1 + \mathcal { C } \right ) \right ]. \quad \left ( 2 9 \right )$$

To arrive from Eq. (28) to Eq. (29) we came across the integral of the form,

$$\mathcal { Z } ( a, T ) = \int _ { - \infty } ^ { \infty } \frac { e ^ { i a t } } { \sinh ^ { 3 } ( \pi T t ) } \ d t. \quad \quad ( 3 0 )$$

The singularity of the integral in Eq. (30) is removed by shifting the time contour by iη , η → 0 in the complex plane. The parameter η is chosen such that ηD glyph[greatermuch] 1 and ( ηT, η ∆ T, η ∆ ) V glyph[lessmuch] 1 with D the band cutoff. We chose the rectangular contour enclosing the singularity at t =0 and use the Cauchy's residue theorem to arrive at the result,

$$\mathcal { Z } ( a, T ) = - i \pi \, \frac { a ^ { 2 } + ( \pi T ) ^ { 2 } } { ( \pi T ) ^ { 2 } } \frac { 1 } { \exp ( a / T ) + 1 }. \quad ( 3 1 )$$

Equation (29) contains all possible terms up to the linear response in ∆ T and ∆ V . Therefore plugging in Eq. (29) into Eq. (27) provides interaction correction up to the quadratic order in ∆ T and ∆ V . To make interaction contributions to the charge current more symmetrical with that of elastic effects, we write

$$\partial _ { \Omega } & = & \frac { N e ( 1 - \mathcal { C } ^ { 2 } ) } { h } \frac { 1 } { 2 } \frac { 1 } { N - 1 } \cos 2 \delta _ { 0 } \, \alpha _ { 1 } ^ { 2 } \left [ \mathcal { K } _ { 2 } \\ \partial _ { \Omega } & \quad + \left ( \frac { \Delta T } { T } ( \pi T ) ^ { 2 } \, ( 1 + \mathcal { C } ) + ( \pi T ) ^ { 2 } \right ) \mathcal { K } _ { 0 } \right ]. \quad ( 3 2 ) \\ \partial _ { \Omega }$$

This equation correctly reproduces the interaction correction up to the quadratic response with the coefficients K 0 2 , given in Eq. (19). Using Eq. (23) and (32), the charge current is given by

$$I _ { c } = I _ { \text{el} } + \delta I _ { \text{int} }.$$

## IV. RESULTS AND DISCUSSION

The non-linear Seebeck effect is quantified by the Seebeck coefficient defined as the ratio of thermo-voltage developed under the condition of zero charge current, ∆ V th ≡ ∆ V | I c =0 , to the applied temperature gradient [54, 55],

$$\mathcal { S } \equiv - \frac { \Delta V _ { \text{th} } } { T _ { L } - T _ { R } } \Big | _ { I _ { \text{c} } = 0 } \,. \quad \quad \quad ( 3 4 ) \quad \text{$\frac{ \text{im}_{\cdot}$} { \text{as}$}$$

In fact, the Seebeck coefficient Eq. (34) contains additional information than the electrical and thermal conductance measurements [56]. While the electrical conductance depends merely on the density of states at the Fermi level, Seebeck coefficient reveals its slope [57]. In addition, the Seebeck coefficient provides the useful informations related to the average energy of charge carriers contributing to the transport processes [58]. We characterize the Seebeck coefficient of a SU( N ) Kondo impurity by defining the dimensionless form of charge current accounting upto the quadratic responses,

$$\mathcal { J } ( N, m ) \equiv \frac { I _ { \text{c} } ( N, m ) } { G _ { 0 } ( N ) T _ { K } ^ { \text{SU} ( N ) } } \quad \quad \quad ( 3 5 ) \quad \text{The $\i$} \quad \text{ takes}$$

$$= \mathcal { L } _ { 1 } ^ { 1 } \Delta \overline { V } + \mathcal { L } _ { 2 } ^ { 1 } \Delta \overline { T } + \mathcal { L } _ { 1 } ^ { 2 } \Delta \overline { V } ^ { 2 } + \mathcal { L } _ { 2 } ^ { 2 } \Delta \overline { T } ^ { 2 } + \mathcal { L } _ { 1 2 } ^ { 1 1 } \Delta \overline { V } \Delta \overline { T }.$$

$$& \mathbb { m } ( N ) \mathcal { I } _ { K } ^ { \widetilde { \mathcal { M } } \cdots \cdots \cdots } \\ & \mathbb { m } ^ { 1 } _ { 2 } \Delta \overline { T } + \mathcal { L } _ { 1 } ^ { 2 } \Delta \overline { V } ^ { 2 } + \mathcal { L } _ { 2 } ^ { 2 } \Delta \overline { T } ^ { 2 } + \mathcal { L } _ { 1 2 } ^ { 1 1 } \Delta \overline { V } \Delta \overline { T }.$$

The maximum conductance of SU( N ) Kondo impurity in the presence of asymmetry is expressed by the relation G 0 ( N ) = ( 1 -C 2 ) Ne /h 2 . From now we use specify the electronic charge e = -1 and the convention ∆ V &gt; 0 and ∆ T &gt; 0. The quantities written in over-line letters represent that they are normalized with corresponding Kondo temperature: T ≡ T/T SU(N) K , ∆ T ≡ ∆ T/T SU(N) K and ∆ V ≡ ∆ V/T SU(N) K . From Eq. (23) and (32) we obtained the transport coefficients L i j and L 11 12 , i, j = 1 2 , for the SU( N ) Kondo impurity,

$$\mathcal { L }$$

$$& \text{for the SU} ( N ) \, Kondo \, \text{impurity}, \quad \\ & \mathcal { L } _ { 1 } ^ { 1 } = \left [ \sin ^ { 2 } \left ( \frac { \pi m } { N } \right ) + \frac { 1 } { 3 } \frac { N + 1 } { N - 1 } \cos \left ( \frac { 2 \pi m } { N } \right ) ( \pi \overline { T } ) ^ { 2 } \right ], \quad & \text{Seebé parar} \\ & \mathcal { L } _ { 2 } ^ { 1 } = - \, \frac { \pi ^ { 2 } } { 3 } \overline { T } \sin \left ( \frac { 2 \pi m } { N } \right ), \ \mathcal { L } _ { 1 } ^ { 2 } = \frac { 1 } { 2 } \mathcal { C } \sin \left ( \frac { 2 \pi m } { N } \right ), \quad & \text{etur} \ \text{ditor to ex;} } \\ & \mathcal { L } _ { 2 } ^ { 2 } = - \, \frac { \pi ^ { 2 } } { 6 } \sin \left ( \frac { 2 \pi m } { N } \right ), \quad & \text{linear} \\ & \mathcal { L } _ { 1 2 } ^ { 1 1 } = - \, \frac { \pi ^ { 2 } } { 3 } \overline { T } \left [ \mathcal { B } \cos \left ( \frac { 2 \pi m } { N } \right ) + 2 \mathcal { C } \mathcal { A } \sin \left ( \frac { 2 \pi m } { N } \right ) \right ]. \\ & \text{The coefficient $A$ is defined in Eq. (8) and $\mathcal{B}$ stands for $\mathcal{B}$ to $text$}$$

The coefficient A is defined in Eq. (8) and B stands for

$$\mathcal { B } \equiv \frac { \mathcal { C } ( N - 2 ) - N - 1 } { N - 1 }. \quad \quad ( 3 7 ) \quad \text{hal} \\ \text{dis} \cdot$$

From Eq. (36) it is seen that the transport coefficients accounting for the linear and quadratic correction in temperature gradient are connected by the relation L 1 2 = 2 T L 2 2 . It is apparent that, merely the asymmetry of the junction is responsible to have the quadratic correction in voltage bias. For half-filled SU( N ) Kondo effects, we observed that L 1 2 = L 2 1 = L 2 2 = 0, therefore, corresponding thermoelectric properties are governed by only two coefficients L 1 1 and L 11 12 . This fact explains that the half-filled SU( N ) Kondo impurity do not offers finite thermo-power even in quadratic-response level of calculations. Another important conclusion can be drawn form Eq. (36) is as follows; for the perfectly symmetrical quarter-filled SU( N ) Kondo correlated systems, the combine effects of temperature gradient and voltage bias tend to vanish L 11 12 ∣ ∣ C =0 ( N,N/ 4)=0. Furthermore, the coefficients characterizing the voltage response do not acquire the temperature correction. These facts should make the non-linear thermoelectric measurement of beyond halffilled SU(4) systems as a trivial procedure. To have more insights of the thermoelectric production in SU( N ) Kondo systems, we solve the zero current condition of the Eq. (35) to get the thermo-voltage upto the quadratic terms in ∆ T ,

$$- \Delta \overline { V } _ { t h } = \mathcal { S } ^ { L R } \Delta \overline { T } + \delta \mathcal { S } ( \Delta \overline { T } ) ^ { 2 } + \mathcal { O } ( \Delta \overline { T } ) ^ { 3 }. \quad ( 3 8 )$$

The Seebeck coefficient S as defined in Eq. (34) then takes the form,

$$\mathcal { S } = \mathcal { S } ^ { L R } + \delta \mathcal { S } \Delta \overline { T } + \mathcal { O } ( \Delta \overline { T } ) ^ { 2 }. \quad \ \ ( 3 9 )$$

Here S LR is the linear response Seebeck coefficient and its first order ∆ T correction is defined by δ S ,

$$\mathcal { S } ^ { L R } \equiv \frac { \mathcal { L } _ { 2 } ^ { 1 } } { \mathcal { L } _ { 1 } ^ { 1 } }, \quad \ \ \ \ \ \ ( 4 0 )$$

$$\delta \mathcal { S } \equiv \left [ \stackrel { \ell } { \frac { \mathcal { L } _ { 2 } ^ { 2 } } { \mathcal { L } _ { 1 } ^ { 1 } } } - \frac { \mathcal { L } _ { 2 } ^ { 1 } \mathcal { L } _ { 1 2 } ^ { 1 1 } } { \left ( \mathcal { L } _ { 1 } ^ { 1 } \right ) ^ { 2 } } + \frac { \left ( \mathcal { L } _ { 2 } ^ { 1 } \right ) ^ { 2 } \mathcal { L } _ { 1 } ^ { 2 } } { \left ( \mathcal { L } _ { 1 } ^ { 1 } \right ) ^ { 3 } } \right ]. \quad ( 4 1 )$$

The transport coefficients defining the linear response Seebeck coefficient S LR are independent of asymmetry parameter C . However, the first order correction δ S bears the strong dependences on the asymmetry parameter via the transport coefficients L 2 1 and L 11 12 . In addition, for the symmetrical setups, we use the Eq. (36) to express the correction factor δ S entirely in terms of linear-response coefficients,

$$\delta \mathcal { S } | _ { \mathcal { C } = 0 } = \frac { \mathcal { S } ^ { L R } } { \overline { T } } \left [ \frac { \sin ^ { 2 } \left ( \frac { \pi m } { N } \right ) } { \mathcal { L } _ { 1 } ^ { 1 } } - \frac { 1 } { 2 } \right ]. \quad \ \ ( 4 2 )$$

To study the effects of coupling asymmetry on the thermoelectric transport properties, we categorize the SU( N ) Kondo impurity into two broad classes, namely, half-filled (PH symmetric) and beyond half-filled and discuss them separately.

## 1. PH symmetric SU( N ) Kondo effects

As we anticipated earlier that for the half-filled SU( N ) Kondo effects the transport coefficients satisfy the relation L 1 2 = L 2 1 = L 2 2 =0, therefore, corresponding

FIG. 2. Linear (LR) and non linear (BLR) Seebeck coefficients with PH symmetric SU(2) and SU(4) Kondo effects for fixed value of the potential scattering δ P .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00415v1-figure-9f892cfe-3aa1-422d-96d6-0c93dcb531ea.jpg)

thermoelectric properties are derived solely from the coefficients L 1 1 and L 11 12 . The non-zero transport coefficients of PH symmetric SU( N ) Kondo effects are summarized below,

$$\mathcal { L } _ { 1 } ^ { 1 } ( N, N / 2 ) = \left [ 1 - \frac { 1 } { 3 } \frac { N + 1 } { N - 1 } ( \pi \overline { T } ) ^ { 2 } \right ], \quad \quad ( 4 3 ) \quad \overset { \text{ital s} } { \Delta \overline { T } } =$$

$$\mathcal { L } _ { 1 2 } ^ { 1 1 } ( N, N / 2 ) = \bar { \frac { 1 } { 3 \bar { T } } } \left [ \frac { \mathcal { C } ( N - 2 ) - \bar { N } - 1 } { N - 1 } \right ] ( \pi \overline { T } ) ^ { 2 }. \ ( 4 4 ) \quad \text{control} \\ \text{erenc}.$$

While for the conventional SU(2) Kondo effects the parameter C does not affects the cross-coefficient, the corresponding measurement in SU( N &gt; 2) PH symmetric systems depends on the coupling asymmetry. However, the PH symmetry of the Kondo impurity realized in QDs is exact only if the dot is tuned to the middle of Coulomb valley [59]. This indicates the possibility of breaking the underlying PH symmetry. This weakly broken PH symmetry of Kondo correlated systems is accounted for by re-normalizing the reference phase shifts such that [60-62],

$$\delta _ { 0 } \rightarrow \tilde { \delta _ { 0 } } = \delta _ { 0 } + \delta _ { P }, \ \delta _ { 0 } \gg \delta _ { P }. \quad \ \ ( 4 5 ) \quad \text{Th} \,.$$

This potential scattering provides the repulsive interactions which breaks the Kondo singlet and contributes to inelastic processes [63]. The first order transport coefficients in Eq. (36) for PH symmetric Kondo correlated systems with an account of the potential scattering effects are then given by

$$& \dots \dots \dots \dots \to \cdots \to \cdots \to \cdots \\ & \mathcal { L } _ { 1 } ^ { 1 } ( N, N / 2 ) \Big | _ { P } = \cos ^ { 2 } \delta _ { P } \left [ 1 - \frac { ( \pi \overline { T } ) ^ { 2 } } { 3 } \frac { N + 1 } { N - 1 } \frac { 2 \cos 2 \delta _ { P } } { 1 + \cos 2 \delta _ { P } } \right ], \quad \mathcal { L } _ { 1 } ^ { 2 } \\ & \mathcal { L } _ { 2 } ^ { 1 } ( N, N / 2 ) \Big | _ { P } = \cos ^ { 2 } \delta _ { P } \left [ \frac { ( \pi \overline { T } ) ^ { 2 } } { 3 \overline { T } } \frac { 2 \sin 2 \delta _ { P } } { 1 + \cos 2 \delta _ { P } } \right ]. \quad ( 4 6 ) \quad \mathcal { L } _ { 1 2 } ^ { 1 1 } \\ & \mathcal { r } = \omega \cdots \cdots \cdots \cdots \cdots \cdots$$

Eq. (46) allows us to compute the linear response Seebeck coefficient of PH symmetric SU( N ) Kondo effects with small potential scattering,

$$\mathcal { S } ^ { \text{LR} } ( N, N / 2 ) \Big | _ { P } = \frac { 2 } { 3 } \, \frac { 1 } { \overline { T } } \, \frac { ( \pi \overline { T } ) ^ { 2 } } { 1 - \frac { ( \pi \overline { T } ) ^ { 2 } } { 3 } \frac { N + 1 } { N - 1 } } } \delta _ { P } + \mathcal { O } ( \delta _ { P } ) ^ { 3 }. \quad \text{ ificen} \\ \text{No to that } \text{do to the normal or for } \text{$N\perp1\}/\int\,/\int\,$N\end{array}$$

Note that due to the numerical factor ( N +1) ( / N -1) in the denominator of Eq. (47), among PH symmetric generalizations of SU( N ) the SU(2) Kondo correlated systems offer highest value of the linear response Seebeck coefficient in the presence of finite potential scattering. Plugging in the Eq. (45) into the transport coefficients Eq. (36) and using them into Eq. (41), we get the first order correction to the Seebeck coefficient upto the linear order in δ P ,

$$\frac { 1 } { 0. 1 4 } \Big | _ { P } & = \frac { \pi ^ { 2 } } { 3 } \, \frac { 1 - \frac { ( \pi \overline { T } ) ^ { 2 } } { 3 } \left ( \frac { N + 1 } { N - 1 } + 2 \mathcal { B } \right ) } { \left [ 1 - \frac { ( \pi \overline { T } ) ^ { 2 } } { 3 } \frac { N + 1 } { N - 1 } \right ] ^ { 2 } } \delta _ { P } + \mathcal { O } ( \delta _ { P } ) ^ { 3 }. \\ \text{oeff} \quad & \text{For SU} (2) Kondo effects the correction $Ea$}. (48) is indeq.$$

For SU(2) Kondo effects the correction Eq. (48) is independent of the asymmetry parameter as can be inferred from Eq. (37). However for SU(4) and other PH symmetric version of SU( N ), the first order correction to the Seebeck effect is weakly asymmetry dependent via the coefficient B ( C ). The linear and non linear Seebeck coefficient with PH symmetric SU(2) and SU(4) Kondo effects are shown in Fig. 2 with the choice of potential scattering term δ P = 0 1 and temperature gradient . ∆ T = 0 05. These significant enhancement of BLR See-. beck coefficients with respect to the corresponding LR contribution get further improve at relatively high reference temperature and large temperature drop across the junction.

## 2. Beyond half-filled SU( N ) Kondo effects

The Kondo correlated systems with N&gt; 2 provide the realization of paradigmatic PH-asymmetric setups. First we start form the SU(3) Kondo effects. The SU(3) Kondo effect can occur either with single electron or two electrons. Furthermore, the SU(3) Kondo systems do not offer the PH symmetric analog. The physics of SU(3) Kondo effect with one and two electrons is related with each other by PH symmetry transformation. Therefore, we discuss the single electron SU(3) Kondo systems, which will ultimately provide the corresponding informations of two electron case. For single electron SU(3) Kondo effects the transport coefficients in Eq. (36) are simplified as

$$\mathcal { L }$$

$$\text{ted} _ { \text{eff} } ^ { \text{ } \colon \infty \nmid \infty \, \infty } \, \text{$\infty$} \\ \mathcal { L } _ { 1 } ^ { 1 } ( 3, 1 ) & = \frac { 3 } { 4 } \left [ 1 - \frac { 4 } { 9 } ( \pi \overline { T } ) ^ { 2 } \right ], \, \mathcal { L } _ { 2 } ^ { 1 } ( 3, 1 ) = - \frac { \pi ^ { 2 } \overline { T } } { 2 \sqrt { 3 } }, \\ \cdot \right ], \, \quad \mathcal { L } _ { 1 } ^ { 2 } ( 3, 1 ) & = \mathcal { C } \frac { \sqrt { 3 } } { 4 }, \, \mathcal { L } _ { 2 } ^ { 2 } ( 3, 1 ) = - \frac { \pi ^ { 2 } } { 4 \sqrt { 3 } }, \\ \mathcal { L } _ { 1 2 } ^ { 1 1 } ( 3, 1 ) & = - \frac { \pi ^ { 2 } } { 3 } \overline { T } \left [ \frac { \mathcal { C } } { 4 } \left ( 1 - 2 \sqrt { \frac { 3 } { \pi } } \frac { \Gamma [ 1 / 3 ] } { \Gamma [ 5 / 6 ] } \right ) - 1 \right ]. \, ( 4 9 ) \\ \text{ed} _ { \text{ } } \, \quad \, \quad \, \quad \, \, \,$$

Therefore while the cross coefficient L 11 12 (3 , 1) is weakly asymmetry dependent, the coefficient L 2 1 (3 1) , is strongly influenced by C . Since all the transport coefficients in Eq. (49) are non-zero, one can solve the zerocurrent equation to get the thermo-voltage developed in SU(3) Kondo effects.

Now we turn to the discussion of SU(4) Kondo effects out of PH symmetric situation. The SU(4) Kondo effects can accommodate upto three electrons. While the two electron case is suffers from the PH symmetry, the

FIG. 3. Left panel: Plot of asymmetry dependent zero current lines in PH asymmetric SU(4) Kondo effects within the qudratic response level of calculations as a function of applied voltage bias and temperature gradient at reference temperature T =0 2. . Right panel: Corresponding Seebeck coefficients for given asymmetry parameter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00415v1-figure-c89696dc-cf40-4b86-aa48-7881e75abf06.jpg)

single and three electron SU(4) systems are regarded to have good thermoelectric performance. Furthermore, the single electron and three electron systems are related to each other by the PH symmetry transformation. Therefore we discuss in details about the thermoelectric of single electron SU(4) Kondo effects. The corresponding transport coefficients are obtained as

$$\mathcal { L } _ { 1 } ^ { 1 } ( 4, 1 ) & = \frac { 1 } { 2 }, \ \mathcal { L } _ { 2 } ^ { 1 } ( 4, 1 ) = - \frac { \pi ^ { 2 } } { 3 } \overline { T }, \ \mathcal { L } _ { 1 } ^ { 2 } ( 4, 1 ) = \frac { \mathcal { C } } { 2 }, \\ \mathcal { L } _ { 2 } ^ { 2 } ( 4, 1 ) & = - \frac { \pi ^ { 2 } } { 6 }, \ \mathcal { L } _ { 1 2 } ^ { 1 1 } ( 4, 1 ) = - \frac { 4 \pi ^ { 2 } \mathcal { C } \overline { T } } { 9 \sqrt { \pi } } \frac { \Gamma [ 1 / 4 ] } { \Gamma [ 3 / 4 ] }. \ \ ( 5 0 ) \quad \text{These} \quad \text{single} \quad \text{age}$$

The cross coefficient L 11 12 (4 , 1) glyph[similarequal]7 32 . C T is very large as compared to other coefficients for relatively large asymmetry parameter. In addition the other coefficient L 2 1 (4 1) is also strongly asymmetry dependent. , Presence of these coefficients is solely manifested by the finite asymmetry of the junction. Therefore we argue that measuring this cross coefficient would be useful while identifying the asymmetry of the junction in addition to its physical implications. Just form the structure of Eq. (50), it is seen that the thermoelectric transport properties of beyond half-filled SU(4) Kondo effects can be easily manipulated by tuning the junction asymmetry. It appears that the effect of asymmetry becomes more pronounce in relatively high temperature gradient regime. The asymmetry parameter C mainly causes to shift the zero-current line either upward or downward with respect to the perfectly symmetric setup. As shown in Fig. 3 the positive value of the asymmetry parameter increases the thermo-voltage, while the opposite effects are apparent for the corresponding negative values. In addition, the beyond linear response contribution always overshoots the corresponding linear response value irrespective of the coupling asymmetry.

## 3. Paradigmatic SU(4) Kondo effects

The cosine factor cos 2 δ 0 in front of the expression of the inelastic current dramatically modifies the low energy transport behavior of SU( N ) Kondo effects. In case of the SU( N ) systems with m electrons satisfying the specific combination such that m/N =(2n + 1) / 4 for n = 0 and 1, the cosine factor cos 2 δ 0 in Eq. (32) amounts to nullify the whole expression. For these specific systems, the beyond Hartree contribution to the self-energy becomes zero being corresponding Hartree contribution finite. In addition, for PH symmetric SU(N) Kondo effects the Hartree contribution vanishes and beyond Hartree contribution becomes finite. Interestingly, the PH asymmetric SU(4) Kondo correlated systems offer vanishing non-Hartree contribution to the self-energy. Since the Hartree contributions can be straightforwardly accounted for by including it in phase shift, the beyond-half filled SU(4) systems can be exactly solved within cubic response and even beyond. This paradigmatic simplication is also applicable for some SU(12) generalizations. From Eq. (23) we obtained two non-zero cubic response coefficients L 3 1 (4 , 1) and L 12 12 (4 , 1) contributing to the charge current of beyond-half filled SU(4) Kondo impurity as

$$\intertext { i o n. } _ { \text{:tric} } \mathcal { J } ( 4, 1 ) | _ { \text{cubic} } = \mathcal { L } _ { 1 } ^ { 3 } ( 4, 1 ) ( \Delta \overline { V } ) ^ { 3 } + \mathcal { L } _ { 1 2 } ^ { 1 2 } ( 4, 1 ) \Delta \overline { V } ( \Delta \overline { T } ) ^ { 2 }.$$

Here the transport coefficients are

$$\mathcal { L } _ { 1 } ^ { 3 } ( 4, 1 ) = - \frac { 1 - 3 \mathcal { C } ^ { 2 } } { 9 \sqrt { \pi } } \frac { \Gamma [ 1 / 4 ] } { \Gamma [ 3 / 4 ] }, \, \mathcal { L } _ { 1 2 } ^ { 1 2 } ( 4, 1 ) = - \mathcal { C } \frac { 2 \pi ^ { 2 } } { 9 \sqrt { \pi } } \frac { \Gamma [ 1 / 4 ] } { \Gamma [ 3 / 4 ] }.$$

These equations show that for the perfectly symmetrical single electron SU(4) Kondo setups, the effects of voltage bias and temperature gradient are not correlated even in cubic response level of calculations. Therefore, only the asymmetry can derives these systems to have combine interplay of voltage bias and temperature gradient. The effects of asymmetry parameter on the Seebeck coefficient in cubic response level of calculations has been presented in Fig. 4 with an example of single electron SU(4) Kondo effects. From Fig. 4 it is seen that with the proper choice (positive value) of asymmetry parameter C the non-linear Seebeck coefficient gets significantly enhanced over the corresponding perfectly symmetrical coupling. This effect is associated with strong asymmetry of the beyond linear response transmission coefficient Eq. (18).

Finally we want to mention that the non-linearly has been also studied by generalizing the definition of Seebeck coefficient with constant current condition [64-66] such that

$$\overline { \mathcal { S } } ( N, m ) = \frac { \partial \mathcal { J } ( N, m ) } { \partial \Delta \overline { T } } \Big / \frac { \partial \mathcal { J } ( N, m ) } { \partial \Delta \overline { V } }. \quad ( 5 1 )$$

In the linear response level of calculation the response coefficient defined in Eq. (51) coincides with the Seebeck coefficient given by Eq. (34). Though their behaviors in non-linear regime is quite different, it has been argued that the coefficient S is indeed experimentally accessible [64] and can provide an important ingredient for the propose of temperature sensing. These effects have been already studied in conventional SU(2) Kondo regime accounting for the linear response of temperature gradient

FIG. 4. Left panel: Lines of zero charge currents in a single electron SU(4) Kondo impurity within cubic response level of calculations. The temperatures of left and right reservoirs (normalized with corresponding Kondo temperature) are varied for given asymmetry parameter at fixed voltage drop ∆ =0 1. V . Right panel: The Seebeck coefficients as a function of asymmetry parameter with beyond half-filled SU(4) Kondo effects for fixed voltage drop ∆ V =0 1. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00415v1-figure-a8692a87-a8e3-493e-bced-b84a139e9b69.jpg)

and finite voltage bias [64, 65]. The central result of the paper expressed in Eq. (36) paved a straightforward way of extending their study with an account of strong nonlinearity in more exotic Kondo correlated system.

## V. CONCLUSIONS

We developed a theoretical framework based on a local Fermi-liquid theory in combination with the out of equilibrium Keldysh approach to study the influences of coupling asymmetry on the thermoelectric transport

- [1] M. S. Dresselhaus, G. Dresselhaus, X. Sun, Z. Zhang, S. B. Cronin, and T. Koga, Physics of the Solid State 41 , 679 (1999).
- [2] S. Datta, Electronic Transport in Mesoscopic Systems , Cambridge Studies in Semiconductor Physics and Microelectronic Engineering (Cambridge University Press, 1995).
- [3] M. Dresselhaus, G. Chen, M. Tang, R. Yang, H. Lee, D. Wang, Z. Ren, J.-P. Fleurial, and P. Gogna, Advanced Materials 19 , 1043 (2007).
- [4] G. Benenti, G. Casati, K. Saito, and R. Whitney, Physics Reports 694 , 1 (2017), fundamental aspects of steady-state conversion of heat to work at the nanoscale.
- [5] Y. M. Blanter and Y. V. Nazarov, Quantum Transport: Introduction to Nanoscience (Cambridge University Press, Cambridge, England, 2009).
- [6] X. Zhang and L.-D. Zhao, Journal of Materiomics 1 , 92 (2015).
- [7] T. A. Costi and V. Zlatic, Phys. Rev. B 81 , 235127 (2010).
- [8] J. Kondo, Progress of Theoretical Physics 32 , 37 (1964).
- [9] R. Scheibner, H. Buhmann, D. Reuter, M. N. Kiselev, and L. W. Molenkamp, Phys. Rev. Lett. 95 , 176602 (2005).
- [10] S. Jezouin, F. D. Parmentier, A. Anthore, U. Gennser, A. Cavanna, Y. Jin, and F. Pierre, Science 342 , 601

of a strongly coupled SU( N ) Kondo impurity. While the linear response Seebeck coefficient is independent of coupling asymmetry, the fundamental role of nonlinearity towards the enhancement of the Seebeck coefficient with a SU( N ) Kondo setup is explored. In addition, we reported the great enhancement of Seebeck coefficient of Kondo impurities by proper tailoring the coupling asymmetry. We explore the importance of potential scattering on the thermoelectric characterization of PH symmetric SU( N ) Kondo effects. The presented analytical expressions of asymmetry dependent transport coefficients for general SU( N ) Kondo effects allow us to make a close connection of our findings with the experimentally studied SU(2) and SU(4) Kondo effects in complex QDs nano structures. Application of developed theoretical framework for the investigation of thermoelectric properties of more exotic Kondo problems such as multi-stage and multi-terminal Kondo screening appears to be a valid avenues for future research.

## ACKNOWLEDGEMENT

We thank Eva Andrei, Yigal Meir, Yuval Oreg, Giulio Casati and Giuliano Benenti for the illuminating discussions. This work was performed in part at Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1607611. This work was partially supported (MK) by a grant from the Simons Foundation.

(2013).

- [11] Z. Iftikhar, S. Jezouin, A. Anthore, U. Gennser, F. D. Parmentier, A. Cavanna, and F. Pierre, Nature 526 , 233 (2015).
- [12] S. Jezouin, Z. Iftikhar, A. Anthore, F. D. Parmentier, U. Gennser, A. Cavanna, A. Ouerghi, I. P. Levkivskyi, E. Idrisov, E. V. Sukhorukov, L. I. Glazman, and F. Pierre, Nature 536 , 60 (2016).
- [13] M. Ferrier, T. Arakawa, T. Hata, R. Fujiwara, R. Delagrange, R. Weil, R. Deblock, R. Sakano, A. Oguri, and K. Kobayashi, Nature Physics 12 , 230 (2016).
- [14] A. Svilans, M. Josefsson, A. M. Burke, S. Fahlvik, C. Thelander, H. Linke, and M. Leijnse, Phys. Rev. Lett. 121 , 206801 (2018).
- [15] B. Dutta, D. Majidi, A. Garcia Corral, P. A. Erdman, S. Florens, T. A. Costi, H. Courtois, and C. B. Winkelmann, Nano Letters 19 , 506 (2019).
- [16] D. B. Karki and M. N. Kiselev, Phys. Rev. B 96 , 121403(R) (2017).
- [17] J. Azema, A.-M. Dar´ e, S. Sch¨fer, a and P. Lombardo, Phys. Rev. B 86 , 075303 (2012).
- [18] D. B. Karki and M. N. Kiselev, arXiv e-prints , arXiv:1906.00724 (2019).
- [19] P. Jarillo-Herrero, J. Kong, H. S. van der Zant, C. Dekker, L. P. Kouwenhoven, and S. D. Franceschi, Nature 434 , 484 (2005).

- [20] A. Makarovski, J. Liu, and G. Finkelstein, Phys. Rev. Lett. 99 , 066801 (2007).
- [21] A. Makarovski, A. Zhukov, J. Liu, and G. Finkelstein, Phys. Rev. B 75 , 241407 (2007).
- [22] M. Ferrier, T. Arakawa, T. Hata, R. Fujiwara, R. Delagrange, R. Deblock, Y. Teratani, R. Sakano, A. Oguri, and K. Kobayashi, Phys. Rev. Lett. 118 , 196803 (2017).
- [23] T. Hata, R. Delagrange, T. Arakawa, S. Lee, R. Deblock, H. Bouchiat, K. Kobayashi, and M. Ferrier, Phys. Rev. Lett. 121 , 247703 (2018).
- [24] A. J. Keller, S. Amasha, I. Weymann, C. P. Moca, I. G. Rau, J. A. Katine, H. Shtrikman, G. Zar´nd, a and D. Goldhaber-Gordon, Nature Physics 10 , 145 (2014).
- [25] G. C. Tettamanzi, J. Verduijn, G. P. Lansbergen, M. Blaauboer, M. J. Calder´n, o R. Aguado, and S. Rogge, Phys. Rev. Lett. 108 , 046803 (2012).
- [26] K. Le Hur, P. Simon, and D. Loss, Phys. Rev. B 75 , 035332 (2007).
- [27] M.-S. Choi, R. L´pez, o and R. Aguado, Phys. Rev. Lett. 95 , 067204 (2005).
- [28] M. Eto, Journal of the Physical Society of Japan 74 , 95 (2005).
- [29] J. S. Lim, M.-S. Choi, M. Y. Choi, R. L´ opez, and R. Aguado, Phys. Rev. B 74 , 205119 (2006).
- [30] J. S. Lim, R. L´pez, o and D. S´nchez, New Journal of a Physics 16 , 015003 (2014).
- [31] Y. Kleeorin and Y. Meir, Phys. Rev. B 96 , 045118 (2017).
- [32] A. Carmi, Y. Oreg, and M. Berkooz, Phys. Rev. Lett. 106 , 106401 (2011).
- [33] R. L´ opez, T. c. v. Rejec, J. Martinek, and R. Zitko, ˇ Phys. Rev. B 87 , 035135 (2013).
- [34] T. Kita, R. Sakano, T. Ohashi, and S.-i. Suga, Journal of the Physical Society of Japan 77 , 094707 (2008).
- [35] I. Kuzmenko and Y. Avishai, Phys. Rev. B 89 , 195110 (2014).
- [36] Y. Nishida, Phys. Rev. Lett. 111 , 135301 (2013).
- [37] J. Bauer, C. Salomon, and E. Demler, Phys. Rev. Lett. 111 , 215304 (2013).
- [38] Y. Nishida, Phys. Rev. A 93 , 011606 (2016).
- [39] I. Kuzmenko, T. Kuzmenko, Y. Avishai, and G.-B. Jo, Phys. Rev. B 93 , 115143 (2016).
- [40] C. Mora, P. Vitushinsky, X. Leyronas, A. A. Clerk, and K. Le Hur, Phys. Rev. B 80 , 155322 (2009).
- [41] R. Delagrange, J. Basset, H. Bouchiat, and R. Deblock, Phys. Rev. B 97 , 041412 (2018).
- [42] P. W. Anderson, Phys. Rev. 124 , 41 (1961).
- [43] H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. Rev. B 21 , 1003 (1980).
- [44] L. I. Glazman and M. E. Raikh, J. Exp. Theor. Phys. 27 , 452 (1988).
- [45] J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149 , 491 (1966).
- [46] O. Parcollet, A. Georges, G. Kotliar, and A. Sengupta, Phys. Rev. B 58 , 3794 (1998).
- [47] C. Mora, Phys. Rev. B 80 , 125304 (2009).
- [48] P. Nozi´res, J. Low Temp. Phys. e 17 , 31 (1974).
- [49] I. Affleck and A. W. W. Ludwig, Phys. Rev. B 48 , 7297 (1993).
- [50] D. L. Cox and A. Zawadowski, Advances in Physics 47 , 599 (1998).
- [51] D. B. Karki, C. Mora, J. von Delft, and M. N. Kiselev, Phys. Rev. B 97 , 195403 (2018).
- [52] D. B. Karki and M. N. Kiselev, Phys. Rev. B 98 , 165443 (2018).
- [53] L. V. Keldysh, Sov. Phys. JETP 20 , 1018 (1965).
- [54] T.-S. Kim and S. Hershfield, Phys. Rev. Lett. 88 , 136601 (2002).
- [55] T.-S. Kim and S. Hershfield, Phys. Rev. B 67 , 165313 (2003).
- [56] C. W. J. Beenakker and A. A. M. Staring, Phys. Rev. B 46 , 9667 (1992).
- [57] V. Zlatic and R. Monnier, Modern Theory of Thermoelectricity (Oxford University Press, 2014).
- [58] K. A. Matveev and A. V. Andreev, Phys. Rev. B 66 , 045301 (2002).
- [59] A. Carmi, Y. Oreg, M. Berkooz, and D. GoldhaberGordon, Phys. Rev. B 86 , 115129 (2012).
- [60] M. Pustilnik and L. I. Glazman, Phys. Rev. Lett. 87 , 216601 (2001).
- [61] M. Pustilnik and L. Glazman, Journal of Physics: Condensed Matter 16 , R513 (2004).
- [62] M. Pustilnik, L. Borda, L. I. Glazman, and J. von Delft, Phys. Rev. B 69 , 115316 (2004).
- [63] J. Park, S.-S. B. Lee, Y. Oreg, and H.-S. Sim, Phys. Rev. Lett. 110 , 246603 (2013).
- [64] A. Dorda, M. Ganahl, S. Andergassen, W. von der Linden, and E. Arrigoni, Phys. Rev. B 94 , 245125 (2016).
- [65] D. P´ erez Daroca, P. Roura-Bas, and A. A. Aligia, Phys. Rev. B 97 , 165433 (2018).
- [66] U. Eckern and K. I. Wysoki´ nski, arXiv e-prints , arXiv:1904.05064 (2019).
D. B. Karki Mikhail N. Kiselev
Nonlinear Seebeck effect of SU($N$) Kondo impurity
We develop a theoretical framework to study the influences of coupling
asymmetry on the thermoelectrics of a strongly coupled SU($N$) Kondo impurity
based on a local Fermi liquid theory. Applying non-equilibrium Keldysh
formalism, we investigate charge current driven by the voltage bias and
temperature gradient in the strong coupling regime of an asymmetrically coupled
SU($N$) quantum impurity. The thermoelectric characterizations are made via
non-linear Seebeck effects. We demonstrate that the beyond particle-hole (PH)
symmetric SU($N$) Kondo variants are highly desirable with respect to the
corresponding PH symmetric setups in order to have significantly improved
thermoelectric performance. The greatly enhanced Seebeck coefficients by
tailoring the coupling asymmetry of beyond PH symmetric SU($N$) Kondo effects
are explored. Apart from presenting the analytical expressions of asymmetry
dependent transport coefficients for general SU($N$) Kondo effects, we make a
close connection of our findings with the experimentally studied SU(2) and
SU(4) Kondo effects in quantum dot nano structures. Seebeck effects associated
with the theoretically proposed SU(3) Kondo effects are discussed in detail.