Project 6: Difference between revisions
Line 22: | Line 22: | ||
To build the BTW "sandpile" that self-organizes into a critical state, start from a <math> 100 \times 100 </math> lattice ''far from equilibrium'' where each site has a random value <math> z(i,j) >3 </math>, and update the whole lattice according to the above rules until the sandpile self-organizes into a statistically stationary (time-independent) critical state with average <math> \langle z \rangle_c < 3 </math>. | To build the BTW "sandpile" that self-organizes into a critical state, start from a <math> 100 \times 100 </math> lattice ''far from equilibrium'' where each site has a random value <math> z(i,j) >3 </math>, and update the whole lattice according to the above rules until the sandpile self-organizes into a statistically stationary (time-independent) critical state with average <math> \langle z \rangle_c < 3 </math>. | ||
a) Find the average height | a) Find explicitly the average height <math> \langle z \rangle_c = \frac{1}{10000} \sum_{i,j=1}^{100} z(i,j) </math> in the SOC state. | ||
b) Add a single grain to this configuration at a randomly chosen site. This will trigger an avalanche of any size (up to the system size), which is why these state is critical. Record the size of the avalanche (number of sites that participate in the toppling rules) and time it takes for avalanche to end. If you collect many avalanches you should find power laws for the distribution function of avalanche sizes and times it take for the avalanche to cease: and . Plot these distributions on a log-log scale and from the slope of the line extract the critical exponents and . Although you can check all sites for stability, you can speed up the update by keeping two arrays which contain and coordinates of those sites that need to be checked for stability (each time a site topples, add its neighbors' positions to these arrays; then remove the last site from the arrays and check its stability; continue this process of removing sites and checking their stability, and adding neighbors to the array until there are no more sites to check). | b) Add a single grain to this configuration at a randomly chosen site. This will trigger an avalanche of any size (up to the system size), which is why these state is critical. Record the size of the avalanche (number of sites that participate in the toppling rules) and time it takes for avalanche to end. If you collect many avalanches you should find power laws for the distribution function of avalanche sizes and times it take for the avalanche to cease: and . Plot these distributions on a log-log scale and from the slope of the line extract the critical exponents and . Although you can check all sites for stability, you can speed up the update by keeping two arrays which contain and coordinates of those sites that need to be checked for stability (each time a site topples, add its neighbors' positions to these arrays; then remove the last site from the arrays and check its stability; continue this process of removing sites and checking their stability, and adding neighbors to the array until there are no more sites to check). |
Revision as of 15:30, 30 April 2012
Self-Organized Criticality
Introduction
Sandpile cellular automata were introduced almost twenty years ago as a paradigmatic example of self-organized criticality (SOC) [1,2], the tendency of slowly driven dissipative systems to display a scale free avalanche response. Such ideas have had an enormous impact in different fields, ranging from magnetic systems, superconductors and mechanics, to geophysics and plasma physics, including in particular the magnetosphere.
The influence also extends beyond physics, to for example biology, human (heart) physiology and cognitive processes or neuroscience. The reason for this success lies in the wide variety of non-equilibrium systems displaying an avalanche response to an external driving. One of the primary aims of SOC was, originally [1], to explain the wide occurrence of noise in natural phenomena, through a direct relation between avalanche scaling and spectral properties.
The criticality of this state, obtained without fine tuning of parameters (such as the tuning of temperature in the case of conventional thermal critical phenomena, as explored in Project 4), is manifested as the absence of characteristic spatial and temporal scales of these avalanches. In the less fancy language: large avalanches occur rather often (there is no exponential decay of avalanche sizes, which would result in a characteristic avalanche size), and there is a variety of power laws without cutoffs in various properties of the system.
The paradigm model for this type of behavior is the celebrated sandpile cellular automaton also know as the Bak-Tang-Wiesenfeld (BTW) model. The local rules for the so-called Bak-Tang-Wiesenfeld (BTW) cellular automaton update the integer field defined on a square lattice synchronously according to:
Thus, if there are more than 3 "grains of sand" on a particular lattice site , then that site relaxes sending one grain of sand to each of four neighbors. If the unstable site happens to be at the boundary, the grains of sand simply leave the system (imagine that they fall off the edge of the table and we are not concerned with them any longer).
Part I for both PHYS460 and PHYS660 students
To build the BTW "sandpile" that self-organizes into a critical state, start from a lattice far from equilibrium where each site has a random value , and update the whole lattice according to the above rules until the sandpile self-organizes into a statistically stationary (time-independent) critical state with average .
a) Find explicitly the average height in the SOC state.
b) Add a single grain to this configuration at a randomly chosen site. This will trigger an avalanche of any size (up to the system size), which is why these state is critical. Record the size of the avalanche (number of sites that participate in the toppling rules) and time it takes for avalanche to end. If you collect many avalanches you should find power laws for the distribution function of avalanche sizes and times it take for the avalanche to cease: and . Plot these distributions on a log-log scale and from the slope of the line extract the critical exponents and . Although you can check all sites for stability, you can speed up the update by keeping two arrays which contain and coordinates of those sites that need to be checked for stability (each time a site topples, add its neighbors' positions to these arrays; then remove the last site from the arrays and check its stability; continue this process of removing sites and checking their stability, and adding neighbors to the array until there are no more sites to check).
c) Select few interesting avalanches and show graphically the shape of the avalanche (i.e., the sites which participate in the avalanche). Mark the sites which have changed their as well as those which after the avalanche is over have the same .
Part II for PHYS660 students only
References
- [1] P. Bak, C. Tang, and K. Wiesenfeld, Self-organized criticality, Phys. Rev. A 38, 364 (1988). [PDF]
- [2] S. K. Grumbacher, K. M. McEwen, D. A. Halverson, D. T. Jacobs, and J. Lindner, Self‐organized criticality: An experiment with sandpiles, Am. J. Phys. 61, 329 (1993). [PDF]
- [3] L. Laurson, M. J. Alava, and S. Zapperi, Power spectra of self-organized critical
sandpiles, J. Stat. Mech. L11001, (2005).