Chad M. Topaz

A Primer of Swarm Equilibria

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 in nature are observed to consist of a concentrated population on the ground separated from an airborne group. Our model can reproduce this configuration; quasi-two-dimensionality of the model plays a critical role.

A model for rolling swarms of locusts

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, the group approaches a crystalline lattice of fixed density if it is H-stable, and in contrast becomes ever denser 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.

Locust dynamics: behavioral phase change and swarming.

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 aggregations structure and reveal transiently traveling clumps of gregarious insects. Our predictions of aggregation and mass gregarization suggest several possible future biological experiments.

Gender Representation on Journal Editorial Boards in the Mathematical Sciences

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.

Flipped Calculus: A Study of Student Performance and Perceptions

Abstract Flipping the classroom refers to moving lectures outside of the classroom to incorporate other activities into a class during its standard meeting time. This pedagogical modality has recently gained traction as a way to center the learning on students in mathematics classrooms. In an effort to better understand the efficacy of this approach, we implemented a controlled study at a small liberal arts college. We compared two sections of the entry-level course applied multivariable calculus I, with one section taught in a traditional lecture-based format and the other taught as a flipped classroom. During our study, we collected and analyzed data related to student performance, as well as perceptions of the approach and attitude toward mathematics in general. Students in both classes scored similarly on graded components of the course, and the majority of students were comfortable with the format of each section. However, some student perceptions and study habits differed.

Social Aggregation in Pea Aphids: Experiment and Random Walk Modeling

From bird flocks to fish schools and ungulate herds to insect swarms, social biological aggregations are found across the natural world. An ongoing challenge in the mathematical modeling of aggregations is to strengthen the connection between models and biological data by quantifying the rules that individuals follow. We model aggregation of the pea aphid, Acyrthosiphon pisum. Specifically, we conduct experiments to track the motion of aphids walking in a featureless circular arena in order to deduce individual-level rules. We observe that each aphid transitions stochastically between a moving and a stationary state. Moving aphids follow a correlated random walk. The probabilities of motion state transitions, as well as the random walk parameters, depend strongly on distance to an aphids nearest neighbor. For large nearest neighbor distances, when an aphid is essentially isolated, its motion is ballistic with aphids moving faster, turning less, and being less likely to stop. In contrast, for short nearest neighbor distances, aphids move more slowly, turn more, and are more likely to become stationary; this behavior constitutes an aggregation mechanism. From the experimental data, we estimate the state transition probabilities and correlated random walk parameters as a function of nearest neighbor distance. With the individual-level model established, we assess whether it reproduces the macroscopic patterns of movement at the group level. To do so, we consider three distributions, namely distance to nearest neighbor, angle to nearest neighbor, and percentage of population moving at any given time. For each of these three distributions, we compare our experimental data to the output of numerical simulations of our nearest neighbor model, and of a control model in which aphids do not interact socially. Our stochastic, social nearest neighbor model reproduces salient features of the experimental data that are not captured by the control.

Forced patterns near a Turing-Hopf bifurcation.

We study time-periodic forcing of spatially extended patterns near a Turing-Hopf bifurcation point. A symmetry-based normal form analysis yields several predictions, including that (i) weak forcing near the intrinsic Hopf frequency enhances or suppresses the Turing amplitude by an amount that scales quadratically with the forcing strength, and (ii) the strongest effect is seen for forcing that is detuned from the Hopf frequency. To apply our results to specific models, we perform a perturbation analysis on general two-component reaction-diffusion systems, which reveals whether the forcing suppresses or enhances the spatial pattern. For the suppressing case, our results are consistent with features of previous experiments on the chlorine dioxide-iodine-malonic acid chemical reaction. However, we also find examples of the enhancing case, which has not yet been observed in experiment. Numerical simulations verify the predicted dependence on the forcing parameters.

Biological Aggregation Driven by Social and Environmental Factors: A Nonlocal Model and Its Degenerate Cahn--Hilliard Approximation

Biological aggregations such as insect swarms and bird flocks may arise from a combination of social interactions and environmental cues. We focus on nonlocal continuum equations, which are often used to model aggregations, and yet which pose significant analytical and computational challenges. Beginning with a particular nonlocal aggregation model [Topaz et al., Bull. Math. Bio., 2006], we derive the minimal well-posed long-wave approximation, which is a degenerate Cahn-Hilliard equation. Energy minimizers of this reduced, local model retain many salient features of those of the nonlocal model, especially for large populations and away from an aggregations boundaries. Using the Cahn-Hilliard model as a testbed, we investigate the degree to which an external potential modeling food sources can be used to suppress peak population density, which is essential for controlling locust outbreaks. A random distribution of food sources tends to increase peak density above its intrinsic value, while a periodic pattern of food sources can decrease it.