Catastrophe theory

Last updated

In mathematics, catastrophe theory is a branch of bifurcation theory in the study of dynamical systems; it is also a particular special case of more general singularity theory in geometry.

Contents

Bifurcation theory studies and classifies phenomena characterized by sudden shifts in behavior arising from small changes in circumstances, analysing how the qualitative nature of equation solutions depends on the parameters that appear in the equation. This may lead to sudden and dramatic changes, for example the unpredictable timing and magnitude of a landslide.

Catastrophe theory originated with the work of the French mathematician René Thom in the 1960s, and became very popular due to the efforts of Christopher Zeeman in the 1970s. It considers the special case where the long-run stable equilibrium can be identified as the minimum of a smooth, well-defined potential function (Lyapunov function). Small changes in certain parameters of a nonlinear system can cause equilibria to appear or disappear, or to change from attracting to repelling and vice versa, leading to large and sudden changes of the behaviour of the system. However, examined in a larger parameter space, catastrophe theory reveals that such bifurcation points tend to occur as part of well-defined qualitative geometrical structures.

In the late 1970s, applications of catastrophe theory to areas outside its scope began to be criticized, especially in biology and social sciences. [1] [2] Zahler and Sussmann, in a 1977 article in Nature , referred to such applications as being "characterised by incorrect reasoning, far-fetched assumptions, erroneous consequences, and exaggerated claims". [3] As a result, catastrophe theory has become less popular in applications. [4]

Elementary catastrophes

Catastrophe theory analyzes degenerate critical points of the potential function — points where not just the first derivative, but one or more higher derivatives of the potential function are also zero. These are called the germs of the catastrophe geometries. The degeneracy of these critical points can be unfolded by expanding the potential function as a Taylor series in small perturbations of the parameters.

When the degenerate points are not merely accidental, but are structurally stable, the degenerate points exist as organising centres for particular geometric structures of lower degeneracy, with critical features in the parameter space around them. If the potential function depends on two or fewer active variables, and four or fewer active parameters, then there are only seven generic structures for these bifurcation geometries, with corresponding standard forms into which the Taylor series around the catastrophe germs can be transformed by diffeomorphism (a smooth transformation whose inverse is also smooth).[ citation needed ] These seven fundamental types are now presented, with the names that Thom gave them.

Potential functions of one active variable

Catastrophe theory studies dynamical systems that describe the evolution [5] of a state variable over time :

In the above equation, is referred to as the potential function, and is often a vector or a scalar which parameterise the potential function. The value of may change over time, and it can also be referred to as the control variable. In the following examples, parameters like are such controls.

Fold catastrophe

Fold catastrophe, with surface
z
=
-
y
2
-
0.1
x
4
{\displaystyle z=-y^{2}-0.1x^{4}}
. Fold catastrophe animation.gif
Fold catastrophe, with surface .
Stable and unstable pair of extrema disappear at a fold bifurcation Fold bifurcation.svg
Stable and unstable pair of extrema disappear at a fold bifurcation

When a< 0, the potential V has two extrema - one stable, and one unstable. If the parameter a is slowly increased, the system can follow the stable minimum point. But at a = 0 the stable and unstable extrema meet, and annihilate. This is the bifurcation point. At a> 0 there is no longer a stable solution. If a physical system is followed through a fold bifurcation, one therefore finds that as a reaches 0, the stability of the a< 0 solution is suddenly lost, and the system will make a sudden transition to a new, very different behaviour. This bifurcation value of the parameter a is sometimes called the "tipping point".

Cusp catastrophe

Diagram of cusp catastrophe, showing curves (brown, red) of x satisfying dV/dx = 0 for parameters (a,b), drawn for parameter b continuously varied, for several values of parameter a.

Outside the cusp locus of bifurcations (blue), for each point (a,b) in parameter space there is only one extremising value of x. Inside the cusp, there are two different values of x giving local minima of V(x) for each (a,b), separated by a value of x giving a local maximum. Cusp catastrophe.svg
Diagram of cusp catastrophe, showing curves (brown, red) of x satisfying dV/dx = 0 for parameters (a,b), drawn for parameter b continuously varied, for several values of parameter a.

Outside the cusp locus of bifurcations (blue), for each point (a,b) in parameter space there is only one extremising value of x. Inside the cusp, there are two different values of x giving local minima of V(x) for each (a,b), separated by a value of x giving a local maximum.
Cusp shape in parameter space (a,b) near the catastrophe point, showing the locus of fold bifurcations separating the region with two stable solutions from the region with one. Cusp shape.svg
Cusp shape in parameter space (a,b) near the catastrophe point, showing the locus of fold bifurcations separating the region with two stable solutions from the region with one.
Pitchfork bifurcation at a = 0 on the surface b = 0 Pitchfork bifurcation left.svg
Pitchfork bifurcation at a = 0 on the surface b = 0
Cusp catastrophe, with surface
z
=
-
0.1
(
x
4
+
y
4
)
-
y
3
+
x
y
{\displaystyle z=-0.1(x^{4}+y^{4})-y^{3}+xy}
. Cusp catastrophe animation gif.gif
Cusp catastrophe, with surface .

The cusp geometry is very common when one explores what happens to a fold bifurcation if a second parameter, b, is added to the control space. Varying the parameters, one finds that there is now a curve (blue) of points in (a,b) space where stability is lost, where the stable solution will suddenly jump to an alternate outcome.

But in a cusp geometry the bifurcation curve loops back on itself, giving a second branch where this alternate solution itself loses stability, and will make a jump back to the original solution set. By repeatedly increasing b and then decreasing it, one can therefore observe hysteresis loops, as the system alternately follows one solution, jumps to the other, follows the other back, and then jumps back to the first.

However, this is only possible in the region of parameter space a< 0. As a is increased, the hysteresis loops become smaller and smaller, until above a = 0 they disappear altogether (the cusp catastrophe), and there is only one stable solution.

One can also consider what happens if one holds b constant and varies a. In the symmetrical case b = 0, one observes a pitchfork bifurcation as a is reduced, with one stable solution suddenly splitting into two stable solutions and one unstable solution as the physical system passes to a< 0 through the cusp point (0,0) (an example of spontaneous symmetry breaking). Away from the cusp point, there is no sudden change in a physical solution being followed: when passing through the curve of fold bifurcations, all that happens is an alternate second solution becomes available.

A famous suggestion is that the cusp catastrophe can be used to model the behaviour of a stressed dog, which may respond by becoming cowed or becoming angry. [6] The suggestion is that at moderate stress (a> 0), the dog will exhibit a smooth transition of response from cowed to angry, depending on how it is provoked. But higher stress levels correspond to moving to the region (a< 0). Then, if the dog starts cowed, it will remain cowed as it is irritated more and more, until it reaches the 'fold' point, when it will suddenly, discontinuously snap through to angry mode. Once in 'angry' mode, it will remain angry, even if the direct irritation parameter is considerably reduced.

A simple mechanical system, the "Zeeman Catastrophe Machine", nicely illustrates a cusp catastrophe. In this device, smooth variations in the position of the end of a spring can cause sudden changes in the rotational position of an attached wheel. [7]

Catastrophic failure of a complex system with parallel redundancy can be evaluated based on the relationship between local and external stresses. The model of the structural fracture mechanics is similar to the cusp catastrophe behavior. The model predicts reserve ability of a complex system.

Other applications include the outer sphere electron transfer frequently encountered in chemical and biological systems, [8] modelling the dynamics of cloud condensation nuclei in the atmosphere, [9] and modelling real estate prices. [10]

Fold bifurcations and the cusp geometry are by far the most important practical consequences of catastrophe theory. They are patterns which reoccur again and again in physics, engineering and mathematical modelling. They produce the strong gravitational lensing events and provide astronomers with one of the methods used for detecting black holes and the dark matter of the universe, via the phenomenon of gravitational lensing producing multiple images of distant quasars. [11]

The remaining simple catastrophe geometries are very specialised in comparison, and presented here only for curiosity value.

Swallowtail catastrophe

Swallowtail catastrophe, with surface
z
=
-
y
4
-
0.1
x
4
+
(
1
-
x
2
)
y
2
+
0.4
x
y
{\displaystyle z=-y^{4}-0.1x^{4}+(1-x^{2})y^{2}+0.4xy} Swallowtail catastrophe animation gif.gif
Swallowtail catastrophe, with surface
Swallowtail catastrophe surface Smallow tail.jpg
Swallowtail catastrophe surface

The control parameter space is three-dimensional. The bifurcation set in parameter space is made up of three surfaces of fold bifurcations, which meet in two lines of cusp bifurcations, which in turn meet at a single swallowtail bifurcation point.

As the parameters go through the surface of fold bifurcations, one minimum and one maximum of the potential function disappear. At the cusp bifurcations, two minima and one maximum are replaced by one minimum; beyond them the fold bifurcations disappear. At the swallowtail point, two minima and two maxima all meet at a single value of x. For values of a> 0, beyond the swallowtail, there is either one maximum-minimum pair, or none at all, depending on the values of b and c. Two of the surfaces of fold bifurcations, and the two lines of cusp bifurcations where they meet for a< 0, therefore disappear at the swallowtail point, to be replaced with only a single surface of fold bifurcations remaining. Salvador Dalí's last painting, The Swallow's Tail , was based on this catastrophe.

Butterfly catastrophe

Butterfly catastrophe, with surface
z
=
-
20
x
5
+
4
x
3
-
2
x
y
-
0.5
(
x
4
+
y
4
)
{\displaystyle z=-20x^{5}+4x^{3}-2xy-0.5(x^{4}+y^{4})}
. Butterfly catastrophe animation gif.gif
Butterfly catastrophe, with surface .

Depending on the parameter values, the potential function may have three, two, or one different local minima, separated by the loci of fold bifurcations. At the butterfly point, the different 3-surfaces of fold bifurcations, the 2-surfaces of cusp bifurcations, and the lines of swallowtail bifurcations all meet up and disappear, leaving a single cusp structure remaining when a> 0.

Potential functions of two active variables

A surface with a hyperbolic umbilic and its focal surface. The hyperbolic umbilic catastrophe is just the upper part of this image. Hyperbolic umbilic focal surface.png
A surface with a hyperbolic umbilic and its focal surface. The hyperbolic umbilic catastrophe is just the upper part of this image.
A surface with an elliptical umbilic, and its focal surface. The elliptic umbilic catastrophe is just the upper part of this image. Elliptical umbilic focal surface.png
A surface with an elliptical umbilic, and its focal surface. The elliptic umbilic catastrophe is just the upper part of this image.

Umbilic catastrophes are examples of corank 2 catastrophes. They can be observed in optics in the focal surfaces created by light reflecting off a surface in three dimensions and are intimately connected with the geometry of nearly spherical surfaces: umbilical point. Thom proposed that the hyperbolic umbilic catastrophe modeled the breaking of a wave and the elliptical umbilic modeled the creation of hair-like structures.

Hyperbolic umbilic catastrophe

Elliptic umbilic catastrophe

Parabolic umbilic catastrophe

Arnold's notation

Vladimir Arnold gave the catastrophes the ADE classification, due to a deep connection with simple Lie groups. [ citation needed ]

There are objects in singularity theory which correspond to most of the other simple Lie groups.

Optics

Caustics underwater are catastrophes, and their singular points are the same even as the surface of the water is everchanging.
The edge of the rainbow is a fold catastrophe, and so it has a diffraction pattern described by the Airy function. This is generic and stable for any shape of the water droplet. Supernumerary-rainbows-jb.jpg
The edge of the rainbow is a fold catastrophe, and so it has a diffraction pattern described by the Airy function. This is generic and stable for any shape of the water droplet.

As predicted by catastrophe theory, singularities are generic, and stable under perturbation. This explains how the bright lines and surfaces are stable under perturbation. The caustics one sees at the bottom of a swimming pool, for example, have a distinctive texture and only has a few types of singular points, even though the surface of the water is ever changing. [12]

The edge of the rainbow, for example, has a fold catastrophe. Due to the wave nature of light, the catastrophe has fine diffraction details described by the Airy function. This is a generic result and does not depend on the precise shape of the water droplet, and so the edge of the rainbow always has the shape of an Airy function. [13] [14] The same Airy function fold catastrophe can be seen in nuclear-nuclear scattering ("nuclear rainbow"). [15]

The cusp catastrophe is the next-simplest to observe. Due to the wave nature of light, the catastrophe has fine diffraction details described by the Pearcey function. [16] Higher-order catastrophes, such as the swallowtail and the butterfly, have also been observed. [17]

Cusp catastrophe caustic, generated from a circle and parallel rays. Circle caustic.png
Cusp catastrophe caustic, generated from a circle and parallel rays.
A photograph of a cusp caustic produced by illuminating a flat surface with a laser beam through a droplet of water. This is a Pearcey function and stable under perturbations. A photograph of a cusp caustic.png
A photograph of a cusp caustic produced by illuminating a flat surface with a laser beam through a droplet of water. This is a Pearcey function and stable under perturbations.

See also

Related Research Articles

<span class="mw-page-title-main">Dynamical system</span> Mathematical model of the time dependence of a point in space

In mathematics, a dynamical system is a system in which a function describes the time dependence of a point in an ambient space, such as in a parametric curve. Examples include the mathematical models that describe the swinging of a clock pendulum, the flow of water in a pipe, the random motion of particles in the air, and the number of fish each springtime in a lake. The most general definition unifies several concepts in mathematics such as ordinary differential equations and ergodic theory by allowing different choices of the space and how time is measured. Time can be measured by integers, by real or complex numbers or can be a more general algebraic object, losing the memory of its physical origin, and the space may be a manifold or simply a set, without the need of a smooth space-time structure defined on it.

In mathematics, an equation is a mathematical formula that expresses the equality of two expressions, by connecting them with the equals sign =. The word equation and its cognates in other languages may have subtly different meanings; for example, in French an équation is defined as containing one or more variables, while in English, any well-formed formula consisting of two expressions related with an equals sign is an equation.

In mathematics, a coefficient is a multiplicative factor involved in some term of a polynomial, a series, or an expression. It may be a number (dimensionless), in which case it is known as a numerical factor. It may also be a constant with units of measurement, in which it is known as a constant multiplier. In general, coefficients may be any expression. When the combination of variables and constants is not necessarily involved in a product, it may be called a parameter.

<span class="mw-page-title-main">Quadratic formula</span> Formula that provides the solutions to a quadratic equation

In elementary algebra, the quadratic formula is a closed-form expression describing the solutions of a quadratic equation. Other ways of solving quadratic equations, such as completing the square, yield the same solutions.

<span class="mw-page-title-main">Algebraic curve</span> Curve defined as zeros of polynomials

In mathematics, an affine algebraic plane curve is the zero set of a polynomial in two variables. A projective algebraic plane curve is the zero set in a projective plane of a homogeneous polynomial in three variables. An affine algebraic plane curve can be completed in a projective algebraic plane curve by homogenizing its defining polynomial. Conversely, a projective algebraic plane curve of homogeneous equation h(x, y, t) = 0 can be restricted to the affine algebraic plane curve of equation h(x, y, 1) = 0. These two operations are each inverse to the other; therefore, the phrase algebraic plane curve is often used without specifying explicitly whether it is the affine or the projective case that is considered.

In mathematics, singularity theory studies spaces that are almost manifolds, but not quite. A string can serve as an example of a one-dimensional manifold, if one neglects its thickness. A singularity can be made by balling it up, dropping it on the floor, and flattening it. In some places the flat string will cross itself in an approximate "X" shape. The points on the floor where it does this are one kind of singularity, the double point: one bit of the floor corresponds to more than one bit of string. Perhaps the string will also touch itself without crossing, like an underlined "U". This is another kind of singularity. Unlike the double point, it is not stable, in the sense that a small push will lift the bottom of the "U" away from the "underline".

<span class="mw-page-title-main">Squeeze mapping</span> Linear mapping permuting rectangles of the same area

In linear algebra, a squeeze mapping, also called a squeeze transformation, is a type of linear map that preserves Euclidean area of regions in the Cartesian plane, but is not a rotation or shear mapping.

<i>The Swallows Tail</i> 1983 painting by Salvador Dalí

The Swallow's Tail — Series of Catastrophes was Salvador Dalí's last painting. It was completed in May 1983, as the final part of a series based on the mathematical catastrophe theory of René Thom.

<span class="mw-page-title-main">Bifurcation theory</span> Study of sudden qualitative behavior changes caused by small parameter changes

Bifurcation theory is the mathematical study of changes in the qualitative or topological structure of a given family of curves, such as the integral curves of a family of vector fields, and the solutions of a family of differential equations. Most commonly applied to the mathematical study of dynamical systems, a bifurcation occurs when a small smooth change made to the parameter values of a system causes a sudden 'qualitative' or topological change in its behavior. Bifurcations occur in both continuous systems and discrete systems.

In the mathematical area of bifurcation theory a saddle-node bifurcation, tangential bifurcation or fold bifurcation is a local bifurcation in which two fixed points of a dynamical system collide and annihilate each other. The term 'saddle-node bifurcation' is most often used in reference to continuous dynamical systems. In discrete dynamical systems, the same bifurcation is often instead called a fold bifurcation. Another name is blue sky bifurcation in reference to the sudden creation of two fixed points.

In mathematics, an orientation of a curve is the choice of one of the two possible directions for travelling on the curve. For example, for Cartesian coordinates, the x-axis is traditionally oriented toward the right, and the y-axis is upward oriented.

In mathematics, a variable is a symbol that represents a mathematical object. A variable may represent a number, a vector, a matrix, a function, the argument of a function, a set, or an element of a set.

In the mathematics of evolving systems, the concept of a center manifold was originally developed to determine stability of degenerate equilibria. Subsequently, the concept of center manifolds was realised to be fundamental to mathematical modelling.

<span class="mw-page-title-main">Umbilic torus</span>

The umbilic torus or umbilic bracelet is a single-edged 3-dimensional shape. The lone edge goes three times around the ring before returning to the starting point. The shape also has a single external face. A cross section of the surface forms a deltoid.

Numerical continuation is a method of computing approximate solutions of a system of parameterized nonlinear equations,

<span class="mw-page-title-main">Umbilical point</span>

In the differential geometry of surfaces in three dimensions, umbilics or umbilical points are points on a surface that are locally spherical. At such points the normal curvatures in all directions are equal, hence, both principal curvatures are equal, and every tangent vector is a principal direction. The name "umbilic" comes from the Latin umbilicus (navel).

In numerical linear algebra, the alternating-direction implicit (ADI) method is an iterative method used to solve Sylvester matrix equations. It is a popular method for solving the large matrix equations that arise in systems theory and control, and can be formulated to construct solutions in a memory-efficient, factored form. It is also used to numerically solve parabolic and elliptic partial differential equations, and is a classic method used for modeling heat conduction and solving the diffusion equation in two or more dimensions. It is an example of an operator splitting method.

In mathematics, and especially affine differential geometry, the affine focal set of a smooth submanifold M embedded in a smooth manifold N is the caustic generated by the affine normal lines. It can be realised as the bifurcation set of a certain family of functions. The bifurcation set is the set of parameter values of the family which yield functions with degenerate singularities. This is not the same as the bifurcation diagram in dynamical systems.

In mathematics, the slow manifold of an equilibrium point of a dynamical system occurs as the most common example of a center manifold. One of the main methods of simplifying dynamical systems, is to reduce the dimension of the system to that of the slow manifold—center manifold theory rigorously justifies the modelling. For example, some global and regional models of the atmosphere or oceans resolve the so-called quasi-geostrophic flow dynamics on the slow manifold of the atmosphere/oceanic dynamics, and is thus crucial to forecasting with a climate model.

In probability theory, a Cauchy process is a type of stochastic process. There are symmetric and asymmetric forms of the Cauchy process. The unspecified term "Cauchy process" is often used to refer to the symmetric Cauchy process.

References

  1. Murray, Stacey R. "The Rise and Fall of Catastrophe Theory". Encyclopedia.com . Retrieved 2 November 2021.
  2. Horgan, John (2015). The End of Science: Facing the Limits of Knowledge in the Twilight of the Scientific Age. New York: Basic Books. p. 213. ISBN   978-0-465-05085-7.
  3. Zahler, Raphael S.; Sussmann, Hector J. (1977). "Claims and accomplishments of applied catastrophe theory". Nature. 269 (5631): 759–763. Bibcode:1977Natur.269..759Z. doi:10.1038/269759a0. ISSN   1476-4687. S2CID   4205198 . Retrieved 2021-11-02.
  4. Rosser, J. Barkley (October 2007). "The rise and fall of catastrophe theory applications in economics: Was the baby thrown out with the bathwater?". Journal of Economic Dynamics and Control. 31 (10): 3255–3280. doi:10.1016/j.jedc.2006.09.013.
  5. Wagenmakers, E. J.; van der Maas, H. L. J.; Molenaar, P. C. M. (2005). "Fitting the cusp catastrophe model". Encyclopedia of Statistics in Behavioral Science.
  6. E.C. Zeeman, Catastrophe Theory, Scientific American , April 1976; pp. 65–70, 75–83
  7. Cross, Daniel J., Interactive rendering of Zeeman's Catastrophe Machine
  8. Xu, F (1990). "Application of catastrophe theory to the ∆G to -∆G relationship in electron transfer reactions". Zeitschrift für Physikalische Chemie. Neue Folge. 166: 79–91. doi:10.1524/zpch.1990.166.Part_1.079. S2CID   101078817.
  9. Arabas, S; Shima, S. (2017). "On the CCN (de)activation nonlinearities". Nonlinear Processes in Geophysics. 24 (3): 535–542. arXiv: 1608.08187 . Bibcode:2017NPGeo..24..535A. doi: 10.5194/npg-24-535-2017 . S2CID   24669360.
  10. Bełej, Mirosław; Kulesza, Sławomir (2012). "Modeling the Real Estate Prices in Olsztyn under Instability Conditions". Folia Oeconomica Stetinensia. 11 (1): 61–72. doi: 10.2478/v10031-012-0008-7 .
  11. A.O. Petters, H. Levine and J. Wambsganss, Singularity Theory and Gravitational Lensing", Birkhäuser Boston (2001)
  12. Berry, M. V.; Upstill, C. (1980-01-01), Wolf, E. (ed.), "IV Catastrophe Optics: Morphologies of Caustics and Their Diffraction Patterns", Progress in Optics, vol. 18, Elsevier, pp. 257–346, retrieved 2024-04-25
  13. Adam, John A. (2002-01-01). "The mathematical physics of rainbows and glories". Physics Reports. 356 (4): 229–365. doi:10.1016/S0370-1573(01)00076-X. ISSN   0370-1573.
  14. "AMS :: Feature Column :: The Mathematics of the Rainbow, Part II". American Mathematical Society. Retrieved 2024-04-25.
  15. Khoa, Dao T; Oertzen, W von; Bohlen, H G; Ohkubo, S (2007-03-01). "Nuclear rainbow scattering and nucleus–nucleus potential". Journal of Physics G: Nuclear and Particle Physics. 34 (3): R111–R164. doi:10.1088/0954-3899/34/3/R01. ISSN   0954-3899.
  16. MacBeath, Darryl (November 2016). The Pearcey function and the cusp catastrophe (Thesis thesis).
  17. Zannotti, Alessandro; Diebel, Falko; Boguslawski, Martin; Denz, Cornelia (May 2017). "Optical catastrophes of the swallowtail and butterfly beams". New Journal of Physics. 19 (5): 053004. doi:10.1088/1367-2630/aa6ecd. ISSN   1367-2630.

Bibliography