Pith. sign in

REVIEW 3 major objections 3 minor 39 references

Spin-Waves without Spin-Waves: A Case for Soliton Propagation in Starling Flocks

T0 review · 3 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read Starling turns may be solitons, not spin waves.

desk verdict A real experimental puzzle and a plausible nonlinear mechanism, but the claimed parameter window is not verified at the measured wavevectors and requires implausibly high polarization if tightened. read the letter →

arxiv 2505.19665 v2 pith:TBN6GH3Z submitted 2025-05-26 cond-mat.stat-mech cond-mat.softphysics.bio-phq-bio.QM

classification cond-mat.stat-mechcond-mat.softphysics.bio-phq-bio.QM
keywords starlingflocksInertialSpinModelFermi-Pasta-Ulam-TsingoucollectiveturnswavessolitonsLorentziancorrelationeffectivestiffness
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Starling flocks turn collectively in a way that looks like a wave: the change of heading sweeps across the group with little damping and a speed that grows with polarization. But new long recordings, made with panning cameras that follow flocks for up to 12.5 seconds, show that the birds' spontaneous velocity fluctuations are overdamped, with a Lorentzian correlation function and no spin-wave peaks. The paper argues that this apparent paradox disappears once the alignment interaction contains a quartic Fermi–Pasta–Ulam–Tsingou term, so the effective stiffness grows with the local phase distortion. Small spontaneous fluctuations feel the small quadratic stiffness and stay overdamped, while a strong turn activates the large quartic stiffness and propagates linearly. If the argument is right, collective turns in starling flocks are solitary waves.

What carries the argument

The central object is the effective stiffness $J_{\mathrm{eff}} = J + J_4\,\delta\varphi^2$, a named identity (equation (7) of the paper) that replaces the constant alignment coupling of the Inertial Spin Model. It carries the argument because it lets weak fluctuations and strong perturbations see different stiffnesses, thereby evading the linear-response tie between propagation and correlation. The planar phase dynamics, once noise and dissipation are dropped, become the Fermi–Pasta–Ulam–Tsingou equations, providing the soliton mechanism for the traveling turns.

What would settle it

Turn one bird by a large angle and measure the restoring torque on its neighbours: the quartic mechanism requires a cubic nonlinearity in the force as the phase distortion approaches order one, so a purely linear restoring force at all amplitudes would close the parameter window.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that the two experimental facts—underdamped linear propagation of collective turns and overdamped Lorentzian spontaneous correlation—are not contradictory once response and unperturbed correlation are governed by different effective stiffnesses. Adding the quartic term $\frac{J_4}{4v_0^4}\sum_{ij} n_{ij}(\mathbf{v}_i-\mathbf{v}_j)^4$ to the ISM Hamiltonian yields an effective stiffness $J_{\mathrm{eff}} = J + J_4\,\delta\varphi^2$, where $\delta\varphi$ is the local phase distortion. In the parameter window $J \lesssim \chi\gamma^2(L/a)^2 \ll J_4 \lesssim J/(1-\Phi)$, spontaneous fluctuations (small $\delta\varphi$) are overdamped, producing the Lorentzian $C(k,\omega)$ with half-width $\omega_L \sim k^2$, while a perturbation with $\delta\varphi \sim O(1)$ becomes underdamped and propagates with speed $c_s \sim 1/\sqrt{1-\Phi}$, matching the observed turns. Numerically, the three-dimensional ISM+FPUT model reproduces both features, and the planar phase dynamics reduces to the Fermi–Pasta–Ulam–Tsingou equations, whose solitary waves the paper proposes as the mechanism of collective turns.

Load-bearing premise

The whole resolution rests on real flocks having a hidden quartic stiffness $J_4$ that is simultaneously large enough to dominate strong turns and small enough to stay dormant under spontaneous fluctuations, yet $J_4$ is never measured directly in the paper.

Editorial extensions

If this is right

  • The new panning-camera data from nine flocks show that the Lorentzian correlation with $\omega_L \sim k^2$ is a robust feature, not an artifact of short recordings.
  • The ISM+FPUT model reproduces in simulation both the overdamped spontaneous correlation and the underdamped linear propagation of a collective turn.
  • The predicted relation $c_s \sim 1/\sqrt{1-\Phi}$ between turn speed and polarization survives the addition of the quartic term, so the new model keeps the agreement with earlier measurements.
  • The quartic term is formally similar to the speed-control mechanism used to explain scale-free speed fluctuations, suggesting a common nonlinear response to large stimuli in bird flocks.
  • The traveling waves observed in flocks are interpreted as FPUT solitons propagating on an adiabatically evolving three-dimensional interaction network.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the FPUT soliton picture is right, the propagation speed of a turn should depend on the amplitude of the heading change—a signature that can be tested on existing recordings without new experiments.
  • The stiffness-switch mechanism is generic: any system whose interaction strength grows with local misalignment will decouple small fluctuations from large perturbations, so similar 'wave propagation without spin-wave correlation' behavior might appear in other collective animal groups or active matter.
  • The paper assumes the interaction network evolves adiabatically during a turn; a sharper test would check whether the turning wave's shape and speed distort measurably as the flock's topology changes.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The paper reports new panning-camera observations of starling flocks showing that collective turns propagate linearly with little damping while the unperturbed velocity correlation function is Lorentzian with half-width scaling as omega_L ~ k^2 and no spin-wave peaks. The authors argue that the Inertial Spin Model (ISM) cannot account for this combination, because its underdamped regime implies spin-wave peaks in the unperturbed correlation. They propose adding a quartic alignment term to the pseudo-Hamiltonian, defining an effective stiffness J_eff = J + J4 (delta_phi)^2, and derive the parameter window in inequality (8): J less than or comparable to chi gamma^2 (L/a)^2 much less than J4 less than or comparable to J/(1-Phi). The claim is that spontaneous fluctuations see J and are overdamped, while strong perturbations activate J4 and become underdamped. Simulations with J=20, J4=10^5, eta=7, T=10^-6 show both a Lorentzian correlation and linear propagation after a 90-degree initiator turn. Finally, in the planar, noiseless, dissipation-free fixed-network limit the equations reduce to the Fermi-Pasta-Ulam-Tsingou (FPUT) form, and the authors suggest that the observed traveling turns are FPUT solitons.

Significance. If the proposed mechanism holds, it resolves a genuine paradox in a natural active system: the separation of the linear response to strong perturbations from the spontaneous fluctuation spectrum through a quartic stiffness. The new experimental dataset, with tracking times up to 12.5 s, is itself a valuable contribution, and the explicit connection between the model and the FPUT equation is elegant. However, the central reconciliation currently rests on a single simulation with hand-picked parameters, a parameter window validated only at k ~ 1/L, and an interpretive leap to solitons that the authors themselves state requires further work. No public code or experimental error bars accompany the analysis, and the quartic stiffness J4 is never estimated from the data. The paper is best read as a promising proof of principle rather than a validated reconciliation.

major comments (3)
  1. The overdamped side of inequality (8) is imposed only at the scale k ~ 1/L, but the Lorentzian correlation functions are measured at k = 1.5-5 m^-1 in the experiments (Figs. S1-S2) and displayed at k = 0.63, 1.26, and 1.88 in the ISM+FPUT simulation (Fig. 3b). For the linearized ISM, the spectral function of transverse fluctuations changes from a central Lorentzian to a spectrum with resolved side peaks when (J/chi)a^2 k^2 exceeds 2 gamma^2, not merely gamma^2. With the quoted simulation parameters (J=20, chi=1.25, eta=7, hence gamma=2.8 and gamma^2=7.84), the mode at k=1.88 has (J/chi)a^2 k^2 between roughly 42 (using the exact cubic-lattice eigenvalue 2J(1-cos k)/chi) and 56 (using the continuum approximation), well above 2 gamma^2=15.7; this mode should therefore show resolvable spin-wave peaks, contradicting the claimed Lorentzian correlation. The paper must specify the frequency range and the peak-resolution criterion used in the analysis, and either correct the quoted parameters or demonstrate that the full anisotropic active correlation remains Lorentzian at these wavevectors.
  2. The soliton interpretation is not supported by any calculation or simulation in the manuscript. Equation (9) is the FPUT equation only in the planar, noise-free, dissipation-free, fixed-network limit; the paper itself states that the problem 'requires to study FPUT solitons on an adiabatically evolving 3d interaction network' and to include weak noise and dissipation, but no such study is presented. The simulations show a propagating disturbance after a 90-degree turn, but no evidence is given that this disturbance is a soliton (for example, an amplitude-width-speed relation, shape persistence, or soliton collision tests) rather than a weakly damped linear wave in the J4-dominated regime. Since the title and abstract advertise 'a case for soliton propagation', this load-bearing interpretive step either needs to be substantiated or the claims should be softened to nonlinear propagating pulses.
  3. The reconciliation is not tested against the new experimental data. Inequality (8) involves J, J4, gamma, chi, L, and Phi, but the paper reports no estimates of these parameters from the starling measurements and does not show that measured polarizations and turn speeds fall inside the window. In particular, if the overdamped condition is imposed at the measured wavevectors rather than at k ~ 1/L, the combination of J4 much greater than chi gamma^2 (L/a)^2, J4 less than or comparable to J/(1-Phi), and J less than or comparable to chi gamma^2/(a^2 k_max^2) implies (1-Phi) less than or approximately 1/(k_max L)^2. For L/a ~ 10 and k_max a ~ 5 this requires Phi greater than about 0.9996, a much stronger constraint than the paper's invocation of the margin 1/(1-Phi) much greater than 1. No evidence is provided that real flocks satisfy this k-resolved version of the window, so the central claim that the FPUT term reconciles theory with experiment remains a proof of principle rather than a validated explanation.
minor comments (3)
  1. The expansion of (v_i - v_j)^2 contains the term (pi_i^2 - pi_j^2)^2/(4 v_0^2), which is subsequently dropped using J << J4; the paper should state why this term is negligible relative to the retained quartic term, since the two terms have different dependence on the phase differences.
  2. The normalization of C(k, omega) and the frequency binning are not defined; please specify the normalization convention and the effective frequency resolution so that the absence of spin-wave peaks can be assessed quantitatively.
  3. The text uses 'Lorentzian' and 'quasi-Lorentzian' interchangeably; it would help to state the fitted functional form and the criterion used to discriminate Lorentzian from spin-wave spectra (for example, the ratio of the spectral value at the expected spin-wave frequency to the central value).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the ISM+FPUT reconciliation is a parameter-window model demonstration, not a hidden reduction of output to input.

full rationale

The paper does not exhibit a circular derivation. The central mechanism is a two-stiffness Hamiltonian: a small quadratic stiffness J governs weak spontaneous fluctuations, and a large quartic stiffness J4 is activated by finite phase distortions. Inequality (8) is explicitly derived from the desired operating conditions (overdamped spontaneous correlation, dormant quartic term, underdamped strong perturbation), which is parameter-window construction rather than a hidden reduction of the output to the input; the paper labels these as 'we want/require/need' conditions. The Lorentzian correlation and linear propagation are then demonstrated by simulations in that window, a proof-of-principle rather than a prediction from fitted data. The cs ~ 1/sqrt(1-Phi) relation is explicitly stated to be 'the same prediction as in the standard ISM' and is inherited, not presented as a new consequence of the FPUT term. The numerous self-citations (e.g., [1,3,7,13-15]) supply background experimental facts and prior model development but are not load-bearing for the new reconciliation: the quartic FPUT step is introduced without appeal to a self-cited uniqueness theorem, and the FPUT soliton interpretation is explicitly flagged as 'tempting' rather than derived. Concerns about whether inequality (8) is sufficient at the experimentally probed wavevectors, or about the identifiability of J and J4 from real flocks, are correctness/validation risks, not circularity.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The central claim rests on the ISM baseline plus a new nonlinear response assumption. The quartic constant J4 and the simulation parameters J, eta, chi, T are not measured from flocks; they are chosen to satisfy the target inequalities. The soliton interpretation uses the known FPUT structure but assumes it survives noise, dissipation, and a time-evolving 3D network. No independent falsifiable prediction is made beyond the already-verified ISM speed-polarization relation.

free parameters (5)
  • J (quadratic alignment stiffness) = J=20 in ISM+FPUT simulations
    Chosen small so that spontaneous fluctuations are overdamped and C(k,omega) is Lorentzian; central to inequality (8).
  • J4 (quartic stiffness) = J4=10^5 in ISM+FPUT simulations
    Chosen large so that finite perturbations become underdamped; the ratio J4/J=5000 is tuned to make the coexistence work, not measured from flocks.
  • eta (spin friction) = eta=7 in overdamped/FPUT; eta=0.7 in underdamped ISM
    Sets the damping scale gamma=eta/(2 chi) that determines the overdamped/underdamped crossover in Eq. (4).
  • chi (generalized inertia) = chi=1.25
    Inertia parameter; along with eta sets gamma; no experimental value is used.
  • T (noise temperature) = T=10^-6
    Controls polarization through (1-Phi) ~ T/J; the chosen value keeps spontaneous fluctuations tiny and helps satisfy inequality (8).
assumptions (6)
  • domain assumption ISM equations of motion (1) describe the unperturbed flock dynamics
    The whole analysis builds on the Inertial Spin Model of refs [1,3]; if the ISM is not the right baseline, the proposed FPUT extension inherits the problem.
  • domain assumption The adjacency matrix n_ij(t) varies slowly compared with the velocity relaxation time
    Invoked before Eq. (4) to justify the quasi-equilibrium dispersion relation; the SM later argues density fluctuations do not change the spin-wave picture.
  • domain assumption Fixed speed constraint |v_i|=v0
    Used to project dynamics on transverse fluctuations and to derive the second-order equation (13).
  • ad hoc to paper Planar velocity approximation for the parameter window
    Inequalities (8) are derived in a planar scalar-phase framework (SM before Eq. (24)); the actual 3D simulations have vector fluctuations, so the quantitative window is approximate.
  • ad hoc to paper FPUT solitons survive noise, dissipation and a time-evolving 3D interaction network
    The identification of collective turns with FPUT solitons is stated as an interpretation and is not derived for the active, noisy, dissipative case.
  • ad hoc to paper Nonlinear response: individuals respond weakly to small misalignments and strongly to large ones
    This behavioral assumption motivates the quartic term; it is plausible but not directly measured.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Spin-Waves without Spin-Waves: A Case for Soliton Propagation in Starling Flocks." pith.science (2026). https://pith.science/paper/TBN6GH3Z

@misc{pith2026250519665,
  author       = {Pith},
  title        = {Pith review of: Spin-Waves without Spin-Waves: A Case for Soliton Propagation in Starling Flocks},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TBN6GH3Z}},
  note         = {Machine review of arXiv:2505.19665}
}
read the original abstract

Collective turns in starling flocks propagate linearly with negligible attenuation, indicating the existence of an underdamped sector in the dispersion relation. Beside granting linear propagation of the phase perturbations, the real part of the frequency should also yield a spin-wave form of the unperturbed correlation function. However, new high-resolution experiments on real flocks show that underdamped traveling waves coexist with an overdamped Lorentzian correlation. Theory and experiments are reconciled once we add to the dynamics a Fermi-Pasta-Ulam-Tsingou term.

Figures

Figures reproduced from arXiv: 2505.19665 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

39 extracted references · 37 canonical work pages

  1. [1]

    Attanasi, A

    A. Attanasi, A. Cavagna, L. Del Castello, I. Giardina, T. S. Grigera, A. Jeli´ c, S. Melillo, L. Parisi, O. Pohl, E. Shen, et al. , Nature Physics 10, 691 (2014)

  2. [2]

    Pomeroy and F

    H. Pomeroy and F. Heppner, The Auk , 256 (1992)

  3. [3]

    Cavagna, L

    A. Cavagna, L. Del Castello, I. Giardina, T. Grigera, A. Jelic, S. Melillo, T. Mora, L. Parisi, E. Silvestri, M. Viale, et al. , Journal of Statistical Physics 158, 601 (2015)

  4. [4]

    P. C. Hohenberg and B. I. Halperin, Reviews of Modern Physics 49, 435 (1977)

  5. [5]

    Vicsek, A

    T. Vicsek, A. Czir´ ok, E. Ben-Jacob, I. Cohen, and O. Shochet, Phys Rev Lett 75, 1226 (1995)

  6. [6]

    Dyson, Physical review 102, 1217 (1956)

    F. Dyson, Physical review 102, 1217 (1956)

  7. [7]

    T. Mora, A. M. Walczak, L. Del Castello, F. Ginelli, S. Melillo, L. Parisi, M. Viale, A. Cavagna, and I. Giar- dina, Nature Physics 12, 1153 (2016)

  8. [8]

    Bialek, A

    W. Bialek, A. Cavagna, I. Giardina, T. Mora, E. Silvestri, M. Viale, and A. M. Walczak, Proc Natl Acad Sci USA 109, 4786 (2012)

Show all 39 references
  1. [9]

    Cavagna, D

    A. Cavagna, D. Conti, I. Giardina, T. S. Grigera, S. Melillo, and M. Viale, Physical Biology 13, 065001 (2016)

  2. [10]

    Toner and Y

    J. Toner and Y. Tu, Phys Rev Lett 75, 4326 (1995)

  3. [11]

    Y. Tu, J. Toner, and M. Ulm, Physical review letters 80, 4819 (1998)

  4. [12]

    Toner and Y

    J. Toner and Y. Tu, Physical review E 58, 4828 (1998)

  5. [13]

    Cavagna, A

    A. Cavagna, A. Culla, L. Di Carlo, I. Giardina, and T. S. Grigera, Comptes Rendus. Physique 20, 319 (2019)

  6. [14]

    Cavagna, A

    A. Cavagna, A. Culla, X. Feng, I. Giardina, T. S. Grigera, W. Kion-Crosby, S. Melillo, G. Pisegna, L. Postiglione, and P. Villegas, Nature Communications13, 2315 (2022)

  7. [15]

    Cavagna, A

    A. Cavagna, A. Culla, and T. S. Grigera, Physical Re- view E 106, 054136 (2022)

  8. [16]

    Fermi, P

    E. Fermi, P. Pasta, S. Ulam, and M. Tsingou, Studies of the nonlinear problems, Tech. Rep. (Los Alamos National Laboratory (LANL), Los Alamos, NM (United States), 1955)

  9. [17]

    Gallavotti, The Fermi-Pasta-Ulam problem: a status report, Vol

    G. Gallavotti, The Fermi-Pasta-Ulam problem: a status report, Vol. 728 (Springer, 2008)

  10. [18]

    Friesecke and J

    G. Friesecke and J. A. Wattis, Communications in math- ematical physics 161, 391 (1994)

  11. [19]

    Friesecke and K

    G. Friesecke and K. Matthies, Discrete and Continuous Dynamical Systems - B 3, 105 (2003)

  12. [20]

    Chen and M

    F. Chen and M. Herrmann, Discrete and Continuous Dy- namical Systems 38, 2305 (2018)

  13. [21]

    Orlowski and K

    A. Orlowski and K. Sobczyk, Reports on Mathematical Physics 27, 59 (1989)

  14. [22]

    Purwins, H

    H.-G. Purwins, H. B¨ odeker, and S. A. and, Advances in Physics 59, 485 (2010)

  15. [23]

    Cavagna, X

    A. Cavagna, X. Feng, S. Melillo, L. Parisi, L. Postiglione, 5 and P. Villegas, IEEE Transactions on Instrumentation and Measurement 70, 1 (2021)

  16. [24]

    Attanasi, A

    A. Attanasi, A. Cavagna, L. Del Castello, I. Giardina, A. Jeli´ c, S. Melillo, L. Parisi, F. Pellacini, E. Shen, E. Sil- vestri, et al. , IEEE transactions on pattern analysis and machine intelligence 37, 2451 (2015)

  17. [25]

    Cavagna, I

    A. Cavagna, I. Giardina, and T. S. Grigera, Physics Reports 728, 1 (2018)

  18. [26]

    Cavagna, J

    A. Cavagna, J. Crist ´ ın, I. Giardina, T. S. Grigera, and M. Veca, Journal of Physics A: Mathematical and Theo- retical 57, 415002 (2024)

  19. [27]

    M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Clarendon Press, Oxford, 1987)

  20. [28]

    W. C. Swope, H. C. Andersen, P. H. Berens, and K. R. Wilson, The Journal of Chemical Physics 76, 637 (1982)

  21. [29]

    Cavagna, L

    A. Cavagna, L. Di Carlo, I. Giardina, T. S. Grigera, and G. Pisegna, Physical Review Research 3, 013210 (2021)

  22. [30]

    Cavagna, L

    A. Cavagna, L. Di Carlo, I. Giardina, T. S. Grigera, S. Melillo, L. Parisi, G. Pisegna, and M. Scandolo, Na- ture Physics 19, 1043 (2023)

  23. [31]

    Quantum Field Theory and Critical Phe- nomena,

    J. Zinn-Justin, “Quantum Field Theory and Critical Phe- nomena,” (2002), international Series of Monographs on Physics 113, 1054 pp. (2002), Fourth Edition

  24. [32]

    Cavagna, I

    A. Cavagna, I. Giardina, T. S. Grigera, A. Jelic, D. Levine, S. Ramaswamy, and M. Viale, Physical Re- view Letters 114, 218101 (2015)

  25. [33]

    Cavagna, L

    A. Cavagna, L. Di Carlo, I. Giardina, L. Grandinetti, T. S. Grigera, and G. Pisegna, Physical Review E 100, 062130 (2019)

  26. [34]

    Cavagna, L

    A. Cavagna, L. Di Carlo, I. Giardina, L. Grandinetti, T. S. Grigera, and G. Pisegna, Physical Review Letters 123, 268001 (2019)

  27. [35]

    Y. Tu, J. Toner, and M. Ulm, Phys. Rev. Lett. 80, 4819 (1998)

  28. [36]

    L. H. Ryder, Quantum field theory (Cambridge university press, 1996)

  29. [37]

    A. Z. Patashinskii and V. L. Pokrovskii, Zh. Eksp. Teor. Fiz. 64, 1445 (1973)

  30. [38]

    A. Z. Patashinskii and V. L. Pokrovskii, Fluctuation The- ory of Phase Transitions (Pergamon Press, 1979). VIDEO CAPTIONS Video S1: A real flock turning. A flock of star- lings turning over their roosting site, at Piazza dei Cin- quecento, in Rome. Images are taken at 10 frame...

  31. [39]

    Depend- ing on the strength of this friction, a crossover between an underdamped and overdamped dynamics is known to arise [33]

    The resulting theory describes the dynamics of ψ and s on a fixed network, and is given by Model G from [4] with the addition of the spin friction η. Depend- ing on the strength of this friction, a crossover between an underdamped and overdamped dynamics is known to arise [33]...

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.