How optimal allocation of limited testing capacity changes epidemic dynamics.

Insufficient testing capacity has been a critical bottleneck in the worldwide fight against COVID-19. Optimizing the deployment of limited testing resources has therefore emerged as a keystone problem in pandemic response planning. Here, we use a modified SEIR model to optimize testing strategies under a constraint of limited testing capacity. We define pre-symptomatic, asymptomatic, and symptomatic infected classes, and assume that positively tested individuals are immediately moved into quarantine. We further define two types of testing. Clinical testing focuses only on the symptomatic class. Non-clinical testing detects pre- and asymptomatic individuals from the general population, and a concentration parameter governs the degree to which such testing can be focused on high infection risk individuals. We then solve for the optimal mix of clinical and non-clinical testing as a function of both testing capacity and the concentration parameter. We find that purely clinical testing is optimal at very low testing capacities, supporting early guidance to ration tests for the sickest patients. Additionally, we find that a mix of clinical and non-clinical testing becomes optimal as testing capacity increases. At high but empirically observed testing capacities, a mix of clinical testing and non-clinical testing, even if extremely unfocused, becomes optimal. We further highlight the advantages of early implementation of testing programs, and of combining optimized testing with contact reduction interventions such as lockdowns, social distancing, and masking.

Implicit versus explicit vector management strategies in models for vector-borne disease epidemiology.

Throughout the vector-borne disease modeling literature, there exist two general frameworks for incorporating vector management strategies (e.g. area-wide adulticide spraying and larval source reduction campaigns) into vector population models, namely, the "implicit" and "explicit" control frameworks. The more simplistic "implicit" framework facilitates derivation of mathematically rigorous results on disease suppression and optimal control, but the biological connection of these results to real-world "explicit" control actions that could guide specific management actions is vague at best. Here, we formally define a biological and mathematical relationship between implicit and explicit control, and we provide mathematical expressions relating the strength of implicit control to management-relevant properties of explicit control for four common intervention strategies. These expressions allow the optimal control and basic reproduction number analyses typically utilized in implicit control modeling to be interpreted directly in terms of real-world actions and real-world monetary costs. Our methods reveal that only certain sub-classes of explicit control protocols are able to be represented as implicit controls, and that implicit control is a meaningful approximation of explicit control only when resonance-like synergistic effects between multiple explicit controls have negligible effects on population reduction. When non-negligible synergy exists, implicit control results, despite their mathematical tidiness, fail to provide accurate predictions regarding vector control and disease spread. Collectively, these elements build an effective bridge between analytically interesting and mathematically tractable implicit control and the challenging, action-oriented explicit control.

Managing disease outbreaks: The importance of vector mobility and spatially heterogeneous control.

Management strategies for control of vector-borne diseases, for example Zika or dengue, include using larvicide and/or adulticide, either through large-scale application by truck or plane or through door-to-door efforts that require obtaining permission to access private property and spray yards. The efficacy of the latter strategy is highly dependent on the compliance of local residents. Here we develop a model for vector-borne disease transmission between mosquitoes and humans in a neighborhood setting, considering a network of houses connected via nearest-neighbor mosquito movement. We incorporate large-scale application of adulticide via aerial spraying through a uniform increase in vector death rates in all sites, and door-to-door application of larval source reduction and adulticide through a decrease in vector emergence rates and an increase in vector death rates in compliant sites only, where control efficacies are directly connected to real-world experimentally measurable control parameters, application frequencies, and control costs. To develop mechanistic insight into the influence of vector motion and compliance clustering on disease controllability, we determine the basic reproduction number R0 for the system, provide analytic results for the extreme cases of no mosquito movement, infinite hopping rates, and utilize degenerate perturbation theory for the case of slow but non-zero hopping rates. We then determine the application frequencies required for each strategy (alone and combined) in order to reduce R0 to unity, along with the associated costs. Cost-optimal strategies are found to depend strongly on mosquito hopping rates, levels of door-to-door compliance, and spatial clustering of compliant houses, and can include aerial spray alone, door-to-door treatment alone, or a combination of both. The optimization scheme developed here provides a flexible tool for disease management planners which translates modeling results into actionable control advice adaptable to system-specific details.

The relationship between controllability, optimal testing resource allocation, and incubation-latent period mismatch as revealed by COVID-19.

The severe shortfall in testing supplies during the initial COVID-19 outbreak and ensuing struggle to manage the pandemic have affirmed the critical importance of optimal supply-constrained resource allocation strategies for controlling novel disease epidemics. To address the challenge of constrained resource optimization for managing diseases with complications like pre- and asymptomatic transmission, we develop an integro partial differential equation compartmental disease model which incorporates realistic latent, incubation, and infectious period distributions along with limited testing supplies for identifying and quarantining infected individuals. Our model overcomes the limitations of typical ordinary differential equation compartmental models by decoupling symptom status from model compartments to allow a more realistic representation of symptom onset and presymptomatic transmission. To analyze the influence of these realistic features on disease controllability, we find optimal strategies for reducing total infection sizes that allocate limited testing resources between 'clinical' testing, which targets symptomatic individuals, and 'non-clinical' testing, which targets non-symptomatic individuals. We apply our model not only to the original, delta, and omicron COVID-19 variants, but also to generically parameterized disease systems with varying mismatches between latent and incubation period distributions, which permit varying degrees of presymptomatic transmission or symptom onset before infectiousness. We find that factors that decrease controllability generally call for reduced levels of non-clinical testing in optimal strategies, while the relationship between incubation-latent mismatch, controllability, and optimal strategies is complicated. In particular, though greater degrees of presymptomatic transmission reduce disease controllability, they may increase or decrease the role of non-clinical testing in optimal strategies depending on other disease factors like transmissibility and latent period length. Importantly, our model allows a spectrum of diseases to be compared within a consistent framework such that lessons learned from COVID-19 can be transferred to resource constrained scenarios in future emerging epidemics and analyzed for optimality.

Dynamic modelling of personal protection control strategies for vector-borne disease limits the role of diversity amplification.

Personal protection measures, such as bed nets and repellents, are important tools for the suppression of vector-borne diseases like malaria and Zika, and the ability of health agencies to distribute protection and encourage its use plays an important role in the efficacy of community-wide disease management strategies. Recent modelling studies have shown that a counterintuitive diversity-driven amplification in community-wide disease levels can result from a population's partial adoption of personal protection measures, potentially to the detriment of disease management efforts. This finding, however, may overestimate the negative impact of partial personal protection as a result of implicit restrictive model assumptions regarding host compliance, access to and longevity of protection measures. We establish a new modelling methodology for incorporating community-wide personal protection distribution programmes in vector-borne disease systems which flexibly accounts for compliance, access, longevity and control strategies by way of a flow between protected and unprotected populations. Our methodology yields large reductions in the severity and occurrence of amplification effects as compared to existing models.

Universal energy diffusion in a quivering billiard.

We introduce and study a model of time-dependent billiard systems with billiard boundaries undergoing infinitesimal wiggling motions. The so-called quivering billiard is simple to simulate, straightforward to analyze, and is a faithful representation of time-dependent billiards in the limit of small boundary displacements. We assert that when a billiard's wall motion approaches the quivering motion, deterministic particle dynamics become inherently stochastic. Particle ensembles in a quivering billiard are shown to evolve to a universal energy distribution through an energy diffusion process, regardless of the billiard's shape or dimensionality, and as a consequence universally display Fermi acceleration. Our model resolves a known discrepancy between the one-dimensional Fermi-Ulam model and the simplified static wall approximation. We argue that the quivering limit is the true fixed wall limit of the Fermi-Ulam model.