Efficiency statistics at all times: Carnot limit at finite power
Abstract
We derive the statistics of the efficiency under the assumption that thermodynamic fluxes fluctuate with normal law, parametrizing it in terms of time, macroscopic efficiency, and a coupling parameter . It has a peculiar behavior: No moments, one sub and one superCarnot maxima corresponding to reverse operating regimes (engine/pump), the most probable efficiency decreasing in time. The limit where the Carnot bound can be saturated gives rise to two extreme situations, one where the machine works at its macroscopic efficiency, with Carnot limit corresponding to no entropy production, and one where for a transient time scaling like microscopic fluctuations are enhanced in such a way that the most probable efficiency approaches Carnot at finite entropy production.
pacs:
05.70.Ln, 05.70.Fh, 88.05.BcEfficiency quantifies how worth a local gain at the expense of a global loss is. In thermodynamics, “losses” are measured by the rate at which entropy is externalized to the environment in the form of a degraded form of energy, while “gain” is the rate at which entropy is expelled from a system to upgrade its own state. Globally, entropy is produced at rate , and the second law of thermodynamics conveys that locally one cannot earn more of what is globally lost. Then, the efficiency is bounded by the (scaled) Carnot efficiency . Alas, the craving of this limit is deluded by the fact it occurs at zero power, which is useless for any activity to be accomplished in a reasonable time.
This picture is only tenable for macroscopic systems. For microscopic systems subject to random fluctuations, the concept of a stochastic efficiency has been recently introduced by Verley et al. gatien1 (); gatien2 (). The first notion one has to revise is that a fluctuating efficiency can indeed exceed the Carnot limit, when in a machine designed to convert in average a form of input power into a form of output power (e.g. an engine producing work at the expense of a heat flow), for a rare event the input and output are reversed (e.g. a pump that employs mechanical work to absorb heat). Moreover, it has been observed that for timesymmetric protocols in the long time limit the Carnot efficiency becomes the least probable in a “large deviation” sense ld () – a very counterintuitive and fascinating result that, in its timeasymmetric variant gatien2 (); gingrich (), is already subject to experimental inquiry roldan (). Corrections at long finite times have been estimated in Ref. gingrich ().
In this work we derive the full probabilit density function (p.d.f.) of the efficiency, under the assumption that thermodynamic fluxes are distributed with multivariate Gaussian with cumulants growing linearly in time. The efficiency p.d.f. displays quite peculiar features. In particular, it does not afford moments of any order, so that there is no average efficiency and meansquare error. Experimentally, this implies that any data analysis should focus on most probable values. About the latter, after an initial transient the distribution becomes bimodal, as observed numerically in Ref. single (). As time elapses, the more pronounced maximum drifts towards the always smaller macroscopic value of the efficiency, while a less pronounced maximum at higher efficiency moves in the superCarnot region towards infinity. We provide a clear physical interpretation of these two peaks. Finally, we argue that the macroscopic framework fails to capture another way of approaching Carnot at finite entropy production, at finite time, when microscopic fluctuations are enhanced so to affect the macroscopic behavior.
Macroscopic nonequilibrium thermodynamics degroot () is rooted on two assumptions, both of which are today being challenged in the framework of the stochastic theory of nonequilibrium thermodynamics st1 (); st3 (): Certain fluxes , with units of an extensive physical quantity per time, take definite values ; Fluxes are linearly related to their conjugate thermodynamic forces via , where the linear response matrix is assumed to be positivesemidefinite and symmetric by virtue of the Onsager reciprocity relations, yielding a nonnegative macroscopic entropy production rate .
We relax the first assumption, by supposing that at a given time fluxes are distributed with law . Each current produces entropy at rate , for a total entropy production rate , with units of per time. Then, the adimensional efficiency
(1) 
is a stochastic variable distributed with p.d.f.
(2)  
where can be assumed to be positive. A remarkable fact one immediately encounters is that the efficiency can fluctuate beyond the Carnot limit. The probability of an efficiency higher than Carnot coincides with the probability of negative entropy production rate,
(3a)  
(3b) 
where is Heaviside’s step function. The rightmost equations follow from the fluctuation theorem ft (); polespo ()
(4) 
which states that processes producing negative entropy are exponentially disfavored with respect to those producing positive entropy. Therefore, that superCarnot efficiencies are unlikely compared to subCarnot efficiencies is an incarnation of the fluctuation theorem.
Exact results can be obtained by assuming that fluxes are distributed with normal multivariate density function
(5) 
where is the determinant. That (onehalf) the correlation matrix should be identified with the linear response matrix is corroborated by the GreenKubo relations
(6) 
another wellknown consequence of the fluctuation theorem andrieux (). The time dependence in Eq. (5) is due to the fact that the timeintegrated fluxes increase linearly in time, and correspondingly so do their cumulants. Under these assumptions the efficiency p.d.f. Eq. (2) can be exactly calculated supp (). It only depends on four adimensional parameters: The macroscopic efficiency , the coupling parameter that for thermoelectric devices benenti () is related to the socalled figure of merit , the average entropy production rate , which sets the time scale and can be reabsorbed by a time reparametrization , and . Being the only extensive parameter, large stands both for large times and the macroscopic limit. We obtain supp ()
(7) 
where is the error function and
(8a)  
(8b) 
Here, is the determinant of the matrix with dimensionless entries . It can be expressed in terms of our parameters as
(9) 
where accounts for the existence of two probability distributions corresponding to given parameters. For to be real, the known bound
(10) 
must hold benenti (). Importantly, is positive semidefinite.
Let us study the efficiency p.d.f. in detail. First, it is a powerlaw distribution with tails
(11) 
which, after submission of this letter, has been proven to be a universal property of efficiency distributions proesmans (). As a consequence, it does not afford finite moments of any order. Hence, the macroscopic efficiency is not the average efficiency , which is not finite.
In Fig. 1 the efficiency distribution is plotted as the bold curve. Remarkably, for a large class of parameters it displays two maxima at and a minimum, the latter slightly off the Carnot efficiency. Hence, not only superCarnot efficiencies are possible, but indeed there appears a local maximum with efficiency higher than Carnot. To understand its physical origin, we distinguish four operational regimes of the machine, according to the signs of the two contributions and to the entropy production rate. The two regimes contributing to positive efficiencies are the machine that employs process 2 flowing along its spontaneous tendency, to drive process 1 against its spontaneous tendency (e.g. heat engine) and the dual machine where the system’s spontaneous tendency is used to drive the environment against its tendency (e.g. the heat pump). Correspondingly we have where
(12) 
and similarly for . Shaded plots are provided in Fig. 1, showing that each of the two maxima is almost exclusively determined by one of the two modes of the machine, the second of which by inversion of input/output has typical efficiency . Regimes and contribute to the tail of the distribution at .
Let us now study the behavior of in scaled time, depicted in Fig. 2. At we obtain a Cauchy distribution, , with maximum at . We have , and equality can only occur for . This implies that the most probable efficiency decreases in time towards . Furthermore, at thermodynamic equilibrium where all the forces vanish, at finite , it can be shown that , which means that systems at equilibrium do not evolve. As time elapses a transition to a bimodal distribution occurs, with the superCarnot maximum drifting to infinity. We can define a critical time at which there appears an inflection point in . Numerical plots of in terms of and show that the critical time is higher the closer to the maximal efficiency and to the loose coupling condition supp (). Finally, in the long time limit one has sloane () and
(13) 
The largetime behavior is captured by the large deviation rate function , which was first calculated and thoroughly analyzed by Verley et al. gatien1 (); gatien2 (). The rate function has only two extrema, a minimum and a maximum , and asymptotically . Then, the more pronounced maximum tends to the macroscopic efficiency , while the minimum tends to the Carnot efficiency. The second maximum does not appear in the large deviation rate function because at infinite time it moves to infinity, since it belongs to a subdominant decay mode. This proves the existence of a critical time , as there must exist another maximum for the distribution to converge.
The quest for Carnot is very subtle. By Eq. (10) the Carnot bound can be saturated in the limit , giving rise to two extreme situations related to the spectrum and eigenvalues of the response matrix . For (tight coupling), by Eq. (9) the correlation matrix becomes degenerate,
(14) 
where are terms of order . For (singular coupling), tends to the inverse of a degenerate matrix, i.e. , with .
To reach Carnot, a second independent condition (selfduality) must hold: attains value , which affords an interesting interpretation in terms of the probability of the inverse efficiency supp (). When , this condition makes either the null eigenvector of relative to its null eigenvalue, or of relative to its finite eigenvalue. In the tightcoupling regime, this condition is known as the stall force prost ().
Expressing the efficiency in terms of the adimensional parameters and (for ) as entin ()
(15) 
one finds that the two limits towards selfduality and towards tight/singular coupling do not commute,
(16) 
Then, a macroscopic Carnot efficiency is “fragile”, as the selfdual forces needed to attain it are those that slightly out of give a dud machine that dissipates to obtain nothing, with macroscopic efficiency .
Nevertheless, the probabilistic level is richer. At tight coupling the bivariate Gaussian Eq. (5) becomes univariate with support along , and the efficiency p.d.f. a Dirac delta centered at the macroscopic efficiency. Then tightlycoupled machines work macroscopically at all scaled times.
More interesting is the singular coupling. Fig. 3 shows that in this limit all extrema tend to accumulate towards the Carnot efficiency, where the density concentrates. Despite the fact that the two peaks survive, convergence to a Dirac delta can be proven by the following argument argu (): From Eqs. (7,8), , and the efficiency p.d.f. converges to a distribution with support in , which is then necessarily a finite combination of derivatives of the Dirac delta, hormander (). Since for all positive test functions , then necessarily but for . Then, singular coupling pushes the most probable efficiency towards Carnot at fixed ; the shadings in Fig. 3 suggest that in this limit the distribution is fairly insensitive to . Moreover, the contour plot in Fig. 3 supports that the most probable efficiency stays at the same value for probability densities evaluated at a fixed time , showing that convergence to is more and more delayed. However, it must be remembered that the physical time scale is set by the entropy production rate. Necessarily the matrix entries of diverge; then in general also diverges. Still, admits a finite eigenvalue. Picking the forces along the relative eigenvector, , one obtains a finite entropy production rate. Oddly, as discussed above, these conditions are met when the macroscopic machine is dud.
To resume: At singular coupling, the effect of fluctuations is macroscopically visible and permits to work close to Carnot efficiency at finite entropy production rate for sufficiently long physical times. The conditions for which the entropy production rate can be held finite are those under which the machine eventually evolves towards a dud fate. Notice that in this regime the system might flip randomly across the close sharp peaks of the p.d.f.. However, the inset in Fig. 3 suggests that at intermediate times reasonably high typical efficiencies will be favored and that a large separation between such peaks (the dark region of zero probability) occurs. Hence, to put it with a motto, a singular machine doomed to be useless might be efficiently useful for some time due to fluctuations; the better in the short run, the worse in the long. By the GreenKubo relation Eq. (6) the singular coupling limit is approached when correlations between the currents diverge and the inverse correlation matrix becomes degenerate. It is tempting to parallel this behavior to the paradigm of criticality at phase transitions, where fluctuations become macroscopic, correlations diverge and the covariance matrix degenerates zanardi (); pole2 (). In practical applications, the best figure of merit reached is (), in fact quite low. Then, we suggest that this insight might indicate a strategy to look for devices with higher figure of merit.
An important observation here to be made is that singular coupling pushes the system far from equilibrium. The framework of stochastic thermodynamics encompasses such systems by assuming that they are subtended by an underlying Markovian dynamics, giving rise to nonGaussian current statistics. Gaussianity is only recovered in the linear regime at large times by the central limit theorem speck (); engel (). While the model of a Brownian particle in a tilted plane studied in Ref. gatien1 () has exact Gaussian propagators as those studied in this paper, in general Markov processes have a more complex behavior in time, in particular the average flux varies as the system evolves, depending on the initial ensemble. Then, the exact short and largetime behavior of the efficiency distribution might become modeldependent. For asymmetric protocols, a signature of nonGaussian behavior is the offCarnot leastprobable efficiency gatien2 (); gingrich (); roldan ().
Nevertheless, our study points out that in the simplest Gaussian scenario the efficiency p.d.f. manifests peculiar features that might possibly be universal: Powerlaw tails, no finite moments, a naturally occurring transition to a bimodal distribution due to reverse working regimes, etc. Particularly intriguing is the limit of a degenerate or singular covariance matrix. While the former case is intrinsically macroscopic and broadly studied benenti (); entin (), we obtain a clear indication that the singular coupling regime displays an interesting behavior that could lead to the enhancement of the efficiency above its macroscopic value. More light is to be shed on these issues by future inquiry on the finitetime statistics of the efficiency in stochastic models proesmans (); prost () in their rich phenomenology, including maximum power generation espositomax (); schmiedl (), multiterminal machines mazza (), broken timereversal symmetry saito (), the insurgence of phase transitions, and in relation to the issue of efficiency enhancement by noise hanggi () or by decoherence plenio (). Experimental setups that could test these predictions are already available exp1 (); exp2 (); exp3 (); exp4 (); exp5 (). The full statistics of the efficiency close to equilibrium has recently been sampled for a Carnot engine realized with a Brownian particle, in the quasistatic limit where the currents’ statistics is Gaussian roldan (), and data analysis farther away from equilibrium might soon be available.
Aknowledgments.
The research was supported by the National Research Fund Luxembourg in the frame of project FNR/A11/02 and of Postdoc Grant 5856127.
References
 (1) G. Verley, T. Willaert, C. Van den Broeck and M. Esposito, Nat. Commun. 5, 4721 (2014).
 (2) G. Verley, T. Willaert, C. Van den Broeck and M. Esposito, Phys. Rev. E 90, 052145 (2014).
 (3) T. R. Gingrich, G. M. Rotskoff, S. Vaikuntanathan, P. L. Geissler, New J. Phys. 16, 102003 (2014).
 (4) H. Touchette, Phys. Rep. 478, 1, (2009).
 (5) I. A. Martinez, E. Roldan, L. Dinis, D. Petrov, J. M. R. Parrondo, R. Rica, arXiv:1412.1282.
 (6) S. Rana, P. S. Pal, A. Saha and A. M. Jayannavar, Phys. Rev. E 90, 042146 (2014).
 (7) S. R. De Groot and P. Mazur, Nonequilibrium thermodynamics (Courier Dover Publications, New York, 2013).
 (8) U. Seifert, Rep. Progr. Phys. 75, 126001 (2012).
 (9) C. Van den Broeck and M. Esposito, Physica A 418, 6 (2015).
 (10) G. N. Bochkov and Y. E. Kuzovlev, Physica A 106, 443 (1981); ibid. 480 (1981).
 (11) M. Polettini and M. Esposito, J. Stat. Mech. P10033 (2014).
 (12) D. Andrieux and P. Gaspard, J. Chem. Phys. 121, 6167 (2004).
 (13) G. Benenti, K. Saito and G. Casati, Phys. Rev. Lett. 106, 230602 (2011).
 (14) K. Proesmans, B. Cleuren, C. Van Den Broeck, arXiv:1411.3531.
 (15) M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover, New York, 1965).
 (16) O. EntinWohlman, J.H. Jiang and Y. Imry, Phys. Rev. E 89, 012123 (2014).

(17)
Discussion with user Kostya_I on mathoverflow.net, http://mathoverflow.net/questions/178859/
powerlawdistributionwithsupportinx0  (18) L. Hörmander, The analysis of linear partial differential operators I (Springer, Berlin, 1983).
 (19) P. Zanardi, P. Giorda and M. Cozzini, Phys. Rev. Lett. 99, 100603 (2007);
 (20) M. Polettini, Eur. Phys. J. B 87, 215 (2014).
 (21) T. Speck and U. Seifert, Phys. Rev. E 70, 066112 (2004).
 (22) J. Hoppenau, and A. Engel, J. Stat. Mech. 06, P06004 (2013).
 (23) F. Jülicher, A. Ajdar and J. Prost, Rev. Mod. Phys. 69, 1269 (1997).
 (24) M. Esposito, K. Lindenberg and C. Van den Broeck, Phys. Rev. Lett. 102, 130602 (2009).
 (25) T. Schmiedl and U. Seifert, Eur. Phys. Lett. 81, 20003 (2008).
 (26) F. Mazza, R. Bosisio, G. Benenti, V. Giovannetti, R. Fazio and F. Taddei, New J. Phys. 16, 085001 (2014).
 (27) K. Saito, G. Benenti, G. Casati, and T. Prosen, Phys. Rev. B 84, 201306 (2011).
 (28) J. Spiechowicz, P. Han̈ggi, J.Łuczka, Phys. Rev. E 90, 032104 (2014).
 (29) F. Caruso, A. W. Chin, A. Datta, S. F. Huelga and M. B. Plenio, J. Chem. Phys. 131, 105106 (2009).
 (30) D. Collin, F. Ritort, C. Jarzynski, S. B. Smith, I. Tinoco and C. Bustamante, Nature 437 231234 (2005).
 (31) A. Bérut, A. Arakelyan, A. Petrosyan, S. Ciliberto, R. Dillenschneider, and E. Lutz, Nature, 483, 187189 (2012).
 (32) J. V. Koski, T. Sagawa, OP. Saira, Y. Yoon, A. Kutvonen, P. Solinas, M. Möttönen, T. AlaNissila and J. P. Pekola, Nature Phys. Lett. 9, 644 (2013).
 (33) S. Ciliberto, A. Imparato, A. Naert and M. Tanase, Phys. Rev. Lett. 110, 180601 (2013).
 (34) C. Tietz, S. Schuler, T. Speck, U. Seifert and J. Wrachtrup, Phys. Rev. Lett. 97, 050602 (2006).
 (35) See the Supplementary Material for details regarding the derivation of the efficiency PDF, the parametrization, and the various limiting situations.
Appendix A Supplementary material
a.1 Efficiency p.d.f.
In this section we derive the probability density function of the efficiency. Without possibility of confusion, we will denote stochastic variables by the values they take. Let us consider two normally distributed stochastic variables (the fluxes)
(17) 
Defining (the forces), we want to calculate the p.d.f. of the efficiency
(18) 
It is convenient to define (the entropy production rates) with averages . Let . We have
(19) 
yielding
(20) 
Letting be the entries of , then the matrix has entries and it can be expressed as
(21) 
Notice that only depends on three parameters, and or, equivalently, , and . Then the efficiency p.d.f. is given by
(22)  
We have
(23) 
where we introduced
(24a)  
(24b)  
(24c) 
Notice that is adimensional, has dimensions of an entropy rate, and of a squared entropy rate. In fact, simple but tedious calculations show that
(25a)  
(25b)  
(25c) 
Notice that and equality can only hold if , a case that we hereby exclude. We now separate , where
(26) 
Performing a change of variables in , we have
(27) 
We can calculate the first and then change sign to to obtain the second. We obtain
Recognizing the complementary error function and defining
(29) 
we obtain
(30) 
Introducing , the full probability distribution reads
(31) 
where we employed the fact that the error function is odd.
Furthermore, defining
(32a)  
(32b)  
(32c)  
(32d) 
we obtain
(33a)  
(33b)  
(33c)  
(33d) 
a.2 Reparametrization
The task is to express in terms of the parameters and . Notice that . Then
(34a)  
(34b) 
From the first we obtain
(35) 
and letting from the second we get
(36) 
yielding
(37a)  
(37b) 
Given
(38a)  
(38b) 
we get
(39a)  
(39b) 
We then obtain
a.3 Equilibrium case
We consider the case at fixed . We have to evaluate Eq. (2) in the main text given that is distributed with
(41) 
We have
(42) 
where we introduced
(43) 
All follows as for the derivation of the general efficiency p.d.f., but for . Then the p.d.f. reads
(44) 
a.4 Degenerate case (Tight coupling)
Let . Under our conditions , we have
(45) 
The eigenvectors of are: relative to eigenvalue , and relative to eigenvalue . Let
(46) 
be the orthogonal matrix that performs the change of coordinate into the diagonal matrix . Then
(47) 
Furthermore, we introduce the MoorePenrose pseudoinverse of which is obtained by inverting all nonvanishing eigenvalues:
(48) 
Letting , it is known that a degenerate normal distribution is supported along the direction
(49) 
and that it has probability density
(50) 
Moreover, the average currents are not arbitrary but we have , which implies
(51) 
Then, using a test function , we have
(52)  
By Eq. (16) in the main text at , .
a.5 Fluctuation relation for selfdual p.d.f.
The two maxima of are due to converse regimes. It is then natural to define the stochastic variable and look at its p.d.f.
(53) 
Selfduality is the condition upon which , implying the following fluctuation relation for the efficiency
(54) 
Let us show that the choice yields this fluctuation relation. First, it simple to show that implies the symmetry of the normal multivariate for the currents
(55) 
In fact, the Jacobian of the transformation is , and the identity between the quadratic polynomials at exponent can be checked by direct substitution:
(56) 
In fact, equating order by order one can also show that the choice are the only ones yielding Eq. (55).
Now let us consider the efficiency p.d.f.:
where on the third line we performed the change of coordinates and we kept into account the order of the extremes of integration by including a suitable absolute value.
a.6 Efficiency at singular coupling
We say the covariance matrix tends to become singular when its inverse tends to become degenerate. Let be small of order . An almost degenerate inverse takes the form