{"id":"5308c409-b1bc-4578-8dc8-42662f9f039c","arxiv_id":"2502.04316","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Broken fore-aft symmetry in a swimming cell produces an asymmetric rotation rate in nonuniform shear, driving net cross-stream migration toward channel walls.","lead":"A model of flagellated microswimmers with head-tail asymmetric shape shows that activity and shape asymmetry make them cross flow streamlines in a channel, drifting toward the walls instead of staying on one streamline. This explains why swimming bacteria can move toward channel walls, a behavior seen in experiments but not captured by simpler symmetric-swimmer models, and suggests flexibility can be used to sort swimmers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Straight-rod flagellum may be the sole source of the HTA rotation-rate asymmetry; a helical-flagellum test is needed to validate the generic mechanism.","rationale":"Good-faith reading: the paper is a mechanistic modeling study. Its central claim is not just that this specific model migrates; it is that broken fore-aft symmetry fundamentally drives active cross-streaming. What would have to be true is that the asymmetric rotation rate is a generic property of HTA swimmers in nonuniform shear, not a special property of a point-attached straight rod. The manuscript identifies the rod as a minimal representation, and the conclusion extends beyond it; the Supplemental Material is referenced for the coefficients, so the derivation cannot be checked here. The straight-rod simplification is therefore the weakest load-bearing element. The reduced oscillator theory is a helpful secondary support, but its comparison to simulation uses a normalization factor y_LC^π taken from the simulation (Appendix A), so it cannot independently validate the amplitude of the mechanism. A helical-flagellum computational test would settle whether the mechanism is generic. This does not amount to a rejection; the model is self-consistent and the qualitative argument about asymmetric time spent facing the wall is plausible, so we keep the reader's CONDITIONAL verdict.","tokens_in":13573,"tokens_out":4372,"duration_ms":46028,"concrete_test":"Run the same Stokes-flow model with the flagellum replaced by a rotating helical filament (e.g., N point forces on a helical centerline with the same total active force and the same cell body), recompute the coefficients β1, β3 and r_sh in Eq. (6), and track trajectories for k=20 at yi/R=0.7. If migration toward the wall vanishes or reverses, the straight-rod representation is load-bearing; if migration persists with comparable β1, β3, the mechanism is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the flow-induced asymmetric rotation rate in Eq. (6), where fore-aft asymmetry enters exclusively through the coefficients β1, β3 and the hydrodynamic-center shift r_sh. These coefficients are computed for a flagellum modeled as a single straight slender rod attached at the sphere surface and hinged at one bend angle ϕ with linear elastic torque. For ℓf→0 the model reduces to an HTS swimmer with β1=β3=r_sh=0, so the entire predicted migration is generated by the straight-rod ansatz. A real flagellar bundle is a rotating helix with distributed flexibility and chirality; the straight rod imposes a particular relation between body orientation and the shear-gradient direction, and the non-uniform rotation rate may be an artifact of that simplified geometry. The paper's generalized conclusion that HTA shape fundamentally drives cross-streaming would be falsified if a helical, still-HTA flagellum produced β1, β3 of different sign, or coefficients dominated by chirality instead of fore-aft asymmetry. Because the Supplemental derivation of β1, β3 is not visible in this version, the robustness of the asymmetry to rod geometry is the key unverified step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a numerical and analytical study of a model flagellated microswimmer—a spherical body with a slender-rod flagellum that can bend at a single hinge—in a pressure-driven Poiseuille flow. Solving the Stokes-flow force and torque balance with image singularities, the authors find that head-tail asymmetric swimmers exhibit three dynamical classes: migration to the walls, escape from the mid-channel region, and indefinite trapping near the mid-channel, depending on initial conditions. A reduced two-variable model (Eqs. (5)–(6)) introduces coefficients β1, β3, and r_sh that parameterize the fore-aft asymmetry, and a van der Pol-type amplitude equation (Eq. (7)) predicts an unstable limit cycle whose stable counterpart traps a population near the center. The authors argue that the asymmetry-induced nonuniform rotation rate, combined with self-propulsion, drives cross-stream migration.","tokens_in":13868,"tokens_out":6244,"duration_ms":61751,"significance":"If the mechanism is generic, this is a significant advance: it identifies a purely viscous, inertialess cross-stream migration mechanism for active particles, distinct from inertia, viscoelasticity, or wall repulsion, and it suggests a practical way to sort swimmers by flagellar flexibility. The paper's strengths are the fully resolved hydrodynamics, the explicit reduced theory with a phase-space explanation, and the comparison with existing experimental parameters in Appendix C showing where HTA effects should have been visible. The central predictions—wall-directed migration, a center-trapped population, and a flexibility-dependent limit cycle—are falsifiable.","major_comments":[{"comment":"The manuscript's generalized conclusion that 'head-tail shape-asymmetry fundamentally drives active cross-streaming' is tested only for a single geometric model: a straight slender rod hinged at one angle. Since all fore-aft asymmetry enters through β1, β3, and r_sh in Eqs. (5)–(6), and these vanish in the ℓf→0 limit, the straight-rod ansatz is the exclusive source of the predicted migration. To rule out that the rotation-rate asymmetry is an artifact of the straight-rod representation, the authors should either compute the same coefficients for a helical (chiral, distributed-flexibility) flagellum or for a different HTA shape (e.g., a two-sphere swimmer) and show that the sign of β1/β3 is shape-independent, or soften the claim to 'this minimal rod model'.","section":"§2, after Eq. (4); Eq. (6)"},{"comment":"The quantitative validation of the reduced theory is weakened by the use of simulation-derived normalizations. The analytical trajectory is plotted as y/y_LC^π with y_LC^π taken from the simulation, so the absolute amplitude of the oscillation is not predicted but matched. Similarly, the 10% 'trapped fraction' estimate in Appendix B uses the limit-cycle dimensions a_max and b_max measured from Fig. A1(c). The paper should either derive y_LC^π from the theory (e.g., via the LC ellipse formula) or clearly state that the theory predicts only the relative dynamics.","section":"Fig. A1(a) and Appendix B"},{"comment":"The coefficients β1, β3 and the reduced coefficients μ1, μ3, k1, k3 are essential to the mechanism, but their derivations are relegated to Supplemental Material [51] and were not available for review. The sign of β1 determines the direction of migration, and Eq. (7) is asserted to follow from eliminating y from Eq. (6) without showing the algebra. For the central claim to be checkable, these derivations must be included in the supplement and made available.","section":"Eqs. (5)–(7) and Supplemental Material"}],"minor_comments":[{"comment":"The sentence 'the shape-HT symmetry is purely broken by the β1 and β3 terms alone' is confusing because r_sh in Eq. (5) is also zero for HTS swimmers; please clarify that the statement refers to the θ̇ equation.","section":"Eq. (6)"},{"comment":"The expression for ψ_a(t) is missing a closing parenthesis in the denominator; please correct the typesetting.","section":"Appendix A, Eq. (A4)"},{"comment":"The text refers to 'the other coefficients are listed in [51]', but the Supplemental Material was not provided with this version; please ensure it is submitted with the revision.","section":"References, item [51]"},{"comment":"The unstable limit cycle (orange) is difficult to distinguish from the separatrix (gray) in the printed figure; using different line styles or colors would improve readability.","section":"Fig. 2(c)"},{"comment":"The inset labels (Panel-1, Panel-2, Panel-3) are not defined in the caption; please add a short description of each class.","section":"Fig. 2(b) caption"}],"recommendation":"major_revision","confidential_remarks":"This is a promising manuscript with a well-posed research question and a physically plausible mechanism. The main concerns are (i) the lack of access to the Supplemental Material, which contains the coefficient derivations, and (ii) the absence of a test of the robustness of the asymmetry to flagellar geometry. Both are fixable within a revision. The paper fits the scope of cond-mat.soft and is likely to be of interest to the active-matter community. I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one if you care about active particles in channel flows. The paper shows that a minimal flagellated microswimmer with broken fore-aft symmetry migrates across streamlines in a Poiseuille flow at low Reynolds number, purely through the asymmetric rotation rate the swimmer experiences in nonuniform shear. That mechanism is new—previous explanations needed inertia, viscoelasticity, walls, field alignment, or unsteady beating. The model is built from Stokes-flow force balance with no fitted parameters, and the reduced two-equation theory captures the phase-space structure convincingly: a stable fixed point at upstream orientation, an unstable limit cycle, and a van der Pol-like amplitude equation. The analytic trajectories match the full simulations well, which gives real confidence. \n\nThe main soft spot is exactly the one the stress-test flags. The entire asymmetry lives in the coefficients β1, β3, and the hydrodynamic-center shift r_sh, all computed for a single straight slender rod hinged at the cell body. The model reduces to an HTS swimmer with zero migration when ℓf→0, so the mechanism is generated entirely by that rod ansatz. A helical, rotating flagellum with distributed flexibility could plausibly change the sign or magnitude of these coefficients, and the paper's conclusion that 'head-tail shape-asymmetry fundamentally drives' cross-streaming is broader than the evidence supports. That is not a fatal flaw—the paper contributes a clear, testable mechanism in a well-defined model—but the generality claim needs either tempering or a test with a more realistic flagellum. \n\nSmaller issues: the full derivation and solver details are in the supplemental material and not visible in the preprint; the analytic comparison in Fig. A1(a) normalizes by the simulation's limit cycle value, which slightly masks the predictive power; and the 10% flushed estimate is an area ratio of an assumed ellipse, so it should be read as a rough estimate. None of these change the central result. \n\nThis paper deserves a serious referee. The revision should make the supplemental derivation accessible and ask the authors to discuss how sensitive the β coefficients are to the straight-rod idealization. If the mechanism survives a helical-flagellum test, it is an important result for bacterial rheotaxis and microfluidic sorting. I'd bring it to reading group.","headline":"A clean mechanistic result on how head-tail asymmetry in a model flagellated swimmer produces cross-stream migration; the main open question is whether the straight-rod flagellum is generic.","tokens_in":14375,"tokens_out":2790,"would_cite":true,"duration_ms":25747,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A head-tail asymmetric microswimmer in a purely viscous channel flow migrates toward the walls because its shape makes it spend more time pointing toward the wall than toward the center.","keywords":["microswimmer","cross-stream migration","head-tail asymmetry","flagellar flexibility","Poiseuille flow","low Reynolds number","active matter","limit cycle"],"falsifier":"In a microchannel Poiseuille flow at low Reynolds number and far from walls, track the orientation angle $\\theta(t)$ of a head-tail-asymmetric swimmer (e.g., an engineered colloid with a rigid asymmetric appendage) over many oscillation periods; if the time it spends pointing toward the wall does not exceed the time it spends pointing toward the channel center, the predicted wall-directed drift cannot occur, and the central mechanism is refuted.","tokens_in":13374,"feed_emoji":"🦠","tokens_out":9177,"duration_ms":82897,"temperature":0.7,"pith_summary":"The paper establishes that a microswimmer with broken head-tail symmetry migrates across streamlines in a purely viscous channel flow, without needing inertia, viscoelasticity, or walls. The mechanism is that the asymmetric shape makes the swimmer rotate faster when pointing toward the channel center than toward the wall, so it spends more time pointing toward the wall and is carried there by its own propulsion. Using a minimal model of a spherical body with a flexible flagellar rod, the authors solve the low-Reynolds-number hydrodynamics self-consistently and show that most such swimmers reach the channel walls, while a small central population near an unstable limit cycle stays trapped at the midsection and is flushed out. Flagellar flexibility accelerates the wall migration and, past a critical stiffness, eliminates the trapped population entirely. The result matters because most real microswimmers are head-tail asymmetric, so cross-stream migration should be common in microchannels rather than an exotic effect.","feed_headline":"Head-tail shape alone steers microswimmers to walls","feed_subtitle":"Broken fore-aft symmetry makes swimmers rotate unevenly, carrying them across streamlines.","key_machinery":"The key object is the rotation-rate equation for the orientation angle $\\theta$: $\\dot{\\theta} = (v_f/R^2)(\\alpha y + 2\\beta_2 \\sin^2\\theta \\, y + \\beta_1 \\sin\\theta + \\beta_3 \\sin^3\\theta)$, where $y$ is the swimmer's vertical position in the channel and the constants $\\beta_1$, $\\beta_3$ and the hydrodynamic-center shift $r_{\\text{sh}}$ are nonzero only when head-tail symmetry is broken. The $\\sin\\theta$ and $\\sin^3\\theta$ terms make the swimmer rotate at different average rates when it points toward the wall than toward the center, and this asymmetry, combined with self-propulsion at speed $v_{\\text{sp}}$, produces the vertical drift in Eq. (5). The second central object is the van der Pol-like amplitude equation $\\ddot{\\psi} + \\mu_1(1-\\zeta\\psi^2)\\dot{\\psi} + k_1\\psi = -\\mu_3\\psi\\dot{\\psi}^2 - k_3\\psi^3$ obtained by linearizing around the fixed point $(0,\\pi)$; its unstable fixed point at amplitude $\\psi^a_{\\text{LC}}$ corresponds to the unstable limit cycle that bounds the center-trapped region.","core_discovery":"The central claim is that head-tail shape-asymmetry is a fundamental driver of active cross-streaming in channel flows. In the model, a pusher with a rigid flagellum in Poiseuille flow shows three trajectory classes: swimmers starting far from the center oscillate toward the nearest wall; those starting at intermediate heights are temporarily trapped at the midsection and then escape to a wall; and those below a height $y_{\\max}$ oscillate with decaying amplitude and remain trapped at the channel center indefinitely. The phase space is organized by an unstable limit cycle that separates the trapped and escaping populations. The underlying mechanism is the asymmetric rotation rate: the swimmer spends a larger fraction of each oscillation period pointing toward the wall, so its self-propulsion produces a net drift toward the wall. The simplified dynamics reduce to a van der Pol-like oscillator around the stable fixed point $(y,\\theta)=(0,\\pi)$, with the asymmetry entering through coefficients $\\beta_1$, $\\beta_3$ and the hydrodynamic-center shift $r_{\\text{sh}}$. Flagellar flexibility increases the average rotation rate, and at a critical rigidity $k_c$ the limit cycle collapses in a subcritical Hopf bifurcation.","pith_inferences":["We infer that head-tail asymmetry may be the dominant cross-stream migration mechanism for bacteria-sized swimmers in typical microchannels, since it operates in a purely Newtonian fluid and requires no external fields or wall contact.","A testable extension would be to measure the rotation-rate asymmetry directly: in a linear shear flow, an asymmetric swimmer should show unequal average $\\dot{\\theta}$ for wall-facing versus center-facing orientations, and the difference should scale with the shear rate.","The model's minimal flagellum suggests that in real flagellated bacteria, the helical shape and distributed flexibility could shift the critical stiffness $k_c$ and the size of the trapped population, but the qualitative phase-space structure (unstable limit cycle, center attractor) should survive."],"forward_implications":["In a mixed population with varying flagellar stiffness, more flexible swimmers reach the walls faster while rigid ones remain center-trapped, so the channel flow can sort swimmers by flexibility.","Tuning the maximum flow speed $v_f$ adjusts the size of the unstable limit cycle and thereby the fraction of the population flushed through the channel center, up to about 10% for a rigid-swimmer population.","For pullers, the limit cycle is stable instead of unstable, so trajectories inside it grow in amplitude and escape the center, the opposite of the pusher case.","The dynamics depend on the channel geometry only through the ratio $v_f/R^2$; trajectories overlap in the $y$-$\\theta$ plane for different $R$ and $v_f$ with the same ratio, indicating an invariant limit-cycle size."],"supporting_citations":[{"why":"Establishes the HTS baseline: a symmetric swimmer in Poiseuille flow oscillates about streamlines with no net migration, the behavior this paper contrasts.","marker":"[10]"},{"why":"Supplies the Jeffery rotation coefficient $\\beta_2$ and the HTS limit that the new model recovers when the flagellar rod shrinks to zero length.","marker":"[12]"},{"why":"Shows that passive bacterial-sized flexible flagella do not cross streamlines strongly, providing the comparison that isolates active HTA-driven migration.","marker":"[44]"},{"why":"Provides the experimental channel parameters and observation window (~60 s) used to argue that HTA migration is slow enough to be missed in that experiment.","marker":"[53]"},{"why":"Provides experimental bacterial trajectories in Poiseuille flow that the model reproduces, with the same wall-ward slopes.","marker":"[55]"},{"why":"Slender-body hydrodynamics used to compute the flagellar drag forces and torques in Eqs. (3)-(4).","marker":"[48]"},{"why":"Slender-rod force-velocity relations used in the flagellar model, including boundary effects relevant to the image-field calculations.","marker":"[49]"},{"why":"Supplies the resistivity tensor and singularity method for the Stokes flow field of the model swimmer.","marker":"[50]"},{"why":"Averaging theory used to derive the slow amplitude dynamics of the van der Pol-like oscillator and locate the unstable limit cycle.","marker":"[57]"}],"fun_headline_variants":["Asymmetric shape drives microswimmers to walls","Lopsided swimmers cross streamlines toward walls","Flex and asymmetry reroute microswimmers in flow","Head-tail mismatch steers swimmers off centerline"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The flagellar bundle is modeled as a single slender rod pinned at one point, bending by one angle $\\phi$ with a linear elastic torque and constant active force, so the predicted drift and limit-cycle structure all rest on the assumption that this minimal rod captures the essential shape asymmetry of a real flagellum.","fun_headline_variants_meta":{"raw":{"variants":["Asymmetric shape drives microswimmers to walls","Lopsided swimmers cross streamlines toward walls","Flex and asymmetry reroute microswimmers in flow","Head-tail mismatch steers swimmers off centerline"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000286,"raw_usage":{"total_tokens":1637,"prompt_tokens":855,"completion_tokens":782,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":722}},"tokens_in":471,"tokens_out":782,"duration_ms":7532,"temperature":1.0,"reasoning_tokens":722,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T22:47:43.722713+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a microchannel Poiseuille flow at low Reynolds number and far from walls, track the orientation angle $\\theta(t)$ of a head-tail-asymmetric swimmer (e.g., an engineered colloid with a rigid asymmetric appendage) over many oscillation periods; if the time it spends pointing toward the wall does not exceed the time it spends pointing toward the channel center, the predicted wall-directed drift cannot occur, and the central mechanism is refuted.","supporting_citations":[{"cited_title":"Z¨ ottl and H","cited_arxiv_id":null,"evidence_quote":"Establishes the HTS baseline: a symmetric swimmer in Poiseuille flow oscillates about streamlines with no net migration, the behavior this paper contrasts."},{"cited_title":"Z¨ ottl and H","cited_arxiv_id":null,"evidence_quote":"Supplies the Jeffery rotation coefficient $\\beta_2$ and the HTS limit that the new model recovers when the flagellar rod shrinks to zero length."},{"cited_title":"Tournus, A","cited_arxiv_id":null,"evidence_quote":"Shows that passive bacterial-sized flexible flagella do not cross streamlines strongly, providing the comparison that isolates active HTA-driven migration."},{"cited_title":"Rusconi, J","cited_arxiv_id":null,"evidence_quote":"Provides the experimental channel parameters and observation window (~60 s) used to argue that HTA migration is slow enough to be missed in that experiment."},{"cited_title":"Junot, N","cited_arxiv_id":null,"evidence_quote":"Provides experimental bacterial trajectories in Poiseuille flow that the model reproduces, with the same wall-ward slopes."},{"cited_title":"Pozrikidis, Introduction to Theoretical and Computa- tional Fluid Dynamics, EBSCO ebook academic collec- tion (OUP USA, 2011)","cited_arxiv_id":null,"evidence_quote":"Slender-body hydrodynamics used to compute the flagellar drag forces and torques in Eqs. (3)-(4)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Slender-rod force-velocity relations used in the flagellar model, including boundary effects relevant to the image-field calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the resistivity tensor and singularity method for the Stokes flow field of the model swimmer."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Averaging theory used to derive the slow amplitude dynamics of the van der Pol-like oscillator and locate the unstable limit cycle."}],"review_version":1}