Quantum Epidemic Models

• Quantum epidemic models are mathematical frameworks that extend classical disease dynamics by mapping epidemic states to quantum systems using methods like spin-½ representations.
• They employ techniques such as Ising Hamiltonians, Lindblad equations, and quantum walks to capture both stochastic effects and quantum coherence, enabling analysis of non-Markovian memory and spatial heterogeneity.
• They facilitate innovative epidemic control strategies by leveraging quantum algorithms like quantum annealing and QAOA to optimize intervention policies on complex networks.

Quantum epidemic models are a class of mathematical frameworks that extend classical compartmental models of infectious disease dynamics by leveraging quantum computing, open quantum systems theory, quantum statistical mechanics, and quantum algorithmic paradigms. These models map stochastic, networked, and spatially heterogeneous epidemic processes onto quantum formalisms including Ising Hamiltonians, Lindblad master equations, quantum walks, and quantum optimization frameworks. This unification enables the study of both quantum-coherent and classical stochastic effects in epidemic dynamics, the incorporation of memory and non-Markovian transition kernels, efficient simulation or control of large-scale network epidemics using quantum resources, and the exploration of phenomena inaccessible to classical models.

1. Quantum Representations of Classical Epidemic Processes

Quantum epidemic models originate from correspondences between classical stochastic epidemic dynamics and quantum mechanical systems. For a basic Susceptible-Infectious (SI) model on a network, each node's state (Susceptible or Infectious) can be identified with the state of a spin-½ system. The full epidemic configuration space maps to a Hilbert space of dimension $2^N$ for $N$ nodes, with system evolution generated by an operator akin to an Ising Hamiltonian. Classical SI transitions (e.g., a node becoming infected at rate $\lambda_i = \beta \sum_{j\in I} \kappa_{ji}$) are encoded as raising/lowering operators or as classical–quantum Markov chains via parameterized Hamiltonian dynamics and projective measurements (Wang et al., 2022).

Second-quantized approaches extend the operator formalism further, introducing bosonic or fermionic creation/annihilation operators for each compartment, constructing a Liouvillian (Hamiltonian) generator $H$, and representing the time-dependent state as a Fock space vector $|\Psi(t)\rangle$. This allows one to capture multi-compartment epidemics (e.g., SIR, SEIR, or models with vital dynamics) and derive time-evolution equations for moments and correlation functions in a Schrödinger-like equation (Mondaini, 2015, Mondaini et al., 2021).

2. Quantum Epidemic Dynamics: Coherence, Dissipation, and Memory

Quantum extensions of compartmental models explore the role of coherent quantum effects, environmental dissipation, and non-Markovian memory kernels. Open quantum system approaches describe the population as a subsystem coupled to a large reservoir (healthy population or environmental bath). The time evolution is governed by Lindblad equations, whose structure can accommodate both time-local (Markovian) and convolution (memory) parts. By assigning transition events (infection, recovery, death) to quantum jump operators with possibly nontrivial memory kernels (e.g., Gaussian), one encodes realistic, non-exponential dwell time distributions for compartments (Bagarello et al., 18 Mar 2025). The Lindblad formalism preserves trace and positivity, with physical admissibility conditions enforced via survival probabilities $g_n(t)$.

Quantum many-body versions, as realized in Rydberg atom arrays, model site/spin-resolved epidemic processes with coherent driving (infection-like transitions), spontaneous decay (immunization or recovery), and dephasing (environment-induced noise). Depending on regime, these systems feature phase transitions in the universality class of dynamic percolation (classical limit, strong dephasing) or display qualitatively new cascades of discontinuous transitions associated with multiple coherent wavefronts in the fully quantum regime (Pérez-Espigares et al., 2017).

Quantum stochastic and path-integral formulations, such as the Doi–Peliti and operator-valued master equation frameworks, provide access to the full probability distribution of epidemic observables and thus enable quantification of intrinsic uncertainty in epidemic trajectories as well as systematic closure and uncertainty analysis (Rojas-Venegas et al., 2024).

3. Quantum Epidemic Spreading on Networks and in Space

Quantum epidemic models generalize classical network epidemics by incorporating quantum channel properties, spatial constraints, and nonclassical propagation mechanisms.

In quantum networks, channel-dependent connectivity is encoded via probabilistically weighted adjacency matrices, where link success rates reflect photon attenuation, distance, and quantum channel transmission probabilities. Epidemic threshold calculations use the spectral radius of these weighted matrices, revealing that quantum networks can exhibit higher raw epidemic thresholds than classical graphs due to their inherent sparsity and local connection patterns. However, when controlling for average degree, the fundamental threshold remains unchanged, underscoring that apparent quantum robustness stems from effective sparsity, not quantum correlations in epidemic events (Hou et al., 20 Oct 2025).

Quantum random walk formalisms model spatial transmission using coherently evolving walkers, where interference and ballistic propagation enable super-diffusive and heavy-tailed spatial infection clusters. QRW-based epidemics interpolate between classical diffusive regimes and rapid, long-range superspreading (high $R_0$, non-Gaussian cluster profiles), suggesting applicability to outbreaks exhibiting anomalous spatial spread (Manna et al., 6 Sep 2025).

Quantum field-theoretic and statistical mechanical approaches map SIR/SIS contagion to Doi–Peliti coherent-state path integrals and gauge-mediated field theories. Loop corrections (vacuum polarization) lead to effective space-time fractional epidemic equations governed by parabolic Riesz potentials and Riemann-Liouville time convolutions, capturing Lévy-flight–like transmission and burst dynamics. The effective reproduction number $R_\text{eff}(k,\omega)$ becomes a spectral function, directly bounded by the microscopic spatial cutoff, providing a rigorous foundation for observed heavy-tailed transmission patterns (Bernal-Alvarado et al., 17 Mar 2026).

4. Quantum-inspired and Quantum-native Epidemic Control

Quantum-compliant epidemic control leverages quantum algorithms—primarily quantum annealing (QA) and the Quantum Approximate Optimization Algorithm (QAOA)—to solve combinatorial optimization problems arising in non-pharmaceutical intervention strategies (e.g., mobility restrictions, isolation policies) on epidemics evolving over networks. The epidemic control problem (either SIS or SIR) is encoded as a Quadratic Unconstrained Binary Optimization (QUBO) task by recasting the intervention decisions as binary variables, allowing mapping to an Ising Hamiltonian for solution on quantum hardware.

Benchmarks indicate that quantum solvers reach near-instantaneous solution times with solution quality comparable to advanced classical metaheuristics, enabling real-time rolling-horizon control for large graphs (e.g., up to $M=107$ nodes for Italian provinces) (Volpe et al., 17 Jul 2025, Zino et al., 30 Aug 2025). The main limitations include short control horizons imposed by quadraticity, embedding and noise constraints of quantum devices, and the need for horizon-extension techniques for multi-step interventions.

5. Analytic Solutions, Exact Results, and Fractional Dynamics

Certain quantum epidemic models furnish exact analytic solutions. For example, an SIR model defined on a nonuniform quantum time grid (via $q$-difference calculus) admits closed product-form expressions for $N$0 at each $N$1-time step, preserving positivity, conservation, and convergence to the continuous classical limit as $N$2; the model retains the same threshold at $N$3 (Lemos-Silva et al., 9 Dec 2025).

In one spatial dimension, quantum spin operator formalisms allow for the exact solution of master equations with enhanced connectivity, yielding time-dependent compartment densities and exact $N$4-point correlation functions. There, dual-neighbor infection leads to sharper, earlier peaks and residuals not predicted by mean-field theory (Williams et al., 2011).

Fractional epidemic dynamics emerge rigorously by integrating out quantum host “vacua,” with the resulting space-time fractional integro-differential equations governed by Riesz and Riemann–Liouville operators. These fractional models explain bursty outbreak statistics, long-range correlations, and the breakdown of local screening in real-world epidemics exhibiting anomalous scaling (Bernal-Alvarado et al., 17 Mar 2026).

6. Limits, Physical Interpretation, and Theoretical Insights

Quantum epidemic models provide a bridge between classical stochastic epidemic theory and the machinery of quantum many-body systems, statistical mechanics, and quantum computation.

A central result is the demonstration of structural equivalence between classical epidemic models and quantum (tight-binding, spin chain, or open quantum system) models under specific mappings. In many cases, quantum phenomena—including Rabi oscillations, von Neumann entropy, and even classical analogues of entanglement—can be reproduced within classical statistical frameworks under linear stochastic dynamics. Non-classical features (e.g., Bell inequality violations) require additional structure beyond stochastic linear models (Pomorski, 2020).

Open quantum systems and Doi–Peliti approaches highlight the intrinsic uncertainty in epidemic forecasting due to demographic stochasticity, establishing that maximal forecast variance aligns with epidemic peak and is irreducible even with perfect model and parameter knowledge (Rojas-Venegas et al., 2024).

Quantum deformations of classical epidemic models via Poisson–Hopf algebra embedding introduce systematic corrections modeling environmental effects, memory, and correlation beyond mean-field or Markovian assumptions, resulting in modified drift, altered fluctuations, and richer dynamic regimes (Esen et al., 2020).

7. Extensions, Scalability, and Applications

Quantum epidemic modeling frameworks are extensible to multi-compartmental models (e.g., SIR, SEIR), can encode time-dependent heterogeneities (e.g., lockdowns, contact network rewiring), and can be applied not only to biological epidemics but also to viral propagation in quantum communication networks, financial contagion, and opinion dynamics.

Resource requirements for full quantum simulation scale linearly for sparse networks, enabling practical simulation on near-term quantum devices for modest sizes (Wang et al., 2022). Quantum-resilient communication network design can be quantitatively informed by epidemic-threshold analysis of weighted network models (Hou et al., 20 Oct 2025). Quantum annealers and hybrid quantum-classical solvers are already deployed for real-time intervention policy optimization in large-scale, dynamically evolving systems (Volpe et al., 17 Jul 2025).

Quantum epidemic models thus constitute a mature and multi-faceted research area, unifying stochastic process theory, network science, quantum information, and computational epidemic control into a comprehensive theoretical and algorithmic toolkit.