What $R_0$ Deletes: Eigenvectors, Non-Normality, and the Social Content of the Basic Reproduction Number

Abstract: The basic reproduction number is the spectral radius of a matrix, $R_0=ρ(K)$. Taking that definition literally, we ask what $K\mapstoρ(K)$ discards. A matrix carries three kinds of information: its dominant eigenvalue, its dominant eigenvectors, and its departure from normality. $R_0$ keeps only the first; the other two are where the epidemic's social structure lives. The right eigenvector is the burden distribution, the left the source distribution; they coincide when the system is normal and diverge under heterogeneity. Across the $177$ national contact matrices of Prem et al., the operator is \emph{never} normal, and once age-specific susceptibility is included, its source and burden eigenvectors are misaligned by a median of $26<sup>{\circ}$, exceeding $40<sup>{\circ}$ in some countries: the groups that drive transmission are systematically not those that bear it. We prove that under reciprocal contact this misalignment obeys a Kantorovich bound set by the susceptibility contrast $q_{\max}/q_{\min}$ alone, and zero when susceptibility is uniform, with the excess in real, non-reciprocal matrices contributed by contact asymmetry. Transient amplification, by contrast, stays small, so the operative social content is the misalignment, not transient blow-up. The omission also has teeth: because minimizing $R_0$ protects those who \emph{spread} infection, while minimizing deaths protects those who \emph{die} from it, the two target different age groups; the former sometimes raises average infection fatality even as it lowers the scalar. When contact is strongly structured and susceptibility is heterogeneous, we suggest reporting $R_0$ along with its eigenvectors rather than reporting it alone.

• The paper demonstrates how summarizing next-generation matrices to a scalar R0 omits critical eigenvector and non-normality details.
• It introduces a matrix factorization (K = QCD) that reveals the misalignment between groups driving transmission and those bearing its burden.
• Empirical analyses using COVID-19 data from 177 countries illustrate that relying solely on R0 can lead to inequitable and suboptimal intervention policies.

Eigenvectors, Non-Normality, and the Limits of $R_0$: An Analytical Perspective

Overview and Motivation

The scalar basic reproduction number $R_0$ remains central to mathematical epidemiology, widely interpreted as a threshold parameter for population invasibility by infectious diseases. Conventionally, in group-structured or heterogeneous models, $R_0$ is defined as the spectral radius $\rho(K)$ of a next-generation matrix $K$, which encodes expected transmissions across groups. However, this approach compresses the operator $K$ to a single eigenvalue, potentially eliding important structural features. This paper interrogates the "deletions" implicit in $K \mapsto \rho(K)$: specifically, it examines the epidemiological meaning and policy ramifications of the dominant eigenvectors and the non-normality of $K$, both of which are erased in scalar $R_0$ summaries.

Mathematical Structure: Decomposition and Non-Normality

The paper presents a matrix factorization of the next-generation operator,

$K = Q C D,$

where $R_0$0 is the contact (mixing) matrix, $R_0$1 encodes group-specific biological susceptibility, and $R_0$2 captures group-specific infectious durations. Consistent with canonical epidemiological frameworks, the dominant right eigenvector $R_0$3 of $R_0$4 represents the asymptotic incidence distribution, while the left eigenvector $R_0$5 captures per capita reproductive contribution. A critical finding is that almost all real-world $R_0$6 operators are non-normal, even when inputs like $R_0$7 are symmetric, due to heterogeneity in $R_0$8 and $R_0$9. This non-normality is the mathematical signature of heterogeneity itself and provokes divergence between source (transmission-driving) and burden (receiving) groups.

Figure 1: Non-normality and eigenvector misalignment in real transmission operators—visualizing $R_0$0 for Italy, distributions of $R_0$1 and $R_0$2 by age, and decomposition of misalignments across 177 countries.

Eigenvector Misalignment: Quantifying Social Content

The divergence between $R_0$3 (burden) and $R_0$4 (source) eigenvectors is formalized by a condition number

$R_0$5

and its geometric interpretation, the source–burden angle $R_0$6. When $R_0$7 is normal and homogeneous, $R_0$8 and $R_0$9; in heterogeneous, non-normal contexts, $\rho(K)$0. Empirically, the paper demonstrates—using 177 national contact matrices and COVID-19 age-susceptibility profiles—that eigenvector misalignment is substantial: the median angle in real data is $\rho(K)$1 (with maxima exceeding $\rho(K)$2). Crucially, the groups that drive transmission are systematically distinct from those bearing its burden, a phenomenon invisible to $\rho(K)$3 alone.

As shown in Figure 1b, the misalignment produces scenarios where, for example, older age groups incur disproportionately high infection burden relative to their reproductive contribution, evidencing a fundamental equity issue in scalar summaries.

Theoretical Bounds: The Kantorovich Inequality

A principal theoretical result is the derivation of a Kantorovich-type upper bound on the source–burden angle in terms of susceptibility contrast $\rho(K)$4. For symmetrized (reciprocal) contact matrices,

$\rho(K)$5

with equality when heterogeneity is maximized as a two-point contrast. Empirical misalignments from real data all respect this bound under symmetrized contacts; the additional observed misalignment is attributed to asymmetry in $\rho(K)$6. This decomposition allows for an attribution of misalignment to susceptibility heterogeneity versus contact asymmetry.

Non-Normality and Transient Dynamics

Non-normality also predicates the possibility of transient amplification—the scenario where subcritical $\rho(K)$7 yields short-lived outbreaks that briefly exceed expected growth rates. While theoretically possible (as demonstrated by artificial matrix constructions in the text), empirical results show that typical human contact structures exhibit limited reactivity (reactivity ratios $\rho(K)$8 across all countries), indicating that transient amplification is not the dominant empirical consequence of non-normality. Rather, it is the eigenvector misalignment that expresses the operative non-normal effect in actual populations.

Policy Implications: Intervention Targeting and Tradeoffs

A striking practical implication is that optimal intervention targeting diverges depending on the criterion—transmission minimization ($\rho(K)$9-optimal) versus mortality minimization. Analytically, the sensitivity of $K$0 to interventions decreasing susceptibility in group $K$1 is proportional to $K$2, reflecting both the reproductive value and incidence in that group. By contrast, mortality minimization is governed by $K$3 weighted by infection fatality rate, often maximized in entirely different demographic groups. Across all tested countries, the age group to shield is always younger when minimizing $K$4 (typically $K$5–$K$6) and older when minimizing mortality risk (consistently $K$7). The $K$8-optimal intervention frequently results in a 1.5–2.0× increase in average per-infection fatality, sharply illustrating the real social cost of interventions based solely on the spectral radius.

Reporting Recommendations and Theoretical Implications

The core theoretical recommendation is that reporting $K$9 in isolation conceals significant structural and social content. The paper recommends the routine joint reporting of:

• $K$0 (invasion threshold),
• $K$1 and $K$2, or $K$3 (burden/source distributions and their misalignment),
• $K$4 (maximum transient amplification risk).

Since these quantities are straightforward to extract in the computation of $K$5, such reporting incurs essentially zero additional computational cost. Practically, this would support more informed policy decisions in heterogeneously structured populations, where homogeneity-based $K$6 loses vital information.

Conclusions and Future Directions

$K$7 is a meaningful but coarse summary statistic that preserves only threshold information while deleting the core details about the social allocation of harm and the heterogeneity of transmission structure. The observed empirical misalignments and documented policy pitfalls demonstrate that detailed spectral analysis is essential for structured populations or pathogens with age- or group-specific differential susceptibility. The analytical framework of this paper, coupled with large-scale empirical validation, opens further lines of research: extending to non-age group structure, integrating more dimensions of heterogeneity (occupation, geography, socioeconomics), and generalizing the Kantorovich control for broader operator classes. The unambiguous implication is that future models and public health reports should pair $K$8 with its eigenstructure for both theoretical clarity and practical equity.