{"id":"69391f3f-b1c7-48e5-85b0-f95df0c216c5","arxiv_id":"2509.07272","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"The paper claims that for an underdamped Brownian particle in a spatial temperature gradient, entropy production and extraction rates vanish at zero force while their time-integrated totals stay finite, so a zero rate does not imply equilibrium.","lead":"This paper studies a tiny particle in a fluid whose temperature changes from hot to cold along one direction, and asks whether a zero rate of entropy production means the particle is at equilibrium. The author claims the answer is no: the particle keeps exchanging heat through its kinetic energy, so it stays out of equilibrium even when the production rate is zero.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"At f=0 the reported 'total entropy production' is a spatial difference of -T(x)/2, not a time integral of the vanishing rate ep_dot; the central claim rests on a category error.","rationale":"In good faith, the paper's intended contribution is to show that in a non-isothermal underdamped system without forces, instantaneous entropy production/extraction rates vanish while cumulative measures remain finite, so a zero rate is not a certificate of equilibrium. The most load-bearing requirement for this claim is that the cumulative measures are actually time integrals of the rates. The paper defines them as such, but every finite total it derives is a spatial difference of a state function, not a time integral of a rate. This is a category error that invalidates the central conclusion independent of the controversial small-mass approximation. The reader's verdict identified the same problem in its rationale, but its 'weakest assumption' was the neglect of the spatial derivative term in Eq. (2); that is a serious mathematical issue but it is not the single most direct failure of the central claim. If the time-integral check is run, it will show that the claimed finite totals do not follow from the stated dynamics, so the conclusion 'zero entropy-production rate does not signify equilibrium' is not established. I therefore find no reason to change the reader's REJECT verdict.","tokens_in":18432,"tokens_out":11164,"duration_ms":127656,"concrete_test":"Using the paper's own definition, evaluate the time integral Delta_ep(T) = int_0^T ep_dot(t) dt in the f->0 limit with ep_dot given by Eq. (24) or (27). The integrand is identically zero, so Delta_ep(T) = 0 for all T. Compare this with the claimed total (Th-Tc)/2 from Eq. (36). If the time integral vanishes while the reported total is positive, the finite 'total entropy production' is a spatial boundary term, not an integrated rate, and the central claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires the time-integrated entropy production/extraction to remain finite while the instantaneous rates vanish. The paper defines these totals in Sec. II as time integrals, e.g. Delta_ep(t) = int_0^t ep_dot(t') dt' and Delta_hd(t) = int_0^t hd_dot(t') dt'. However, in the f->0 limit the paper's own expressions for the rates, Eqs. (24) and (27), give ep_dot = hd_dot = 0 identically. Consequently, any time integral of these rates is exactly zero, not (Th-Tc)/2 > 0. The finite quantities reported in Eqs. (31)-(33), (36)-(38), and (48)-(50) are instead boundary differences of the state function Hd(x) = -<mv^2/2 + U> = -T(x)/2, obtained by integrating dHd/dx over x. These are not time integrals and do not represent accumulated entropy production. The abstract explicitly says 'time-integrated values remain finite,' but no such time integral is ever evaluated; Delta_Hd = (Th-Tc)/2 is a spatial state-function difference. This is internally inconsistent with the paper's own entropy balance (Sec. II), which states that a nonequilibrium steady state requires ep_dot = hd_dot > 0. If ep_dot = 0, the system satisfies the paper's criterion for equilibrium, contradicting the headline conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies underdamped Langevin dynamics with position-dependent temperature, deriving entropy production and extraction rates for quadratic, linear, and piecewise-constant temperature profiles, and for a ratchet potential. The central claim is that for a free particle (no external force) in a spatial temperature gradient, the instantaneous entropy production and extraction rates vanish, yet the time-integrated totals remain finite, implying that a zero entropy production rate does not certify equilibrium. The paper also discusses ratchet transport, multiplicative noise, and temperature-dependent friction.","tokens_in":18907,"tokens_out":9433,"duration_ms":119736,"significance":"If the central claim were correct, it would challenge the standard stochastic-thermodynamics identification of vanishing entropy production rate with equilibrium and would have broad implications for non-isothermal systems. The paper also attempts to provide exact analytical formulas for several temperature profiles, which is useful in principle. However, the central result rests on an invalid stationary ansatz and on a category error in which spatial differences of a state function are presented as time integrals of vanishing rates. The ratchet section contains additional unverified ansätze. Thus the main conclusions are not established, and the paper in its current form does not provide a sound basis for the advertised conceptual claim.","major_comments":[{"comment":"The steady-state distribution in Eq. (22) is not a solution of the Fokker-Planck equation Eq. (2) for general γ. Substituting P ∝ exp[-mγ(v-f/γ)^2/(2T)] into Eq. (2) gives a v^2 coefficient proportional to γ-1 in the stationary equation, and the exponent is not dimensionless for a friction coefficient with physical units. The v-integral is 1/√γ, so the claimed normalization to unity also fails unless γ=1. Because Eqs. (24), (27), (41), and (56) all follow from Eq. (22), the central rates are not predictions of the stated model. In addition, with Eq. (22), J'=-f P/m, so Eq. (12) gives ˙ep<0, contradicting the positive values reported in Eq. (24).","section":"Sec. III, Eq. (22)"},{"comment":"The paper defines Δep(t)=∫0^t ˙ep dt and Δhd(t)=∫0^t ˙hd dt in Sec. II. For f=0, Eqs. (24), (27), and (41) give ˙ep=˙hd=0 identically, so any time integral is zero. The finite quantities reported in Eqs. (31)-(33), (36)-(38), and (48)-(50) are spatial differences of the state function Hd(x)=-T(x)/2, not time integrals. The abstract and Sec. VII explicitly claim 'time-integrated values remain finite,' but no such time integral is ever evaluated. This is internally inconsistent with the paper's own definitions.","section":"Sec. II vs Secs. III-IV"},{"comment":"The paper states that a nonequilibrium steady state requires ˙ep=˙hd>0, while equilibrium corresponds to ˙ep=˙hd=0. At f=0 the derived rates vanish, so by the paper's own criterion the system is at equilibrium, not in a nonequilibrium steady state. The proposed alternative indicator, a finite Hd, is not part of the entropy balance and is not shown to be a valid measure of irreversibility; a state-function difference does not by itself imply dissipative dynamics. The conclusion 'zero entropy production rate does not signify equilibrium' is therefore not supported by the framework used.","section":"Sec. II, NESS criterion"},{"comment":"The 'more rigorous' distribution Eq. (39) is presented without verification that it satisfies Eq. (2); it appears to be another local ansatz obtained by neglecting ∂(vP)/∂x. The sentence 'We now integrate the entropy production and extraction rates over x and t' is misleading: Eq. (41) is a rate, and the subsequent totals, e.g., Eq. (49), are again state-function differences, not time integrals. No time integration is actually performed.","section":"Sec. IV, Eq. (39)"}],"minor_comments":[{"comment":"The Fokker-Planck equation is rewritten as ∂P/∂t = k + ∂J'/∂v, but from Eq. (2) the k-term enters with a minus sign. This sign error should be corrected, although it does not affect the later results if k=0 is imposed.","section":"Sec. II, Eq. (9)"},{"comment":"The notation ˙hd(x)=˙ep(x,t) is confusing: the right-hand side contains t but no t-dependence, and the left-hand side is a function of x only.","section":"Sec. III, Eq. (24)"},{"comment":"The caption says 'as a function of the rescaled temperature τ and mass U0'; U0 is the barrier height, not a mass. Please correct the caption and the corresponding axis labels.","section":"Sec. III, Fig. 3 caption"},{"comment":"The multiplicative-noise temperature T(x)=√D |x|^{-z/2} is singular at x=0 and is not dimensionally a temperature. No regularization or domain restriction is given, so integrals over x may diverge.","section":"Sec. VI A, Eq. (67)"},{"comment":"For the piecewise-constant profile Eq. (35), the derivation of ΔHd=(Th-Tc)/2 is not shown. Interface contributions at x=L0/2 need to be addressed explicitly.","section":"Sec. III, piecewise profile"},{"comment":"The text contains numerous typos and OCR-like artifacts (e.g., '⣨' for angle brackets, 'Tome et. at.'), duplicated sentences in Sec. II, and inconsistent notation such as ˙Hd vs Hd and ˙Ep vs ˙ep. A careful copyedit is needed.","section":"Throughout"}],"recommendation":"reject","confidential_remarks":"The central conclusion is not supported by the manuscript's own equations: the 'time-integrated' totals are spatial state-function differences, and the stationary ansatz Eq. (22) does not solve the stated Fokker-Planck equation. These are load-bearing errors that cannot be repaired by local revision. The ratchet section, while possibly interesting, is disconnected from the central claim and rests on similar unverified ansätze. I do not see a path to acceptance within the manuscript's current scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: the central claim is a category error. The finite \"total entropy production\" at f=0 is a spatial difference of the state function Hd(x) = -T(x)/2, not the time integral of the vanishing rates. The paper defines Delta_ep(t) and Delta_hd(t) as time integrals, so if ep_dot = hd_dot = 0 those integrals are exactly zero. The positive numbers in Eqs. (31)-(33), (36)-(38), and (48)-(50) come from integrating dHd/dx over x, which is a different quantity.\n\nThe paper does a few things worth acknowledging. It asks a legitimate question about whether a vanishing rate certifies equilibrium, and it engages seriously with the Tome/Seifert/Sekimoto framework. The ratchet section reproduces standard Brownian motor behavior, and the closed-form expressions for the temperature profiles are derived clearly from the chosen ansatz.\n\nThe soft spots are load-bearing, not cosmetic. First, Eq. (22) is not a solution of the Fokker-Planck equation Eq. (2). The exponent contains an extra factor of gamma, making it dimensionally inconsistent unless gamma=1, and even then it only works because the spatial derivative term partial_x(vP) is dropped. That term is exactly where the temperature gradient enters. Dropping it removes the mechanism the paper credits with irreversibility. Second, the move from H_dot = -d<E>/dt to Hd = -<E> is fine as an energy balance, but it is not a time integral of the entropy production rate. At steady state d<E>/dt = 0, so the global heat dissipation rate is zero. The finite Delta_Hd is just the change in internal energy as the particle moves across the temperature profile; it is not accumulated entropy production. Third, Eq. (12) as written has a minus sign that would make all the rates negative; the results flip the sign without comment. And the paper's own Section II says a nonequilibrium steady state requires ep_dot = hd_dot > 0, which contradicts the claim that ep_dot = 0 with the system still out of equilibrium.\n\nThe conclusion that zero entropy production rate does not imply equilibrium for a free particle in a temperature gradient is likely false for the exact Kramers equation. The true steady state has a nonzero velocity-space current due to the spatial variation of T(x), and the entropy production rate should be positive. The paper's claim is an artifact of the local-Maxwellian approximation.\n\nBottom line: the manuscript is not sound as a contribution to stochastic thermodynamics. The ratchet part may have some value, but the abstract and Sections III-IV rest on an invalid ansatz and a category error. I would not cite it, and I would not send it to peer review in its current form. If the author wants to make the point about cumulative measures, they need to do the honest calculation of the full Kramers equation.","headline":"The central claim is a category error: the finite 'totals' are state-function differences, not time integrals of the vanishing rates, and the steady-state ansatz drops the spatial gradient that is supposed to drive irreversibility.","tokens_in":19287,"tokens_out":8916,"would_cite":false,"duration_ms":111730,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82C31"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a force-free underdamped Brownian particle in a spatial temperature gradient, the entropy production and extraction rates decay to zero at long times, yet their time-integrated totals stay finite—so a vanishing rate does not certify equ","keywords":["underdamped Brownian motion","entropy production rate","entropy extraction","spatial temperature gradient","nonequilibrium steady state","kinetic energy heat transfer","Brownian ratchet","thermophoretic drift"],"falsifier":"Numerically solve the full underdamped Fokker–Planck equation for a linear temperature gradient at zero external force without dropping the v∂P/∂x term, and check whether the integrated entropy production and extraction stay positive; if they vanish, the local-Maxwellian approximation, not the kinetic-energy heat-flow mechanism, is producing the claim.","tokens_in":18359,"feed_emoji":"🌡️","tokens_out":9011,"duration_ms":91894,"temperature":0.7,"pith_summary":"This paper asks whether a vanishing entropy-production rate guarantees thermodynamic equilibrium, and answers no. For an underdamped (inertia-retaining) Brownian particle moving in a non-uniform temperature landscape with no external force, the instantaneous entropy production and entropy extraction rates decay to zero at long times even though heat continues to flow from hot to cold regions. The time-integrated amounts of entropy produced and extracted remain finite and positive, and the paper attributes this persistent irreversibility to heat exchange carried by kinetic energy rather than by external driving. If the claim is right, rate-based diagnostics of equilibrium are incomplete in non-isothermal underdamped systems; cumulative heat transfer is the quantity that reveals the system's nonequilibrium character. The same framework also yields a Brownian ratchet whose velocity depends on particle mass, so thermal gradients alone can produce directed, sortable motion.","feed_headline":"Entropy rate hits zero, but system stays irreversible","feed_subtitle":"In a thermal gradient, a force-free particle's entropy rates vanish while their cumulative totals stay positive.","key_machinery":"The load-bearing object is the steady-state probability distribution P(x,v), approximated as a local Maxwellian (a Gaussian velocity distribution whose temperature varies with position) while the spatial-derivative term v∂P/∂x in the Fokker–Planck equation (the evolution equation for P) is dropped. From this ansatz the paper derives closed-form expressions for the entropy production rate and entropy extraction rate, both equal to f^2/(γT(x)), and for the integrated heat exchange Hd = −⟨(mv^2/2 + U(x))⟩, which stays nonzero in the force-free limit. The mechanism credited with irreversibility is heat exchange via kinetic energy: a particle crossing the temperature landscape carries kinetic ene","core_discovery":"The central claim is that in the zero-force limit of the underdamped Langevin dynamics with a spatially varying temperature, the local entropy production rate and entropy extraction rate vanish—scaling as f^2/(γT(x))—because both are computed from a current that disappears with the external force. But the time- or space-integrated totals of entropy production and extraction remain strictly positive: for linear, quadratic, and piecewise-constant temperature profiles the integrated heat exchange reduces to the same ΔHd = (Th − Tc)/2, with equal positive total entropy production. The paper concludes that entropy conservation on this model does not settle the equilibrium question; the system kee","pith_inferences":["The paper does not test whether the finite integrated totals survive a full solution of the Fokker–Planck equation that retains the v∂P/∂x term; a numerical check of that term would separate the physical conclusion from the local-Maxwellian approximation.","The same rate-versus-integral distinction may apply to overdamped systems with temperature-dependent mobility, where a noise-induced drift can produce nonzero cumulative dissipation with zero applied force.","The predicted mass dependence of ratchet velocity suggests a concrete microfluidic separation experiment: a binary colloid mixture in a periodic thermal landscape should split by mass, with direction and speed controlled by load and barrier height.","An experimental test could measure integrated heat exchange of a Brownian particle cyclically driven through a temperature gradient; the claim predicts positive totals even at instants when the entropy production rate crosses zero."],"forward_implications":["In a force-free underdamped system with a sustained temperature gradient, checking instantaneous entropy production or extraction rates alone would falsely indicate equilibrium; integrated measures are the meaningful diagnostic.","The total heat exchange over one spatial period takes the same value, ΔHd = (Th − Tc)/2, for linear, quadratic, and piecewise-constant temperature profiles, so the result is insensitive to the profile shape.","Adding a periodic ratchet potential to the thermal gradient produces unidirectional motion even at zero load, with velocity dependent on mass, barrier height, and noise intensity—supporting particle sorting along a reaction coordinate.","Spatially varying viscous friction of exponential temperature-dependent form increases entropy production and extraction rates, making inhomogeneous friction an additional driver of nonequilibrium behavior.","The distinction between vanishing rates and nonzero integrated quantities implies that systems can sit in a nonequilibrium steady state while standard rate-based order parameters vanish."],"supporting_citations":[{"why":"Prior underdamped model that this paper extends to quadratic, linear, and piecewise temperature profiles.","marker":"[35]"},{"why":"Author's recent underdamped study whose thermodynamic analysis is extended here.","marker":"[6]"},{"why":"Supplies the trajectory-level methodology for computing ensemble entropy production and extraction rates.","marker":"[7]"},{"why":"Provides the method for deriving entropy production and dissipation rates from the time derivative of Gibbs entropy.","marker":"[8]"},{"why":"Gives the isothermal formulas for entropy production and entropy extraction rates in terms of the probability current.","marker":"[14]"},{"why":"Defines the stochastic-energetics heat dissipation rate used to derive H_d.","marker":"[46]"},{"why":"Companion stochastic-energetics formulation underlying the heat dissipation and free-energy relations.","marker":"[47]"},{"why":"Earlier model results for entropy production and extraction in Brownian ratchets that the ratchet section builds on.","marker":"[19]"},{"why":"Earlier Brownian-motor results for velocity and current reversal used to benchmark the ratchet findings.","marker":"[20]"},{"why":"Prior derivation of entropy extraction rates that informs the rate expressions used here.","marker":"[27]"}],"fun_headline_variants":["Zero entropy rate, yet irreversibility persists","Entropy rate vanishes; irreversibility remains","No force, zero entropy rate, still irreversible","Vanishing entropy rate fails to prove equilibrium","Heat flow leaves lasting irreversibility"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The derivation assumes that the Fokker–Planck term describing how the probability distribution changes along the temperature gradient can be dropped because the particle mass is small and the boundaries are periodic; if that term is not negligible, the predicted vanishing rates and finite integrated totals no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["Zero entropy rate, yet irreversibility persists","Entropy rate vanishes; irreversibility remains","No force, zero entropy rate, still irreversible","Vanishing entropy rate fails to prove equilibrium","Heat flow leaves lasting irreversibility"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000121,"raw_usage":{"total_tokens":880,"prompt_tokens":642,"completion_tokens":238,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":386,"completion_tokens_details":{"reasoning_tokens":170}},"tokens_in":386,"tokens_out":238,"duration_ms":4367,"temperature":1.0,"reasoning_tokens":170,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T22:31:07.946155+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the full underdamped Fokker–Planck equation for a linear temperature gradient at zero external force without dropping the v∂P/∂x term, and check whether the integrated entropy production and extraction stay positive; if they vanish, the local-Maxwellian approximation, not the kinetic-energy heat-flow mechanism, is producing the claim.","supporting_citations":[{"cited_title":"Defaveri, C","cited_arxiv_id":null,"evidence_quote":"Prior underdamped model that this paper extends to quadratic, linear, and piecewise temperature profiles."},{"cited_title":"Zia and B","cited_arxiv_id":null,"evidence_quote":"Author's recent underdamped study whose thermodynamic analysis is extended here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the trajectory-level methodology for computing ensemble entropy production and extraction rates."},{"cited_title":"Seifert, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the method for deriving entropy production and dissipation rates from the time derivative of Gibbs entropy."},{"cited_title":"Harris and G.M","cited_arxiv_id":null,"evidence_quote":"Gives the isothermal formulas for entropy production and entropy extraction rates in terms of the probability current."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the stochastic-energetics heat dissipation rate used to derive H_d."},{"cited_title":"Sekimoto, J","cited_arxiv_id":null,"evidence_quote":"Companion stochastic-energetics formulation underlying the heat dissipation and free-energy relations."},{"cited_title":"Crochik and T","cited_arxiv_id":null,"evidence_quote":"Earlier model results for entropy production and extraction in Brownian ratchets that the ratchet section builds on."},{"cited_title":"Asfaw, Phys","cited_arxiv_id":null,"evidence_quote":"Earlier Brownian-motor results for velocity and current reversal used to benchmark the ratchet findings."},{"cited_title":"Campa, Y","cited_arxiv_id":null,"evidence_quote":"Prior derivation of entropy extraction rates that informs the rate expressions used here."}],"review_version":1}