{"id":"60d33179-e5d6-4d3b-9ddd-c9a1b8de4705","arxiv_id":"1908.09731","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For a driven harmonic oscillator and an expanding piston, the characteristic function of work is derived with path integrals and proven equivalent to Schrödinger-based results.","lead":"This paper uses Feynman's path integrals to derive analytical formulas for the work statistics of two driven quantum systems, a harmonic oscillator and an expanding piston, and shows the formulas match older Schrödinger-based results. It also applies the same method to the classical versions of both systems, giving a unified route to work fluctuations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Piston result rests on unproved interchange of infinite sums and Fourier transform in Eqs. (47)-(50); a numerical partial-sum check against TDSE would resolve it.","rationale":"The reader's weakest-assumption diagnosis is essentially correct: the load-bearing gap in the piston derivation is the unproved interchange of infinite sums with integration and Fourier transformation. I set agreement to 'partial' because the reader's stated reason that the Erfi terms grow exponentially with n is not quite right; their asymptotic form is oscillatory with 1/n decay, so the issue is conditional, not exponential, convergence. That technical correction does not weaken the concern. The numerical test I propose targets exactly the step from Eq. (47) to Eqs. (50)-(51), and an independent TDSE computation would settle whether the formal series is only conditionally valid or actually incorrect. The harmonic-oscillator part of the paper is supported by an external match to Ref. [44], and Appendix B gives a separate Schroedinger derivation of the transition probability, so I do not see grounds to reject the paper. Keeping the reader's CONDITIONAL verdict is appropriate.","tokens_in":18907,"tokens_out":8869,"duration_ms":99079,"concrete_test":"For the Fig. 1 parameters (m=1, l0=0.01, lf=0.02, u=0.01, beta=1), compute the partial sums S_N(n1,n3)=sum_{n2,n4<=N} A1 A2 A3 A4 in Eq. (51) for a fixed small transition, say n1=n3=1, and compare with |<E^tau_n3|U|E^0_n1>|^2 obtained by numerically integrating the time-dependent Schroedinger equation (e.g., split-step in scaled coordinates). If S_N converges to the TDSE value as N increases, the formal series is validated for that case; if S_N oscillates or converges to a different value, Eqs. (50)-(51) need a regularization argument. Also check whether the series of absolute values diverges; if so, the term-by-term Fourier inversion remains unproved even if a summation method happens to work.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The expanding-piston derivation is carried by four manipulations that are never justified: term-by-term integration of the infinite propagator sums to obtain Eq. (39), the further integrations producing the four-fold sum in Eq. (47), the Fourier inversion from Eq. (47) to the delta-comb in Eq. (50), and the identification of the n2,n4 sub-sums as the transition probabilities in Eq. (51). Every step commutes an infinite sum with an integral or with a distribution. The Erfi factors are not, as one might fear, exponentially growing: for large argument the combinations in A1...A4 oscillate and decay only like n^{-1} e^{-i n^2/2}, so the series are at best conditionally convergent. Conditional convergence is exactly the regime where term-by-term Fourier inversion can fail: the distributional limit need not equal the formal sum of deltas unless uniform convergence or an explicit regularization is established. Appendix B is a genuine independent check of Eq. (51), but it uses the same infinite series for the transition probability, so it does not certify the interchange used to get Eqs. (47) and (50). The harmonic-oscillator result is externally verified against Ref. [44], so the concern is specific to the piston claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":19137,"tokens_out":5859,"duration_ms":63278,"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":[{"comment":"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":"Section III.B, Eqs. (47)-(50)"},{"comment":"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":"Section III.B, Eqs. (37)-(39) and (39)-(47)"},{"comment":"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.","section":"Section IV, Eq. (54) and Fig. 1"}],"minor_comments":[{"comment":"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.","section":"Eqs. (29)-(30)"},{"comment":"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":"Appendix A"},{"comment":"The factor (-1)^{5/4} requires a branch specification, since different choices change the overall phase of the propagator.","section":"Section III.B, Eq. (39)"},{"comment":"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.","section":"Fig. 1"},{"comment":"There are minor typos, including 'protocal' for 'protocol' and 'forth coming' for 'forthcoming'.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely publishable after the convergence issues are addressed. The harmonic-oscillator result is solid and externally verified; the piston result needs either a rigorous justification of the infinite-sum/Fourier interchanges or a numerical convergence check against the Schrödinger solution. I would not reject: the issues are localized and appear fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is the expanding-piston calculation. The harmonic oscillator section re-derives a known result, so the novelty is the method rather than the formula, but the piston characteristic function in Eq. (47) and the transition probabilities in Eq. (51) are new, and Appendix B independently checks them against the Schrödinger-based results of Doescher–Rice and Quan–Jarzynski. That cross-check is the strongest part of the paper: no free parameters, no fitting, and the final probabilities match an independent formalism. The classical harmonic oscillator limit also agrees with Ref. [44], and the numerical quantum-to-classical plot for the piston is a reasonable sanity check.\n\nThe soft spots are real but not fatal. The step from Eq. (47) to the delta-comb Eq. (50) involves Fourier inverting a four-fold infinite sum whose terms are only conditionally convergent; the authors never justify the interchange. Appendix B proves the final probabilities are correct, which actually covers much of the concern, because the TPM work distribution is a delta-comb by definition once the transition probabilities are known. What remains unproven is the path-integral evaluation itself, which also commutes sums with integrals to reach Eq. (47). Since the final result is externally verified, I read this as a rigor gap rather than a correctness failure, but the authors should add a convergence argument or at least a numerical partial-sum test against the TDSE results. Appendix A is also an algebra dump — many steps are just 'after some simplification' — though the final expression checks out against a published result. The classical-limit discussion in Sec. IV is more asserted than proved; for the piston they only have numerics, which is honest but weaker than the quantum part.\n\nThe citation pattern is fine. The paper leans on Ref. [52] for the path-integral representation and on the authors' own earlier work, but the relation to Ref. [44] and the independent check in Appendix B are handled fairly.\n\nWho is this for? People working on quantum work statistics and path-integral methods in stochastic thermodynamics. It is a competent, useful toolbox paper that extends an existing formalism to a nontrivial model. It deserves serious peer review; I would send it out and ask the authors to tighten the convergence issue and expand Appendix A. I'd probably cite the piston result if I worked in this area.","headline":"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.","tokens_in":19678,"tokens_out":1977,"would_cite":true,"duration_ms":25111,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Path integrals reproduce the two-point-measurement work statistics exactly for the driven oscillator and the expanding piston.","keywords":["work statistics","path integral","characteristic function","two-point measurement","quantum thermodynamics","time-dependent harmonic oscillator","expanding piston","quantum-classical correspondence"],"falsifier":"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.","tokens_in":18694,"feed_emoji":"⚛️","tokens_out":11583,"duration_ms":109062,"temperature":0.7,"pith_summary":"The paper aims to show that work statistics in quantum thermodynamics can be computed directly from path integrals, rather than by solving wave equations for the system's states. It derives the characteristic function of work under the two-point measurement scheme for two solvable models: a harmonic oscillator with a time-dependent frequency and a free particle in a rigid box with one uniformly moving wall. For the oscillator, the path-integral result is closed-form and identical to the previously known wave-mechanics result; for the piston, it is an infinite sum over collision classes of classical paths, and its Fourier transform gives the work distribution and transition probabilities, which an appendix proves agree with the wave-mechanics answer. The same path-integral route also gives the classical characteristic function of the driven oscillator analytically, and numerical evaluation of the piston sums shows convergence to the known classical distribution as $\\hbar\\to 0$.","feed_headline":"Path integrals yield exact work statistics for two quantum models","feed_subtitle":"The driven oscillator and expanding piston match two-point-measurement results, with explicit transition probabilities.","key_machinery":"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.","core_discovery":"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).","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the path-integral expression for the characteristic function of work, Eq. (11), from which both model calculations start.","marker":"[52]"},{"why":"Provides the oscillator work-statistics result (Eq. (17)) that Eq. (33) is shown to equal, and the classical counterpart (Eq. (25)) matched by Eq. (60).","marker":"[44]"},{"why":"Gives the exact semiclassical propagator and action formula for the harmonic oscillator used in Eqs. (20)-(21).","marker":"[53]"},{"why":"Supplies the four-class collision enumeration, actions, and sign factors for the free-particle-in-a-box propagator used for the piston.","marker":"[54]"},{"why":"Provides the expanding-piston model and its exact time-dependent wave-mechanics solution, which Appendix B uses to prove Eq. (51).","marker":"[46]"},{"why":"Gives the exact wave functions of the uniformly moving box used in the Appendix B derivation of the transition probabilities.","marker":"[56]"},{"why":"Supplies the classical work distribution of the expanding piston used in Fig. 1 as the $\\hbar\\to0$ comparison.","marker":"[58]"},{"why":"Provides the integral identities, including the imaginary error function, used to evaluate Eqs. (26) and (40)-(49).","marker":"[55]"}],"fun_headline_variants":["Path integrals give exact work statistics for two quantum systems","Quantum work statistics via path integrals: oscillator and piston","Path-integral work distributions match wave-equation results in two models","Exact work characteristic function from path integrals in two systems"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Path integrals give exact work statistics for two quantum systems","Quantum work statistics via path integrals: oscillator and piston","Path-integral work distributions match wave-equation results in two models","Exact work characteristic function from path integrals in two systems"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000604,"raw_usage":{"total_tokens":2760,"prompt_tokens":829,"completion_tokens":1931,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":1864}},"tokens_in":445,"tokens_out":1931,"duration_ms":13005,"temperature":1.0,"reasoning_tokens":1864,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:03:22.198505+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Liu, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the path-integral expression for the characteristic function of work, Eq. (11), from which both model calculations start."},{"cited_title":"D´ ora, A","cited_arxiv_id":null,"evidence_quote":"Provides the oscillator work-statistics result (Eq. (17)) that Eq. (33) is shown to equal, and the classical counterpart (Eq. (25)) matched by Eq. (60)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the exact semiclassical propagator and action formula for the harmonic oscillator used in Eqs. (20)-(21)."},{"cited_title":"Deffner and E","cited_arxiv_id":null,"evidence_quote":"Provides the expanding-piston model and its exact time-dependent wave-mechanics solution, which Appendix B uses to prove Eq. (51)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the exact wave functions of the uniformly moving box used in the Appendix B derivation of the transition probabilities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical work distribution of the expanding piston used in Fig. 1 as the $\\hbar\\to0$ comparison."}],"review_version":1}