This work presents a control strategy for stochastic compartmental epidemic models that operates indirectly by controlling an associated Fokker-Planck partial differential equation. Instead of manipulating individual stochastic sample paths, the framework targets the time evolution of the probability distribution of system states. Steering the distribution allows the designer to shape the ensemble of possible epidemic outcomes and yields robustness to uncertainty in both model dynamics and initial conditions.
The authors position the approach as suitable for settings where uncertainty is intrinsic to transmission dynamics and to initial state estimates. By optimizing at the distributional level, policy objectives may be expressed in terms of statistical or risk-sensitive measures of epidemic outcomes rather than single-run trajectories.
The central formulation is a PDE-constrained optimization problem in which the state equation is the Fokker-Planck equation that governs the probability density of the stochastic compartmental model. Controls act through parameters that enter the underlying stochastic dynamics and, through those dynamics, modify the Fokker-Planck evolution. Cost functionals encode policy goals; different choices represent alternative priorities that a decision maker might have.
This formulation explicitly allows uncertainty in initial data by considering the Fokker-Planck initial condition as a probability distribution. The optimization problem therefore seeks control functions that produce desirable distributional outcomes across a range of initial uncertainties.
The article provides a full mathematical analysis of the PDE-constrained optimization problem. Key theoretical steps include defining the control-to-state mapping that associates each admissible control with the solution of the Fokker-Planck PDE, and analyzing properties of that map to support existence results.
Using this analysis, the authors prove existence of optimal controls for the formulated problem. The paper focuses on rigorous justification of solvability of the constrained optimization, ensuring that optimal control solutions exist under the assumptions laid out in the development of the model and cost functional.
To characterize optimal controls, the study applies the Pontryagin minimum principle in the context of the PDE-constrained problem. This yields necessary conditions for optimality expressed via an optimality system involving the Fokker-Planck state equation and associated adjoint conditions. The characterization provides a theoretical basis for computing and interpreting optimal control laws in practice.
For numerical approximation of optimal controls, the authors describe the application of the sequential quadratic Hamiltonian method to their PDE-constrained optimal control problem. This iterative numerical technique generates approximations of optimal control maps by solving a sequence of quadratic subproblems derived from the Hamiltonian formulation of the control problem.
The manuscript emphasizes that this approach produces numerical approximations of the optimal control maps suitable for implementation in the context of the modeled stochastic epidemic.
To demonstrate the framework, the authors apply it to a minimal stochastic susceptible-infected-recovered (SIR) model. They explore different choices of cost functional that reflect distinct policy-maker concerns. The examples illustrate how varying the cost objective modifies the optimal control strategy computed via the Fokker-Planck PDE approach.
The demonstration is presented as a proof of concept showing the framework’s practical application to a canonical epidemiological model. The source describes these illustrative computations but does not report specific numerical outcomes or comparative performance metrics beyond the methodological demonstration.
Controlling the Fokker-Planck equation shifts the control problem from individual stochastic realizations to distributions, which is a principled way to handle uncertainty in both dynamics and initial conditions. Strengths highlighted by the formulation include robustness to variability and a direct route to optimizing distributional targets important in public-health decision making.
Limitations inherent to the approach—such as the mathematical and computational complexity of PDE-constrained optimization and the need to specify appropriate cost functionals that capture policy objectives—are implicit in the development. The source emphasizes the theoretical existence and characterization of optimal controls and the availability of a numerical method, but does not provide exhaustive performance benchmarks or comparisons with alternative control strategies.
The article "A Fokker-Planck framework for control of epidemics" is authored by Christian Parkinson and Souvik Roy and appears in J Math Biol (2026). Identifiers reported in the source include PMID 42747468 and DOI 10.1007/s00285-026-02462-7. Affiliations listed are the Department of Mathematics & Department of Computational Mathematics, Science and Engineering at Michigan State University and the Department of Mathematics at the University of Texas at Arlington.
The authors declare no conflict of interest. Grant support noted in the source includes Division of Mathematical Sciences awards (2309491 and 2230790). Keywords provided in the publication include Epidemiology, Fokker-Planck equation, Optimal control, PDE-constrained optimization, and the sequential quadratic Hamiltonian method.