Complex quadratic polynomial

Since ${\displaystyle f_{c}(x)\,}$ is affine conjugate to the general form of the quadratic polynomial it is often used to study complex dynamics and to create images of Mandelbrot, Julia and Fatou sets.

When one wants change from ${\displaystyle \theta \,}$ to ${\displaystyle c\,}$:[5]

${\displaystyle c=c(\theta )={\frac {e^{2\pi \theta i}}{2}}\left(1-{\frac {e^{2\pi \theta i}}{2}}\right).}$

When one wants change from ${\displaystyle r\,}$ to ${\displaystyle c,\,}$ the parameter transformation is[6]

${\displaystyle c=c(r)\,=\,{\frac {1-(r-1)^{2}}{4}}=-{\frac {r}{2}}\left({\frac {r-2}{2}}\right)}$

and the transformation between the variables in ${\displaystyle z_{t+1}=z_{t}^{2}+c}$ and ${\displaystyle x_{t+1}=rx_{t}(1-x_{t})}$ is

${\displaystyle z=r\left({\frac {1}{2}}-x\right).}$

There is semi-conjugacy between the dyadic transformation (the doubling map) and the quadratic polynomial case of c = –2.

Here ${\displaystyle f^{n}\,}$ denotes the n-th iteration of the function ${\displaystyle f\,}$ (and not exponentiation of the function):

${\displaystyle f_{c}^{n}(z)=f_{c}^{1}(f_{c}^{n-1}(z))\,}$

so

${\displaystyle z_{n}=f_{c}^{n}(z_{0}).\,}$

Because of the possible confusion with exponentiation, some authors write ${\displaystyle f^{\circ n}\,}$ for the nth iterate of the function ${\displaystyle f.\,}$

The monic and centered form ${\displaystyle f_{c}(x)=x^{2}+c\,}$ can be marked by:

• the parameter ${\displaystyle c}$
• the external angle ${\displaystyle \theta }$ of the ray that lands:
  • at c in M on the parameter plane
  • at z = c in J(f) on the dynamic plane

so :

${\displaystyle f_{c}=f_{\theta }\,}$

${\displaystyle c=c({\theta })}$

The monic and centered form, sometimes called the Douady-Hubbard family of quadratic polynomials,[7] is typically used with variable ${\displaystyle z\,}$ and parameter ${\displaystyle c\,}$:

${\displaystyle f_{c}(z)=z^{2}+c.\,}$

When it is used as an evolution function of the discrete nonlinear dynamical system

${\displaystyle z_{n+1}=f_{c}(z_{n})\,}$

it is named the quadratic map:[8]

${\displaystyle f_{c}:z\to z^{2}+c.\,}$

The Mandelbrot set is the set of values of the parameter c for which the initial condition z₀ = 0 does not cause the iterates to diverge to infinity.

A critical point of ${\displaystyle f_{c}\,}$ is a point ${\displaystyle z_{cr}\,}$ in the dynamical plane such that the derivative vanishes:

${\displaystyle f_{c}'(z_{cr})=0.\,}$

Since

${\displaystyle f_{c}'(z)={\frac {d}{dz}}f_{c}(z)=2z}$

implies

${\displaystyle z_{cr}=0\,}$

we see that the only (finite) critical point of ${\displaystyle f_{c}\,}$ is the point ${\displaystyle z_{cr}=0\,}$.

${\displaystyle z_{0}}$ is an initial point for Mandelbrot set iteration.[9]

A critical value ${\displaystyle z_{cv}\ }$ of ${\displaystyle f_{c}\,}$ is the image of a critical point:

${\displaystyle z_{cv}=f_{c}(z_{cr})\,}$

Since

${\displaystyle z_{cr}=0\,}$

we have

${\displaystyle z_{cv}=c.\,}$

So the parameter ${\displaystyle c\,}$ is the critical value of ${\displaystyle f_{c}(z).\,}$

Dynamical plane with critical orbit falling into 3-period cycle

Dynamical plane with Julia set and critical orbit.

Dynamical plane : changes of critical orbit along internal ray of main cardioid for angle 1/6

Critical orbit tending to weakly attracting fixed point with abs(multiplier)=0.99993612384259

The forward orbit of a critical point is called a critical orbit. Critical orbits are very important because every attracting periodic orbit attracts a critical point, so studying the critical orbits helps us understand the dynamics in the Fatou set.[10][11][12]

${\displaystyle z_{0}=z_{cr}=0\,}$

${\displaystyle z_{1}=f_{c}(z_{0})=c\,}$

${\displaystyle z_{2}=f_{c}(z_{1})=c^{2}+c\,}$

${\displaystyle z_{3}=f_{c}(z_{2})=(c^{2}+c)^{2}+c\,}$

${\displaystyle ...\,}$

This orbit falls into an attracting periodic cycle if one exists.

Wikimedia Commons has media related to Critical orbits.

The critical sector is a sector of the dynamical plane containing the critical point.

${\displaystyle P_{n}(c)=f_{c}^{n}(z_{cr})=f_{c}^{n}(0)\,}$

so

${\displaystyle P_{0}(c)=0\,}$

${\displaystyle P_{1}(c)=c\,}$

${\displaystyle P_{2}(c)=c^{2}+c\,}$

${\displaystyle P_{3}(c)=(c^{2}+c)^{2}+c\,}$

These polynomials are used for:

• finding centers of these Mandelbrot set components of period n. Centers are roots of n-th critical polynomials

${\displaystyle centers=\{c:P_{n}(c)=0\}\,}$

• finding roots of Mandelbrot set components of period n (local minimum of ${\displaystyle P_{n}(c)\,}$)
• Misiurewicz points

${\displaystyle M_{n,k}=\{c:P_{k}(c)=P_{k+n}(c)\}\,}$

Critical curves

Diagrams of critical polynomials are called critical curves.[13]

These curves create the skeleton (the dark lines) of a bifurcation diagram.[14][15]

One can use the Julia-Mandelbrot 4-dimensional ( 4D) space for a global analysis of this dynamical system.[16]

w-plane and c-plane

In this space there are 2 basic types of 2-D planes:

• the dynamical (dynamic) plane, ${\displaystyle f_{c}\,}$-plane or c-plane
• the parameter plane or z-plane

There is also another plane used to analyze such dynamical systems w-plane:

• the conjugation plane[17]
• model plane[18]

Gamma parameter plane for complex logistic map ${\displaystyle z_{n+1}=\gamma z_{n}\left(1-z_{n}\right),}$

Multiplier map

The phase space of a quadratic map is called its parameter plane. Here:

${\displaystyle z_{0}=z_{cr}\,}$ is constant and ${\displaystyle c\,}$ is variable.

There is no dynamics here. It is only a set of parameter values. There are no orbits on the parameter plane.

The parameter plane consists of:

• The Mandelbrot set
  • The bifurcation locus = boundary of Mandelbrot set with
    • root points
  • Bounded hyperbolic components of the Mandelbrot set = interior of Mandelbrot set[19] with internal rays
• exterior of Mandelbrot set with
  • external rays
  • equipotential lines

There are many different subtypes of the parameter plane.[20][21]

See also :

• Boettcher map which maps exterior of mandelbrot set to the exterior of unit disc
• multiplier map which maps interior of hyperbolic component of Mandelbrot set to the interior of unit disc

 "The polynomial Pc maps each dynamical ray to another ray doubling the angle (which we measure in full turns, i.e. 0 = 1 = 2π rad = 360◦), and the dynamical rays of any polynomial “look like straight rays” near infinity. This allows us to study the Mandelbrot and Julia sets combinatorially, replacing the dynamical plane by the unit circle, rays by angles, and the quadratic polynomial by the doubling modulo one map." Virpi K a u k o[22]

On the dynamical plane one can find:

• The Julia set
• The Filled Julia set
• The Fatou set
• Orbits

The dynamical plane consists of:

• Fatou set
• Julia set

Here, ${\displaystyle c\,}$ is a constant and ${\displaystyle z\,}$ is a variable.

The two-dimensional dynamical plane can be treated as a Poincaré cross-section of three-dimensional space of continuous dynamical system.[23][24]

Dynamical z-planes can be divided in two groups :

• ${\displaystyle f_{0}}$ plane for ${\displaystyle c=0}$ ( see complex squaring map )
• ${\displaystyle f_{c}}$ planes ( all other planes for ${\displaystyle c\neq 0}$ )

The extended complex plane plus a point at infinity

• the Riemann sphere

On the parameter plane:

• ${\displaystyle c}$ is a variable
• ${\displaystyle z_{0}=0}$ is constant

The first derivative of ${\displaystyle f_{c}^{n}(z_{0})}$ with respect to c is

${\displaystyle z_{n}'={\frac {d}{dc}}f_{c}^{n}(z_{0}).}$

This derivative can be found by iteration starting with

${\displaystyle z_{0}'={\frac {d}{dc}}f_{c}^{0}(z_{0})=1}$

and then replacing at every consecutive step

${\displaystyle z_{n+1}'={\frac {d}{dc}}f_{c}^{n+1}(z_{0})=2\cdot {}f_{c}^{n}(z)\cdot {\frac {d}{dc}}f_{c}^{n}(z_{0})+1=2\cdot z_{n}\cdot z_{n}'+1.}$

This can easily be verified by using the chain rule for the derivative.

This derivative is used in the distance estimation method for drawing a Mandelbrot set.

On the dynamical plane:

• ${\displaystyle z}$ is a variable;
• ${\displaystyle c}$ is a constant.

At a fixed point ${\displaystyle z_{0}\,,}$

${\displaystyle f_{c}'(z_{0})={\frac {d}{dz}}f_{c}(z_{0})=2z_{0}.}$

At a periodic point z₀ of period p the first derivative of a function

${\displaystyle (f_{c}^{p})'(z_{0})={\frac {d}{dz}}f_{c}^{p}(z_{0})=\prod _{i=0}^{p-1}f_{c}'(z_{i})=2^{p}\prod _{i=0}^{p-1}z_{i}=\lambda }$

is often represented by ${\displaystyle \lambda }$ and referred to as the multiplier or the Lyapunov characteristic number. Its logarithm is known as the Lyapunov exponent. It used to check the stability of periodic (also fixed) points.

At a nonperiodic point, the derivative, denoted by ${\displaystyle z'_{n},\,}$ can be found by iteration starting with

${\displaystyle z'_{0}=1,\,}$

and then using

${\displaystyle z'_{n}=2*z_{n-1}*z'_{n-1}.\,}$

This derivative is used for computing the external distance to the Julia set.

The Schwarzian derivative (SD for short) of f is:[25]

${\displaystyle (Sf)(z)={\frac {f'''(z)}{f'(z)}}-{\frac {3}{2}}\left({\frac {f''(z)}{f'(z)}}\right)^{2}}$.