Cited by 40Mathematical and Theoretical Epidemiology and Ecology ModelsComplex Network Analysis TechniquesTheoretical and Computational Physics

The DOI matches its Crossref record; 21 of the 21 references it cites were found in a registry or an index; none of the 5 sources its claims cite is checked for support yet.

You can run, you can hide: The epidemiology and statistical mechanics of zombies

Open access

Abstract

We use a popular fictional disease, zombies, in order to introduce techniques used in modern epidemiology modeling, and ideas and techniques used in the numerical study of critical phenomena. We consider variants of zombie models, from fully connected continuous time dynamics to a full scale exact stochastic dynamic simulation of a zombie outbreak on the continental United States. Along the way, we offer a closed form analytical expression for the fully connected differential equation, and demonstrate that the single person per site two dimensional square lattice version of zombies lies in the percolation universality class. We end with a quantitative study of the full scale US outbreak, including the average susceptibility of different geographical regions.

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper uses zombie outbreaks as a fictional model system to illustrate epidemiological modeling and critical phenomena, from a solvable well-mixed SZR model to stochastic lattice simulations and a continental US-scale outbreak. It shows the 2D single-person-per-site zombie model lies in the percolation universality class and provides a quantitative US zombie-risk simulation.

Key findings

  • The well-mixed SZR model admits a closed-form analytical solution; the final state depends on the sign of \(P = Z_0 + (1-\alpha)S_0\), with zombies winning for \(\alpha < 1\) and humans winning for \(\alpha > 1\).
  • In stochastic simulations, the extinction probability for a supercritical outbreak is \(P_{\text{ext}} = \alpha\); for \(\alpha = 0.6\), population 200, and one initial infected, zombies won only 40% of runs.
  • On a 2D square lattice with one person per site, the zombie model has critical point \(\alpha_c = 0.43734613(57)\), follows percolation scaling with \(\tau = 187/91\), and has fractal dimension \(d_f = 1.89(4)\), compared with the exact percolation value \(91/48 = 1.895833\).
  • In the US-scale simulation using 2010 Census data, with \(\alpha = 0.8\), \(\beta = 3.6 \times 10^{-3}\)/h, and \(\mu = 0.0914\)/h, most of the US population becomes zombies within the first week; after 28 days the greatest risk is in northeastern Pennsylvania, while after 7 days risk concentrates around populous metropolitan areas.

Methods

  • Deterministic SZR compartmental ODE model with susceptible, zombie, and removed states, nondimensionalized and solved analytically; compared with a density-dependent SIR model.
  • Stochastic simulations using the Gillespie algorithm and, for the US scale, the Next Reaction Method; branching-process analysis for extinction probability.
  • 2D square-lattice simulations with periodic boundary conditions, finite-size scaling, power-law cluster-size distributions, and box-counting fractal dimension; US simulation on a 1500 × 900 grid built from 2010 Census block data, with an added exposed state and zombie diffusion.

Limitations

  • The paper does not include a formal limitations section. It states that real-world SZR-like dynamics would require refinement because medical personnel are only a small fraction of susceptible hosts, and it notes late-time computational inefficiency from simulating zombie diffusion.
  • Evident from the text: parameters are chosen from popular-culture depictions and simple physical estimates rather than empirical zombie data; the US lattice uses binary closing to add artificial connectivity; humans are assumed not to move; and the outbreak model is explicitly fictional.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1103/physreve.92.052801

Cites

21 of 21 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Verification

Claims the paper makes

0 of 5 cited sources checked · 0 supported · 0 flagged

  • Findingreasoning

    If zombies win in the SZR model, the final zombie population is Z∞ = Z0 + (1 - α)S0.

    “If zombies win we have Z ∞ = Z 0 + (1 -α)S 0 .”
    Quote in the paperSubstantial claim
  • Findingreasoning

    P = Z + (1 - α)S is a constant of the SZR dynamics.

    “so P is a constant of the dynamics.”
    Quote in the paperSubstantial claim
  • Findingown data

    For α = 0.6 with a population of 200 and one infected individual to start, the zombies win only 40% of the time.

    “For the parameter settings used in Figs. 1 and 2, namely α = 0.6 with a population of 200 and one infected individual to start, the zombies win only 40% of the time.”
    Quote in the paperNot a substantial claim
  • Findingown data

    In stochastic dynamics, it is not always guaranteed that a supercritical (α < 1) outbreak will take over the entire susceptible population.

    “Another implication of stochastic dynamics is that it is not always guaranteed that a supercritical (α < 1) outbreak will take over the entire susceptible population.”
    Quote in the paperSubstantial claim
  • Findingreasoning

    If α < 1, the SZR system always flows to a final state composed entirely of zombies and no humans, where P denotes the number of zombies that survive.

    “and the system will always flow to a final state composed of entirely zombies and no humans, where P denotes the number of zombies that survive.”
    Quote in the paperSubstantial claim
  • The quote only partly states this claim. 1
  • It may state two separate claims as one. 15
  • Probably not a substantial claim. 6
  • The quote only partly states this claim.

    Findingreasoning

    For an α < 1 simulation in the limit of large populations, the probability of extinction is P_ext = κ/(β + κ) = α.

    “The probability of extinction is just given by our dimensionless inverse virulence α.”
    Quote in the paperSubstantial claim
  • It may state two separate claims as one.

    Findingreasoning

    Unlike the SZR model, the SIR epidemic is self-limiting and there will always be some susceptibles left at the end of the outbreak, no matter how virulent the disease is.

    “Unlike the SZR model, here we see that no matter how virulent the disease is, the epidemic will be self-limiting, and there will always be some susceptibles left at the end of the outbreak.”
    Quote in the paperSubstantial claim
  • It may state two separate claims as one.

    Findingreasoning

    The SIR model does not admit a closed form analytical solution, but a parametric solution can be found.

    “The SIR model does not admit a closed form analytical solution, but we can find a parametric solution by dividing the first equation by the third, revealing”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Findingown data

    When α is large on a fixed lattice, zombies do not spread very far and are always defeated by their neighboring humans.

    “When α is large, the zombies do not spread very far, always being defeated by their neighboring humans.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim. The paper’s conclusions may not depend on it.

    Findingnone given

    The mean population of the occupied grid sites is 420.

    “The mean population of the occupied grid sites is 420; the median population of an occupied site is 13.”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Findingreasoning

    The relevant virulence parameter μ for the density-dependent SIR model has a population dependence, unlike the SZR α parameter which depends only on disease constants β and κ.

    “Unlike the α parameter for SZR which depends only on our disease constants β,κ, the relevant virulence for the density dependent SIR model (μ) has a population dependence.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim. The paper’s conclusions may not depend on it. It may state two separate claims as one.

    Findingown data

    The paper's fractional fluctuations are less than 0.1%, while the experimental results vary over several orders of magnitude.

    “Notice the small vertical scale: our fractional fluctuations are less than 0.1%, while our experimental results vary over several orders of magnitude.”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Findingown data

    In the early stages of the outbreak, while the population is roughly homogeneous, the zombie plague spreads out in roughly uniform circles, where the speed of infection is tied to local population density.

    “In the early stages of the outbreak, while the population is roughly homogeneous, the zombie plague spreads out in roughly uniform circles, where the speed of the infection is tied to the local population density.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: Humans do not move during the zombie outbreak.
    • The authors state this as an assumption: Zombies behave like a random gas inside a cell, so the probability that a zombie will cross a cell boundary is roughly one-fourth of the zombies within v t of the edge moving across that edge in a small amount of time.
    • The authors state this as an assumption: Zombies undergo random motion inside a cell at 1 ft/s and interact with a human anytime they come within 100 ft.
  • It may state two separate claims as one.

    Findingown data

    After four weeks, much of the United States has fallen, but it takes a very long time for the zombies to diffuse and capture the remaining portions of the United States.

    “After four weeks, much of the United States has fallen, but it takes a very long time for the zombies to diffuse and capture the remaining portions of the United States.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: Humans do not move during the zombie outbreak.
    • The authors state this as an assumption: Zombies behave like a random gas inside a cell, so the probability that a zombie will cross a cell boundary is roughly one-fourth of the zombies within v t of the edge moving across that edge in a small amount of time.
    • The authors state this as an assumption: Zombies undergo random motion inside a cell at 1 ft/s and interact with a human anytime they come within 100 ft.
  • Probably not a substantial claim. The paper’s conclusions may not depend on it.

    Findingown data

    The most populated grid site is downtown New York City, with 299,616 individuals.

    “The most populated grid site is downtown New York City, with 299,616 individuals.”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Findingown data

    There is a transition in the susceptible population near α = 1, corresponding to where the continuous dynamics would show a sharp boundary; in the stochastic model the boundary is blurred due to stochasticity.

    “We can clearly see a transition in the susceptible population near α = 1 corresponding to where our continuous dynamics would show a sharp boundary.Here the boundary is blurred, again due to the stochasticity.”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Findingown data

    For α > α_c (when the zombies lose), the large outbreaks are suppressed below the power law, and for α < α_c (when the zombies often win), the largest outbreaks bend up.

    “Notice that for α > α c (when the zombies lose) the large outbreaks are suppressed below the power law, and for α < α c (when the zombies often win) the largest outbreaks bend up.”
    Quote in the paperSubstantial claim
  • It may state two separate claims as one.

    Interpretationreasoning

    This greater risk is partly because those regions have lots of individuals who could potentially serve as patient zero, and partly due to the rapid spread of zombies in those areas.

    “This is partly because those regions have lots of individuals who could potentially serve as patient zero, and partly due to the rapid spread of zombies in those areas.”
    Quote in the paperSubstantial claim
    • The authors state this as an assumption: Humans do not move during the zombie outbreak.
    • The authors state this as an assumption: Zombies behave like a random gas inside a cell, so the probability that a zombie will cross a cell boundary is roughly one-fourth of the zombies within v t of the edge moving across that edge in a small amount of time.
    • The authors state this as an assumption: Zombies undergo random motion inside a cell at 1 ft/s and interact with a human anytime they come within 100 ft.
  • Probably not a substantial claim. The paper may give no reason why it holds. It may state two separate claims as one.

    Interpretationreasoning

    Medical personnel represent only a small fraction of all susceptible hosts, so a refinement to an SZR-type model would be required to account for their role.

    “(Medical personnel, however, represent only a small fraction of all susceptible hosts, so a refinement to an SZR-type model would be required to account for this.)”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Interpretationnone given

    The SZR model requires that susceptible hosts directly participate in the removal of zombie hosts from the infectious population, leading to runaway outbreaks as susceptible hosts are depleted.

    “As noted, the SZR model requires that susceptible hosts directly participate in the removal of zombie hosts from the infectious population, leading to runaway outbreaks as susceptible hosts are depleted.”
    Quote in the paperNot a substantial claim
  • It may state two separate claims as one.

    Premisereasoning

    For zombies, it is appropriate to model the bite parameter β as a rate per person per unit time, so that a zombie in a larger population is proportionally more effective at infecting and is killed proportionally faster.

    “By claiming that zombies can be modeled by a single bite parameter β that itself is a rate per person per unit time, we are claiming that a zombie in a block with 5000 people would be one hundred times as effective at infecting new zombies as a zombie in a block with fifty people; similarly the zombie in question would be killed one hundred times faster.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: Zombies directly seek out hosts to infect, at which point the human and zombie engage in a duel to the death, making density-dependent interactions appropriate.
  • It may state two separate claims as one.

    Premisecitation

    The SIR model on a two-dimensional lattice with a single person per site falls into the percolation universality class, though details of its cluster growth can differ.

    “The SIR model on a two-dimensional lattice with a single person per site falls into the percolation universality class [19], though details of its cluster growth can differ [20].”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [19]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

    • cites [20]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • It may state two separate claims as one.

    Interpretationreasoning

    The holes of all sizes in the critical cluster reflect proximity to the threshold: the battle between zombies and humans is so evenly matched that one gets an emergent scale invariance in the survival patterns.

    “This reflects the proximity to the threshold: the battle between zombies and humans is so evenly matched, that one gets an emergent scale invariance in the survival patterns.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim.

    Premisecitation

    It is often the case that there are several universality classes with critical exponents close to one another, as in the case of percolation.

    “However, it is often the case that there are several universality classes with critical exponents close to one another as in the case of percolation [22].”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [22]

      Couldn’t be matched to a paper in Crossref, so it can’t be checked.

  • Probably not a substantial claim.

    Premisenone given

    The exact percolation value of the fractal dimension is d_f = 91/48 = 1.895833.

    “compared to the exact percolation value of d f = 91/48 = 1.895833.”
    Quote in the paperSubstantial claim
20 of 20
  • The paper relies on this without support.

    The SZR model's population is fully connected and well-mixed, so any infected individual can infect any susceptible individual with equal probability.

    “Until now, we have considered fully connected, well-mixed populations, where any infected individual can infect any susceptible individual with equal probability.”
  • The paper relies on this without support.

    On the two-dimensional square lattice, infection spreads through nearest neighbor bonds only.

    “The infection spreads through nearest neighbor bonds only.”
  • The paper relies on this without support.

    Zombies directly seek out hosts to infect, at which point the human and zombie engage in a duel to the death, making density-dependent interactions appropriate.

    “Zombies directly seek out hosts to infect, at which point the human and zombie engage in a duel to the (un)death.”
  • The paper relies on this without support.

    Zombies do not fly airplanes.

    “Zombies do not fly airplanes, so our model is closer in spirit to the spread of certain agricultural infestations, where the disease spreads across a lattice of sites along the two-dimensional surface of the Earth (although not in those cases where pathogens are transported long distances by atmospheric currents).”
  • The authors state this as an assumption.

    At criticality, an outbreak of size Bs is formed by what would have been B smaller outbreaks of size s which happened to trigger one another, and these in turn are formed by what would have been B outbreaks of size s/B.

    “Let us oversimplify to say that at criticality an outbreak of size Bs is formed by what would have been B smaller outbreaks of size s which happened to trigger one another, and these in turn are formed by what would have been B outbreaks of size s/B.”
  • The authors state this as an assumption.

    If the probabilities and form of mutual triggering are the same at each scale, then many properties of outbreaks would be the same after rescaling sizes by a factor of B.

    “If the probabilities and form of this mutual triggering is the same at each scale, then it would not surprise us that many properties of the outbreaks would be the same, after rescaling the sizes by a factor of B.”
  • The authors state this as an assumption.

    Humans do not move during the zombie outbreak.

    “We assume that the humans do not move, not only for computational efficiency, but because, as we will see, the zombie outbreaks tend to happen rather quickly, and we expect large transportation networks to shut down in the first days, pinning most people to their homes.”
  • The authors state this as an assumption.

    Zombies behave like a random gas inside a cell, so the probability that a zombie will cross a cell boundary is roughly one-fourth of the zombies within v t of the edge moving across that edge in a small amount of time.

    “we assume that the zombies behave like a random gas inside the cell, so that the probability that a zombie will cross a cell boundary is roughly 1 4 Z L 2 Lv t, that is, one-fourth of the zombies within v t of the edge will move across that edge in a small amount of time.”
  • The authors state this as an assumption.

    Zombies undergo random motion inside a cell at 1 ft/s and interact with a human anytime they come within 100 ft.

    “we similarly assume that the zombies are undergoing random motion inside the cell at 1 ft/s, and they interact with a human anytime they come within 100 ft.”
9 of 9

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

  • Zombie movement equation omits three outgoing directions

    Major

    Reaction (45) moves a zombie at rate μZi to each nearest neighbor. On the stated four-neighbor square lattice, total outgoing movement is 4μZi, whereas equation (41) subtracts only μZi.

  • Figure caption gives the wrong power-law slope

    Minor

    The text derives a cumulative cluster-size law proportional to s^(2−τ), whose log-log slope is 2−τ. Figure 7's caption instead calls the slope 2τ.

  • Exposed-state equation omits new infections

    Minor

    Equation (38) gives the exposed population only a loss term, but reaction (42) moves bitten humans into that state. The equation therefore cannot describe the stated infection reaction.

3 of 3

Citations

Cites

Flagged first, then by citations

Research works printed but not linked: 2

  1. P. Munz, 2009, Infect. Dis. Model. Res. Prog., 4, 133
  2. B. Mota, 2014, Bull. Am. Phys. Soc., 59, 1252P

Cited by

Sorted by citations

Metadata

Authors
Alemi, A. A., Bierbaum, M., Myers, C. R., & Sethna, J. P.
Affiliations
Cornell University (US)
Venue
Physical Review E
Published
2015-11-02
Volume
92
Issue
5
Pages
052801-052801
ISSN
  • 1539-3755 (ISSN-L)
  • 1550-2376
OpenAlex
W2111668503
Cited by
40
References
15
Citations by year
2026: 1 · 2025: 2 · 2024: 5 · 2023: 3 · 2021: 19 · 2020: 4 · 2019: 2 · 2018: 1 · 2017: 1 · 2016: 1 · 2015: 1
FWCI
1.30
Type
journal article
Language
en
Open access
Locations
  • Physical Review E (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • arXiv (Cornell University) (submittedVersion)
  • arXiv (Cornell University)
Funders
  • National Science Foundation (IIS–1247696)
  • Cornell University
Topics
  • Physical Sciences › Mathematics › Applied Mathematics › Mathematical and Theoretical Epidemiology and Ecology Models
  • Physical Sciences › Physics and Astronomy › Statistical and Nonlinear Physics › Complex Network Analysis Techniques
  • Physical Sciences › Physics and Astronomy › Condensed Matter Physics › Theoretical and Computational Physics
Keywords
epidemic modeling, critical phenomena, stochastic epidemic simulation, percolation, susceptibility, square lattice, differential equation model, statistical mechanics
MeSH terms
Disease Outbreaks, Humans, Stochastic Processes, Models, Statistical, Epidemiologic Studies

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-27, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-27 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review