How to measure the entropy of a mesoscopic system via thermoelectric transport

TL;DR

This paper introduces a method to measure the entropy of mesoscopic systems through thermoelectric transport, providing analytical proof, numerical validation, and experimental application to quantum dots, revealing entropy changes and resolving longstanding puzzles.

Contribution

It proposes a novel technique to extract entropy from transport measurements in mesoscopic systems, applicable to systems with arbitrary spectra and degeneracies.

Findings

Method successfully applied to quantum dot thermoelectric data

Able to determine entropy changes across Coulomb-blockade valleys

Resolves the puzzle of finite thermoelectric response at particle-hole symmetry

Abstract

Entropy is a fundamental thermodynamic quantity indicative of the accessible degrees of freedom in a system. While it has been suggested that the entropy of a mesoscopic system can yield nontrivial information on emergence of exotic states, its measurement in such small electron-number system is a daunting task. Here we propose a method to extract the entropy of a Coulomb-blockaded mesoscopic system from transport measurements. We prove analytically and demonstrate numerically the applicability of the method to such a mesoscopic system of arbitrary spectrum and degeneracies. We then apply our procedure to measurements of thermoelectric response of a single quantum dot, and demonstrate how it can be used to deduce the entropy change across Coulomb-blockade valleys, resolving, along the way, a long-standing puzzle of the experimentally observed finite thermoelectric response at the apparent particle-hole symmetric point.

Figures

Equations

$$G(\mu ,T)=\sum _{ij}{G}_{ij}(\mu ,T)=\sum _{ij}{\mathcal{T}}_{ij}^{(0)}\times \left[(P_i^{(N+1)}(\mu ,T)+P_j^{(N)}(\mu ,T)\right]\frac{df(E_i^{(N+1)}-E_j^{(N)}-\mu ,T)}{d\mu }$$

$$G_{ij}(\mu ,T)=C(T)G_{ij}^{NI}(\mu +{\Delta }_{ij}(T),T)$$

$$\Delta_{ij}(T)=\frac{E^{(N+1)}_{j}-E^{(N)}_{i}}{2}+\frac{T}{2}\log\big{[}\frac{\sum_{j}g^{(N+1)}_{j}e^{-E^{(N+1)}_{j}/T}}{\sum_{i}g^{(N)}_{i}e^{-E^{(N)}_{i}/T}}\big{]}$$

$${\mathrm{TR}}_{ij}(\mu ,T)=C(T){\mathrm{TR}}_{ij}^{NI}(\mu +{\Delta }_{ij}(T),T)+G_{ij}(\mu ,T){\Delta }_{ij}(T)/T$$

$${\mathrm{TR}}^{NI}(\mu ,T)={\gamma }_{1}T\frac{d{G}^{NI}(\mu ,{\gamma }_{2}T)}{d\mu },$$

Full text

Introduction

The entropy of a mesoscopic system can yield nontrivial information on emergence of exotic states, such as two-channel Kondo impurity [4], non-abelian anyons in the $\nu=5/2$ regime [1, 2] or Majorana modes in topological superconductors [3]. Nevertheless, the measurement of entropy in such small electron number systems is highly nontrivial. Refs.  5 and 6 used the asymmetry of the in and out tunneling processes in a quantum dot (QD) to deternmine the degeneracy of the QD states, while recent elegant experiments [7] have employed the thermodynamic Maxwell relation between entropy evolution and chemical potential, $(\partial\mu/\partial T)_{n}=-(\partial S/\partial n)_{T}$, in order to directly measure entropy transitions in semiconductor QDs. This latter experiment required measurements of another thermodynamic quantity - the charge of the system as a function of gate voltage, for different temperatures, and hence a specially designed device. Here we propose a different approach to this problem: can one extract information about the entropy from transport measurements ? Obviously, this requires a measurement of both particle and thermal (entropy/heat) transport. This question has been addressed in the context of bulk solids [8, 9, 10, 11], with sometimes debated points of view. A general relation exists between the low-temperature thermopower and specific-heat (entropy) of a free electron gas, and this relation appears to apply in a number of materials[9, 10]. However, thermopower is, quite generally, a transport coefficient and its relation to entropy has been shown to be questionable in systems with strongly anisotropic transport for instance[11]. In the opposite high-temperature limit, where temperature is the largest energy scale in the system, general relations between the thermopower and derivatives of the entropy can be derived, embodied in the Heikes [12, 8, 13] and Kelvin [14, 11] formulas.

Consider an arbitrary mesoscopic system in the Coulomb-blockade regime (where only $N$ and $N+1$-particle states are energetically relevant), whose entropy one wishes to measure. The method we propose here is based on a general observation, which is also an important result of our work: if one weakly couples this system to leads, the conductance of such an interacting system can be put in the form of a non-interacting conductance formula, provided one takes into account a temperature-dependent shift of the chemical potential (gate voltage). The thermal response (TR), in turn, can be written in a similar manner, where the temperature-dependent shift in the chemical potential produces an extra contribution. We show that this extra term, which can be determined by comparing the actual thermal response of the system to that of the related non-interacting system (which can be estimated using a newly introduced high-temperature version of the original Mott formula[15]), can be used to extract the entropy of such a mesoscopic system even in the case of arbitrary spectrum and degeneracies, and then demonstrate the usefulness of the approach by applying it to several model systems. One big advantage of our formulation is that one can apply it to any such mesosopic system where measurements of both electrical conductance and thermopower are available. This allows us to apply our procedure to existing data of thermoelectric response of a single QD, and demonstrate how it can be used to deduce the entropy change and the QD’s degeneracy. In the process we explain the long standing puzzle of the observation of a non-zero thermopower at the apparent electron-hole symmetry point in the Coulomb Blockade (CB) valley [16, 17].

Results

General Formulation

Consider a general mesoscopic system with many-body eigenstates $\Psi^{(N)}_{i}$, where $N$ is the number of electrons in that state, with energies $E^{(N)}_{i}$ (with $g^{(N)}_{i}$ the degeneracy of the energy $E^{(N)}_{i}$), whose entropy one wishes to measure. In order to perturb the system as little as possible, we weakly couple the mesoscopic system to two reservoirs (with coupling $V_{i}$ for each state $i$). In this weak-coupling limit $\Gamma_{ij}=2\pi V_{i}V_{j}\rho$, the characteristic level broadening, with $\rho$ the density of states in the reservoirs, obeys $\Gamma_{ij}\ll T$, where $T$ is the temperature. In this limit the conductance $G$ through the mesoscopic system can be written as the sum of individual transitions from state $i$ with $N$ electrons to state $j$ with $N+1$ electrons [18]

[TABLE]

where $\mathcal{T}^{(0)}_{ij}$ is equal to $\Gamma_{ij}$ times the overlap of the $N+1$-particle many-body wave function $\Psi^{(N+1)}_{j}$ with the $N-$particle wavefunction $\Psi^{(N)}_{i}$ , with the addition of the electron tunneling in from the leads (or the reverse process) (see Supplementary Information, Eq. 1). In the above $f(E,T)$ is the equilibrium Fermi function, $\mu$ the chemical potential, and $P_{i}^{(N)}(\mu,T)=e^{-(E_{i}^{(N)}-\mu N)/T}/Z$ is the equilibrium probability of the system to be in the $N$-particle many-body state $i$, with $Z$ the partition function (except for the experimental part, we use $k_{B}=1$ throughout the paper, where $k_{B}$ is the Boltzman coefficient, so that temperature has units of energy and entropy is dimensionless). A similar expression can be written for the TR, defined as $dI/dT$, the change in the linear-response current due to temperature difference between the leads, in analogy to conductance, with $df/d\mu$ being replaced by $df/dT$. We assume that the Coulomb energy is significantly larger than $T$ and $\Gamma$ so that for a given chemical potential, $G$ involves transitions between states with only $N$ or $N+1$ particles. A crucial step in our formulation is the demonstration that the above general expressions for the conductance and the thermal response for an arbitrary interacting system can be accurately written, in the vicinity of each $N\rightarrow N+1$ transition, as those for a non-interacting system, but with a temperature-dependent effective chemical potential (see Supplementary Note 1):

[TABLE]

where $G^{NI}_{ij}$ is the conductance for a non-interacting system with same spectrum and couplings, and $C(T)$ is some temperature dependent prefactor, that will drop out when the relation between G and TR is derived. This temperature dependent shift in the chemical potential is given by

[TABLE]

In the simple case of a transition from an empty state into a single level, with degeneracy $g$, this shift reduces to $\frac{1}{2}T\log g$, which has been noticed before [19, 2], and has been measured experimentally [5]. In that case this shift was attributed to the fact the chemical potential has to shift in order to compensate for the fact there are $g$ ways for an electron to tunnel into the QD, while having a single channel for tunneling out, an asymmetry that has been verified experimentally [5, 20]. In contrast, our expression indicates that in the case of many levels, which has not been discussed before, the temperature-dependent part of the shift does not depend on which level the electron tunnels through, and what its degeneracy is. This part of the shift is identical for all transitions, and is equal one half of the difference of the canonical free energies between the CB valleys corresponding to $N$ and $N+1$ electrons.

The explicit dependence of $\Delta_{ij}$ on $T$ allows us to write, in a similar manner to Eq. 2, an explicit expression for the TR of a general interacting system in terms of its conductance and the TR of the related non-interacting system,

[TABLE]

In order to derive an equation for $\mathrm{TR}^{NI}$, the thermal response of a non-interacting system with same spectrum and couplings, we generalize the Mott formula [15], valid for $T\ll\Gamma$, to the regime $T\gg\Gamma$ (see Eq. 5 in the Methods section and Supplementary Note 2 for derivation). Thus, the deviation of the $\mathrm{TR}$ from $\mathrm{TR}^{NI}$ (calculated from the conductance) allows us to estimate $\Delta_{i,j}(T)$, and consequently the entropy difference between the consecutive CB valleys: $\Delta S_{N\rightarrow N+1}=2d\Delta_{ij}(T)/dT$.

So, given the experimentally or numerically obtained $G(\mu,T)$ and $TR(\mu,T)$, the procedure we propose for finding the entropy difference between consecutive CB valleys is the following: 1. Given $G(\mu,T)$, one can use our variant of the Mott formula (Eq. 5 in the Methods section) to evaluate the first term on the right-hand side (RHS) of Eq. 4. 2. For a given temperature, the difference between this term and the actual TR, which is a function of the chemical potential, is proportional to $G(\mu,T)$. We denote this proportionality constant A(T) (note that $A(T)$ is the only fitting number required, for a given temperature, to map the two functions on top of each other). 3. Given the obtained $A(T)$, the difference in entropy between the valleys is then given by $\Delta S_{N\rightarrow N+1}=2d\left[T\times A(T)\right]/dT$. (A step by step description of the fitting process is detailed in Supplementary Note 3).

In the following we demonstrate the usefulness of this formalism in model systems, where one can compare the entropy obtained using the above relation to that calculated directly from thermodynamic considerations, and finally we apply our formalism to available experimental data.

Comparison to numerical calculations

Let us start with a simple example where in each $N$-electron subspace there are $g^{(N)}$ degenerate $N$-particle states of energy $E^{(N)}$, and all other states can be ignored (i.e. the level spacing is much higher than temperature). In this case the entropy $S_{N}$ in each valley is equal to $\log g^{(N)}$, and is temperature independent. In this case, one indeed finds that the proportionality constant is temperature independent, $A(T)=\log(g^{(N+1)}/g^{(N)})/2$. Fig. 1b illustrates the correspondence between the TR obtained directly, using Eq. 1, and that obtained by the RHS of Eq. 4 (red circles), for a four-fold degenerate interacting QD, relevant, for example, to a carbon nanotube QD (see also experimental section below). The conductance used in evaluating both terms in the RHS of Eq. 4 was also obtained via Eq. 1 (and is shown in Fig. 1a). In this case there are 4 CB peaks, separating valleys with degeneracies $g^{(N)}=1,4,6,4$ and $1$ for $N=0,1,2,3$ and $4$. In order to construct the estimate for the TR in Fig. 1b we have used the above fitting procedure separately for each peak, as the entropy difference between consecutive valleys is different for each peak. The figure displays an almost perfect agreement between the direct calculation of the TR and that obtained by our Ansatz.

In this case, as the entropy change $\Delta S$ between the valleys is temperature independent, the estimate of $A$ at a single temperature is directly proportional to the entropy change through $\Delta S=2A$. In particular, the entropy change across the first CB peak is a direct measure of the degeneracy of the QD ($4$ in the above example). We have repeated the procedure for QDs of arbitrary degeneracy. Fig. 1c depicts the entropy change deduced using our procedure (red circles), compared to the expected change in entropy($\log g^{(N+1)}$). We see a perfect agreement even up to large degeneracies. As mentioned above, some aspects of this simple case of a single degenerate level have been addressed before, and it has been suggested that the thermopower through a single-level QD can be used, e.g., to deduce the nature of the neutral modes in the fractional quantum Hall regime [2].

The advantage of our procedure lies in its application to a multi-level mesoscopic system, such as a multi-level QD, or to a multi-dot system, where the entropy is temperature dependent. As an example, let us consider the case of two singly degenerate levels, with level spacing $\Delta\epsilon$ (describing, for example, a single-level QD in a magnetic field). One expects that when $T\ll\Delta\epsilon$ the entropy of the single-electron system will be equal to zero, while for higher temperature, larger than $\Delta\epsilon$, it will increase to $\log 2$. As the entropy is temperature dependent, one has to perform the procedure for all $T$ in order to extract $A(T)$, its derivative, and consequently the entropy. For simplicity, we assume that the transition through one of the levels dominates the transport, so Eq. 4, which corresponds to a transition between specific states, will also reflect the full transport coefficient of the system. As we will demonstrate, even though a single transition dominates the transport, the resulting procedure yields the full entropy change in the system.

Fig. 1d and e depict, respectively, the calculated conductance and TR, again using Eq. 1, for a specific temperature, $T=\Delta\epsilon$. Fig. 1e also shows the TR derived from our procedure - the fitting leads to $A(T=\Delta\epsilon)$ for this temperature. Repeating the same procedure for many temperatures, one is able to produce the whole curve $A(T)$, and then the entropy change, $\Delta S=2d\left[TA(T)\right]/dT$. The resulting estimate for the entropy change is plotted in Fig. 1f along with the thermodynamic calculation of the entropy change: $\Delta S_{N\rightarrow N+1}=-\partial\left[F_{N+1}(T,\mu)-F_{N}(T,\mu)\right]/\partial T$ with $F_{N}(T,\mu)$ the free energy of the $N$-electron system. Again we observe excellent agreement between the entropy deduced in our procedure and the direct calculation. In Supplementary Note 4 we discuss our procedure for the case when several transitions are relevant to the total transport.

Interestingly, while this formalism was derived for the weak coupling ($\Gamma\ll T$) regime, empirically its validity extends outside this strict regime. Since Eq. 1 does not apply to the regime $\Gamma\gtrsim T$, we have employed here the numerical-renormalization-group (NRG) method (see Methods), which is accurate down to zero temperature. Fig. 2 demonstrates the validity of our formalism and shows that the estimates of the entropy, using our procedure for the cases of a two-fold ($SU(2)$) and four-fold ($SU(4)$) degenerate single-level QD, agree with expected values ($\log 2$ and $\log 4$, respectively), down to $T\simeq 0.1\Gamma$. The fitting procedure that corresponds to Eq. 4 remains accurate throughout the presented region of temperatures with coefficient of determination ($R^{2}$) values of close to unity (crosses in Fig. 2c,d). Thus, at least for these two models, our approach extends to couplings to the leads $\Gamma$ which are of the order or even larger than temperature.

Application to Experiments

One of the main advantages of our approach, compared, e.g. to that of ref. 7, is that it can be readily applied to any previous transport experiment in a mesoscopic system, for which conductance and TR data are available. As an example of the usefulness of the suggested procedure, we have analyzed recent thermoelectric measurement results[21] through a QD device, formed in a two dimensional electron system of a GaAs/AlGaAs heterostructure using split-gate technology. This technology allows for a high degree of control over system parameters such as QD energy and tunnel coupling $\Gamma$ between the QD and the reservoirs, by adjusting the voltages applied to the split gates. The sample is shown in the inset to Fig 3b. Gates B1, B2 and B3 are used to form the QD (yellow dot). The tunnel coupling between the QD and the reservoirs H and C can be controlled symmetrically adjusting the gate voltage applied to gate B1. Gate P, the so-called plunger gate, is used to continuously tune the electrochemical potential of the QD, and consequently the number of electrons on the QD. Gate G is not used in these experiments and is kept at ground at all times.

The sample is cooled down in a dilution refrigerator, with an electron base temperature of $\approx 230$ mK, in the presence of a small perpendicular magnetic field (B = 0.6 T) [22]. In order to establish a temperature difference $\Delta T$ across the QD, a small heating current was applied to reservoir H (see Methods section and Supplementary Note 5), thereby mainly enhancing the electron temperature in that reservoir. The thermovoltage $V_{th}$ is then obtained by recording the voltage drop across the QD as a response to the temperature increase in reservoir H under open circuit conditions (see methods section and Supplementary Note 5 for further details), thus $V_{th}=\mathrm{TR}\times\Delta T/G$.

Fig. 3a and b depict the experimental data for $G$ and $V_{th}$, respectively, for a pair of CB peaks. Interestingly, the data show that at points of apparent particle-hole symmetry in the conductance (e.g. arrow in Fig. 3b and crossing point in Fig. 4b), $V_{th}$ does not vanish as would be expected from the usual, spin-degenerate QD, described by the standard single-impurity Anderson model [23]. This experimental observation (see also Refs.  16, 17) is to this day an unresolved puzzle in the field (see Ref.  31 for an attempt to resolve this puzzle).

In the following we detail our analysis of these CB peaks. It has been noted before [17] that under the condition of heating one reservoir, the actual temperature of the QD can differ greatly from the fridge’s temperature. Since in the present case where $T<\Gamma\simeq 550\mu eV$, the actual temperature cannot be deduced from the width of the CB peaks, we use the temperature as an additional fitting parameter. In addition, since the x-axis relation between the conductance measurement (Fig. 3a) and the thermovoltage measurements (Fig. 3b) were not experimentally established, another fitting parameter is introduced: the x-axis shift in the measured conductance relative to the measured thermovoltage. The results of fitting the TR to Eq. 4 are depicted in Fig. 3c. As can be seen in the figure, there is a good agreement between the fit and the observed TR in the vicinity of each peak, again using only a few fitting parameters to fit the whole curve (see Supplementary Note 3 for a detailed step-by-step of the analysis of the experimental data using our procedure), illustrating the experimental validity of our approach. Due to the limited availability of the data we used $G(\mu,T)$ instead of $G(\mu,\gamma_{2}T)$ to estimate $\mathrm{TR}^{NI}$. However, this should make a little difference when $T<\Gamma$.

In applying our method to the experiment, one needs to translate the measured $V_{th}$ to the thermo-electric response TR by dividing by $\Delta T$. This value, however, is not easily and accurately determined in an experiment and thus leads to uncertainties in the absolute values of the entropy changes across the peaks. On the other hand, the ratio of these entropy changes across consecutive peaks is independent of $\Delta T$, and is found to be $-2.07\pm 0.13$ for the two peaks depicted in Fig. 3 (the errors estimate is due to variation in possible fitting region around the peaks, see Supplementary Note 3). The simplest scenario giving rise to such a ratio, is that the entropy change across the first peak is $\log 4$ while the second is $-\log 2$. This means that the first peak signals a transition into a four-fold degenerate state, while the second peak may either correspond to a transition from a four-fold degenerate to a two-fold degenerate state, or from a two-fold degenerate state to a non-degenerate state. This suggests a deviation from the naive picture of consecutive filling of a four-fold degenerate state. Including this scenario into our fit, $\Delta T$ is found to be $\approx 20mK$, which is close to the experimental estimate of being of the order of 30 mK (see methods and Supplementary Note 5).

While the degeneracy of these two levels seems fortuitous, such a model, in fact, has been claimed to be generic for transport through QDs [24, 25, 26], and has been invoked to explain the repeating phase jumps in the transmission phase through such a dot[27, 28]. In these works this is caused by two overlapping levels with different tunneling widths. At each conductance valley the narrow level is filled by an additional electron, shifting the energies of the narrow and the wide level differently, thus leading naturally, due to the degeneracy, to the entropy change of $\log 4$ across the first peak. In this scenario, after the second conductance peak the narrow level is doubly occupied, and does not play an additional role in transport, while the wide level is shifted up to overlap with another narrow level, and the process repeats itself. This explained the repeated phase change across consecutive conductance peaks[27, 28], and is, in fact, consistent with the observation that the upshift of the TR from zero at the apparent particle-hole symmetric point happens in consecutive pairs of conductance peaks [16].

Experimentally, one can easily change the tunneling rates $\Gamma$ between the QD and the leads through the split gate technique. These data, depicted in Fig. 4a and b, can then be used to differentiate between these possible scenarios. We found that the model that best reproduces the experimental findings, is that of a QD with two spinful states with an energy difference $\Delta\epsilon$ that depends on gate voltage (in the model we used the same $\Gamma$ for both levels to avoid additional parameters). Similar evolution of the degeneracy as a function of chemical potential has already been observed in quantum nano-tubes[29].

In this model, around the gate voltage corresponding to the first peak (QD energy $\epsilon\sim-0.75$meV), the two levels are almost degenerate yielding a net four-fold degeneracy($g_{N}=0,g_{N+1}=4$) which is lifted as the gate voltage is tuned toward the second peak, around QD energy $\epsilon\sim 0.75$meV ($g_{N}=2,g_{N+1}=1$) (as illustrated in the inset of Fig. 4d). This interpretation leads to the observed values of entropy change.

Fig. 4c and d depict NRG calculation of a specific model for various values of $\Gamma$, where the energy difference between the levels changes linearly with chemical potential, $\Delta\epsilon=a+b(\mu-\epsilon)$, with $a=-0.01D,b=0.13$ ($D$ the bandwidth of the leads). The model reproduces the essential experimental features and those captured by varying $\Gamma$. Some features in the experimental data, such as small side peaks in the lower two values of $\Gamma$, attributed to excited states [30], are not captured within the current simple model. Interestingly, this model naturally reproduces the non-zero value of the TR at the seemingly particle-hole symmetric point which is also visible in the experimental data (crossing point in Fig. 4b, marked by an arrow). This anomalous increase of the TR around the middle of the valley is attributed to a non-trivial degeneracy, thus providing a natural explanation that this value of gate voltage does not correspond, in fact, to a particle-hole symmetric point. (An alternative explanation, based on non-linear effects, was suggested in recent work[31].)

Discussion

In this work, we have derived a theoretical connection between the entropy and transport coefficients in mesoscopic junctions. This connection relates the TR of a Coulomb-blockaded mesoscopic system with arbitrary many-body levels to the conductance and the entropy change between adjacent CB valleys. While the derivation was introduced for weak coupling $\Gamma$ between the system and the leads (in comparison with temperature), we have demonstrated numerically that, for the case of 2-fold and 4-fold degenerate QD, the method is accurate also for temperatures well below $\Gamma$. This allowed us to apply the method to experimental data in that regime, which yielded non-trivial, and in fact unexpected information about the entropy in each CB valley. The deduced theoretical model, which described the experimental QD, reproduced the measured thermopower and resolved the long-standing puzzle of a finite TR in the apparent particle-hole symmetric point.

The success of this procedure suggests possible venues to extend this analysis especially towards the study of entropy of exotic states. One direction would be to extend the method to low temperatures, thus enabling the determination the degeneracy of the ground state of the full system. This, for example, is particularly relevant to exotic phases, such as the two-channel Kondo system, where the zero temperature entropy is non zero. If the TR of this system can be utilized to deduce the entropy of the ground state, this can be a smoking gun for the observation of the two channel Kondo ground state [32] or other such non Fermi liquid ground states. Such an extension has also been suggested in parallel by Sela et al. [37] to measure the fractional entropy of Majorana zero modes.

Methods

High Temperature Mott Relation

In relating the non-interacting conductance and TR we use a high temperature adaptation of the Mott relation[15].

[TABLE]

where the superscript $NI$ denotes a non-interacting system, and $\gamma_{2}=2/\sqrt{3},\gamma_{1}=2\gamma_{2}^{3}$ are universal values related to properties of the Fermi function (for derivation see Supplementary Note 2).

Numerical Renormalization Group

for the density-matrix numerical renormalization group (DM-NRG) results we used the open-access Budapest Flexible DM-NRG code[34, 35]. The expectation values and the transmission spectral function, required for the evaluation of the conductance through the double dot device [18], were calculated, assuming, for simplicity, equal couplings to the left and right leads, $\Gamma=\pi\rho V^{2}$, and equal and constant density of states $\rho=1/2D$ in the two leads, with a symmetric band of bandwidth $2D$, around the Fermi energy. The NRG simulation is able to output the many body discreet energy states that the system can occupy and their respective spectral weight, $\epsilon_{i},w_{i}$. Transport coefficient are then calculated using $G(\mu,T)=\Gamma\pi\sum w_{i}df(\epsilon_{i}-\mu,T)/d\mu$ and $\mathrm{TR}(\mu,T)=\Gamma\pi\sum w_{i}df(\epsilon_{i}-\mu,T)/dT$.

Experiment

Our sample is designed similar to the one used by Scheibner et al. [16]. The electron reservoir H which serves as a hot lead for the quantum dot in our thermopower experiments is shaped into a channel of width $w=2\mu m$ and length $l=20\mu m$ (see Supplementary Fig.  5). The QD is situated on one side of the channel, delimited by gates B1 and B2,while the opposite side of the channel is delimited by the two gates Q1 and Q2, forming a quantum point contact (QPC) which is positioned exactly opposite to the quantum dot. The QPC is adjusted to the conductance plateau at G = 10 $e^{2}/h$. It separates the heating channel H from the reservoir REF which is kept at ground potential. At the two ends of the heating channel (separated by the distance $l=20\mu m$) the 2DES opens up quickly into large reservoirs. The channel can be contacted electrically through two Ohmic contacts $I_{1}$ and $I_{2}$. We apply a heating current $I_{h}=70nA$ to the channel, which is modulated at a low frequency $\omega=13$ Hz. Because at low temperature electron-electron scattering is the dominant scattering mechanism on length scales up to several 10 $\mu$m in our system, the power $P_{h}$ introduced through $I_{h}$ is dissipated inside the channel only into the electron gas while in the larger reservoirs outside the channel, $P_{h}$ is dissipated into the lattice through electron-lattice interaction. From here the heat gets removed efficiently by the dilution refrigerator. In this manner we establish a locally enhanced electronic temperature in the channel while the rest of the 2DES remains approximately at base temperature. Using the thermopower of the QPC as a thermometer [36] we estimate that for the given $I_{h}$, $T_{el}$ in the channel increases by $\Delta T\approx 30$ mK. We note that because $I_{h}$ gets modulated with $\omega$, the temperature in the heating channel oscillates with $2\omega$ since the dissipated power $P_{h}\propto I_{h}^{2}\propto sin^{2}(\omega t)\propto cos(2\omega t)$. This provides all temperature driven effects with a clear signature of an oscillation frequency of $2\omega$. The thermovoltage $V_{th}$ of the QD is obtained by measuring the potential difference between the contacts of the two cold reservoirs $V_{ref}$ and $V_{C}$ using a Lock-In amplifier operating at $2\omega=26Hz$. Since the QPC is adjusted to a conductance plateau its contribution to the $V_{th}$ is zero. Hence the measured signal can be attributed fully to the QD. In order to suppress any potential fluctuations at $\omega$ in close vicinity to the QD structure, which may occur due to unwanted capacitive coupling inside the sample, we let the excitation voltage for the heating current at both contacts of the heating channel oscillate symmetrically with respect to ground. Since reservoir REF is kept grounded, this suppresses oscillations of the electrical potential at $\omega$ around the QD structure.

Data Availability

The datasets generated and analysed in the study are available upon request from the corresponding authors.

Acknowledgements

We thank L. Maier for sample fabrication and D. Reuter and A.D. Wieck (Ruhr-University Bochum) for providing the hetero structure. The work in Würzburg has been funded by the Deutsche Forschungsgemeinschaft DFG (SPP1386). We thank P. Moca and G. Zarand for scientific discussions and help with the DM-NRG code. Y.M. acknowledges support from ISF grant 292/15. HT acknowledges funding through the European research council (Grant No. 339306, METIQUM). AG and YM acknowledge the KITP at the University of California at Santa Barbara, where this project has initiated, supported in part by NSF Grant No. PHY17-48958, NIH Grant No. R25GM067110, and the Gordon and Betty Moore Foundation Grant No. 2919.01.

Author contributions statement

AG, YM and LM initiated the project. AG, YK and YM developed the theory. YK performed the numerical calculations, and the experiments were performed by HT, HB and LM. All authors contributed to the writing of the paper.

Competing interests

The authors declare no competing interests.