We construct a continuum model for biological aggregations in which individuals experience long-range social attraction and short range dispersal. For the case of one spatial dimension, we study the steady states analytically and numerically. There exist strongly nonlinear states with compact support and steep edges that correspond to localized biological aggregations, or clumps. These steady state clumps are approached through a dynamic coarsening process. In the limit of large population size, the clumps approach a constant density swarm with abrupt edges. We use energy arguments to understand the nonlinear selection of clump solutions, and to predict the internal density in the large population limit. The energy result holds in higher dimensions, as well, and is demonstrated via numerical simulations in two dimensions.

No abstract

We study equilibrium configurations of swarming biological organisms subject to exogenous and pairwise endogenous forces. Beginning with a discrete dynamical model, we derive a variational description of the corresponding continuum population density. Equilibrium solutions are extrema of an energy functional, and satisfy a Fredholm integral equation. We find conditions for the extrema to be local minimizers, global minimizers, and minimizers with respect to infinitesimal Lagrangian displacements of mass. In one spatial dimension, for a variety of exogenous forces, endogenous forces, and domain configurations, we find exact analytical expressions for the equilibria. These agree closely with numerical simulations of the underlying discrete model.The exact solutions provide a sampling of the wide variety of equilibrium configurations possible within our general swarm modeling framework. The equilibria typically are compactly supported and may contain δ-concentrations or jump discontinuities at the edge of the support. We apply our methods to a model of locust swarms, which are observed in nature to consist of a concentrated population on the ground separated from an airborne group. Our model can reproduce this configuration; quasi-twodimensionality of the model plays a critical role.

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.

We classify and predict the asymptotic dynamics of a class of swarming models. The model consists of a conservation equation in one dimension describing the movement of a population density field. The velocity is found by convolving the density with a kernel describing attractive-repulsive social interactions. The kernel's first moment and its limiting behavior at the origin determine whether the population asymptotically spreads, contracts, or reaches steady-state. For the spreading case, the dynamics approach those of the porous medium equation. The widening, compactly-supported population has edges that behave like traveling waves whose speed, density and slope we calculate. For the contracting case, the dynamics of the cumulative density approach those of Burgers' equation. We derive an analytical upper bound for the finite blow-up time after which the solution forms one or more δ-functions.

We study gender representation on the editorial boards of 435 journals in the mathematical sciences. Women are known to comprise approximately 15% of tenure-stream faculty positions in doctoral-granting mathematical sciences departments in the United States. Compared to this group, we find that 8.9% of the 13067 editorships in our study are held by women. We describe group variations within the editorships by identifying specific journals, subfields, publishers, and countries that significantly exceed or fall short of this average. To enable our study, we develop a semi-automated method for inferring gender that has an estimated accuracy of 97.5%. Our findings provide the first measure of gender distribution on editorial boards in the mathematical sciences, offer insights that suggest future studies in the mathematical sciences, and introduce new methods that enable large-scale studies of gender distribution in other fields.

Locusts exhibit two interconvertible behavioral phases, solitarious and gregarious. While solitarious individuals are repelled from other locusts, gregarious insects are attracted to conspecifics and can form large aggregations such as marching hopper bands. Numerous biological experiments at the individual level have shown how crowding biases conversion towards the gregarious form. To understand the formation of marching locust hopper bands, we study phase change at the collective level, and in a quantitative framework. Specifically, we construct a partial integrodifferential equation model incorporating the interplay between phase change and spatial movement at the individual level in order to predict the dynamics of hopper band formation at the population level. Stability analysis of our model reveals conditions for an outbreak, characterized by a large scale transition to the gregarious phase. A model reduction enables quantification of the temporal dynamics of each phase, of the proportion of the population that will eventually gregarize, and of the time scale for this to occur. Numerical simulations provide descriptions of the aggregation's structure and reveal transiently traveling clumps of gregarious insects. Our predictions of aggregation and mass gregarization suggest several possible future biological experiments.

Abstract. We construct an individual-based kinematic model of rolling migratory locust swarms. The model incorporates social interactions, gravity, wind, and the effect of the impenetrable boundary formed by the ground. We study the model using numerical simulations and tools from statistical mechanics, namely the notion of H-stability. For a free-space swarm (no wind and gravity), as the number of locusts increases, it approaches a crystalline lattice of fixed density if it is H-stable, and in contrast becomes ever more dense if it is catastrophic. Numerical simulations suggest that whether or not a swarm rolls depends on the statistical mechanical properties of the corresponding free-space swarm. For a swarm that is H-stable in free space, gravity causes the group to land and form a crystalline lattice. Wind, in turn, smears the swarm out along the ground until all individuals are stationary. In contrast, for a swarm that is catastrophic in free space, gravity causes the group to land and form a bubble-like shape. In the presence of wind, the swarm migrates with a rolling motion similar to natural locust swarms. The rolling structure is similar to that observed by biologists, and includes a takeoff zone, a landing zone, and a stationary zone where grounded locusts can rest and feed.