Schelling's model of segregation

Simulation of the model. Agents will move at each step until the fraction of neighbors that are from their own group is greater than or equal to ${\displaystyle B_{\textrm {a}}}$. For equal sized populations, ${\displaystyle B_{\textrm {a}}\geq B_{\textrm {seg}}\approx 1/3}$ leads to the groups segregating themselves.

The original model is set in an ${\displaystyle N\times N}$ grid. Agents are split into two groups and occupy the spaces of the grid and only one agent can occupy a space at a time. Agents desire a fraction ${\displaystyle B_{\textrm {a}}}$ of their neighborhood (in this case defined to be the eight adjacent agents around them) to be from the same group. Increasing ${\displaystyle B_{\textrm {a}}}$ corresponds to increasing the agent's intolerance of outsiders.

Each round consists of agents checking their neighborhood to see if the fraction of neighbors ${\displaystyle B}$ that matches their group—ignoring empty spaces—is greater than or equal ${\displaystyle B_{\textrm {a}}}$. If ${\displaystyle B<B_{\textrm {a}}}$ then the agent will choose to relocate to a vacant spot where ${\displaystyle B\geq B_{\textrm {a}}}$. This continues until every agent is satisfied. Every agent is not guaranteed to be satisfied and in these cases it is of interest to study the patterns (if any) of the agent dynamics.

While studying populations dynamics of two groups of equal size, Schelling found a threshold ${\displaystyle B_{\textrm {seg}}}$ such that ${\displaystyle B_{\textrm {a}}<B_{\textrm {seg}}}$ leads to a random population configuration and ${\displaystyle B_{\textrm {a}}\geq B_{\textrm {seg}}}$ leads to a segregated population. The value of ${\displaystyle B_{\textrm {seg}}}$ was approximately ${\displaystyle {\frac {1}{3}}}$. This points to how individuals with even a small amount of in-group preference can form segregated societies. There are different parameterizations and variants of the model and a 'unified' approach is presented in [6] allowing the simulations to explore the thresholds for different segregation events to occur.

There have been observations that the fundamental dynamics of the agents resemble the mechanics used in the Ising model of ferromagnetism.[7][8][9][10] This primarily relies on the similar nature in which each occupied grid location calculates an aggregate measure based upon the similarities of the adjacent grid cells. If each agent produces a satisfaction based upon their homophilic satisfaction threshold as ${\textstyle [0,1]}$ then the summation of those values can provide an indication for the segregation of the state that is analogous to the clustering of the aligned spins in a magnetic material. If each cell is a member of a group ${\textstyle m_{n}\in {m_{1},m_{2},m_{empty}}}$, then the local homogeneity can be found via

${\textstyle l(m_{n})=\sum _{i=-1}^{1}\sum _{j=-1}^{1}\left(\delta _{m_{(ni,nj)},m_{(ni+i,nj+j)}}:i,j\neq 0\right)}$ where the 1-d position of ${\displaystyle n}$ can be translated into i,j coordinates of ni,nj. Then the state of whether the agent ${\displaystyle m_{n}}$ moves to a randomly empty grid cell position or 'remains' is defined by:

Agents that have their local homogeneity constraint satisfied 'remain' ${\displaystyle R}$ in that position between iterations. The total ${\displaystyle R}$ for the grid is plotted as an average over 500 simulations.

${\displaystyle r\left(m_{n}\right)={\begin{cases}\left(l\left(m_{n}\right)\geq B_{a}\right),&{\text{if}}:m_{n}\notin {m_{empty}}\\0,&{\text{if}}:m_{n}\in {m_{empty}}\end{cases}}}$

Each agent produces a binary value, so that for each grid configuration of agents of both groups, a vector can be produced of the remain due to satisfaction or not. The overall satisfaction from the remain states of all the agents can be computed;${\textstyle R=\sum _{n=1}^{N}r(m_{n})}$.

${\displaystyle R}$ then provides a measure for the amount of homogeneity (segregation) on the grid and can be used with the maximum possible value (total sum of agents) as a 'density' of segregation over the simulation of movements as is performed in.[6][11] Following the approach of [9] ${\displaystyle R}$ can be interpreted as a macrostate whose density ${\displaystyle \Omega }$ can be estimated by sampling via the Monte Carlo method the grid space from the random initialisations of the grid to produce a calculation of the entropy; ${\textstyle S=k_{B}{\text{ln}}\Omega (R).}$ This allows a trace of the entropy to be computed over the iterations of the simulation as is done with other physical systems.

The canonical Schelling model does not consider variables which may affect the agent's ability to relocate positions in the grid. The work of [3] investigates a model extension where the utility available to agents to move governs this action. It can explain some of the patterns seen where groups do not segregate due to the financial barrier homogeneous zones produce as a result of high demand. The consideration of the financial aspect is also investigated in [12] and.[13] The work of [10] further develops this concept of the importance of the monetary factor in the decision making, and uses it to extend the model with a dual dynamic where agents radiate their income store whenever a movement is made. This also provides a means to produce a more complete model where the trace of the entropy is non-decreasing and adds support that social systems obey the Second law of thermodynamics.[14]

Schelling's model has also been studied from a game-theoretic perspective: In Schelling games, agents strategically strive to maximize their utilities by relocating to a position with the highest fraction of neighboring agents from the same group.[15][16]

• Critical mass
• Domino effect
• Gerrymandering
• Social dynamics
• Tipping point