68
$\begingroup$

The complex power series $$\sum_{n=1}^{\infty}\frac{z^{n^2}}{n^2}$$ has radius $1$ (Ratio Test) and is absolutely convergent along $|z|=1$. Recalling something that my calculus professor (Ray Mayer, emeritus of Reed College) showed me 15 years ago, I started looking at a "graph" of this function. More precisely, here are plots of the images of $z$ with constant magnitude under this power series.

graph of the fractal

(The mapped curves are the images of $|z|=1$, $|z|=0.9$, and $|z|=0.8$. At the far right you can see $\sum\limits_{n=1}^{\infty}\frac1{n^2}=\frac{\pi^2}{6}$.)

So... what the heck is going on with the fractal behavior of the image of the boundary? Is there any understanding of this kind of behavior from power series? For instance, rotating $z$ by some angles might leave you with roughly the original series after a bit of rotation, scaling and translation. But I haven't been able to see how that would all come together.

I have a hunch that the "berries" along the inside of the leaf happen around values of $z$ with interesting arguments, but I haven't sat down to map out what those arguments might be.


EDIT: Indeed, the "leaftips" and "berries" seem to happen at regular $z$ values. Starting at the largest leaftip in quadrant I and moving clockwise, the leaftips are the images of $\exp\left(\frac{\pi}{2}i\right), \exp\left(\frac{\pi}{4}i\right), \exp\left(\frac{\pi}{6}i\right)$,... Similarly, starting from the largest berry and moving clockwise the "seeds" of the berries are the images of $\exp\left(\frac{\pi}{3}i\right), \exp\left(\frac{\pi}{5}i\right), \exp\left(\frac{\pi}{7}i\right)$,...

It appears that the "arc" from the large berry in quadrant IV to the large berry in quadrant I is parametrized by $\exp(it)$ with $t\in\left(-\frac{\pi}{3},\frac{\pi}{3}\right)$. Further, that the large leaf in quadrant I that crosses into quadrant II is parametrized by $t\in\left(\frac{3\pi}{7},\frac{3\pi}{5}\right)$. These two sections (and indeed any of the "leaves", "subleaves", etc.) ought to be similar in the fractal sense.

  • 4
    I should note that this is just one example. This behavior is showing up for similar graphs of $\sum_{n=n_0}^{\infty}\frac{z^{f(n)}}{f(n)}$, where $f$ is other increasing functions from $\mathbb{N}$ to $\mathbb{N}$. So far, this one and the one with $f(n)=\binom{2n}{n}$ are my favorites.2011-11-23
  • 0
    Very intriguing! A random observation: your series is an antiderivative of $\sum_n z^{f(n)-1}$, an (almost) [lacunary](http://en.wikipedia.org/wiki/Lacunary_function) series. This sort of suggests trying to replace the "nice" $f(n)$ by $f(n)+1$. Also: what happens if you integrate once more?2011-11-23
  • 0
    How do you compute the series?2011-11-23
  • 0
    @lhf By approximation - just summing over the first $N$ terms. Although I think that smaller choices for N and m would still give smooth pictures in shorter time, here are the Maple commands that I used: `with(plots); m := 2000; N := 200; L := seq(evalf(sum(exp(I*n^2*t)/n^2, n = 1 .. N)), t = 0 .. 3.1415926, 3.1415926/m): display(complexplot([L]), complexplot(map(conjugate, [L])));`2011-11-23
  • 0
    @Henning That's interesting - I wonder if experts on lacunary series would have anything to say about this then. I tried plotting antiderivatives in a few cases, and they don't exhibit the same fractal behavior. In hindsight that makes sense, since these functions are locally close enough to constant that their antiderivatives would be locally close to linear.2011-11-23
  • 0
    @Henning: And summing $\frac{z^{n^2+1}}{n^2+1}$ still produces fractal behavior. If the graph above has a maple leaf shape and wraps once around the origin, I'd say the graph for $f(n)=n^2+1$ has more of an arrowhead shape and it wraps around twice.2011-11-23
  • 0
    If you separate out real and imaginary parts, note that the components look awfully like those canonical "continuous-but-nowhere-differentiable" functions, like Riemann's or Weierstrass's...2011-11-23
  • 4
    ...and speaking of which, the imaginary part of your sum essentially already is [Weierstrass's function](http://mathworld.wolfram.com/WeierstrassFunction.html).2011-11-23
  • 0
    BTW, you might want to see [this paper](http://dx.doi.org/10.1090/S0002-9939-98-04387-1).2011-11-23
  • 0
    In case my comments have been too obtuse: the fractal behavior observed in your curves is already apparent in the Weierstrass nowhere-differentiable function; it stands to reason that your curves will inherit that behavior.2011-12-05
  • 0
    alex.jordan: I suggest you erase the first term of the series, i.e z, and then instead of the all the series use the 2 nd and 3 rd terms. Do as you did above. You will get pointwise sum o two small circles with different radii, but "pointwise" sum will be according to actions of z^2 and z^9 on the circle.2011-12-05
  • 0
    I meant z^4. If you add z and increase the number of terms used in the approximation of the series, you will get a bigger and better approximation to the above fractal.2011-12-05
  • 1
    @J.M. I agree that the known information about Weierstrauss's function will almost certainly help (if not fully) explain this image. In another week I might finally have time to read up on it and I might ask you to link that article as an answer. But I'm also interested in the similar fractal behavior that I see for any $\sum \frac{x^{f(n)}}{f(n)}$ where $f$ grows fast enough for the series to be absolutely convergent when $|x|=1$. $\sum \frac{x^{n!}}{n!}$, $\sum \frac{x^{\binom{2n}{n}}}{\binom{2n}{n}}$, and $\sum \frac{x^{2^n}}{2^n}$ are some examples.2011-12-05
  • 1
    Note that to get the Weierstrass function you actually want to use $f(n)=b^n$ (maybe with a factor of $\pi$). In your base example, with $f(n)=n^2$, you actually get Riemann's "nowhere" differentiable function $$R(\theta)=\sum_{n\ge 0} {\sin\big(n^2\theta\big)\over n^2}$$ This is actually differentiable at points of the form $\pi{2p+1\over 2q+1}$, where $p$ and $q$ are integers. I'm not sure what goes on for other $f(n)$ (e.g. $f(n)=n^{k}$ for $k>2$), though $f(n)=n!$ nearly gives Darboux's function. See this paper for references: http://epubl.ltu.se/1402-1617/2003/320/LTU-EX-03320-SE.pdf2011-12-08
  • 0
    If anyone would like to summarize what has been said in the comments, I'll gladly award the bounty to them.2011-12-09
  • 0
    Some putzing around with the search terms `"lacunary Fourier series" fractal` turned up a few... interesting results.2011-12-10
  • 1
    $f(n)=n$ is actually also... exciting.2011-12-10
  • 1
    It gets even more funny when you take different functions $f(n)$ and $g(n)$ for the powers and the factors resp., i.e. $\sum_{n=1}z^{f(n)}/g(n)$.2011-12-10
  • 0
    Take $f(n)=2.5\cdot 2^n+1$ and $g(n)=(1.62)^n$.2011-12-10
  • 0
    I wish to add the following: The derivative of $f(z)$ is the well-known theta function:(http://functions.wolfram.com/EllipticFunctions/EllipticTheta1/)2011-12-13
  • 0
    A related question is now [on MO](http://mathoverflow.net/questions/83311/)2011-12-13
  • 1
    A Google search gave me this... diagrams like the ones here http://www.uam.es/personal_pdi/ciencias/acordoba/documentos/articulos/DifferentiabilityandDimension.pdf , exact calculations here http://mathworld.wolfram.com/WeierstrassFunction.html .2011-12-13
  • 0
    This generates the graph:http://www.wolframalpha.com/input/?i=ParametricPlot%5B+%7B%5Csum%5E%7B500%7D_%7Bn%3D1%7D%5Cfrac%7Bcos%5Bn%5E%7B2%7Dx%5D%7D%7Bn%5E%7B2%7D%7D%2C%5Csum%5E%7B500%7D_%7Bn%3D1%7D%5Cfrac%7BSin%5Bn%5E%7B2%7Dx%5D%7D%7Bn%5E%7B2%7D%7D%7D%2C%7Bx%2C0%2C2%5Cpi%7D%5D2011-12-14
  • 0
    I'm not sure I'm adding anything here, but there is a result - the Hadamard Gap theorem: http://en.wikipedia.org/wiki/Ostrowski%E2%80%93Hadamard_gap_theorem - that if the powers in a power series grow sufficiently quickly, the series is lacunary. One would then expect some sort of highly nondifferentiable behavior on the boundary of the disc of convergence, because otherwise what's the obstruction to analytic continuation? So perhaps it's not so surprising that we get these crazy fractals for a wide class of power series where the powers grow fast?2011-12-22

2 Answers 2