- Research
- Open access
- Published:

# Modeling the movement of *Oecophylla smaragdina* on short-length scales in an unfamiliar environment

*Movement Ecology*
**volume 11**, Article number: 64 (2023)

## Abstract

The movement of individual weaver ants, of *Oecophylla smaragdina*, was previously tracked within an unfamiliar arena. We develop an empirical model, based on Brownian motion with a linear drag and constant driving force, to explain the observed distribution of ants over position and velocity. Parameters are fixed according to the isotropic, homogeneous distribution observed near the middle of the arena. Then, with no adjustable parameters, the model accounts for all features of the measured population distribution. The tendency of ants to remain near arena edges is largely explained as a statistical property of bounded stochastic motion though evidence for active wall-following behavior appears in individual ant trajectories. Members of this ant species are capable of impressive feats of collective action and long-range navigation. But we argue that they use a simplistic algorithm, captured semi-quantitatively by the model provided, to navigate within the confined region.

## Introduction

The social coherence of ant colonies has fascinated naturalists for centuries [1, 2]. Members of a colony work together to complete elaborate tasks. It requires the combined effort of many weaver ants, of *Oecophylla smaragdina* [3, 4] from tropical Asia and Australia, to bend and glue tree leaves while building their nest. This species, along with one other in the genus, is also notable for its colonies use of ‘living bridges’, made of hundreds of individuals, to negotiate gaps. Such collective action is only possible when individuals move in a precise manner and communicate, on both short and long length scales [5,6,7], to influence the motion of others. The weaver ants are known to employ a sophisticated system of communication, combining chemical signals with body motion to deliver distinct directives [3, 4, 8,9,10].

There is a large volume of literature on how the motion of one ant is influenced by messages received from another [2, 8, 11] and on the tools that an individual ant uses to navigate [12,13,14,15,16,17,18,19,20,21,22,23]. Desert ants [24, 25] in an unfamiliar environment use their view of the panorama [26] and path integration of their vector displacement [27] to navigate back to the nest. When foraging, they may carry out systematic searches of the local territory, moving in widening circles like other insects [28, 29], or merge navigational tactics with searching algorithms when they have incomplete knowledge of their whereabouts [30,31,32].

The movement of insects that do not have a certain destination is often modeled as a generalised random walk [33,34,35,36] similar to Brownian motion. Observed by Brown in 1827 and explained theoretically by Einstein and Langevin in the early 1900s [37,38,39,40,41], Brownian motion is exemplified by a small grain immersed in water that moves erratically because it receives impulses from collisions with neighboring molecules. Its velocity changes, like those of a crawling ant, are difficult to predict so they are treated as random quantities. While having this stochastic component, ant movement is also partly deterministic because the ant responds predictably to features in the local environment [42].

Ants are among many animals, including other insects [43], fish [44], rodents [45] and even humans during an evacuation [46] that exhibit wall-following. This is the tendency to move along a one-dimensional inhomogeneity within a two-dimensional space [47,48,49] (when the animal maintains bodily contact with the wall, it is termed thigmotaxis). The behavior has plausible adaptive value since a physical wall offers partial shelter and any 1D feature can be used as a directional guide [48,49,50]. Faced with a choice between two otherwise identical bridges, black garden ants select the one with a wall along its edge [51]. Rock ants tends to move parallel to distant walls [52] while desert ants follow the mid-lines between rows of shrubs [53]. Models of the Brownian motion of cockroaches [54] and harvester ants [55] in a confined space include position-dependent terms that represent an affinity for arena boundaries [56].

We previously carried out an experimental study [57] of the motion of a single weaver ant confined to a featureless, unfamiliar arena. An individual ant was snatched from a nest, placed on a square ceramic tile bounded by water, and its motion tracked. The only macroscopic inhomogeneity of the arena is its boundary, to which the ant tended to remain close. We chose the weaver ant because of its propensity for correlated action and its availability in our area.

In this paper we develop a model, based on Langevin theory, of the movement algorithm exhibited by the weaver ants and compare the results of numerical simulations and analytical calculations to experiment. Random velocity changes of model ants are governed by a fixed probability distribution augmented by deterministic drag and driving forces along the direction of motion. After it encounters the arena boundary and is forced to stop, the model ant immediately resumes its fixed movement algorithm. The model reproduces qualitative features of the observed distribution of ants over position and velocity without adjustable parameters.

We use this model to investigate whether an enhanced ant density at arena edges is evidence that ants follow walls actively, by modifying their movement algorithm when they are close to the edge and/or preferentially moving towards this edge, or whether it can be explained as a result of passive stochastic motion bounded by the arena. Our model includes no active wall-following behavior: model ants use the same movement algorithm near edges as in the open arena and are no more likely to turn towards a boundary as away from it. They nonetheless exhibit a high probability to be found within a short distance from the arena edge, in agreement with experiment. Going beyond this description of the population distribution, we consider below the residence times of individuals in the near-edge region, for both real and model ants, to isolate active wall-following behavior and evaluate its significance.

The paper is organized in the following way. In section “Experimental Results” we summarize the results of our previous experiment and discuss a new analysis of the ant distribution data. In section “Model of ant motion” we develop a model based on the distribution of ants in the arena interior, away from the boundaries. In section “Analytic approximation” we study the model analytically in order to better understand its key properties. In section “Results and discussion” we present the results of numerical simulations of the model and compare them to the experimental data in detail. We discuss the implications of the work in section “Results and discussion” before concluding in section “Conclusions”.

## Experimental results

Previously we described an experiment (see Ref. [57] for details) in which individual weaver ants, from *Oecophylla smaragdina*, were captured from one of many nests located on campus at Suranaree University of Technology, Thailand. All ants were the aggressive [58] female worker ants that comprise the majority of the colony. Ants were transported to the laboratory and one was released into an arena, a dry ceramic tile surrounded by water, where its motion was tracked via a stationary camera for five minutes. The time between an individual being captured and its observation completed was between ten and forty minutes, after which it was returned to its nest. The ant is unlikely to be motivated by hunger during the trial, but we do not speculate about its purpose for moving about the arena. Its motion was nearly continuous, without apparent change in its nature, throughout the trial.

A single ant was tracked as it moved on a square tile of length \(L=30\) cm centered on the origin and surrounded by a channel of water that it never attempted to cross. The tile and camera were located at a fixed position within a small room. The arena is defined by \(|x|<L/2\) and \(|y|<L/2\) where the orientation of the *x* and *y* axes was parallel to the edges of the tile. The ant position \({{\textbf {x}}}(t)=(x[t],y[t])\) was measured at time intervals \(\Delta t=1/15\) s over a duration \(T=4500\Delta t=300\)s, and the velocity \({{\textbf {v}}}(t)=({{\textbf {x}}}[t+\Delta t]-{{\textbf {x}}}[t])/\Delta t\) and change in velocity \(\Delta {{\textbf {v}}}(t)={{\textbf {v}}}(t+\Delta t)-{{\textbf {v}}}(t)\) were determined. The experiment was repeated on \(N=60\) individual ants, using many similar tiles, each washed and re-used multiple times.

Sample data, showing the trajectory of \({{\textbf {x}}}(t)\) for one ant over the full duration *T*, is shown in the top left panel of Fig. 1. The ant was released into the arena by inverting a plastic container near the middle of the arena, waiting for the ant to crawl out, then removing the container and starting the camera. For the trial shown, the first position measurement was \({{\textbf {x}}}(0)\approx (8,8)\) cm and the ant mainly stays close to the arena boundary, which was typical.

All data were combined into a single set, treated as a statistical distribution. By studying single ant motion using this distribution, we are averaging over the peculiarities of individuals. The *j*th ant has a position \({{\textbf {x}}}_j(t)\) and corresponding velocity \({{\textbf {v}}}_j(t)\), where \(t=0,\Delta t,2\Delta t,..T\). We have an experimental probability distribution defined by

The time dependence of \(\Pi ({{\textbf {x}}},{{\textbf {v}}},t)\) was weak for all \(t>>\Delta t\), so the time-averaged quantity \(\Pi ({{\textbf {x}}},{{\textbf {v}}})=T^{-1}\int dt \Pi ({{\textbf {x}}},{{\textbf {v}}},t)\) can be interpreted as an equilibrium value. Equation 1 is equivalent to a Boltzmann distribution function [59] in many-particle physics, but is here used to describe the motion of many individual ants moving in the arena at different times (i.e. different trials). The model below cannot be applied to a population of ants in the arena together because it includes no ant-ant interactions.

The density \(n({{\textbf {x}}})\) and velocity distribution \(P({{\textbf {v}}})\) are defined by

The limits of the velocity integrals are \(\pm \infty\), those of the position integrals are \(\pm L/2\), and \(\Pi ({{\textbf {x}}},{{\textbf {v}}})\), \(n({{\textbf {x}}})\), and \(P({{\textbf {v}}})\) are all normalized probability distributions.

A two-dimensional map of \(n({{\textbf {x}}})\) is shown at the top right of Fig. 1. The density decreases rapidly with distance from the boundary, approaching a value \(n({{\textbf {x}}})\approx n_{\textrm{int}}\) that is constant throughout the interior. Ants spent only 20 percent of their time at positions further than 3 cm from the boundary, with \(n_{\textrm{int}}\approx 4\cdot 10^{-4}\)cm\(^{-2}\). They spent most of their time near the boundaries: within 3 cm of \(x=\pm L/2\) and \(y=-L/2\) the average density was \(5.3 n_{\textrm{int}}\) while within 3 cm of \(y=L/2\) it was \(8.2 n_{\textrm{int}}\). The preferred \(y=L/2\) boundary was positioned closest to the window in the laboratory room.

The map of the velocity distribution \(P({{\textbf {v}}})\) is shown at the bottom of Fig. 1. The data has been broken up into two subsets: (i) ants in the interior and (ii) within 3 cm from a boundary. The former distribution is isotropic and depends only on speed \(v=|{{\textbf {v}}}|\). *P*(*v*) initially decreases to a local minimum then increases, exhibiting a circular shoulder, before falling off at high speed. For ants near the boundary, \(P({{\textbf {v}}})\) has four arms that correspond to members moving parallel to the square edges.

The plots in Fig. 1 provide an overview of ant behavior during the experiment. We develop a model of this behavior using the measured distribution of velocity changes within the arena interior. Data for the interior is the minority, but is simplest to analyze because motion is not affected by boundaries. So, we consider \(\Delta {{\textbf {v}}}\) measured when ants were more than 3 cm from the boundary. In coordinates relative to velocity \(\hat{{{\textbf {v}}}}\), it is given by

where \(\Delta v_T\) affects direction, with a positive value corresponding to a left turn, and \(\Delta v_L\) changes speed (the vector \(\hat{{\textbf {z}}}\) points up from the arena).

We use the same function *p*(*q*) to fit the distribution over \(\Delta v_T=q\) or \(\Delta v_L=q\), it is

where the mean \(\bar{q}\) and standard deviation \(\sigma\) are the fitting parameters. The distribution over \(\Delta {{\textbf {v}}}\) depends weakly on position within the arena interior. It does vary according to how fast the ant was moving when the velocity change occurred, so we sub-divided data for the interior according to speed *v*. The results are shown Fig. 2.

For velocity changes normal to motion \(\Delta v_T\), shown on the left side of Fig. 2, Eq. 4 appears to provide a good fit to the distribution. In the insets we plot the best-fit values for \(\sigma\) and \(\bar{q}\) versus speed. The standard deviation has an average value of \(0.996\pm 0.007\) cm/s, where the uncertainty is statistical. It can be approximated by a constant \(\sigma =\sigma _T\approx 1.00\) cm/s for, while it deviates from this value at low speeds, the dominant fraction of the data falls within the region where \(\sigma\) is speed-independent. The best-fit value for the mean \(\bar{q}\) is always zero within error.

We checked whether \(\Delta v_T(t)\) values for different times are independent by calculating the Pearson correlation coefficient [60] of \(\Delta v_T(t)\) with \(\Delta v_T(t+3\Delta t)\) and found a value of 0.01, indicating weak correlations. (We could not consider \(\Delta v(t)\) values more closely-spaced in time because, according to their definition, \(\Delta t\Delta v(t)\equiv x(t+2\Delta t)-2x(t+\Delta t)+x(t)\), they share *x*(*t*) data points.) So \(\Delta v_T(t)\) for each time *t* can be regarded as an independent variable with a distribution *p*(*q*), a mean \(\bar{q}=0\) and root-mean-square \(\sigma _T\).

The distribution of changes parallel to velocity \(\Delta v_L\) is shown on the right in Fig. 2. The fit with Eq. 4 appears reasonable, and the standard deviation \(\sigma =\sigma _L\approx 1.25\)cm/s can again be treated as a constant (its average value was \(1.251\pm 0.008\) cm/s). But the mean \(\bar{q}\ne 0\) has a characteristic *v* dependence, plotted in the upper-right inset. It increases with *v* to a small positive maximum before decreasing linearly into negative values. The majority of the data set falls within the range of speed where \(\bar{q}\) is linearly decreasing and the linear fit shown is relevant. Below, we will assume that \(\Delta v_L\) is composed of two terms, one random and the other deterministic. The latter results in the non-zero mean of \(\Delta v_L\) and some positive correlation between velocity changes at nearby times.

## Model of ant motion

The distribution of velocity changes in the arena interior suggests that a simple model, based on Brownian motion, may be applicable. We assume that the observed motion of the ant in the interior is characteristic of its movement algorithm. The arena boundaries will be treated as constraints that do not otherwise modify movement. A model ant has a position and velocity that can be updated in simulations according to

where the second equation corresponds to Newton’s law for a particle of unit mass. The term \({\textbf {F}}_j(t)={\textbf {F}}({\textbf {v}}_j[t])\) is a deterministic force while \({\textbf {q}}_j(t)\) is a random impulse occurring each time step. Equation 5 is the basic rate equation of the Langevin theory of Brownian motion.

The force is written as

where the first term is a linear drag force, with a time constant \(\tau _D\), and the second is a constant driving force in the forward direction, with a speed constant \(v_0\). The impulse \({\textbf {q}}_j(t)={\textbf {q}}=q_T\hat{{\textbf {z}}}\times \hat{{\textbf {v}}}+q_L\hat{{\textbf {v}}}\) where \(q_L\) is a random number governed by Eq. 4 with \(\sigma =\sigma _L\) and \(\bar{q}=0\) while \(q_T\) is a random number governed by Eq. 4 with \(\sigma =\sigma _T\) and \(\bar{q}=0\). Since the mean impulse \({\textbf {q}}\) is zero, an average over all model ants with a given velocity \({\textbf {v}}_j(t)={\textbf {v}}\) results in

where angular brackets denote this average. We can use the linear fit to the mean \(\Delta {\textbf {v}}\), shown in the insets of Fig. 2, to obtain \(\tau _D=0.55\) s and \(v_0=6.3\) cm/s. The effects of the deterministic force are small during a single time step, since \(|{\textbf {F}}|\Delta t<<\sigma _L,\sigma _T\), but are important on long time scales since they add constructively.

While we have based the movement algorithm on ant motion in the interior, we have to model the response of an ant to the boundary in order to compare simulations with data. We used the following protocol for the ant-boundary interaction. In simulations, if \({{\textbf { x}}}(t+\Delta t)\) is found outside the arena then we replace it with the nearest position on the boundary and set the velocity component normal to the boundary equal to zero. The position and velocity parallel to the boundary are not modified. With the subsequent time step, position and velocity are updated via Eq. 5 with the only difference being that the random velocity change cannot take the ant immediately back out of the arena: if the boundary is on the ant’s right side then \(q_T(t+\Delta t)\) must be positive. This latter constraint, a detail that does not change any qualitative results, was used to prevent model ants from briefly sliding along the boundary, which seemed artificial.

The response of the model ant to the boundary is a key feature of our study, so we briefly consider a few alternative approaches. We could have analysed in detail the observed motion of the ant near the boundary and attempted to model it. This would require numerous empirical parameters and reduce the entire exercise to a fit of measured motion within our particular experimental configuration. We sought instead a minimal model that might have some predictive value. That said, while we updated the model ant velocity at the \(x=L/2\) boundary according to \((v_x,v_y)\rightarrow (0,v_y)\), there are other simple procedures that could have been used. One would have model ants reflecting elastically off boundaries like billiard balls, \((v_x,v_y)\rightarrow (-v_x,v_y)\). This was dismissed because it cannot give rise to the measured distribution and seems to take the analogy with Newtonian particles to an absurd extreme. Another, to have ant speed vanish at the boundary \((v_x,v_y)\rightarrow (0,0)\), was rejected because it does not allow rapid motion along the boundary, which is seen at a first glance of the experiment. Any more sophisticated response of the ant to a boundary encounter, such as pausing for some time or changing direction in a particular way, would employ more fitting parameters.

Finally, one could introduce wall-following to the model by hand using a position-dependent deterministic force \({\textbf {F}}({\textbf {x}},{\textbf {v}})\) that attracts ants to the boundary at short distances. We avoided this for two reasons. First, such a force would introduce new fitting parameters that characterize its strength and range. Second, we wanted to investigate the possibility that the observed distribution of ants could be realized by a model that did not include any active wall-following behavior. Here we define *active* wall-following as a position-dependent algorithm that favors motion towards or along the boundary. In contrast, a model like ours that predicts an increased ant density at boundaries without employing a position-dependent movement algorithm could be said to exhibit *passive* wall-following behavior. Overall, we view the protocol we adopted as the simplest one that had any chance of reproducing the observed distributions. The model ant is forced to stop at the boundary but then resumes its motion, using the same algorithm it employs in the open arena, as if nothing happened.

The simulation was carried out many times and the distribution extracted from the numerical data and time-averaged. We usually started each simulated trial with the initial conditions \({\textbf {x}}_j(0)={\textbf {v}}_j(0)=0\) and the ant’s body orientation equally likely to be in any direction, but the equilibrium distribution was independent of initial conditions. The values of the model parameters: \(\sigma _T\), \(\sigma _L\), \(\tau _D\), \(v_0\) were all obtained experimentally from the distribution of \(\Delta {\textbf {v}}\) in the arena interior. The model has no free parameters and the fixed parameters can only be properties of the movement algorithm for an ant in an unbounded arena. A comparison between the simulated and measured distributions in the bounded arena is thus meaningful, i.e. the model is falsifiable.

## Analytic approximation

### The rate equation

An approximate rate equation for the model detailed above is presented in this section. The goal is to obtain some analytic understanding of its properties. Particularly, we explain the length and speed scales of the model, which can be associated with those observed in the data. Using Eqs. 1, 5 and 6, we write the distribution at time \(t+\Delta t\) as

The probability that a given member receives an impulse \({\textbf {q}}\) within a time step is \(d{\textbf {q}}p({\textbf {q}})\). We have

where \(q_L\) and \(q_T\) are the components of \({\textbf {q}}\) that are parallel and perpendicular to \({\textbf {v}}\), respectively. The limits of the \(q_L\) and \(q_T\) integrals in Eq. 9 are \(\pm \infty\).

If we drop terms of order \(\Delta t^2\) and \(\Delta t|{\textbf {q}}|\) then the first Dirac delta function in Eq. 8 can be replaced by \(\delta ({\textbf {x}}_j+{\textbf {v}}\Delta t-{\textbf {x}})\). The second can be rewritten as

where we used the property \(\delta (f[x])=\delta (x-x_0)/|df/dx|_{x=x_0}\) with \(f(x_0)=0\). Substituting these expressions and rearranging, we have

Now carrying out an expansion to first order in \(\Delta t\) and second order in \({\textbf {q}}\) we get

The last term, the average over squared random impulses, is \(\Delta \Pi _{\textrm{R}}/\Delta t= \Lambda _{xx}+\Lambda _{yy}+\Lambda _{xy}\) with

where \({\textbf {v}}=(v_x,v_y)=v(\cos \theta ,\sin \theta )\).

It is convenient to express this equation in dimensionless coordinates. We have parameters with the units of speed \(v_\infty\) and distance \(\ell\), and a dimensionless measure of anisotropy \(\alpha\) that are

Dimensionless coordinates for velocity \({\textbf {u}}={\textbf {v}}/v_\infty\), position \({\textbf {r}}={\textbf {x}}/\ell\), time \(\tau =t/\tau _D\), and force \({\textbf {f}}={\textbf {F}}\tau _D/v_\infty =-{\textbf {u}}+u_0{\hat{{\textbf {u}}}}\), with \(u_0=v_0/v_\infty\), are introduced. The size of the arena is \(R=L/\ell\). The dimensions of \(\Pi ({\textbf {x}},{\textbf {v}},t)\) can be removed by multiplying it by \(v_\infty ^2\ell ^2\). The dimensionless rate equation for the distribution \(\Pi ({\textbf {r}},{\textbf {u}},\tau )\) is given by

The equilibrium distribution \(\Pi ({\textbf {r}},{\textbf {u}})\) solves Eq. 17 with the left side equal to zero.

The model equilibrium distribution has two speed scales, \(v_0\) and \(v_\infty \approx 2.5\)cm/s and one length scale \(\ell \approx 1.4\)cm. The value \(v_0\), at which the driving and drag forces balance, is a mean speed for members moving in a given direction. The effect of random impulses are contained within \(v_\infty\), the mean-square deviation of velocity from this mean. An ant remains approximately localized within a region of width equal \(\ell\) because it is unlikely to travel this far without turning around. The anisotropy factor is \(\alpha \approx 0.36\).

### Equilibrium distribution for a homogeneous, unbounded arena

We can study the model in the simple case of an unbounded arena with area \(R^2\). (The actual, more difficult, case of a bounded arena is treated approximately in Appendix B). This is done by imposing periodic boundary conditions and can be imagined as the 2D arena stretched over a closed 3D shape. The equilibrium distribution \(\Pi ({\textbf {r}},{\textbf {u}})\) is proportional to that for an arena of infinite size: every position and direction is symmetry-equivalent so the distribution, independent of \({\textbf {r}}\) and \(\theta\), depends only on speed *u*.

For isotropic impulses, \(\alpha =0\), the distribution satisfies

which is solved to give \(\Pi ({\textbf {r}},{\textbf {u}})=\Pi (u)\propto \textrm{e}^{-u^2/2+u_0u}\). (The second order differential equation has two independent solutions but we must choose the one that vanishes at large speed.) As a function of speed, in variables with units, this distribution has shoulder peaks at \(v=\pm v_0\), each with a half-width \(v_\infty\).

For anisotropic impulses, the equilibrium distribution is still independent of \(\theta\). This is ensured by symmetry if the initial state of the ensemble is uniformly distributed over all ant-body orientations. A normalized solution to the rate equation for \(\alpha \ne 0\) is

with the normalization constant

The anisotropy \(\alpha\) of the random impulses, while not having a huge effect on the distribution, shifts the position of the shoulders to slightly higher speed. They are centered on \(u=\pm v_s/v_\infty\) where

So the shoulder peak position \(v_s\approx 6.7\)cm/s

## Results and discussion

We compare the equilibrium distributions obtained from model simulations and experiment. First, consider the one-dimensional density and velocity distributions:

and their counterparts *n*(*y*) and \(P(v_y)\).

In the upper left main panel of Fig. 3 we show the model densities *n*(*x*) versus *x* and *n*(*y*) versus *y*. The inset shows the experimental results for the same quantities. Note the enhanced value for the measured density *n*(*y*) near \(y=L/2\), the deviation from square symmetry that we ascribed to the presence of the laboratory window. The model density respects square symmetry exactly.

The model density has the same qualitative behavior as the observed value: it is large at the boundary, then decreases over a length scale \(\ell\) to a value that is constant throughout the arena interior. The analytic calculation (see Appendix B) of the 1D version of the model suggests that *n*(*x*) has an integrable divergence at the boundary and decreases exponentially. That is, at a small distance \(d=L/2-x\) from the right boundary(or \(d=-L/2+x\) from the left) we have \(n(x)\sim d^{-2/3}\) for \(d<<\ell\) while \(n(x)-n(0)\approx \textrm{e}^{-d/\ell }\) for \(d>>\ell\).

This behavior of *n*(*x*) is simply understood. Any ant that arrives at a boundary tends to remain nearby because, after it stops, it is improbable that a sequence of impulses in the same direction (away from the boundary) will carry it deep into the interior. More probably, it wanders a short distance from the boundary before returning and stopping again. This gives an enhancement of *n*(*x*) within distance \(\ell\) of boundary. Members that are much further than \(\ell\) have little chance of wandering far enough in one direction to encounter the boundary, so they behave as if it did not exist.

The velocity distributions \(P(v_x)\) and \(P(v_y)\) are plotted in the upper right of Fig. 3, with the model result in the main panel and experimental values in the inset. The overall distribution \(P({\textbf {v}})\) exhibits square symmetry, so \(P(v_x)=P(v_y)=P(-v_x)\). The function \(P(v_x)\) is sharply peaked at \(v_x=0\) and initially decreases with \(|v_x|\). This peak is dominated by ants moving along \(y=\pm L/2\) boundaries that frequently have the normal component of their velocity reset to zero. After reaching a local minimum, \(P(v_x)\) recovers to exhibit shoulder peaks that are characteristic of the interior population. They are centered near \(|v_x|\approx v_s\) with a half-width \(v_\infty\). The shape of the simulated distribution resembles that of the experiment, at least all features are represented. From Appendix A, in which we allow model parameters to vary, we get an indication that the observed \(P(v_x)\) function might be fully reproduced within this model framework.

A 2D representation of \(P({\textbf {v}})\) is illustrated in the lower panels of Fig. 3 using contour plots for the model (main panels) and experiment (insets). On the left we show \(P({\textbf {v}})\) for ants in the arena interior (more than 3 cm from the boundary). The distribution is isotropic, with a valley and ring-shaped plateau seen at finite speed (this is the 2D realization of the shoulders). The experimental result, while more ragged, has the same qualitative character. On the right we plot \(P({\textbf {v}})\) for ants in the 3 cm wide boundary strip. Since they remain close to the boundary, they move rapidly along it, giving the characteristic four-lobed shape.

It should be emphasized that model parameters were obtained from the homogeneous, isotropic distribution of \(\Delta {\textbf {v}}\) within the arena interior. The variation of \(n({\textbf {x}})\) near the boundary, and features of \(P({\textbf {v}})\), including any direction-dependence, were not assumed within the model but emerged as calculated properties of it. The basis for our claim that the model is in semi-quantitative agreement with the data is that the model length scale \(\ell\) and speed scales \(v_\infty\), \(v_0\) are consistent with observed values.

The density \(n({\textbf {x}})\), being strongly peaked at the arena boundaries, exhibits the signature of wall-following, a well-known property of ant motion [51]. An animal that exhibits active wall-following behavior should show a marked preference to remain near the wall and move differently when close to a wall than it does in open space. In Ref. [54] the motion of a cockroach was parameterized separately for phases when the animal was in the arena interior and when it was within antenna-length of the boundary. The harvester ants studied in Ref. [55] were observed to spend most of their time in the boundary region when confined to an arena similar to ours. Their corresponding model included parameters characterising diffusive motion in the arena interior and independent parameters associated with the probability of an ant leaving the boundary region. Note that the open-arena motion of Ref. [55] differs somewhat from our picture since it took velocity direction to change randomly while speed remained constant whereas, based on the measured distribution \(p(\Delta {\textbf {v}})\) for weaver ants, we varied both direction and speed.

Our model ants cannot walk beyond an arena boundary so they stop when they reach the edge. But when they resume moving they are just as likely to turn away from the boundary as towards it. In fact, we slightly biased the model by forcing the ant at the boundary to take its first step into the arena interior. The model ants have no preference to remain close to the boundary, it is simply that the diffusive motion they undergo in the open arena does not effectively lead them away from its edge. As such, the fact that the simulated density \(n({\textbf {r}})\) is peaked at the boundary results from passive wall-following behavior exhibited by model ants. But we can look further for evidence of active wall-following behavior in real weaver ants, i.e for a position-sensitive movement algorithm, by examining the distribution of resident times for individuals in the boundary strip.

We extracted the distribution \(P(t_b)\) for the time \(t_b\) an ant remains continuously in the 3 cm boundary strip after visiting the arena edge. Specifically, we logged all times \(t=t_1\) when the ant displacement achieved a local maximum |*x*(*t*)| within 0.5 cm of the \(|x|=L/2\) boundary and the time \(t=t_2=t_1+t_b\) when it subsequently moved beyond the 3 cm boundary strip. Accordingly, \(t=t_2\) is the earliest time \(t>t_1\) for which \(|x(t)<\)12 cm. Note that an ant whose center of mass comes within 0.5 cm of the arena edge, half the length of its body, can reasonably be said to have contacted the boundary [61]. We consider only the *x* coordinate for this distribution in order to avoid having to make special provisions for arena-corner effects.

The normalized distribution \(P(t_b)\) is shown in the left panel of Fig. 4 along with corresponding model prediction. The latter was obtained from simulation data using the same procedure followed for the experimental data. The complementary distributions \(P(v_x,t_1)\) of velocity perpendicular to the boundary at \(t=t_1\) are included in the right panel. Note that, since \(x(t_1)\) is a maximum, all ants have a velocity \(v_x(t_1)\) that is negative, i.e. towards the interior, or zero. The model distribution \(P(t_b)\) decays with time \(t_b\) on the scale of approximately \(\tau _D\). The experimental distribution also drops significantly over this time scale, confirming that the high density \(n({\textbf {r}})\) at the boundary can largely be explained by the passive wall-following behavior present in the model. But there is a long tail in the measured distribution, giving a significantly higher probability \(P(t_b)\) for large \(t_b\) values in the experiment than in the model simulation.

The tail of the \(P(t_b)\) distribution can be attributed to an additional tendency for active wall-following present among real weaver ants. As noted above, such an effect could be incorporated into the model using a position-dependent force term in Eq. 6. While our current interest is the population distribution shown in Fig. 3, we note that in performing fits of a parameterized force \({\textbf {F}}({\textbf {r}},{\textbf {v}})\) to data like that of Fig. 4, our model ants provide a control group that establishes the baseline distribution \(P(t_b)\) in the absence of active wall-following. The force parameters obtained for different experimental groups, perhaps those with ants following a tactile wall instead of the arena boundary used here, could thus be more meaningfully compared.

The movement algorithm that the ants employ during the conditions of the experiment appears to be simplistic. The ants can be, at most, weakly influenced by long-range sensory input, such as their view of the panorama, measurement of light polarization or other sources. We saw that \(P({\textbf {v}})\) and, to a lesser extent, \(n({\textbf {x}})\) exhibited square symmetry. (The enhanced \(n({\textbf {x}})\) for the arena edge nearest the window was a quantitative, not qualitative, symmetry violation.) This symmetry is an indication that ants were using local information to navigate, rather than looking beyond the arena for guidance. For, while the arena itself (an average over square tiles) was square-symmetric, the panorama of the laboratory was not.

There is also no indication that our ants were searching the arena in a systematic way, such as turning in widening circles. First, they do not favor one turning direction over another, as you would expect for those turning circles with a certain sense of rotation. More importantly, the correlation between consecutive impulses normal to velocity is negligible. If an ant was following a smooth trajectory then its velocity change during a given time step would be strongly correlated with that of the previous step. Instead, it appears that the ant is changing its velocity haphazardly.

The experiment subjected the weaver ant to unnatural conditions, and it is possible that the ant was unable to utilize its usual navigational tools within the confined arena. This might explain its apparently confused motion. The results above would then have little to say about the dispersion of weaver ants in nature. We are currently [62] using the same experimental configuration with a pair of ants, and studying whether the frequent interactions between the partners influence their individual motion. So, in this worst case, the model could still prove useful for studying communication among weaver ants.

On the other hand, if the observed motion is representative of weaver ant movement in natural situations then it raises questions about how motion described by the model above could be of any tactical advantage. In a homogeneous environment where it is difficult to keep track of one’s position, a systematic search pattern might be difficult to follow. An individual searching for a faint chemical trail on a barren landscape might behave like a Brownian particle, which tends to remain localized near a starting point but does not repeatedly follow the same path [63,64,65]. Such randomized motion would, eventually, leave no stone unturned.

## Conclusions

We considered the motion of an individual weaver ant moving in a square arena. Experimental data on the ant distribution over position and velocity were presented. We developed a model to explain the observed results. The distributions were non-trivial: the ant showed a strong tendency to remain near the arena boundaries and its distribution over speed was non-monotonic with several characteristic features.

The model assumed that the ant undergoes random velocity change according to a fixed probability distribution and that it has no preference whatsoever for position within the arena or direction of motion. With these minimal assumptions and without adjustable parameters, the model captures the data semi-quantitatively. The work provides insight into the movement algorithm used by ants when confined in an unfamiliar region.

## Availability of data and materials

Experimental data and numerical codes available upon request from authors.

## Change history

### 28 February 2024

A Correction to this paper has been published: https://doi.org/10.1186/s40462-024-00456-y

## References

Holldobler B, Wilson EO. The Ants. Cambridge: Harvard University Press; 1994.

Holldobler B, Wilson EO. Journey to the ants. Cambridge: Bellknap Press of Harvard University Press; 1994.

Cole AC, Jones JW. A study of the weaver ant,

*Oecophylla smaragdina*. Am Midl Nat. 1948;39:641.Holldobler B. Territorial behavior in the green tree ant. Biotropica. 1983;15:241.

McCreery H, Breed M. Cooperative transport in ants: a review of proximate mechanisms. Insectes sociaux. 2014;61:99.

Berman S. Sensors for gesture recognition systems. Proc IEEE. 2011;99:1470.

Feinerman O, et al. The physics of cooperative transport in groups of ants. Nat Phys. 2018;14(7):683.

Vander Meer RK, et al. Pheromone communication in social insects: ants, wasps, bees, and termites. Boca Raton: CRC Press; 2019.

Stepejns DW, Krebs JR. Foraging theory. Press: Princeton University; 1986.

Kempraj V, Park SJ, De Faveri S, Taylor PW. Overlooked scents: chemical profile of soma, volatile emissions and trails of the green tree ant,

*Oecophylla smaragdina*. Molecules. 2020;25:2112.Traniello J. Foraging strategies in ants. Ann Rev Entomol. 1989;34(1):191.

Collett M, Collett T, Bisch S, et al. Local and global vectors in desert ant navigation. Nature. 1998;394:269.

Wehner R. Desert ant navigation: how miniature brains solve complex tasks. J Comp Physiol. 2003;189:579.

Wehner R, Boyer M, Loertscher F, Sommer S, Menzi U. Ant navigation: one-way routes rather than maps. Curr Biol. 2006;16:75.

Knaden M, Wehner R. Uncertainty about nest position influences systematic search strategies in desert ants. J Exp Biol. 2006;209:26.

Muller M, Wehner R. Wind and sky compass cues in desert ant navigation. Naturwissenschaften. 2007;94:589.

Baddeley B, Graham P, Husbands P, Philippides A. A model of ant route navigation driven by scene familiarity. PLoS Comput Biol. 2012;8: e1002336.

Moller R. A model of ant navigation based on visual prediction. J Theor Biol. 2012;305:118.

Knaden M, Graham P. The sensory ecology of ant navigation: from natural environments to neural mechanisms. Ann Rev Entomol. 2015;61:63.

Fleischmann PN, et al. The geomagnetic field is a compass cue in

*Cataglyphis*ant navigation. Curr Biol. 2018;28:1440.Jander R, Jander U. The light and magnetic compass of the weaver ant,

*Oecophylla smaragdina*(Hymenoptera: formicidae). Ethology. 1998;104:743.Ogawa Y, et al. Physiological properties of the visual system in the green weaver Ant,

*Oecophylla smaragdina*. 2022. https://www.researchsquare.com/article/rs-1958013/v1.Amorim P, Goudon T, Peruani F. An ant navigation model based on Weber’s law. J Math Biol. 2019;78:943.

Schwartz S, Wystrach A, Cheng K. Ants’ navigation in an unfamiliar environment is influenced by their experience of a familiar route. Sci Rep. 2017;7:14161.

Wehner R, Hoinville T, Cruse H, Cheng K. Steering intermediate courses: desert ants combine information from various navigational routines. J Comp Phys A. 2016;202:459.

Wystrach A, Beugnon G, Cheng K. Landmarks or panoramas: What do navigating ants attend to for guidance? Front Zool. 2011;8(1):1.

Müller M, Wehner R. Path integration in desert ants,

*Cataglyphis fortis*. R Proc Nat Acad Sci. 1988;85(14):5287.Schultheiss P, Cheng K, Reynolds AM. Searching behavior in social hymenoptera. Learn Motiv. 2015;50:59.

Cheng K, Schultheiss P, Schwarz S, Wystrach A, Wehner R. Beginnings of a synthetic approach to desert ant navigation. Behav Process. 2014;102:51.

Collett M. Spatial Memory in insect navigation Current. Biology. 2012;22:927.

Legge ELG, Wystrach A, Spetch ML, Cheng KJ. Combining sky and earth: desert ants (

*Melophorus bagoti*) show weighted integration of celestial and terrestrial cues. Exp Biol. 2014;217:4159.Wystrach A, Mangan M, Webb B. Optimal cue integration in ants. Proc R Soc B-Bio Sci. 2015;282:20151484.

Bell WJ. Searching behaviour: The behavioural ecology of finding resources. London: Champman and Hall; 1990.

Benichou O, Loverdo C, Moreau M, Voituriez R. Intermittent search strategies. Rev Mod Phys. 2011;83:81.

Li L, Peng H, Kurths J, Yang Y, Schellnhuber HJ. Chaos-order transition in foraging behavior of ants. Proc Natl Acad Sci. 2014;111:8392.

Paiva L, et al. Scale-free movement patterns in termites emerge from social interactions and preferential attachments. Proc Nat Acad Sci. 2021;118(20): e2004369118.

Einstein A. On the motion-required by the molecular kinetic theory of heat-of small particles suspended in a stationary liquid. Ann d Physik 1V. 1905;17:549–60.

Langevin P. On the theory of Brownian motion. C R Acad Sci Paris. 1908;146:530.

Einstein A. Investigations on the theory of Brownian motion. London: Methue & Co. Ltd.; 1966.

Van Kampen NG. Stochastic processes in physics and chemistry. Amsterdam: Elsevier; 2007.

Wagner RJ, Such K, Hobbs E, Vernerey FJ. Treadmilling and dynamic protrusions in fire ant rafts. Soc Interface. 2021;18:20210213.

Turchin P. Quantitative analysis of movement. Sunderland: Sinauer Assoc; 1998.

Salazar MO, Planas-Sitja I, Sempo G, Deneubourg JL. Individual thigmotactic preference affects the fleeing behavior of the American cockroach (Blattodea: Blattidae). J Insect Sci. 2018;18:9.

Patton P, Windsor S, Coombs J. Active wall following by Mexican blind cavefish (

*Astyanax mexicanus*). Comp Physiol A. 2010;196:853.Lamprea MR, Cardenas FP, Setem J, Morato S. Thigmotactic responses in an open-field. Braz J Med Biol Res. 2008;41:135.

Xue S, Jiang R, Wong SC, Feliciani C, Shi X, Jia B. Wall-following behaviour during evacuation under limited visibility: experiment and modelling. Transp A Transp Sci. 2020;16(3):626.

Fraenkel G, Gunn DL. The orientation of animals. New York: Dover Publications; 1961.

Holldobler B, Taylor RW. A behavioral study of the primitive ant Nothomyrmecia macrops Clark. Insectes Soc. 1983;30:384–401.

Collett TS, Collett M, Wehner R. The guidance of desert ants by extended landmarks. J Exp Biol. 2001;204:1635.

Klotz JH, Reid BL. Nocturnal orientation in the black carpenter ant

*Camponotus pennsylvanicus*(DeGeer) (Hymenoptera: Formicidae). Insectes Soc. 1993;40:95.Dussutour A, Deneuborg JL, Fourcassie V. Amplification of individual preferences in a social context: the case of wall-following in ants. Proc Biol Sci. 2005;272(1564):705.

Pratt SC, Brooks SE, Franks NR. The use of edges in visual navigation by the ant

*Leptothorax albipennis*. Ethology. 2001;107:1125.Heusser D, Wehner R. The visual centring response in desert ants,

*Cataglyphis fortis*. J Exp Biol. 2002;205:585.Jeanson R, Blanco S, Fournier R, Deneubourg JL, Fourcassié V, Theraulaz G. A model of animal movements in a bounded space. J Theor Biol. 2003;225(4):443.

Casellas E, Gautrais J, Fournier R, Blanco S, Combe M, Fourcassié V, Theraulaz G, Jost CJ. From individual to collective displacements in heterogeneous environments. Theor Biol. 2008;250:424.

Weitz S, Blanco S, Fournier R, Gautrais J, Jost J, Theralaux G. Residence times and boundary-following behavior in animals. Phys Rev E. 2014;89: 052715.

Thiwatwaranikul T, et al. Modeling navigation by weaver ants in an unfamiliar, featureless environment. Phys Rev E. 2020;101: 052404.

Van Mele P, Cuc NT. Weaver ants,

*Oecophylla smaragdina*(Hymenoptera: Formicidae), as biocontrol agents on African mahogany trees. J Pest Manag. 2010;4:295.Abrikosov AA. Introduction to the theory of normal metals. New York: Academic Press; 1972.

Denoting by angular brackets a time average: \(<f(t)>\equiv T^{-1}\int _0^T dt f(t)\), the unitless Pearson correlation coefficient for the correlation between \(a\equiv \Delta v_T(t)\) and \(b\equiv \Delta v_T(t+3\Delta t)\) is \(r=\langle AB\rangle /\sqrt{\langle A^2\rangle \langle B^2\rangle }\) where \(A=a-\langle a\rangle\), \(B=b-\langle a\rangle\). Note that \(\langle a \rangle =\langle b \rangle =0\) and \(\langle A^2\rangle =\langle B^2\rangle =\sigma _T^2\). For the correlation between \(\Delta v_T(t)\) and itself we get \(r=1\). A value \(|r|<<1\) suggests the variables are independent.

The typical weaver ant in the experiment has a body length slightly less than 1 cm. Here we define an ant to have reached the arena boundary if its position comes to within 0.5 cm from \(|x|=15 cm\). It is difficut to be more precise: the tile edge is not perfect and the ant position, determined from its optical center of mass, often went slightly beyond the boundary as the ant swung its body around.

Paisanpan P, et al. Unpublished.

Reynolds AM. Scale-free animal movement patterns: Lévy walks outperform fractional Brownian motions and fractional Lévy motions in random search scenarios. J Phys A Math Theor. 2009;42: 434006.

Palyulin V, Chechkin A, Klages R, Metzler R. Search reliability and search efficiency of combined Lévy-Brownian motion: long relocations mingled with thorough local exploration. J Phys A Math Theor. 2016;49: 394002.

Khadem S, Klapp S, Klages R. Search efficiency of discrete fractional Brownian motion in a random distribution of targets. Phys Rev Res. 2021;3: 023169.

## Acknowledgements

This work was supported by Suranaree University of Technology (SUT).

## Funding

This work was supported by the Research Center for Theoretical Physics of Suranaree University of Technology; PP and LC received additional funding from DPST Scholarships granted by the Institute for the Promotion of Teaching Science and Technology, Thailand.

## Author information

### Authors and Affiliations

### Contributions

TT, PP, and SS did the original experiment using a tracking code written by SS; they shared their data and analysis throughout the project and helped develop the theoretical model. LC helped develop the theoretical model, wrote the simulation codes and did all numerical calculations and plots. MS helped develop the theoretical model, did the analytical calculations and wrote the paper.

### Corresponding author

## Ethics declarations

### Ethical approval and consent to participate

Not applicable.

### Competing interests

The authors declare that they have no Competing interests.

## Additional information

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

The original version of this article was revised: a spelling mistake in the title and in the Abstract, Introduction and Experimental results sections was corrected.

## Appendices

### Appendix A: variation of distribution with model parameters

To better understand the range of behavior possible for the model, we relax the requirement that parameters must be fixed to measured values. As described above, the density in the interior is constant by symmetry. It formally diverges at the boundary since the current of ants incident on the \(x=L/2\) boundary during a given time step will all be placed at \(x=L/2\) for the next time step: a finite number of members occupying a single point. The value of the basic distance scale \(\ell\), over which *n*(*x*) varies, changes as we modify parameters: larger mean-square impulses \(\sigma _L^2\) and a smaller drag force coefficient \(1/\tau _D\) increase the region that a given ant explores, extending the range over which boundaries influence the distribution. But the qualitative behavior of the *n*(*x*) does not change dramatically.

The distinctive shape of \(P(v_x)\) with round shoulder peaks coming from the interior and a sharp central peak from the boundaries, is sensitive to parameter values. In Fig. 5 we show the distribution for various parameter choices. In the upper panels we vary \(\tau _D\) and \(v_0\), which control the deterministic force, and in the lower panels vary \(\sigma _L\) and \(\sigma _T\), the size of the random impulses.

The shoulder peaks are centered on \(\pm v_0\) and have a width \(v_\infty\). Increasing the driving force \(v_0\) pushes them outward without changing their height. Upon decreasing \(v_0\), the shoulders can merge with the central peak, so the width of the latter increases significantly. If the drag coefficient \(1/\tau _D\) increases then accumulated random impulses cannot achieve large speeds so \(v_\infty\) is decreased. The result is the shoulder peaks become narrow, tall and more pronounced. On the other hand, members become more narrowly localized in space so \(\ell\) decreases and the boundaries affect a smaller region. The height of the central peak slightly decreases as a result. If \(1/\tau _D\) decreases and \(v_\infty\) gets large then the shoulders become wider and shorter, merging into a smooth distribution.

With the deterministic force fixed we can vary the size of the root-mean-square random impulses \(\sigma _L\) and \(\sigma _T\). Variation of \(\sigma _L\) has a similar effect to changing \(\tau _D\), as is clear from the two plots on the left of Fig. 5. This is because the distribution results from a competition between the deterministic force, which scales as \(1/\tau _D\), and the random impulses that scale as \(\sigma _L^2\). So, increasing one of these is similar to decreasing the other. Recall, in Eq. 19, the anisotropy of the impulses produced only logarithmic corrections to the distribution in the interior. As a result, when we change \(\sigma _T\) the shoulder features look the same, as seen in the lower right of Fig. 5. But the central peak is more sensitive to \(\sigma _T\). After a model ant hits the boundary its body is aligned parallel with it (because the normal component of its velocity is reset to zero). If \(\sigma _T\) is small, it cannot change direction effectively so it moves along the edge and remains in contact with the boundary. The height of the central peak, which indicates the weight of the boundary population, increases as a result.

### Appendix B: Equilibrium distribution for a 1D arena

Since the rate equation for the bounded 2D model is more difficult to solve, we consider a simpler 1D model in this section. The solution for the 1D equilibrium distribution reproduces the key properties of the density *n*(*x*) and velocity distribution \(P(v_x)\) found above in the 2D simulations and experiment. The 1D model can be viewed as the limiting case of extremely anisotropic impulses for the 2D model. We set \(\sigma _T=0\) so an ant initially oriented along the *x* axis remains so.

The arena is one-dimensional with length \(R=L/\ell\) centered on the origin \(x=0\), with parameters for speed \(v_\infty\) and distance \(\ell\) defined in Eq. 16. The position and velocity are *x* and *v* have dimensionless counterparts \(r=x/\ell\) and \(u=v/v_\infty\). The force *F* and its dimensionless version *f* are

and the rate equation that determines the equilibrium distribution \(\Pi (r,u)\) is

Following the protocol for the simulations described above, we know that the boundary distribution \(\Pi (R/2,u)=\Pi (-R/2,-u)\) should be similar to \(\delta (u)\) because all members that reach a boundary are restarted with zero velocity. However, at a boundary, the distribution is not an even function of *u*.

To solve Eq. 24 it is helpful first to find eigenfunctions \(g_\lambda (u)\) of the velocity operator appearing in it, in the sense that

for a given eigenvalue \(\lambda\). (The primes denote derivatives with respect to *u*.) Having found such functions, we can write the general solution for the equilibrium distribution as a sum over \(g_\lambda (u)\) functions weighted by coefficients \(c_\lambda (r)\) that vary exponentially with position,

Also, Eq. 25 implies \(g_{-\lambda }(u)=g_\lambda (-u)\) and the equilibrium distribution must satisfy reflection symmetry so \(\Pi (r,u)=\Pi (-r,-u)\). This means we can restrict the sum to non-negative eigenvalues and write

The \(\lambda =0\) term is position-independent: it gives the equilibrium distribution for an unbounded arena. The \(\lambda >0\) terms vary with position: the first term on the right of Eq. 27 dominates near the right boundary but is negligible at the left boundary. The second term has the opposite property.

The eigenvalue equation, Eq. 25, for \(\lambda =0\) requires that

where the last equality is required for a (normalizable) physical distribution. It is solved by

Clearly, \(g_0(u)\) is the 1D analogue of the equilibrium distribution found above, in Eq. 19, for the 2D unbounded arena. According to Eq. 27, the distribution approaches \(g_0(u)\) as one moves further from a boundary.

The eigenfunction \(g_\lambda (u)\) for \(\lambda \ne 0\) can be expressed exactly as a Hermite function, with an order that increases with \(\lambda\), multiplied by a Gaussian factor. It oscillates increasingly rapidly as \(\lambda\) increases, so a \(g_\lambda (u)\) function is analogous to a Fourier component. Near the boundary, many \(g_\lambda (u)\) contribute and their sum gives a velocity distribution that is sharply peaked near \(u=0\).

To make it easier to sum over \(g_\lambda (u)\) with \(\lambda >0\) we approximate these functions. Consider Eq. 25 when \(\lambda>>1\) and the terms expected to dominate it satisfy \(g''_\lambda (u)-\lambda u g_\lambda (u)=0\). This is solved by an Airy function of the first kind: \(g_\lambda (u)\propto \textrm{Ai}(\lambda ^{1/3}u)\). (An Airy function of the second kind also solves this equation but would give a physically incorrect distribution.) We substitute a trial form \(g_\lambda (u)=\textrm{Ai}(\lambda ^{1/3}u)\exp (\chi _1[u])\), with \(\chi _1(u)\) an unknown function, into Eq. 25 and obtain

For sufficiently large \(\lambda\) the first expression in parentheses, which is multiplied by \(\lambda ^{1/3}\), must vanish. It does if \(\chi _1(u)=\chi _0(u)/2\). The result

agrees reasonably well with the true eigenfunction for \(\lambda \approx 1\) and becomes exact in the limit of large \(\lambda\). It is convenient because it is well behaved for any eigenvalue so we can treat \(\lambda >0\) as a continuous parameter. Note that \(\textrm{Ai}(0)=(3^{2/3}\Gamma [2/3])^{-1}\approx 0.36\) and the Airy function is positive and exponentially decaying when its argument is positive but oscillates in sign for negative argument.

We replace the sum over non-zero \(\lambda\) with an integral using an interpolating function \(c(\lambda )\equiv c_\lambda (R/2)\). A minimum value of \(\lambda =1\) will be imposed on the continuous spectrum–the precise minimum value will not affect important results below. The distribution becomes:

The cutoff of the \(\lambda\) integral can be chosen according to the desired position resolution \(\Delta r\) and velocity resolution \(\Delta u\) of the distribution. The terms with \(\lambda \gg 1/\Delta r\) are zero at any measurable distance from the boundary, so they can be ignored. Also, since \(\textrm{Ai}(\pm \infty )=0\), terms with \(\lambda>>1/\Delta u\) contribute nothing to the distribution. For our, qualitative, purpose we set \(\Lambda =\infty\) whenever the integral converges. Also, if the \(\lambda =0\) term does not dominate the sum then we can extend the lower integration limit to \(\lambda =0\) (it does not matter if we use Eq. 31 to poorly approximate \(g_0(u)\) in this case).

The \(c(\lambda )\) coefficients should be determined by matching the \(r=R/2\) distribution

to a known boundary value. Current conservation at the boundary determines the rate that members are replaced at \(r=R/2\) with zero velocity. But as a rough approximation to this procedure, and in analogy with the corresponding situation for a Fourier transform, we assume that \(u=0\) members at \(r=R/2\) are uniformly distributed among all \(\lambda\) components. This means \(c(\lambda )=c_1\approx c_0\) with \(c_1\) a constant.

The border distribution is now

where we used the fact that the integral is strongly peaked about \(u=0\) to set \(\chi _0(u)\approx \chi _0(0)=0\) and then used the defining differential equation of the Airy function to do the \(\lambda\) integral. At \(u=0\) the function has a value \(\Pi (R/2,0)\approx \Lambda\) and is peaked at the slightly negative velocity of \(u=u_{\max }\approx -1.47\Lambda ^{-1/3}\) with \(\Pi (R/2,u_{\max })\approx 1.46\Lambda\). For positive velocities it drops to zero on a scale of \(u_{\max }\) while for negative velocities it oscillates on this same scale, rapidly averaging to zero. In Fig. 6 we plot it using a modest value \(\Lambda ^{1/3}=20\) so the qualitative behavior can be clearly seen.

The overall velocity distribution is similarly obtained

where

The function *B*(*u*) comes from the population near the boundaries and determines the velocity distribution at small *u*. It has a \(u=0\) peak of height \(B(0)=2\ln \Lambda\) and width (FWHM) of approximately \(4/\ln \Lambda\). While *B*(*u*) depends on \(\Lambda\), this dependence is weak. In Fig. 6 we plot *B*(*u*) and *P*(*u*) with \(\Lambda =10^6\) and \(u_0=1.5\), \(R=4\) and \(c_0=c_1\). (The vertical axis for plots of *P*(*u*) and *n*(*r*) is arbitrary in this figure: the amplitude of these distributions should be fixed by normalization.) The cusp in *P*(*u*) at \(u=0\) comes from that of the *B*(*u*) function while the shoulders, coming from \(g_0(u)\), correspond to \(u=\pm u_0\).

To obtain the density

we integrate Eq. 32 numerically. The \(\lambda\) integral converges so we can take \(\Lambda \rightarrow \infty\). The result for \(n(r)-n(0)\) is plotted in Fig. 6 versus the distance from either boundary.

Insight can be gained by considering the case of no driving force, i.e. setting \(u_0=0\), so that analytic results for *n*(*r*) are possible. Integrating Eq. 32 over *u*, using the large \(\lambda\) approximation of \(2\sqrt{\pi }\textrm{e}^{2\lambda ^2/3}\textrm{Ai}(\lambda ^{4/3})\approx \lambda ^{-1/3}\), and taking \(\Lambda \rightarrow \infty\) we find

Writing the positive distance from the left boundary as \(d=R/2+r\) or from the right boundary as \(d=R/2-r\), we have

If we had extended the lower limit of the integral to \(\lambda =0\) then the density, for all distances *d*, would be given by the \(d<<1\) expression above (so the lower cutoff in the \(\lambda\) integral affects the long-distance decay of the boundary population). In Fig. 6 we compare the \(d^{-2/3}\) divergence with the numerical calculation of *n*(*r*). To within a constant multiplicative factor, the \(u_0=0\) approximation of the density is similar as the result for finite \(u_0\).

The results above are consistent with what we saw in the experiment and the model simulations. The density *n*(*x*) is large at the boundary because so many slow-moving ants remain near it. The density decreases with distance *d* and ants further than \(\ell\) from the boundary act as though they are in an infinite arena. The sharp \(v=0\) peak in the velocity distribution *P*(*v*) is due to ants at the boundary that are often abruptly stopped. The shoulder features in *P*(*v*) are properties of the distribution in the arena interior: they are local maxima in the distribution centered on \(v=v_0\) with a width of \(v_\infty\).

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

## About this article

### Cite this article

Charoonratana, L., Thiwatwaranikul, T., Paisanpan, P. *et al.* Modeling the movement of *Oecophylla smaragdina* on short-length scales in an unfamiliar environment.
*Mov Ecol* **11**, 64 (2023). https://doi.org/10.1186/s40462-023-00426-w

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/s40462-023-00426-w