Interacting particle system

From HandWiki
Revision as of 21:47, 6 February 2024 by Nautica (talk | contribs) (linkage)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Short description: Type of stochastic process

In probability theory, an interacting particle system (IPS) is a stochastic process [math]\displaystyle{ (X(t))_{t \in \mathbb R^+} }[/math] on some configuration space [math]\displaystyle{ \Omega= S^G }[/math] given by a site space, a countably-infinite-order graph [math]\displaystyle{ G }[/math] and a local state space, a compact metric space [math]\displaystyle{ S }[/math]. More precisely IPS are continuous-time Markov jump processes describing the collective behavior of stochastically interacting components. IPS are the continuous-time analogue of stochastic cellular automata.

Among the main examples are the voter model, the contact process, the asymmetric simple exclusion process (ASEP), the Glauber dynamics and in particular the stochastic Ising model.

IPS are usually defined via their Markov generator giving rise to a unique Markov process using Markov semigroups and the Hille-Yosida theorem. The generator again is given via so-called transition rates [math]\displaystyle{ c_\Lambda(\eta,\xi)\gt 0 }[/math] where [math]\displaystyle{ \Lambda\subset G }[/math] is a finite set of sites and [math]\displaystyle{ \eta,\xi\in\Omega }[/math] with [math]\displaystyle{ \eta_i=\xi_i }[/math] for all [math]\displaystyle{ i\notin\Lambda }[/math]. The rates describe exponential waiting times of the process to jump from configuration [math]\displaystyle{ \eta }[/math] into configuration [math]\displaystyle{ \xi }[/math]. More generally the transition rates are given in form of a finite measure [math]\displaystyle{ c_\Lambda(\eta,d\xi) }[/math] on [math]\displaystyle{ S^\Lambda }[/math].

The generator [math]\displaystyle{ L }[/math] of an IPS has the following form. First, the domain of [math]\displaystyle{ L }[/math] is a subset of the space of "observables", that is, the set of real valued continuous functions on the configuration space [math]\displaystyle{ \Omega }[/math]. Then for any observable [math]\displaystyle{ f }[/math] in the domain of [math]\displaystyle{ L }[/math], one has

[math]\displaystyle{ Lf(\eta)=\sum_\Lambda\int_{\xi:\xi_{\Lambda^c}=\eta_{\Lambda^c}}c_\Lambda(\eta,d\xi)[f(\xi)-f(\eta)] }[/math].

For example, for the stochastic Ising model we have [math]\displaystyle{ G=\mathbb Z^d }[/math], [math]\displaystyle{ S=\{-1,+1\} }[/math], [math]\displaystyle{ c_\Lambda=0 }[/math] if [math]\displaystyle{ \Lambda\neq\{i\} }[/math] for some [math]\displaystyle{ i\in G }[/math] and

[math]\displaystyle{ c_i(\eta,\eta^i)=\exp[-\beta\sum_{j:|j-i|=1}\eta_i\eta_j] }[/math]

where [math]\displaystyle{ \eta^i }[/math] is the configuration equal to [math]\displaystyle{ \eta }[/math] except it is flipped at site [math]\displaystyle{ i }[/math]. [math]\displaystyle{ \beta }[/math] is a new parameter modeling the inverse temperature.

The Voter model

Main page: Voter model

The voter model (usually in continuous time, but there are discrete versions as well) is a process similar to the contact process. In this process [math]\displaystyle{ \eta(x) }[/math] is taken to represent a voter's attitude on a particular topic. Voters reconsider their opinions at times distributed according to independent exponential random variables (this gives a Poisson process locally – note that there are in general infinitely many voters so no global Poisson process can be used). At times of reconsideration, a voter chooses one neighbor uniformly from amongst all neighbors and takes that neighbor's opinion. One can generalize the process by allowing the picking of neighbors to be something other than uniform.

Discrete time process

In the discrete time voter model in one dimension, [math]\displaystyle{ \xi_t(x): \mathbb{Z} \to \{0,1\} }[/math] represents the state of particle [math]\displaystyle{ x }[/math] at time [math]\displaystyle{ t }[/math]. Informally each individual is arranged on a line and can "see" other individuals that are within a radius, [math]\displaystyle{ r }[/math]. If more than a certain proportion, [math]\displaystyle{ \theta }[/math] of these people disagree then the individual changes her attitude, otherwise she keeps it the same. Durrett and Steif (1993) and Steif (1994) show that for large radii there is a critical value [math]\displaystyle{ \theta_c }[/math] such that if [math]\displaystyle{ \theta \gt \theta_c }[/math] most individuals never change, and for [math]\displaystyle{ \theta \in (1/2, \theta_c) }[/math] in the limit most sites agree. (Both of these results assume the probability of [math]\displaystyle{ \xi_0(x) = 1 }[/math] is one half.)

This process has a natural generalization to more dimensions, some results for this are discussed in Durrett and Steif (1993).

Continuous time process

The continuous time process is similar in that it imagines each individual has a belief at a time and changes it based on the attitudes of its neighbors. The process is described informally by Liggett (1985, 226), "Periodically (i.e., at independent exponential times), an individual reassesses his view in a rather simple way: he chooses a 'friend' at random with certain probabilities and adopts his position." A model was constructed with this interpretation by Holley and Liggett (1975).

This process is equivalent to a process first suggested by Clifford and Sudbury (1973) where animals are in conflict over territory and are equally matched. A site is selected to be invaded by a neighbor at a given time.

References