Cannibalism is often an extreme interaction in the animal species to quell competition for limited resources. To model this critical factor, we improve the predator–prey model with nonlocal competition effect by incorporating the cannibalism term, and different kernels for competition are considered in this model numerically. We give the critical conditions leading to the double Hopf bifurcation, in which the gestation time delay and the diffusion coefficient were selected as the bifurcation parameters. The innovation of the work lies near the double Hopf bifurcation point, and the stable homogeneous and inhomogeneous periodic solutions can coexist. The theoretical results of the extended centre manifold reduction and normal form method are in good agreement with the numerical simulation.

Conventional preytaxis systems assume that prey-tactic velocity is proportional to the prey density gradient. However, many experiments exploring the predator–prey interactions show that it is the predator’s acceleration instead of velocity that is proportional to the prey density gradient in the prey-tactic movement, which we call preytaxis with prey-induced acceleration. Mathematical models of preytaxis with prey-induced acceleration were proposed by Arditi et al. ((2001) Theor. Popul. Biol. 59(3), 207–221) and Sapoukhina et al. ((2003) Am. Nat. 162(1), 61–76) to interpret the spatial heterogeneity of predators and prey observed in experiments. This paper is devoted to exploring the qualitative behaviour of such preytaxis systems with prey-induced acceleration and establishing the global existence of classical solutions with uniform-in-time bounds in all spatial dimensions. Moreover, we prove the global stability of spatially homogeneous prey-only and coexistence steady states with decay rates under certain conditions on system parameters. For the parameters outside the stability regime, we perform linear stability analysis to find the possible patterning regimes and use numerical simulations to demonstrate that spatially inhomogeneous time-periodic patterns will typically arise from the preytaxis system with prey-induced acceleration. Noticing that conventional preytaxis systems are unable to produce spatial patterns, our results imply that the preytaxis with prey-induced acceleration is indeed more appropriate than conventional preytaxis to interpret the spatial heterogeneity resulting from predator–prey interactions.

This is the second part of our study on the spatially heterogeneous predator–prey model where the interaction is governed by a Crowley–Martin type functional response. In part I, we have proved that when the predator competition is strong (i.e. k is large), the model has at most one positive steady-state solution for any $\mu \in \mathbb {R}$, moreover it is globally asymptotically stable for any $\mu >0$. This part is denoted to study the effect of saturation. Our result shows that the large saturation coefficient (i.e. large m) can not only lead to the uniqueness of positive solutions, but also lead to the multiplicity of positive solutions, moreover the stability of the corresponding positive solutions is also completely obtained. This work can be regarded as a supplement of Ref. [10].

Mealybug is an important pest of cassava plant in Thailand and tropical countries, leading to severe damage of crop yield. One of the most successful controls of mealybug spread is using its natural enemies such as green lacewings, where the development of mathematical models forecasting mealybug population dynamics improves implementation of biological control. In this work, the Sharpe–Lotka–McKendrick equation is extended and combined with an integro-differential equation to study population dynamics of mealybugs (prey) and released green lacewings (predator). Here, an age-dependent formula is employed for mealybug population. The solutions and the stability of the system are considered. The steady age distributions and their bifurcation diagrams are presented. Finally, the threshold of the rate of released green lacewings for mealybug extermination is investigated.

In this paper we consider a diffusive stochastic predator–prey model with a nonlinear functional response and the randomness is assumed to be of Gaussian nature. A large deviation principle is established for solution processes of the considered model by implementing the weak convergence technique.

We present a stability analysis of steady-state solutions of a continuous-time predator–prey population dynamics model subject to Allee effects on the prey population which occur at low population density. Numerical simulations show that the system subject to an Allee effect takes a much longer time to reach its stable steady-state solution. This result differs from that obtained for the discrete-time version of the same model.

We have developed a chemical kinetics simulation that can be used as both an educationaland research tool. The simulator is designed as an accessible, open-source project thatcan be run on a laptop with a student-friendly interface. The application can potentiallybe scaled to run in parallel for large simulations. The simulation has been successfullyused in a classroom setting for teaching basic electrochemical properties. We have shownthat this can be used for simulating fundamental molecular and chemical processes and evensimplified models of predator–prey interactions. By giving the simulated entities spatialextent in the lattice, the particles do not interpenetrate, and clusters of particles canspatially exclude one another. Our simulation demonstrates that spatial inhomogeneityleads to different results than those that are obtained by using standard ordinarydifferential equation models, as previously reported.