Recently i read that the space of completely holomorphic (also at the cusps) modular forms $M_k(\Gamma(N))$ possesses a basis having Fourier coefficients in $\mathbb{Z}[\zeta_N]$ where $\zeta_N = e^{2 \pi i / N}$.
Can somebody point out a reference for this?
I already know the following things: At least for $k \geq 2$, $S_k(\Gamma(N))$ -- the subspace of cusp forms -- possesses a basis having Fourier coefficients in $\mathbb{Z}$ (see Shimura, Thm 3.52). What is missing is the Eisenstein series $G^{v}$ (see Diamond/Shurman, Thm 4.2.3). All the Fourier coiefficients except the first one do indeed lie inside $\mathbb{Z}[\zeta_N]$ (up to a constant in $\mathbb{Q}$) but the constant term of the Eienstein series is (in the case that $v_1 \equiv 0 \mod N$) the term
$\sum_{n \in \mathbb{Z} \setminus \{0\}, n \equiv v_1 \mod N} \frac{1}{n^k}$
This is the Hurwitz Zeta Function up to the term $N^{-k}$. The question here is: is this value in $\mathbb{Z}[\zeta_N]$ (up to some denominator) or is there a completely different way to see that such a basis with Fourier coeffs in $\mathbb{Z}[\zeta_N]$ exists?
Note that i am aware of this post: Is there a Miller basis for M_k(N)? but i was not able to locate the result in these books.
best and thanks!
Fabian Werner