{"id":"23c27f63-12cc-41a3-8a62-745811f31178","arxiv_id":"2607.05239","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Longer bacterial chains in nematic liquid crystals swim faster because relative-speed sorting preferentially packs faster cells into longer chains, not because force balance itself depends on length.","lead":"Motile E. coli in nematic liquid crystals form long-lived swimming chains whose average speed rises with chain length. The effect is explained as swimming-limited aggregation: faster cells catch slower ones first and are preferentially sorted into longer chains.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"The length–speed correlation is partly a one-parameter fit of the dry model at a matched dynamical age; the pre-aggregation speed distribution is not independently measured.","rationale":"The reader correctly isolates the dry mechanical coupling model and the post-hoc matching of dynamical age as the weakest links. The qualitative sorting picture is robust: the minimal three-body model (Eqs. 1–2, Fig. 2) already shows that faster cells meet neighbours sooner, and the agent-based runs produce a positive length–speed correlation whenever speed variance is appreciable. The load-bearing quantitative claim, however, rests on a one-parameter fit performed after aggregation has begun, using a residual rather than a true initial speed distribution. That is precisely the concern the reader flags. No stronger internal inconsistency appears; hydrodynamics, elasticity and reversals are acknowledged omissions that would mainly shift the effective age or slope, not reverse the sign of the correlation. Hence the verdict remains CONDITIONAL, with the same confidence level, pending an independent pre-aggregation speed measurement.","tokens_in":23994,"tokens_out":606,"duration_ms":5233,"concrete_test":"Measure the speed distribution of solitary bacteria in the same DSCG chambers at the earliest accessible time after mixing (or in a dilute control that suppresses aggregation) and re-run the agent-based simulations with that fixed distribution (no free σ_V). If the predicted V_n(n) at the experimental ⟨n⟩ no longer matches Fig. 1e within error, the sorting explanation of the length–speed correlation is under-constrained.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that longer chains swim faster because of relative-velocity self-sorting, not because chain speed itself depends on length. Under the dry model (Appendix E, Eqs. 3–5) chain speed is exactly the arithmetic mean of constituents, so any positive length–speed correlation must come from non-random composition. The simulations reproduce Fig. 1e only after (i) estimating V from the already-aggregated experimental dataset and (ii) fitting the single free parameter σ_V so that the simulated speed–length curve at the time t* when ⟨n⟩ ≈ 2.2 matches experiment. The paper itself notes (§IV) that the solitary-bacterium speed histogram in Fig. 1f is the residual population after faster cells have already joined chains, and that experimental t = 0 is ill-defined. Thus the quantitative match is not an independent prediction of the sorting mechanism; it is a consistency check under the dry-model assumptions. If the true pre-aggregation speed distribution were narrower, or if hydrodynamics/elasticity systematically raised chain speeds with n, the same curve could be obtained without the claimed self-sorting.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript reports that motile E. coli in nematic liquid crystals form long-lived 1D chains whose average swimming speed increases with chain length n. A single merging event is consistent with a dry force-balance model in which chain speed equals the arithmetic mean of precursor speeds (Fig. 1d; Appendix E, Eqs. 3–5). The authors attribute the positive length–speed correlation to a dynamic self-sorting process driven by cell-to-cell speed variability: faster bacteria encounter neighbours sooner and are preferentially incorporated into longer chains. This is supported by a minimal three-body encounter model (Eqs. 1–2, Fig. 2) and by agent-based simulations of irreversible 1D aggregation that, after estimating mean speed V from the experimental weighted average and fitting one free parameter σ_V at the simulation time t* where ⟨n⟩ matches experiment (≈2.2), reproduce both the mean speed–length curve and the shape of the speed distributions (Fig. 1e–f). The authors term the process “swimming-limited aggregation” and contrast its time scale with passive LC-mediated aggregation.","tokens_in":24321,"tokens_out":1249,"duration_ms":9307,"significance":"If the self-sorting interpretation holds, the work identifies a clean, experimentally accessible regime of active aggregation whose kinetics and emergent length–speed correlation are controlled by the intrinsic speed distribution of the agents rather than by reaction rates or elastic attractions. The combination of large trajectory statistics, a transparent dry kinematic rule, a minimal analytic encounter model, and publicly released simulation code is a genuine strength and makes the mechanism falsifiable in related confined active systems. The result is of clear interest to soft-matter and active-matter communities studying motility-induced organization in anisotropic media.","major_comments":[{"comment":"§IV and Fig. 1e–f: The quantitative match is obtained by estimating V from the already-aggregated experimental dataset (conserved under the dry model) and fitting the single free parameter σ_V so that the simulated speed–length curve at the dynamical age t* (defined by matching experimental ⟨n⟩≈2.2) reproduces experiment. The paper itself notes that the solitary-bacterium histogram in Fig. 1f is residual after faster cells have joined chains and that experimental t=0 is ill-defined. The agreement is therefore a consistency check under the dry-model assumptions rather than an independent prediction of the sorting mechanism. An independent measurement (or at least a robust bound) of the pre-aggregation speed distribution, or an explicit demonstration that the qualitative positive correlation survives without fitting σ_V to the same curve, is needed to secure the central claim.","section":null},{"comment":"Appendix E, Eqs. (3)–(5) and §III.B: The dry mechanical coupling model forces chain speed to equal the arithmetic mean of constituents by assuming identical drag coefficients, instantaneous flagellar rearrangement after head-on collisions, and complete neglect of hydrodynamics and LC elasticity once a chain has formed. These assumptions are used both to interpret the single merging event (Fig. 1d) and to drive the simulations that are then fitted to experiment. The manuscript acknowledges that hydrodynamics would likely increase the slope of the speed–length curve (and therefore shift the matching time t*). Because any systematic n-dependent enhancement of chain speed would produce a positive correlation even without self-sorting, the paper should either (i) provide a quantitative estimate of the hydrodynamic/elastic correction or (ii) show that the observed correlation cannot be explain","section":null}],"minor_comments":[{"comment":"Fig. 1e inset and §II: The observation count drops rapidly with n; error bars on the longest chains are large. A brief statement of the minimum sample size used for the mean-speed points would help the reader assess robustness.","section":null},{"comment":"§III.A, Eq. (1): The notation for the cumulative distribution F_V and the absolute-value arguments is dense; a short sentence defining the two bracketed factors would improve readability.","section":null},{"comment":"Fig. 3(d) inset: The non-monotonic speed–length relation reported for σ_V/V=1/3 is interesting but only briefly discussed; a sentence on whether this feature is expected to appear in experiment would be useful.","section":null},{"comment":"Appendix B: The passive-aggregation estimate relies on an anchoring energy taken from B. subtilis; a short caveat that W has not been measured for E. coli would be appropriate.","section":null},{"comment":"Code availability is excellent (GitHub links for both image analysis and simulations). Adding a short README note on the exact random-seed protocol used for the 50-run averages would further aid reproducibility.","section":null}],"recommendation":"major_revision","confidential_remarks":"The central physical idea (relative-velocity self-sorting under a dry averaging rule) is attractive and the experimental statistics are solid. The main risk is that the quantitative agreement is partly circular (V taken from the aggregated data, σ_V fitted at matched ⟨n⟩). I would accept after the authors either supply an independent pre-aggregation speed distribution or clearly reframe the theory as a consistency check and strengthen the argument that hydrodynamics/elasticity alone cannot produce the observed correlation. Scope is appropriate for a soft-matter journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new piece is the positive length–speed correlation for E. coli chains in nematic DSCG, plus a transparent mechanism: faster cells meet neighbours sooner and get locked into longer chains. Mushenheim already saw reversible chains; this paper quantifies the speed trend and shows why it appears even though a newly merged chain must swim at the arithmetic mean of its parts.\n\nWhat works. Thousands of trajectories give a clean mean-speed-vs-n curve and right-skewed speed histograms. The single merging event in Fig. 1d matches the dry force-balance prediction. The nearest-neighbour encounter calculation is short and useful: median waiting time falls with central speed and rises with speed variance. Agent-based runs with log-normal speeds reproduce the qualitative sorting, the ageing of the slope, and the residual solitary-cell histogram once you match mean chain length. Code and processed data are public.\n\nSoft spots, in proportion. The quantitative overlay in Fig. 1e is not free. V is taken from the already-aggregated experimental weighted average; σ_V is then fitted so that the simulated curve at the time when ⟨n⟩≈2.2 matches experiment. The paper itself notes that the solitary-cell distribution is the residual after fast cells have already joined chains and that experimental t=0 is ill-defined. So the match is a consistency check under the dry assumptions (identical drag, instantaneous flagellar reorientation, no hydrodynamics or elasticity once the chain forms), not an independent forecast of the pre-aggregation speed distribution. If hydrodynamics systematically lowered drag with n, or if the true initial σ_V were narrower, you could get a similar curve without pure self-sorting. Those omissions are acknowledged and do not kill the qualitative claim, but they keep the quantitative claim modest.\n\nMath and citations look solid; the dry model is derived cleanly in the appendix and the literature on LC bacteria and ballistic aggregation is properly placed. This is for people who work on active matter in anisotropic media or bacterial self-organization. It deserves a serious referee. I would engage with it and expect to cite the sorting idea.","headline":"Clear experimental observation of faster longer chains, explained by relative-velocity self-sorting; quantitative match is a one-parameter consistency check under a dry model, not an independent prediction.","tokens_in":24907,"tokens_out":515,"would_cite":true,"duration_ms":5486,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Longer bacterial chains swim faster because faster cells preferentially join them, a self-sorting process the authors call swimming-limited aggregation.","keywords":["swimming-limited aggregation","bacterial chains","nematic liquid crystals","self-sorting","speed variance","E. coli","one-dimensional aggregation","active matter"],"falsifier":"Measure the length–speed curve of chains formed from a genetically or environmentally prepared population whose swimming-speed variance is much narrower than the wild-type distribution used here; if the positive correlation vanishes or becomes non-monotonic while mean speed remains unchanged, the self-sorting claim fails.","tokens_in":24883,"feed_emoji":"🦠","tokens_out":543,"duration_ms":5072,"temperature":0.7,"pith_summary":"Escherichia coli swimming in one-dimensional lanes inside a nematic liquid crystal spontaneously form long-lived chains. Force balance alone would make every chain swim at the average speed of its members, yet experiments show that longer chains are systematically faster. The authors resolve the paradox by showing that chains are not random assemblies: cell-to-cell speed differences make faster bacteria catch slower ones first, so the fastest cells are preferentially incorporated into longer chains. A minimal three-body encounter model and agent-based simulations of thousands of particles with log-normally distributed speeds reproduce the measured length–speed curve and its dependence on speed variance. The same framework predicts that the correlation weakens as the system ages. The result identifies a distinct aggregation regime whose kinetics and statistics are controlled by the intrinsic speed distribution of the swimmers rather than by diffusion or elastic attraction.","feed_headline":"Faster bacteria join longer chains, making them swim faster","feed_subtitle":"Speed differences, not random collisions, sort cells into chains whose length tracks average speed.","key_machinery":"Swimming-limited aggregation: a one-dimensional process in which irreversible chain formation is driven solely by relative swimming speeds (sampled from a lognormal distribution) under a dry force-balance rule that sets every new chain’s speed to the arithmetic mean of its constituents.","core_discovery":"Longer chains of E. coli swimming along the director of a nematic liquid crystal move systematically faster than shorter ones because of a dynamic self-sorting process: relative velocities cause faster bacteria to meet and join chains earlier, so longer chains become enriched in high-speed cells. This swimming-limited aggregation is constrained by the variance of the individual speed distribution and becomes weaker with time.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Longer E. coli chains swim faster via speed-based self-sorting","Faster bacteria join early, enriching long chains that outpace short ones","Speed variance drives bacterial chains to grow faster as they lengthen","Relative velocities sort swimmers so longer chains gain higher average speed","Self-assembled bacterial chains speed up with length due to dynamic sorting"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The model assumes every bacterium has the same drag coefficient and that flagella rearrange instantly after collisions, so chain speed is always exactly the average of the members’ intrinsic speeds.","fun_headline_variants_meta":{"raw":{"variants":["Longer E. coli chains swim faster via speed-based self-sorting","Faster bacteria join early, enriching long chains that outpace short ones","Speed variance drives bacterial chains to grow faster as they lengthen","Relative velocities sort swimmers so longer chains gain higher average speed","Self-assembled bacterial chains speed up with length due to dynamic sorting"]},"model":"grok-4.5","effort":"low","cost_usd":0.004172,"raw_usage":{"total_tokens":1280,"prompt_tokens":781,"num_sources_used":0,"completion_tokens":95,"cost_in_usd_ticks":41720000,"prompt_tokens_details":{"text_tokens":781,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":404,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":781,"tokens_out":95,"duration_ms":3928,"temperature":1.0,"reasoning_tokens":404,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T07:33:10.197052+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Measure the length–speed curve of chains formed from a genetically or environmentally prepared population whose swimming-speed variance is much narrower than the wild-type distribution used here; if the positive correlation vanishes or becomes non-monotonic while mean speed remains unchanged, the self-sorting claim fails.","supporting_citations":[],"review_version":2}