REVIEW 3 major objections 5 minor 63 references
Path integral approach to the calculation of the characteristic function of work
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Path integrals reproduce the two-point-measurement work statistics exactly for the driven oscillator and the expanding piston.
desk verdict A solid derivation paper: the harmonic oscillator work statistics reproduce Ref. [44] and the expanding-piston characteristic function is new, with the main soft spot being unproved interchange of infinite sums and Fourier inversion that Appendix B only partially mitigates. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the path-integral representation of the characteristic function, Eq. (11): $\chi_W(\nu)=\int e^{\frac{i}{\hbar}(S_1[x]-S_2[y])}\rho(x_i,y_i)\delta(x_f-y_f)$, with forward and backward actions $S_1$, $S_2$ that each contain a segment of duration $\hbar\nu$ at fixed work parameter. This split turns the problem into two propagators: for the oscillator the semiclassical propagator is exact because only one classical path contributes, while for the piston the exact propagator is a sum over infinitely many classical paths classified into four collision classes with sign factors from half-wave loss. Gaussian integrations over the intermediate and boundary positions reduce the oscillator case to Eq. (33); the piston case introduces the imaginary error function through the collision-class sums and produces the four-fold sum Eq. (47). The same forward/backward split is then reduced to the classical work functional in the $\hbar\to0$ limit.
What would settle it
Take the right-hand side of Eq. (47) with upper summation limits $N$ on all four indices, Fourier-transform the truncated sum numerically, and compare the resulting weights with the explicit wave-mechanics transition probabilities in Appendix B for modest $\hbar$; if the truncated sums do not converge as $N$ grows, or converge to different weights, the claimed equality fails.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that the characteristic function of work for the two-point measurement scheme admits an exact path-integral evaluation for both prototype systems. For the time-dependent harmonic oscillator, the expression reduces to Eq. (33), which is exactly Eq. (17) of Ref. [44]. For the free particle in an expanding piston, the characteristic function is the four-fold infinite sum in Eq. (47), built from products of imaginary error functions and the four collision classes of classical paths; Fourier inversion gives the delta-weighted work distribution Eq. (50), and the individual coefficients are the transition probabilities Eq. (51). Appendix B proves that these transition probabilities coincide with those obtained from the time-dependent wave-equation solution. The paper also obtains the classical characteristic function of the driven oscillator directly by path integration over classical trajectories, Eq. (60).
Load-bearing premise
The piston result depends on being allowed to take an infinite four-fold sum of terms that grow exponentially with the summation indices, move it inside a Fourier integral, and reorder it term by term; if those sums need regularization or do not converge, Eqs. (50) and (51) are not established.
Editorial extensions
If this is right
- For the driven harmonic oscillator, work statistics can be obtained by evaluating classical-path actions and Gaussian integrals, without expanding in instantaneous eigenstates.
- For the expanding piston, Eqs. (50) and (51) give explicit transition probabilities between instantaneous energy levels of the moving box, an object otherwise available only through numerical wave-packet propagation.
- The classical work distribution of the driven oscillator follows from the same path-integral framework directly, rather than by taking the $\hbar\to0$ limit of the quantum characteristic function.
- The numerical accumulation of Eq. (50) approaches the known classical expanding-piston distribution as $\hbar\to0$, supporting a trajectory-level quantum-to-classical correspondence.
- The paper states that the forward/backward action split is in principle applicable to open quantum systems and quantum fields as well.
Reading between the lines
- One can read Eq. (47) as defining the work distribution by analytic continuation of the sums; a natural step not taken in the paper would be to study the convergence radius or to regularize the sums and check that the result is independent of the regularization.
- The same collision-class enumeration would apply to a box whose wall motion is not uniform, or to higher-dimensional pistons, as long as the classical paths remain piecewise linear; whether the imaginary-error-function integrals would still close is a testable question.
- The paper's equivalence between path-integral and wave-mechanics results suggests that the trajectory work defined by the forward/backward action difference can be assigned a well-defined probability at finite $\hbar$ in these models; proving that for general potentials would require an additional argument beyond what the paper supplies.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives the characteristic function of work for two isolated quantum systems—a harmonic oscillator with a time-dependent frequency and a free particle in a uniformly expanding piston—using Feynman path integrals, and compares the results with Schrödinger-based calculations. It also obtains classical work statistics by taking the classical limit of the path-integral expression and illustrates the quantum-to-classical crossover numerically for the piston. The main claimed achievements are exact analytical work distributions for these two models and a demonstration that the path-integral approach reproduces the two-point-measurement work statistics.
Significance. If the derivations are valid, the paper offers a useful demonstration that Feynman's path integral can serve as a unified tool for quantum and classical work statistics. The harmonic-oscillator characteristic function exactly matches the known result of Deffner and Lutz (Ref. [44]), and the piston transition probabilities are cross-checked against time-dependent Schrödinger solutions in Appendix B. The manuscript contains no fitted parameters and the central quantum results are checked against independent references, which is a genuine strength. The remaining risk is concentrated in the piston derivation, where infinite sums and Fourier transforms are interchanged without a rigorous justification; because of this, the quantitative status of Eq. (50) and the classical-limit demonstration is not yet fully established.
major comments (3)
- [Section III.B, Eqs. (47)-(50)] The step from the characteristic function to the work distribution is a term-by-term Fourier inversion of an infinite four-fold sum. The Erfi-containing coefficients do not decay absolutely (the relevant combinations oscillate and decay only like n^{-1} e^{-i n^2/2} in the large-argument regime), so the distributional interchange of the sum and the Fourier integral is not automatic. Since Eq. (50) is the central piston work distribution and is used in Fig. 1, the authors need to justify the interchange by proving convergence in a suitable test-function space or by introducing a regulator and taking the limit, or they should verify Eq. (50) numerically against a direct Schrödinger evaluation of P(W) for finite truncations. Appendix B checks Eq. (51), but it does not certify the interchange that produces Eqs. (47) and (50).
- [Section III.B, Eqs. (37)-(39) and (39)-(47)] The derivation also integrates the infinite propagator sums term by term over the intermediate position x_b and over the four coordinate variables. No convergence theorem is cited for these operations, and the sums are only conditionally convergent. If any of these interchanges fails, Eq. (47) is not an established characteristic function, even though Eq. (51) is later confirmed by the Schrödinger calculation. Please provide a justification of the interchanges or a numerical check of Eq. (47), for instance by comparing partial sums with the exact TPM characteristic function obtained from the Schrödinger transition probabilities.
- [Section IV, Eq. (54) and Fig. 1] The classical limit is introduced through a stationary-phase argument, but for the expanding piston the authors state that the usual classical work functional is not applicable, and the quantum-to-classical transition is concluded from a single numerical figure with no quantitative convergence test. This is load-bearing for the paper's claim that the path-integral approach works in both quantum and classical thermodynamics; either provide a quantitative analysis of the ℏ→0 limit or explicitly restrict the conclusion to the harmonic-oscillator case.
minor comments (5)
- [Eqs. (29)-(30)] The harmonic-oscillator eigenstate sums start at n=1 while the spectrum E_n^0 = ℏω0(n+1/2) begins at n=0; the normalization of the density matrix appears inconsistent as written. Please correct the index convention or the prefactor.
- [Appendix A] The passages 'After some simplification' and 'After further simplifications' hide the most involved algebra, from Eqs. (A1) through (A8). Since Eq. (33) is a central exact result, please include a more detailed derivation or provide a symbolic-checkable supplementary file.
- [Section III.B, Eq. (39)] The factor (-1)^{5/4} requires a branch specification, since different choices change the overall phase of the propagator.
- [Fig. 1] The caption states 'Color online', but the text does not identify which curve corresponds to which value of ℏ; please add a legend or explicit description. Also define the quantum accumulated work distribution more precisely, including the range of W′ and the meaning of W_min.
- [Throughout] There are minor typos, including 'protocal' for 'protocol' and 'forth coming' for 'forthcoming'.
Circularity Check
No circularity: path-integral derivation is self-contained; the sole self-citation [52] is an exact representation of the TPM characteristic function and is not used to force the target results.
full rationale
The paper's central results are not obtained by fitting a parameter or by defining the output into the input. The characteristic function starts from the standard two-point-measurement expression Eq. (9). The path-integral expression Eq. (11) is attributed to Ref. [52], which shares author H. T. Quan, but the paper immediately states the explicit actions S1 and S2 in Eq. (10), and the rewrite is an exact identity, not an ansatz fitted to the answers. This self-citation supplies a mathematical representation, not a uniqueness argument and not a fitted quantity, so under the review rules it is not load-bearing circularity. For the harmonic oscillator, Eq. (33) is explicitly compared with Eq. (17) of the independent Ref. [44], and the classical analogue Eq. (60) with Eq. (25) of the same Ref. [44]; these are external benchmarks, not renamings of the paper's inputs. For the expanding piston, Eq. (47) is followed by Eq. (50), and Eq. (51) is independently rederived in Appendix B from the time-dependent Schrödinger equation using the exact solutions of Refs. [46,56]. Although Ref. [46] also has an author overlap, Appendix B contains the derivation and the exact solutions originate with the external Ref. [56], so this is a cross-check rather than a circular reduction. The potential lack of justification for interchanging infinite sums and Fourier transforms between Eqs. (47) and (50) is a mathematical-rigor and convergence concern, not a circularity: it does not make the claimed result identical to an input by construction. No fitted parameters, no input-dependent predictions, no uniqueness theorems imported from the authors, and no renamed known results masquerading as derivations were found.
Assumptions & free parameters
assumptions (6)
- domain assumption Semiclassical propagators are exact for the two models.
- domain assumption The action sum in Eq. (35) with four collision classes gives the exact propagator for the expanding piston.
- domain assumption The initial density matrix is thermal and the work is defined by two projective measurements.
- ad hoc to paper The Fourier transform and the infinite sums in Eq. (47) may be interchanged.
- domain assumption In the classical limit the two quantum paths collapse to the classical trajectory satisfying Newton's equation.
- standard math Gaussian and Fresnel integrations are valid for complex parameters.
Cite this review
Pith. "Pith review of Path integral approach to the calculation of the characteristic function of work." pith.science (2026). https://pith.science/paper/P6F4HMWR
@misc{pith2026190809731,
author = {Pith},
title = {Pith review of: Path integral approach to the calculation of the characteristic function of work},
year = {2026},
howpublished = {\url{https://pith.science/paper/P6F4HMWR}},
note = {Machine review of arXiv:1908.09731}
}
read the original abstract
Work statistics characterizes important features of a non-equilibrium thermodynamic process. But the calculation of the work statistics in an arbitrary non-equilibrium process is usually a cumbersome task. In this work, we study the work statistics in quantum systems by employing Feynman's path-integral approach. We derive the analytical work distributions of two prototype quantum systems. The results are proved to be equivalent to the results obtained based on Schr\"{o}dinger's formalism. We also calculate the work distributions in their classical counterparts by employing the path-integral approach. Our study demonstrates the effectiveness of the path-integral approach to the calculation of work statistics in both quantum and classical thermodynamics, and brings important insights to the understanding of the trajectory work in quantum systems.
Figures
Reference graph
Works this paper leans on
- [44]
- [52]
-
[1]
Sekimoto, Stochastic Energetics, Lecture Notes in Physics, V ol
K. Sekimoto, Stochastic Energetics, Lecture Notes in Physics, V ol. 799 (Springer-V erlag, Berlin, 2010)
work page 2010
-
[2]
This result is identical to Eq. (25) in Ref. [44]. We would like to emphasize that we obtain the classical char- acteristic function of work (Eq. (60)) by doing path integra l in the classical trajectory space. But in Ref. [44], this resul t was obtained by taking the classical limit ( ℏ → 0) of the quantum characteristic function of work (33). We believe...
- [3]
-
[4]
and the Crooks fluctuation theorem [7]. By analogy with the partition function, which characterize s completely the properties of a thermal equilibrium state, t he work statistics is an essential function which encodes impo r- tant information about the non-equilibrium thermodynamic process [42, 43]. Nevertheless, the calculation of the work statistics in ...
-
[5]
Seifert, Rep
U. Seifert, Rep. Prog. Phys. 75, 126001 (2012)
2012
-
[6]
Jarzynski, Phys
C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997)
1997
Show all 63 references
-
[7]
Klages, W
Nonequilibrium Statistical Physics of Small Systems: Fluc tua- tion Relations and Beyond, edited by R. Klages, W. Just, and C. Jarzynski (Wiley-VCH, New Y ork, 2013)
2013
-
[8]
Jarzynski, Phys
C. Jarzynski, Phys. Rev. E 56, 5018 (1997)
1997
-
[9]
G. E. Crooks, Phys. Rev. E 60, 2721 (1999)
1999
-
[10]
G. E. Crooks, Phys. Rev. E 61, 2361 (2000)
2000
-
[11]
Hummer and A
G. Hummer and A. Szabo, Proc. Natl. Acad. Sci. U.S.A. 98, 3658 (2001)
2001
-
[12]
Liphardt, S
J. Liphardt, S. Dumont, S. B. Smith, I. Tinoco Jr, and C. B usta- mante, Science 296, 1832 (2002)
2002
-
[13]
Collin, F
D. Collin, F. Ritort, C. Jarzynski, S. Smith, I. Tinoco, and C. Bustamante, Nature (London) 437, 231 (2005)
2005
-
[14]
G. M. Wang, E. M. Sevick, E. Mittag, D. J. Searles, and D. J . Evans, Phys. Rev. Lett. 89, 050601 (2002)
2002
-
[15]
Ciliberto, R
S. Ciliberto, R. Gomez-Solano, and A. Petrosyan, Annu. Rev. Cond. Matt. 4, 235 (2013)
2013
-
[16]
J. P . Pekola and I. M. Khaymovich, Annu. Rev. Cond. Matt. 10, 193 (2019)
2019
-
[17]
Junier, A
I. Junier, A. Mossa, M. Manosas, and F. Ritort, Phys. Rev . Lett. 102, 070602 (2009)
2009
-
[18]
Blickle, T
V . Blickle, T. Speck, L. Helden, U. Seifert, and C. Bechi nger, Phys. Rev. Lett. 96, 070603 (2006)
2006
-
[19]
S. An, J. N. Zhang, M. Um, D. Lv, Y . Lu, J. Zhang, Z. Q. Yin, H. T. Quan, K. Kim, Nat. Phys. 11, 193 (2015)
2015
-
[20]
T. M. Hoang, R. Pan, J. Ahn, J. Bang, H. T. Quan, and T. Li, Phys. Rev. Lett. 120, 080602 (2018)
2018
-
[21]
Talkner, E
P . Talkner, E. Lutz, and P . H¨ anggi, Phys. Rev. E 75, 050102 (2007)
2007
-
[22]
Tasaki, arXiv:cond-mat/0009244 [cond-mat.stat-m ech]
H. Tasaki, arXiv:cond-mat/0009244 [cond-mat.stat-m ech]
-
[23]
Kurchan, arXiv:cond-mat/0007360 [cond-mat.stat- mech]
J. Kurchan, arXiv:cond-mat/0007360 [cond-mat.stat- mech]
-
[24]
Talkner and P
P . Talkner and P . H¨ anggi, Phys. Rev. E93, 022131 (2016)
2016
-
[25]
H. K. Yadalam and U. Harbola, Phys. Rev. A 99, 063802 (2019)
2019
-
[26]
Sampaio, S
R. Sampaio, S. Suomela, and T. Ala-Nissila, Phys. Rev. E 94, 062122 (2016)
2016
-
[27]
Brandner and U
K. Brandner and U. Seifert, Phys. Rev. E 93, 062134 (2016)
2016
-
[28]
Perarnau-Llobet, E
M. Perarnau-Llobet, E. B¨ aumer, K. V . Hovhannisyan, M. Hu- ber, and A. Acin, Phys. Rev. Lett. 118, 070601 (2017)
2017
-
[29]
K. Funo, M. Ueda, T. Sagawa, arXiv:1803.04778v2 [cond- mat.stat-mech]
- [30]
- [31]
-
[32]
Suomela, J
S. Suomela, J. Salmilehto, I. G. Savenko, T. Ala-Nissil a, and M. M¨ ott¨ onen, Phys. Rev. E91, 022126 (2015)
2015
-
[33]
Engel, Europhys
A. Engel, Europhys. Lett. 79, 10003 (2007)
2007
-
[34]
Subas ¸ı and B
Y . Subas ¸ı and B. L. Hu, Phys. Rev. E85, 011112 (2012)
2012
-
[35]
F. W. J. Hekking and J. P . Pekola, Phys. Rev. Lett. 111, 093602 (2013)
2013
-
[36]
Solinas and S
P . Solinas and S. Gasparinetti, Phys. Rev. E 92, 042150 (2015)
2015
-
[37]
B¨ aumer, M
E. B¨ aumer, M. Lostaglio, M. Perarnau-Llobet, and R. Sampaio, arXiv:1805.10096v2 [quant-ph]
-
[38]
Sampaio, S
R. Sampaio, S. Suomela, T. Ala-Nissila, J. Anders, and T . G. Philbin, Phys. Rev. A 97, 012131 (2018)
2018
-
[39]
Guarnieri, N
G. Guarnieri, N. H. Y . Ng, K. Modi, J. Eisert, M. Paternos tro, and J. Goold, Phys. Rev. E 99, 050101 (2019)
2019
- [40]
-
[41]
B. P . V enkatesh, G. Watanabe, and P . Talkner, New J. Phys. 17, 075018 (2015)
2015
-
[42]
A. E. Allahverdyan, Phys. Rev. E 90, 032137 (2014)
2014
-
[43]
H. J. D. Miller and J. Anders, New J. Phys. 19, 062001 (2017)
2017
- [45]
-
[46]
Deffner and E
S. Deffner and E. Lutz, Phys. Rev. E 77, 021128 (2008)
2008
-
[47]
Deffner, O
S. Deffner, O. Abah, and E. Lutz, Chem. Phys. 375, 200 (2010)
2010
-
[48]
H. T. Quan and C. Jarzynski, Phys. Rev. E 85, 031102 (2011)
2011
-
[49]
Z. Gong, S. Deffner, and H. T. Quan, Phys. Rev. E 90, 062121 (2014)
2014
-
[50]
Talkner, P
P . Talkner, P . S. Burada, and P . H¨ anggi, Phys. Rev. E78, 011115 (2008)
2008
-
[51]
Fei and H
Z. Fei and H. T. Quan, to be published
-
[53]
Z. Fei, H. T. Quan, and F. Liu, Phys. Rev. E 98, 012132 (2018)
2018
-
[54]
Funo and H
K. Funo and H. T. Quan, Phys. Rev. Lett. 121, 040602 (2018)
2018
-
[55]
Husimi, Prog
K. Husimi, Prog. Theor. Phys. 9, 381 (1953)
1953
-
[56]
M. G. E. Da Luz and B. K. Cheng, J. Phys. A: Math. Gen. 25, L1043 (1992)
1992
-
[57]
I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, Cambridge, Massachusetts, 2007)
2007
-
[58]
S. W. Doescher and M. H. Rice, Am. J. Phys. 37, 1246 (1969)
1969
-
[59]
Z. Gong, Y . Lan, H. T. Quan, Phys. Rev. Lett. 117, 180603 (2016)
2016
-
[60]
R. C. Lua and A. Y . Grosberg, J. Phys. Chem. B. 109, 6805 (2005)
2005
-
[61]
Bartolotta and S
A. Bartolotta and S. Deffner, Phys. Rev. X 8, 011033 (2018)
2018
-
[62]
Ortega, E
A. Ortega, E. McKay, ´A. M. Alhambra, and E. Mart´ ın- Mart´ ınez, Phys. Rev. Lett.122, 240604 (2019)
2019
-
[63]
J. J. Dong and Y . F. Yang, Phys. Rev. B 100, 035124 (2019)
2019
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.