Evolutionary foraging in grids NeurIPS 2026

Evolutionary foraging in grids: Intermittent search dynamics emerge in finite, depletable landscapes

Shailendra Bhandari · Alex Szorkovszky · Anis Yazidi · Pedro G. Lind

Accepted at NeurIPS 2026 · arXiv:2609.39239 · DOI

Source code ↗ Paper (PDF) Coming soon Poster Coming soon arXiv ↗

Abstract

How search strategies evolve in finite, depletable landscapes remains a central question in foraging theory. We study this problem with an evolutionary simulation in which agents forage on a two-dimensional toroidal lattice containing non-renewable resources distributed either uniformly or as Lévy dust. Each agent carries a heritable genome encoding step lengths, velocities, and turning angles, and selection acts on a fitness function combining energetic gain, movement cost, and coverage efficiency. By allowing movement traits to evolve without imposing a prescribed power-law step-length distribution, we test whether the evolved trajectories are better described by intermittent-search or Lévy-walk dynamics. Strikingly, our results indicate that, in the finite depletion-driven landscapes considered here, evolved search is more consistent with intermittent dynamics than with strict scale-free Lévy motion. To characterize the effective dynamics of the evolved trajectories, we take the average adjusted coefficient of determination when fitting second- and fourth-order displacement moments to intermittent search and Lévy-walk models. While a Lévy-like random walk fits with the numerical results from the evolutionary search very well (mean adjusted R̄² > 0.9 in most tested conditions), the intermittent search achieves a closer fit, with mean adjusted coefficients of determination (R̄² > 0.99) for all tested resource distributions. This preference holds across the tested grid sizes and resource densities. Five independent evolutionary runs per environment on a 503 × 503 grid at nominal resource density ρ = 0.15 reproduce this preference for the uniform environment and the five Lévy-dust environments. Evolution rapidly reshapes the movement genome toward short displacements while retaining a sparse tail of longer relocations, consistent with local exploitation punctuated by occasional transfer. The framework provides a controlled setting for studying how search rules emerge under resource limitation and may inform resource-constrained exploration in autonomous systems. Project webpage: https://evo-foraging.github.io/

Resource landscapes and best evolved trajectories

Figure
Resource landscapes and best evolved trajectories
Figure: Resource landscapes and best evolved trajectories for the L = 503 torus. Top row: placement-step distributions used to generate the resource fields, shown for the uniform control and for Lévy-dust environments with µ ∈ {1.0, 1.5, 2.0, 2.5, 3.0}. Main panels: resource locations (gray points) overlaid with the trajectory of the best evolved agent in each environment (colored line). Changing µ changes the spatial correlation structure of the resource field by altering the frequency of long placement steps: smaller µ gives heavier-tailed placement steps, whereas larger µ suppresses long placement steps. The inset in the µ = 3.0 panel shows a magnified local patch.
Original full-resolution PNG ↗

Method

Here, we develop an evolutionary simulation framework that combines genome-encoded movement traits, finite lifetimes, explicit resource depletion, and controlled resource clustering. Agents forage on two-dimensional toroidal lattices containing non-renewable resources distributed either as Lévy dust or uniformly. Selection acts on a fitness function that rewards energetic efficiency and low-redundancy spatial coverage. We analyze how movement genomes evolve across resource environments, classify evolved trajectories using second- and fourth-moment fits to LW and IS models, and quantify waiting-time statistics under depletion. This framework allows us to ask whether adaptive search in finite, depletable landscapes converges toward strict scale-free LW dynamics or toward IS, combining local exploitation with occasional relocation.

Figures

Evolution of fitness

Figure
Evolution of fitness
Figure: mean fitness over 1500 generations for each resource environment, with shaded standard errors across agents within a single evolutionary run at each generation.
Original vector PDF ↗

Step-length genome distributions

Figure
Step-length genome distributions
Figure: distributions of step-length genome entries at generations 0, 200, and 1500. Selection rapidly concentrates the encoded step-length pool at short values while retaining a sparse tail of longer entries.
Original vector PDF ↗

Moment-based trajectory comparison

Figure
Moment-based trajectory comparison
Figure: Moment-based comparison of evolved trajectories on the 503 × 503 grid. Empirical second moments m2 and fourth moments m4 are compared with best-fit IS and LW models. Columns correspond to the six resource environments. Across environments, IS gives a closer fit to the empirical moment curves.
Original vector PDF ↗

Movement traits and return intervals

Figure
Movement traits and return intervals
Figure: Movement traits and return-interval statistics of the best evolved agents across grid sizes after 1500 generations. Columns correspond to resource environments. Top row: velocity-genome entries of the best evolved agents for L = 211, 503, and 1009. Middle row: turn-angle increments recorded along the best evolved trajectories. Bottom row: running maximum of return intervals Mjreturn versus return-event index j.
Original vector PDF ↗

Resource maps and trajectories at high density

Figure
Resource maps and trajectories at high density
Figure: best evolved trajectories overlaid on the corresponding resource maps for the highest tested density, ρ = 0.45.
Original full-resolution PNG ↗

Paper and code

Shailendra Bhandari, Alex Szorkovszky, Anis Yazidi, and Pedro G. Lind.

arXiv:2609.39239 · 10.48550/arXiv.2609.39239

Source code ↗ Paper (PDF) Coming soon Poster Coming soon arXiv ↗