Nonequilibrium steady-state thermoelectrics of Kondo-correlated quantum dots

The transport across a Kondo-correlated quantum dot coupled to two leads with independent temperatures and chemical potentials is studied using a controlled nonperturbative, and in this sense numerically exact, treatment based on a hybrid numerical renormalization group combined with time-dependent density matrix renormalization group (NRG-tDMRG). In the Kondo regime, for sufficiently large fixed voltage bias $V\ensuremath{\gtrsim}{T}_{K}$, with ${T}_{K}$ the Kondo temperature, we find a peak in the conductance vs the temperature gradient $\mathrm{\ensuremath{\Delta}}T={T}_{R}\ensuremath{-}{T}_{L}$ across left and right lead. Focusing then on zero voltage bias but finite $\mathrm{\ensuremath{\Delta}}T$ far beyond linear response, we reveal the dependence of the characteristic zero-bias conductance on the individual lead temperatures. We find that the finite-$\mathrm{\ensuremath{\Delta}}T$ data behaves quantitatively similar to linear response with an effective equilibrium temperature derived from the different lead temperatures. The regime of sign changes in the Seebeck coefficient, signaling the presence of Kondo correlations, and its dependence on the individual lead temperatures provide a complete picture of the Kondo regime in the presence of finite-temperature gradients. The results from the zero-bias conductance and Seebeck coefficient studies unveil an approximate ``Kondo circle'' in the ${T}_{L}/{T}_{R}$ plane as the regime within which the Kondo correlations dominate. We also study the heat current and the corresponding heat conductance vs finite $\mathrm{\ensuremath{\Delta}}T$. We provide a polynomial fit for our numerical results for the thermocurrent as a function of the individual lead temperatures, which may be used to fit experimental data in the Kondo regime.