
==== Front
R Soc Open Sci
R Soc Open Sci
RSOS
royopensci
Royal Society Open Science
2054-5703
The Royal Society

rsos240794
10.1098/rsos.240794
100110011004304424Physics and Biophysics
Research Articles
Non-local interaction in discrete Ricci curvature-induced biological aggregation
Non-local interaction in discrete Ricci curvature-induced biological aggregation
https://orcid.org/0000-0002-7171-7960
Beuria Jyotiranjan 1 2 Conceptualization Formal analysis Methodology Writing – original draft Writing – review and editing jbeuria@iitmandi.ac.in

Behera Laxmidhar 1 3 Conceptualization Supervision Writing – review and editing lbehera@iitk.ac.in

1 IKSMHA Center, IIT Mandi , Mandi, India
2 IKS Research Center, ISS Delhi , Delhi, India
3 Department of Electrical Engineering, IIT Kanpur , Kanpur, India
9 2024
04 9 2024 September 4, 2024
04 9 2024 September 4, 2024
11 9 24079414 5 2024 May 14, 2024
20 7 2024 July 20, 2024
22 7 2024 July 22, 2024
© 2024 The Author(s).
2024
https://creativecommons.org/licenses/by/4.0/ Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.

We investigate the collective dynamics of multi-agent systems in two- and three-dimensional environments generated by minimizing discrete Ricci curvature with local and non-local interaction neighbourhoods. We find that even a single effective topological neighbour suffices for significant order in a system with non-local topological interactions. We also explore topological information flow patterns and clustering dynamics using Hodge spectral entropy and mean Forman–Ricci curvature.

biocomplexity
; biophysics
; applied mathematics
==== Body
pmc1. Introduction

Biological aggregation is one of the most ubiquitous and spectacular displays of coordinated behaviour of a collection of self-propelled particles [1–3]. The universality of this phenomenon can be seen in systems of very different sizes and scales, such as the flocking of birds [4,5], schooling of fish [6–8], bacterial colonies [9,10], locust swarms [11], sheep herds [12,13], and even human crowding [14–16] or robot swarming [17]. One of the characteristics of such a coordinated motion is the emergence of order in the system through the exchange of information among the agents. Several models for collective motion have been developed, primarily through local microscopic interaction mechanisms. Also, some recent studies have shown that active agents can undergo flocking without even explicitly aligning with neighbours [18–20]. However, a common feature of biological aggregation is to spread the consensus direction of motion among the constituent agents.

The study of collective motion is an active area of research both theoretically [18] and experimentally [21]. One of the most popular models of biological aggregation is the Vicsek model [22]. In this model, during each time step, the velocity of each agent is adjusted to match the average velocity of its neighbouring agents. Additionally, a random unit vector, scaled by a fixed noise strength, is added to this adjusted velocity. All the agents update their positions simultaneously at every time step based on the orientations. The interaction radius r defines the neighbourhood of interaction. Any two neighbouring agents must be separated by less than or equal to r. In other words, agents can only be influenced by their neighbours up to a radial distance r. In this model, an interesting phenomenon happens that all agents’ headings will converge to the same value or consensus for high group densities and low noise levels. The collective motion reaches a consensus faster with a larger interaction radius r, implying a denser interaction network and faster information flow.

Most of the models for collective motion are primarily based on local or short-range interaction neighbourhoods. Interestingly, some recent studies have examined interaction mechanisms extending beyond strict locality [23–26]. Proper consideration of non-local interactions is essential when the sensory systems of animals or birds are known to support such interactions physiologically. Natural swarms often exhibit long-ranged correlations [27], and many flocking animals or birds rely heavily on their vision [28,29], which can extend well beyond the immediate flock size. A study by Strandburg-Peshkin et al. [30] has demonstrated that biological organisms efficiently determine their movement directions using line-of-sight interactions. Notably, models with non-local interactions can easily produce velocity correlations across large distances [26]. This condition is more stringent for models when interactions are only local.

On top of that, the information flow among the agents is best modelled using a topological neighbourhood. The research conducted by Ballerini et al. [4] on starling flocks revealed that interaction among individuals is topological. Each agent coordinates with a fixed number, approximately seven, of its closest neighbours, regardless of their distances. This finding contradicts the assumptions made by most models of self-organized collective motion, which typically rely on metric interactions. Also, purely metric interactions require a much larger number of interacting neighbours [31,32]. Ballerini et al. have argued that topological interactions provide stronger cohesion than metric interactions, making them more effective from an anti-predatory perspective. Also, comparing metric to topological models in three dimensions, topological models exhibit superior stability [31], especially with spatially balanced neighbour selection, supporting findings from starling flock experiments [4]. This spatially balanced topological model requires fewer interacting neighbours. One of the goals of this work is to demonstrate the efficacy of discrete Ricci curvature-driven topological interaction in reducing the number of neighbours to as low as one.

From the above discussions, we realize that non-locality and the number of interacting neighbours are crucial for biological aggregation. Thus, we present a topological model that leverages the metric-based features and non-locality by constructing the Vietoris–Rips complex. In this model, an agent’s new orientation is not obtained by averaging over all possible nearest metric or topological neighbours. Instead, we prescribe a mechanism for following the direction of the most curved nearest neighbour(s) within an annular region around the agent. This model represents the entire system at a particular time step as a topological space using a Vietoris–Rips complex. The Ricci curvature measures how much the volume of a small ball in space changes compared to Euclidean space. We calculate a variant of discrete Ricci curvature by Forman [33] for the edges and vertices in this simplicial complex.

From a geometric perspective, the Forman–Ricci curvature (FRC) measures the information flow at the ends of edges within a network [34]. A large negative value for FRC of an edge indicates a higher spread of information at the ends of an edge. Such an edge has numerous neighbouring edges at both vertices, resembling a funnel connecting many other vertices. Interestingly, it has been observed that most real-world networks possess negative FRCs [35] for the edges. Thus, we expect that minimizing FRC while selecting interacting neighbours will favour biological aggregation. We also study the Hodge spectral entropy [36,37] that might offer insights into the higher-order information flow in the topological space formed by the agents. However, before we describe the model and the numerical simulations, we delve into a brief mathematical introduction to the rich field of topology and discrete Ricci curvature.

2. Forman–Ricci curvature and topology

2.1. Simplicial complex of the point cloud

Constructing simplicial complexes is central to computational geometry techniques such as topological data analysis. It offers a systematic approach to constructing topological spaces using basic combinatorial units called simplices, facilitating the study of complex geometries through more manageable combinatorial and counting tasks. A k-dimensional simplex arises by taking a convex hull of k+1 points (see figure 1). For example, a 0-simplex corresponds to a vertex, a 1-simplex to an edge, a 2-simplex to a filled triangle and a 3-simplex to a filled tetrahedron. This construction can also be generalized to higher-dimensional polytopes. In a k-dimensional simplex, simplices with dimensions <k constitute its faces. Thus, for a 2-simplex, its edges and vertices constitute the faces.

Figure 1. 0-simplex is a vertex, 1-simplex is an edge, 2-simplex is a filled triangle and 3-simplex is a filled tetrahedron.

0-simplex is a vertex, 1-simplex is an edge, 2-simplex is a filled triangle and 3-simplex is a filled tetrahedron.

Let X={x0,x1,...,xm−1}∈ℝm represent a point cloud data, where each point xi resides in ℝn. Let r denote a fixed radius. The Vietoris–Rips complex of X is an abstract simplicial complex. Its k−simplices consist of (k + 1)-tuples of points in X such that the pairwise distance between them is less than or equal to r. This maximal radial limit r is also referred to as the filtration parameter. Mathematically, the Vietoris–Rips complex, also known as the Rips complex K, is expressed as:

(2.1) K={σ⊆X|d(xi,xj)≤r∀xi≠xj∈σ},

where d(xi,xj) denotes the Euclidean metric between two points.

2.2. Constructing chain complexes

A k−chain, denoted by Ck, represents a formal sum of a subset of k−simplices belonging to the simplicial complex K. It is expressed as:

(2.2) Ck=∑iαiσki,

where σk represents the k−simplex and αi is a real number (αi∈ℝ). It is to be noted that Ck forms an abelian group under component-wise addition. The k−chain group is generated by the k−cycles, which are k−chains without boundaries. The k−th boundary operator ∂k on a k−simplex σk with vertices (v0,v1,...,vk) is defined as:

(2.3) ∂k(σk)=∑i=0k(−1)i(v0,v1,...,v^i,...,vk),

where v^i denotes the vertex removed from σk. This operation maps a k−chain Ck to a (k−1)-chain Ck−1. A chain complex is formed by the sequence of boundary operators acting on the chain groups:

(2.4) Ck→∂kCk−1→∂k−1Ck−2...C1→∂1C0→∂00.

The kernel of the boundary operator ∂k consists of all k-chains Ck that have no boundary, while the image of ∂k is the set of (k−1)-chains Ck−1 that represent the boundaries of k-chains Ck. This can be mathematically expressed as

(2.5) ker(∂k)={c∈Ck|∂kCk=0}im(∂k)={d∈Ck−1|∃c∈Ck:d=∂k(c)}.

We note that k-boundaries are elements of im(∂k+1), whereas elements of ker(∂k) correspond to k-cycles. Moreover, under addition, the sets of k-boundaries Bk and k-cycles Zk form abelian subgroups of Ck. Since k-boundaries are also k-cycles, it is important to observe that ker(∂k)⊂im(∂k+1). With Bk⊂Zk⊂Ck, the groups Bk, Zk and Ck exhibit a nested structure. The quotient group that represents cycles modulo boundaries is called the k-th homology group, or Hk. It can be stated as:

(2.6) Hk=ZkBk=ker(∂k)im(∂k+1).

Here, Hk(K) represents the quotient vector space whose generators are given in terms of k-cycles that are not boundaries of any (k+1)-simplices. The rank of Hk(K) is often referred to as the k-th Betti number βk(K). The Betti number βk(K) indicates the number of k-dimensional holes in the simplicial complex K that are not boundaries of any (k+1)-simplices. For example, β0(K) is the number of connected components in the simplicial complex. Euler characteristic χ is a topological invariant of the topological space. It is expressed as an alternate sum of Betti numbers βk. Mathematically,

(2.7) χ=∑k=0n(−1)kβk.

2.3. Forman Ricci curvature

Let α and α‾ denote k-dimensional simplices within the simplicial complex K. If there exists a simplex β in K such that both β>α and β>α‾, then α and α‾ possess a common co-face β, and they are termed as upper adjacent. Similarly, α and α‾ are considered lower adjacent if they share a common face γ, a (k−1)-simplex with γ<α and γ<α‾. When α and α‾ are either lower or upper adjacent but not both, they are referred to as parallel. The FRC is expressed as [33]

(2.8) Rk(α)=N(Upper adjacent simplices)+N(Lower adjacent simplices)−N(Parallel simplices).

In the context of weighted simplicial complexes with weights denoted as w, Rk(α) is expressed as

(2.9) Rk(α)=wα[∑β>αwβwα+∑γ<αwγwα]−wα∑α¯≠α[∑β>α,α¯wαwα¯wβ−∑γ<α,α¯wγwαwα¯].

Considering up to 1-simplex in simplicial complex K, FRC for an edge reduces to

(2.10) R1(α)=wα(∑γ<αwγwα−∑α¯≠α∑γ<α,α¯wγwαwα¯).

Recently, there has been considerable interest in using FRC to study real networks [35,38–41] we encounter in various fields.

2.4. Hodge spectral entropy

The Hodge Laplacian is computed with respect to boundary operators. The Hodge Laplacian L[n] on a d-dimensional simplicial complex is defined for each n=0,…,d as

(2.11) L[0]=L[0]upforn=0,

(2.12) L[n]=L[n]up+L[n]downforn>0,

where L[n]up and L[n]down encode the diffusion from n-simplices to n-simplices via (n+1)-simplices and (n−1)-simplices, respectively. They can be expressed in terms of the boundary operators as follows:

(2.13) L[n]down=BnBn†,

(2.14) L[n]up=Bn+1†Bn+1.

Using the Hodge Laplacian, we compute the von Neumann entropy for the simplicial complex as follows [36]:

(2.15) Snhs=−Tr(ρnlnρn),

where ρn is expressed in terms of the Hodge Laplacian as follows:

(2.16) ρn=e−βL[n]Tr(e−βL[n]),

where β is the damping factor and is different from the Betti numbers βk discussed earlier.

From a computational point of view, Snhs can also be expressed in terms of the partition function Zn and eigenvalues of L[n] as follows:

(2.17) Snhs=β⟨λn⟩+ln(Zn),

(2.18) Zn=Tr(e−βL[n]),

(2.19) ⟨λn⟩=∑ie−βλi(L[n])λi(L[n])Zn,

where λi(L[n]) is the ith generic eigenvalue of L[n]. With these mathematical preliminaries introduced above, we are ready to explore the model based on discrete Ricci curvature.

3. Ricci curvature-based non-local model

The model consists of N self-propelled agents moving at a constant speed (v0) in a periodic box of length L and dimension d. The orientation of the agents keeps changing randomly. The model we consider in this work is a variant of the standard Vicsek model [2,22]. The Vicsek model is simple yet very powerful in predicting collective behaviour. The velocity of an agent is expressed as follows:

(3.1) v→i(t+1)=1|Ni|∑j∈Niv→j(t)+v0ηr^i(t),

where

— v→i(t+1) is the velocity of agent i at time t+1.

— |Ni| represents the number of nearest neighbours of agent i.

— v→j(t) is the velocity of neighbour j of agent i at time t.

— η is the noise strength.

— r^i(t) is a random unit vector representing noise.

This equation represents the update rule for each agent’s velocity in the Vicsek model, where agents adjust their velocities to align with the average velocity of their neighbours plus a random noise term. The variants of the Vicsek model differ in how an agent’s nearest neighbours are selected. There are two primary variations, one with a topological neighbourhood and the other with a metric neighbourhood.

In figure 2, we suggest a non-local scheme for finding the most effective nearest neighbours. In the standard Vicsek model, the neighbours are found within a cutoff radius around the agent. Following King and Turner [26], we construct an annular disk of inner radius (rmin) and outer radius (rmax) around every agent. Traditional network models find the k−nearest neighbours of an agent and average their orientations to update the agent’s orientation up to a noise factor governed by η. However, we first chose to embed the entire system in a Vietoris–Rips complex. Thus, rather than dealing with fragmented portions of the system, we wish to capture both the local and non-local factors in deciding the new orientation of the agent.

Figure 2. Schematic representation of the simplicial complex (up to edges) structure and the neighbouring agents in an annular region around an agent.

Schematic representation of the simplicial complex (up to edges) structure and the neighbouring agents in an annular region around an agent.

The non-locality enters in two folds. First, we exclude the closest agents by choosing rmin>0. Second, instead of averaging over all the effective agents in the annular disk, we perform a weighted averaging of the orientations of the agents featuring the lowest FRC. We assign weights to the nodes and edges in the simplicial complex formed by the ensemble of agents. This process can also be extended to higher-order simplices for more precise results. Due to computational limitations, we restrict the homology dimension to 1. The weighting scheme is as follows.

Each edge’s weight, denoted as wij, is determined by the sum of the weights of the two vertices it connects, divided by 1+rij, where rij represents the distance between the vertices. For each edge connecting vertices i and j, the weight wij is given by:

(3.2) wij=wi+wj1+rij,

where wi and wj are the weights of vertices i and j given by wi=1+degree(i). This weight assignment scheme ensures that each edge’s weight is influenced by the combined weights of its connected vertices and inversely proportional to their distance. Also, the vertices with a larger number of connecting edges get larger weights.

The vertex FRC Ri for agents is calculated as the sum of Rij over all edges incident to the ith vertex. Mathematically, it can be expressed as [42]:

(3.3) Ri=∑jRij,

where Rij represents the curvature associated with the edge connecting the ith vertex to its neighbouring vertices. This vertex curvature measure aggregates of the curvature contributions from all edges incident to a particular vertex, providing insights into the curvature at individual vertices within the network [34].

To determine the updated direction of movement of an agent, we select the k vertices with the lowest Ri values, where Ri represents the vertex curvature. Let vi denote the velocity vector associated with the i-th vertex. The weighted average velocity v→i(t+1) is then calculated as follows:

(3.4) v→i(t+1)=v0∑i=1kwi⋅v→i(t)‖∑i=1kwi⋅v→i(t)‖+v0ηr^i(t),

where wi represents the weight assigned to the ith vertex, and k is the number of effective interacting vertices. This procedure ensures that the collective direction of movement is influenced by the directions of the most curved agents from the chosen k agents. In case of only one effective nearest neighbour, the expression becomes

(3.5) v→i(t+1)=v→Rmini(t)+v0ηr^i(t),

where v→Rmini is the velocity of the most curved agent for i-th agent. As we will see in the numerical results, following the most curved agent can also bring appreciable order within the system. Thus, it is an economical model in terms of the number of nearest neighbours to be considered.

Before we proceed to analyse the results, it is important to define the order parameter. The order parameter ⟨ϕv⟩ represents the average alignment of velocities of agents in a system. It provides insights into the collective order or alignment in the system. Traditionally, it has been studied for phase transition in Vicsek-like models. Mathematically, it is expressed as:

(3.6) ⟨ϕv⟩=1N‖∑i=1Nv→i‖v→i‖‖

We also define a very useful quantity, radial distribution function g(r), given by:

(3.7) g(r)=1ρN∑i≠j⟨δ(r−|ri−rj|)⟩,

where the average is performed over the shell volume. g(r) gives a quantitative estimate of the density of particles at distance r from an agent. From a practical point of view, one divides the volume into shells of width Δr at a certain radial distance r and counts the number of agents in that shell [43]. Peaks in g(r) are the signature of different kinds of structure formation in the flock due to stronger attractive forces. Similarly, the troughs in g(r) indicate repulsive forces in the flock.

4. Numerical results

We simulate N=100 agents for two-dimensional and N=1000 for the three-dimensional periodic box for 105 time steps such that the number density ρ=1 for both two- and three-dimensional. The distance between two agents is the usual Euclidean distance in a periodic box. The speed of agents v0 is fixed at 0.5. We choose |rmax−rmin|=1 with two representative r values, rmin=0 and rmin=1. rmin=0 represents the so-called local model, and rmin=1 is the non-local model. This is because, for rmin>0, the nearby local agents do not play a role explicitly. We average different features over 10 samples of simulation for each combination of rmin and η to obtain statistical stability. As briefly discussed in the introduction section, the experiment on the starling flock favours the number of effective neighbours as low as seven [4]. Thus, the goal of our discussion will be to see how this model performs with as little as one nearest neighbour. We will also contrast the results with the seven nearest neighbours model. The computation was performed using CuPy: a NumPy-Compatible Library for NVIDIA GPU calculations [44] on a workstation with an Intel Xeon W-3265 CPU and A4000 NVIDIA GPU.

In figures 3 and 4, we present the evolution of the order parameter ⟨ϕv⟩ with time in two- and three-dimensional cases, respectively, for zero and non-zero noise strength (η=0.2). It is interesting to observe that in the absence of noise (η=0.0), ⟨ϕv⟩ approaches complete order (⟨ϕv⟩=1.0) very quickly even with only one effective nearest neighbour. This is true whether rmin=0 or rmin=1 in both two- and three-dimensional cases. Thus, Ricci curvature-induced local interaction is very efficient and economical in generating coherence at a large scale. Also, as expected, with the introduction of noise (η=0.2), the order disrupts, and NN=1 becomes less effective compared to the NN=7 case.

Figure 3. Time evolution of order parameter for the two-dimensional case. NN denotes the number of effective neighbours being considered. It is to be noted that when ⟨ϕv⟩≈1.0, the lines for NN=1 get buried below that for NN=7. (a) for rmin=0.0 and noise strength η=0.0, (b) for rmin=0.0 and noise strength η=0.2, (c) for rmin=1.0 and noise strength η=0.0 and (d) for rmin=1.0 and noise strength η=0.2.

Time evolution of order parameter for the two-dimensional case.

Figure 4. Time evolution of order parameter for the three-dimensional case. NN denotes the number of effective neighbours being considered. It is to be noted that when ⟨ϕv⟩≈1.0, the lines for NN=1 get buried below that for NN=7. (a) for rmin=0.0 and noise strength η=0.0, (b) for rmin=0.0 and noise strength η=0.2, (c) for rmin=1.0 and noise strength η=0.0 and (d) for rmin=1.0 and noise strength η=0.2.

Time evolution of order parameter for the 2D case.

It is to be noted that there is a difference between two- and three-dimensional cases. Unlike the two-dimensional case in figure 3b, the three-dimensional case with η=0.2 shown in figure 4b says that for rmin=0, only one effective agent is insufficient for attaining appreciable order. However, this is not true for the non-local case with rmin=1 as shown in figure 4d. With the introduction of not-so-large noise strength (η=0.2), the disruption of order happens for both small and large NN. Yet, the non-local model (rmin=1) performs much better compared to the local model (rmin=0).

Next, we study the mean Ricci curvature for NN=1 and NN=7 in two- and three-dimensional periodic boxes. Here, we have performed the average over all time steps and all agents. In figure 5a,b, we present the mean FRC of the vertices with the variation of noise strength η in two dimensions for NN=1 and NN=7, respectively. The three-dimensional counterpart is shown in figure 5c,d. We observe that at η=0.0, NN=1 has large negative R for rmin=0 compared to rmin=1 (see figure 5a,c). A cross-over exists between the cases for rmin=0 and rmin=1 at non-zero η. For the three-dimensional case, the crossover occurs at a smaller η than the two-dimensional case. Interestingly, the NN=7 cases have a slightly different pattern (see figure 5b and d). At η=0.0, R for rmin=1 is more negative compared to rmin=0. For larger η, rmin=1 scenarios have lesser R both in two and three dimensions.

Figure 5. The variation of the mean of the Ricci curvature R as the noise strength (η) varies. (a) NN=1 in two dimensions, (b) NN=7 in two dimensions, (c) NN=1 in three dimensions and (d) NN=7 in three dimensions.

The variation of the mean of the Ricci curvature

The negative R is the signature of the large clustering of agents and edge connectivity. This also tells us about a large information exchange between edges and vertices. A large negative value of R for rmin=0 at a lower η for NN=1 and subsequent crossover with the rmin=1 scenario indicates the order in the system is very much sensitive to the strength of the noise. Although η≈0 favours the local model rmin=0 compared to the non-local rmin=1 one, as η increases, the local model quickly becomes ineffective. In the three-dimensional case with NN=1, the crossover occurs at a lower η than the two-dimensional case. The point of crossover is pushed to a larger η for NN=7. This is expected because a larger number of neighbours will decide the correct alignment even in the presence of stronger noise.

In figure 6, we present the often discussed second-order phase transition of the Vicsek models for the two- and three-dimensional periodic boxes. Irrespective of the number of effective neighbours, the phase transition for the non-local model (rmin=1) takes place at a slightly larger noise strength η compared to the local model (rmin=0). Also, for NN=7, the system remains aligned in the presence of larger η compared to the NN=1 cases. These observations corroborate the previous discussion on mean Ricci curvature that the non-local models are preferred in the presence of larger noise strength.

Figure 6. The variation of the mean order parameter ⟨ϕv⟩ as the noise strength (η) varies. (a) NN=1 in two dimensions, (b) NN=7 in two dimensions, (c) NN=1 in three dimensions and (d) NN=7 in three dimensions.

The variation of the mean order parameter

In figure 7, we present the Hodge spectral entropy (Skhs) for the edges and the vertices in two dimensions as a function of time. The k-th Hodge spectral entropy measures information transmission across k-simplices. In figure 7a and c, we observe that S0hs approaches a constant value for NN=7. It is to be noted that the zeroth Betti number β0 is the number of connected components in the simplicial complex. As the number of effective neighbours increases, the system quickly attains appreciable order, even in the presence of noise. Large S0hs indicates large agent-to-agent information flow through the edges connecting each other. We further note that the time evolution of S0hs is more stable for the non-local case than the local one. In the case of a three-dimensional periodic box (see figure 8), the time evolution of S0hs remains similar to that of the two-dimensional periodic box.

Figure 7. Time evolution of Hodge Spectral entropy Skhs for the local and non-local model in two-dimensional periodic box. (a) S0hs with rmin=0, (b) S1hs with rmin=0, (c) S0hs with rmin=1.0 and (d) S1hs with rmin=1.0.

Time evolution of Hodge Spectral entropy

Figure 8. Time evolution of Hodge Spectral entropy Skhs for the local and non-local model in three-dimensional periodic box. (a) S0hs with rmin=0, (b) S1hs with rmin=0, (c) S0hs with rmin=1.0 and (d) S1hs with rmin=1.0.

Time evolution of Hodge Spectral entropy

S1hs encodes the edge-to-edge information flow. In figure 7b,d, we observe that S1hs for NN=1 is higher compared to the NN=7 case. However, the difference is higher for the non-local case, indicating a larger difference in the information flow through the vertices in the non-local case. As mentioned earlier, the discussion for Hodge spectral entropies S1hs and S0hs is very much similar for the three-dimensional periodic box (see figure 8). In the three dimensional case, the difference between S1hs for NN=1 and NN=7 is stronger compared to the two-dimensional case. We also note that Skhs remains almost constant with time, indicating continuous information flow for all orders of homology.

In figure 9, we present the radial distribution function g(r) for both two- and three-dimensional. In figure 9a, we see that for r≈0.0 and NN=1, the local model features larger g(r) compared to the non-local one. The trend is similar for NN=7. Interestingly, the non-local model with NN=1 features the largest g(r) in the three-dimensional case compared to the local one. However, the trend is the opposite for the NN=1 case. Generically, g(r) is almost featureless, as we do not observe any prominent peaks or troughs.

Figure 9. Radial distribution function g(r) with η=0.2. (a) in two dimensions and (b) in three dimensions.

Radial distribution function

5. Conclusion

The study delves into the Ricci curvature-induced dynamics of self-propelled agents within two- and three-dimensional periodic boxes, contrasting local and non-local interaction models. We find that the system can achieve significant order even with a single effective nearest neighbour. Although such a framework is sensitive to noise, the non-local interaction, particularly in the three-dimensional periodic box, is an effective mechanism to create alignment among the agents. Mean Ricci curvature analysis unveils the system’s sensitivity to noise levels, with local models demonstrating better order at lower noise strengths. However, the non-local model features lower mean Ricci curvature for larger noise than the local model. Furthermore, our analysis of Hodge spectral entropy and mean Ricci curvature provides deeper insights into the topological information transmission within the system. Although the radial distribution function is almost featureless, the analysis highlights a nuanced difference between local and non-local interaction models, particularly in the close vicinity of an agent.

Overall, the study underscores the intricate interplay between the nature of the topological interactions, noise strength and dimensionality. The non-local Ricci curvature-induced topological interaction-based model is effective and economical for collective motion in biological systems. These findings contribute to our understanding of complex systems and can help design effective control strategies for collective motion in various real-world applications.

Acknowledgements

We thank our parent organisations for their kind support.

Ethics

This work did not require ethical approval from a human subject or animal welfare committee.

Data accessibility

Data and relevant code for this research work are stored in GitHub: [45] and have been archived within the Zenodo repository: [46].

Declaration of AI use

We have not used AI-assisted technologies in creating this article.

Authors’ contributions

J.B.: conceptualization, formal analysis, methodology, writing—original draft, writing—review and editing; L.B.: conceptualization, supervision, writing—review and editing.

Both authors gave final approval for publication and agreed to be held accountable for the work performed therein.

Conflict of interest declaration

We declare we have no competing interests.

Funding

No funding has been received for this article. APC is paid by IIT Mandi.
==== Refs
References

1. Sumpter DJ. 2010 Collective animal behavior. Princeton, NJ: Princeton University Press.
2. Vicsek T, Zafeiris A. 2012 Collective motion. Phys. Rep. 517 , 71–140. (10.1016/j.physrep.2012.03.004)
3. Popkin G. 2016 The physics of life. Nature 529 , 16–18. (10.1038/529016a)26738578
4. Ballerini M et al . 2008 Interaction ruling animal collective behavior depends on topological rather than metric distance: evidence from a field study. Proc. Natl Acad. Sci. USA 105 , 1232–1237. (10.1073/pnas.0711437105)18227508
5. Nagy M, Ákos Z, Biro D, Vicsek T. 2010 Hierarchical group dynamics in pigeon flocks. Nature 464 , 890–893. (10.1038/nature08891)20376149
6. Lopez U, Gautrais J, Couzin ID, Theraulaz G. 2012 From behavioural analyses to models of collective motion in fish schools. Interface Focus 2 , 693–707. (10.1098/rsfs.2012.0033)24312723
7. Rosenthal SB, Twomey CR, Hartnett AT, Wu HS, Couzin ID. 2015 Revealing the hidden networks of interaction in mobile animal groups allows prediction of complex behavioral contagion. Proc. Natl Acad. Sci. USA 112 , 4690–4695. (10.1073/pnas.1420068112)25825752
8. Miller N, Gerlai R. 2012 From schooling to shoaling: patterns of collective motion in zebrafish (Danio rerio). PLoS One 7 , e48865. (10.1371/journal.pone.0048865)23166599
9. Zhang HP, Be’er A, Florin EL, Swinney HL. 2010 Collective motion and density fluctuations in bacterial colonies. Proc. Natl Acad. Sci. USA 107 , 13626–13630. (10.1073/pnas.1001651107)20643957
10. Rabani A, Ariel G, Be’er A. 2013 Collective motion of spherical bacteria. PLoS One 8 , e83760. (10.1371/journal.pone.0083760)24376741
11. Buhl J, Sword GA, Simpson SJ. 2012 Using field data to test locust migratory band collective movement models. Interface Focus 2 , 757–763. (10.1098/rsfs.2012.0024)24312729
12. Gómez-Nava L, Bon R, Peruani F. 2022 Intermittent collective motion in sheep results from alternating the role of leader and follower. Nat. Phys. 18 , 1494–1501. (10.1038/s41567-022-01769-8)
13. Ginelli F, Peruani F, Pillot MH, Chaté H, Theraulaz G, Bon R. 2015 Intermittent collective dynamics emerge from conflicting imperatives in sheep herds. Proc. Natl Acad. Sci. USA 112 , 12729–12734. (10.1073/pnas.1503749112)26417082
14. Zhou B, Tang X, Wang X. 2013 Measuring crowd collectiveness. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 3049–3056. (10.1109/CVPR.2013.392)
15. Ma Y, Lee EWM, Shi M, Yuen RKK. 2021 Spontaneous synchronization of motion in pedestrian crowds of different densities. Nat. Hum. Behav. 5 , 447–457. (10.1038/s41562-020-00997-3)33398140
16. Keta YE, Jack RL, Berthier L. 2022 Disordered collective motion in dense assemblies of persistent particles. Phys. Rev. Lett. 129 , 048002. (10.1103/PhysRevLett.129.048002)35939008
17. Zhao H, Liu H, Leung YW, Chu X. 2018 Self-adaptive collective motion of swarm robots. IEEE Trans. Autom. Sci. Eng. 15 , 1533–1545. (10.1109/TASE.2018.2840828)
18. Shaebani MR, Wysocki A, Winkler RG, Gompper G, Rieger H. 2020 Computational models for activematter. Nat. Rev. Phys. 2 , 181–199. (10.1038/s42254-020-0152-1)
19. Bär M, Großmann R, Heidenreich S, Peruani F. 2020 Self-propelled rods: insights and perspectives for active matter. Annu. Rev. Condens. Matter Phys. 11 , 441–466. (10.1146/annurev-conmatphys-031119-050611)
20. Chaté H. 2020 Dry aligning dilute active matter. Annu. Rev. Condens. Matter Phys. 11 , 189–212. (10.1146/annurev-conmatphys-031119-050752)
21. Méhes E, Vicsek T. 2014 Collective motion of cells: from experiments to models. Integr. Biol. 6 , 831–854. (10.1039/c4ib00115j)
22. Vicsek T, Czirók A, Ben-Jacob E, Cohen I, Shochet O. 1995 Novel type of phase transition in a system of self-driven particles. Phys. Rev. Lett. 75 , 1226–1229. (10.1103/PhysRevLett.75.1226)10060237
23. Tadmor E, Tan C. 2014 Critical thresholds in flocking hydrodynamics with non-local alignment. Philos. Trans. A Math. Phys. Eng. Sci. 372 , 20130401. (10.1098/rsta.2013.0401)25288813
24. Zumaya M, Larralde H, Aldana M. 2018 Delay in the dispersal of flocks moving in unbounded space using long-range interactions. Sci. Rep. 8 , 15872. (10.1038/s41598-018-34208-x)30367121
25. Mao Z, Li Z, Karniadakis GE. 2019 Nonlocal flocking dynamics: learning the fractional order of pdes from particle simulations. Commun. Appl. Math. Comput. 1 , 597–619. (10.1007/s42967-019-00031-y)
26. King AE, Turner MS. 2021 Non-local interactions in collective motion. R. Soc. Open Sci. 8 , 201536. (10.1098/rsos.201536)33959323
27. Bialek W, Cavagna A, Giardina I, Mora T, Silvestri E, Viale M, Walczak AM. 2012 Statistical mechanics for natural flocks of birds. Proc. Natl Acad. Sci. USA 109 , 4786–4791. (10.1073/pnas.1118633109)22427355
28. Bastien R, Romanczuk P. 2020 A model of collective behavior based purely on vision. Sci. Adv. 6 , eaay0792. (10.1126/sciadv.aay0792)32076645
29. Krongauz DL, Ayali A, Kaminka GA. 2024 Vision-based collective motion: a locust-inspired reductionist model. PLoS Comput. Biol. 20 , e1011796. (10.1371/journal.pcbi.1011796)38285716
30. Strandburg-Peshkin A et al . 2013 Visual sensory networks and effective information transfer in animal groups. Curr. Biol. 23 , R709–R711. (10.1016/j.cub.2013.07.059)24028946
31. Camperi M, Cavagna A, Giardina I, Parisi G, Silvestri E. 2012 Spatially balanced topological interaction grants optimal cohesion in flocking models. Interface Focus 2 , 715–725. (10.1098/rsfs.2012.0026)24312725
32. Kumar V, De R. 2021 Efficient flocking: metric versus topological interactions. R. Soc. Open Sci. 8 , 202158. (10.1098/rsos.202158)34631117
33. Forman R. 2003 Bochner’s method for cell complexes and combinatorial ricci curvature. Discrete Comput. Geom. 29 , 323–374. (10.1007/s00454-002-0743-x)
34. Roy I, Vijayaraghavan S, Ramaia SJ, Samal A. 2020 Forman-ricci curvature and persistent homology of unweighted complex networks. Chaos. Solitons. Fractals. 140 , 110260. (10.1016/j.chaos.2020.110260)
35. Sreejith R, Mohanraj K, Jost J, Saucan E, Samal A. 2016 Forman curvature for complex networks. J. Stat. Mech. 2016 , 063206. (10.1088/1742-5468/2016/06/063206)
36. Baccini F, Geraci F, Bianconi G. 2022 Weighted simplicial complexes and their representation power of higher-order network data and topology. Phys. Rev. E 106 , 034319. (10.1103/PhysRevE.106.034319)36266916
37. Gupta KV, Beuria J, Behera L. 2024 Characterizing EEG signals of meditative states using persistent homology and hodge spectral entropy. Biomed. Signal Process. Control 89 , 105779. (10.1016/j.bspc.2023.105779)
38. Sreejith R, Jost J, Saucan E, Samal A. 2017 Systematic evaluation of a new combinatorial curvature for complex networks. Chaos. Solitons. Fractals. 101 , 50–67. (10.1016/j.chaos.2017.05.021)
39. Samal A, Sreejith R, Gu J, Liu S, Saucan E, Jost J. 2018 Comparative analysis of two discretizations of ricci curvature for complex networks. Sci. Rep. 8 , 8650. (10.1038/s41598-018-27001-3)29872167
40. Saucan E, Sreejith R, Vivek-Ananth R, Jost J, Samal A. 2019 Discrete ricci curvatures for directed networks. Chaos. Solitons. Fractals. 118 , 347–360. (10.1016/j.chaos.2018.11.031)
41. Beuria J. 2024 Intrinsic geometry of collider observations and forman ricci curvature. Phys. Rev. D 110 , 035023. (10.1103/PhysRevD.110.035023)
42. Weber M, Saucan E, Jost J. 2017 Characterizing complex networks with forman-ricci curvature and associated geometric flows. J. Complex Netw. 5 , 527–550. (10.1093/comnet/cnw030)
43. Cavagna A, Cimarelli A, Giardina I, Orlandi A, Parisi G, Procaccini A, Santagati R, Stefanini F. 2008 New statistical tools for analyzing the structure of animal groups. Math. Biosci. 214 , 32–37. (10.1016/j.mbs.2008.05.006)18586280
44. Okuta R, Unno Y, Nishino D, Hido S, Loomis C. 2017 CuPy: a NumPy-compatible library for NVIDIA GPU calculations. In Proceedings of Workshop on Machine Learning Systems (LearningSys) in the Thirty-First Annual Conference on Neural Information Processing Systems (NIPS).
45. Jbeuria. 2024 Vicsek-ricci. GitHub. https://github.com/jbeuria/vicsek-ricci
46. Jbeuria. 2024 Jbeuria/vicsek-ricci: version accepted in RSOS. Zenodo. https://zenodo.org/records/12788392
