{"id":"53396c4a-897b-4cd0-94f5-267b6b181dd7","arxiv_id":"2506.14025","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A flow-interaction model of in-line flapping flyers predicts quantized spacing, spring-like bonds, and resonantly amplified waves that limit passive group size to roughly 4 to 9 members.","lead":"This paper presents a mathematical model of flying formations in which each flapping flyer feels the wake of the flyer ahead. The model reproduces several lab experiments and predicts that passive groups form a fragile 'crystal' that transmits growing waves, limiting group size to about 4 to 9 members.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Erase-and-replace wake assumption is load-bearing; if wakes advect or superpose instead of being overwritten, the quantized lattice and flonon cascade lack support.","rationale":"The reader's weakest assumption correctly identifies the erase-and-replace rule, and I agree that it is the most load-bearing point. The rest of the model--Newtonian balance, the 3/2-power drag law, state-dependent delays, and the steady-theory derivations--is internally consistent, and the comparisons with published two-flyer and five-flyer experiments give real support to the modeling program. The analytical expressions in Eqs. (20)-(21), (25)-(26), and (35) are derived rather than fitted, and the fitted drag coefficient calibrates the isolated-flyer speed curve without directly determining the structural predictions. The problem is that the structural predictions all pass through the assumption that a flyer's wake replaces the upstream wake. That assumption is stated explicitly but is not derived from the fluid mechanics, and it is not strongly constrained by the experiments: two-flyer tests have only one relevant wake, while five-flyer tests can be matched without establishing the transport mechanism. Because the nearest-neighbor reduction is exactly what makes the flonon resonance cascade pairwise and sequential, and because the maximum-group-size predictions in Sec. V are consequences of that cascade, the erase-and-replace rule is load-bearing. If a no-erase superposition variant changes the equilibrium spacing or removes downstream amplification, the central claim about passive flow-mediated matter would need to be weakened. If it does not, the concern is resolved and the conditional verdict can stand or be upgraded. I therefore recommend no change to the CONDITIONAL verdict.","tokens_in":34818,"tokens_out":5665,"duration_ms":67808,"concrete_test":"Simulate an N>=3 array with a no-erase variant of Sec. II.A: replace the wake signal at flyer n by a superposition of upstream wakes, e.g. W_{n-1}(x,t) = sum_{m<n} V_m(t_{m->n}(t)) e^{-(t - t_{m->n})/tau}, with t_{m->n} defined by X_n(t)=X_m(t_{m->n}), and keep all other equations, parameters, and initial conditions identical to Sec. V. If the equilibrium spacing remains S*=j+1/6 and downstream amplitude gain G>1 persists for N=5, the erase-and-replace rule is not uniquely load-bearing. If the spacing shifts or flonon gain vanishes, the central material analogy fails without that rule.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's material analogy rests on the erase-and-replace wake rule in Sec. II.A: each flyer overwrites the wake signal of its upstream neighbor with its own signal, so interactions are nearest-neighbor and one-way downstream. This is what reduces Eq. (3) to a single delayed term V_{n-1}(t_n(t))e^{-(t-t_n)/tau}. From that term come the equilibrium condition (Eq. 20), the quantized spacing S*=j+1/6 (Eq. 21), the Hookean spring constant (Eq. 25), and the pairwise resonance-cascade explanation of flonons in Sec. IV.D, which explicitly treats each pair as isolated because the 2-to-3 bond cannot influence flyer 2. If real wakes are not erased but advect past several downstream members and superpose, a follower should feel a sum of wakes from all upstream flyers, with delays set by when each upstream flyer occupied the follower's current position, rather than only the immediate predecessor's signal. The interaction graph remains one-way, but it becomes long-range, so the nearest-neighbor chain reduction and the sequential-pair gain argument no longer follow. The validations in Secs. III.C-III.E do not isolate this rule: pairwise experiments cannot distinguish erasure from superposition, and the five-flyer experiment is compared after the fact rather than used to test the wake-transport mechanism. The maximum-group-size, flonon, and fragility predictions in Sec. V are therefore conditional on an unvalidated kinematic overwrite rule.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a delay-differential-equation model for in-line formations of flapping flyers, in which each flyer's wake is \"erased and replaced\" by the next downstream flyer, so that interactions are nearest-neighbor and one-way. The model uses a skin-friction drag scaling as U^{3/2} and a state-dependent delay tied to the time an upstream flyer occupied the follower's current position. The paper validates the model against a series of prior robophysical experiments: a single isolated flyer, a self-interacting flyer in a cyclic domain, two flyers with identical and with distinct kinematics, and a five-flyer group. Analytical steady-state analysis yields quantized equilibrium spacings S* = j + 1/6, a Hookean spring constant k, and a resonance frequency f_R; numerical simulations show downstream-amplifying longitudinal waves (\"flonons\") and a maximum cohesive group size of roughly 4 to 9 flyers in the strong-interaction (tau -> infinity) limit. The authors interpret these results as supporting a view of flying formations as flow-mediated matter.","tokens_in":35208,"tokens_out":11036,"duration_ms":118951,"significance":"If the model is accepted as a faithful minimal phenomenology, the paper makes a valuable contribution by providing an analytically tractable framework that connects individual flapping kinematics to collective ordering, elasticity, and instability. The derivation of the quantized spacing and spring constant directly from the delay equations, rather than by fitting to aggregate data, is a genuine strength, as is the reproduction of hysteresis, force-displacement curves, and flonon amplification without adjusting the interaction rule per experiment. The predictions of maximum stable group size and fragility maps are falsifiable and potentially useful for interpreting animal group behavior. However, the central physical mechanism is the erase-and-replace wake rule, which is assumed rather than measured; all subsequent claims inherit this assumption. The paper would be significant for the modeling community even if the wake-transport rule turns out to be an oversimplification, provided the authors clearly delimit the domain of validity and test the sensitivity of the headline predictions to the rule.","major_comments":[{"comment":"The erase-and-replace wake scheme is the load-bearing assumption of the paper: it reduces Eq. (3) to a single delayed term and leads directly to the nearest-neighbor, one-way interaction graph on which the quantized spacing (Eq. 20), the spring constant (Eq. 25), and the pairwise resonance-cascade explanation of flonons (Sec. IV.D) all rest. The manuscript, however, provides no direct experimental or numerical evidence that a follower's wake erases the upstream signal, as opposed to superposing with it after advection past several downstream members. The pairwise validations cannot distinguish these alternatives, and the five-flyer comparison in Sec. III.E is an a posteriori reproduction rather than a designed test of the transport rule. I request either (i) a sensitivity study in which the five-flyer case is simulated with a superposition or finite-range-advection wake rule, (ii) a direct comparison to wake-flow measurements if available, or (iii) an explicit, detailed statement that the erase-and-replace rule is an untested modeling postulate, together with a discussion of which observable would most clearly falsify it.","section":"Sec. II.A and Sec. III.E"},{"comment":"The manuscript reports that the theoretical tensile strength overestimates the numerical values by about 40% and the tensile strains by about 60%, attributing the discrepancy to unsteady terms. Because the material analogy is a central theme, this quantitative mismatch should be addressed rather than stated in passing. Please provide confidence intervals or error bars for the numerical force-displacement data, and either (i) derive a first-order unsteady correction to Eq. (29) that reduces the gap, or (ii) soften the claims of quantitative agreement and explicitly state that only the compressive side is well captured by the steady theory.","section":"Sec. IV.C (tensile strength)"},{"comment":"The text says the 3/2-power drag law is \"crucial\" for recovering the experimentally observed hysteresis loops, but no counterfactual simulation with the previously used quadratic drag U^2 is shown. Since this is an empirical modeling claim that motivates the entire model revision, please provide a direct comparison: run the same cyclic self-interacting setup with a quadratic drag term and show whether the multi-loop hysteresis disappears, or revise the wording to \"consistent with\" rather than \"crucial.\"","section":"Sec. III.B (drag law and hysteresis)"},{"comment":"All experimental comparisons in Secs. III.A-III.E are plotted without error bars, and in Figs. 5(c) and 9 the experimental data points are not shown at all (only cited in the text). To support the claim that the model \"faithfully reproduces\" the experiments, please overlay the experimental measurements with their reported uncertainties, or, if the original papers do not provide uncertainties, state that explicitly and discuss the resulting limitations on the validation.","section":"Sec. III (experimental comparisons)"}],"minor_comments":[{"comment":"The inequality is written as \"2πS* tan(2πS*) ln(1/(2cos(2πS*)) > 1\", which appears to be missing a division sign; it should read [2πS* tan(2πS*)] / ln(1/(2cos(2πS*))) > 1.","section":"Sec. IV.B, Eq. (27)"},{"comment":"The caption contains the typo \"single.yer\" in the dotted cyan line description; it should be \"single flyer.\"","section":"Fig. 4(b) caption"},{"comment":"The notation U1,2(t), V1,2(t), and V2,1(t1,2(t)) in the closed two-flyer equation is confusing; please define the index conventions explicitly (e.g., which subscript refers to the source and which to the receiver).","section":"Sec. II.B, Eq. (5)"},{"comment":"The assumption that the wake-flow speed at a flyer's location equals its flapping speed (W_n(X_n(t), t) = V_n(t)) is introduced without a citation or empirical justification; a brief note on the experimental or theoretical basis for this equality would help the reader judge its validity.","section":"Sec. II.A"},{"comment":"The allowed wave numbers in the nonreciprocal chain are stated as \"q_m = mπ/λ or q_m = (2m±1/2)π/[(N−1)λ]\" without a derivation; please clarify the boundary-condition setup (especially the role of the ghost point and the fixed-free conditions) and show how these values follow from the determinant condition.","section":"Sec. V.B, Eqs. (40)-(42)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a modeling study that heavily cites the authors' own prior experimental and theoretical work (refs. 19, 21, 24, 25, 26, 47). This is not inappropriate for a follow-up modeling paper, but it places extra weight on the five-flyer comparison being a genuine test rather than a post-hoc fit. The erase-and-replace assumption is the key risk; if the authors can demonstrate robustness to relaxation of this rule, the paper would become substantially stronger. The scope is appropriate for physics.flu-dyn, though the material-science framing may appeal to a broader readership."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a solid modeling extension of the Ristroph group's own previous work, and it lands some genuinely new analytic results. The 3/2-power drag revision is a real improvement with observable consequences, and the analytic formulas for the equilibrium spacing (S* = j + 1/6), the spring constant, the resonance frequency, and the stability bounds for dissimilar amplitudes are all cleanly derived. The model reproduces the main qualitative features of the earlier robophysical experiments: hysteresis, quantized spacing, force-displacement curves, and the flonon-like amplification. That is real evidence the model is not just a toy.\n\nThe soft spots are in the assumptions and the validation. The erase-and-replace wake scheme in Sec. II.A is doing load-bearing work: it reduces interactions to nearest-neighbor one-way, gives the equilibrium condition, and makes the flonon cascade a sequential-pair argument. But the paper never tests that physical picture against wake-transport data. Pairwise experiments cannot distinguish erasure from superposition, and the five-flyer comparison is post hoc rather than a designed test. If real wakes advect past several downstream members and superpose, the lattice constant and the gain argument would both change. This is not a fatal flaw, but it is a real gap. The tensile-strength theory overestimates numerics by about 40%, and the paper explains why but does not resolve it. The drag coefficient CD is fitted to single-flyer data, and there are no error bars on the experimental comparisons or any released code, which makes it hard to judge how much of the quantitative agreement is tuning.\n\nStill, the analytic core is internally consistent, and the paper is honest about its limitations. This is a paper for people working on flapping locomotion, flow-mediated interactions, or active matter analogies. It deserves a serious referee: the new formulas are useful enough that a careful treatment of the wake-transport assumption would materially strengthen it. I would send it to peer review.\n\nRecommendation: send it, but the referee should ask for a sensitivity test that relaxes erase-and-replace, even in a simplified way, to see whether the qualitative predictions survive.","headline":"A useful modeling extension with clean analytics, but the erase-and-replace wake rule carries a lot of weight and deserves a sensitivity check before the material analogy is taken as established.","tokens_in":35681,"tokens_out":2807,"would_cite":true,"duration_ms":32254,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76Z10","34K60"],"pacs":["47.63.mf"],"model":"deepseek-v4-flash","headline":"In-line flapping flyers self-organize into a soft crystal with quantized spacing and spring-like wake bonds, while self-amplifying 'flonon' waves limit cohesive groups to roughly 4-9 members.","keywords":["formation flight","flow-mediated interactions","flapping foils","wake memory","delay differential equations","flonons","soft crystals","group cohesion"],"falsifier":"In a three-foil in-line experiment, independently vary the last foil's flapping amplitude or frequency and measure the leader's speed: the model predicts zero upstream response, so any measurable change in leader speed falsifies the one-way nearest-neighbor assumption. Conversely, measure the distribution of gaps in a five-foil array: the model predicts peaks at $S=j+1/6$ (gaps of about $1.2\\lambda$); a gap distribution without a peak at 1.2 times the trajectory wavelength would falsify the lattice constant.","tokens_in":34632,"feed_emoji":"🕊️","tokens_out":7553,"duration_ms":71954,"temperature":0.7,"pith_summary":"This paper argues that the orderly formations seen in flying and swimming groups can be understood as a material assembled by flow forces alone. The authors revise a follower-wake model of flapping flyers and show it reproduces a decade of robophysical experiments on pairs and small groups, including discrete stable spacings, hysteresis, and force-spacing curves that behave like spring-like bonds. In the strong-interaction limit the model predicts that in-line arrays form a soft crystal: members settle at quantized gaps of about 1.2 wake wavelengths, bonds have a measurable stiffness, and disturbances travel down the line as longitudinal waves ('flonons') that grow in amplitude, eventually colliding neighbors and fracturing the formation. Because the interactions are one-way and nearest-neighbor, the amplification cascades pair by pair, and the model gives a maximum cohesive group size of roughly 4 to 9 flyers without feedback control. If right, passive flow-mediated physics sets both the order and the fragility of animal collectives, and active sensing would be needed to stabilize long natural formations.","feed_headline":"Wakes alone order flyers into soft crystals of 4 to 9","feed_subtitle":"Model reproduces lab experiments and predicts quantized spacing, springy bonds, and self-amplifying waves.","key_machinery":"The central object is the follower-wake interaction model (Eqs. 2-4), a system of nonlinear, state-dependent delay differential equations in which each flyer's thrust depends on $(V_n - V_{n-1}(t_n)e^{-(t-t_n)/\\tau})^2$, with memory time $t_n$ defined implicitly by $X_n(t)=X_{n-1}(t_n)$ and obeying its own evolution equation. The erase-and-replace wake scheme makes interactions strictly nearest-neighbor and one-way downstream. The analysis then reduces the pairwise equilibrium to $\\cos(2\\pi S^*) = e^{-S^*/(f\\tau)}$, whose $\\tau\\to\\infty$ limit gives $S^*=j+1/6$; linear stability gives the spring constant $k$, and the driven-damped-oscillator gain factor explains flonon amplification. A dispersion relation $\\Omega(q)=\\sqrt{k/M}(1-e^{iq\\lambda})^{1/2}$ for longitudinal waves in a one-way (diodic) mass-spring chain supplies the group speed. The flapping-speed-squared thrust law, the 3/2-power skin-friction drag, and the exponential wake decay are the three physical ingredients that make the balance produce quantization.","core_discovery":"The central claim is that a minimal delay-differential model of wake interactions is sufficient to produce the crystalline order observed in in-line flapping-flyer formations, and to explain its breakdown. Each flyer's thrust is taken to depend on the square of its flapping speed relative to the 'erase-and-replace' wake signal left by the immediately upstream neighbor, which decays exponentially and is encountered after a memory delay that is itself a state variable. Steady-state analysis of a leader-follower pair shows that equilibrium spacings satisfy $\\cos(2\\pi S^*)=1/2$ in the long-lived-wake limit, quantizing the lattice spacing to $S^*=j+1/6$ (dimensionless gap of about 1.2 wake wavelengths), with $j+5/6$ unstable. Linear perturbation theory yields a spring constant $k$ for the inter-flyer bond, a resonant frequency, and compressive and tensile strengths of roughly 20% and 10% strain. The one-way directionality of the bonds lets each pair be treated as an isolated driven damped oscillator, so an oscillatory disturbance is amplified by each successive pair; nonlinear simulations show spontaneous collisions for groups beyond roughly 4-9 members.","pith_inferences":["If the erase-and-replace scheme is a fair idealization of real high-Reynolds-number wakes, the same framework should apply to fish schools in in-line configurations, predicting the same $j+1/6$ quantization with the tail-beat wavelength replacing the flapping wavelength; this is a testable prediction against existing schooling data.","The flonon mechanism implies that a trailing member can be a passive amplifier of disturbances the leader never directly sensed, which could be exploited in engineered swarm sensing rather than suppressed.","The predicted group-size ceiling gives a concrete benchmark for behavioral studies: any natural columnar formation longer than about nine members likely requires active control or kinematic variability that detunes the resonance cascade.","Measuring flonon propagation times in experiments with different body inertias would let one infer the effective spring constant of the flow bond in vivo, offering a non-intrusive probe of the group's mechanical state."],"forward_implications":["In-line formations of flapping flyers should display quantized equilibrium spacings at $S^*=j+1/6$, with the first stable gap at about $1.2\\lambda$, matching prior two-foil experiments.","Each inter-flyer bond should respond like a spring with stiffness $k$ given by Eq. (25), failing in compression at roughly 20% strain and in tension at roughly 10% strain.","A small oscillatory disturbance applied to the leader should grow in amplitude as it travels downstream, with the per-pair gain given by the driven-damped-oscillator resonance curve.","Without feedback control, cohesive in-line groups should be limited to roughly 4-9 members, with the largest groups at small dimensionless mass and flapping Reynolds number, and instability times of 2-20 flaps.","Information propagates downstream at a speed bounded by the flight speed at low mass and by the non-reciprocal wave speed at high mass."],"supporting_citations":[{"why":"Supplies the two-flyer robophysical experiments whose quantized stable spacings, force-spacing curves, and stability map the model reproduces.","marker":"[19]"},{"why":"Supplies the self-interacting flyer experiments with hysteresis and multiple flight modes that fix the drag law and validate the cyclic-boundary model.","marker":"[21]"},{"why":"Supplies experiments with two differently flapping flyers whose stable positions, collisions, and separations give the model's phase diagram.","marker":"[24]"},{"why":"Supplies the few-flyer experiments showing ordered formations disrupted by self-amplifying waves, the empirical target for 'flonons' and group fragility.","marker":"[25]"},{"why":"Provides the discrete-time lattice theory of hydrodynamically interacting flapping swimmers whose bistability analysis the steady theory extends.","marker":"[26]"},{"why":"Provides the thrust scaling and optimal Strouhal range used to justify the flapping-speed-squared thrust law.","marker":"[30]"},{"why":"Provides measurements fixing the thrust coefficient used in the model.","marker":"[39]"},{"why":"Provides the Blasius boundary-layer drag law whose 3/2-power speed scaling the model adopts to reproduce hysteresis.","marker":"[42]"},{"why":"Supplies the non-reciprocal metamaterial continuum whose dispersion relation is adapted for the diodic wave speed at high mass.","marker":"[64]"}],"fun_headline_variants":["Wake interactions crystallize flyer formations","Flying formations behave like soft crystals","Flow-mediated bonds turn flyers into crystals","Quantized spacing emerges from wake interactions","Flapping flyers form soft crystals via wakes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the erase-and-replace wake rule: a flyer's wake fully overwrites its upstream neighbor's at the position of the downstream flyer, so interactions are only between nearest neighbors, only downstream, and the wake speed equals the flapping speed.","fun_headline_variants_meta":{"raw":{"variants":["Wake interactions crystallize flyer formations","Flying formations behave like soft crystals","Flow-mediated bonds turn flyers into crystals","Quantized spacing emerges from wake interactions","Flapping flyers form soft crystals via wakes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000412,"raw_usage":{"total_tokens":2199,"prompt_tokens":1080,"completion_tokens":1119,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":696,"completion_tokens_details":{"reasoning_tokens":1055}},"tokens_in":696,"tokens_out":1119,"duration_ms":9709,"temperature":1.0,"reasoning_tokens":1055,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:24:58.466171+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a three-foil in-line experiment, independently vary the last foil's flapping amplitude or frequency and measure the leader's speed: the model predicts zero upstream response, so any measurable change in leader speed falsifies the one-way nearest-neighbor assumption. Conversely, measure the distribution of gaps in a five-foil array: the model predicts peaks at $S=j+1/6$ (gaps of about $1.2\\lambda$); a gap distribution without a peak at 1.2 times the trajectory wavelength would falsify the lattice constant.","supporting_citations":[{"cited_title":"Lighthill.Mathematical biofluiddynamics","cited_arxiv_id":null,"evidence_quote":"Supplies the two-flyer robophysical experiments whose quantized stable spacings, force-spacing curves, and stability map the model reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the self-interacting flyer experiments with hysteresis and multiple flight modes that fix the drag law and validate the cyclic-boundary model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies experiments with two differently flapping flyers whose stable positions, collisions, and separations give the model's phase diagram."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the few-flyer experiments showing ordered formations disrupted by self-amplifying waves, the empirical target for 'flonons' and group fragility."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the discrete-time lattice theory of hydrodynamically interacting flapping swimmers whose bistability analysis the steady theory extends."},{"cited_title":"On the stability of an in-line formation of hydrodynamically interacting flapping plates","cited_arxiv_id":"2410.04626","evidence_quote":"Provides the thrust scaling and optimal Strouhal range used to justify the flapping-speed-squared thrust law."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides measurements fixing the thrust coefficient used in the model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Blasius boundary-layer drag law whose 3/2-power speed scaling the model adopts to reproduce hysteresis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the non-reciprocal metamaterial continuum whose dispersion relation is adapted for the diodic wave speed at high mass."}],"review_version":1}