A Lower Bound on Work Extraction Probability Prescribed by Nonequilibrium Work Relation

In nonequilibrium processes, work extraction from a system is subject to random fluctuations associated with the statistical distribution prescribed by its environment. The probability of extracting work above a given arbitrary threshold can be a measure of restriction imposed by experimental circumstances. We present a lower bound for the probability when the work value lies in a finite range. For the case of fixed maximum work, the lower bound gets larger as the free energy difference between initial and final states becomes larger. We point out also that an upper bound previously reported in the literature is a direct consequence of the well-known second mean value theorem for definite integrals.


Introduction
It has been well recognized that the random nature of thermal fluctuation in nonequilibrium processes for microscopic systems entails statistical distribution of work values when we measure them.Thus, we may consider that seeking fundamental constraints or bounds for the probability distribution would represent one of research directions among the vast literature of nonequilibrium science.This article aims at obtaining such bounds.Of course the statement of the 2nd law of thermodynamics bolsters any possible extraction of work; i.e., average work W done on a system invariably exceeds the value of equilibrium free energy difference ∆F between initial and final thermal states (i.e., W ≥ ∆F) with holding equality when the process is reversibly done [1].
We note that in the process of the extraction of work by a mechanically coupled device with a system, its value can be regarded as a random variable with an associated probability distribution.However, the form of the extracted work distribution is not determined from the outset by a first principle behind it.We can only estimate it from measurements.
Up to the present date, a robust constraint on the distribution of the value of extracted work is imposed via the relation exp(−β∆F) = e −βW , (1) which is referred to as the nonequilibrium work relation or the Jarzynski identity [2], in which β is the inverse temperature 1/k B T with the Boltzmann constant k B and the brackets in the right hand side denote an integral of this factor over the observed work distribution P(W), i.e., the statistical mean.Conversely, however, getting information regarding the form of the work distribution itself from this identity is implicit and therefore the precise form must be determined through experimental measurements or numerical simulations based on each model.
The constraint on the work distribution that exceeds a certain value ζ, which is defined by the integral W max ζ P(W)dW is our present primary interest, where we set the maximum work W max .This quantity of course takes on the property of the distribution P(W).A previous study mentions that P(W) is approximately Gaussian according to the central limit theorem when a process proceeds sufficiently slowly [3], and the left tail of it is exponentially suppressed [4].To the contrary, it is reported that in the experiment of single molecule force measurements the distribution of the forward work value shows two distinct peaks, i.e., bimodal distributions [5].The deviation from the Gaussian shape has been reported in several cases.For example, Ref. [6] confirmed that the numerical simulations of a one-dimensional particle model and of soft sphere fluid particles in three dimension show highly non-Gaussian work distribution.In addition, the work value, at which the value of the quantity P(W)e −βW is peaked, is very different from the average work W except for slow switching procedure (thermalization).Moreover, a direct evidence of the profile P(W) by experiments reports the distinct non-Gaussian work distribution [7], where the motion of an overdamped colloidal particle in a time-dependent nonharmonic potential was traced.Also in quantum regime, the observed work distribution in a trapped ion system is found to be different from the Gaussian form [8].
To get useful bounds also for these non-Gaussian cases, we need to invoke mathematical inequalities that any work distribution must satisfy.Indeed, useful physical bounds are often obtained as direct consequences of applying inequalities.For example, the fact that the extracted work from a system on average cannot exceed the decrease in free energy of the system has the strong connection to the Jensen's inequality.That is, the Jensen's inequality ϕ( X ) ≤ ϕ(X) that holds for a convex function ϕ of a variable X readily provides the relation e −β W ≤ e −βW that leads to the above statement, i.e., W ≥ ∆F by combining the identity Equation (1).
With the above situation in mind, the present focus of this paper is how the probability of work extraction that exceeds an arbitrary value ζ denoted as P(W ≥ ζ) can be bounded from below under nonequilibrium processes governed by the equality Equation (1).To achieve this goal, we take advantage of employing the decreasing property of the exponential function.This paper proceeds as follows: In the next section, before seeking the bound, we see an equality that the probability P(W ≥ ζ) must satisfy in terms of the free energy difference and of the end point values of work.In Section 3, a lower bound is presented and some remarks on upper bounds stemming from the random nature of the work variable W are appended.In Section 4, we discuss how the obtained lower bound works for a previous experimental study.The last section summarizes our results and provides a conclusion.

Equality For The Probability P(W ≥ ζ)
Before we derive a lower bound, we first consider an equality that always holds for the probability P(W ≥ ζ).We first restrict the work value in a finite interval, i.e., ζ ∈ [W min , W max ] so that the probability of the work done on a system P(W) takes nonzero value within this range.This is legitimate as a measuring device is mechanically coupled with a system and detects its conceivable finite changes.Recall here that the well-known second mean value theorem for definite integrals (e.g., [9]).Let φ(x) be a bounded and monotonic function on In precise, φ(a) and φ(b) are understood as φ(a) = lim (3) which leads to an equality by rearranging the above This tells us only the existence of the work value ζ somewhere between the interval, and not indicating more than anything that.Indeed, the right hand side does not contain information of where the ζ is located.Recently in the study [10], the r.h.s of Equation ( 4) was derived as an upper bound on the probability P(W ≥ ζ) in a different way from ours and its attainability of this limit is analyzed for two examples: discrete quantum processes for a two-level system and the Szilard-like heat engine composed of a single molecule [11].Note that the sign convention of work W in ref. [10] is opposite from the present consideration; that is, while we set the sign of work W done on a system positive (i.e., the factor is e −βW ), the work done by a system is set to be positive, so that the factor is expressed as e βW .The important indication of this observation is that the upper bound claimed in [10] is nothing but the immediate consequence of the second mean value theorem for definite integrals; that is, the equality always holds.For the unrestricted case of the maximum work (i.e., W max = ∞), we obtain P(W ≥ ζ) = 1 − e −(∆F−W min ) .On the other hand, when the minimum work is unrestricted, the Equation ( 4) does not make sense, indicating that the strategy of separation of the work domain by ζ in the second line of Equation ( 4) is not an effective method in that case.

A Lower Bound On The Probability P(W ≥ ζ)
A way forward to deriving a lower bound on the probability P(W ≥ ζ) can be sought by simply using the decreasing property of the exponential function.With this, we can include the unrestricted case of the minimum work (i.e., W min = −∞) and of the maximum work (i.e., W max = ∞).We first consider the case where ζ lies in the finite interval [W min , W max ].We thus evaluate the exponential factor as The inequality in the second line follows from the decreasing property e −βx ≥ e −βc (x ∈ [−∞, c]) applied for each term in the first line.Therefore, the following lower bound is readily obtained by rearranging terms; Note that contrary to the equality Equation ( 4), the right hand side of the above contains the information of the position of ζ.In addition, the value of the minimum work W min is irrelevant, because it does not appear in the bound.We show in Figure 1 that the lower bound lines are decreasing functions of ζ for a fixed maximum work value.We observe that the bound lines move rightward as ∆F increases and the value of bounds approach unity as the threshold value ζ gets smaller, which is intuitively consistent.The bound value vanishes when ζ equals to the free energy difference.Just to be certain, we stress that the bound Equation ( 6) is relevant only for ζ ≤ ∆F and it is simply replaced by P(W ≥ ζ) ≥ 0 for a range ζ ≥ ∆F.

Some Remarks on the Upper Bounds
As we mentioned in Introduction, the extracted work value is a random variable whose distribution is not a priori known.In this case, probability theory in general gives restrictions for a nonnegative random variable W (e.g., [12]): e βW (exponential Chebyshev's inequality) (7) where Var(W) is the variance of work value, i.e., (W − W ) 2 .The first inequality certainly provides an upper bound on the probability of work value greater or equal to ζ.It is not clear, however, whether it is tight or not due to its universality, that is, it may not be the best possible.For a general purpose, an upper bound would indeed be desired to reflect shape information of the distribution of the work variable W. The direct application of the third one (i.e., the exponential Chebyshev's inequality) gives an upper bound P(W ≥ ζ) ≤ e β(W max −ζ) because the work value does not exceed the maximum, i.e., W ≤ W max .However, this is actually not a useful bound because it always holds with equality when ζ = W max .Furthermore, it gives an exponentially loose bound as the threshold work value ζ departs from the maximum.In sum, if the mean and the variance of the distribution are known, arbitrary distributions can be upper bounded by the Chebyshev inequalities, therefore, how and in what forms of inequalities we should invoke seems an important question to link the free energy difference to work data by real experiments.

Discussion
To obtain the probability P(W ≥ ζ) directly from an experiment of nonequilibrium process is not an easy task.One usually needs the information of the work distribution P(W) and integrates it.Even in this procedure, a reliable bound, if exists, is a helpful guide.We have resorted to the decreasing property of the exponential factor to obtain a lower bound Equation (6).To see how the obtained lower bound works, we examine the work distribution reported in Ref. [7], where the measured distribution of applied work on an overdamped colloidal particle in a time-dependent non-harmonic potential exhibits asymmetry and non-Gaussian.According to Ref. [7], the histogram of the work values is well fitted by the Pearson type III distribution [13], i.e., a gamma distribution with moments: the mean 2.4 k B T and the variance 11.6 k B T. We thus reconstruct the histogram and perform the nonlinear regression fitting to obtain the probability P(W ≥ ζ) by integrating P(W) of the form: The k and θ are referred to as the shape and the rate parameters, respectively.We can deduce the value of W min as −4 k B T from the Figure 4 of Ref. [7], below which the experimentally measured counts negligibly contribute.We determined the parameters as k = 7.32 and θ = 1.16 by the nonlinear least-squares method.The first moment of this distribution, which is given as k/θ − W min , is thus determined as 2.32 [k B T].This value is consistent with the calculated one 2.4 [k B T] in Ref. [7].With these values, the probability P(W ≥ ζ) is thus expressed as where γ(k, x) is the lower incomplete gamma function defined as [13].In Figure 2, we show the curve of the probability P(W ≥ ζ) together with the reproduced work distribution P(W) in the inset.We can read out the value of βW max as 18.0 from the work distribution, which can be regarded as the unrestricted case βW max → ∞, because the value e −βW max almost vanishes in this case.From the reported value e − W k B T 1.03 [7], we can deduce the free energy difference as β∆F = −0.03by applying the Jarzinski equality.Because the free energy difference in the experiment was close to zero, the bound line covers relatively smaller domain of ζ.We think that the bound line extends towards larger ζ if the free energy difference set to larger in the experiment.From this consideration, we think that the present lower bound can be a useful guide over a certain range of the free energy difference.A more refined lower bound to cover a small free energy difference is worth seeking as a future work.
Finally, although the present study considers only the case in which the nonequilibrium work relation Equation ( 1) holds, one should bear in mind the scope of the applicability of this equality.In some single-molecule pulling experiments, where RNA molecules receive time dependent forces, the standard link between work (external force × displacement) and changes in Hamiltonian may not hold [14][15][16].In such cases thermodynamic free-energy changes cannot be estimated accurately and the breakdown of the equality relation can be implied as a corollary.A validity condition for the equality is addressed in Ref. [17], which states that the phase space extension of the system, i.e., the sum over all microscopic states of the system must be kept constant before and after the process.For example, in adiabatic vacuum expansion of ideal gases, the condition cannot be satisfied because no work is done on the gas system during the sudden increase of volume.This leads to vanishing free energy change for any initial microstates in the associated initial phase space extension.Accordingly, such a process is not covered by the Jarzynski's equality and is excluded by the scope of the present bound.It therefore needs to devise a lower bound separately from our approach.

Conclusions
The bound for the probability of observed work above a given arbitrary intermediate value ζ is interesting in its own right in nonequilibrium processes.Since the measured work value is a random variable under a statistical distribution that reflects system's nature, the role of such limitation in experimental settings would be much more restrictive than it has been in the previous studies based on specific models.Although the probability of the work value is amenable to limit by the Chebyshev's inequalities, the presented bound serves as a useful guide.We pointed out also that the previously found upper bound Equation ( 4) is nothing more than a consequence of the application of the mean value theorem for definite integrals in itself and it is thus an equality.The implication of this observation is interesting and important because the probabilistic violation of the second law of thermodynamics can be allowed in general.As long as we limit nonequilibrium processes that proceed under the prescription of the Jarzynski's work relation, our result should also universally holds.

FFigure 1 .
Figure 1.The lower bound lines for the probability P(W ≥ ζ) are provided by Equation (6) as a function of ζ for some fixed values of the free energy difference ∆F.We used a fixed value of W max = 100 and β = 1.0.The bounds vanish when ζ = ∆F.Also with other values of W max and β, the behaviour is similar.

Figure 2 .
Figure 2. The plot data in the inset shows the work distribution constructed from the measured histogram by the experiment of Ref. [7] and the solid line shows the nonlinear least-squares fit to the gamma distribution Equation (8).The solid blue line indicates the probability P(W ≥ ζ) represented by the lower incomplete gamma function Equation (9).The dashed line indicates the lower bound by Equation (6).The values β∆F = −0.03and βW max = 18.0 are used for the calculations of the bound line.