The series $$\sum_{k=1}^\infty \zeta (2k)z^{2k}$$ (where $\zeta$ is the zeta function) has a non-zero radius of convergence (in $z$). This led me to the following question:
Let $$C_k=\sum_{m=1}^\infty\sum_{n=0}^\infty \frac{1}{(m^2+n^2)^{2k}}$$ where $k$ is an integer $\ge 1$. Does $$\sum_{k=1}^\infty C_k z^{4k}$$ have a non-zero radius of convergence as well?
In fact, $C_k$ is "almost" the Eisenstein series (https://en.wikipedia.org/wiki/Eisenstein_series) of weight $4k$ evaluated at $i$ (known to converge absolutely), with the summand replaced by absolute value: $$C_k=\sum_{m=1}^\infty\sum_{n=0}^\infty \frac{1}{|m+ni|^{4k}}.$$
I tried to build on the analogy with the zeta function. Using the fact that $m^2+n^2\ge n^2$, I obtained $$C_k\le \zeta (4k)+\sum_{m=1}^\infty \zeta (4k)$$ but, unfortunately, this doesn't lead anywhere.