Numéricamente obtuve estas raíces con partes imaginarias crecientes :
$-4.3409477472 + i\ 4.4770216332$
$-4.2880579538 + i\ 13.633957460$
$-4.3851693625 + i\ 22.647453867$
$-4.2763847378 + i\ 31.711259005$
$-4.3658043306 + i\ 40.8352457753$
$-4.3147004415 + i\ 49.8019753881$
$-4.3083714889 + i\ 58.9723863592$
...
La regularidad y las expresiones analíticas para las raíces exactas parecen estar comprometidas.
EDIT2: Como se trata de una serie de Dirichlet (con coeficientes $c_n=e^{-n^2}$ ) y puesto que hay convergencia absoluta en todas partes (por la prueba M de Weierstrass, como ha observado Eric Naslund) podemos utilizar el propiedades analíticas de las series de Dirichlet para concluir que la abscisa de convergencia de la serie de Dirichlet es $-\infty$ para que $G$ es analítica en todo el plano. Aplicando la Teorema de factorización de Weierstrass deducimos que $G(z)$ puede escribirse como un producto que incluye sus raíces.
Los valores aproximados fueron muy bien presentados por David Speyer, así que vamos a compararlos con los valores reales. He tabulado las distancias entre las raíces reales y las aproximadas para obtener : $[0.056814, 0.054425, 0.058860, 0.053906, 0.057958, 0.055621, 0.055335]$
En cuanto a los argumentos (de las raíces reales) la diferencia entre los argumentos consecutivos $\mod (2 \pi)$ son : $[2.5436, 2.6434, 2.6083, 2.5665, 2.6760, 2.5343]$
Tanto la distancia como la rotación parecen regulares. Esto podría indicar que las raíces exactas "giran regularmente" alrededor de las raíces aproximadas como $k$ crece. En efecto, una "aproximación de segundo orden" para la raíz de índice k es :
$$\frac{-3}{\ln(2)} + \frac{(2k+1)\pi}{\ln(2)} i - \alpha e^{i(2k+1)\beta} \approx - 4.328 + (2k+1) 4.532 i - 0.05621 e^{-4.9793i(2k+1)}$$ $$\mathrm{with}\ \alpha=\frac{e^{3\frac{ln(3)}{\ln(2)}-8}}{\ln(2)},\ \beta=-\pi \frac{\ln(3)}{\ln(2)}\ \mathrm{and}\ k \in \mathbb{Z}.$$
Para conseguirlo he añadido un término a la aproximación de David: $g_3(z) := e^{-1} + e^{-4}/2^z + e^{-9}/3^z$
multiplicado por $e$ , considerado el $e^{1-9}/3^z$ término una "pequeña perturbación" de $1$ para obtener : $-e^{e^{-8}/3^z}\approx e^{-3-z\ln(2)}$ . Tomando los logaritmos de ambos lados y sustituyendo la z de la izquierda por la solución de "primer orden" se obtiene $z\approx \left(-3+(2k+1)i\pi-e^{-8-\ln(3)(-3+(2k+1)i\pi)/\ln(2)}\right)/\ln(2).$
Con esta nueva aproximación la distancia mediana es ahora $\approx 0.0023$ pero me abstendré de la búsqueda de epiciclos en esta nueva astronomía que todavía requiere un Kepler...
RAÍCES EXACTAS - Utilizando el mismo método (con los términos exactos a la izquierda en lugar de $-e^{e^{-8}/3^z}$ ) podemos encontrar las raíces exactas mediante iteraciones. Definamos : $$f_n(z)=1+e\sum_{m=3}^{n+1} \frac{e^{-m^2}}{m^z}\ \mathrm{and} \ r_1(k)=\frac{-3+(2 k+1) i \pi}{\ln(2)}$$
entonces la raíz correspondiente a $k$ se obtiene mediante iteraciones de : $$\ r_n(k)=r_1(k)-\frac{\ln(f_n(r_{n-1}(k))}{\ln(2)}$$
Ejemplo para $k=6$ :
$r_1(6)= -4.328085123 + i\ 58.92068184$
$r_2(6)= -4.310848945 + i\ 58.97454441$
$r_3(6)= -4.308195489 + i\ 58.97248418$
$r_4(6)= -4.308369505 + i\ 58.97237420$
$r_5(6)= -4.308372230 + i\ 58.97238650$
$r_6(6)= -4.308371466 + i\ 58.97238640$
$r_7(6)= -4.308371487 + i\ 58.97238636$
$r_8(6)= -4.308371489 + i\ 58.97238636$
(la distancia del valor real es superior a $16$ veces menor en cada iteración, lo mismo ocurre con $k=0$ ).
¡Buen problema!
OTRAS CAPAS - Investiguemos la "segunda capa" detectada por deoxygerbe con primeras raíces en :
$-12.3840905498065732 + i\ 7.69101465838042938$
$-12.2523952898036964 + i\ 23.2196214862555024$
$-12.3153823799387642 + i\ 38.8071497678513693$
La "primera capa" fue generada por las soluciones de $0=e^{-1^2}/1^z+e^{-2^2}/2^z$ (primer+segundo término de la serie).
La "segunda capa" parece haber sido generada por las soluciones de $0=e^{-2^2}/2^z+e^{-3^2}/3^z$ (segundo+tercer término de la serie). Las soluciones pueden obtenerse a partir de las anteriores con la sustitución de $-3=1^2-2^2$ por $-5=2^2-3^2$ y de $\ln(2)=\ln(2/1)$ por $\ln(3/2)$ para obtener : $$r^{(2)}(k)=\frac{-5+(2k+1)i\pi}{\ln(\frac32)}$$
He utilizado la notación $r^{(2)}_k$ con el pensamiento (¡confortado por su clara imagen!) de que deben existir otras capas cerca de : $$r^{(L)}(k)=\frac{-(2L+1)+(2k+1)i\pi}{\ln(\frac{L+1}{L})}$$
En efecto, una solución cercana a $r^{(3)}(0)$ (la "tercera capa" en $x\approx -24.33$ ) es : $-24.400057077391736+ i\ 10.878966908988947$ .
Además, la abscisa de la "cuarta capa" cerca de $-40.33$ con ordenadas cerca de $14.0$ y $42.2$ parecen corresponder bastante bien a su imagen.
De todos modos vamos a tratar de manejar todas las capas a la vez y definir : $$f^{(L)}_n(z)=e^{L^2}L^z\left(\sum_{m=1}^{n+1} \frac{e^{-m^2}} {m^z}- \frac{e^{-(L+1)^2}}{(L+1)^z}\right)$$ $$r^{(L)}_1(k)=\frac{-(2 L+1)+(2 k+1) i \pi}{\ln(\frac{L+1}{L})}$$
entonces la raíz correspondiente a $k$ en la capa $L$ debe obtenerse mediante iteraciones de : $$\ r^{(L)}_n(k)=r^{(L)}_1(k)-\frac{\ln(f^{(L)}_n(r^{(L)}_{n-1}(k))}{\ln(\frac{L+1}{L})}$$
Esto parece proporcionar todas las raíces "visibles", pero no excluyo la existencia de otras familias (que implican, por ejemplo, el primer y el tercer término), después de todo, me perdí una infinidad de capas... :-)