Optimization of External Stimuli for Populations of Theta Neurons via Mean-Field Feedback Control

We study the problem of designing 'robust' external excitations for control and synchronization of an assembly of homotypic harmonic oscillators representing so-called theta neurons. The model of theta neuron (Theta model) captures, in main, the bursting behavior of spiking cells in the brain of biological beings, enduring periodic oscillations of the electric potential in their membrane. Our task is to find an external stimulus (control), which steers all neurons of a given population to their desired phases (i.e., excites/slows down its spiking activity) with the highest probability. Our methodology is the following: The optimization problem at hand is formulated as an optimal mean-field control problem for the local continuity equation in the space of probability measures. To solve this problem numerically, we propose an indirect deterministic descent method based on an exact representation of the increment (infinite-order variation) of the objective functional. We illustrate the modus operandi of the proposed method discuss some aspects of its practical realization and provide some results of numerical experiments.