The Black Death Anomaly: A Non-Abelian Field Theory of Epidemiological Safe Zones

Abstract

Classical reaction-diffusion models of the 14th-century Black Death fail to explain the rapid genetic radiation of \textit{Yersinia pestis} and the anomalous emergence of vast, untouched geographic safe zones, such as Central Europe. In this work, we resolve these historical anomalies by embedding macroscopic pathogen dynamics within a non-Abelian gauge theory. Utilizing the Doi-Peliti formalism, we map the stochastic master equation of a multi-strain epidemic into a covariant classical field theory. We introduce an $SU(N)$ environmental gauge field, $\mathbf{A}_μ$, which actively couples geographic displacement to phenotypic mutation, treating evolutionary drift as a spatial transport phenomenon. We demonstrate via linear stability analysis that this covariant advection drives a Differential Flow (Turing-Hopf) instability, spontaneously breaking spatial symmetry to generate traveling waves of mutation. Furthermore, by extending the pathogen multiplet to the large-$N$ ('t Hooft) continuum limit, we prove that historical safe zones are not statistical outliers nor the result of perfect quarantine, but are mathematically necessary topological voids. In this continuous limit, the destructive interference of the mutating wavefronts analytically resolves into a stable, isotropic macroscopic node governed by a zeroth-order Bessel function ($J_0$), precisely mapping onto the historical survival of Poland and Bohemia.

Figures (2)

• Figure 1: a Turing-Hopf instability driven by a uniform SU(2) gauge field. It presents three heatmaps across a 2D space (x and y axes from 0 to 100) showing the spatial distribution and spontaneous pattern formation for three populations. Variant A ($I_A$) (left, in red tones) shows regions of spontaneous wave formation, with density varying between approximately 0.06097 and 0.06101. Variant B ($I_B$) (middle, in blue tones) similarly displays a wave-like pattern, complementary to Variant A, with density values ranging from about 0.06099 to 0.06103. The Susceptible Population ($S$) (right, in viridis tones) exhibits distinct regions of high and low susceptibility, varying narrowly around 0.5 (specifically from $0.5 - 3\times 10^{-6}$ to $0.5 + 3\times 10^{-6}$). The areas of higher susceptibility are labeled as spatially generated "safe zones," emerging from the complex dynamics between the infected variants.
• Figure 2: This map visualizes the precise, steady-state spatial topology of host survival probability, derived analytically in the preceding section for the infinitely diverse mutational limit. The vibrant colormap represents the normalized intensity—or squared macroscopic amplitude—of the zeroth-order Bessel function of the first kind: $S(\mathbf{r}) \propto |J_0(k |\mathbf{r} - \mathbf{r}_c|)|^2$. The concentric diffraction pattern is explicitly centered on the historical haven coordinates of Southern Poland ($\mathbf{r}_c$, marked by a red star). Physically, this translates to an undepleted host population that is topologically protected from infection.Deep purple represents the concentric diffraction 'dead zone' ripples.