Certain convolution formulas for multiple series (Q2441401)

For functions \(\beta,f:\mathbb{Z}^r\setminus\{\mathbf{0}\}\to\mathbb{C}\) and \(\theta\in (0,1)\cap\mathbb{Q}\), let \[ Z(s;\beta,f) =\sum_{\mathbf{m}\in\mathbb{Z}^r}\frac{\beta(\mathbf{m})}{f(\mathbf{m})^s} \] and \[ \begin{aligned} Z_2(&s_1,s_2,s_3,s_4;\beta,f,\theta)\\ &=\sum_{\mathbf{m}\in\mathbb{Z}^r}\sum_{\mathbf{n}\in\mathbb{Z}^r} \frac{\beta(\mathbf{m})\beta(\mathbf{n})} {f(\mathbf{m})^{s_1}f(\mathbf{n})^{s_2}\left\{f(\mathbf{m})^{2s_3}-2\cos(\theta\pi)f(\mathbf{m})^{s_3}f(\mathbf{n})^{s_4}+f(\mathbf{n})^{2s_4}\right\}}. \end{aligned} \] Here, \(\sum^{'}_{\mathbf{m}\in\mathbb{Z}^r}\) means that the summation is taken over \(\mathbb{Z}^r\setminus\{\mathbf{0}\}\). Notice that \(Z(s;\beta,f)\) includes the multiple zeta functions of Euler-Zagier type, of Barnes type and the Epstein zeta functions. In this paper, under some analytic assumptions, the author establishes a functional relation between \(Z(s;\beta,f)\) and \(Z_2(s_1,s_2,s_3,s_4;\beta,f,\theta)\). More precisely, for non-negative integer \(N\), positive integers \(a,b\) satisfying \((a,b)=1\) with \(a<2b\) and \(s\in\mathbb{C}\) with \(\mathrm{Re}(s)>R\), it holds that \[ \begin{aligned} & Z_2\left(s-2d,s+4dbN,d,d;\beta,f,\frac{a}{2b}\right)\\ &\;\;\;-\cos\left(\frac{a}{2b}\pi\right)Z_2\left(s-d,s+4dbN-d,d,d;\beta,f,\frac{a}{2b}\right)\\ &=\frac{1}{2}\sum^{4bN}_{\mu=0}e^{-\pi ia\mu/2b}Z(s+d\mu;\beta,f)Z(s+4dbN-d\mu;\beta,f), \end{aligned} \] where the real number \(R>0\) and the positive integer \(d\) come from the above assumptions. As concrete examples, we have \[ \begin{aligned} \sum_{m,n\geq 1}\sum_{k,l\geq 1}\frac{1}{m^2k^2(m+n)(k+l)^3\left\{(m+n)^2+(k+l)^2\right\}} &=\frac{1}{8}\left(\pi^2\zeta(3)-11\zeta(5)\right)^2,\\ \sum_{m,n\geq 1}\sum_{k,l\geq 1}\frac{1}{(m+ni)^2(k+li)^4\left\{(m+ni)^2+(k+li)^2\right\}} &=\frac{\varpi^8}{7200}-\frac{\varpi^4\pi^4}{5400}+\frac{\pi^8}{16200} \end{aligned} \] and \[ \begin{aligned}\sum_{(m,n)\in\mathbb{Z}^2}\;& \sum_{(k,l)\in\mathbb{Z}^2} \frac{1}{(m^2+mn+n^2)(k^2+kl+l^2)^3\left\{(m^2+mn+n^2)^2+(k^2+kl+l^2)^2\right\}}\\ &=\frac{32\pi^6}{2187}\zeta(3)^2, \end{aligned} \] where \(\zeta(s)\) is the Riemann zeta function and \(\varpi=2.622\ldots\) is the lemniscate constant. The proof is given by considering a convolution type double integral of \(Z(s;\beta,f)\) via a generalization of the method established by \textit{S. Kanemitsu} et al. [J. Math. Soc. Japan 57, 1167--1177 (2005; Zbl 1092.11037)] and its slight modification by the author [Bull. Aust. Math. Soc. 86, 461--472 (2012; Zbl 1277.11088)].