Women have become better represented in business, academia, and government over time, yet a dearth of women at the highest levels of leadership remains. Sociologists have attributed the leaky progression of women through professional hierarchies to various cultural and psychological factors, such as self-segregation and bias. Here, we present a minimal mathematical model that reveals the relative role that bias and homophily (self-seeking) may play in the ascension of women through professional hierarchies. Unlike previous models, our novel model predicts that gender parity is not inevitable, and deliberate intervention may be required to achieve gender balance in several fields. To validate the model, we analyze a new database of gender fractionation over time for 16 professional hierarchies. We quantify the degree of homophily and bias in each professional hierarchy, and we propose specific interventions to achieve gender parity more quickly. *

Motivated by rhythms in the olfactory system of the brain, we investigate the synchronization of all-to-all pulse-coupled neuronal oscillators exhibiting various types of mixed-mode oscillations (MMOs) composed of sub-threshold oscillations (STOs) and action potentials (“spikes”). We focus particularly on the impact of the delay in the interaction. In the weak-coupling regime, we reduce the system to a Kuramoto-type equation with non-sinusoidal phase coupling and the associated Fokker-Planck equation. Its linear stability analysis identifies the appearance of various cluster states. Their type depends sensitively on the delay and the width of the pulses. Interestingly, long delays do not imply slow population rhythms, and the number of emerging clusters only loosely depends on the number of STOs. Direct simulations of the oscillator equations reveal that for quantitative agreement of the weak-coupling theory the coupling strength and the noise have to be extremely small. Even moderate noise leads to significant skipping of STO cycles, which can enhance the diffusion coefficient in the Fokker-Planck equation by two orders of magnitude. Introducing an effective diffusion coefficient extends the range of agreement significantly. Numerical simulations of the Fokker-Planck equation reveal bistability and solutions with oscillatory order parameters that result from nonlinear mode interactions. These are confirmed in simulations of the full spiking model.