“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Fluctuation analysis of rotational speeds of the bacterial flagellar motor.

We measured the dependence of the variance in the rotation rate of tethered cells of Escherichia coli on the mean rotation rate over a regime in which the motor generates constant torque. This dependence was compared with that of broken motors. In either case, motor torque was augmented with externally applied torque. We show that, in contrast to broken motors, functioning motors in this regime do not freely rotationally diffuse and that the variance measurements are consistent with the predicted values of a stepping mechanism with exponentially distributed waiting times (a Poisson stepper) that steps approximately 400 times per revolution.