pith. sign in

arxiv: 2604.20469 · v2 · submitted 2026-04-22 · 🧮 math.AP · q-bio.PE

Short-wave signal versus indirect prey-taxis

Pith reviewed 2026-05-09 23:42 UTC · model grok-4.3

classification 🧮 math.AP q-bio.PE
keywords short-wave asymptoticsprey-taxisPatlak-Keller-Segelpredator-preyexternal signalstabilitycross-diffusionKapitza averaging
0
0 comments X

The pith

A complete asymptotic expansion describes short-wave solutions of predator-prey PDEs with indirect taxis driven by an external signal.

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

The paper builds a full asymptotic series for solutions to quasi-linear reaction-diffusion systems that model predators moving in response to a signal produced by prey. That signal is itself modulated by a rapidly varying external field independent of the population densities. The construction holds in any spatial dimension and for arbitrary inter- and intra-specific reaction kinetics. Readers care because the same expansion supplies explicit criteria for when the external signal stabilizes or destabilizes the underlying species equilibrium.

Core claim

The authors derive a complete asymptotic expansion for short-wave solutions of quasi-linear second-order PDE systems with cross-diffusion of Patlak-Keller-Segel type, where predators respond to a driving signal produced by prey whose production depends on an independent external field. They then apply the expansion to analyze stability or instability of equilibria induced by the signal, following the averaging ideas of Kapitza's inverted-pendulum theory.

What carries the argument

The complete short-wave asymptotic expansion of the solution, which reduces the original system to effective slow-scale equations whose stability can be read off directly.

Load-bearing premise

The external signal must be short-wave and the production of the driving signal must depend on an external field whose intensity is independent of the predator-prey state.

What would settle it

Compute the first correction term in the expansion for a concrete periodic external signal, then compare the resulting slow-scale stability threshold against direct numerical integration of the full PDE system as the wavelength tends to zero.

Figures

Figures reproduced from arXiv: 2604.20469 by Andrey Morgulis, Karrar Malal.

Figure 1
Figure 1. Figure 1: The sketch of a neutral set in the case of N = n = 1, where the parameter is A , and the wave number is k. often get relevant for studying the pattern formation in the dynamics of cells or populations, e.g. [7]-[27]. 2.5.2 Normal modes Let the external signal undergoes no slow modulation – that is, equalities (77) holds true. Then the linearized system (81)-(83) possesses the translational invariance, cf. … view at source ↗
read the original abstract

We address a short-wave asymptotic for one class of quasi-linear second-order PDE systems involving the cross-diffusion described by the so-called Patlak-Keller-Segel law. It is common to employ these equations for modeling the predator-prey community with the prey-taxis that means the interactions of two species of particles or cells or anything else through which the species called "predators" is capable of moving directionally while searching for the other species called "prey." However, we suppose the predators to be sensitive not to the prey density but to a driving signal produced by the prey. Additionally, the production of the driving signal is assumed to be sensitive to the intensity of an external field, which is independent from the community state. This is what we call the external signal. It can be due to the spatiotemporal inhomogeneity of the environment arising from natural or artificial reasons. We assume that the external signal takes a general short-wave form and construct a complete asymptotic expansion for the short-wave solutions with no restrictions on the spatial dimension or kinetics of inter- or intra-specific reactions. Further, we apply the short wave asymptotic to studying the stability or instability induced by the external signal following Kapitza's theory for the upside-down pendulum. Applying the general results to some special classes of external signals, we get examples of suppressing the taxical transport, examples of robustness of the species equilibrium to the signal up to a very strong stabilization or, oppositely, destabilization and somewhat like blurring the borderline in the parametric space between the areas of stability and instability of this equilibrium. These results contribute to filling the gap in the literature, since the theory and techniques for the asymptotic integration of systems described above represent a weakly charted area.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript constructs a complete (all-order) short-wave asymptotic expansion for solutions of a quasi-linear cross-diffusion system of Patlak-Keller-Segel type that models indirect prey-taxis driven by a general external short-wave signal independent of the community state. The derivation uses a two-scale ansatz, collects powers of the large wave number, and solves the resulting transport equations order by order. The expansion is applied, via Kapitza averaging, to analyze stability or instability of equilibria for special classes of signals, yielding examples of suppression of taxical transport, strong stabilization, destabilization, and blurring of stability boundaries, with no a-priori restrictions on spatial dimension or reaction kinetics.

Significance. If the claimed construction is valid, the work supplies a general tool for high-frequency asymptotics in cross-diffusion systems that is largely absent from the literature. The independence of the external field from the densities is used to close the hierarchy at every order, which is a structural strength. The subsequent stability analysis illustrates concrete ecological implications (signal-induced control of taxis), and the absence of dimensional or kinetic restrictions, if substantiated, would broaden applicability beyond the usual one-dimensional or linear-kinetics settings.

major comments (2)
  1. [Derivation of the expansion] The central derivation (substitution of the two-scale ansatz and order-by-order solution of the transport equations) is asserted to close without uncontrolled source terms for arbitrary kinetics, but the quasi-linear cross-diffusion coefficients depend on the densities that are themselves expanded; explicit verification that the resulting linear problems at each order remain solvable (or that secular terms are absent) for general nonlinear reaction terms is not visible and is load-bearing for the 'no restrictions on kinetics' claim.
  2. [Stability analysis] The stability conclusions rest on applying Kapitza averaging to the effective system obtained after the short-wave expansion; it is unclear whether the averaged coefficients remain well-defined and the stability thresholds uniform when the external signal is only assumed to be of general short-wave form rather than the special classes treated in the examples.
minor comments (3)
  1. [Abstract] The abstract uses the informal phrase 'somewhat like blurring the borderline in the parametric space'; replace with a precise statement of how the stability region in parameter space is modified.
  2. [Notation and setup] Notation for the two-scale ansatz and the large wave-number parameter should be introduced once and used consistently; currently the transition from the original PDE system to the expanded equations is abrupt.
  3. [Introduction] Add references to prior short-wave or high-frequency asymptotic analyses of Keller-Segel or cross-diffusion systems to situate the contribution.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, positive assessment, and recommendation of minor revision. We address the two major comments point by point below, indicating the revisions that will be incorporated.

read point-by-point responses
  1. Referee: The central derivation (substitution of the two-scale ansatz and order-by-order solution of the transport equations) is asserted to close without uncontrolled source terms for arbitrary kinetics, but the quasi-linear cross-diffusion coefficients depend on the densities that are themselves expanded; explicit verification that the resulting linear problems at each order remain solvable (or that secular terms are absent) for general nonlinear reaction terms is not visible and is load-bearing for the 'no restrictions on kinetics' claim.

    Authors: We thank the referee for this observation. In the derivation, the external signal enters the quasi-linear coefficients independently of the densities. Substituting the two-scale ansatz therefore produces, at each order, a linear transport problem in the fast variable whose coefficients depend on the known signal and on lower-order terms already determined. Solvability follows from periodicity in the fast scale and the fact that the right-hand side, after averaging, satisfies the necessary compatibility condition; the nonlinear reaction terms contribute only to the slow-scale modulation equations and do not generate secular growth. We will insert a short clarifying paragraph immediately after the recursive construction to make this structure explicit for arbitrary smooth kinetics. revision: yes

  2. Referee: The stability conclusions rest on applying Kapitza averaging to the effective system obtained after the short-wave expansion; it is unclear whether the averaged coefficients remain well-defined and the stability thresholds uniform when the external signal is only assumed to be of general short-wave form rather than the special classes treated in the examples.

    Authors: The referee is correct that the explicit stability thresholds are computed for concrete signal classes. However, the Kapitza averaging step itself is performed on the general effective system obtained from the short-wave expansion. For any short-wave signal satisfying the standing regularity and periodicity assumptions, the averaged coefficients are obtained by integration over the fast period and are therefore well-defined; the resulting stability criteria are expressed directly in terms of these averages and are uniform over the admissible class. The special signals serve only to obtain closed-form expressions. We will add a sentence in the stability section stating that the averaging procedure and the uniformity of the thresholds hold for the general short-wave signals treated in the expansion. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained asymptotic construction

full rationale

The central result is the explicit construction of a complete short-wave asymptotic expansion obtained by substituting a two-scale ansatz into the quasi-linear cross-diffusion PDE system, collecting powers of the large wave number, and solving the resulting transport equations order by order. The assumption that the external signal is independent of the community state closes the hierarchy at each order without introducing uncontrolled terms or fitted quantities. No step reduces by construction to its own inputs, no self-citation is load-bearing for the derivation, and no ansatz is smuggled via prior work. The method is a standard multiple-scale analysis whose validity rests on the PDEs themselves rather than on any target stability result.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard short-wave asymptotic techniques for PDEs together with modeling assumptions that the external field is independent and short-wave. No free parameters, invented entities, or ad-hoc constants are mentioned.

axioms (2)
  • standard math Standard short-wave asymptotic expansion techniques apply to the given class of quasi-linear second-order PDE systems
    Invoked to construct the complete expansion without dimensional or kinetic restrictions.
  • domain assumption The external signal is of general short-wave form and independent of community state
    Required to apply the asymptotics and perform the Kapitza-style stability analysis.

pith-pipeline@v0.9.0 · 5610 in / 1364 out tokens · 65772 ms · 2026-05-09T23:42:49.583787+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

55 extracted references · 16 canonical work pages

  1. [1]

    Kareiva, P., & Odell, G. (1987). Swarms of predators exhibit ‘preytaxis’ if individual predators use area-restricted search. The American Naturalist, 130(2), 233-270

  2. [2]

    Amann, H. (1993). Nonhomogeneous linear and quasilinear elliptic and parabolic boundary value problems. In Function spaces, differential operators and nonlinear analysis (pp. 9-126). Springer Fachmedien Wiesbaden

  3. [3]

    (1994) From the dynamics of population autowaves generated by living cells to neuroinformatics

    Ivanitskii, G.R., Medvinskii, A.B., Tsyganov, M.A. (1994) From the dynamics of population autowaves generated by living cells to neuroinformatics. Physics-Uspekhi, 37(10): 961–989

  4. [4]

    Horstmann, D. (2004). From 1970 until present: the Keller-Segel model in chemotaxis and its consequences. II, Jahresber. Deutsch. Math.-Verein., 106, 51-69

  5. [5]

    & Painter, K.J

    Hillen, T. & Painter, K.J. (2009). A user’s guide to PDE models for chemotaxis. Journal of mathematical biology, 58(1-2), 183-217

  6. [6]

    (2015) Toward a mathematical theory of Keller- Segel models of pattern formation in biological tissues

    Bellomo N.,Bellouquid A., Tao Y., Winkler M. (2015) Toward a mathematical theory of Keller- Segel models of pattern formation in biological tissues. Math. Models Methods Appl. Sci. 25(09): 1663-1763

  7. [7]

    & Karev, G.P

    Berezovskaya, F.S. & Karev, G.P. (1999). Bifurcations of travelling waves in population taxis models. Physics-Uspekhi, 42(9), 917-929

  8. [8]

    (2000) Slow taxis in a predator-prey model

    Govorukhin V., Morgulis A., Tyutyunov Y. (2000) Slow taxis in a predator-prey model. Doklady Mathematics 61(3): 420-422

  9. [9]

    (2001) Directed movement of predators and the emergence of density-dependence in predator-prey models

    Arditi R., Tyutyunov Y., Morgulis A., Govorukhin V., Senina I. (2001) Directed movement of predators and the emergence of density-dependence in predator-prey models. Theoretical Population Biology 59(3): 207-221

  10. [10]

    Sapoukhina, N., Tyutyunov, Y., & Arditi, R. (2003). The role of prey taxis in biological control: a spatial theoretical model. The American Naturalist, 162(1), 61-76

  11. [11]

    (2007) Chemotaxis-induced spatio-temporal heterogeneity in multi-species host-parasitoid systems

    Pearce I.G., Chaplain M.A.J., Schofield P.G., Anderson A.R.A, Hubbard S.F. (2007) Chemotaxis-induced spatio-temporal heterogeneity in multi-species host-parasitoid systems. J. Math. Biol. 55(3): 365-388

  12. [12]

    Chakraborty, A., Singh, M., Lucy, D., & Ridland, P. (2007). Predator–prey model with prey- taxis and diffusion. Mathematical and computer modelling, 46(3-4), 482-498

  13. [13]

    M., Hillen, T., & Lewis, M

    Lee, J. M., Hillen, T., & Lewis, M. A. (2009). Pattern formation in prey-taxis systems. Journal of biological dynamics, 3(6), 551-573

  14. [14]

    J., & Hillen, T

    Painter, K. J., & Hillen, T. (2011). Spatio-temporal chaos in a chemotaxis model. Physica D: Nonlinear Phenomena, 240(4-5), 363-375

  15. [15]

    Banerjee, M., & Petrovskii, S. (2011). Self-organised spatial patterns and chaos in a ratio- dependent predator–prey system. Theoretical Ecology, 4, 37-53

  16. [16]

    (2016) Predator-prey model with diffusion and indirect prey-taxis

    Tello I.J., Wrzosek D. (2016) Predator-prey model with diffusion and indirect prey-taxis. Math. Models Methods Appl. Sci. 26(11): 2129-2162. A. Morgulis & K. Malal.p. 28 of 30

  17. [17]

    (2017) Prey-taxis destabilizes homogeneous stationary state in spatial Gause-Kolmogorov-type model for predator-prey system

    Tyutyunov Y., Titova L., Senina I. (2017) Prey-taxis destabilizes homogeneous stationary state in spatial Gause-Kolmogorov-type model for predator-prey system. Ecological Complexity; 31: 170-180

  18. [18]

    (2017) Time-periodic and stable patterns of a two-competing- species Keller-Segel chemotaxis model: Effect of cellular growth

    Wang Q., Yang J., Zhang L. (2017) Time-periodic and stable patterns of a two-competing- species Keller-Segel chemotaxis model: Effect of cellular growth. Discrete & Continuous Dy- namical Systems - B. ; 22(9):3547-3574

  19. [19]

    (2017) Prey-taxis destabilizes homogeneous stationary state in spatial Gause-Kolmogorov-type model for predator-prey system

    Tyutyunov Y., Titova L., Senina I. (2017) Prey-taxis destabilizes homogeneous stationary state in spatial Gause-Kolmogorov-type model for predator-prey system. Ecological Complexity; 31: 170-180. Mathematics, 14(7), 1165. https://doi.org/10.3390/math14071165

  20. [20]

    V., Sen, D., Titova, L

    Tyutyunov, Y. V., Sen, D., Titova, L. I., & Banerjee, M. (2019). Predator over- comes the Allee effect due to indirect prey–taxis. Ecological complexity, 39, 100772. https://doi.org/10.1016/j.ecocom.2019.100772

  21. [21]

    The Dynamics of a Predator–Prey Model with Diffusion and Indirect Prey- Taxis

    Wang, J., Wang, M. The Dynamics of a Predator–Prey Model with Diffusion and Indirect Prey- Taxis. J. Dyn. Diff. Equat., 32, 1291–1310 (2020). https://doi.org/10.1007/s10884-019-09778-7

  22. [22]

    R., Petrovskii, S., & Banerjee, M

    Chowdhury, P. R., Petrovskii, S., & Banerjee, M. (2021). Oscillations and pattern formation in a slow–fast prey–predator system. Bulletin of mathematical biology, 83, 1-41

  23. [23]

    R., Petrovskii, S., & Banerjee, M

    Chowdhury, P. R., Petrovskii, S., & Banerjee, M. (2022). Effect of slow–fast time scale on transient dynamics in a realistic prey-predator system. Mathematics, 10(5), 699

  24. [24]

    Li, S. (2023). Positive steady-state solutions for a class of prey-predator systems with indirect prey-taxis. SIAM Journal on Mathematical Analysis, 55(6), 6342-6374. https://doi.org/10.1137/22M1529518

  25. [25]

    Roy Chowdhury, P., Banerjee, M., & Petrovskii, S. (2024). A two-timescale model of plank- ton–oxygen dynamics predicts formation of oxygen minimum zones and global anoxia. Journal of Mathematical Biology, 89(1), 8

  26. [26]

    Mu, C., Tao, W., & Wang, Z. A. (2024). Global dynamics and spatiotemporal heterogeneity of a preytaxis model with prey-induced acceleration. European Journal of Applied Mathematics, 1-33

  27. [27]

    (2024) Taxis-Driven Pattern Formation in Tri-Trophic Food Chain Model with Omnivory

    Giricheva, E. (2024) Taxis-Driven Pattern Formation in Tri-Trophic Food Chain Model with Omnivory. Mathematics, 12, 290

  28. [28]

    Wu, S., Wang, Y., & Geng, D. (2025). Dynamic and pattern formation of a diffusive preda- tor–prey model with indirect prey-taxis and indirect predator-taxis. Nonlinear Analysis: Real World Applications, 84, 104299. https://doi.org/10.1016/j.nonrwa.2024.104299

  29. [29]

    Xie, Z., & Li, Y. (2025). Classical solutions to a pursuit-evasion model with prey-taxis and indirect predator-taxis. Mathematical Methods in the Applied Sciences, 48(4), 4117-4143. https://doi.org/10.1002/mma.10536

  30. [30]

    Ignacio Tello, J., & Wrzosek, D. (2026). From indirect to direct taxis by fast re- action limit. Mathematical Models and Methods in Applied Sciences, 36(01), 173- 204.https://doi.org/10.1142/S0218202526500041 A. Morgulis & K. Malal.p. 29 of 30

  31. [31]

    Poggiale, J.-C., & Gauduchon, M. (2026). A General Class of Growth Models in a Patchy Environment Exhibiting Enhanced Production

  32. [32]

    A brief introduction to homogenization and miscellaneous applications

    Allaire, G. A brief introduction to homogenization and miscellaneous applications. IESAIM: Proceedings. In ESAIM: Proceedings (Vol. 37, pp. 1-49). EDP Sciences; 2012

  33. [33]

    Mechanics

    Landau L.D., Lifshitz E.M. Mechanics. 3d edition. Elsivier, 1982. 224 p

  34. [34]

    M., Simonenko, I

    Zen’kovskaya, S. M., Simonenko, I. B. (1966). Effect of high frequency vibration on convection initiation. Fluid Dynamics, 1(5), 35-37

  35. [35]

    Zen’kovskaya, S. M. (1968). Study of convection in a liquid layer with vibrating forces. Fluid Dynamics, 3(1), 35-36

  36. [36]

    Lyubimov; Tatiana P

    Dmitri V. Lyubimov; Tatiana P. Lyubimova; Anatoli A. Tcherepanov; Bernard H. Roux. Vi- bration influence on fluid interfaces. Comptes Rendus. M´ ecanique, Microgravity / La micrope- santeur, Volume 332 (2004) no. 5-6, pp. 467-472. doi: 10.1016/j.crme.2004.01.013

  37. [37]

    V. I. Yudovich, Vibration dynamics of systems with constraints, Dokl. Math., 42:6 (1997), 322–325

  38. [38]

    Yudovich, The dynamics of a particle on a smooth vibrating surface, Jour- nal of Applied Mathematics and Mechanics, Volume 62, Issue 6, 1998,P

    V.I. Yudovich, The dynamics of a particle on a smooth vibrating surface, Jour- nal of Applied Mathematics and Mechanics, Volume 62, Issue 6, 1998,P. 893-900, https://doi.org/10.1016/S0021-8928(98)00114-2

  39. [39]

    A., Morgulis, A

    Vladimirov, V. A., Morgulis, A. B. (2014). Relative equilibria in the Bjerknes problem. Siberian Mathematical Journal, 55(1), 35-48

  40. [40]

    Babenko, A.V., Polyakova, O.R., Tovstik, T.P. (2023). Conceptual Generalizations of the Kapitsa Problem. In: Altenbach, H., Irschik, H., Porubov, A.V. (eds) Progress in Continuum Mechanics. Advanced Structured Materials, vol 196. Springer, Cham. https://doi.org/10.1007/978-3-031-43736-6-4

  41. [41]

    Prey-Taxis vs

    Morgulis, A.; Malal, K.H. Prey-Taxis vs. An External Signal: Short-Wave Asymptotic and Stability Analysis. Mathematics 2025, 13, 261. https://doi.org/10.3390/math13020261

  42. [42]

    Prey-Taxis Vs A Shortwave External Signal In Multiple Dimensions

    Morgulis, A., Malal, K. Prey-Taxis Vs A Shortwave External Signal In Multiple Dimensions. J Math Sci (2025). https://doi.org/10.1007/s10958-025-08063-x

  43. [43]

    Morgulis, A., & Ilin, K. (2020). Indirect taxis on a fluctuating environment. Mathematics, 8(11), 2052

  44. [44]

    (1953) A strong maximum principle for parabolic equations

    Nirenberg, L. (1953) A strong maximum principle for parabolic equations. Comm. Pure Appl. Math. 6(2): 167-177

  45. [45]

    Second order equations of elliptic and parabolic type

    Landis, E.M. Second order equations of elliptic and parabolic type. American Math. Soc., 1997. 203 p

  46. [46]

    Yudovich, V. I. The linearization method in hydrodynamical stability theory American Math- ematical Soc., 1989, 170 p

  47. [47]

    Local bifurcations, center manifolds, and normal forms in infinite- dimensional dynamical systems

    Haragus, M., & Iooss, G. Local bifurcations, center manifolds, and normal forms in infinite- dimensional dynamical systems. Springer Science & Business Media, 2010. A. Morgulis & K. Malal.p. 30 of 30

  48. [48]

    and Titova LI (2021) Ratio-Dependence in Predator-Prey Systems as an Edge and Basic Minimal Model of Predator Interference

    Tyutyunov Yu.V. and Titova LI (2021) Ratio-Dependence in Predator-Prey Systems as an Edge and Basic Minimal Model of Predator Interference. Front. Ecol. Evol. 9:725041. doi: 10.3389/fevo.2021.725041

  49. [49]

    Theory of Matrices

    Gantmacher, F.R. Theory of Matrices. Vol. 2. Chelsea publishing, New York, USA, 1959

  50. [50]

    Levenshtam, V. B. Justification of the averaging method for the convection problem with high- frequency vibrations, Siberian Math. J., 34:2 (1993), 280–296

  51. [51]

    Levenshtam, V.B. Justification of the averaging method for parabolic equations containing rapidly oscillating terms with large amplitudes, // Izvestiya: Mathematics, 2006, Volume 70, Issue 2, Pages 233–263 DOI: https://doi.org/10.1070/IM2006v070n02ABEH002311

  52. [52]

    Justification of the averaging method for a system of equations with the Navier–Stokes operator in the principal part

    Levenshtam, V.B. Justification of the averaging method for a system of equations with the Navier–Stokes operator in the principal part. St. Petersburg Mathematical Journal, 2015, Vol- ume 26, Issue 1, Pages 69–90 DOI: https://doi.org/10.1090/S1061-0022-2014-01331-3

  53. [53]

    Averaging a High-Frequency Hyperbolic System of Quasi- linear Equations with Large Summands

    Levenshtam, V.B. Averaging a High-Frequency Hyperbolic System of Quasi- linear Equations with Large Summands. Sib. Math. J. 65, 934–942 (2024). https://doi.org/10.1134/S0037446624040189

  54. [54]

    (1992) Homogenization and two-scale convergence

    Allaire, G. (1992) Homogenization and two-scale convergence. SIAM Journal on Math. Analysis, 23(6): 1482-1518

  55. [55]

    Homogenization of partial differential equations

    Marchenko V.A., Khruslov E.Y. Homogenization of partial differential equations. Springer Sci- ence & Business Media; 2008 Dec 22