Mandelbrot Sets Generated from the Riemann Zeta Function



Updated: 2013/9/28
[List]
Published: 2013/7/21
This page lets you explore animated Mandelbrot-type sets generated from the Riemann zeta function.
Mandelbrot set generated from the zeta function
Mandelbrot Set Generated from the Zeta Function

The Classical Mandelbrot Set

The classical Mandelbrot set is the set of complex numbers c for which the sequence zn does not diverge to infinity under the recurrence relation

z0 = 0
zk+1 = zk2 + c

Starting from Costas Simeonidis's GPL-licensed Mandelbrot viewer, I modified the code and added animation. Points inside the set are usually shown in black, but here they are colored according to their complex values.

Full view of the Mandelbrot set
Full Mandelbrot View

(Iterations: 64)

Drag the mouse to draw a rectangle and zoom in on that region. The default iteration count, or recursion depth, is 16. Click [Start] to animate the process by gradually increasing the iteration count, so you can watch the set take shape step by step.

Zooming in near the origin gives the image below.

Mandelbrot set near the origin
Near the Origin

(Iterations: 64)

You can view a spiral structure here.

Spiral structure in the Mandelbrot set
Mandelbrot Spiral

(Iterations: 256)

The animation smoothly zooms from the full view down into this spiral.

An Exponential-Map Mandelbrot Set

Replacing the quadratic map with an exponential function also produces a Mandelbrot-like set:

F(z) = ez
zk+1 = F(zk) + c, z0 = 0

The result can be viewed here.

Overview of the exponential-map Mandelbrot set
Exponential-Map Overview

(Iterations: 4)

The default iteration count is 4. Press [Start] to animate the set as the iteration count increases.

A magnified view around the origin is shown below.

Exponential-map Mandelbrot set near the origin
Exponential Map near the Origin

(Iterations: 64)

You can also view the spiral region here, although rendering may take some time.

Spiral in the exponential-map Mandelbrot set
Exponential-Map Spiral

(Iterations: 256)

The Riemann Zeta Function

The Riemann zeta function is defined by the series

$$ \zeta(s)=\frac{1}{1^s}+\frac{1}{2^s}+\frac{1}{3^s}+\cdots $$

You can explore its graph here.

Graph of the Riemann zeta function
Graph of the Zeta Function

(Terms: 16)

Rendering is done with HTML Canvas and Web Workers. The number of terms is used to truncate the series, and the default is 16. Click [Start] to animate the graph as the number of terms increases and to watch the zeros appear.

A Zeta-Function Mandelbrot Set

Replacing F(z) with the zeta function produces another Mandelbrot-type set:

F(z) = ζ(z)
zk+1 = F(zk) + c, z0 = 0

The full view is available here.

Overview of the zeta-function Mandelbrot set
Zeta-Map Overview

(Iterations: 4, Terms: 16)

Use the [Start] buttons for iterations and terms to animate those parameters separately.

This set was introduced in 1998 by Wuon at Cambridge University (arXiv link).

A zoom near the first non-trivial zeros is shown below. Rendering may be slow.

Zeta-function Mandelbrot set near the zeta zeros
Near the Zeta Zeros

(Iterations: 64, Terms: 8)

You can view the spiral here. This view is especially slow to render.

Zeta-function Mandelbrot spiral
Zeta-Map Spiral

(Iterations: 256, Terms: 8)

A flower-like structure appears in this region.

Comparing the Mandelbrot Variants

Although the maps are different, these Mandelbrot-type sets display surprisingly similar patterns.

Comparison of Mandelbrot-Type Sets
F(z) z2 ez ζ(z)
Full View Full Mandelbrot set Full exponential-map Mandelbrot set Full zeta-function Mandelbrot set
Near
Origin/
Zeros
Mandelbrot set near the origin Exponential-map Mandelbrot set near the origin Zeta-function Mandelbrot set near the zeros
Spiral Mandelbrot spiral Exponential-map spiral Zeta-function spiral

The resemblance may be related to the following facts:

$$ e^z = \sum_{k=0}^{\infty} \frac{1}{k!}z^k $$ $$ \zeta(s)=\sum_{k=1}^{\infty}\frac{1}{k^s}=\sum_{k=1}^{\infty}e^{-s\log k} $$

Could There Be a Physical Meaning?

Each point in a Mandelbrot set either remains bounded or diverges to infinity. In a loose analogy, the universe may also either expand forever or eventually collapse. Perhaps each point could be imagined as representing a different universe.

How the Zeta Function Is Computed

Because the naïve series converges slowly, the program switches between several formulas as needed:

Accuracy tests for these formulas are summarized here:

Accuracy Comparison of the Algorithms

User Guide

Basic Controls

Animation Controls

Other Features


Related Articles:
Home >  Quantum Mechanics Articles


© 2013, xseek-qm.net


広告

量子革命: アインシュタインとボーア、偉大なる頭脳の激突 (マンジット・クマール, 青木薫 訳)(2017/1/28) 量子力学の奥深くに隠されているもの コペンハーゲン解釈から多世界理論へ (ショーン キャロル, 塩原通緒 訳)(2020/9/25) 量子力学の多世界解釈 なぜあなたは無数に存在するのか (和田純夫)(2022/12/15)