Next Article in Journal
Entropy Generation Results of Convenience But without Purposeful Analysis and Due Comprehension—Guidelines for Authors
Next Article in Special Issue
Non-Extensive Entropic Distance Based on Diffusion: Restrictions on Parameters in Entropy Formulae
Previous Article in Journal
Information and the Quantum World
Previous Article in Special Issue
Energy Flows in Low-Entropy Complex Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interacting Brownian Swarms: Some Analytical Results

STI/IMT/LPM, Ecole Polytechnique Fédérale de Lausanne, Station 17 (Batiment BM), CH-1015 Lausanne, Switzerland
*
Authors to whom correspondence should be addressed.
Submission received: 18 November 2015 / Revised: 7 January 2016 / Accepted: 11 January 2016 / Published: 14 January 2016

Abstract

:
We consider the dynamics of swarms of scalar Brownian agents subject to local imitation mechanisms implemented using mutual rank-based interactions. For appropriate values of the underlying control parameters, the swarm propagates tightly and the distances separating successive agents are iid exponential random variables. Implicitly, the implementation of rank-based mutual interactions, requires that agents have infinite interaction ranges. Using the probabilistic size of the swarm’s support, we analytically estimate the critical interaction range below that flocked swarms cannot survive. In the second part of the paper, we consider the interactions between two flocked swarms of Brownian agents with finite interaction ranges. Both swarms travel with different barycentric velocities, and agents from both swarms indifferently interact with each other. For appropriate initial configurations, both swarms eventually collide (i.e., all agents interact). Depending on the values of the control parameters, one of the following patterns emerges after collision: (i) Both swarms remain essentially flocked, or (ii) the swarms become ultimately quasi-free and recover their nominal barycentric speeds. We derive a set of analytical flocking conditions based on the generalized rank-based Brownian motion. An extensive set of numerical simulations corroborates our analytical findings.

1. Introduction

Elementary models of multi-agents swarms often rely on the assumption that all individuals obey identical dynamical rules. This hypothesis of agents’ homogeneity significantly simplifies the mathematical modeling and its treatment. In particular, for large swarms, homogeneity enables analytical discussions based on the mean-field (MF) approach. In the MF approach, the effects of mutual interactions are aggregated into an external effective force field, driving the dynamic of a single randomly picked agent, whose behavior reflects the whole society’s evolution. Homogeneity is, however, rarely encountered in actual situations. A natural question is to then ask how heterogeneity between agents affects and/or possibly destroys the potential emergence of collective (i.e., flocking) dynamics. In the presence of heterogeneity, analytical approaches are more difficult as, a priori, the MF approach cannot be directly employed. In multi-agent swarms, heterogeneity can manifest itself in many different forms and degrees. It may concern the individual dynamics of each agent and/or the way agents interact with their fellows. An extreme heterogeneous situation is realized when, due to its specificities, a single agent is actually leading the entire swarm [1,2,3,4,5]. In this paper, we focus on comparatively low heterogeneities encountered when only the individual dynamics differ but not the interacting rules. For this case, a paradigmatic and pioneering example is the Kuramoto synchronization of heterogeneous phase oscillators for which the role of heterogeneity can, to a certain extent, be analytically discussed [6,7]. More recently, models based on rank-based drifted Brownian motion (RBM) [8], and, in particular, the so-called Hybrid Atlas Model (HAM) [9], offer another possibility to analytically discuss the heterogeneous swarm dynamics, without relying on the MF approach. In the sequel, we shall basically use the HAM modeling framework to study a new type of global heterogeneity. This is achieved when two sub-swarms are forced to mutually interact. Each sub-swarm, formed by homogeneous agents, is initially tightly flocked and obeys specific collective dynamics. The global society formed by the union of the two sub-swarms is clearly heterogeneous. This example depicts a situation where two distinct sub-swarms (agents of each sub-swarm are different) interact via a local competition mechanism (here, agents systematically try to catch their leaders in a finite observation range). In this paper, we are able to analytically characterize the post-collisional outputs. Depending on external control parameters (individual sub-swarm dynamics, strength of competition), we observe that the post-collisional situation can either be a single tight swarm or a segregation into the nominal sub-swarm.
The paper is organized as follows: In Section 2, we introduce the microscopic agent’s dynamic and link it to RBMs. In Section 3, relying on the HAM theory, we derive the flocking conditions for single swarms composed of scalar agents interacting via imitation. The mutual interactions are effectively of long-range types (i.e., each agent can observe all leaders independent of their distance). For a finite number of tightly flocked agents, the (scalar) inter-distance between successive fellows is a stationary random variable, and, hence, so is the size of the swarm extension. This enables us to define a characteristic stationary swarm extension. In Section 4, we consider a couple of sub-swarms, initially well separated on the real line and with different average barycentric speeds. As the positive barycentric velocity of the leftmost sub-swarm is chosen, larger than the one of the rightmost sub-swarm, collision is unavoidable. For large but finite observation ranges of the sub-swarm (i.e., range larger than the sub-swarm characteristics of each sub-swarm extension derived in Section 3), we study the asymptotic post-collisional behavior of the whole system. The conditions for either the emergence of a global flocking swarm or quasi-freedom of the sub-swarms are approximatively derived. All analytical findings are tested and corroborated by simulation experiments. Concluding remarks, along with perspectives, are finally presented in Section 5.

2. Flocking of Interacting Brownian Agents

We consider a swarm of N Brownian agents A j , j = 1 , 2 , , N , diffusing on R . Agent A i follows the dynamic:
d X i ( t ) = d U , i ( X ( t ) , t ) d t + σ d W i ( t ) , X i ( 0 ) = x i , 0 ,
where X ( t ) = X 1 ( t ) , X 2 ( t ) , , X N ( t ) and d W i ( t ) ( i = 1 , . . , N ) are independent White Gaussian Noise (WGN) sources. The set of individual drifts d U , i ( X ( t ) , t ) for i = 1 , 2 , , N are rank-dependent interactions that are defined as follows:
(i)
In real-time, agent A i counts the (time-dependent) number N U , i ( t ) of leading fellows located within an observation interval U, namely the number of A k for k i that are located in R U , i ( t ) with:
R U , i ( t ) : = { x R 0 < x - X i ( t ) U } .
(ii)
The A i -rank-based drift in Equation (1) is then determined by:
d U , i ( X ( t ) , t ) : = α N U , i ( t ) N , i = 1 , 2 , , N ,
with α R + a velocity scale factor.
According to Equation (1), agent A i is positively accelerated by its leaders located in R U , i ( t ) but ignores the presence of its followers. These types of dynamics generically model imitation mechanisms; a typical example for the diffusion of innovation is shown in [10].
For very large swarms (i.e., when N ), we may define an agent probability density ρ ( x , t x 0 ) at position x at time t. Adopting a mean-field (MF) approach, ρ ( x , t x 0 ) obeys a (nonlinear) Fokker–Planck parabolic partial differential equation (pde):
t ρ ( x , t x 0 ) = F ρ ( x , t x 0 ) , F ( · ) : = σ 2 2 x , x ( · ) - α x 0 U ρ ( y + x , t x 0 ) d y ( · ) .
For U = , the model can be analytically investigated, and its solution behaves as a flocked swarm of agents that propagate like a soliton wave [10,11,12]. Conversely, for small U, a first order expansion of Equation (4) yields Burgers’ equation, that can be analytically solved via the well known Hopf–Cole (logarithmic) transformation [11]. In this small U regime, ρ ( x , t ) exhibits a diffusive behavior that asymptotically vanishes with time. Qualitatively, when U is very small, the agents mutual interactions are not strong enough to sustain a cooperative soliton dynamic pattern. By introducing a decreasing distance-dependent modulation on rank-based drifts, it has been showed that a bifurcation point, controlled by the distance decay rate, exists that separates both propagating regimes [10,12].
In general, however, the joint nonlocal and nonlinear characteristics of Equation (4) precludes direct analytical characterization of this bifurcation point. Alternative approaches then have to be found. In the present paper, focusing on Equation (1) with N < , we shall use available results for rank-based Brownian motions [8,9] to approximately estimate the observation range required to generate a soliton-like propagation.

Rank-Based Brownian Motions (RBMs)

For calculation convenience, we now consider the dynamics of Equation (1) with respect to the swarm’s barycenter. That is to say, Equation (3) becomes:
d U , i ( X ( t ) , t ) : = α N U , i ( t ) N - γ , γ : = 1 N l = 0 N - 1 α l N = α N - 1 2 N ( average   barycentric   drift ) .
In other words, agents A i , i = 1 , 2 , , N counts, in real time, the number of its leaders A k with k i . It then updates its drift d U , i ( X ( t ) , t ) as given in Equation (5), where γ is the average drift velocity of the swarm. Since the WGN in Equation (1) achieves unbounded realizations, the ability for all A i ’s to exactly count the number of their leaders effectively requires an infinite observation range U.
By construction, one can verify that the stability conditions:
k = 1 l α k - 1 N - γ = α 2 N l · l - N < 0 1 l < N
are verified, meaning that the laggards have a systematic average tendency to catch up with the leaders.
The rank-based drifts effectively implement a stylized gravity model, where the set of conditions for Equation (6) ensures the existence of a stationary probability measure for the RBM. This stationary regime corresponds to a tight propagating swarm that effectively behaves like a soliton propagating wave. Moreover, [8,13] established that in this stationary state, the inter-distances between successive agents are independent, λ l -exponentially distributed random variables with:
λ l = α N σ 2 l · N - l , l = 1 , 2 , , N - 1 .
Since the swarm’s stationary size is the sum of ( N - 1 ) independent exponentially distributed r.v.’s, it defines a hypoexponential probability distribution H ( x ) with density:
d H ( x ) = l = 1 N - 1 λ l e - λ l x l = 1 , k l N - 1 λ k λ k - λ l d x .
Accordingly, the average size of the swarm μ S ( N ) and its corresponding variance σ S ( N ) behave as follows:
μ S ( N ) = 0 x d H ( x ) = l = 1 N - 1 λ l - 1 = N σ 2 α l = 1 N - 1 1 l ( N - l ) N 2 σ 2 α log ( N ) , σ S 2 ( N ) = 0 ( x - μ S ) 2 d H ( x ) = l = 1 N - 1 λ l - 2 N σ 4 α 2 [ π 2 3 + 4 log ( N ) N ] N σ 4 π 2 3 α 2 .
The coefficient of variation c v S 2 : = σ S 2 ( N ) μ S 2 ( N ) , which measures the stochasticity of the swarm’s size, vanishes in the mean-field limit ( N ). Observe that, in the mean-field limit, the agents’ density ρ ( x , t ) of the RBM directly coincides with Burgers’ dynamics for U = . This leads to a tight soliton-like propagation density with infinite support. This is consistent with Equation (9), where the swarm’s size slowly diverges for N .

3. Heuristic Characterization of the Observation Threshold Leading to Cooperative Dynamics

Consider the nominal dynamics Equation (1) and now assume finite observation ranges (i.e., 0 < U < ). Due to the unboundedness of the WGN, an exact ranking procedure is therefore not strictly possible. A non vanishing probability (decreasing with increasing U’s) exists to find outlying leaders (i.e., located beyond the observation range of some of their followers). While the average swarm size μ S ( N ) from Equation (9) is valid only when U = , we postulate that it can still be used to approximately characterize the minimal observation range U c < required to sustain a tight swarm. Indeed, whenever U < U c , it is highly probable that some agents cannot determine their respective ranks. These outlying agents will underestimate their instantaneous drifts (from their incomplete observations), ultimately leading to their escape from the bulk. Outlying agents always escape the swarm from the rear, leading with time to an effective “evaporation”.
For consistency, the initial conditions X i ( 0 ) , to be drawn from a probability measure with finite support Ω 0 , namely X i , 0 Ω 0 with Ω 0 = [ l , r ] R and length ( Ω 0 ) = r - l μ S ( N ) must be restricted. Let us now define a characteristic relaxation time τ relax = μ S ( N ) σ 2 that corresponds to the time needed by the stochastic processes to diffuse over the interval of length μ S . For observation ranges U > μ S and initial conditions X i , 0 Ω 0 , all N agents mutually interact at time t = 0 , and, therefore, correctly determine their relative ranks. We assume that this will remain true with high probability, for times 0 < t < τ relax . In other words, we expect that the initially tight swarm is likely to remain flocked for t τ relax , effectively a meta-stable state. For asymptotic times, however (i.e., t τ relax ), the swarm’s tightness will be destroyed. Indeed, for any finite U, outliers will ultimately escape the attraction generated by the bulk of the remaining swarm. Therefore, our heuristic postulates that the characteristic observation range, separating these two regimes, U c = μ S ( N ) . The actual validity of these simple heuristics have been tested in several numerical experiments that are reported in Figure 1 and Figure 2.
Figure 1. End probability distribution P ( x , T ) at time T = 200 > τ relax = μ s σ 2 = 13 . 58 1 2 of swarms of N = 500 agents, with respect to U. Here, all agents initially start at x 0 = - 15 . (a): tight swarm for U = 3 U c = 3 μ S (red) and “diffusive evaporating” swarm for U = U c 3 (blue). (b): tight swarm for U = 2 U c (red) and “diffusive evaporating” swarm for U = U c 2 (blue). The respective observation ranges U are depicted at the top of each figure.
Figure 1. End probability distribution P ( x , T ) at time T = 200 > τ relax = μ s σ 2 = 13 . 58 1 2 of swarms of N = 500 agents, with respect to U. Here, all agents initially start at x 0 = - 15 . (a): tight swarm for U = 3 U c = 3 μ S (red) and “diffusive evaporating” swarm for U = U c 3 (blue). (b): tight swarm for U = 2 U c (red) and “diffusive evaporating” swarm for U = U c 2 (blue). The respective observation ranges U are depicted at the top of each figure.
Entropy 18 00027 g001
Figure 2. End probability distribution P ( x , T ) at time T = τ relax = 13 . 58 of a swarm of N = 500 agents, initially starting at x 0 = - 15 . The observation range U = U c = μ S is depicted at the top of the figure. Notice that the evaporation (i.e., destruction of swarm tightness) does not start before the relaxation time τ relax .
Figure 2. End probability distribution P ( x , T ) at time T = τ relax = 13 . 58 of a swarm of N = 500 agents, initially starting at x 0 = - 15 . The observation range U = U c = μ S is depicted at the top of the figure. Notice that the evaporation (i.e., destruction of swarm tightness) does not start before the relaxation time τ relax .
Entropy 18 00027 g002

4. Interactions between Collinear Colliding Swarms

In this section, we investigate the mutual interactions between two initially distinct swarms S 1 and S 2 composed of A 1 , k , k = 1 , 2 , , N 1 respectively A 2 , j , j = 1 , 2 , , N 2 agents driven by rank-based interactions. Initially, we let all agents from S 1 be randomly positioned with X 0 , 1 [ l 1 , r 1 ] : = Ω 1 ; similarly, for S 2 , we assume that X 0 , 2 [ l 2 , r 2 ] : = Ω 2 . Assume non-overlapping supports Ω 1 Ω 2 = with r 1 < l 2 (interval Ω 1 is located on the left of Ω 2 ) and finite observation ranges U such that U < l 2 - r 1 ). Accordingly, at time t = 0 , agents from S 1 do not interact with agents from S 2 and vice versa. We further assume that initially the swarm S k with k = 1 , 2 has an average barycentric drift velocity d k > 0 and that d 1 > d 2 . This implies that the S 1 -agents will, on average, catch and ultimately overcome their S 2 -fellows. Therefore, for large times, the S 1 -agents exhibit a net tendency to become S 2 -leaders. Qualitatively, after a S 1 - S 2 collision, the following dynamic features occur:
(i)
S 1 -agents become leaders and thus are less influenced by (or possibly almost independent from) S 2 -agents. This implies that, with time, S 1 -agents will exhibit a net tendency to recover their nominal drifts (i.e., the drifts realized before the collision).
(ii)
S 2 -agents feel the presence of their leaders from S 1 and, therefore, have a net tendency to increase their drifts.
Therefore, depending on the values of the population sizes N 1 and N 2 and drifts d 1 and d 2 , two alternative situations emerge:
(a)
Mutual capture of swarms. After collision, S 1 and S 2 aggregate into a global tight swarm S glob . The average barycentric drift d glob of S glob is larger than d 2 . Thus, in this case, mutual interactions generate an increase in the average velocity of the global population.
(b)
Quasi-asymptotic freedom of swarms. After collision, and for asymptotic times, the swarms S 1 and S 2 evolve almost without interactions (quasi-free swarm evolutions). In the quasi-free regime, S 1 and S 2 recover their respective initial individual barycentric velocities d 1 and d 2 .
In the sequel, the implementation of the velocity condition d 1 > d 2 is realized by considering two distinctive rank-based dynamics: (i) The Hybrid-Atlas model (HAM) introduced by T. Ichiba et al. [9] and (ii) a new modified HAM model (MHAM). For both cases, we analytically estimate the threshold values of the relevant control parameters that lead to either outcome (a) or (b).

4.1. Colliding Swarms Driven by Hybrid Atlas Models

In this first situation, we assume that the swarms S 1 and S 2 are populated by slightly different types of agents. Agents S 1 and S 2 share a common rank-based drift, given in Equation (1), but we endow the S 1 s with an extra systematic constant drift Γ 1 > 0 :
d U , i ( X ( t ) , t ) = α N U , i ( t ) N + 𝟙 { i N 1 } Γ 1 ,
with N = N 1 + N 2 .
This ensures that S 1 -agents, on average, initially travel faster than S 2 -agents, and, therefore, d 1 > d 2 . Since the initial S 1 -support Ω 1 is located on the left of Ω 2 , a S 1 - S 2 -collision is unavoidable. At this point, we emphasize that when U = , the global heterogeneous society S glob = S 1 S 2 , when subject to the dynamics given by Equation (1), is a special case of HAM dynamics [9]. In [9], the authors derive the set of combinatorial stability conditions that ensure the existence of a tight swarm stationary state. Heuristically speaking, one has to explicitly verify that, in all possible ranking configurations, the agents S 1 and S 2 are systematically driven by attractive drifts directed toward the barycenter of S glob .
For U < , the hypotheses in [9] are not strictly realized (since the perfect ranking determination would indeed require U = ). Along the same lines as in Section 3, we assume that for large enough but finite U, the agents’ ranking in S glob during collision remains approximately unaffected. Accordingly, the HAM tightness conditions also remain approximately valid in S glob , during collision. The average barycentric drift of S glob , during collision and when S 1 is ahead of S 2 , reads:
γ H A M = 1 N i = 1 N 1 ( α i - 1 N + Γ 1 ) + i = N 1 + 1 N α i - 1 N = ( - 1 + N ) α + 2 N 1 Γ 1 2 N .
Therefore, the recentered rank-based drift d U , l ( X ( t ) , t ) of agent A l from S glob , in this configuration, becomes:
d U , i ( X ( t ) , t ) = α N U , i ( t ) N + 𝟙 { i N 1 } Γ 1 - γ H A M , l = 1 , 2 , , N - 1 .
According to [9], the set of combinatorial tightness stability conditions for HAM dynamics, under any permutation p Σ N of the N = N 1 + N 2 agents read as:
k = 1 l α k - 1 N + 𝟙 { p ( k ) N 1 } Γ 1 - γ H A M < ! 0 , for   l = 1 , 2 , , N - 1 , p Σ N ,
where < ! expresses the necessary condition.
To check the validity of Equation (13), it is sufficient to focus on the most critical agent configuration. This configuration is realized when all faster agents take the ranks 1 , . . . , N 1 (i.e., when S 1 is ahead of S 2 ). This corresponds to the choice l = N 1 in Equation (13), under the permutation p = ( 1 , 2 , . . . , N ) , where:
j = 0 N 1 - 1 α j N + Γ 1 = - N 1 N 2 ( α - 2 Γ 1 ) 2 N .
To ensure that this sum is negative, one must finally have:
Γ 1 < α 2 .
Remark that the condition in Equation (14) does not depend on the individual swarm sizes N 1 and N 2 , with our specific choice of agent drifts Equation (10).

4.2. Colliding Swarms Driven by Modified Hybrid Atlas (MHAM) Dynamics

Here, we assume that the N 1 agents from S 1 follow the dynamics of Equation (1), with α = α 1 ˜ > 0 and similarly, S 2 is subject to the same dynamics but with α = α 2 ˜ < α 1 ˜ . From now on and without loss of generality, we focus only on the α 2 ˜ = 1 and α 1 ˜ = ρ . As before, we initially configure the system so that swarm S 1 and S 2 do not interact at time t = 0 (i.e., X 0 , 1 = [ l 1 , r 1 ] = Ω 1 , X 0 , 2 = [ l 2 , r 2 ] = Ω 2 , Ω 1 Ω 2 = and the common observation range U is such that U > l 2 - r 1 ). When ρ > 1 , we have d 1 > d 2 and we look for a critical ratio ρ c under which a tight swarm survives a S 1 - S 2 collision. To this aim, let us consider a “post-colliding” configuration realized when S 1 is completely ahead of S 2 . In term of stability, this configuration is the most critical: S 1 -agents can only escape from S 2 via this configuration. Therefore, the S glob tightness has to be checked under this specific configuration.
Accordingly, whenever the stability conditions from [9] are verified under this critical configuration, they will also be verified under any other configuration. Therefore, the tightness of S glob only depends on this critical configuration. Formally, we now rewrite Equation (5) for S glob and verify whether the stability requirements given by [9] are fulfilled for this critical case. In this case, the recentered dynamics read:
d U , k ( X ( t ) , t ) = ρ N U , k ( t ) N - γ M H A M for k N 1 , N U , k ( t ) N - γ M H A M otherwise . γ M H A M = ρ N k = 1 N 1 k - 1 N + 1 N k = N 1 + 1 N k - 1 N = ρ N 1 ( N 1 - 1 ) + N 2 ( 2 N 1 + N 2 - 1 ) 2 N 2 .
From [9], the stability conditions for S glob read:
k = 1 l ( 1 + 𝟙 { k N 1 } ( ρ - 1 ) ) · k - 1 N - γ M H A M < ! 0 , 1 l N - 1 .
The most critical condition is realized for l = N 1 (this condition tests whether the distance between S 1 and S 2 remains stationary). It reads:
k = 1 N 1 ρ k - 1 N - γ M H A M = N 1 N 2 ( N 2 - 1 - N 1 ( ρ - 2 ) + ρ ) 2 N 2 < ! 0 ,
implying that the critical threshold ρ c finally reads:
ρ c = 1 + N N 1 - 1 .
Notice that for N 1 = 1 , ρ c cannot possibly exist (i.e., ρ c = ). The first agent has no drift, since it never observes any leader. Therefore, a society of one leader cannot possibly drive a swarm. Instead a single leader always remains flocked with S 2 . For N 1 = N 2 , Figure 3 shows the critical value ρ c . When N , ρ c asymptotically converges toward ρ c = 3 . Figure 4 shows the behavior of ρ c as a function of the proportion N 1 N (faster agents in S glob ), for N agents. For N and N 1 N 1 , ρ c 2 .
Figure 3. Value of the critical threshold ρ c , with regard to the swarm size N of the HAM, when N 1 = N 2 . The asymptotic value N of ρ c = 3 is depicted in red.
Figure 3. Value of the critical threshold ρ c , with regard to the swarm size N of the HAM, when N 1 = N 2 . The asymptotic value N of ρ c = 3 is depicted in red.
Entropy 18 00027 g003
Figure 4. Value of the critical threshold ρ c , with regard to the proportion N 1 / N of fast agents in the HAM, when N .
Figure 4. Value of the critical threshold ρ c , with regard to the proportion N 1 / N of fast agents in the HAM, when N .
Entropy 18 00027 g004

4.3. Numerical Simulations

We consider two different cases, both with N 1 = N 2 = 250 agents. In both simulations, the society S 1 starts with all agents located at x 0 , 1 = - 15 , while agents from S 2 start at x 0 , 2 = + 15 . The common observation range U is selected large enough ( U > max ( μ S , 1 , μ S , 2 ) ), and σ = 1 . For the experiments, we vary ρ, in order to compare i) a case in which both societies remain flocked together with ii) a case leading to the quasi-freedom behavior.
Figure 5 shows one realization in which ρ < ρ c , allowing both societies to remain flocked. The left plot shows the final distribution of agents from both societies, as S 1 leads S 2 . In the right plot of Figure 5, we show the average barycentric speeds of both societies. Notice how S 1 is initially accelerated as its agents interact with those from S 2 . Later, when S 1 has overcome S 2 , agents from S 2 get an extra drift due to interactions with S 1 . This finally enables both societies to evolve at the same average barycentric speed (and therefore to remain flocked).
In Figure 6, we show the same scenario, for one realization with ρ > ρ c (quasi-freedom regime). In the right plot, we still notice the initial drift gained by S 1 as its agents interact with those from S 2 . Once S 1 has overcome S 2 , agents from S 2 get an extra drift, that is not large enough to keep both societies flocked together. The extra drift gained by S 2 finally vanishes, once S 1 exits the observation range of the agents from S 2 . Finally, S 1 and S 2 recover their nominal speeds.
Figure 5. Final distributions P 1 , 2 ( x , T ) and average barycentric speeds v 1 , 2 ( t ) for one realization, in which both societies that start at x 0 , { 1 , 2 } = 15 achieve flocking. Here, ρ = 2 . 5 < ρ c , U = 20 > U c = 5 . 4 , σ = 1 and N 1 = N 2 = 250 agents. (a): distribution of the agents from S 1 (blue) and S 2 (red), at ending time T = 100 . (b): average barycentric speed of each society, with respect to time.
Figure 5. Final distributions P 1 , 2 ( x , T ) and average barycentric speeds v 1 , 2 ( t ) for one realization, in which both societies that start at x 0 , { 1 , 2 } = 15 achieve flocking. Here, ρ = 2 . 5 < ρ c , U = 20 > U c = 5 . 4 , σ = 1 and N 1 = N 2 = 250 agents. (a): distribution of the agents from S 1 (blue) and S 2 (red), at ending time T = 100 . (b): average barycentric speed of each society, with respect to time.
Entropy 18 00027 g005
Figure 6. Final distributions P 1 , 2 ( x , T ) and average barycentric speeds v 1 , 2 ( t ) for one realization, in which both societies that start at x 0 , { 1 , 2 } = 15 do not achieve flocking. Here, ρ = 6 > ρ c , U = 20 > U c = 2 . 26 , σ = 1 and N 1 = N 2 = 250 agents. (a): distribution of the agents from S 1 (blue) and S 2 (red), at ending time T = 100 . (b): average barycentric speed of each society, with respect to time.
Figure 6. Final distributions P 1 , 2 ( x , T ) and average barycentric speeds v 1 , 2 ( t ) for one realization, in which both societies that start at x 0 , { 1 , 2 } = 15 do not achieve flocking. Here, ρ = 6 > ρ c , U = 20 > U c = 2 . 26 , σ = 1 and N 1 = N 2 = 250 agents. (a): distribution of the agents from S 1 (blue) and S 2 (red), at ending time T = 100 . (b): average barycentric speed of each society, with respect to time.
Entropy 18 00027 g006

5. Conclusions and Perspectives

Tightly flocked swarms of agents are externally perceived as plastic (i.e., deformable) macroscopic bodies with dynamics resulting from the compromise between the individual evolutions, and the nature and range of the mutual interactions. Having macroscopic bodies at hand, basic physics naturally suggests the questions: (i) how such do solid bodies mutually interact; (ii) what emerges after collisions, and (iii) what type of information can we get from studying the results of the collisions? The discussion of collisions presents some similarities with the physics of solitons’ interactions where the underlying nonlinearities precludes the superposition principle to be be invoked. We show here that the rank-based Brownian dynamics is one possible approach to explore, in a partly analytical way, some of the challenging questions related to swarm interactions.

Acknowledgments

Detailed comments from anonymous referees strongly contributed to the presentation and quality of this paper. Guillaume Sartoretti is funded by the Swiss National Science Foundation.

Author Contributions

Max-Olivier Hongler came up with the core concept of the paper; Guillaume Sartoretti and Max-Olivier Hongler made the analytical calculations; Guillaume Sartoretti performed the numerical experiments; Max-Olivier Hongler and Guillaume Sartoretti wrote the paper. Both authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Couzin, I.D.; Krause, J.; Franks, N.R.; Levin, S.A. Effective leadership and decision-making in animal groups on the move. Nature 2005, 433, 513–516. [Google Scholar] [CrossRef] [PubMed]
  2. Sartoretti, G.; Hongler, M.O.; Filliger, R. The Estimation Problem and Heterogenous Swarms of Autonomous Agents. In Proceedings of Stochastic Modeling Techniques and Data Analysis International Conference (SMTDA 2014), Lisbon, Portugal, 11–14 June 2014.
  3. Han, J.; Li, M.; Guo, L. Soft control on collective behavior of a group of autonomous agents by a shill agent. J. Syst. Sci. Complex. 2006, 19, 54–62. [Google Scholar] [CrossRef]
  4. Sartoretti, G.A.; Hongler, M.O. Soft Control of Swarms: Analytical Approach. In Proceedings of the 5th International Conference on Agents and Artificial Intelligence (ICAART 2013), Barcelona, Spain, 15–18 February 2013; pp. 147–153.
  5. Sartoretti, G.; Hongler, M.-O. Soft Control of Self-Organized Locally Interacting Brownian Planar Agents. In Computer Aided Systems Theory - EUROCAST 2013; Moreno-Díaz, R., Pichler, F., Quesada-Arencibia, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; pp. 45–52. [Google Scholar]
  6. Strogatz, S.H. From Kuramoto to Crawford: Exploring the onset of synchronization in populations of coupled oscillators. Physica D 2000, 143, 1–20. [Google Scholar] [CrossRef]
  7. Acebrón, J.A.; Bonilla, L.L.; Vicente, C.J.P.; Ritort, F.; Spigler, R. The Kuramoto model: A simple paradigm for synchronization phenomena. Rev. Mod. Phys. 2005, 77. [Google Scholar] [CrossRef]
  8. Pal, S.; Pitman, J. One dimensional Brownian particle systems with rank-dependent drifts. Ann. Appl. Probab. 2008, 18, 2179–2207. [Google Scholar] [CrossRef]
  9. Ichiba, T.; Papathanakos, V.; Banner, A.; Karatzas, I.; Fernholz, R. Hybrid Atlas Model. Ann. Appl. Probab. 2011, 21, 609–644. [Google Scholar] [CrossRef]
  10. Hongler, M.O.; Gallay, O.; Hashemi, F. Imitation’s Impact on the Dynamics of Long-Wave Growth. Econ. Model. 2015. submited. [Google Scholar]
  11. Jourdain, B.; Malrieu, F. Propagation of chaos and Poincaré inequalities for a system of particles interacting through their CDF. Ann. Appl. Probab. 2008, 18, 1706–1736. [Google Scholar] [CrossRef]
  12. Hongler, M.O.; Filliger, R.; Gallay, O. Local versus nonlocal barycentric interactions in 1D agent dynamics. Math. Biosci. Eng. 2014, 11, 303–315. [Google Scholar] [PubMed]
  13. Banner, A.D.; Fernholz, R.; Karatzas, I. Atlas Models of Equity Markets. Ann. Appl. Prob. 2005, 15, 2296–2330. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Sartoretti, G.; Hongler, M.-O. Interacting Brownian Swarms: Some Analytical Results. Entropy 2016, 18, 27. https://0-doi-org.brum.beds.ac.uk/10.3390/e18010027

AMA Style

Sartoretti G, Hongler M-O. Interacting Brownian Swarms: Some Analytical Results. Entropy. 2016; 18(1):27. https://0-doi-org.brum.beds.ac.uk/10.3390/e18010027

Chicago/Turabian Style

Sartoretti, Guillaume, and Max-Olivier Hongler. 2016. "Interacting Brownian Swarms: Some Analytical Results" Entropy 18, no. 1: 27. https://0-doi-org.brum.beds.ac.uk/10.3390/e18010027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop