What is the Riemann surface of $y=\sqrt{z+z^2+z^4+\cdots +z^{2^n}+\cdots}$? The following appears as the second-to-last problem of Stewart's Complex Analysis:

Describe the Riemann surface of the function $y=\sqrt{z+z^2+z^4+\cdots +z^{2^n}+\cdots}$.

This problem intimidated me when I first saw it as an undergrad, as the series under the root isn't the expansion of any function I know (then or now). That renders the Riemann surface more subtle than usual, since any info about poles and zeroes must be found from this series alone. For instance, it has a zero at $z=0$ and therefore a branch point. But there should be other zeroes — indeed, infinitely many zeroes— since the number associated with any truncation of the series grows as $2^N$!. This suggests that the Riemann surface of this function must be exotic.
I'll give a sketch of an argument below, taking the pedestrian approach of hunting poles and zeroes of the function. However, I'd be glad to see an elegant and advanced perspective of the problem, particularly if it uses some ideas which I wouldn't have encountered as an undergraduate. The more food for thought, the better!
EDIT: As noted by O.L.'s answer below, the series above is the canonical example of a lacunary function with the unit circle as natural boundary. What that leaves open is the nature of the resulting Riemann surface, which would seem to hinge upon the quantity of the function's zeros within the unit circle. David Speyer has given an argument in his answer that this number is at least 10. Is there a more definitive result on the zeroes of this lacuuary function?
 A: This answer is closer to a sketch than a proof. If you notice gaps/mistakes/nonsense in this presentation, please let me know...
We want to examine the behavior of the function $f(z)$ inside the square root. It evidently vanishes linearly at the origin. Indeed, for all $|z|<1$ this series is convergent by comparison with the geometric series. If $|z|>1$, however, the series diverges since the terms do not tend to zero. Hence the unit disk is the disk of convergence.
Compared to this straightforward picture, the behavior for $|z|=1$ is more exotic. Consider first a rational rotation, i.e. $z=e^{i\phi}$ for some $\phi\in 2\pi\mathbb{Q}$. Then $2^n\phi$ is either $0$ or $\phi$ modulo $2\pi$ for some sufficiently large integer $N$. (This is equivalent to the binary representation of any rational number eventually terminating or becoming periodic.) In the former case, it follows that $z^{2^n}=1$ for $n\geq N$. For the latter, $z^{2^n}$ is a certain periodic sequence of roots of unity; since this sequence does not include 1, the sum of each period is nonzero and identical. In both cases we may conclude that the series diverges. So every rational rotation gives a singularity.
On the other hand, if $\phi/2\pi$ is irrational, then the sequence $2^n\phi$ modulo $2\pi$ is neither periodic nor terminating: therefore the orbit of $2^n\phi$ will be dense in $[0,2\pi]$. As a consequence, the series sums to zero. So every irrational rotation corresponds to a zero.
Therefore it would appear that the unit circle is interwoven with (countably many) singularities and (uncountably many) zeroes. What does this entail for $y=\sqrt{f(z)}?$ Here my understanding goes from loose to shaky. In the standard case, each zero of $f(z)$ would be a branch point of the corresponding Riemann surface. Consequently one would seem to have a Riemann surface of (uncountably) infinite genus.
Under the (generous) assumption that this is all sound, there are questions which linger for me:


*

*What is the nature of the singularities at commensurate phases? (Are they poles of some multiplicity, or are they essential singularities?)

*While there are uncountably many zeroes, there is a symmetry under complex conjugation for all roots save $z=0$. Consequently the 'number' of such zeroes is even and we expect to be able to draw branch cuts between all of those on the unit circle. This leaves the branch point at zero; usually I would conclude that there will be a branch point at infinity, but I'm not sure that holds here.

*If one truncates the series expansion at some $2^n$, then the number of zeroes is finite. Most of them are in the neighborhood of the unit circle as would be expected. However, there is also a zero for some $z\in(0,1)$ (as follows from an intermediate value argument.) What happens to this negative real zero when we go to the infinite series?

A: This is an extended comment to the answer of Semiclassical.

The series
$$f(z)=\sum_{k=0}^{\infty}z^{2^k}$$
is a canonical example of a function having natural boundary (the so-called lacunary series). Here the natural boundary is given by the unit circle $|z|=1$, inside which the series is absolutely convergent. 
A more general statement is the Ostrowski–Hadamard gap theorem:

Let $\{p_k\}$ be a sequence of positive integers satisfying $\lambda p_k<p_{k+1}$ for any $k\in\mathbb{N}$ with some $\lambda>1$, and let $\{a_k\}$ be a sequence of complex numbers such that the radius of convergence of the series $f(z)=\displaystyle\sum_{k\in\mathbb{N}}a_kz^{p_k}$ is equal to $1$. Then no point on the unit circle $|z|=1$ is a regular point of $f(z)$.

The answers to these MO questions (1, 2) may also be of interest.
