Contact process (mathematics) explained
The contact process is a stochastic process used to model population growth on the set of sites
of a
graph in which occupied sites become vacant at a constant rate, while vacant sites become occupied at a rate proportional to the number of occupied neighboring sites. Therefore, if we denote by
the proportionality constant, each site remains occupied for a random time period which is
exponentially distributed parameter 1 and places descendants at every vacant neighboring site at times of events of a
Poisson process parameter
during this period. All processes are
independent of one another and of the random period of time sites remains occupied. The contact process can also be interpreted as a model for the spread of an infection by thinking of particles as a bacterium spreading over individuals that are positioned at the sites of
, occupied sites correspond to infected individuals, whereas vacant correspond to healthy ones.
The main quantity of interest is the number of particles in the process, say
, in the first interpretation, which corresponds to the number of infected sites in the second one. Therefore, the process
survives whenever the number of particles is positive for all times, which corresponds to the case that there are always infected individuals in the second one. For any infinite graph
there exists a positive and finite critical value
so that if
then survival of the process starting from a finite number of particles occurs with positive probability, while if
their extinction is almost certain. Note that by Latin:
reductio ad absurdum and the
infinite monkey theorem, survival of the process is equivalent to
, as
, whereas extinction is equivalent to
, as
, and therefore, it is natural to ask about the rate at which
when the process survives.
Mathematical definition
If the state of the process at time
is
, then a site
in
is occupied, say by a particle, if
and vacant if
. The contact process is a continuous-time
Markov process with state space
, where
is a finite or countable
graph, usually
, and a special case of an
interacting particle system.More specifically, the dynamics of the basic contact process is defined by the following transition rates: at site
,
0 → 1 atrateλ\sumy:y\simx\xit(y),
where the sum is over all the neighbors
of
in
. This means that each site waits an exponential time with the corresponding rate, and then flips (so 0 becomes 1 and vice versa).
Connection to percolation
The contact process is a stochastic process that is closely connected to percolation theory. Ted Harris (1974) noted that the contact process on
when infections and recoveries can occur only in discrete times
corresponds to one-step-at-a-time bond percolation on the graph obtained by orienting each edge of
in the direction of increasing coordinate-value.
The law of large numbers on the integers
A law of large numbers for the number of particles in the process on the integers informally means that for all large
,
is approximately equal to
for some positive constant
. Harris (1974) proved that, if the process survives, then the rate of growth of
is at most and at least linear in time. A weak law of large numbers (that the process converges in probability) was shown by
Durrett (1980). A few years later, Durrett and Griffeath (1983) improved this to a strong law of large numbers, giving almost sure convergence of the process.
Die out at criticality
Contact processes on all integer lattices die out almost surely at the critical value.[1]
Durrett's conjecture and the central limit theorem
Durrett conjectured in survey papers and lecture notes during the 1980s and early 1990s regarding the central limit theorem for the Harris contact process, viz. that, if the process survives, then for all large
,
equals
and the error equals
multiplied by a (random) error distributed according to a standard
Gaussian distribution.
[2] [3] [4] Durrett's conjecture turned out to be correct for a different value of
as
proved in 2018.
[5] Further reading
- Durrett . Richard. Rick Durrett . 1980 . On the Growth of One Dimensional Contact Processes . The Annals of Probability . 8 . 5 . 890–907 . 10.1214/aop/1176994619. free .
- Durrett . Richard. Rick Durrett . David Griffeath . 1983 . Supercritical Contact Processes on Z . The Annals of Probability . 11 . 1 . 1–15 . 10.1214/aop/1176993655 . free .
- Book: Liggett
, Thomas M.
. Thomas M. Liggett
. Thomas M. Liggett. Interacting Particle Systems . registration. 1985 . Springer Verlag . New York . 978-0-387-96069-2 .
- Thomas M. Liggett, "Stochastic Interacting Systems: Contact, Voter and Exclusion Processes", Springer-Verlag, 1999.
Notes and References
- Bezuidenhout . Carol . Grimmett . Geoffrey . 10.1214/aop/1176990627 . 4 . . 2244329 . 1071804 . 1462–1482 . The critical contact process dies out . 18 . 1990.
- Durrett. Richard. Oriented Percolation in Two Dimensions Number. The Annals of Probability. 12. 4. 1984. 999–1040. 10.1214/aop/1176993140. free.
- Durrett. Richard. Lecture Notes on Particle Systems and Percolation. Wadsworth.
- .Durrett. Richard. The contact process, 1974–1989. Cornell University, Mathematical Sciences Institute..
- Tzioufas. Achillefs. The Central Limit Theorem for Supercritical Oriented Percolation in Two Dimensions. Journal of Statistical Physics. 171. 5. 2018. 802–821. 10.1007/s10955-018-2040-y. 1411.4543. 2018JSP...171..802T. 119174423.