I was reading Irresistible Integrals by Victor H. Moll, where I encountered the following Taylor series expansion of $\ln\Gamma(1+x)$ $$ \ln\Gamma(1+x)=-\gamma x + \sum_{k=2}^{\infty}\dfrac{(-1)^k \zeta (k)}{k}x^k. $$
I am looking for the Taylor series for $\ln^2\Gamma(1+x)$. I guess that maybe I can just square the result I already have, like this $$ \begin{align} \ln^2\Gamma(1+x) &= \left(-\gamma x +\sum_{k=2}^{\infty}\dfrac{(-1)^k \zeta (k)}{k}x^k\right)^2\\ &= \gamma^2x^2 -2\gamma \sum_{k=2}^\infty \dfrac{(-1)^k\zeta(k)x^{k+1}}{k}+ \left(\sum_{k=2}^{\infty}\dfrac{(-1)^k \zeta (k)}{k}x^k\right)^2. \end{align} $$
If this is true, how can I square the sum? I think Cauchy product can help, but that will be very complicated and requires more manipulations which are way too much for me as I am new to the manipulations of zeta function.
Any help?