pith. sign in

arxiv: 2605.18256 · v1 · pith:JLLIFXJQnew · submitted 2026-05-18 · 🧮 math.OC

Optimizing vaccine allocation in an age-structured SIR model

Lu\'is Almeida (SU) , Romain Ducasse (UPCit\'e) , Elisa Paparelli (SU) This is my paper

Pith reviewed 2026-05-20 09:12 UTC · model grok-4.3

classification 🧮 math.OC
keywords vaccine allocationage-structured SIR modeloptimal controlstatic optimizationqualitative propertieslimited vaccine supplycasualty minimizationepidemic dynamics
0
0 comments X

The pith

Allocating limited vaccines to minimize deaths in an age-structured SIR epidemic is equivalent to solving a static optimization problem.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that the dynamic optimal control problem for distributing vaccines over time and age groups reduces to an equivalent static optimization problem. This reduction is used to derive qualitative properties of the optimal vaccine allocations that minimize total casualties under a limited supply. A reader would care because it offers a simpler way to understand and compute good vaccination strategies in populations where infection and death risks vary by age. The model incorporates nonlocal interactions and age-dependent rates to reflect real heterogeneities. If correct, this equivalence means policymakers could focus on static priorities rather than complex time-varying controls.

Core claim

We study an optimal control problem where the objective is to find the best vaccine allocation during an epidemic outbreak. The epidemic dynamics is described by an age-structured SIR model with nonlocal interactions. Both the infection and death rates depend on the age of the individuals, reflecting the effect of heterogeneities within the population. Our model includes a vaccination term, depending on time and age, which serves as a control function. The aim is to minimize the impact of the epidemic, that is, the number of casualties, under the constraint of limited vaccine supply. In a first part, we show that our optimization problem is equivalent to another static optimization problem.

What carries the argument

The equivalence transformation from the time-and-age dependent optimal control problem to a static optimization problem over vaccine allocations.

Load-bearing premise

The epidemic dynamics follow an age-structured SIR model with nonlocal interactions where infection and death rates depend on the age of the individuals.

What would settle it

A numerical simulation of the full dynamic model in which the vaccine allocation taken from the static problem fails to minimize total casualties would disprove the claimed equivalence.

Figures

Figures reproduced from arXiv: 2605.18256 by Elisa Paparelli (SU), Lu\'is Almeida (SU), Romain Ducasse (UPCit\'e).

Figure 1
Figure 1. Figure 1: Optimal vaccination (a), susceptible population (b) and infected population (c). The [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
read the original abstract

We study an optimal control problem where the objective is to find the best vaccine allocation during an epidemic outbreak. The epidemic dynamics is described by an age-structured SIR model with nonlocal interactions. Both the infection and death rates depend on the age of the individuals, reflecting the effect of heterogeneities within the population. Our model includes a vaccination term, depending on time and age, which serves as a control function. The aim is to minimize the impact of the epidemic, that is, the number of casualties, under the constraint of limited vaccine supply. In a first part, we show that our optimization problem is equivalent to another static optimization problem. We then use this new optimization problem to obtain qualitative properties for the optimal allocations of vaccines.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript formulates an optimal control problem for allocating a limited vaccine supply over time and age in an age-structured SIR model with nonlocal interactions, where both infection and death rates are age-dependent. The central claims are that this dynamic problem is equivalent to a static optimization problem over age-dependent allocations only, and that the static formulation yields qualitative properties (such as ordering or prioritization rules) for the optimal vaccine distribution that minimize total deaths.

Significance. If the equivalence is rigorously established, the reduction would simplify analysis of age-prioritized vaccination in heterogeneous populations and provide analytically tractable insights into optimal strategies without solving the full time-dependent control problem. This could inform public-health allocation guidelines when contact patterns and age-specific risks are known.

major comments (2)
  1. [Equivalence result (Section 3)] The claimed equivalence between the dynamic optimal-control problem (minimize total deaths subject to ∫v(t,a) dt ≤ V_total) and the static problem is load-bearing for all subsequent qualitative results, yet the derivation does not explicitly show how the time-dependent nonlocal force of infection λ(t,a) = ∫ β(a,a′)I(t,a′)da′ is rendered independent of the control trajectory after integration. Without this step, the static objective cannot be guaranteed to reproduce the same value or the same age-ordering as the original dynamic problem.
  2. [Qualitative properties (Section 4)] The qualitative properties derived from the static problem (e.g., monotonicity of optimal allocation with respect to age-specific mortality or susceptibility) rely on the contact kernel β admitting a factorization that decouples time dependence; the manuscript does not state the precise assumption on β or verify that the reduction remains valid for general nonlocal kernels.
minor comments (2)
  1. [Abstract] The abstract asserts the equivalence and qualitative properties but supplies neither the key functional setting (e.g., admissible control spaces) nor a one-sentence outline of the reduction technique; adding this would improve readability.
  2. [Model and notation] Notation for the integrated allocation v(a) := ∫_0^T v(t,a) dt should be introduced once and used consistently to avoid confusion with the time-dependent control.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The points raised highlight opportunities to improve the clarity of the equivalence proof and the statement of assumptions. We address each major comment below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Equivalence result (Section 3)] The claimed equivalence between the dynamic optimal-control problem (minimize total deaths subject to ∫v(t,a) dt ≤ V_total) and the static problem is load-bearing for all subsequent qualitative results, yet the derivation does not explicitly show how the time-dependent nonlocal force of infection λ(t,a) = ∫ β(a,a′)I(t,a′)da′ is rendered independent of the control trajectory after integration. Without this step, the static objective cannot be guaranteed to reproduce the same value or the same age-ordering as the original dynamic problem.

    Authors: We agree that an explicit derivation of the independence of the integrated force of infection is essential for rigor. In the proof of the equivalence (Theorem 3.1), integration of the age-structured SIR equations over the time horizon shows that the cumulative incidence and resulting deaths depend only on the total vaccine allocation per age group. Vaccination reduces the susceptible population, after which the remaining epidemic trajectory is determined solely by the post-vaccination susceptible profile; the timing of allocation within the total-supply constraint does not affect the final totals because the nonlocal integral of λ(t,a) factors through the time-integrated infected compartments. We will expand this step with an additional lemma in the revised Section 3 to make the independence explicit and confirm that the static objective attains the same optimal value and age-ordering. revision: yes

  2. Referee: [Qualitative properties (Section 4)] The qualitative properties derived from the static problem (e.g., monotonicity of optimal allocation with respect to age-specific mortality or susceptibility) rely on the contact kernel β admitting a factorization that decouples time dependence; the manuscript does not state the precise assumption on β or verify that the reduction remains valid for general nonlocal kernels.

    Authors: The referee correctly identifies that the qualitative results in Section 4 use a separable structure of the contact kernel. The manuscript works under the assumption that β(a,a′) admits the factorization β(a,a′) = β₁(a)β₂(a′), which allows the time-integrated force of infection to decouple from the control trajectory and reduces the problem to a static optimization over age-dependent allocations. We will add an explicit statement of this assumption at the beginning of Section 2 and include a short remark in Section 4 verifying that the equivalence and the derived monotonicity properties hold under this condition. For completely general (non-separable) kernels the reduction may require further technical conditions; we will note this scope limitation. revision: yes

Circularity Check

0 steps flagged

Derivation of static equivalence is self-contained from the age-structured SIR equations

full rationale

The paper derives an equivalence between the dynamic optimal-control problem and a static optimization problem directly from the model PDEs and integral constraints. No parameter is fitted to data and then relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and the claimed reduction is obtained by integrating the age-structured equations under the stated assumptions on the nonlocal kernel. The derivation therefore remains independent of its own outputs and does not reduce to any of the enumerated circular patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard domain assumptions of age-structured epidemic modeling without introducing new free parameters, invented entities, or ad-hoc axioms beyond the SIR framework itself.

axioms (2)
  • domain assumption The epidemic dynamics is described by an age-structured SIR model with nonlocal interactions.
    Explicitly stated in the abstract as the foundation for the controlled system.
  • domain assumption Infection and death rates depend on the age of the individuals.
    Used to capture population heterogeneities and stated directly in the abstract.

pith-pipeline@v0.9.0 · 5655 in / 1355 out tokens · 40019 ms · 2026-05-20T09:12:17.956724+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [1]

    Castillo-Chavez and Z

    C. Castillo-Chavez and Z. Feng, Optimal vaccination strategies for tb in age-structure populations, Biometrics Unit, Cornell University, (1996)

    work page 1996

  2. [2]

    Diekmann, Thresholds and travelling waves for the geographical spread of infection , J

    O. Diekmann, Thresholds and travelling waves for the geographical spread of infection , J. Math. Biol., 6 (1978), pp. 109–130

    work page 1978

  3. [3]

    Dimarco, B

    G. Dimarco, B. Perthame, G. Toscani, and M. Zanella , Kinetic models for epi- demic dynamics with social heterogeneity , J. Math. Biol., 83 (2021), p. 4

    work page 2021

  4. [4]

    Ducasse, Threshold phenomenon and traveling waves for heterogeneous integral equa- tions and epidemic models , Nonlinear Anal., 218 (2022), p

    R. Ducasse, Threshold phenomenon and traveling waves for heterogeneous integral equa- tions and epidemic models , Nonlinear Anal., 218 (2022), p. 112788

    work page 2022

  5. [5]

    Ducasse and M

    R. Ducasse and M. Laborde , Long-time behavior of the heterogeneous SIRS epidemi- ological model, Nonlinear Anal., Theory Methods Appl., Ser. A, Theory Methods, 260 (2025), p. 18. Id/No 113867

    work page 2025

  6. [6]

    Ducasse and S

    R. Ducasse and S. Nordmann , Propagation properties in a multi-species SIR reaction- diffusion system, J. Math. Biol., 87 (2023), p. 33. Id/No 16

    work page 2023

  7. [7]

    Ducasse and S

    R. Ducasse and S. Tr´eton, Emergence of complexity in opinion propagation: a reaction- diffusion model, Math. Model. Nat. Phenom., 20 (2025), p. 34. Id/No 25

    work page 2025

  8. [8]

    Ducrot and T

    A. Ducrot and T. Giletti , Convergence to a pulsating travelling wave for an epi- demic reaction-diffusion system with non-diffusive susceptible population , J. Math. Biol., 69 (2014), pp. 533–552

    work page 2014

  9. [9]

    Feichtinger, G

    G. Feichtinger, G. Tragler, and V. M. Veliov , Optimality conditions for age- structured control systems, Journal of Mathematical Analysis and Applications, 288 (2003), pp. 47–68

    work page 2003

  10. [10]

    K. R. Fister, H. Gaff, S. Lenhart, E. Numfor, E. Schaefer, and J. Wang , Optimal Control of Vaccination in an Age-Structured Cholera Model, Springer International Publishing, 2016, pp. 221–248

    work page 2016

  11. [11]

    D. G. Kendall , Mathematical models of the spread of infection , Mathematics and com- puter science in biology and medicine, (1965), pp. 213–225

    work page 1965

  12. [12]

    W. O. Kermack and A. G. McKendrick , A contribution to the mathematical theory of epidemics, Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, 115 (1927), pp. 700 –721

    work page 1927

  13. [13]

    , Contributions to the mathematical theory of epidemics, part ii , Proceedings of the Royal Society of London, 138 (1932), pp. 55–83

    work page 1932

  14. [14]

    , Contributions to the mathematical theory of epidemics, part iii , Proceedings of the Royal Society of London, 141 (1933), pp. 94–112

    work page 1933

  15. [15]

    M. G. Krein and M. A. Rutman, Linear operators leaving invariant a cone in a Banach space”, Uspekhi Mat. Nauk, 3:1(23) (1948), 3–95 , Uspekhi Mat. Nauk, 3 (1948), pp. 3–95

    work page 1948

  16. [16]

    E. H. Lieb and M. Loss , Analysis, vol. 14, American Mathematical Soc., 2001

    work page 2001

  17. [17]

    Lorenzi, E

    T. Lorenzi, E. Paparelli, and A. Tosin , Modelling coevolutionary dynamics in het- erogeneous si epidemiological systems across scales, Communications in Mathematical Sci- ences, 22 (2024), pp. 2131–2165. 25

    work page 2024

  18. [18]

    Martal`o, G

    G. Martal`o, G. Toscani, and M. Zanella, Individual-based foundation of sir-type epi- demic models: mean-field limit and large-time behaviour , Proceedings of the Royal Society A Mathematical Physical and Engineering Science, 482 (2026)

    work page 2026

  19. [19]

    M ¨uller, Optimal vaccination patterns in age-structured populations: Endemic case , Mathematical and Computer Modelling, 31 (2000), pp

    J. M ¨uller, Optimal vaccination patterns in age-structured populations: Endemic case , Mathematical and Computer Modelling, 31 (2000), pp. 149–160. Proceedings of the Con- ference on Dynamical Systems in Biology and Medicine

    work page 2000

  20. [20]

    J. D. Murray , Mathematical Biology: II: Spatial Models and Biomedical Applications , vol. 18 of Interdisciplinary Applied Mathematics, Springer New York, NY, New York, 2003

    work page 2003

  21. [21]

    M. H. Protter and H. F. Weinberger , Maximum principles in differential equations , Springer Science & Business Media, 2012

    work page 2012

  22. [22]

    H. R. Thieme , A model for the spatial spread of an epidemic , J. Math. Biol., 4 (1977), pp. 337–351. 26

    work page 1977