Stochastic epidemic models with varying infectivity and waning immunity

We study an individual-based stochastic epidemic model in which infected individuals become susceptible again following each infection. In contrast to classical compartment models, after each infection, the infectivity is a random function of the time elapsed since one's infection. Similarly, recovered individuals become gradually susceptible after some time according to a random susceptibility function. We study the large population asymptotic behaviour of the model, by proving a functional law of large numbers (FLLN) and investigating the endemic equilibria properties of the limit. The limit depends on the law of the susceptibility random functions but only on the mean infectivity functions. The FLLN is proved by constructing a sequence of i.i.d. auxiliary processes and adapting the approach from the theory of propagation of chaos. The limit is a generalisation of a PDE model introduced by Kermack and McKendrick, and we show how this PDE model can be obtained as a special case of our FLLN limit.% for a particular set of infectivity and susceptibility random functions and initial conditions. For the endemic equilibria, if $ R_0 $ is lower than (or equal to) some threshold, the epidemic does not last forever and eventually disappears from the population, while if $ R_0 $ is larger than this threshold, the epidemic will not disappear and there exists an endemic equilibrium. The value of this threshold turns out to depend on the harmonic mean of the susceptibility a long time after an infection, a fact which was not previously known.