Periodic points of complex quadratic mappings explained

This article describes periodic points of some complex quadratic maps. A map is a formula for computing a value of a variable based on its own previous value or values; a quadratic map is one that involves the previous value raised to the powers one and two; and a complex map is one in which the variable and the parameters are complex numbers. A periodic point of a map is a value of the variable that occurs repeatedly after intervals of a fixed length.

These periodic points play a role in the theories of Fatou and Julia sets.

Definitions

Let

f _{c(z)} = z ^{2+c}

be the complex quadric mapping, where

z and c are complex numbers.

Notationally,

f ^{(k)} _{c} (z) is the k -fold composition of f _{c} with itself (not to be confused with the k th derivative of f _{c} )—that is, the value after the k-th iteration of the function f _{c.} Thusf ^{(k)} _{c} (z) = f ^{(k -1 )} _{c(f} _{c} (z)). p are points z of the dynamical plane such that f ^{(p)} _{c} (z) = z,

where

p is the smallest positive integer for which the equation holds at that z.

We can introduce a new function:

F _{p(z,f)} = f ^{(p)} _{c} (z) - z,

so periodic points are zeros of function

F _{p(z,f)} : points z satisfyingF _{p(z,f)} = 0 , 2 ^{p.}

Number of periodic points

The degree of the polynomial

F _{p(z,f)} describing periodic points is d = 2 ^{p} so it has exactly d = 2 ^{p} complex roots (= periodic points), counted with multiplicity.

Stability of periodic points (orbit) - multiplier

The multiplier (or eigenvalue, derivative)

m(f ^{p,z} _{0 ) =λ} of a rational map f iterated p times at cyclic point z _{0} is defined as:m(f ^{p,z} _{0 )} = λ = \begin{cases} f ^{p} (z _{0),} & if z _{0} \ne inf ty \\ (1)/(f ^{p} (z _{0)}) , & if z _{0} = inf ty \end{cases}

where

f ^{p\prime} (z _{0)} is the first derivative of f ^{p} with respect to z at z _{0} .

Because the multiplier is the same at all periodic points on a given orbit, it is called a multiplier of the periodic orbit.

The multiplier is:

• a complex number;
• invariant under conjugation of any rational map at its fixed point;^[1]
• used to check stability of periodic (also fixed) points with stability index

abs( λ ).

A periodic point is^[2]

• attracting when

abs( λ ) < 1 ;

• • super-attracting when

abs( λ ) = 0 ;

• • attracting but not super-attracting when

0 < abs( λ ) < 1 ;

• indifferent when

abs( λ ) = 1 ;

• • rationally indifferent or parabolic if

λ is a root of unity;

• • irrationally indifferent if

abs( λ ) =1 but multiplier is not a root of unity;

• repelling when

abs( λ ) > 1 .

Periodic points

• that are attracting are always in the Fatou set;
• that are repelling are in the Julia set;
• that are indifferent fixed points may be in one or the other.^[3] A parabolic periodic point is in the Julia set.

Period-1 points (fixed points)

Finite fixed points

Let us begin by finding all finite points left unchanged by one application of

f . These are the points that satisfy f _{c(z)=z} . That is, we wish to solvez ^{2+c=z,}

which can be rewritten as

z ^{2-z+c=0.}

Since this is an ordinary quadratic equation in one unknown, we can apply the standard quadratic solution formula:

\alpha _{1} = (1- \sqrt{ 1-4 c)/() } and \alpha _{2} = (1+ \sqrt{ 1-4 c)/() }.So for c \in C \setminus \{ 1 / 4 \} we have two finite fixed points \alpha _{1} and \alpha _{2} .

Since

\alpha _{1} = (1)/(2) - m and \alpha _{2} = (1)/(2) + m where m = (\sqrt{ 1-4 c)/() },

we have

\alpha _{1} + \alpha _{2} = 1 .

Thus fixed points are symmetrical about

z = 1 / 2 .

Complex dynamics

Here different notation is commonly used:^[4]

\alpha _{c} = (1- \sqrt{ 1-4 c)/() } with multiplier λ _{\alpha _{c}} = 1- \sqrt{ 1-4 c}

and

\beta _{c} = (1+ \sqrt{ 1-4 c)/() } with multiplier λ _{\beta _{c}} = 1+ \sqrt{ 1-4 c}.

Again we have

\alpha _{c} + \beta _{c} = 1 .

Since the derivative with respect to z is

P _{c'(z)} = (d)/(dz) P _{c(z)} = 2 z ,

we have

P _{c'(\alpha} _{c)} + P _{c'(\beta} _{c)=} 2 \alpha _{c} + 2 \beta _{c} = 2 (\alpha _{c} + \beta _{c)} = 2 .

This implies that

P _{c} can have at most one attractive fixed point.

These points are distinguished by the facts that:

\beta _{c} is:

• • the landing point of the external ray for angle=0 for

c \in M \setminus \left\{ 1 / 4 \right\}

• • the most repelling fixed point of the Julia set
  • the one on the right (whenever fixed point are not symmetrical around the real axis), it is the extreme right point for connected Julia sets (except for cauliflower).^[5]

\alpha _{c} is:

• • the landing point of several rays
  • attracting when

c is in the main cardioid of the Mandelbrot set, in which case it is in the interior of a filled-in Julia set, and therefore belongs to the Fatou set (strictly to the basin of attraction of finite fixed point)

• • parabolic at the root point of the limb of the Mandelbrot set
  • repelling for other values of

c

Special cases

An important case of the quadratic mapping is

c =0 . In this case, we get \alpha _{1} = 0 and \alpha _{2=1} . In this case, 0 is a superattractive fixed point, and 1 belongs to the Julia set.

Only one fixed point

We have

\alpha _{1=\alpha} _{2} exactly when 1-4 c =0 . This equation has one solution, c =1 / 4 , in which case \alpha _{1=\alpha} _{2=1/2} . In fact c =1 / 4 is the largest positive, purely real value for which a finite attractor exists.

Infinite fixed point

C to the Riemann sphere (extended complex plane) \hat{C } by adding infinity:\hat{C } = \mathbb \cup \

and extend

f _{c} such that f _{c(infty)=infty.}

Then infinity is:

• superattracting
• a fixed point of

f _{c} :^[6] $$f\_c(\backslash infty)=\backslash infty=f^\_c(\backslash infty).$$

Period-2 cycles

Period-2 cycles are two distinct points

\beta _{1} and \beta _{2} such that f _{c(\beta} _{1)} = \beta _{2} and f _{c(\beta} _{2)} = \beta _{1} , and hencef _{c(f} _{c(\beta} _{n))} = \beta _{n}

for

n \in \{ 1 , 2 \} :f _{c(f} _{c(z))} = (z ^{2+c)} ^{2+c} = z ^{4} + 2 cz ^{2} + c ^{2} + c.

Equating this to z, we obtain

z ^{4} + 2 cz ^{2} - z + c ^{2} + c = 0 .

This equation is a polynomial of degree 4, and so has four (possibly non-distinct) solutions. However, we already know two of the solutions. They are

\alpha _{1} and \alpha _{2} , computed above, since if these points are left unchanged by one application of f , then clearly they will be unchanged by more than one application of f .

Our 4th-order polynomial can therefore be factored in 2 ways:

First method of factorization

(z - \alpha _{1)(z-\alpha} _{2)(z-\beta} _{1)(z-\beta} _{2)} = 0 .

This expands directly as

x ^{4} - Ax ^{3} + Bx ^{2} - Cx + D = 0 (note the alternating signs), whereD = \alpha _{1} \alpha _{2} \beta _{1} \beta _{2,} C = \alpha _{1} \alpha _{2} \beta _{1} + \alpha _{1} \alpha _{2} \beta _{2} + \alpha _{1} \beta _{1} \beta _{2} + \alpha _{2} \beta _{1} \beta _{2,} B = \alpha _{1} \alpha _{2} + \alpha _{1} \beta _{1} + \alpha _{1} \beta _{2} + \alpha _{2} \beta _{1} + \alpha _{2} \beta _{2} + \beta _{1} \beta _{2,} A = \alpha _{1} + \alpha _{2} + \beta _{1} + \beta _{2.}

We already have two solutions, and only need the other two. Hence the problem is equivalent to solving a quadratic polynomial. In particular, note that

\alpha _{1} + \alpha _{2} = (1- \sqrt{ 1-4 c)/() } + \frac = \frac = 1

and

\alpha _{1} \alpha _{2} = (( 1- \sqrt{ 1-4 c)/()( 1+ \sqrt{ 1-4 c})}{ 4 }) = (1 ^{2} - (\sqrt{ 1-4 c)/()) ^{2}{4}=} (1-1+4 c)/(4) = (4 c)/(4) = c.

Adding these to the above, we get

D = c \beta _{1} \beta _{2} and A = 1 + \beta _{1} + \beta _{2} . Matching these against the coefficients from expanding f , we getD = c \beta _{1} \beta _{2} = c ^{2} + c and A = 1 + \beta _{1} + \beta _{2} = 0 .

From this, we easily get

\beta _{1} \beta _{2} = c + 1 and \beta _{1} + \beta _{2} = -1 .

From here, we construct a quadratic equation with

A' = 1 , B = 1 , C = c +1 and apply the standard solution formula to get\beta _{1} = (-1- \sqrt{ -3-4 c)/() } and \beta _{2} = (-1+ \sqrt{ -3-4 c)/() }.

Closer examination shows that:

f _{c(\beta} _{1)} = \beta _{2} and f _{c(\beta} _{2)} = \beta _{1,}

meaning these two points are the two points on a single period-2 cycle.

Second method of factorization

We can factor the quartic by using polynomial long division to divide out the factors

(z - \alpha _{1)} and (z - \alpha _{2),} which account for the two fixed points \alpha _{1} and \alpha _{2} (whose values were given earlier and which still remain at the fixed point after two iterations):(z ^{2+c)} ^{2} + c - z = (z ^{2} + c - z)(z ^{2} + z + c +1 ).

The roots of the first factor are the two fixed points. They are repelling outside the main cardioid.

The second factor has the two roots

(-1 \pm \sqrt{ -3-4 c)/() }. \,

These two roots, which are the same as those found by the first method, form the period-2 orbit.^[7]

Special cases

Again, let us look at

c =0 . Then\beta _{1} = (-1- i\sqrt{ 3)/() } and \beta _{2} = (-1+ i\sqrt{ 3)/() },

both of which are complex numbers. We have

| \beta _{1} | = | \beta _{2} | = 1 . Thus, both these points are "hiding" in the Julia set.Another special case is c =-1 , which gives \beta _{1} = 0 and \beta _{2} = -1 . This gives the well-known superattractive cycle found in the largest period-2 lobe of the quadratic Mandelbrot set.

Cycles for period greater than 2

The degree of the equation

f ^{(n)} (z) = z is 2ⁿ; thus for example, to find the points on a 3-cycle we would need to solve an equation of degree 8. After factoring out the factors giving the two fixed points, we would have a sixth degree equation.

There is no general solution in radicals to polynomial equations of degree five or higher, so the points on a cycle of period greater than 2 must in general be computed using numerical methods. However, in the specific case of period 4 the cyclical points have lengthy expressions in radicals.^[8]

In the case c = –2, trigonometric solutions exist for the periodic points of all periods. The case

z _{n+1} = z ^{2-2} _{n} is equivalent to the logistic map case r = 4: x _{n+1} =4 x _{n(1-x} _{n).} Here the equivalence is given by z =2-4 x. One of the k-cycles of the logistic variable x (all of which cycles are repelling) is \sin ^{2 \left( (2 \pi)/(2 ^{k-1})} \right), \sin ^{2 \left( 2 ⋅ (2 \pi)/(2 ^{k-1})} \right), \sin ^{2\left(2} ^{2 ⋅ (2 \pi)/(2 ^{k-1})} \right), \sin ^{2\left(2} ^{3 ⋅ (2 \pi)/(2 ^{k-1})} \right), ... , \sin ^{2\left(2} ^{k-1} (2 \pi)/(2 ^{k-1}) \right).

Further reading

• Geometrical properties of polynomial roots
• Alan F. Beardon, Iteration of Rational Functions, Springer 1991,
• Michael F. Barnsley (Author), Stephen G. Demko (Editor), Chaotic Dynamics and Fractals (Notes and Reports in Mathematics in Science and Engineering Series) Academic Pr (April 1986),
• Wolf Jung : Homeomorphisms on Edges of the Mandelbrot Set. Ph.D. thesis of 2002
• The permutations of periodic points in quadratic polynominials by J Leahy

External links

• Algebraic solution of Mandelbrot orbital boundaries by Donald D. Cross
• Brown Method by Robert P. Munafo
• arXiv:hep-th/0501235v2 V.Dolotin, A.Morozov: Algebraic Geometry of Discrete Dynamics. The case of one variable.

Notes and References

1. Alan F. Beardon, Iteration of Rational Functions, Springer 1991,, p. 41
2. Alan F. Beardon, Iteration of Rational Functions, Springer 1991,, page 99
3. http://www.ijon.de/mathe/julia/some_julia_sets_1_en.html Some Julia sets by Michael Becker
4. http://www.math.nagoya-u.ac.jp/~kawahira/works/cauliflower.pdf On the regular leaf space of the cauliflower by Tomoki Kawahira Source: Kodai Math. J. Volume 26, Number 2 (2003), 167-178.
5. http://www.ibiblio.org/e-notes/MSet/Attractor.htm Periodic attractor by Evgeny Demidov
6. [Robert L. Devaney|R L Devaney]
7. http://www.ibiblio.org/e-notes/MSet/Attractor.htm Period 2 orbit by Evgeny Demidov
8. https://arxiv.org/abs/0802.2565 Gvozden Rukavina : Quadratic recurrence equations - exact explicit solution of period four fixed points functions in bifurcation diagram