Topological data analysis of biological aggregation models.
Bottom Line:
To interpret our results, we introduce a visualization that displays Betti numbers over simulation time and topological persistence scale.We compare our topological results to order parameters typically used to quantify the global behavior of aggregations, such as polarization and angular momentum.The topological calculations reveal events and structure not captured by the order parameters.
View Article:
PubMed Central - PubMed
Affiliation: Department of Mathematics, Statistics, and Computer Science, Macalester College, Saint Paul, Minnesota, United States of America.
ABSTRACT
We apply tools from topological data analysis to two mathematical models inspired by biological aggregations such as bird flocks, fish schools, and insect swarms. Our data consists of numerical simulation output from the models of Vicsek and D'Orsogna. These models are dynamical systems describing the movement of agents who interact via alignment, attraction, and/or repulsion. Each simulation time frame is a point cloud in position-velocity space. We analyze the topological structure of these point clouds, interpreting the persistent homology by calculating the first few Betti numbers. These Betti numbers count connected components, topological circles, and trapped volumes present in the data. To interpret our results, we introduce a visualization that displays Betti numbers over simulation time and topological persistence scale. We compare our topological results to order parameters typically used to quantify the global behavior of aggregations, such as polarization and angular momentum. The topological calculations reveal events and structure not captured by the order parameters. No MeSH data available. |
Related In:
Results -
Collection
License getmorefigures.php?uid=PMC4430537&req=5
Mentions: Fig 12(A) shows time series of three order parameter metrics used in [12] to characterize the global behavior of the system. Polarization P, defined in (Eq 1), is similar to φ(t) for the Vicsek model; it measures the degree to which agents are aligned. Because the group is traveling around a circle, P (red curve) remains low for the duration of the simulation. Angular momentum M, also defined in (Eq 1), helps quantify the rotation of the group. For a perfect mill structure, M = 1. For our simulation, we see an evolution from M ≈ 0 early in time to M ≈ 0.93 for later times (green curve). However, as mentioned previously, the metric M cannot distinguish between single and double mills, and so [12] introduces the absolute angular momentum Mabs, defined in (Eq 2). The fact that Mabs approaches unity signals that the asymptotic behavior of the group is rotational. The fact that M approaches a number close to but less than unity signals that a small minority of the group members are rotating counter to the majority. |
View Article: PubMed Central - PubMed
Affiliation: Department of Mathematics, Statistics, and Computer Science, Macalester College, Saint Paul, Minnesota, United States of America.
No MeSH data available.