The food seeking behavior of slime mold: a macroscopic approach

Starting from a particle model we derive a macroscopic aggregation-diffusion equation for the evolution of slime mold under the assumption of propagation of chaos in the large particle limit. We analyze properties of the macroscopic model in the stationary case and study the behavior of the slime mold between food sources. The efficient numerical simulation of the aggregation-diffusion equation allows for a detailed analysis of the interplay between the different regimes drift, interaction and diffusion.


Introduction and Background
Physarum Polycephalum, or a slime mold, is an amoeboid organism that is notable for its ability to perform complex tasks despite its relatively simple biological strucutre and lack of a brain. Surprisingly, slime molds are known to be able to solve mazes, construct robust networks, and solve shortest path and spanning tree problems [1-3, 5, 6, 16, 19, 23, 24]. The simplicity of the slime mold organism and the complexity of its emergent behaviors have inspired research into Physarum Polycephalum as a model for collective behavior [13][14][15]25]. A hallmark of research in collective dynamics is complex global structure emerging as a result of relatively simple local interactions between agents. Indeed, cellular automata models inspired by physarum polycephalum food seeking behavior have been found to be good models of other collective phenomena such as the formation of transportation networks in a country or the "cosmic web" of stellar material between galaxies [4,7].
When seeking food, physarum polycephalum moves to concentrate most of its mass on food sources -however it maintains a network of mass between food sources. In the absence of food it begins to retract to a central mass. Past research into models of physarum polycephalum mainly attempt to model this food-seeking behavior and have taken several different directions. Many are discrete cellular automata models -the slime mold is represented by a collection of agents that evolve discretely in time and space. Different interaction rules and sometimes even different "types" of agents are employed to model different slime mold behaviors however almost all models include a rule that causes slime mold agents to move towards food sources [12,17,18,28,29]. Other approaches note that physarum polycephalum moves and forms networks by moving fluid through tubules within its mass -much like an electrical network [25][26][27]. These approaches focus on modeling physarum polycephalum as a network of ordinary differential equations and attempt to model the network that physarum polycephalum builds between food sources once they are found. In particular they do not attempt to model the food seeking behavior of the organism.
In this manuscript we propose a model hierarchy aimed at modeling the food seeking behavior of Physarum Polycephalum. The main behaviors we attempt to model are aggregation of mass on food sources while maintaining connected paths between food sources and retraction in the absence of food sources. We first define an agent based model that is continuous in time and space through a system of stochastic differential equations. The evolution of each agent is governed by three terms -an interaction term with other agents that causes agents to aggregate, a drift term that causes agents to move towards food sources and a noise term intended to model slime mold foraging behavior in the absence of food sources. We then rigorously show under the assumption of propagation of chaos that as the number of agents approaches infinity, the agent based model converges in a sense to a deterministic density based model in the form of an aggregation-diffusion equation. This macroscopic model includes terms that correspond to the drift, interaction and noise terms present in the microscopic model. We then embark on an analysis of the macroscopic equation from the point of view of modeling Physarum food seeking behavior. We first discuss a framework for choosing an appropriate interaction kernel through an analysis of the steady states of the macroscopic model in the case that there are no food sources. Then, we rigorously show that there is a parameter regime in which the only possible steady state of the macroscopic model is zero and thus this regime is not suitable for modeling Physarum food seeking behavior. Finally, we simulate the macroscopic model using the blob method for aggregation-diffusion equations presented in [8] and find that it is necessary to scale the three main terms of the model in order to reproduce aggregation on food sources. We then simulate the evolution of the model in different scaling regimes and with different choices of interaction kernel and qualitatively examine which regimes result in the main behaviors we attempt to model. We find that there are several regimes that result in the qualitative desired behaviors.

Model Derivation
The main features of slime mold food seeking behavior are aggregation of mass on food sources while maintaining a connected mass and retraction in the absence of food sources. To model these behaviors we first introduce an agent-based model. Here, the slime mold is represented by a swarm of N agents who's trajectories, X i (t) ∈ R d , are modeled as a system of stochastic differential equations. In the following we will refer to this model as the microscopic model in order to reinforce the intuition that the agents represent slime mold "particles".
Model 1 (Microscopic Model) Consider a collection of N agents with (random) positions (X 1 (t), ..., X n (t)) ∈ (R d ) N . The microscopic model is given by the evolution: We note that if we assume that ∇V and ∇W are globally Lipschitz and further if there exists then it can be shown that (1) has a unique strong solution. Intuitively, V represents the density of chemo-attractants given off by food sources. For this reason, we usually think of V as a sum of radially symmetric, positive functions. Each term represents the attractants emanating from a single food source. A prototypical example would be a sum of Gaussians centered at different positions. W represents the interaction between agents. Notice that, roughly, if ∇W (X j (t) − X i (t)) ≈ X j (t) − X i (t), then the interaction causes agents to attract each other. lLikewise if ∇W (X j (t) − X i (t)) ≈ X i (t) − X j (t), then the agents will repel each other. As we have noted, Physarum tends to maintain a connected mass, even when aggregating around disparate food sources. Therefore, this interaction term should cause agents to remain locally close, i.e. agents should not repel each other at long ranges. The noise term, CdB i (t), models physarum foraging behavior. In the absence of chemo-attractants, Physarum forages by diffusing outward from its initial position while still maintaining a connected mass. Here, in the absence of a food source (V = 0) the noise term should cause particles to spread from their initial positions while the interaction term causes them to still have some propensity to remain locally "together".
We include the parameters A, B and C for the case that simulations of (1) demonstrate the need to scale the terms' influence on each particles trajectory to produce evolutions that resemble Physarum. However, any simulation of (1) that would give an accurate picture of how (1) models Physarum would have to use a very large number of particles; this becomes computationally intractable. Therefore, we must find a way to examine the behavior or (1), for a large number of particles, without computing the trajectory of every particle explicitly. We will show that in the large particle limit that the evolution of the marginal distribution of any particle is given by an aggregation-diffusion equation. There exist computationally efficient methods for simulating such equations, therefore we can examine the behavior of (1) by simulating the evolution of the marginal density. The technique we will employ to derive an appropriate model is the assumption of propagation of chaos; i.e. if the number of agents agents is large, the trajectories of any two agents can be assumed to be independent. This is true in many cases, see for example [9,10]. Theorem 1 In the N → ∞ limit and for a given, A, B, C > 0, W, V ∈ C 1 (R d ) and ρ 0 the evolution of the marginal distribution of any agent evolving according to (1) is a solution of the aggregationdiffusion equation under the assumption of propagation of chaos.

Proof
Consider a finite group of N agents, (X i (t), t ≥ 0)i N ⊆ R d , i = 1, ..., N evolving according to (1). Let µ N be the empirical distribution of the configuration of agents. That is: where δ X i (t) is the Dirac measure with unit mass at X i (t). Let F be a test function in C ∞ c ((R d ) N ), therefore by definition of the Dirac measure we deduce: we first derive an expression for the time evolution of the expectation of this quantity. Let X(t) = (X 1 (t), ..., X N (t)). By conditioning we find where ρ N (x, t) is the joint distribution of the collection of agents. Using Ito's Lemma (see e.g. [30]) we can further expand the conditioned term in the above We first focus on the third term of (3). Notice that as each X i evolves according to (1), since the noise terms for disparate agents are independent, we have that: Therefore, the third term of (3) simplifies as: (4) Now, turning to the second term of (3) we find by appyling (1) that: As the stochastic integral is a martingale, we have that the last term of (5) is equal to 0. Therefore, combining (4) and (5), we have that: (6) Therefore, by (2) and the continuity of the trajectories, (6) implies P almost-surely that: Therefore, ρ N is a weak solution to the following initial value problem: Using (8), we now compute the time evolution of the marginal distribution of one agent in the N → ∞ limit under the assumption of propagation of chaos, that is for any finite collection of M agents we have that their joint distribution satisfies: Without loss of generality we will consider the marginal distribution of the first agent which is given by: and therefore, by (8), we have that ρ(x 1 , t) satisfies (weakly): We analyze (10) term by term. First, by applying (9) we find that Next, if we again apply (9) and integrate by parts in each term of the sum where j = 1 we find that Similarly we deduce We can simplify (13) further by integrating under the divergence operator and applying (9), notice that: Therefore: Therefore, combining (11), (12) and (15) we have by (10) that the marginal of any agent satisfies: as desired.
As the equation derived in Theorem 1 describes how the entire collection of particles evolves from a density standpoint, in the following we will refer to it as the macroscopic model.

Model 2 (Macroscopic model)
The macroscopic model is given by the evolution of the aggregationdiffusion equation Here, V and W continue to represent the chemoattractant of food sources and the slime mold's propensity to aggregate respectively. In the following we will study (17) from the viewpoint of using it to model Physarum food-seeking behavior.

Model Properties
Many questions remain about how to utilize the macroscopic equation derived in the previous section to model a slime mold. In this section we present some analysis of stationary states of the equation in order to inform some of these modeling choices. First, we investigate stationary states in the case where there is no food present as a way to gain some information about what choice of kernel function might be reasonable. Inspired by [20], we believe that Gaussians represent a reasonable model for stationary states of a slime mold -however we find that Gaussian stationary states are only possible in the case that the kernel function, W , is quadratic. Next, we investigate stationary states of the scaled equation. We find through a fixed point argument that if the diffusion scaling parameter, C, is sufficiently higher than A and B that the only stationary state of the equation is 0 (even in the presence of a food source). This suggests that when tuning the scaling parameters that there cannot be "too much" diffusion.

Zero-food stationary states
In the case that there are no food sources a slime mold will "retract" to a more compact configuration [21][22][23]. Often, this configuration is roughly radially symmetric. Given that, a reasonable model for the configuration of the slime mold in the case that there is no food could be a Gaussian. We now investigate if there are conditions on the kernel function that are imposed by the assumption that stationary states are Gaussian. For simplicity we work in one dimension however the calculations are analagous in higher dimensions. In the case that there are no food sources (17) becomes: Therefore a stationary state satisfies (in the following we will ignore the scaling parameters, A and B, as they do not change the computation): and finally that: (19) represents a general condition that zero food stationary states must satisfy. If we assume that stationary states are Gaussian, i.e. that: then by (19) we have that: The above can only hold if W (y) = ay for some a ∈ R. Therefore, Gaussians are zero food stationary states only if the interaction kernel is quadratic. Further, by plugging in ay for W (y) in the above we can see that the following must hold: So, for a given quadratic kernel there are a family of Gaussians that are potential zero food stationary states (depending on the mass of the initial profile). For an illustration of this fact see Figure 1.
Here, we simulate (17) in the case of no drift term using the "blob method" for aggregation-diffusion equations introduced in [8]. We will discuss the numerical method in more detail in the Numerics section of the paper. We choose an initial condition representative of a slime mold agregated around two food sources at x = 1 and x = −1. We find in accordance with our calculation above that the mass profile appears to converge to a Gaussian configuration. This also illustrates a modeling property of (17); in the case that food sources "run out" the profile will "retract" to a central configuration.
Effectively, by starting from information about the zero food stationary state, we have "solved" for the correct kernel. This analysis suggests an empirical method for determining the correct kernel function from experimental data. Instead of our assumption that zero food stationary states are Gaussian, statistical analysis of slime mold configurations in the presence of zero food sources could provide information about the stationary state configuration. This information could then be used to "solve" for the kernel that results in this "correct" stationary state.

High diffusion stationary states
We now turn to considering stationary states of the scaled equation: We will see in a later section that the "fair" regime (A = B = C) does not result in food-source aggregation. Therefore we are motivated to study different scalings of the terms in (17). Here, we will rule out a large regime of scaling parameters by employing a fixed point argument to show that if the diffusion parameter, C, is sufficiently larger then the drift and interaction parameters then the only stationary state of (17) is 0 on bounded domains.

Theorem 2
Let Ω ⊆ R d be compact and connected. For a given A, B < 0 and V, W inH 2 (Ω) ∩ L ∞ (Ω) then for any sufficiently large C > 0, 0 is the unique stationary state of:

Proof
Recall that ρ is a stationary state of (17) if: We define: Given ρ ∈ H 1 0 (Ω) we know that G(ρ) ∈ L 2 (Ω) and that the equation: has a unique solution in H 1 0 (Ω). Given K > 0, define H 1 0 (Ω, K) := { ρ ∈ H 1 0 (Ω) : | |ρ|| H 1 0 ≤ K }. It can be shown (see, for example, Theorem 4 in section 6.3 of [11]) that ifρ is the unique solution of (21) then the following estimate holds for some D > 0: Therefore, by applying standard estimates we can see that for a sufficiently large choice of C if ρ ∈ H 1 0 (Ω, K) then we must also have that the solution to (21),ρ, belongs to ρ ∈ H 1 0 (Ω, K) as well. DefineÃ : H 1 0 (Ω, K) → H 1 0 (Ω, K) via: We will show thatÃ is a contraction and thus has a unique fixed point by Banach's fixed point theorem. Since 0 is trivially a fixed point ofÃ and K is arbitrary this implies that 0 is the only fixed point ofÃ and thus the only stationary state of (17) .
We will show that there exists a 0 < γ(A, B, C) < 1 such that: Since Ω is bounded by assumption we have by Poincare's inequality that: for some α > 0. Therefore sinceρ 1 andρ 2 are both solutions of (21) we must have by another application of the Poincare inequality that: Therefore we have that: and by (23) that: We claim that for some γ > 0 that depends only on A, B and C for a given W and V . Expanding we find that: We will continue our analysis term by term. Starting with the first term and applying the Cauchy-Schwarz inequality we find that (27) Now, turning to the second term of (26) and applying the Cauchy-Schwarz inequality we find: We will further expand the last term of (28). Notice that: Now, combining (28) and (29) we find Finally, we turn to the third term of (26). Notice that: Now, combining our results in (27), (30) and (31) we find that: where γ ≥ 0 and is independent of ρ 1 and ρ 2 as we have assumed a uniform bound on their norms (they are members of H(Ω, K)). Therefore, if W and V are given, γ can be made arbitrarily small if C is sufficiently larger than A and B. Additionally we have that for some β > 0 and therefore Combining (24) and (33) we obtain: . Therefore for C sufficiently larger than A and B, the mappingÃ is a contraction as desired.
From a modeling point of view, Theorem 2 shows that too strong a contribution from the diffusion term will destroy any of the qualitative behaviors we wish to capture. We will see this explicitly in our numerical studies of (17) in the next section.

Numerics -Qualitative features of the macroscopic model
We now present some numerical investigations of the macroscopic model derived in the previous section.
These investigations are aimed at replicating qualitative properties of slime mold movement seen in [21][22][23], specifically aggregation around food sources while still maintaining a "connected" mass. Recall that here, the function V models the density of chemoattractants dispersed by food sources. Therefore, V will always take the form of a sum of radially symmetric positive functions where each term is centered on a food source. We simulate (34) via the "blob method for aggregation-diffusion equations" introduced in [8] which approximates (34) by solving the ODE: for a collection of N particles, (x 1 (t), ..., x N (t), t ≥ 0), who's initial positions are a regular grid on the domain on which we'd like to approximate (34) and who's masses are given by m i = ρ(0, X i (0)). To recover an approximation of ρ from the positions of the particles we convolve the particle solution with a mollifier ϕ :ρ We assume mollifiers are always of the form: where d is the dimension of the domain. It is shown in [8] that under some regularity and growth conditions on W and V thatρ converges to the solution of (34). As a first example we will consider (34) in one dimension. We will take the food source to be Therefore we can think of the food sources as being located at x = −1 and x = 1. For the parameters of the ODE we take: • initial mass profile given by ρ 0 (x) = 1 √ 2πσ 2 e − x 2 2σ 2 , σ 2 = 0.0625 as in [8].
Note that as the variance of the initial profile is very small that, at the particle level, slime mold particles have a very low probability of being found outside of [−2.1, 2.1] (initially). In Figure 2 we perform a first simulation of (17) in one dimension via the blob method in the so called "fair regime" (A = B = C = 1). Here, we choose a quadratic interaction kernel -specifically we impose that W (x) = x. Qualitatively, we do not observe aggregation of slime mold mass around food sources to the degree present in [23]. Therefore, motivated by this example, in the following we will attempt to modify (34) in order to better model slime mold food seeking behavior. In Figure 2 we observe that diffusion appears to be the dominating effect. Therefore, in order to examine the interplay of the diffusion term against the other two terms in (34) we will employ the scaling parameters, A, B, C > 0, in order to control the contributions from the drift term, interaction term and diffusion term respectively. Roughly we examine two main regimes; "drift dominated" (where A > B, C) and "interaction dominated" (where B > A, C). We also examine so called "competition regimes" where two scaling parameters are equal but dominate the third. We do not examine the diffusion dominated regime as we've previously seen that if A = B = C (we will refer to this as the "fair regime") then diffusion is the dominating effect. In each regime we examine three different interaction kernelsquadratic, polynomial and Gaussian. We choose the same ODE parameters as in the previous fair regime simulation.

Food source dominated regime
We will first examine simulations of (17) in the regime where the drift term (which models attraction to the food sources) is the dominating factor. We first examine the case where the kernel is quadratic Here, the kernel is purely attractive. At the agent based level all agents exert a "pull" on other agents with a force that is proportional to their spacial seperation. The results can be seen in Figure 3. Here, we qualitatively observe that in the case where A = 10, B = C = 1 that we have strong food source aggregation however, the mass profile splits into two distinct bumps -the strong contribution from the drift term prevents a connected mass from being maintained. In the other three cases a connected mass is maintained and food source aggregation is observed. However, in the two cases where C = 5 the degree of aggregation is relatively weak likely due to the stronger contribution from the diffusion term.
Next, we examine the case where the kernel is polynomial and given by W (x) = x 4 4 − x 2 2 . Intuitively, at the particle level, this means that this kernel is attractive at long-ranges but includes a short-range repulsion to prevent particles from becoming too close. The results can be seen in Figure 4. Here, we observe results that are qualitatively very similar to the results obtained in the quadratic case. This is not surprising given that the polynomial kernel is largely attractive and at long ranges the degree of attraction is proportional to the separation between particles -analagous to the quadratic kernel. However in the case where A = 10, B = 5, C = 1 the evolution with a polynomial kernel does not maintain a connected mass in contrast to the evolution witha quadratic kernel. This is likely due to the short-range repulsion that is present in the polynomial kernel. At the particle level, particles in the middle of the domain are repelled away from eachother while simultaneously being pulled by a food-source -encouraging a separation.
Finally, we examine the case where the kernel is Gaussian and given by Here, similar to the quadratic kernel the interaction is purely attractive. However, unlike the polynomial the strength of attraction between two particles is not proportional to their separation. Instead, attraction becomes stronger (but not in an unbounded fashion) as two particles become closer and weaker as they become farther apart. This encodes the modeling assumption that interactions should be "local" -particles who are physically distant from eachother should not interact strongly. The results can be seen in Figure 5.
Here, we observe food source aggregation in every case. Similar to the polynomial and wuadratic   Food source aggregation occurs in all cases, however in two cases a connected mass is not maintained. kernels the aggregation is stronger when the contribution from the diffusion term is small. However, similar to the polynomial kernel a connected mass is not maintained in either of these cases. Additionally, in the cases where a connected mass is maintained we still observe a sharp drop off in the middle of the profile -indicating formation of a separation. We beleive that the tendency to separate in this case is caused by a combination of the strong local interaction and the strong contribution from the food sources; particles that start close to a food source are "trapped" and remain close to the food source as the contribution from the food source and the local interaction reinforce eachother. Then, particles towards the center of the domain feel a pull towards food sources form both the food source itself and the "trapped" particles.

Interaction dominated regime
We now examine simulations of (17) in the regime where the interaction term is the strongest contribution. Again, we start with the case where the kernel function W is quadratic and given by W (x) = x 2 2 . The results can be seen in Figure 6. Here, we do not see any food source aggregation. This can possibly be explained by the fact that in these simulations the initial profile is Gaussian. We have previously seen that Gaussians are stationary states for (17) in the case of a quadratic kernel and no drift term.
Next we examine the case where the kernel function is polynomial, of the same form as in the food source dominated regime. The results can be seen in Figure 7. Here, analagous to the food source dominated regime we see more pronounced aggregation in the cases where the diffusion term has the weakest contribution. However the aggregation is less pronounced than the respective examples in the food source dominated regime due to the stronger contribution from the interaction term. Unlike the analagous cases in the drift dominated regime we do not observe food source aggregation in the cases where C = 5, this is possibly due to the stronger contribution from the short range repulsion in the interaction term.   Finally, we examine the evolution in the case of a Gaussian kernel function. The results can be seen in Figure 8. Here, in contrast to the respective cases in the drift dominated regime, we do not observe food source aggregation in the cases where the diffusion term has the second largest contribution. Addtionally, unlike the analagous cases in the drift dominated regime and the polynomial kernel in the interaction dominated regime, we do not observe aggregation in the case where B = 10, C = 1, A = 1. We do observe aggregation in the case where A = 5, B = 10, C = 1 however similar to the polynomial kernel the effect is less pronounced than in the analagous case in the drift dominated regime. Both of these effects are consistent with a stronger contribution from the attractive kernel.

Competition regime
We now simulate cases where two parameters dominate in competition. The results can be seen in Figure 9, Figure 10, and Figure 11. Interestingly, across all three interaction kernels, pronounced food source aggregation is only observed in cases where A = B = 10. All of these cases except one exhibit the qualitative features we wish to capture as they do maintain a connected mass. All other cases do not result in food source aggregation. This can be explained by the fact that in all of these cases the diffusion term has a significant contribution. The observation that strong contributions from the diffusion term dampen food source agregation is a general pattern that we have also observed in the other regimes and is additionally supported by the statement of Theorem 2.
The cases in the competition regime where we do see food source aggregation also reinforce another pattern observed across all regimes; examples that have the qualitative features we wish to capture have strong contributions from the drift term, the interaction term, or both. However, we also find that cases in which a connected mass was not maintained only occured where the drift term was dominant. Likewise, cases in which food source aggregation did not occur in the absence of strong diffusion only occured in the interaction dominated regime. These observations suggest a general modeling principal; contributions from the drift and diffusion terms should be scaled similarly and larger than the contribution from the diffusion term.

Conclusion
In this manuscript we presented a model hierarchy aimed at modeling the food seeking behavior of the slime mold Physarum Polycephalum. We were principally interested in capturing the slime mold's ability to aggregate around disparate food sources while maintaining a connected mass. We first presented a particle based model which includes three main features -a drift term to model a gradient of chemoattractant produced by food sources, an interaction term to model the slime mold's propensity to maintain a connected mass, and a diffusion term to model slime mold foraging behavior. Additionally, we included scaling parameters for each of these three terms in order to study how varying their contributions could aid in modeling Physarum. However simulating enough particles to realistically model the evolution of slime mold would be computationally intractable. Therefore, under the assumption of propagation of chaos we showed that in the large particle limit the evolution of the particle model can be described by a macroscopic aggregation-diffusion equation which can be efficiently simulated. Before embarking on simulating the equation we made some analytical observations via analysis of the equation's stationary states. First, we found that in the case of no drift term that an assumption of Gaussian stationary states implies that the interaction kernel must be quadratic and discussed how this could be used as a heuristic to determine a kernel that realistically models physarum behavior. Then, we showed that if the contribution from the diffusion term is sufficiently larger than the interaction and drift terms that the only possible stationary state is 0, deeming this regime as unsuitable for modeling Physarum.
We then ran a series of simulations of the macroscopic model in order to investigate how to scale the three contributing terms and choose the interaction kernel in order to model Physarum food   (17) in the competition regime with kernel given by W (x) = seeking behavior. We simulated three different parameter regimes; a regime in which the drift term dominated, a regime in which the interaction term dominated and a regime in which the dominating terms were in competition. In each regime we examined three different interaction kernels, an attractive kernel with a local interaction, an attractive kernel with long range interaction, and a kernel that included a short range repulsion. We observed cases that reproduced connected mass food source aggregation across all three regimes with all kernels however several general patterns did emerge. We found that strong contributions from the diffusion term tend to dampen food source aggregation even in the presence of strong contributions from the other terms. Additionally, we found that in the absence of strong diffusion, dominance of the drift term can result in mass separation while dominance of the interaction term can prevent food aggregation. This suggests a general modeling strategy of scaling contributions from the drift and interaction terms similarly and larger than the diffusion term.
There are many ways in which this study could be extended. Aggregation-diffusion equations have been widely studied in the case of no drift term [9]. A main finding is that there is a critical mass above which the interaction term becomes dominant and causes a "blow up" in the mass profile. From the particle view this can be thought of intuitvely as all particles converging on one point. Below the critical mass, the tendency of the diffusion term to cause the mass profile to spread prevents a blow up from occuring. We did not observe this phenomena in any of our simulations (even in the interaction dominated regime). It would be interesting to consider whether the presence of a drift term can prevent blow up from occuring at all or mitigate it in the sense that the critical mass at which blow up occurs becomes higher. Additionally, as we did observe the desired qualitative behavior in several cases it would be useful to compare the evolution of our model to data from actual slime mold growth in a simple situation such as the two food source case we consider. This would be an important first step in the quantitative validation of this model.