Power, Efficiency and Fluctuations in a Quantum Point Contact as Steady-State Thermoelectric Heat Engine

The trade-off between large power output, high efficiency and small fluctuations in the operation of heat engines has recently received interest in the context of thermodynamic uncertainty relations (TURs). Here we provide a concrete illustration of this trade-off by theoretically investigating the operation of a quantum point contact (QPC) with an energy-dependent transmission function as a steady-state thermoelectric heat engine. As a starting point, we review and extend previous analysis of the power production and efficiency. Thereafter the power fluctuations and the bound jointly imposed on the power, efficiency, and fluctuations by the TURs are analyzed as additional performance quantifiers. We allow for arbitrary smoothness of the transmission probability of the QPC, which exhibits a close to step-like dependence in energy, and consider both the linear and the non-linear regime of operation. It is found that for a broad range of parameters, the power production reaches nearly its theoretical maximum value, with efficiencies more than half of the Carnot efficiency and at the same time with rather small fluctuations. Moreover, we show that by demanding a non-zero power production, in the linear regime a stronger TUR can be formulated in terms of the thermoelectric figure of merit. Interestingly, this bound holds also in a wide parameter regime beyond linear response for our QPC device.


Introduction
Nanoscale thermodynamics has attracted considerable attention during the last three decades. Key motivations are the prospect of on-chip cooling and power production as well as an enhanced thermoelectric performance arising from unique properties of nanoscale systems, such as quantum size effects and strongly energy-dependent transport properties [1][2][3][4][5][6][7][8][9]. Among various nanoscale systems, quantum point contacts (QPC) [10] are arguably the simplest devices which show a thermoelectric response [11]. A requirement of such a response is an energy-dependent transmission probability [12,13], which breaks the electron-hole symmetry. Within non-interacting scattering theory, the transmission probability fully determines the thermoelectric response of a two-terminal device. The QPC and similar devices provide a particularly interesting thermoelectric platform as their transmission probability may approximate a step function, maximizing the power generation [14,15]. This feature is in contrast to the case of a quantum dot, where the transmission probability may approximate a Dirac delta distribution, maximizing the efficiency of heat-to-power conversion [16][17][18][19][20].

Model System and Transport Theory
We consider the two-terminal setup shown in Figure 1, with a single-mode QPC connected to a left (L) and a right (R) electronic reservoir, characterized by electrochemical potentials µ L = µ 0 − eV L and µ R = µ 0 − eV R , and kept at temperatures T L = T 0 (cold reservoir) and T R = T 0 + ∆T (hot reservoir), respectively. Here, V L and V R are externally applied voltages, µ 0 denotes the electrochemical potential in the absence of voltage bias, T 0 corresponds to the background temperature and ∆T ≥ 0 stands for the temperature difference due to heating of the right reservoir. In the following, we always set µ 0 as the reference energy.

Quantum Point Contact
We employ the established model for a QPC [12] and describe the energy-dependent transmission probability as This is a step-like function of the energy E, see Figure 1b, where E 0 and γ denote the position and width in energy of the step, respectively. For a vanishingly small width, γ → 0, the transmission probability reduces to a step function, D(E) → θ(E − E 0 ). In experiments with 2DEGs, the width or smoothness of the QPC barrier γ, typically takes values of the order of 1 meV (corresponding to temperatures of the order of 10 K) [11,32,48,49].
The results presented in this paper are equally valid for different types of conductors, where the transmission function has a (smooth) step-like behavior, such as quantum wires with interfaces or controlled by finger gates. Here, smoothness parameters γ of values down to several µeV are expected (corresponding to temperatures of the order of 10-100 mK) [29].  (1), with a step positioned at energy E 0 and energy-smearing width, or smoothness, γ. The solid line shows the limit of vanishing width γ → 0.

Non-Linear-Transport Theory
The transport properties of the system are described by scattering theory [21]. In the non-lineartransport regime, the scattering properties of the QPC become dependent on the applied voltages V L and V R of the reservoirs and V g of the QPC-defining split gate [50] and possibly also on the temperature bias [5,35,51]. Since the details of this dependence will not be of importance for our analysis, we for simplicity consider a basic model with the QPC-potential capacitively coupled with equal strength, C L = C R = C, to the two terminals L and R. This leads to a modification of the transmission probability as where C g is the split gate-QPC capacitance. In the following, we absorb the gate potential dependence into the step energy E 0 + eC g V g /(2C + C g ) → E 0 (V g ) ≡ E 0 . This modification of the transmission probability guarantees a gauge-invariant formulation of the problem with observable quantities only dependent on the potential differences V = V L − V R , V L − V g and V R − V g . We here refrain from including the effect of a large temperature difference in the treatment of the transmission probability D(E), which is not required by fundamental principles such as gauge invariance and which has been little addressed so far, and postpone its study to future work.
For the study of the average currents of interest, namely charge current, I α , and heat current, J α , we now consider a symmetric biasing V L = −V R = V/2. We can then write the average currents that are flowing out of reservoir α as and Here, α should be understood as follows: L = R and R = L, whereas τ L = 1 and τ R = −1.
In Equations (2) and (3), we have introduced the Fermi distribution functions f α (E), While current conservation ensures I L = −I R ≡ I, energy conservation results in J L = −J R − IV.
To analyze the fluctuations in the system we also need the zero-frequency charge-current noise, given by [52] In addition to the study of the noise, it is often convenient to analyze the Fano factor being a measure of how much the noise deviates from the one of Poissonian statistics (for which F = 1).

Thermodynamic Laws and Performance Quantifiers
The laws of thermodynamics set very general constraints on the quantities introduced above and on the performance quantifiers, which we are going to study in this paper. We describe these quantities within scattering theory, known to correctly reproduce the laws of thermodynamics [5]. The first law of thermodynamics guarantees energy conservation and can be written as Here, we have introduced the electrical power produced, where −V I > 0 if the current flows against the applied bias. Please note that throughout the work, we limit our analysis of performance quantifiers to the relevant regime of positive power production. The second law of thermodynamics states that the entropy production σ is non-negative. In our two-terminal geometry, it can be written as This expression determines the direction of energy flows through the system. It equals zero in case that a process is reversible.
To determine the performance of the QPC as a heat engine, we now consider three independent quantities and combine them with each other. The first performance quantifier is given by the electrical power, Equation (8), which following the first law, Equation (7), is fully produced from heat.
The second performance quantifier we consider is given by the efficiency where J R is the heat current that flows out of the hot reservoir. As long as power is positive, the efficiency is bounded by the second law of thermodynamics, Equation (9), where η C denotes the Carnot efficiency. The dissipation arising from an inefficient heat to work conversion is quantified by the entropy, which thereby relates efficiency and produced electrical power to each other It is desirable to have a thermoelectric heat engine which not only produces large power, at high efficiency, but also minimizes fluctuations. The third independent performance quantifier is therefore provided by the low-frequency power fluctuations Interestingly, a trade-off between these quantities in the form of a TUR usually exists, as discussed in more detail in Section 6. This trade-off is typically written in the form of [38,39] where we have introduced the coefficient Q TUR . While this inequality is not always fulfilled for systems well described by scattering theory, see e.g., the discussion in [45,53], we find it to be respected in our system for all parameter values. Importantly this coefficient can be cast into the form [42] where we used Equations (12) and (13). Thus, under the constraint of positive power production and efficiency, we identify Q TUR as a convenient combined performance quantifier, accounting for power production, efficiency and power fluctuations together, where 1/2 sets the optimum value.

Linear-Response Regime
To compare with the much more studied linear-transport regime, we here present the relevant transport properties in this limit. Specifically, with small applied voltage and thermal bias, we can write the heat and charge current in the convenient matrix form [21], where (only due to linear response!) J = J L = −J R , and the matrix elements are defined as with Here f 0 (E) is the Fermi-Dirac distribution in Equation (4) with V α = 0 and T α = T 0 . In the same limit, the charge-current noise reduces to the equilibrium noise, given by S I = 2k B T 0 G, in accordance with the fluctuation-dissipation relation. Another performance quantifier, which is often used in the linear response, is the figure of merit ZT. It is given by [5] in terms of the response coefficients given above.

Power Production
To characterize the performance of the engine, we first consider the power P produced. The power as a function of applied bias V, for different values of the step energy E 0 and temperature difference ∆T, is shown in Figure 2 for both sharp (γ → 0) and smooth (γ = k B T 0 ) transmission step.
Sharp barrier (γ = 0) As is seen from the figure, a common feature for all P-vs-V curves is that they first increase monotonically from P = 0 at V = 0 with increasing negative voltage. At some voltage V max the power reaches its maximal value, P V max , and then decreases monotonically to zero, reached at the stopping voltage V s . The maximum power with respect to voltage is a function of E 0 /(k B T 0 ), ∆T/T 0 and γ/(k B T 0 ), i.e., In addition, we note that in the linear-response regime we have Figure 2, it is clear that the power as a function of voltage depends strongly on all parameters E 0 /(k B T 0 ), ∆T/T 0 and γ/(k B T 0 ). In particular, going from the linear to the non-linear regime by increasing ∆T/T 0 , the maximum power P V max might increase or decrease depending on the step properties γ and E 0 .

Maximum Power
To further analyze the properties of P V max , we first recall from the seminal work of Whitney [14,15,54] that the power is bounded from above by quantum mechanical constraints. It was shown that the upper bound is reached for a QPC with a sharp step, γ → 0, for which, using Equations (2) and (8), the power becomes Maximizing this expression with respect to eV/(k B T 0 ) and E 0 /(k B T 0 ) we find that the maximizing voltage is given by eV max = −ξk B ∆T where ξ ≈ 1.14 is the solution of ln(1 + e −ξ ) = −ξe −ξ /(1 + e −ξ ) [5]. Moreover, the maximizing step energy E 0,max and temperature difference ∆T max are related via [14,28] Inserting this expression, together with the relation for the maximizing voltage, into Equation (21) we reach the upper bound for the power established by Whitney [54] and related to the Pendry bound [55], which, we emphasize, holds in the linear as well as in the non-linear regime. To relate to this upper bound, in Figure 3a-c we present a set of density plots of P V max as a function of E 0 /(k B T 0 ) and ∆T/T 0 for different values of step smoothness parameters γ. From the figure it is clear that for a sharp step, γ → 0, there is a broad range of E 0 /(k B T 0 ) and ∆T/T 0 around the dashed line in the (E 0 , ∆T)-space, given by Equation (22), for which P V max is close to the theoretical maximum value P W . For a step smoothness up to γ ∼ k B T 0 , the situation changes only noticeably for small ∆T/T 0 . This is illustrated clearly in the close-ups in Figure 3d,e. Increasing the smoothness further, the region with maximum power close to P W shifts to higher values E 0 and ∆T, although still largely centered around Equation (22), as is shown in Figure 3f.
To provide a more quantitative analysis of this behavior, below we investigate two limiting cases for γ in further detail.

Small Smoothness Parameter
In the limit, where the value of the smoothness parameter γ is small, γ/(k B T 0 ) 1, the expression for the transmission probability in Equation (1) can be expanded to leading order in γ as [56] Inserting this into the expression for the charge current, Equation (2), and performing a partial integration for the delta function derivative, we get the power with P sharp given in Equation (21). To estimate how the overall maximum power is modified due to finite smoothness we insert into Equation (25) the values for eV/(k B T 0 ), E 0 /(k B T 0 ) and ∆T/T 0 along the line in the (E 0 , ∆T)-space, see Figure 3, which gives the bounded power for the sharp barrier. We find noting that E 0,max and ∆T are related via Equation (22). This expression quantifies the effect of the barrier smoothness visible in Figure 3, namely that the maximum power P V max in the region along the line in the (E 0 , ∆T) plane defined by Equation (22) is mainly affected for small ∆T/T 0 , and approaches P W in the strongly non-linear regime, ∆T/T 0 1.

Smoothness
Also in the case where the barrier gets smoother, such that γ equals the base temperature, γ = k B T 0 , we can perform an analytical investigation. Focusing on the linear-response regime, ∆T/T 0 1, where the effect of the smoothness is most pronounced, we can write the power in a compact form as As discussed above, in the linear-response regime the maximizing voltage is V max = V s /2, where the stopping voltage V s is the voltage that makes the expression in the curly bracket in Equation (27) vanish. Further maximizing over E 0 then gives E 0,max = 1.6 k B T 0 , which inserted into the power expression gives From Figure 3 it is clear that both E 0,max and P V,E 0 max are in good agreement with the numerical result.

Efficiency
Taking into account the aspect of limited resources, the power output is often not the most significant performance quantifier. A more relevant quantity is then the efficiency of a device. For a heat engine, it is defined as the power output divided by the heat absorbed from the hot bath, Equation (10).
We show the efficiency of the QPC as a steady-state thermoelectric heat engine in Figure 4. Panels (a) to (d) show the efficiency for the sharp barrier as function of voltage eV/(k B ∆T) for different temperature differences ∆T/T 0 and step energies E 0 /(k B T 0 ). For small absolute values of the step energies, see panels (a) and (b) for two examples with E 0 /(k B T 0 ) = −1, 0, the efficiency is rather small with respect to the Carnot efficiency, η/η C 0.25 and its overall shape only weakly depends on the temperature difference. This is radically different for larger values of E 0 : panels (c) and (d) of Figure 4 show a strong increase of the efficiency, which for E 0 /(k B T 0 ) = 15 and large temperature differences can reach about 90% of the Carnot efficiency. Also, the stopping voltage V s , at which the efficiency is zero and the device stops working as a thermoelectric, is strongly increased, depending on the temperature difference.  For large E 0 , see panel (d) of Figure 4, and small temperature differences, where large maximum efficiency values are reached, the efficiency-voltage relation takes a close-to-triangular shape. In this regime, we have that E ± eV/2 T 0 , T 0 + ∆T for all energies above the step energy E 0 . Therefore, only the tails of the Fermi functions contribute in Equations (2) and (3) and the efficiency in linear response in ∆T can be approximated as This formula describes well the triangular shape of the curves in panel (d), including the stopping voltage at small ∆T and large E 0 , given by eV s /k B ∆T ≈ −E 0 /k B T 0 , from the argument of the Heaviside function θ in Equation (29). We note that for V → V s the efficiency η → ∆T/T 0 ≈ η C , i.e., the efficiency approaches the Carnot limit, see Equation (11). The mechanism for this is the same as described in Ref. [20]; transport effectively takes place in a very narrow energy interval around E 0 , where the distribution functions f L (E 0 ) ≈ f R (E 0 ). Panels (e) to (h) of Figure 4 show results for the changes in the efficiency for a smooth barrier, γ = k B T 0 . At temperature differences that are much larger than the smoothness-here the case for k B ∆T/γ = 5, 15-the results for the efficiency are very similar to the case of the sharp barrier. This agrees with the discussion on the power production in the previous section, Section 3. At small temperature differences, however, the efficiency gets strongly reduced by the effect of the smoothness. This is particularly striking for large step energies, see panels (g)-(h) for E 0 /k B T 0 = 5, 15, respectively. Here, efficiencies that were close to Carnot efficiency for a sharp barrier get reduced by a factor three due to the barrier smoothness. The reason is that increasing smoothness leads to a broadening of the energy interval where the transport takes place, and hence a breakdown of the mechanism for Carnot efficiency discussed in Ref. [20].

Maximum Efficiency
We now focus our study on the maximum value of the efficiency that can be reached over the whole range of voltages, η V max , a function of E 0 /(k B T 0 ), ∆T/T 0 , γ/(k B T 0 ). The results of this maximization procedure are shown in Figure 5, where panel (a) corresponds to a sharp barrier (γ = 0) while smooth barriers with γ = k B T 0 and γ = 3k B T 0 are presented in panels (b) and (c), respectively. Two important results can be immediately seen from these density plots of the efficiency η V max as a function of temperature difference ∆T/T 0 and step energy E 0 /(k B T 0 ). First, we confirm the observations about the response to small temperature differences made from Figure 4. While for a sharp barrier, efficiencies close to Carnot efficiency are reached in the linear response (close to the stopping voltage, as we know from Figure 4), for even only slightly smoothed barriers this is not the case anymore. For γ/k B T 0 = 3, the maximum efficiency in the linear response is even suppressed down towards zero. This clearly shows that whenever the barrier step is not truly sharp, non-linear response is required to get a thermoelectric response with large power output and with high efficiency. Second, panels (b) and (c) of Figure 5 show that for temperature differences much larger than the smoothness-or, in other words, with one of the reservoir temperatures being much larger than the smoothness-almost the same (large) efficiency as in the sharp-barrier case is found, as long as the step energy is sufficiently large. Note, however, that these large-efficiency regions are far from those regions, which were previously identified as the ones of large power output, and are furthermore limited to regions with very large temperature differences and step energies.

Power-Efficiency Relations
The relation between power and efficiency for a sharp barrier, γ = 0, was investigated in detail in Refs. [14,15,28,36]. A convenient way to present the efficiency at a given power output, and vice versa, is in the form of lasso diagrams, as shown in Figure 6. At small step energies, E 0 /(k B T 0 ) = −1, 0, the maximum power as well as the maximum efficiency are relatively small. However, maximum efficiency and maximum output power basically happen at the same parameter values. This is advantageous for operation of a thermoelectric device, where one typically must decide whether to optimize the engine operation with respect to efficiency or power output.
This trend continues also for larger step energies, see panels (c) and (d) of Figure 6, as long as the temperature difference is larger or of the order of the step energy, k B ∆T E 0 (meaning that T 0 + ∆T > E 0 /k B > T 0 ). In this case, the power output is close to its maximum value P ≈ P W , while the efficiency still takes values of up to the order of η ≈ 0.6 η C , in agreement with the bounds discussed in Refs. [14,15,28]. These results clearly show the promising opportunities of step-shaped energy-dependent transmissions, as they can possibly be realized in QPCs, for thermoelectric power production.
Please note that the impressively large values for the efficiency at maximum output power do not, however, violate the Curzon-Ahlbohrn [25] bound, η CA , which relates to the Carnot efficiency as This predicts a bound on the efficiency at maximum power of η CA = 0.5 η C in linear response in ∆T. That this bound is respected, can for example be verified by noting that the efficiency at maximum power of the grey solid line for ∆T/T 0 = 0.01 in panel (c) is only slightly above 0.4η C . Equally, one can check from the green dashed-dotted line in the same panel that the efficiency at maximum power does not exceed the bound for ∆T/T 0 = 5 given by η CA = 0.7 η C . For step energies that are large with respect to the temperature of both reservoirs, T 0 , T 0 + ∆T < E 0 /k B , the power output is reduced, the maximum efficiency, however, increases. In the limit of linear response in the temperature difference, efficiencies close to Carnot efficiency are reached at the expense of close-to-zero power output.

Power Fluctuations and Inverse Fano Factor
During recent years it has become clear that in addition to the power and the efficiency as performance indicators of a heat engine, the fluctuations of the power output, S P , should also be considered [40]. A reliable operation of the heat engine, i.e., where fluctuations are limited, is desirable. This is particularly relevant for nanoscale devices, where fluctuations are always a sizable effect. To analyze the effect of power fluctuations, we note that the relevant fluctuations in this QPC steady-state thermoelectric heat engine are the charge-current fluctuations, since S P = V 2 S I . Therefore, we shift the analysis of power fluctuations to the more straightforward analysis of the Fano factor, see Equation (6).
In Figure 7, we plot the inverse Fano factor 1/F as a function of voltage eV/(k B ∆T) for different barrier smoothness γ, thermal gradients ∆T/T 0 , and step energies E 0 /(k B T 0 ). Please note that we set the inverse Fano factor to zero outside the parameter range where power is produced, to be able to use it as a performance quantifier. This performance quantifier 1/F is desired to be large, meaning that current fluctuations are small with respect to the average. For all parameters, we find that increasing the (negative) voltage decreases the inverse Fano factor 1/F (meaning that the Fano factor F increases). This behavior is attributed to the decrease in charge current as the voltage is moved closer to the stopping voltage V s , while the total noise is less affected. For small voltages, as well as small and negative step energies, increasing the thermal gradient generally increases the inverse Fano factor. These results can be understood from the linear-response expression for the currents and noise, Equations (16)- (18) and below, giving the inverse Fano factor where the absolute value can be omitted when focusing on the voltage window in which power is produced. This expression increases with ∆T and decreases as V goes to more negative values. Increasing ∆T thus increases the current without an accompanied increase in fluctuations because S I = 2k B T 0 G is independent of the bias in the linear response. For large step energies E 0 , the inverse Fano factor no longer increases monotonically in ∆T but a non-monotonic behavior is observed, indicating a more subtle interplay between the fluctuations and the mean value of the current. We note that for almost all parameter values, the inverse Fano factor is substantially smaller than one which can be attributed to the relatively large thermal noise in the present system. In Figure 8, the inverse Fano factor maximized over the voltage, (1/F) V max , is shown for the same parameters as used in Figures 3 and 5. Please note that the maximization only includes the voltage window where positive electrical power is produced. We find that for all three values of smoothness, the maximum inverse Fano factor increases monotonically with increasing ∆T, saturating at values a bit above unity. The Fano factor is thus slightly below unity, a signature of almost uncorrelated, close-to Poissonian, charge transfer (for Poissonian statistics, F = 1). At small ∆T T 0 , close to equilibrium, the noise is large even though the average electrical current is small. As noted above, this is purely due to thermal fluctuations, resulting in a small inverse Fano factor.

Thermodynamic Uncertainty Relation
We now turn to the investigation of the TUR, cf. Equation (14), which provides a combined performance quantifier accounting for power output, efficiency and power fluctuations. We first consider the TUR-coefficient Q TUR in the linear-response regime. Together with Equation (14), we therefore use the relations for power, power fluctuations, entropy and efficiencies, given in Equations (8)- (13). The linear-response expressions for the charge and heat currents occurring in these relations are given in Equation (16) and we furthermore use S I = 2k B T 0 G. With this, we find Maximizing this expression with respect to voltage we find V max → ±∞, resulting in Q TUR = 1/2, and hence, the inequality becoming an equality. However, this voltage is within a voltage regime where power is dissipated (P < 0) and not produced; power production (P ≥ 0) would instead require V s ≤ V ≤ 0. Thus, this is not of practical relevance for the engine performance. Adding the extra condition that P ≥ 0 we instead find V max → 0. The corresponding value of Q TUR on the left-hand side of Equation (32) then becomes L 2 T 0 /(2GK). Expressing this in terms of the figure of merit ZT, given in Equation (19), we can write the bound on the operationally meaningful TUR-coefficient in the linear-response regime as This shows that in the linear response, the parameters of the steady-state thermoelectric heat engine are actually subjected to a tighter bound than given by Equation (32). Please note that this bound is saturated in the limit V = 0, where the power production, the power fluctuations, as well as the efficiency all vanish. Also, only for ideal thermoelectrics, with ZT → ∞, does the bound become 1/2. As seen in Figure 9d, this maximal bound is actually reached for large step energies E 0 . . We note that we set the Q TUR to zero outside the parameter regime where power is produced. The black-dashed line in all panels corresponds to the bound that is given by Equation (33).
The full TUR-coefficient beyond linear response is illustrated in Figures 9 and 10. We find that the inequality Q TUR ≤ 1/2 is always respected, even though this is not guaranteed by scattering theory [45]. Interestingly, we find the tighter bound in Equation (33) to be respected for most parameters, even though the inequality is only proven to hold in the linear-response regime and the bound is expressed in terms of linear-response quantities (given by Equation (17)), only. Violations of the bound given in Equation (33) beyond linear response are observed for sufficiently low E 0 and when the temperature difference is of the order of the magnitude of E 0 (cf. Figure 9a for a sharp barrier). The regimes where a violation can occur are extended when the barrier is smooth (cf. Figure 9e,f). These violations agree with the general notion that dissipation increases when moving away from the linear response [39]. Furthermore, from Figure 9, we find that in the linear response, as well as for small and negative E 0 , Q TUR decreases monotonically as the (negative) voltage is increased. This reflects the behavior of the inverse Fano factor in Figure 7. Importantly, for sharp step energies E 0 , and beyond the linear response, Q TUR is a non-monotonic function of the voltage and takes on its maximum at a point where power production is finite. This non-monotonicity is a consequence of the interplay between the monotonically decreasing inverse Fano factor and the strongly increasing efficiency and power (cf. Figures 2 and 4), as the voltage is changed to more negative values. Figure 10 shows the TUR-coefficient maximized over voltage, Q V TUR,max , as a function of the thermal gradient ∆T and the step energy E 0 . As for the inverse Fano factor, the maximization only includes the voltage window where power is non-negative. For all values of the barrier smoothness, we find that Q V TUR,max generally decreases as a function of ∆T, and a closer inspection reveals small non-monotonic features related to the small violations of Equation (33). This contrasts with the maximized inverse Fano factor, which shows the opposite behavior, cf. Figure 8. The decrease of the fluctuations with ∆T is thus overcompensated by an increase in dissipation which results in the highest values for Q V TUR,max being reached in the linear-response regime. This shows that Q V TUR,max is maximal in regimes, where the η V max is large. Note however that the maximal Q TUR is reached at zero voltage, the maximized efficiency η is reached close to the stopping voltage V S = 0. Furthermore, no features of the line of optimal power production close to P W can be identified in the panels of Figure 10. Smooth barrier (γ = 3k B T 0 ) Figure 10. TUR-coefficient maximized over the bias V, Q V TUR,max , as a function of temperature difference ∆T and step energy E 0 , for three different values of barrier smoothness, γ/k B T 0 = 0, 1, 3 (a-c).

Conclusions and Outlook
In this paper, we have reviewed and extended the analysis of a QPC (or a QPC-like) device, with a transmission probability with a smoothed step-like energy-dependence, as a steady-state thermoelectric heat engine. The interest in a QPC for heat-to-power conversion derives from its optimal performance with respect to the output power, which goes along with rather large efficiencies. We have analyzed the influence of the barrier smoothness on this behavior and found that strong non-linear-response conditions are required to recover a comparable performance.
In addition to the typically studied performance quantifiers-output power and efficiency-we have broadened the analysis by adding the power fluctuations as an independent quantification of performance. The bound on the combination of these three quantities set by the recently identified thermodynamic uncertainty relation, suggests investigating this as a combined performance quantifier.
We have shown that the bound of the thermodynamic uncertainty relation is further restricted if one adds the practical constraint of finite (positive!) output power. In the linear response, we quantify this restriction by the figure of merit ZT. Interestingly, we have found that this combined performance quantifier maximized over the voltage has large values in those parameter regions in which the maximized efficiency is large, while regions of maximal output power are not distinguished. However, while efficiencies take their maximal value in regions close to the stopping voltage in which finite power is produced, accounting for fluctuations shifts the optimal performance value to the limit of zero voltage and zero power production.
Whether this result is unique to the QPC as steady-state heat engine or can be generalized for other thermoelectric devices is a topic of further studies. Our analysis also naturally raises the question of how to quantify the performance of the QPC when operated as a refrigerator [14,15]. Given that QPCs are standard components in many mesoscopic experiments and both the currents and noise are experimentally accessible, we anticipate that our results could be tested in experiments in the near future.
Author Contributions: The project was suggested and conceived by J.S. Numerical maximization of different quantities were mostly performed by S.K. and P.S., analytical estimates were mainly performed by N.D., M.M., and P.S. Thermodynamic interpretation of the TUR aspects was provided by P.P.P. All authors (S.K., N.D., M.M., P.P.P., J.S., P.S.) contributed to the discussion of results and the preparation of the manuscript.