REVIEW 2 major objections 4 minor
Floquet Driving of Enzymatic Reactions: Counting Statistics and Long-Time Currents
T0 review · 2 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read Stepwise periodic switching of reaction rates in a Michaelis–Menten enzyme produces a net product current: (1/126)(2π/Ω) at high frequency and (1/100)(Ω/2π) at low frequency, with the sign set by the cycling direction of the rate steps.
desk verdict The Floquet current formulas in this paper are correct as far as I can verify; the only real gaps are rigor and reproducibility, not substance. 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 central machinery is the counting-field-dependent Markov generator W(χ;t), whose one-period propagator is decomposed by Floquet theory into an effective Floquet generator W_eff(χ) and a kicked state |Ψ(0;t)⟩; the current is obtained from the χ-derivative of W_eff evaluated in that state. At high frequency, W_eff is constructed by a van Vleck expansion in powers of 1/Ω, with the commutators of Fourier components summed into closed forms involving the Hurwitz zeta function. At low frequency, repeated projection onto the instantaneous dominant eigenmode yields a geometric phase φ(χ), whose χ-derivative gives a current linear in Ω. The stepwise protocol is an ordered sequence γ = {W_0, …, W_
What would settle it
Numerically integrate the two-state master equation for the three-step protocol γ_ac over many periods at Ω = 10π and measure the long-time slope of N_P(t); if it does not converge to (1/126)(2π/Ω) = 1/630 within the stated O(Ω⁻²) error, the high-frequency formula is wrong. At Ω = 0.2π, the low-frequency prediction is (1/100)(Ω/2π) = 1/1000, and both slopes should reverse sign when the step order is reversed to γ_c.
Extended reading notes
Core claim
The paper's central claim is that the long-time current in a periodically driven Markov process is controlled by a counting-field Floquet generator: J∞(Ω) = ⟨S|Ĵ_eff(0)|ψ0⟩, the derivative of the effective generator with respect to the counting field, evaluated in the kicked state. For the three-step protocol with generators (1,1,1,1) → (2,1,1,1) → (1,1,1,2), the authors derive explicit asymptotics: J∞(Ω) = (1/126)(2π/Ω) + O(Ω⁻²) in the high-frequency regime and J∞(Ω) = (1/100)(Ω/2π) in the low-frequency regime, and they show that reversing the cyclic order γ_ac → γ_c flips the sign. In the high-frequency limit the driven reaction is equivalent to a static enzyme with renormalized rates; in
Load-bearing premise
The load-bearing premise is that the van Vleck high-frequency expansion, imported from quantum Floquet theory, remains a valid asymptotic series for the non-Hermitian, abruptly switched Markov generator; the paper offers numerical agreement rather than a proven error bound.
Editorial extensions
If this is right
- Periodic switching alone, with no bias in the time-averaged rates, generates a net product current whose direction is set by the cyclic order of the switching steps.
- In the high-frequency limit, the driven enzyme behaves like a static enzyme with renormalized reaction rates; for the three-step example the effective rate k1 increases and k_-2 decreases by an amount proportional to 1/Ω.
- The long-time current scales as 1/Ω for fast switching and as Ω for slow switching, so the same protocol exhibits opposite frequency dependencies in the two regimes, with a crossover in between.
- Reversing the driving sequence γ_ac to γ_c reverses the sign of the current in both the high- and low-frequency regimes.
- The geometric part of the low-frequency current is independent of how the total period is split among the steps, as long as each step is long enough for relaxation; only the dynamical part depends on the dwell-time weights.
- The current can vanish when averaged rates are balanced, yet appear purely from the ordering of the rate steps, establishing a stochastic-pumping mechanism under strong, abrupt driving.
Reading between the lines
- A direct extension left implicit in the paper: higher derivatives of the cumulant generating function at χ = 0 should yield the full counting statistics — diffusion coefficient, skewness, and beyond — for product counts under periodic driving, so the same Floquet machinery could describe noise and rare-event statistics.
- Because the N-step protocol has a continuous-driving limit as N → ∞, the derived high- and low-frequency formulas should reproduce earlier smooth-driving results in that limit; comparing the coefficients as N grows would test whether the asymptotic structure is universal across driving shapes.
- The effective-rate picture at high Ω suggests a practical control scheme: by choosing Ω, one can tune the renormalized Michaelis–Menten parameters continuously, which could allow optogenetically switched enzymes to be programmed as tunable rate elements in a biochemical circuit.
- The robustness of the geometric current to dwell-time redistribution implies that even irregular gating of an enzyme — where step intervals are not equal — would pump at the same per-period rate as long as each interval reaches quasi-steady state; this may be relevant to in vivo G-protein-driven cAMP signaling, though the paper does not make that claim.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper develops a Floquet formalism for periodically driven continuous-time Markov processes in the presence of a counting field, with the goal of computing long-time currents and counting statistics. The central formal object is Eq. (25), which expresses the net current at finite time through the effective Floquet generator and the kicked state. The paper then derives asymptotic descriptions in two limits: a high-frequency van Vleck expansion of the effective generator (Sec. III.B, Eqs. (33), (62)–(64)) and a low-frequency geometric-phase formulation (Sec. III.C, Eqs. (44)–(46)). These are applied to a two-state Michaelis–Menten enzyme with a three-step discrete driving protocol (Eqs. (51)–(53)). The main quantitative claims are J∞(Ω) = (1/126)(2π/Ω) + O(Ω^{-2}) at high frequency (Eq. (75)) and J∞(Ω) = (1/100)(Ω/2π) at low frequency (Eq. (83)), with sign reversal when the cyclic order of generators is reversed (Eqs. (76), (84)); both are compared with numerical integration in Fig. 5. An effective-reaction-rate interpretation of the high-frequency regime is given in Eqs. (68)–(70).
Significance. If the claims hold, the paper provides a non-perturbative extension of earlier perturbative treatments of stochastic pumping and geometric currents in chemical reaction networks. The high-frequency effective-generator picture is intuitively appealing and directly usable, while the low-frequency geometric decomposition extends known results to discrete, abruptly switched protocols. The paper has no fitted parameters: the asymptotic coefficients are determined analytically and compared with numerical integration. The explicit criterion for the low-frequency regime (Sec. IV.C.3) and the numerical convergence study in App. A are useful additions. The work is likely to be of interest to researchers in stochastic thermodynamics, Floquet engineering, and biophysical modeling of optogenetically controlled enzymes.
major comments (2)
- [§IV.A.2 (Eq. (83)) and App. D] The low-frequency asymptotic coefficient J_dyn = 0 and J_geom = (1/100)(Ω/2π) are stated as 'explicitly calculated', but no calculation is shown. Eqs. (78)–(80) require the dominant left and right eigenvectors of the three 2×2 generators and their χ-derivatives. This is a finite but nontrivial computation, and it is load-bearing for the low-frequency central claim and for the numerical comparison in Fig. 5. I independently checked the coefficient, but the manuscript should make the derivation traceable—either in the main text or an appendix—by giving the eigenvectors and the χ-derivative overlaps that produce the factor 1/100 and the vanishing dynamical contribution.
- [App. C / §III.B] The van Vleck high-frequency expansion, Eqs. (27), (62)–(64), is imported from Hermitian quantum Floquet theory without an error bound or convergence argument for the non-Hermitian, piecewise-constant Markov generator W(χ;t). The validity criterion in §IV.C.3 is a heuristic magnitude estimate (Ω/π ≫ 1) rather than a proof. For finite-dimensional bounded generators, the one-period propagator is analytic near T = 0, so (1/T) log U(T) has a convergent expansion and the gap is likely fillable. I ask the authors to add a short justification (or explicitly state that the expansion is formal and used in the asymptotic sense) so that the central 1/Ω formula is not left as an unproved import.
minor comments (4)
- [§IV.A.1, Eq. (73)] Please state the normalization of |ψ0⟩ used in the high-frequency calculation. Eq. (E15) uses the normalization ⟨S|ψ0⟩ = 1; this is not explicit when the vector (3/7, 4/7)^T is introduced, and it is needed to reproduce Eq. (75) from Eq. (33).
- [Fig. 5 / numerical methods] The numerical simulations are described only qualitatively. Please add a brief numerical-methods paragraph (ODE solver, time-step control, tolerances, number of periods) and, ideally, a data/code availability statement or repository, since the paper claims numerical validation of the two asymptotics.
- [Eq. (64) and App. D] The Fourier-mode notation W_m, with negative indices, is used heavily. A one-sentence definition of the Fourier decomposition and the index convention (as in Eq. (D6)) before Eq. (64) would improve readability and avoid confusion about the sign conventions in the commutators.
- [Eq. (D21)] In the Ω^{-2} coefficient, the expression '8π^4/37' appears arithmetically surprising; it may be a typographical artifact. Please check whether the intended denominator is 729 (or another value) and correct it.
Circularity Check
No significant circularity: the current coefficients are computed from the specified generators, not fitted from the data they are compared with.
full rationale
I walked the claimed derivation chain. The formal current expression Eq. (25) follows from the Floquet decomposition (17) and the parameter-derivative identity (G1), and the high-frequency reduction to Eq. (33) is carried out in App. E using the van Vleck expansion. In the three-step example, every input is fixed: the generator structure (5), the rate tuples (51)-(53), and the dwell times tau=T/3. The high-frequency coefficient 1/126 in Eq. (75) is the explicit matrix element <S|J_eff|psi0> computed from Eqs. (66)-(73); the low-frequency coefficient 1/100 in Eq. (83) is a direct overlap/geometric-phase calculation from the dominant eigenvectors of the three step generators. No parameter is fitted to the numerical long-time current of Fig. 5; the numerics serve only as a consistency check. The choice J_dyn=0 is disclosed in the text as a protocol selection designed to isolate the commutator/geometric term, not as a fitted input that produces the quoted coefficients. The only author-overlapping reference used at a substantive point is [45] for the van Vleck formula, but App. C re-derives the formula from the standard Magnus/log-derivative identity and gives the coefficient-matching steps, so the high-frequency result does not reduce to a self-citation. No authors' uniqueness theorem is invoked to force a choice. The lack of a rigorous non-Hermitian error bound and the omission of the explicit low-frequency eigenvector algebra are accuracy/reproducibility caveats, not circular reductions.
Assumptions & free parameters
free parameters (2)
- rate constants of three-step protocol =
k1=2 in I1; k-2=2 in I2; all others =1
- step count N and dwell time =
N=3, τ=T/3
assumptions (5)
- domain assumption Two-state continuous-time Markov chain with S and P as reservoirs (chemostat)
- domain assumption For χ near 0 the dominant eigenvalue λ0(χ;t) is simple and analytic, and all other eigenvalues have negative real parts
- domain assumption Van Vleck expansion (Eq. (27)) is asymptotic for the non-Hermitian Markov generator at O(Ω^{-1})
- domain assumption Time-scale separation T ≫ δt ≫ τ_rel in the low-frequency regime
- standard math Duhamel/differentiation identities and Drazin pseudoinverse machinery
Cite this review
Pith. "Pith review of Floquet Driving of Enzymatic Reactions: Counting Statistics and Long-Time Currents." pith.science (2026). https://pith.science/paper/LRDDSGTU
@misc{pith2026260717072,
author = {Pith},
title = {Pith review of: Floquet Driving of Enzymatic Reactions: Counting Statistics and Long-Time Currents},
year = {2026},
howpublished = {\url{https://pith.science/paper/LRDDSGTU}},
note = {Machine review of arXiv:2607.17072}
}
read the original abstract
Technologies for artificially controlling chemical reaction systems, such as optogenetics, are rapidly advancing, making it increasingly important to understand reaction dynamics under time-dependent control. When the modulation of reaction rates is periodic in time, the Floquet formalism provides a systematic framework. We develop a Floquet theory for classical stochastic processes that enables the calculation of the current and its counting statistics under such periodic modulation. In particular, we formulate the theory in terms of a counting field and derive general expressions for the first cumulant and the corresponding current. The current is expressed using the effective Floquet generator and the kicked state, and we further obtain general asymptotic expressions for the current in both the high- and low-frequency regimes. As a concrete example to test our analytical expressions, we then apply the results to discrete Floquet driving -- a non-perturbative, stepwise protocol. The setup is motivated by a biochemical system known as cyclic adenosine monophosphate (cAMP) production, which is an enzymatic reaction activated and inhibited by G-proteins. This is formulated as a discretely driven Michaelis--Menten-type reaction model, in which the catalytic activity is switched on and off abruptly in time, and we obtain analytical expressions and numerical results showing how periodic switching of reaction rates generates a long-time product current. In particular, in the high-frequency limit, we show that the effect of the periodic driving can be interpreted through an effective modification of the chemical reaction rates. These results provide a basis for Floquet analysis of periodically driven chemical reactions.
Figures
Figures from the paper (4 more)
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.