Turing-Hopf bifurcation of a delayed predator-prey system with prey-taxis incorporating a Beddington-DeAngelis functional response

When Alan Turing published his 1952 work on morphogenesis, he revealed something startling: identical cells, governed by identical rules, can spontaneously organize themselves into stripes, spots, and other patterns. All it takes is diffusion—the seemingly uniform, random mixing of molecules—acting on a system poised at the edge of instability. Seven decades later, that insight has migrated from chemistry into ecology, where biologists and mathematicians now ask whether animal populations, too, can self-organize into spatial structures without any external blueprint. A new theoretical study from the University of Science and Technology Beijing suggests they can, and that the recipe for such ecological pattern-making is far richer than anyone expected. The work, published in Frontiers in Physics by Yan Meng, Bingran Guo, and Jiaxin Xiao, dissects a predator-prey model in which two powerful forces—the deliberate movement of predators toward prey, and the inevitable delays in how ecosystems process energy—collide to generate an astonishing gallery of spatial and temporal patterns.

The mathematical skeleton of the study is a Leslie-Gower predator-prey system, a workhorse of population ecology that improves on the classic Lotka-Volterra equations by acknowledging that predators face environmental carrying capacity and logistic limits on their own growth. But the Beijing team layered on three crucial biological refinements. First, they replaced the naive linear feeding assumption with a Beddington-DeAngelis functional response, a predation law that captures how feeding rates saturate at high prey densities and, critically, how predators interfere with one another—competing for space, wasting time in confrontations, and thereby reducing each other’s hunting efficiency. Second, they incorporated prey-taxis: the directed, purposeful movement of predators up the gradient of prey density, a foraging strategy first formalized by Kareiva and Odell in 1987. Third, they added a time delay representing the lag between a predator consuming prey and that energy actually converting into new predator biomass—a delay that in real ecosystems reflects digestion, gestation, and resource regeneration.

Each of these ingredients, on its own, is known to destabilize uniform population distributions. Prey-taxis can either stabilize a homogeneous ecosystem—when predators’ attraction to prey smooths out spatial variations—or shatter it into heterogeneous fragments, depending on the direction and strength of the movement. Time delays are notorious for inducing Hopf bifurcations, the mathematical threshold at which a steady state loses stability and gives way to oscillations. What makes the new study distinctive is its systematic exploration of what happens when both mechanisms operate simultaneously. The researchers used the prey-taxis sensitivity coefficient and the time delay as bifurcation parameters, deriving precise conditions under which the system undergoes a Turing-Hopf bifurcation—a rare and mathematically intricate event where Turing instability (the pattern-forming mechanism Turing identified) and Hopf instability (the oscillation-triggering mechanism) occur at the same parameter point.

To understand the significance, it helps to recall what a Turing bifurcation does in an ecological context. In a reaction-diffusion system, Turing instability arises when diffusion, which normally erases differences, paradoxically amplifies them—creating stable spatial patterns from an initially uniform state. The Beijing team derived the conditions for Turing bifurcation in their prey-taxis model, identifying the critical wave number at which the first unstable mode emerges. They then went further, employing weakly nonlinear analysis to derive an amplitude equation that classifies the bifurcation as supercritical or subcritical. In a supercritical bifurcation, the emerging pattern is stable and grows continuously from the homogeneous state; in a subcritical bifurcation, the pattern is unstable and the system may jump abruptly to a different regime. This distinction matters enormously for predicting how real ecosystems will respond to changing conditions.

The theoretical predictions were validated by numerical simulations on a two-dimensional grid, using an explicit Euler scheme with finite-difference discretization of the Laplacian and careful treatment of the nonlinear taxis term. The results reveal a remarkable spectrum of spatiotemporal patterns. When the prey-taxis sensitivity is below the Turing threshold, the system settles into homogeneous steady patterns—uniform distributions with no spatial structure. Crossing that threshold triggers non-homogeneous steady patterns: stripes and spots that persist over time. As the taxis parameter increases further, regular stripes break down into spot patterns, which the authors interpret as temporally and spatially periodic oscillations. The simulations show that repulsive prey-taxis, where predators move away from high prey density, is particularly effective at inducing spatial heterogeneity, sometimes forming isolated refuge structures within specific parameter ranges.

But the most dramatic patterns emerge when time delay joins the picture. Below the critical delay, the system exhibits spatially heterogeneous but temporally stable patterns—stripes that persist indefinitely. As the delay crosses its Hopf threshold, those stable stripes destabilize and give way to non-homogeneous periodic patterns: structures that oscillate both in space and time. Pushing the delay further produces mixed wave patterns, a chaotic interplay of traveling waves and standing oscillations. The spatiotemporal series reveal population densities at fixed spatial points oscillating periodically or irregularly, while spatial snapshots show the ecosystem transitioning from uniform distributions through regular stripes to irregular spot patterns. This progression—from homogeneous stability through Turing-driven spatial patterns to Turing-Hopf-driven spatiotemporal chaos—offers a mathematical narrative of how ecological complexity can arise from simple local rules.

One finding stands out for its ecological implications: increased predator interference in the Beddington-DeAngelis response leads to a more even spatial distribution of populations. In other words, when predators spend more time interfering with one another—competing, confronting, wasting energy in mutual antagonism—the resulting patterns become smoother and more uniform. This suggests that predator social behavior, often overlooked in population models, plays a crucial role in shaping the spatial architecture of ecosystems. It also hints at a potential stabilizing mechanism: interference among predators may act as a buffer against the fragmentation of populations into isolated patches, with consequences for biodiversity, refuge formation, and the persistence of species in heterogeneous landscapes.

The study’s framework extends beyond predator-prey dynamics. Prey-taxis is mathematically analogous to chemotaxis, the directed movement of cells toward chemical signals that governs bacterial foraging, immune cell migration, and tumor growth. Time delays appear in neural systems, epidemiological models, and gene regulatory networks. The Turing-Hopf bifurcation, as a mathematical object, describes the collision of spatial and temporal instabilities in any reaction-diffusion system with delayed feedback. By deriving explicit conditions for this bifurcation in a biologically motivated model, the Beijing team has contributed a tool that pattern-formation researchers across disciplines can adapt. The saturation form of the prey-taxis coefficient—which limits predator movement speed at high densities—adds another layer of realism, capturing the finite perceptual and locomotor abilities of real organisms.

What makes this work compelling is its insistence that neither prey-taxis nor time delay alone can account for the full richness of ecological pattern formation. Comparing their results with models that consider only taxis-diffusion or rely solely on time delay, the authors demonstrate that the combined effect of these two mechanisms is necessary to generate the full spectrum of observed patterns—from homogeneous steady states through non-homogeneous periodic oscillations to mixed wave chaos. The interplay between directed movement and delayed energy conversion, they argue, profoundly enriches the dynamical complexity of predator-prey systems, offering a valuable framework for studying pattern formation in other ecological models that incorporate taxis mechanisms and delays.

As ecologists grapple with real-world questions—how climate change alters the spatial distribution of species, how habitat fragmentation affects predator-prey interactions, how to predict regime shifts in ecosystems—mathematical models like this one provide a crucial bridge between abstract theory and observable pattern. The Beijing team’s work reminds us that the patterns we see in nature, from the stripes of vegetation in arid landscapes to the patchy distributions of plankton in the ocean, may arise not from external forces but from the internal logic of movement, delay, and interaction. Turing’s legacy, it seems, is still unfolding—one bifurcation at a time.

Subject of Research: Science News Article

Article Title: Turing-Hopf bifurcation of a delayed predator-prey system with prey-taxis incorporating a Beddington-DeAngelis functional response

Article References: Meng, Y., Guo, B., & Xiao, J. (2026). Turing-Hopf bifurcation of a delayed predator-prey system with prey-taxis incorporating a Beddington-DeAngelis functional response. Frontiers in Physics, 14, Article 1919666. https://doi.org/10.3389/fphy.2026.1919666

Image Credits: AI Generated

DOI: 10.3389/fphy.2026.1919666

Keywords: predator-prey dynamics, Turing-Hopf bifurcation, prey-taxis, time delay, pattern formation, Beddington-DeAngelis functional response, reaction-diffusion, ecological modeling, Turing-Hopf, bifurcation, delayed, predator-prey