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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- 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.
- 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.
- 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)
- 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.
- 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.
- 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
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
free parameters (5)
- J (quadratic alignment stiffness) =
J=20 in ISM+FPUT simulations
- J4 (quartic stiffness) =
J4=10^5 in ISM+FPUT simulations
- eta (spin friction) =
eta=7 in overdamped/FPUT; eta=0.7 in underdamped ISM
- chi (generalized inertia) =
chi=1.25
- T (noise temperature) =
T=10^-6
assumptions (6)
- domain assumption ISM equations of motion (1) describe the unperturbed flock dynamics
- domain assumption The adjacency matrix n_ij(t) varies slowly compared with the velocity relaxation time
- domain assumption Fixed speed constraint |v_i|=v0
- ad hoc to paper Planar velocity approximation for the parameter window
- ad hoc to paper FPUT solitons survive noise, dissipation and a time-evolving 3D interaction network
- ad hoc to paper Nonlinear response: individuals respond weakly to small misalignments and strongly to large ones
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
Reference graph
Works this paper leans on
-
[1]
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)
work page 2014
- [2]
-
[3]
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)
work page 2015
-
[4]
P. C. Hohenberg and B. I. Halperin, Reviews of Modern Physics 49, 435 (1977)
1977
- [5]
- [6]
-
[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)
work page 2016
- [8]
Show all 39 references
-
[9]
Cavagna, D
A. Cavagna, D. Conti, I. Giardina, T. S. Grigera, S. Melillo, and M. Viale, Physical Biology 13, 065001 (2016)
2016
-
[10]
Toner and Y
J. Toner and Y. Tu, Phys Rev Lett 75, 4326 (1995)
1995
-
[11]
Y. Tu, J. Toner, and M. Ulm, Physical review letters 80, 4819 (1998)
1998
-
[12]
Toner and Y
J. Toner and Y. Tu, Physical review E 58, 4828 (1998)
1998
-
[13]
Cavagna, A
A. Cavagna, A. Culla, L. Di Carlo, I. Giardina, and T. S. Grigera, Comptes Rendus. Physique 20, 319 (2019)
2019
-
[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)
2022
-
[15]
Cavagna, A
A. Cavagna, A. Culla, and T. S. Grigera, Physical Re- view E 106, 054136 (2022)
2022
-
[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)
1955
-
[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)
2008
-
[18]
Friesecke and J
G. Friesecke and J. A. Wattis, Communications in math- ematical physics 161, 391 (1994)
1994
-
[19]
Friesecke and K
G. Friesecke and K. Matthies, Discrete and Continuous Dynamical Systems - B 3, 105 (2003)
2003
-
[20]
Chen and M
F. Chen and M. Herrmann, Discrete and Continuous Dy- namical Systems 38, 2305 (2018)
2018
-
[21]
Orlowski and K
A. Orlowski and K. Sobczyk, Reports on Mathematical Physics 27, 59 (1989)
1989
-
[22]
Purwins, H
H.-G. Purwins, H. B¨ odeker, and S. A. and, Advances in Physics 59, 485 (2010)
2010
-
[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)
2021
-
[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)
2015
-
[25]
Cavagna, I
A. Cavagna, I. Giardina, and T. S. Grigera, Physics Reports 728, 1 (2018)
2018
-
[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)
2024
-
[27]
M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Clarendon Press, Oxford, 1987)
1987
-
[28]
W. C. Swope, H. C. Andersen, P. H. Berens, and K. R. Wilson, The Journal of Chemical Physics 76, 637 (1982)
1982
-
[29]
Cavagna, L
A. Cavagna, L. Di Carlo, I. Giardina, T. S. Grigera, and G. Pisegna, Physical Review Research 3, 013210 (2021)
2021
-
[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)
2023
-
[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
2002
-
[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)
2015
-
[33]
Cavagna, L
A. Cavagna, L. Di Carlo, I. Giardina, L. Grandinetti, T. S. Grigera, and G. Pisegna, Physical Review E 100, 062130 (2019)
2019
-
[34]
Cavagna, L
A. Cavagna, L. Di Carlo, I. Giardina, L. Grandinetti, T. S. Grigera, and G. Pisegna, Physical Review Letters 123, 268001 (2019)
2019
-
[35]
Y. Tu, J. Toner, and M. Ulm, Phys. Rev. Lett. 80, 4819 (1998)
1998
-
[36]
L. H. Ryder, Quantum field theory (Cambridge university press, 1996)
1996
-
[37]
A. Z. Patashinskii and V. L. Pokrovskii, Zh. Eksp. Teor. Fiz. 64, 1445 (1973)
1973
-
[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...
1979
-
[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]...
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.