Douady rabbit

Last updated

A Douady rabbit is a fractal derived from the Julia set of the function , when parameter is near the center of one of the period three bulbs of the Mandelbrot set for a complex quadratic map. It is named after French mathematician Adrien Douady.

Contents

An example of a Douady rabbit. The colors show the number of iterations required to escape. Douady rabbit, plotted with Matplotlb.svg
An example of a Douady rabbit. The colors show the number of iterations required to escape.

Background

The Douady rabbit is generated by iterating the Mandelbrot set map on the complex plane, where parameter is fixed to lie in one of the two period three bulb off the main cardioid and ranging over the plane. The resulting image can be colored by corresponding each pixel with a starting value and calculating the amount of iterations required before the value of escapes a bounded region, after which it will diverge toward infinity.

It can also be described using the logistic map form of the complex quadratic map, specifically

This can be seen to be equivalent to the more common form via the affine transformation

which leads to

Setting produces the well-known iteration .

Irrespective of the specific iteration used, the filled Julia set associated with a given value of (or ) consists of all starting points (or ) for which the iteration remains bounded. Then, the Mandelbrot set consists of those values of (or ) for which the associated filled Julia set is connected. The Mandelbrot set can be viewed with respect to either or .

MandelbrotLambda.jpg
The Mandelbrot set in the plane
MandelbrotMuDouadyRabbit.jpg
The Mandelbrot set in the plane

Noting that is invariant under the substitution , the Mandelbrot set with respect to has additional horizontal symmetry. Since and are affine transformations of one another, or more specifically a similarity transformation, consisting of only scaling, rotation and translation, the filled Julia sets look similar for either form of the iteration given above.

Detailed description

Douady rabbit in an exponential family DouadyRabbitInExponentialFamily.jpg
Douady rabbit in an exponential family
Lamination of the rabbit Julia set Rabbit-lamination.png
Lamination of the rabbit Julia set
Representation of the dynamics inside the rabbit Douday rabbit rough dynamics.png
Representation of the dynamics inside the rabbit

We now continue to describe the Douady rabbit utilising the Mandelbrot set with respect to as shown in Figure 2. In this figure, the Mandelbrot set superficially appears as two back-to-back unit disks with sprouts or buds, such as the sprouts at the one- and five-o'clock positions on the right disk or the sprouts at the seven- and eleven-o'clock positions on the left disk. When is within one of these four sprouts, the associated filled Julia set in the mapping plane is said to be a Douady rabbit. For these values of , it can be shown that has and one other point as unstable (repelling) fixed points, and as an attracting fixed point. Moreover, the map has three attracting fixed points. A Douady rabbit consists of the three attracting fixed points , , and and their basins of attraction.

For example, Figure 4 shows the Douady rabbit in the plane when , a point in the five-o'clock sprout of the right disk. For this value of , the map has the repelling fixed points and . The three attracting fixed points of (also called period-three fixed points) have the locations

The red, green, and yellow points lie in the basins , , and of , respectively. The white points lie in the basin of .

The action of on these fixed points is given by the relations , , and .

Corresponding to these relations there are the results

Figure 4: Douady rabbit for
g
=
2.55268
-
0.959456
i
{\displaystyle \gamma =2.55268-0.959456i}
or
m
=
0.122565
-
0.744864
i
{\displaystyle \mu =0.122565-0.744864i} ColorDouadyRabbit1.jpg
Figure 4: Douady rabbit for or

As a second example, Figure 5 shows a Douady rabbit when , a point in the eleven-o'clock sprout on the left disk ( is invariant under this transformation). This rabbit is more symmetrical in the plane. The period-three fixed points then are located at

The repelling fixed points of itself are located at and . The three major lobes on the left, which contain the period-three fixed points ,, and , meet at the fixed point , and their counterparts on the right meet at the point . It can be shown that the effect of on points near the origin consists of a counterclockwise rotation about the origin of , or very nearly , followed by scaling (dilation) by a factor of .

Figure 5: Douady rabbit for
g
=
-
0.55268
+
0.959456
i
{\displaystyle \gamma =-0.55268+0.959456i}
or
m
=
0.122565
-
0.744864
i
{\displaystyle \mu =0.122565-0.744864i} ColorDouadyRabbit2.jpg
Figure 5: Douady rabbit for or

Variants

A twisted rabbit [1] is the composition of a rabbit polynomial with powers of Dehn twists about its ears. [2]

The corabbit is the symmetrical image of the rabbit. Here parameter . It is one of 2 other polynomials inducing the same permutation of their post-critical set are the rabbit.

3D

The Julia set has no direct analog in three dimensions.

4D

A quaternion Julia set with parameters and a cross-section in the plane. The Douady rabbit is visible in the cross-section.

Quaternion Julia Douady rabbit.jpg

Embedded

Parabolic julia set c = -1.125 + 0.21650635094611*i.png

A small embedded homeomorphic copy of rabbit in the center of a Julia set [3]

Fat

The fat rabbit or chubby rabbit has c at the root of the 1/3-limb of the Mandelbrot set. It has a parabolic fixed point with 3 petals. [4]

n-th eared

In general, the rabbit for the th bulb of the main cardioid will have ears [5] For example, a period four bulb rabbit has three ears.

Perturbed

Perturbed rabbit [6]

Twisted rabbit problem

In the early 1980s, Hubbard posed the so-called twisted rabbit problem, a polynomial classification problem. The goal is to determine Thurston equivalence types[ definition needed ] of functions of complex numbers that usually are not given by a formula (these are called topological polynomials): [7]

The problem was originally solved by Laurent Bartholdi and Volodymyr Nekrashevych [8] using iterated monodromic groups. The generalization of the problem to the case where the number of post-critical points is arbitrarily large has been solved as well. [9]

See also

Related Research Articles

<span class="mw-page-title-main">Electroweak interaction</span> Unified description of electromagnetism and the weak interaction

In particle physics, the electroweak interaction or electroweak force is the unified description of two of the four known fundamental interactions of nature: electromagnetism (electromagnetic interaction) and the weak interaction. Although these two forces appear very different at everyday low energies, the theory models them as two different aspects of the same force. Above the unification energy, on the order of 246 GeV, they would merge into a single force. Thus, if the temperature is high enough – approximately 1015 K – then the electromagnetic force and weak force merge into a combined electroweak force. During the quark epoch (shortly after the Big Bang), the electroweak force split into the electromagnetic and weak force. It is thought that the required temperature of 1015 K has not been seen widely throughout the universe since before the quark epoch, and currently the highest human-made temperature in thermal equilibrium is around 5.5x1012 K (from the Large Hadron Collider).

<span class="mw-page-title-main">Lorentz transformation</span> Family of linear transformations

In physics, the Lorentz transformations are a six-parameter family of linear transformations from a coordinate frame in spacetime to another frame that moves at a constant velocity relative to the former. The respective inverse transformation is then parameterized by the negative of this velocity. The transformations are named after the Dutch physicist Hendrik Lorentz.

<span class="mw-page-title-main">Pauli matrices</span> Matrices important in quantum mechanics and the study of spin

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices that are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

<span class="mw-page-title-main">Stress–energy tensor</span> Tensor describing energy momentum density in spacetime

The stress–energy tensor, sometimes called the stress–energy–momentum tensor or the energy–momentum tensor, is a tensor physical quantity that describes the density and flux of energy and momentum in spacetime, generalizing the stress tensor of Newtonian physics. It is an attribute of matter, radiation, and non-gravitational force fields. This density and flux of energy and momentum are the sources of the gravitational field in the Einstein field equations of general relativity, just as mass density is the source of such a field in Newtonian gravity.

In quantum field theory, Wilson loops are gauge invariant operators arising from the parallel transport of gauge variables around closed loops. They encode all gauge information of the theory, allowing for the construction of loop representations which fully describe gauge theories in terms of these loops. In pure gauge theory they play the role of order operators for confinement, where they satisfy what is known as the area law. Originally formulated by Kenneth G. Wilson in 1974, they were used to construct links and plaquettes which are the fundamental parameters in lattice gauge theory. Wilson loops fall into the broader class of loop operators, with some other notable examples being 't Hooft loops, which are magnetic duals to Wilson loops, and Polyakov loops, which are the thermal version of Wilson loops.

Variational Bayesian methods are a family of techniques for approximating intractable integrals arising in Bayesian inference and machine learning. They are typically used in complex statistical models consisting of observed variables as well as unknown parameters and latent variables, with various sorts of relationships among the three types of random variables, as might be described by a graphical model. As typical in Bayesian inference, the parameters and latent variables are grouped together as "unobserved variables". Variational Bayesian methods are primarily used for two purposes:

  1. To provide an analytical approximation to the posterior probability of the unobserved variables, in order to do statistical inference over these variables.
  2. To derive a lower bound for the marginal likelihood of the observed data. This is typically used for performing model selection, the general idea being that a higher marginal likelihood for a given model indicates a better fit of the data by that model and hence a greater probability that the model in question was the one that generated the data.

In rotordynamics, the rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane.

<span class="mw-page-title-main">Electromagnetic tensor</span> Mathematical object that describes the electromagnetic field in spacetime

In electromagnetism, the electromagnetic tensor or electromagnetic field tensor is a mathematical object that describes the electromagnetic field in spacetime. The field tensor was first used after the four-dimensional tensor formulation of special relativity was introduced by Hermann Minkowski. The tensor allows related physical laws to be written very concisely, and allows for the quantization of the electromagnetic field by Lagrangian formulation described below.

<span class="mw-page-title-main">Mathematical formulation of the Standard Model</span> Mathematics of a particle physics model

This article describes the mathematics of the Standard Model of particle physics, a gauge quantum field theory containing the internal symmetries of the unitary product group SU(3) × SU(2) × U(1). The theory is commonly viewed as describing the fundamental set of particles – the leptons, quarks, gauge bosons and the Higgs boson.

In general relativity, the Gibbons–Hawking–York boundary term is a term that needs to be added to the Einstein–Hilbert action when the underlying spacetime manifold has a boundary.

In mathematical physics, the gamma matrices, also called the Dirac matrices, are a set of conventional matrices with specific anticommutation relations that ensure they generate a matrix representation of the Clifford algebra It is also possible to define higher-dimensional gamma matrices. When interpreted as the matrices of the action of a set of orthogonal basis vectors for contravariant vectors in Minkowski space, the column vectors on which the matrices act become a space of spinors, on which the Clifford algebra of spacetime acts. This in turn makes it possible to represent infinitesimal spatial rotations and Lorentz boosts. Spinors facilitate spacetime computations in general, and in particular are fundamental to the Dirac equation for relativistic spin particles. Gamma matrices were introduced by Dirac in 1928.

The circles of Apollonius are any of several sets of circles associated with Apollonius of Perga, a renowned Greek geometer. Most of these circles are found in planar Euclidean geometry, but analogs have been defined on other surfaces; for example, counterparts on the surface of a sphere can be defined through stereographic projection.

In mathematics, more precisely in measure theory, an atom is a measurable set which has positive measure and contains no set of smaller positive measure. A measure which has no atoms is called non-atomic or atomless.

<span class="mw-page-title-main">Møller scattering</span> Electron-electron scattering

Møller scattering is the name given to electron-electron scattering in quantum field theory, named after the Danish physicist Christian Møller. The electron interaction that is idealized in Møller scattering forms the theoretical basis of many familiar phenomena such as the repulsion of electrons in the helium atom. While formerly many particle colliders were designed specifically for electron-electron collisions, more recently electron-positron colliders have become more common. Nevertheless, Møller scattering remains a paradigmatic process within the theory of particle interactions.

An external ray is a curve that runs from infinity toward a Julia or Mandelbrot set. Although this curve is only rarely a half-line (ray) it is called a ray because it is an image of a ray.

In mathematics, the Harish-Chandra isomorphism, introduced by Harish-Chandra (1951), is an isomorphism of commutative rings constructed in the theory of Lie algebras. The isomorphism maps the center of the universal enveloping algebra of a reductive Lie algebra to the elements of the symmetric algebra of a Cartan subalgebra that are invariant under the Weyl group .

In 3-dimensional topology, a part of the mathematical field of geometric topology, the Casson invariant is an integer-valued invariant of oriented integral homology 3-spheres, introduced by Andrew Casson.

The filled-in Julia set of a polynomial is a Julia set and its interior, non-escaping set

A hydrogen-like atom (or hydrogenic atom) is any atom or ion with a single valence electron. These atoms are isoelectronic with hydrogen. Examples of hydrogen-like atoms include, but are not limited to, hydrogen itself, all alkali metals such as Rb and Cs, singly ionized alkaline earth metals such as Ca+ and Sr+ and other ions such as He+, Li2+, and Be3+ and isotopes of any of the above. A hydrogen-like atom includes a positively charged core consisting of the atomic nucleus and any core electrons as well as a single valence electron. Because helium is common in the universe, the spectroscopy of singly ionized helium is important in EUV astronomy, for example, of DO white dwarf stars.

In probability theory, the family of complex normal distributions, denoted or , characterizes complex random variables whose real and imaginary parts are jointly normal. The complex normal family has three parameters: location parameter μ, covariance matrix , and the relation matrix . The standard complex normal is the univariate distribution with , , and .

References

  1. "A Geometric Solution to the Twisted Rabbit Problem by Jim Belk, University of St Andrews" (PDF). Archived (PDF) from the original on 2022-11-01. Retrieved 2022-05-03.
  2. Laurent Bartholdi; Volodymyr Nekrashevych (2006). "Thurston equivalence of topological polynomials". Acta Mathematica. 197: 1–51. arXiv: math/0510082 . doi:10.1007/s11511-006-0007-3.
  3. "Period-n Rabbit Renormalization. 'Rabbit's show' by Evgeny Demidov". Archived from the original on 2022-05-03. Retrieved 2022-05-03.
  4. Note on dynamically stable perturbations of parabolics by Tomoki Kawahira Archived October 2, 2006, at the Wayback Machine
  5. "Twisted Three-Eared Rabbits: Identifying Topological Quadratics Up To Thurston Equivalence by Adam Chodof" (PDF). Archived (PDF) from the original on 2022-05-03. Retrieved 2022-05-03.
  6. "Recent Research Papers (Only since 1999) Robert L. Devaney: Rabbits, Basilicas, and Other Julia Sets Wrapped in Sierpinski Carpets". Archived from the original on 2019-10-23. Retrieved 2020-04-07.
  7. "Polynomials, dynamics, and trees by Becca Winarski" (PDF). Archived (PDF) from the original on 2022-11-01. Retrieved 2022-05-08.
  8. Laurent Bartholdi; Volodymyr Nekrashevych (2005). "Thurston equivalence of topological polynomials". arXiv: math/0510082v3 .
  9. James Belk; Justin Lanier; Dan Margalit; Rebecca R. Winarski (2019). "Recognizing Topological Polynomials by Lifting Trees". arXiv: 1906.07680v1 [math.DS].

This article incorporates material from Douady Rabbit on PlanetMath, which is licensed under the Creative Commons Attribution/Share-Alike License.