On the Numerical Computation of R0 in Periodic Environments

We propose a novel approach to approximate the basic reproduction number R0 as the spectral radius of the Next-Generation Operator in time-periodic population models by characterizing the latter via evolution semigroups. Once birth/infection and transition operators are identified, we discretize them via either Fourier or Chebyshev collocation methods. Then R0 is obtained by solving a generalized matrix eigenvalue problem. The order of convergence of the approximating reproduction numbers to the true one is shown to depend on the regularity of the model coefficients, and spectral accuracy is proved. We validate the theoretical results by discussing applications to epidemiology, viz. a large-size multigroup epidemic model with periodic contact rates, and a vector-borne disease model with seasonal vector recruitment. We illustrate how the method facilitates implementation compared to existing approaches and how it can be easily adapted to also compute type-reproduction numbers.