Schelling's model of segregation

Last updated December 02, 2023

Schelling's model of segregation is an agent-based model developed by economist Thomas Schelling.^[1]^[2] Schelling's model does not include outside factors that place pressure on agents to segregate such as Jim Crow laws in the United States, but Schelling's work does demonstrate that having people with "mild" in-group preference towards their own group could still lead to a highly segregated society via de facto segregation.^[3]^[4]^[5]

Model

The original model is set in an $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 $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 $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 $B$ that matches their group—ignoring empty spaces—is greater than or equal $B_{\textrm {a}}$ . If $B<B_{\textrm {a}}$ then the agent will choose to relocate to a vacant spot where $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 $B_{\textrm {seg}}$ such that $B_{\textrm {a}}<B_{\textrm {seg}}$ leads to a random population configuration and $B_{\textrm {a}}\geq B_{\textrm {seg}}$ leads to a segregated population. The value of $B_{\textrm {seg}}$ was approximately ${\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.

Physical model analogies

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 $n$ can be translated into i,j coordinates of ni,nj. Then the state of whether the agent $m_{n}$ moves to a randomly empty grid cell position or 'remains' is defined by:

$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})}$ .

$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] $R$ can be interpreted as a macrostate whose density $\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.

Broader model considerations

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]

Related Research Articles

<span class="mw-page-title-main">Molecular dynamics</span> Computer simulations to discover and understand chemical properties

Molecular dynamics (MD) is a computer simulation method for analyzing the physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a fixed period of time, giving a view of the dynamic "evolution" of the system. In the most common version, the trajectories of atoms and molecules are determined by numerically solving Newton's equations of motion for a system of interacting particles, where forces between the particles and their potential energies are often calculated using interatomic potentials or molecular mechanical force fields. The method is applied mostly in chemical physics, materials science, and biophysics.

<span class="mw-page-title-main">Crowd simulation</span> Model of movement

Crowd simulation is the process of simulating the movement of a large number of entities or characters. It is commonly used to create virtual scenes for visual media like films and video games, and is also used in crisis training, architecture and urban planning, and evacuation simulation.

In chemical kinetics, a reaction rate constant or reaction rate coefficient is a proportionality constant which quantifies the rate and direction of a chemical reaction by relating it with the concentration of reactants.

In probability and statistics, the Kumaraswamy's double bounded distribution is a family of continuous probability distributions defined on the interval (0,1). It is similar to the Beta distribution, but much simpler to use especially in simulation studies since its probability density function, cumulative distribution function and quantile functions can be expressed in closed form. This distribution was originally proposed by Poondi Kumaraswamy for variables that are lower and upper bounded with a zero-inflation. This was extended to inflations at both extremes [0,1] in later work with S. G. Fletcher.

A dynamical billiard is a dynamical system in which a particle alternates between free motion and specular reflections from a boundary. When the particle hits the boundary it reflects from it without loss of speed. Billiards are Hamiltonian idealizations of the game of billiards, but where the region contained by the boundary can have shapes other than rectangular and even be multidimensional. Dynamical billiards may also be studied on non-Euclidean geometries; indeed, the first studies of billiards established their ergodic motion on surfaces of constant negative curvature. The study of billiards which are kept out of a region, rather than being kept in a region, is known as outer billiard theory.

A satellite galaxy is a smaller companion galaxy that travels on bound orbits within the gravitational potential of a more massive and luminous host galaxy. Satellite galaxies and their constituents are bound to their host galaxy, in the same way that planets within our own solar system are gravitationally bound to the Sun. While most satellite galaxies are dwarf galaxies, satellite galaxies of large galaxy clusters can be much more massive. The Milky Way is orbited by about fifty satellite galaxies, the largest of which is the Large Magellanic Cloud.

<span class="mw-page-title-main">Bond graph</span> Graphical representation of a dynamic system

A bond graph is a graphical representation of a physical dynamic system. It allows the conversion of the system into a state-space representation. It is similar to a block diagram or signal-flow graph, with the major difference that the arcs in bond graphs represent bi-directional exchange of physical energy, while those in block diagrams and signal-flow graphs represent uni-directional flow of information. Bond graphs are multi-energy domain and domain neutral. This means a bond graph can incorporate multiple domains seamlessly.

The Bose–Hubbard model gives a description of the physics of interacting spinless bosons on a lattice. It is closely related to the Hubbard model that originated in solid-state physics as an approximate description of superconducting systems and the motion of electrons between the atoms of a crystalline solid. The model was introduced by Gersch and Knollman in 1963 in the context of granular superconductors. The model rose to prominence in the 1980s after it was found to capture the essence of the superfluid-insulator transition in a way that was much more mathematically tractable than fermionic metal-insulator models.

The material point method (MPM) is a numerical technique used to simulate the behavior of solids, liquids, gases, and any other continuum material. Especially, it is a robust spatial discretization method for simulating multi-phase (solid-fluid-gas) interactions. In the MPM, a continuum body is described by a number of small Lagrangian elements referred to as 'material points'. These material points are surrounded by a background mesh/grid that is used to calculate terms such as the deformation gradient. Unlike other mesh-based methods like the finite element method, finite volume method or finite difference method, the MPM is not a mesh based method and is instead categorized as a meshless/meshfree or continuum-based particle method, examples of which are smoothed particle hydrodynamics and peridynamics. Despite the presence of a background mesh, the MPM does not encounter the drawbacks of mesh-based methods which makes it a promising and powerful tool in computational mechanics.

An important question in statistical mechanics is the dependence of model behaviour on the dimension of the system. The shortcut model was introduced in the course of studying this dependence. The model interpolates between discrete regular lattices of integer dimension.

A phase-field model is a mathematical model for solving interfacial problems. It has mainly been applied to solidification dynamics, but it has also been applied to other situations such as viscous fingering, fracture mechanics, hydrogen embrittlement, and vesicle dynamics.

Sugarscape is a model for artificially intelligent agent-based social simulation following some or all rules presented by Joshua M. Epstein & Robert Axtell in their book Growing Artificial Societies.

The kicked rotator, also spelled as kicked rotor, is a paradigmatic model for both Hamiltonian chaos and quantum chaos. It describes a free rotating stick in an inhomogeneous "gravitation like" field that is periodically switched on in short pulses. The model is described by the Hamiltonian

<span class="mw-page-title-main">Metadynamics</span> Scientific computer simulation method

Metadynamics is a computer simulation method in computational physics, chemistry and biology. It is used to estimate the free energy and other state functions of a system, where ergodicity is hindered by the form of the system's energy landscape. It was first suggested by Alessandro Laio and Michele Parrinello in 2002 and is usually applied within molecular dynamics simulations. MTD closely resembles a number of newer methods such as adaptively biased molecular dynamics, adaptive reaction coordinate forces and local elevation umbrella sampling. More recently, both the original and well-tempered metadynamics were derived in the context of importance sampling and shown to be a special case of the adaptive biasing potential setting. MTD is related to the Wang–Landau sampling.

Interatomic potentials are mathematical functions to calculate the potential energy of a system of atoms with given positions in space. Interatomic potentials are widely used as the physical basis of molecular mechanics and molecular dynamics simulations in computational chemistry, computational physics and computational materials science to explain and predict materials properties. Examples of quantitative properties and qualitative phenomena that are explored with interatomic potentials include lattice parameters, surface energies, interfacial energies, adsorption, cohesion, thermal expansion, and elastic and plastic material behavior, as well as chemical reactions.

Quantum thermodynamics is the study of the relations between two independent physical theories: thermodynamics and quantum mechanics. The two independent theories address the physical phenomena of light and matter. In 1905, Albert Einstein argued that the requirement of consistency between thermodynamics and electromagnetism leads to the conclusion that light is quantized, obtaining the relation $. This paper is the dawn of quantum theory. In a few decades quantum theory became established with an independent set of rules. Currently quantum thermodynamics addresses the emergence of thermodynamic laws from quantum mechanics. It differs from quantum statistical mechanics in the emphasis on dynamical processes out of equilibrium. In addition, there is a quest for the theory to be relevant for a single individual quantum system.$

Coarse-grained modeling, coarse-grained models, aim at simulating the behaviour of complex systems using their coarse-grained (simplified) representation. Coarse-grained models are widely used for molecular modeling of biomolecules at various granularity levels.

Supersymmetric theory of stochastic dynamics or stochastics (STS) is an exact theory of stochastic (partial) differential equations (SDEs), the class of mathematical models with the widest applicability covering, in particular, all continuous time dynamical systems, with and without noise. The main utility of the theory from the physical point of view is a rigorous theoretical explanation of the ubiquitous spontaneous long-range dynamical behavior that manifests itself across disciplines via such phenomena as 1/f, flicker, and crackling noises and the power-law statistics, or Zipf's law, of instantonic processes like earthquakes and neuroavalanches. From the mathematical point of view, STS is interesting because it bridges the two major parts of mathematical physics – the dynamical systems theory and topological field theories. Besides these and related disciplines such as algebraic topology and supersymmetric field theories, STS is also connected with the traditional theory of stochastic differential equations and the theory of pseudo-Hermitian operators.

The harmonic mean p-value(HMP) is a statistical technique for addressing the multiple comparisons problem that controls the strong-sense family-wise error rate (this claim has been disputed). It improves on the power of Bonferroni correction by performing combined tests, i.e. by testing whether groups of p-values are statistically significant, like Fisher's method. However, it avoids the restrictive assumption that the p-values are independent, unlike Fisher's method. Consequently, it controls the false positive rate when tests are dependent, at the expense of less power (i.e. a higher false negative rate) when tests are independent. Besides providing an alternative to approaches such as Bonferroni correction that controls the stringent family-wise error rate, it also provides an alternative to the widely-used Benjamini-Hochberg procedure (BH) for controlling the less-stringent false discovery rate. This is because the power of the HMP to detect significant groups of hypotheses is greater than the power of BH to detect significant individual hypotheses.

In many-body physics, the problem of analytic continuation is that of numerically extracting the spectral density of a Green function given its values on the imaginary axis. It is a necessary post-processing step for calculating dynamical properties of physical systems from quantum Monte Carlo simulations, which often compute Green function values only at imaginary-times or Matsubara frequencies.

References

↑ Thomas C. Schelling (1978) Micromotives and Macrobehavior, Norton. Description, preview.
↑ Schelling, Thomas C. (1971). "Dynamic models of segregation". The Journal of Mathematical Sociology. Informa UK Limited. 1 (2): 143–186. doi:10.1080/0022250x.1971.9989794. ISSN 0022-250X.
1 2 Hatna, Erez; Benenson, Itzhak (2012). "The Schelling Model of Ethnic Residential Dynamics: Beyond the Integrated - Segregated Dichotomy of Patterns". Journal of Artificial Societies and Social Simulation. 15 (1): 6. doi: 10.18564/jasss.1873 . ISSN 1460-7425.
↑ Vinkovic, D.; Kirman, A. (2006-12-06). "A physical analogue of the Schelling model". Proceedings of the National Academy of Sciences. 103 (51): 19261–19265. Bibcode:2006PNAS..10319261V. doi: 10.1073/pnas.0609371103 . ISSN 0027-8424. PMC 1748214 . PMID 17151197.
↑ Zhang, Junfu (2011). "Tipping and residential segregation: a unified Schelling model". Journal of Regional Science. Wiley. 51 (1): 167–193. doi:10.1111/j.1467-9787.2010.00671.x. hdl: 10419/36217 . ISSN 0022-4146. S2CID 17624822.
1 2 Rogers, Tim; McKane, Alan J (2011-07-12). "A unified framework for Schelling's model of segregation". Journal of Statistical Mechanics: Theory and Experiment. IOP Publishing. 2011 (7): P07006. arXiv: 1104.1971 . Bibcode:2011JSMTE..07..006R. doi:10.1088/1742-5468/2011/07/p07006. ISSN 1742-5468. S2CID 14701950.
↑ Stauffer, D.; Solomon, S. (2007). "Ising, Schelling and self-organising segregation". The European Physical Journal B. Springer Science and Business Media LLC. 57 (4): 473–479. arXiv: physics/0701051 . Bibcode:2007EPJB...57..473S. doi:10.1140/epjb/e2007-00181-8. ISSN 1434-6028. S2CID 119519603.
↑ Ódor, Géza (2008). "Self-organizing, two-temperature Ising model describing human segregation". International Journal of Modern Physics C. 19 (3): 393–398. arXiv: 0710.1496 . Bibcode:2008IJMPC..19..393O. doi:10.1142/s0129183108012212. ISSN 0129-1831. S2CID 18016863.
1 2 Mantzaris, Alexander; Marich, John; Halfman, Tristin (2018-08-21). "Examining the Schelling Model Simulation through an Estimation of Its Entropy". Entropy. MDPI AG. 20 (9): 623. Bibcode:2018Entrp..20..623M. doi: 10.3390/e20090623 . ISSN 1099-4300. PMC 7513139 . PMID 33265712.
1 2 Mantzaris, Alexander V. (2020-10-12). "Incorporating a monetary variable into the Schelling model addresses the issue of a decreasing entropy trace". Scientific Reports. Springer Science and Business Media LLC. 10 (1): 17005. Bibcode:2020NatSR..1017005M. doi:10.1038/s41598-020-74125-6. ISSN 2045-2322. PMC 7552411 . PMID 33046767.
↑ Nielsen, Alexander Valentin; Gade, Annika Lund; Juul, Jeppe; Strandkvist, Charlotte (2015-11-04). "Schelling model of cell segregation based only on local information". Physical Review E. American Physical Society (APS). 92 (5): 052705. Bibcode:2015PhRvE..92e2705N. doi:10.1103/physreve.92.052705. ISSN 1539-3755. PMID 26651721.
↑ Hatna, Erez; Benenson, Itzhak, eds. (2012-03-12). "Geosimulation of Income-Based Urban Residential Patterns". Advanced Geo-simulation Models. BENTHAM SCIENCE PUBLISHERS. pp. 111–125. doi:10.2174/978160805222611101010111. ISBN 978-1-60805-222-6.
↑ Benenson, Itzhak; Hatna, Erez; Or, Ehud (2009-04-16). "From Schelling to Spatially Explicit Modeling of Urban Ethnic and Economic Residential Dynamics". Sociological Methods & Research. SAGE Publications. 37 (4): 463–497. doi:10.1177/0049124109334792. ISSN 0049-1241. S2CID 120002044.
↑ Bailey, Kenneth D. (1997). "System entropy analysis". Kybernetes. Emerald. 26 (6/7): 674–688. doi:10.1108/03684929710169852. ISSN 0368-492X.
↑ Chauhan, Ankit; Lenzner, Pascal; Molitor, Louise (2018). Schelling Segregation with Strategic Agents. Proceedings of the 11th International Symposium on Algorithmic Game Theory. pp. 137–149. arXiv: 1806.08713 . doi:10.1007/978-3-319-99660-8_13.
↑ Elkind, Edith; Gan, Jiarui; Igarashi, Ayumi; Suksompong, Warut; Voudouris, Alexandros A. (2019). Schelling Games on Graphs. Proceedings of the 28th International Joint Conference on Artificial Intelligence. pp. 266–272.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Thomas C. Schelling (1978) Micromotives and Macrobehavior, Norton. Description, preview.

[2] Schelling, Thomas C. (1971). "Dynamic models of segregation". The Journal of Mathematical Sociology. Informa UK Limited. 1 (2): 143–186. doi:10.1080/0022250x.1971.9989794. ISSN 0022-250X.

[Hatna2012-3] 1 2 Hatna, Erez; Benenson, Itzhak (2012). "The Schelling Model of Ethnic Residential Dynamics: Beyond the Integrated - Segregated Dichotomy of Patterns". Journal of Artificial Societies and Social Simulation. 15 (1): 6. doi: 10.18564/jasss.1873 . ISSN 1460-7425.

[4] Vinkovic, D.; Kirman, A. (2006-12-06). "A physical analogue of the Schelling model". Proceedings of the National Academy of Sciences. 103 (51): 19261–19265. Bibcode:2006PNAS..10319261V. doi: 10.1073/pnas.0609371103 . ISSN 0027-8424. PMC 1748214 . PMID 17151197.

[5] Zhang, Junfu (2011). "Tipping and residential segregation: a unified Schelling model". Journal of Regional Science. Wiley. 51 (1): 167–193. doi:10.1111/j.1467-9787.2010.00671.x. hdl: 10419/36217 . ISSN 0022-4146. S2CID 17624822.

[Rogers2011-6] 1 2 Rogers, Tim; McKane, Alan J (2011-07-12). "A unified framework for Schelling's model of segregation". Journal of Statistical Mechanics: Theory and Experiment. IOP Publishing. 2011 (7): P07006. arXiv: 1104.1971 . Bibcode:2011JSMTE..07..006R. doi:10.1088/1742-5468/2011/07/p07006. ISSN 1742-5468. S2CID 14701950.

[7] Stauffer, D.; Solomon, S. (2007). "Ising, Schelling and self-organising segregation". The European Physical Journal B. Springer Science and Business Media LLC. 57 (4): 473–479. arXiv: physics/0701051 . Bibcode:2007EPJB...57..473S. doi:10.1140/epjb/e2007-00181-8. ISSN 1434-6028. S2CID 119519603.

[8] Ódor, Géza (2008). "Self-organizing, two-temperature Ising model describing human segregation". International Journal of Modern Physics C. 19 (3): 393–398. arXiv: 0710.1496 . Bibcode:2008IJMPC..19..393O. doi:10.1142/s0129183108012212. ISSN 0129-1831. S2CID 18016863.

[Mantzaris2018-9] 1 2 Mantzaris, Alexander; Marich, John; Halfman, Tristin (2018-08-21). "Examining the Schelling Model Simulation through an Estimation of Its Entropy". Entropy. MDPI AG. 20 (9): 623. Bibcode:2018Entrp..20..623M. doi: 10.3390/e20090623 . ISSN 1099-4300. PMC 7513139 . PMID 33265712.

[Mantzaris2020-10] 1 2 Mantzaris, Alexander V. (2020-10-12). "Incorporating a monetary variable into the Schelling model addresses the issue of a decreasing entropy trace". Scientific Reports. Springer Science and Business Media LLC. 10 (1): 17005. Bibcode:2020NatSR..1017005M. doi:10.1038/s41598-020-74125-6. ISSN 2045-2322. PMC 7552411 . PMID 33046767.

[11] Nielsen, Alexander Valentin; Gade, Annika Lund; Juul, Jeppe; Strandkvist, Charlotte (2015-11-04). "Schelling model of cell segregation based only on local information". Physical Review E. American Physical Society (APS). 92 (5): 052705. Bibcode:2015PhRvE..92e2705N. doi:10.1103/physreve.92.052705. ISSN 1539-3755. PMID 26651721.

[12] Hatna, Erez; Benenson, Itzhak, eds. (2012-03-12). "Geosimulation of Income-Based Urban Residential Patterns". Advanced Geo-simulation Models. BENTHAM SCIENCE PUBLISHERS. pp. 111–125. doi:10.2174/978160805222611101010111. ISBN 978-1-60805-222-6.

[13] Benenson, Itzhak; Hatna, Erez; Or, Ehud (2009-04-16). "From Schelling to Spatially Explicit Modeling of Urban Ethnic and Economic Residential Dynamics". Sociological Methods & Research. SAGE Publications. 37 (4): 463–497. doi:10.1177/0049124109334792. ISSN 0049-1241. S2CID 120002044.

[14] Bailey, Kenneth D. (1997). "System entropy analysis". Kybernetes. Emerald. 26 (6/7): 674–688. doi:10.1108/03684929710169852. ISSN 0368-492X.

[15] Chauhan, Ankit; Lenzner, Pascal; Molitor, Louise (2018). Schelling Segregation with Strategic Agents. Proceedings of the 11th International Symposium on Algorithmic Game Theory. pp. 137–149. arXiv: 1806.08713 . doi:10.1007/978-3-319-99660-8_13.

[16] Elkind, Edith; Gan, Jiarui; Igarashi, Ayumi; Suksompong, Warut; Voudouris, Alexandros A. (2019). Schelling Games on Graphs. Proceedings of the 28th International Joint Conference on Artificial Intelligence. pp. 266–272.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]