{"id":"d87dbae7-1979-4f68-8e81-2c37b3c14b1f","arxiv_id":"2411.18103","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For a two-rotor bacterial flagella model with white noise, the phase difference is an Ornstein-Uhlenbeck process and the total heat separates into linear, periodic, and Brownian terms under slow-variable and small-phase-difference approximations.","lead":"This paper derives approximate equations describing how the phase difference and heat dissipation of two hydrodynamically coupled flagella fluctuate in individual noisy runs. The formulas, which combine a linear trend, an oscillation, and Brownian motion, match simulations when synchronization is strong and noise is small.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central formulas survive, but the printed derivation is internally inconsistent: Eq. (23) has the wrong sign, so Eqs. (24)-(26) cannot be obtained from it as written.","rationale":"I verified that substituting the proposed phase solution (25) into the exact heat expression (18) reproduces the deterministic and Brownian parts of (26) at leading order, and that the phase-difference Ornstein-Uhlenbeck solution and its variance are consistent with the stated small-phase-difference condition. The reader's weakest assumption about the sin delta approximation is real and is explicitly acknowledged by the paper's own Fig. 4, so it does not by itself invalidate the conditional claim. The more immediately load-bearing issue is that the printed derivation of the central explicit formulas contains a sign error in Eq. (23) and an incorrect averaging period in Section III: as written, Eq. (24) does not follow from Eq. (23), and averaging over pi/omega does not remove the fast term oscillating at 2omega_phi. These are presentation-and-derivation defects rather than defects in the final formulas, so the paper needs revision but the conditional scientific claim can stand.","tokens_in":8246,"tokens_out":29600,"duration_ms":264468,"concrete_test":"Independently integrate the printed Eq. (23) with delta approximated by 0 and compare the result term-by-term with Eq. (24); if the signs of the Komega_phi t and periodic terms differ, the central derivation is internally inconsistent as printed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The explicit solution for the phase sum is the load-bearing step connecting the SDEs (3)-(4) to the heat formula (26). Combining the exact equations for dphi1 and dphi2 gives dSigma/dt = -Komega sin(Sigma+beta+2omega_phi t) + Komega sin beta cos delta + 2 sqrt(Domega) xi_+. The manuscript's Eq. (23) prints the negative of this right-hand side: +Komega sin(...) - Komega cos delta sin beta. Direct integration of the printed Eq. (23) with delta approximately 0 yields Sigma(t) = -Komega_phi t - (Komega/(2omega_phi))[cos(2omega_phi t+beta)-cos beta] + noise, whereas Eq. (24) has +Komega_phi t and the opposite periodic sign. Since Eqs. (25) and (26) are built on Eq. (24), a reader cannot verify the central claim from the printed derivation without silently correcting the sign. A related inconsistency appears in Section III, where the fast term is said to average out over the period pi/omega; the actual fast frequency is 2omega_phi = 2omega sin beta, so the stated period does not zero the term. Both issues are checkable and fixable, and the final formulas appear correct after the intended corrections, but the derivation as printed is not reproducible.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two hydrodynamically coupled bacterial flagella modeled as active rotors driven by white noise. Using slow-variable and small-phase-difference approximations, it derives approximate stochastic-process descriptions: the phase difference becomes an Ornstein-Uhlenbeck process, and the individual phases and total heat dissipation are expressed as explicit functions of time and Wiener processes. The analytical formulas are compared with Euler-Maruyama simulations for three values of the force-tilt angle beta. The central claim is that single-trajectory heat fluctuations can be characterized without ensemble averaging.","tokens_in":8607,"tokens_out":10550,"duration_ms":88996,"significance":"If the derivation is corrected, the paper provides an analytically tractable stochastic-thermodynamic description of synchronization in a concrete flagellar model, going beyond ensemble-averaged quantities. The formulas are parameter-free in the sense that K, D, and beta are inputs from the model, and the numerical verification for beta = 0.25 pi is a genuine test. The method is model-specific but potentially transferable to other mesoscale synchronizing systems. The main weakness is that several printed equations are internally inconsistent as written, so the central derivation is not reproducible without silent corrections.","major_comments":[{"comment":"Equations (1) and (3) are inconsistent as printed. With the stated definitions F_phi = -F sin(beta) and F_r = F cos(beta), the hydrodynamic term in Eq. (1), after division by zeta b, becomes -K omega cos(phi_2 - beta) sin(phi_1), whereas Eq. (3) contains -K omega cos(phi_2 + beta) sin(phi_1). Since Eq. (3) is the starting point for all subsequent analytic and numerical work, Eq. (1) must be corrected or the sign convention for F_r must be stated consistently.","section":"Section II, Eqs. (1) and (3)"},{"comment":"The fast oscillatory term in Eq. (7) is printed as sin(Phi_i + Phi_j + 2 omega t + beta), and the text states that it averages out over the period pi/omega. From the definition Phi_i = phi_i - omega_phi t, the fast phase is actually 2 omega_phi t, so the averaging period should be pi/omega_phi, not pi/omega (unless sin(beta)=1). As printed, the averaging step is not justified and should be corrected; the final stationary distribution may survive the correction, but the derivation must be internally consistent.","section":"Section III, Eq. (7)"},{"comment":"The stationary solution of Eq. (9) is P proportional to exp[(K cos(beta)/(2D)) cos(Phi_2 - Phi_1)], with normalization 16 pi^2 I0(K cos(beta)/(2D)). The printed argument K/(2D cos(beta)) has the cosine in the wrong place and is inconsistent with the variance 2D/(K cos(beta)) quoted later in the same section. This must be corrected, as it also affects the stochastic entropy expression in Eq. (13).","section":"Section III, Eqs. (10)-(12)"},{"comment":"Equation (23) is printed with the wrong signs relative to the correct reduction of Eq. (22) to the sum variable. Using sigma = Sigma + 2 omega_phi t, the correct equation is dSigma/dt = -K omega sin(Sigma + beta + 2 omega_phi t) + K omega sin(beta) cos(delta) + 2 sqrt(D omega) xi_+, not the printed +K omega sin(...) - K omega cos(delta) sin(beta). Direct integration of the printed equation would produce the opposite signs in Eq. (24), and since Eqs. (25) and (26) are built on Eq. (24), the central formulas cannot be derived from the printed text without a silent sign correction. Please fix the sign and also clarify the K omega_phi notation in Eq. (22).","section":"Section IV, Eqs. (23)-(24)"}],"minor_comments":[{"comment":"The sign convention for delta is inconsistent: the text before Eq. (19) defines delta = phi_1 - phi_2, while the caption of Fig. 3 defines V_delta = <delta^2> = <(phi_2 - phi_1)^2>. Please make the convention uniform, since the sign matters for expressions such as Eq. (19).","section":"Section IV and Fig. 3"},{"comment":"The term printed as 'K omega phi cos delta' is ambiguous. If it means K omega_phi cos(delta), then because omega_phi = omega sin(beta) it equals K omega sin(beta) cos(delta); please write this explicitly to avoid dimensional confusion.","section":"Section IV, Eq. (22)"},{"comment":"Only the beta = 0.25 pi case provides a full quantitative test of Eqs. (20) and (26); for beta = 0.01 pi the slow-variable approximation is not valid, and for beta = 0.49 pi the small-phase-difference approximation breaks down at long times. The paper should state this limitation explicitly rather than implying uniform numerical verification across the three cases.","section":"Section V"},{"comment":"There are several typographical and wording issues: 'under the the slow variable' in the Discussion, 'the coefficient before that is also from the addition of two noise' near Eq. (25), and the term 'dimensionless heat' for Q' = Q/(zeta b^2), which is not dimensionless as written. These should be cleaned up.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The errors identified in the major comments are substantial but appear to be correctable typographical and sign mistakes rather than fundamental flaws. The core formulas seem to be consistent with the simulations for beta = 0.25 pi once the intended corrections are made. I recommend major revision and a careful re-derivation of Eqs. (7), (10), and (23) to ensure the printed derivation is reproducible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is Eq. (26), a closed-form single-trajectory expression for heat dissipation in a two-rotor flagellar model, plus the Ornstein-Uhlenbeck phase-difference process Eq. (20). That goes beyond the ensemble-averaged energetics in Refs. [27,28] and is a legitimate extension. The steady-state distribution is a standard von Mises result and the model itself comes from Ref. [34], so the novelty is concentrated in the trajectory-level formulas.\n\nThe paper does a fair job of checking itself. The derivation is parameter-free: K, D, and beta are set by the earlier model, with no fitting to the simulation. That is a real strength. The Euler-Maruyama checks in the valid regime (beta = 0.25 pi) match the corrected formulas, and the authors honestly show the failure at beta = 0.49 pi, where the small-phase approximation breaks down and the heat variance grows as t^3 rather than linearly. That honesty counts.\n\nBut the derivation as printed is not reproducible. The stress-test note is right: Eq. (23) has the wrong sign. Integrating what is printed gives Sigma(t) with -K omega_phi t and the opposite periodic sign, not Eq. (24). Equations (25) and (26) are built on Eq. (24), so a reader cannot verify the central claim without silently fixing the sign. There is also a period inconsistency in Section III: the fast term is said to average out over the period pi/omega, but the actual fast frequency is 2 omega_phi = 2 omega sin beta, so the stated period does not zero the term. The reader's note about Eq. (10) putting cos beta in the wrong place and Eq. (22)/Eq. (23) notation inconsistencies is also accurate. These look like fixable typos rather than conceptual errors, but as printed they block verification.\n\nThe title overstates scope. This is a specific simplified rotor model with white noise and no experimental data; the validity window is narrow, and the paper itself demonstrates the small-phase approximation failing at long times for beta close to pi/2. The Discussion is appropriately cautious about the 'slow variable and strong coupling approximations,' but the title and abstract suggest broader applicability than the model delivers.\n\nWho is this for? Stochastic thermodynamics practitioners and bead-assay experimentalists who want single-trajectory predictions. They will find Eq. (26) useful, and the paper is written accessibly. The citation pattern is fine; leaning on Ref. [34] is legitimate self-citation because the model comes from there.\n\nI would send this to a serious referee. The core result is checkable and probably correct, and the errors are correctable. The referee should verify the sign corrections and ask the authors to state the approximation windows explicitly. With those fixes, this is a solid extension worth publishing.","headline":"Genuinely new single-trajectory heat formula for a two-rotor flagellar model, but the printed derivation has sign and period errors that must be fixed before the central claim is verifiable.","tokens_in":9048,"tokens_out":2483,"would_cite":true,"duration_ms":21743,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60H10","82C31"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives explicit formulas for the phase difference and heat dissipation along individual trajectories of two synchronizing flagella under white noise.","keywords":["flagellar synchronization","stochastic thermodynamics","phase oscillators","heat dissipation","Langevin equation","Ornstein-Uhlenbeck process","hydrodynamic coupling","single-trajectory energetics"],"falsifier":"Track two bead-rotor flagella (or simulate the full Langevin equations) with parameters such that $K\\cos\\beta$ is comparable to or smaller than $D$, for example $\\beta$ near $\\pi/2$ with weak coupling, and measure the single-trajectory heat variance at long times: if it grows approximately as $t^3$ rather than linearly, the closed-form formula for $Q'(t)$ is not a valid description in that regime.","tokens_in":8057,"feed_emoji":"🦠","tokens_out":5229,"duration_ms":46361,"temperature":0.7,"pith_summary":"This paper asks whether the energetics of synchronizing flagella can be described trajectory by trajectory, without averaging over many noise realizations. For two hydrodynamically coupled rotors driven by white noise, the authors derive explicit stochastic-process formulas for the phase difference and the total heat dissipation under two approximations: the slow-variable approximation and the small-phase-difference approximation. They show that the phase difference is an Ornstein-Uhlenbeck process whose variance approaches $2D/(K\\cos\\beta)$, and that heat dissipation contains a linear drift, a periodic term, and a Brownian term. A sympathetic reader would care because it makes single-trajectory heat fluctuations analytically tractable in a biologically motivated setting, connecting synchronization to stochastic thermodynamics.","feed_headline":"Explicit heat formula for individual noisy flagella trajectories","feed_subtitle":"Closed-form phase and heat for each noise realization, verified by simulation, replace ensemble averages with single-shot predictions.","key_machinery":"The central object is the pair of stochastic phase equations for two rotors, analyzed through the slow variable $\\Phi_i=\\varphi_i-\\omega_\\varphi t$: separating the fast intrinsic rotation leaves a slow, noise-driven phase dynamics. The second load-bearing element is the small-phase-difference linearization $\\sin\\delta\\approx\\delta$, which turns the phase-difference equation into an exactly solvable Ornstein-Uhlenbeck process; the same linearization, inserted into the Stratonovich definition of heat, yields the closed-form formulas for phase and heat. The slow-variable approximation is what allows the rapidly oscillating coupling term to be averaged away.","core_discovery":"Starting from a Langevin description of two bead rotors modeling tethered flagella coupled by Blake-tensor hydrodynamic interactions, the paper claims the following explicit solution. When $2D/(K\\cos\\beta)\\ll 1$ and the slow variable approximation holds, the phase difference $\\delta(t)=\\varphi_1-\\varphi_2$ is a zero-mean Ornstein-Uhlenbeck process with variance $V(t)=\\frac{2D}{K\\cos\\beta}\\left(1-e^{-2K\\omega\\cos\\beta\\, t}\\right)$; the individual phases are $\\varphi_i(t)=\\omega_\\varphi\\left(1+\\frac{K}{2}\\right)t+\\frac{K\\omega}{4\\omega_\\varphi}[\\cos(2\\omega_\\varphi t+\\beta)-\\cos\\beta]+\\sqrt{2D\\omega}\\,W_i(t)$ up to first order in $K$ and $\\sqrt{D}$; and the normalized total heat is $Q'(t)=2\\omega_\\varphi^2(1+K)t-2K\\omega\\sin(\\omega_\\varphi t+\\beta)\\sin(\\omega_\\varphi t)+2\\omega_\\varphi\\sqrt{D\\omega}\\,W_t$. These formulas are validated by numerical simulations when both approximations hold; for $\\beta=0.49\\pi$, where the small-phase-difference assumption fails at long times, the heat variance grows as $t^3$ rather than linearly.","pith_inferences":["Implicit in the results is that the $t^3$ growth of heat variance for $\\beta$ near $\\pi/2$ is a regime marker: when phase diffusion is no longer confined by coupling, heat fluctuations acquire superdiffusive character, so the linear variance formula should be read as conditional on synchronization strength.","One could test the method's scope by measuring, in the same setup, the distribution of single-trajectory heat at fixed time; the theory predicts a Gaussian heat fluctuation driven by the Brownian term, whereas beyond the small-angle regime the distribution should broaden and skew.","The slow-variable scheme suggests a direct extension to $N$ rotors: each pair contributes an Ornstein-Uhlenbeck phase difference, and total heat would be a sum of independent such contributions only in the weak-coupling limit; otherwise pair correlations enter."],"forward_implications":["For parameters satisfying both approximations, single-trajectory heat can be characterized without ensemble averaging; the variance of heat grows linearly in time in the regime where the formulas apply.","The steady-state phase-difference distribution is a von Mises distribution, so the tendency to synchronize is controlled by the ratio $K/(2D\\cos\\beta)$, and no synchronization tendency exists at $\\beta=\\pi/2$.","The explicit phase formula gives a parameter-free prediction for bead-assay trajectories, allowing the work done by a bead on the surrounding fluid to be compared with measured time series.","Because the model is of a generic form, the same approximation scheme can be carried over to other mesoscale synchronizing systems driven by white noise, and extended to additive active noise."],"supporting_citations":[{"why":"Supplies the two-rotor model of tethered flagella as active rotors that pump the surrounding fluid.","marker":"[34]"},{"why":"Provides the Blake tensor Green function with the no-slip boundary condition used to derive the hydrodynamic coupling between rotors.","marker":"[35]"},{"why":"Supplies the stochastic thermodynamics definition of heat as a Stratonovich product used in the dissipation calculation.","marker":"[36]"},{"why":"Gives the stochastic energetics framework that defines heat for overdamped Langevin dynamics.","marker":"[37]"},{"why":"Describes the bead assay with which single flagellar rotation trajectories can be measured experimentally for comparison.","marker":"[31]"}],"fun_headline_variants":["Exact phase and heat formulas for single noisy flagella pairs","Beyond averages: stochastic heat for each flagella trajectory","Closed-form noise: phase and heat for individual flagella","Flagella sync: explicit stochastic phase and heat formulas","Single-run thermodynamics: exact flagella phase and heat"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the phase difference between the two rotors stays small enough that $\\sin\\delta$ can be replaced by $\\delta$; if noise is strong relative to the synchronization tendency, this fails and the explicit phase and heat formulas no longer describe the late-time behavior.","fun_headline_variants_meta":{"raw":{"variants":["Exact phase and heat formulas for single noisy flagella pairs","Beyond averages: stochastic heat for each flagella trajectory","Closed-form noise: phase and heat for individual flagella","Flagella sync: explicit stochastic phase and heat formulas","Single-run thermodynamics: exact flagella phase and heat"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000466,"raw_usage":{"total_tokens":2316,"prompt_tokens":926,"completion_tokens":1390,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":542,"completion_tokens_details":{"reasoning_tokens":1311}},"tokens_in":542,"tokens_out":1390,"duration_ms":9404,"temperature":1.0,"reasoning_tokens":1311,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:32:30.270549+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track two bead-rotor flagella (or simulate the full Langevin equations) with parameters such that $K\\cos\\beta$ is comparable to or smaller than $D$, for example $\\beta$ near $\\pi/2$ with weak coupling, and measure the single-trajectory heat variance at long times: if it grows approximately as $t^3$ rather than linearly, the closed-form formula for $Q'(t)$ is not a valid description in that regime.","supporting_citations":[{"cited_title":"Silverman and M","cited_arxiv_id":null,"evidence_quote":"Supplies the two-rotor model of tethered flagella as active rotors that pump the surrounding fluid."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Blake tensor Green function with the no-slip boundary condition used to derive the hydrodynamic coupling between rotors."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the stochastic thermodynamics definition of heat as a Stratonovich product used in the dissipation calculation."},{"cited_title":"Nakamura, Y","cited_arxiv_id":null,"evidence_quote":"Gives the stochastic energetics framework that defines heat for overdamped Langevin dynamics."},{"cited_title":"Izumida, H","cited_arxiv_id":null,"evidence_quote":"Describes the bead assay with which single flagellar rotation trajectories can be measured experimentally for comparison."}],"review_version":1}