{"id":"6d400bbc-513a-4991-825b-850eef37a435","arxiv_id":"1908.04316","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":5,"one_line_summary":"Cells with memory matching the environmental correlation scale extract extra drift from spatial fluctuations, exceeding Keller-Segel gradient-sensing predictions.","lead":"This paper shows that E. coli chemotaxis models with a cellular memory can drift faster through rugged chemoattractant landscapes than predicted by instantaneous gradient sensing. The authors derive a formula for this enhanced drift, linking it to the memory time scale relative to the spatial correlation length of the environment.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 25's correction is orders of magnitude too small under reported parameters to explain the 10–30% enhancement shown in Figs. 3–4.","rationale":"In good faith, the paper's qualitative message—that memory can use spatial correlations to enhance drift—may still be true, and the AB simulations may display such an enhancement. My concern is specifically about the quantitative bridge between the numerics and the analytical formula, which is the paper's central theoretical contribution. The reader's weakest assumption concerned the second-order expansion and kernel truncation, and the acknowledged overprediction in one regime. I agree that approximation is imperfect, but the more immediate issue is that Eq. (25), as written with the reported parameters, predicts an effect roughly five orders of magnitude smaller than the enhancement plotted in the figures. That is not a matter of approximation quality; it is a mismatch between the stated parameter values and the reported results. The test I propose is straightforward because the code and data are publicly available, so the authors or a referee can settle it quickly. If the test shows that Eq. (25) does reproduce the curves because the figure parameters were misreported, then the concern is resolved and the paper's verdict could stand. If the test confirms the mismatch, then Eq. (25) needs correction, the figure captions need correction, or the central claim that Eq. (25) explains the enhancement must be weakened. Therefore the appropriate verdict is CONDITIONAL: accept only after this quantitative consistency check is resolved.","tokens_in":122,"tokens_out":10663,"duration_ms":133486,"concrete_test":"Using the public code, evaluate Eq. (25) with the Fig. 3 parameters (βα = 0.05, βση = 10⁻³) on the same γ, μ grid as Fig. 3c and compute Vμ/vKS*. If the maximum deviation from 1 is below 0.1%, the analytic curve cannot be the red line shown in the paper; then check whether the AB circles are instead consistent with a much larger effective ση or with an additional term not present in Eq. (25).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central analytical result is Eq. (25), which gives Δvμ ∝ β²ση². With the parameters stated for Fig. 3 (βα = 0.05, βση = 10⁻³), β²ση² = 10⁻⁶, and the dimensionless prefactor in Eq. (25) is O(10⁻²) at most over the plotted γ, μ range. Thus Δvμ ≈ 10⁻⁸–10⁻⁷, while vKS ≈ 0.006–0.007. The predicted ratio Vμ/vKS* therefore differs from 1 by less than 10⁻⁴. Yet Fig. 3c and Fig. 4a show AB drift velocities 1.05–1.3 times vKS*, with the red analytical curve claimed to describe them. This is an internal consistency problem: the displayed magnitude of the enhancement cannot originate from the second-order correlation term that is the paper's mechanism, unless Eq. (25) is missing a prefactor or the simulations actually used a much larger ση than reported. The reader's identified failure in the short-memory, mild-ruggedness regime is separate; even outside that regime, the stated small-noise values make Δvμ negligible, so the analytic explanation does not quantitatively account for the numerics as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies how cellular memory affects E. coli chemotaxis in one-dimensional rugged attractant landscapes. It introduces an agent-based (AB) model with a bi-lobed response kernel, validates it against the Keller-Segel (KS) model in constant shallow gradients, and reports that AB cells drift faster than the KS prediction when the spatial correlation length of the landscape is comparable to the run length and the cellular memory is of comparable magnitude. To explain this, the authors extend de Gennes' run-time derivation to second order, obtaining an analytic drift velocity Vmu = vKS + Delta v_mu (Eqs. 23-25). They also report that short-memory cells produce long-lived multimodal population distributions in rugged landscapes. Code and data are deposited with DOIs.","tokens_in":22640,"tokens_out":18366,"duration_ms":196998,"significance":"If the central claim held, the paper would be important: it proposes a parameter-free, analytically tractable mechanism by which cellular memory extracts information from spatial correlations, with falsifiable predictions for the dependence of drift enhancement on the memory gamma and correlation length mu. The manuscript has clear strengths: the analytic derivation is explicit in Supplementary Notes 3-5, the AB model is carefully checked against KS in constant gradients, the small-response condition is verified in the simulated regimes, and code and data are publicly available. However, the quantitative inconsistency described in the major comments means that the central claim, as stated, is not supported by the reported equations and parameters.","major_comments":[{"comment":"The central result is internally inconsistent with the stated simulation parameters. For Fig. 3, the caption reports perceived gradient alpha beta = 0.05 and beta sigma_eta = 10^-3. With these values, v*KS = 0.148 alpha beta approximately 7.4 x 10^-3, while the dimensionless prefactor multiplying beta^2 sigma_eta^2 in Eq. (25) is at most about 10^-2 over the plotted gamma and mu ranges. Since beta^2 sigma_eta^2 = 10^-6, the maximum correction satisfies Delta v_mu / v*KS less than or similar to 10^-6; near gamma*KS it is even smaller. Yet Figs. 3c and 4a display enhancements of 5-30% above vKS and show a red Vmu curve tracking those data. Equation (25) as written therefore cannot produce the displayed enhancement. The authors should correct the prefactor or scaling of Eq. (25), state explicitly the parameters used to draw the red curves in Figs. 3-4, and re-verify the AB results against the corrected formula; if the large enhancement persists, an additional mechanism beyond the second-order correlation term must be identified.","section":"Derivation of drift speed, Eq. (25); Figs. 3-4"},{"comment":"The paper's own stated limitation reinforces the quantitative concern. In the regime of short memory and mildly rugged landscapes (0.05 < mu < 0.5, gamma < gamma*KS; Fig. 3c), the authors state that Vmu overpredicts the AB drift and attribute the discrepancy to population heterogeneity not captured by the second-order moment. However, at the parameters of Fig. 3, the bracket in Eq. (25) in this regime is of order 10^-3 or smaller and can change sign; after multiplication by beta^2 sigma_eta^2 = 10^-6, the resulting correction is at most about 10^-9, far too small to be visible on the scale of the figure or to produce the claimed overprediction. This indicates that the second-order correlation term is not the controlling contribution to the observed deviations from KS, and the qualitative attribution to population heterogeneity is not quantitatively supported by the stated magnitudes.","section":"The effect of memory on the drift speed, Fig. 3c"}],"minor_comments":[{"comment":"Please state the values of beta, alpha, sigma_eta, mu, and m used for the curves and error bars in Fig. 4, so that Eq. (25) can be evaluated directly; currently only the Fig. 3 caption gives these parameters.","section":"Fig. 4 caption"},{"comment":"The definitions of tau, t-hat, and the integration order in Eq. (22) are only clear after reading Supplementary Note 5; add a sentence or pointer to Supplementary Figure 6 in the main text to guide the reader.","section":"Eq. (22)"},{"comment":"The notation beta sigma_eta in the figure captions is ambiguous because it is the product beta times sigma_eta; please write it as beta sigma_eta consistently to avoid reading it as a single variable.","section":"Figure captions, notation"},{"comment":"The y-axis of Supplementary Figure 4 is labelled sigma^2_eta, but the quantity plotted is the normalised response variance sigma^2_Lambda; please correct the axis label to match the text.","section":"Supplementary Figure 4"},{"comment":"The statement that 'the KS model accurately predicts the behaviour of the AB population' should be qualified by the parameter range of the validation (here, alpha beta less than or similar to 0.1 and the small-response condition).","section":"Discussion, first paragraph"}],"recommendation":"major_revision","confidential_remarks":"I disagree with the reader's acceptance recommendation. The mismatch between Eq. (25) and the reported enhancement is an order-of-magnitude internal inconsistency that directly affects the abstract's central claim. I recommend major revision rather than rejection, because a corrected prefactor or corrected parameter reporting may resolve the issue; however, if the authors cannot make Eq. (25) quantitatively agree with the AB data after re-derivation and re-checking, the manuscript should not be accepted, since the claimed explanation of the enhanced drift would then be absent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nI read Gosztolai and Barahona on cellular memory in chemotaxis. The idea is appealing: in rugged environments, cells with memory could exploit spatial correlations to drift faster than Keller-Segel gradient-sensing predicts. The paper builds an agent-based model, simulates run-and-tumble cells in a linear gradient plus correlated noise, and derives an analytical correction Δvμ using a second-order expansion of de Gennes' run-time formula. That derivation is explicit, parameter-free, and the code/data are available. The recovery of KS in the limits μ→0, ∞ and γ→0, ∞ is a nice sanity check.\n\nBut there is a serious internal consistency problem. With the stated parameters for Fig. 3, βα = 0.05 and βση = 10^-3, the correction Eq. (25) is Δvμ ≈ (βση)^2 × O(0.01) ≈ 10^-8, while vKS ≈ 6×10^-3. So Vμ/vKS differs from 1 by less than 10^-4. Yet the AB simulations show drift enhancement of 10-30%, and the authors claim the red curve (their Vμ) describes it. That cannot come from the second-order term as written. Either Eq. (25) is missing a prefactor, the simulations used a much larger ση than reported, or the plotted analytical curve is not actually Eq. (25). This is load-bearing: the paper's title claim is that memory enhances navigation and that the enhancement is explained by the correlation term. As written, the explanation does not account for the numerics.\n\nThe reader's earlier note flagged overprediction for short memory in mildly rugged landscapes; that is real but separate. The magnitude mismatch is more basic.\n\nWhat's good: the problem is well-posed, the AB model is standard, the multimodality/bet-hedging observation (Fig. 5) is interesting even if the mechanism is not fully resolved, and the authors' honest discussion of the failed regime is rare and welcome.\n\nMy take: this needs major revision. The authors should either correct the prefactor, use parameters consistent with the observed enhancement, or produce the actual numbers used in the simulations. As it stands, the central analytical result is not validated by the paper's own figures. I would still send it to review, because the question matters and the authors have the machinery to fix it, but I would not cite it yet.","headline":"The analytical correction in Eq. (25) is orders of magnitude too small under the reported parameters to explain the simulated enhancement, so the paper's central claim is not supported as written.","tokens_in":23144,"tokens_out":5883,"would_cite":false,"duration_ms":55979,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92C17"],"pacs":[],"model":"deepseek-v4-flash","headline":"In rugged chemoattractant landscapes, bacteria with memory drift faster than the Keller–Segel gradient-sensing prediction by exploiting spatial correlations.","keywords":["chemotaxis","cellular memory","run-and-tumble","Keller-Segel","rugged landscape","spatial correlations","drift velocity","Escherichia coli"],"falsifier":"Track individual E. coli in a microfluidic channel that superimposes spatial noise of known correlation length µ on a linear attractant gradient; if the measured mean drift does not exceed the Keller–Segel prediction for memories γ comparable to the correlation time, the claimed correlation-exploitation mechanism fails, and the region where the formula overpredicts (0.05<µ<0.5, γ<1/4) is a discriminating test.","tokens_in":22132,"feed_emoji":"🦠","tokens_out":7568,"duration_ms":71119,"temperature":0.7,"pith_summary":"The paper argues that the intrinsic memory of the bacterial chemotaxis pathway is not just a noise filter but a navigation asset: in environments where the attractant has spatial correlations on the scale of a run, cells with memory can drift faster than the classical Keller–Segel (KS) gradient-sensing prediction. The authors establish this with an agent-based run-and-tumble model of E. coli and a new analytical formula for the average drift velocity in rugged landscapes, expressed as the KS velocity plus a positive correction that couples the memory kernel to the landscape's correlation structure. The correction vanishes for constant gradients and white-noise landscapes, and grows when the memory time is comparable to the time scale of perceived fluctuations, with the largest gain near a correlation length of about half a run length. If correct, the result provides a concrete mechanism by which bacteria extract information from environmental structure beyond the local gradient, and it predicts when the standard KS equation must be corrected.","feed_headline":"Memory lets bacteria outrun gradient-sensing in rugged landscapes","feed_subtitle":"E. coli drift faster than Keller-Segel predicts when memory matches the correlation scale of the landscape.","key_machinery":"The central mechanism is the bi-lobed chemotactic memory kernel K(t) = (β/γ)$e^{{-t/γ}}$(t/γ − $t^{2}$/($2γ^{2}$)), which gives positive weight to recent attractant samples and negative weight to older ones, convolved with the perceived signal to set the tumbling rate λ(t) = 1 − Λ(t). The argument extends de Gennes' drift-velocity derivation by keeping the second-order term in the small-response expansion of the run-time average; that term couples the kernel to the spatial autocovariance Cη of the Ornstein–Uhlenbeck noise and produces the positive correction Δvµ. The load-bearing identity is the closed-form integral (Eq. 25) for Δvµ in terms of γ and µ, which turns the abstract idea of using correlations into a quantitative prediction and explains why the optimum sits at memory comparable to perceived fluctuation time.","core_discovery":"The central discovery is that cells with a bi-lobed memory kernel K(t) and constant run speed v0 achieve an average drift velocity ⟨vAB(Sη)⟩ξ ≈ vKS + Δvµ in a landscape Sη(x) = αx + η(x) with additive Ornstein–Uhlenbeck noise of correlation length µ, where vKS = 2βαγ/(1+2γ)^3 is the Keller–Segel drift and Δvµ is a positive correction (Eq. 25) arising from the overlap of the kernel with the spatial autocovariance of the noise. This means the cell uses the correlations it encounters while swimming, not just the instantaneous gradient, and the effect is largest when the memory time γ is commensurate with the perceived fluctuation time. The formula recovers KS in the limits µ→0, µ→∞, γ→0, γ→∞, and α→0, and it matches agent-based simulations across a broad range of memory and correlation length, with optimal memory γ*_µ ≥ 1/4 and maximum speed-up at µ ≈ 1/2. The paper also shows that short-memory cells in mildly rugged landscapes split into long-lived multimodal subpopulations, while long memory keeps the population unimodal Gaussian.","pith_inferences":["If the correction term holds in two and three dimensions, memory tuning could act as a population-level bet-hedging mechanism, with subpopulations of different memory values scanning rugged territories at different rates—a testable prediction for heterogeneous microbial habitats.","The same kernel–covariance coupling could be ported to other bi-lobed sensing systems, such as visual neurons in saccadic search, where the memory is the temporal response of the receptive field and the landscape is the image statistics.","A direct experiment could validate the formula by placing E. coli in microfluidic gradients with engineered spatial noise of known correlation length and comparing measured drift to the KS baseline; the discrepancy should peak near correlation lengths of about 0.1–1 run lengths and vanish for smooth gradients.","The authors' overprediction for short memories hints that higher-order response statistics or population heterogeneity must enter the closure; adding a third-order or variance term might resolve the discrepancy and yield a more accurate effective-drift equation."],"forward_implications":["In constant shallow gradients, agent-based cells match KS drift exactly, so memory gives no advantage when the environment is smooth.","In rugged landscapes with correlation length near one run length, the KS model measurably underestimates population drift, by an amount the formula quantifies.","The optimal memory for navigation is always at least the KS optimum γ*_KS = 1/4, and the largest speed-up over KS occurs when the correlation length is about half a run length.","The enhancement disappears in the white-noise and constant-gradient limits, so the effect is specifically due to spatial correlations rather than noise alone.","Short-memory cells in mildly rugged landscapes lose population coherence, forming long-lived multimodal distributions, while long-memory cells stay unimodal."],"supporting_citations":[{"why":"Basis for the run-and-tumble drift-velocity derivation that this paper extends to second order.","marker":"de Gennes36"},{"why":"Justifies the Keller–Segel limit under the small-response condition |Λ|≪1, which licenses the analytical expansion.","marker":"Xue25"},{"why":"Source of the agent-based model and of the optimal-filtering interpretation of the memory kernel.","marker":"Becker et al.32"},{"why":"Supplies the empirical run-and-tumble parameters and the run-length scale used in non-dimensionalisation.","marker":"Berg and Brown19"},{"why":"Provides the pathway model of tumbling-rate modulation that motivates the bi-lobed response kernel.","marker":"Tu et al.29"},{"why":"Derives the Keller–Segel drift from individual-based run-and-tumble motion, the baseline being corrected.","marker":"Rousset and Samaey41"},{"why":"Gives the diffusion limit and shallow-gradient condition framing the validity of the KS approximation.","marker":"Erban and Othmer42"}],"fun_headline_variants":["Memory lets bacteria outsmart rugged chemical landscapes","Bacteria with memory drift faster in rough chemoattractant fields","Cellular memory gives E. coli a navigation edge on bumpy terrain","Memory-matched bacteria beat pure gradient sensing in rugged mazes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The drift-velocity formula rests on treating the tumbling-rate response as a small perturbation (|Λ|≪1) and ignoring correlations of order higher than two and inputs preceding the last tumble.","fun_headline_variants_meta":{"raw":{"variants":["Memory lets bacteria outsmart rugged chemical landscapes","Bacteria with memory drift faster in rough chemoattractant fields","Cellular memory gives E. coli a navigation edge on bumpy terrain","Memory-matched bacteria beat pure gradient sensing in rugged mazes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1361,"prompt_tokens":939,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":353}},"tokens_in":555,"tokens_out":422,"duration_ms":5432,"temperature":1.0,"reasoning_tokens":353,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:45:27.953297+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track individual E. coli in a microfluidic channel that superimposes spatial noise of known correlation length µ on a linear attractant gradient; if the measured mean drift does not exceed the Keller–Segel prediction for memories γ comparable to the correlation time, the claimed correlation-exploitation mechanism fails, and the region where the formula overpredicts (0.05<µ<0.5, γ<1/4) is a discriminating test.","supporting_citations":[],"review_version":1}