REVIEW 3 major objections 3 minor 2 cited by
Mathematical Models for Fish Schooling
T0 review · 3 major / 3 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read This note reviews two stochastic models of fish schooling and reports a simulation in which a 25-fish school, using attraction, repulsion, velocity matching, obstacle reflection, and scent-gradient sensing, reaches a food source while keepi
desk verdict A clear review of the authors' own fish-schooling SDE models, but the one new simulation example is not reproducible and quietly simulates a modified model. 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 machinery is a five-component force balance in which the environment force carries all information about obstacles and food. In free space, the pairwise terms are the generalized inverse-square attraction $-\alpha r^p (x_i-x_j)/\|x_i-x_j\|^p$, the repulsion $\alpha r^q (x_i-x_j)/\|x_i-x_j\|^q$, and the distance-weighted velocity matching $-\beta(r^p/\|x_i-x_j\|^p + r^q/\|x_i-x_j\|^q)(v_i-v_j)$, with $r$ the critical distance separating attraction from repulsion. In the foraging model, the environment force adds an obstacle-reflection term $-\gamma(R^P/\|x_i-x_i^*\|^P + R^Q/\|x_i-x_i^*\|^Q)(v_i-v_i^*)$, where $x_i^*\in\partial\Omega$ is the predicted collision point and $v_i^*$ is the ref
What would settle it
Run model (3.3) with the Section 4 parameters but with several independent Brownian-motion realizations and a stated integrator and time step; if the school fails to reach the food in a majority of runs, or if any run blows up when two fish nearly coincide, then the displayed pattern is a single-trial artifact rather than a property of the model.
Extended reading notes
Core claim
The central claim is that model (3.3) reproduces collective foraging. The model writes each fish's position as an SDE with Brownian noise and its velocity as a deterministic balance of four pairwise forces—attraction, repulsion, velocity matching—plus an obstacle-reflection force and a food-scent gradient $\alpha_3\nabla X(x_i)$. The scent density $X$ is not prescribed by hand but solves $-\alpha_1\Delta X + \alpha_2 X = F$ with homogeneous Neumann boundary conditions, where $F$ is the food density. In the Section 4 example, with $N=25$ fish starting in $[0,2]\times[3.5,4]$ and parameters $\alpha=1$, $\beta=0.5$, $\gamma=1$, $\alpha_3=5.5$, $\sigma=0.001$, and obstacle rectangles $[2,2.5]\ti
Load-bearing premise
The example rests on a single random realization with hand-picked parameter values and an unspecified numerical integration scheme; if that one trajectory is not representative of model (3.3), or if the singular pairwise forces at $\|x_i-x_j\|\to 0$ break the simulation, then the claimed 'the school reaches to the food source' does not follow from the model.
Editorial extensions
If this is right
- Leaderless fish can keep a school coherent from pairwise local rules alone, since model (2.1) with the five components has no global information.
- Adding a scent-gradient force to the environment term is sufficient to route a school to a resource around obstacles in model (3.3).
- The elliptic equation for scent turns any food-density function $F$ into a smooth potential field, so the same model can handle multiple food sources or moving resources by changing $F$.
- The free-space model (2.1) is the special case of (3.3) in which obstacle and food forces vanish, so the reviewed SDE results for free space carry over to the foraging model.
Reading between the lines
- A robustness check not reported in the paper: run the same parameters with many independent noise realizations and measure the fraction of runs that reach the food and keep the school shape.
- The scent potential $X$ solves a screened Poisson equation; the same force structure could represent chemical gradients, temperature plumes, or light fields, so the model may transfer to other taxis behaviors without changing the equations.
- The singular $1/\|x_i-x_j\|^p$ forces require a numerical scheme that keeps fish separated or regularizes the potential; the paper does not state its integrator or time step, so a reproduction would need to fix those choices.
- All parameters are stated in dimensionless units (e.g., $r=0.1$, $R=0.2$ on a $7\times4$ domain); applying the model to real fish would require rescaling distances, speeds, and noise to experimental data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This short note reviews two stochastic models for fish schooling developed by the authors. Section 2 presents an SDE model in free space, Eq. (2.1), combining attraction, repulsion, velocity alignment, environmental force, and noise. Section 3 extends it to domains with obstacles and a fixed food resource by coupling the fish SDEs to an elliptic PDE for scent density, Eqs. (3.1)-(3.3); fish are attracted along the scent gradient. Section 4 gives a two-dimensional simulation of 25 fish in an obstacle domain, with a food source at (6,0.1), and Figure 1 shows snapshots at t=0,50,100,150. The caption claims the school reaches the food source while maintaining school structure. The paper cites several earlier works, mostly by the same authors, for existence and more detailed results.
Significance. If the model and the numerical example were substantiated, the paper would contribute a useful stochastic agent-based formulation of collective foraging, with the interesting feature that food information is transmitted through a scalar scent field solving an elliptic PDE rather than through direct attraction to the food. The equations are stated explicitly and the potential-function approach is transparent, and the review of the two model families is concise. However, the central demonstration in Section 4 is a single unreproducible trajectory, and the SDEs have singular interaction terms whose well-posedness is not established. The paper therefore does not, in its current form, provide a reliable computational or mathematical validation of the claimed foraging behavior.
major comments (3)
- [Section 4, Eq. (3.3) / speed-cap paragraph] The numerical example is not a simulation of the stated model (3.3). The text adds a speed-limiting projection: 'If the magnitude of vi exceeds ||v||max, our program would reset vi to a vector of magnitude ||v||max and same direction.' This reset is not part of the SDE (3.3), so Figure 1 illustrates a hybrid model. The captions and text claim the school reaches the food source 'for the model (3.3)', but no argument is given that solutions of (3.3) are unaffected by the reset, or that the reset is a negligible regularization. This is load-bearing because the plotted trajectories may owe their stability and obstacle avoidance precisely to this unmodeled projection.
- [Section 3, Eq. (3.3), interaction terms] The interaction terms in (3.3) are singular at xi = xj: they contain (xi - xj)/||xi - xj||^p and (vi - vj)/||xi - xj||^q with p=3, q=4 in the example. At xi = xj the right-hand side is undefined, and for nearby particles the force diverges as a high power of the inverse distance. The paper gives no regularization, no collision-avoidance exclusion, no invariant preventing collisions, and no well-posedness theorem for (3.3) with these coefficients. A naive numerical integrator may encounter enormous forces during near-collisions, and the plotted trajectory may depend on an implicit or ad hoc regularization. This is a core mathematical gap, not a presentation issue, for any claim that (3.3) itself produces the shown behavior.
- [Section 4, Figure 1 / numerical protocol] The example is a single random realization with no numerical scheme, time step, random seed, or precise initial configuration ('taken randomly'). Consequently the figure cannot be reproduped or checked, and a single trajectory cannot support a claim that the model produces collective foraging. With N=25 and the strong singular forces, different seeds or integrator choices may give qualitatively different outcomes. The lack of ensemble statistics, error bars, or comparison with a baseline (e.g., without scent gradient, or without velocity alignment) means the central statement 'The school reaches to the food source' is not a demonstrated property of (3.3) but an anecdote.
minor comments (3)
- [Section 3, Eq. (3.3), x_i^* and v_i^*] The obstacle-avoidance term depends on a point x_i^* on ∂Ω and a reflection vector v_i^*, but the text does not say how x_i^* is determined from xi and vi for a general obstacle domain. For the non-convex domain in Section 4 the reflection law is not unique. A precise definition would be needed to make the model reproducible.
- [Section 4, parameters] The list of parameters is given, but no rationale or sensitivity study is provided. In particular, alpha3=5.5, sigma=0.001, and the speed cap 0.8 are hand-picked. A short sensitivity discussion or at least a second realization would strengthen the presentation.
- [General] There are several typographical slips and awkward phrasings, e.g., 'mathching' in Section 3, 'The sc hool reaches' in Figure 1 caption, and 'stochastic differentiation' where 'differential' seems intended. The paper would also benefit from a statement of which results are new here versus reviewed from [8], [10], [16].
Circularity Check
No significant circularity: §4 is an illustrative simulation of a fully stated model; food-seeking is a modeling input, not a hidden prediction.
full rationale
The paper is a review of the authors' own stochastic models and presents a single illustrative simulation in §4. The model equations (3.3) are stated in full, so the simulation is self-contained rather than relying on an unstated cited theorem. The food-seeking term α3∇X(xi) is an explicit modeling assumption, and the §4 caption's observation that the school reaches the food is a consistency check of that assumption, not an independent prediction derived from something else. No parameter is fitted to data, and no quantity is predicted that was used as an input. The numerous self-citations identify the origins of the models but are not load-bearing for the example's claim; the example is presented directly with its own equations and figure. The speed reset and missing numerical details are reproducibility/correctness concerns, not circularity. Hence no circular step is identified.
Assumptions & free parameters
free parameters (10)
- alpha (attraction coefficient) =
1
- beta (velocity alignment coefficient) =
0.5
- gamma (obstacle avoidance strength) =
1
- r (attraction/repulsion critical distance) =
0.1
- R (obstacle avoidance distance) =
0.2
- alpha_3 (scent sensitivity) =
5.5
- sigma (noise intensity) =
0.001
- v_max (speed cap) =
0.8
- alpha_1 (scent diffusion) =
0.1
- alpha_2 (scent decline rate) =
0.2
assumptions (4)
- domain assumption Behavioral rules (R1)-(R3) of Camazine et al. are a valid basis for fish schooling models.
- domain assumption The scent field X satisfies the elliptic Neumann problem (3.1), and fish move up its gradient (R5).
- standard math The SDE systems (2.1) and (3.3) have unique strong solutions.
- ad hoc to paper The numerical simulation in Section 4 converges to the true solution of (3.3).
Cite this review
Pith. "Pith review of Mathematical Models for Fish Schooling." pith.science (2026). https://pith.science/paper/3NYEJTT4
@misc{pith2026250808310,
author = {Pith},
title = {Pith review of: Mathematical Models for Fish Schooling},
year = {2026},
howpublished = {\url{https://pith.science/paper/3NYEJTT4}},
note = {Machine review of arXiv:2508.08310}
}
read the original abstract
This note reviews our mathematical models for fish schooling, considered in free space, and in space with obstacle and food resource. These models are performed by stochastic differential equations or stochastic partial differential equations. We then present an example for the model in the last case.
Forward citations
Cited by 2 Pith papers
-
Towards High-Fidelity and Controllable Bioacoustic Generation via Enhanced Diffusion Learning
BirdDiff combines a multi-band enhancement stage with a DiffWave-based diffusion generator and multimodal conditioning, reporting substantially better bird-call synthesis metrics than DiffWave on a 12-species propriet...
-
Deep Learning for Automated Identification of Vietnamese Timber Species: A Tool for Ecological Monitoring and Conservation
ShuffleNetV2 achieves 99.29% accuracy and 99.35% F1-score in classifying ten Vietnamese timber species from images, outperforming heavier models in speed-accuracy trade-off.
Reference graph
Works this paper leans on
-
[1]
Swarming behavior of a fish school consisting of a large number of individuals often surprises us
Introduction. Swarming behavior of a fish school consisting of a large number of individuals often surprises us. They swim co herently match- ing their velocity without collision and maintaining a cons tant scale of school, even though they have only moderate ability of infor mation pro- cessing and of execution of programming. Several mathematical models ...
arXiv 2010
-
[2]
In [ 16], we introduced five components constructing a SDE model for fish schooling in free space
Fish schooling model based on five components in free spaces. In [ 16], we introduced five components constructing a SDE model for fish schooling in free space. These components are based on the fo llowing three local rules of Camazine-Deneubourg-Franks-Sneyd-Theraulaz-Bonabeau ([2]): (R1) The school has no leaders and each fish follows the same be havioral r...
-
[3]
When two particles i and j are far from each other, both would move toward each other
Attraction force. When two particles i and j are far from each other, both would move toward each other. In our model, this force is a generalization of the inverse-square law of universal grav itation: −αrp(xi − xj) ∥xi − xj∥p , where 1 < p < ∞ and r > 0 are constants, α a coefficient of attraction among individuals
-
[4]
When two particles i and j are close enough, both would move far from each other
Repulsion force. When two particles i and j are close enough, both would move far from each other. This force is a generalizatio n of the Van der Waals forces: αrq(xi − xj) ∥xi − xj∥q , where 1 < p < q < ∞ is a fixed exponent
-
[5]
Alignment or velocity matching. The velocity matching of the particle i to the particle j also has a similar weight depending on the distance ∥xi − xj∥: −β ( rp ∥xi − xj∥p + rq ∥xi − xj∥q ) (vi − vj), where β is a coefficient of velocity matching among individuals
-
[6]
The individual i react to the environment a force Fi(xi, vi)
Reaction to the environment. The individual i react to the environment a force Fi(xi, vi). FISH SCHOOLING 3
-
[7]
All particles are subject to random factors or nois e
Noise. All particles are subject to random factors or nois e. These five components form our SDE model of the form (2.1) dxi(t) = vidt + σidwi(t), i = 1, 2, . . . N, dvi(t) = { − α N∑ j=1,j̸=i ( rp ∥xi−xj ∥p − rq ∥xi−xj ∥q ) (xi − xj) −β N∑ j=1,j̸=i ( rp ∥xi−xj∥p + rq ∥xi−xj∥q ) (vi − vj) +Fi(xi, vi) } dt, i = 1, 2, . ...
-
[8]
In [ 10], we introduced a fish schooling model in spaces with ob- stacle
Fish schooling model in spaces with obstacle and food re- source. In [ 10], we introduced a fish schooling model in spaces with ob- stacle. For this model, in addition to the three local rules ( a)–(c) in the previous section, we newly presented a local rule of obstacl e avoidance for individual fish: (R4) Each fish executes an action for avoiding obstacle a...
Show all 34 references
-
[9]
In this section, we give an example for the model ( 3.3) in two-dimensional space
An example. In this section, we give an example for the model ( 3.3) in two-dimensional space. The reader can find more results on (3.3) in [ 8]. Put a food resource at a small circle of radius 0 .04 and center (6 , 0.1). More precisely, the function F of food resource in ( 3.1...
-
[10]
Aoki, A simulation study on the schooling mechanism in fish, Bull
I. Aoki, A simulation study on the schooling mechanism in fish, Bull. Japan. Soc. Sci. Fish, 48 (1982), 1081–1088
1982
-
[11]
Camazine, J
S. Camazine, J. L. Deneubourg, N. R. Franks, J. Sneyd, G. T heraulaz, E. Bonabeau, Self-organization in Biological System , Princeton University Press, 2001
2001
-
[12]
Cucker, E
F. Cucker, E. Mordecki, Flocking in noisy environments, J. Math. Pures Appl. , 89 (2008), 278–296
2008
-
[13]
Cucker, S
F. Cucker, S. Smale, On the mathematics of emergence, Japan. J. Math. , 2 (2007), 197–227. FISH SCHOOLING 7
2007
-
[14]
Cucker, S
F. Cucker, S. Smale, Emergence behavior in flocks, IEEE Trans. Automat. Control 52 (2007), 852–862
2007
-
[15]
M. R. D’Orsogna, Y. L. Chuang, A. L. Bertozzi, L. S. Chayes , Self-propelled particles with soft-core interactions: Patterns, stability, and col lapse, Phys. Rev. Lett. , 104302 (2006), 1–4
2006
-
[16]
A. Huth, C. Wissel, The simulation of the movement of fish s chool, J. Theor. Biol. , 156 (1992), 365–385
1992
-
[17]
T. V. Ta . , L. T. H. Nguyen, A stochastic differential equation model fo r foraging behavior of fish schools, Phys. Biol. , 15 (2018), 036007. arXiv:1509.00063
2018 arXiv
-
[18]
L. T. H. Nguyen, T. V. Ta . , A. Yagi, Quantitative investigations for ODE model describing fish schooling, Sci. Math. Jpn. , 77 (2014), 403–413
2014
-
[19]
L. T. H. Nguyen, T. V. Ta . , A. Yagi, Obstacle avoiding patterns and cohesiveness of fish school, J. Theor. Biol. , 406 (2016), 116–123
2016
-
[20]
Oboshi, S
T. Oboshi, S. Kato, A. Mutoh, H. Itoh, Collective or scat tering: evolving schooling behaviors to escape from predator, Artif. Life , VIII (2002), 386–389
2002
-
[21]
Olfati-Saber, Flocking for multi-agent dynamic sys tems: Algorithms and Theory, IEEE Trans
R. Olfati-Saber, Flocking for multi-agent dynamic sys tems: Algorithms and Theory, IEEE Trans. Automat. Control , 51 (2006), 401–420
2006
-
[22]
C. W. Reynolds, Flocks, herds, and schools: a distribut ed behavioral model, Computer Graphics, 21 (1987), 25–34
1987
-
[23]
Shklarsh, G
A. Shklarsh, G. Ariel, E. Schneidman, E. Ben-Jacob, Sma rt swarms of bacteria- inspired agents with performance adaptable interactions, PLoS Comput. Biol. , 7 (2011), e1002177
2011
-
[24]
T. V. Ta. , L. T. H. Nguyen, A. Yagi, Flocking and non-flocking behavior in a stochastic Cucker-Smale system, Anal. Appl. , 12 (2014), 63–73
2014
-
[25]
Uchitane, T
T. Uchitane, T. V. Ta . , A. Yagi, An ordinary differential equation model for fish schooling, Sci. Math. Jpn. , 75 (2012), 339–350
2012
-
[26]
Vicsek, A
T. Vicsek, A. Czirok, E. Ben-Jacob, I. Cohen, O. Shochet , Novel type of phase tran- sition in a system of self-driven particles, Phys. Rev. Lett. , 75 (1995), 1226–1229
1995
-
[27]
Zienkiewicz, D
A. Zienkiewicz, D. A. W. Barton, M. Porfiri, M. di Bernard o, Data-driven stochastic modelling of zebrafish locomotion, J. Math. Biol. , 71 (2015), 1081–1105
2015
-
[28]
T. V. Ta, L. T. H. Nguyen, A. Yagi, A sustainability condi tion for stochastic forest model, Commun. Pure Appl. Anal. , 16 (2017), 699–718. doi:10.3934/cpaa.2017034
2017 doi
-
[29]
A. D. Hartono, L. T. H. Nguyen, T. V. Ta, A stochastic diffe rential equation model for predator-avoidance fish schooling, Math. Biosci. , 367 (2024), 109112
2024
-
[30]
T. V. Ta, Dynamics of species in a non-autonomous Lotka- Volterra system, Acta Math. Acad. Paedagog. Nyh´ azi. (N.S.), 25 (2009), 45–54
2009
-
[31]
Math-for-Industry
A. D. Hartono, T. V. Ta, L. T. H. Nguyen, A geometrical str ucture for predator- avoidance fish schooling, Proc. Forum “Math-for-Industry” 2022 - Mathematics of Public Health and Sustainability , 75–89, ISSN:2188-1200 (2023)
2022
-
[32]
T. V. Ta, Survival of three species in a non-autonomous L otka–Volterra system, J. Math. Anal. Appl. , 362 (2010), 427–437
2010
-
[33]
Math-for-Industry
L. T. H. Nguyen, T. V. Ta, A. Yagi, A brief review of some sw arming models using stochastic differential equations, In: Cheng, J., Dinghua, X., Saeki, O., Shirai, T. (eds) Proc. Forum “Math-for-Industry” 2018. Mathematics for Ind ustry, vol 35, Springer, Singapore, 2021
2018
-
[34]
J. Qi, T. Casse, M. Harada, L. T. H. Nguyen, T. V. Ta, Quant ifying fish school frag- mentation under predation using stochastic differential eq uations, arXiv:2508.00953 (2025). 8 Linh Thi Hoai Nguyen Department of Immunobiology and Neuroscience Medical Institute of Bioregulati...
2025 arXiv
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.