Nonlinear dynamical systems describe neural activity at various scales and are frequently used to study brain functions and the impact of external perturbations. Here, we explore methods from optimal control theory (OCT) to study efficient, stimulating “control” signals designed to make the neural activity match desired targets. Efficiency is quantified by a cost functional, which trades control strength against closeness to the target activity. Pontryagin’s principle then enables to compute the cost-minimizing control signal. We then apply OCT to a Wilson–Cowan model of coupled excitatory and inhibitory neural populations. The model exhibits an oscillatory regime, low- and high-activity fixed points, and a bistable regime where low- and high-activity states coexist. We compute an optimal control for a state-switching (bistable regime) and a phase-shifting task (oscillatory regime) and allow for a finite transition period before penalizing the deviation from the target state. For the state-switching task, pulses of limited input strength push the activity minimally into the target basin of attraction. Pulse shapes do not change qualitatively when varying the duration of the transition period. For the phase-shifting task, periodic control signals cover the whole transition period. Amplitudes decrease when transition periods are extended, and their shapes are related to the phase sensitivity profile of the model to pulsed perturbations. Penalizing control strength via the integrated 1-norm yields control inputs targeting only one population for both tasks. Whether control inputs drive the excitatory or inhibitory population depends on the state-space location.

We apply the framework of nonlinear optimal control to a biophysically realistic neural mass model, which consists of two mutually coupled populations of deterministic excitatory and inhibitory neurons. External control signals are realized by time-dependent inputs to both populations. Optimality is defined by two alternative cost functions that trade the deviation of the controlled variable from its target value against the “strength” of the control, which is quantified by the integrated 1- and 2-norms of the control signal. We focus on a bistable region in state space where one low- (“down state”) and one high-activity (“up state”) stable fixed points coexist. With methods of nonlinear optimal control, we search for the most cost-efficient control function to switch between both activity states. For a broad range of parameters, we find that cost-efficient control strategies consist of a pulse of finite duration to push the state variables only minimally into the basin of attraction of the target state. This strategy only breaks down once we impose time constraints that force the system to switch on a time scale comparable to the duration of the control pulse. Penalizing control strength via the integrated 1-norm (2-norm) yields control inputs targeting one or both populations. However, whether control inputs to the excitatory or the inhibitory population dominate, depends on the location in state space relative to the bifurcation lines. Our study highlights the applicability of nonlinear optimal control to understand neuronal processing under constraints better.