Is it possible to simplify the below sum? or maybe find a nice closed form (or asymptotics)?
$$2\sum_{n = 0}^{m - 1}(-1)^n\zeta(4m - 2n + 1)\zeta(2n + 2)$$
I am very much aware about even values of Riemann zeta (Euler's formula in terms of Bernoulli numbers) and that would definitely give us a nice closed form for the term $\zeta(2n + 2)$, but what about the term $\zeta(4m - 2n + 1)$? This seems too complicated to simplify simply because we know very little about the odd values of Riemann zeta. However I'm not familiar with other closed forms for $\zeta(2n + 1)$, so I cannot say for sure that the sum can't be simplified more.
Out of curiosity: Can we simplify or find a closed form for the product of even zeta and odd zeta that is $\zeta(2n + 1)\zeta(2n)$? Or can we get a nice closed form for $\zeta(4m + 3)$?
Thanks.
EDIT $\textbf{1}$: @Claude Leibovici, in his answer claims that
If $m$ is even and large then
$$ 2\sum_{n = 0}^{m - 1}(-1)^n\,\zeta(4m - 2n + 1)\,\zeta(2n + 2) \approx \pi \coth (\pi )-2 \tag{1}\label{1}$$
and if $m$ is odd then $$2\sum_{n = 0}^{m - 1}(-1)^n\,\zeta(4m - 2n + 1)\,\zeta(2n + 2) \approx \pi \coth (\pi ) \tag{2}\label{2}$$
However, I have no idea on how to prove $\eqref{1}$ and $\eqref{2}$, so proving $\eqref{1}$ and $\eqref{2}$ (probably using some elementary tools) is sufficient for answering the question completely. Thanks.
As is shown in this answer $$ \sum_{k=1}^\infty\zeta(2k)x^{2k}=\frac12(1-\pi x\cot(\pi x))\tag1 $$ substituting $x\mapsto ix$, we get $$ \sum_{k=1}^\infty(-1)^{k-1}\zeta(2k)x^{2k}=\frac12(\pi x\coth(\pi x)-1)\tag2 $$ and therefore, $$ \sum_{k=1}^\infty(-1)^{k-1}(\zeta(2k)-1)x^{2k}=\frac12(\pi x\coth(\pi x)-1)-\frac1{1+x^2}\tag3 $$ Taking the limit of $(3)$ as $x\to1$, we get $$ \sum_{k=1}^\infty(-1)^{k-1}(\zeta(2k)-1)=\frac12(\pi\coth(\pi)-2)\tag4 $$ Since $n\le m-1$, we have that $4m-2n+1\ge2m+3$. This means that $\zeta(4m-2n+1)=1+O\!\left(4^{-m}\right)$ $$ \begin{align} &2\sum_{n=0}^{m-1}(-1)^n\color{#C00}{\zeta(4m-2n+1)}\color{#090}{\zeta(2n+2)}\tag{5a}\\ &=2\color{#C00}{\left(1+O\!\left(4^{-m}\right)\right)}\sum_{n=1}^m(-1)^{n-1}(\color{#090}{\zeta(2n)}-1)+2(m\bmod2)\tag{5b}\\ &=2\left(1+O\!\left(4^{-m}\right)\right)\frac12(\pi\coth(\pi)-2)+2(m\bmod2)\tag{5c}\\[6pt] &=\bbox[5px,border:2px solid #C0A000]{\pi\coth(\pi)-2+2(m\bmod2)+O\!\left(4^{-m}\right)}\tag{5d} \end{align} $$ Explanation:
$\text{(5b)}$: apply $\zeta(4m-2n+1)=1+O\!\left(4^{-m}\right)$
$\phantom{\text{(5b):}}$ then substitute $n\mapsto n-1$
$\phantom{\text{(5b):}}$ then compensate for $2\sum_{n=1}^m(-1)^n$ with $2(m\bmod2)$
$\text{(5c)}$: apply $(4)$ noting that $\zeta(2n)-1=O\!\left(4^{-n}\right)$
$\text{(5d)}$: rearrange to make things look nicer
Note that since $n\le m-1$, $$ \zeta(4m-2n+1)-\zeta(4m-2n+2)=O\!\left(4^{-m}\right)\tag6 $$ and so $\text{(5d)}$ also works for $\zeta(4m-2n+2)$ in place of $\zeta(4m-2n+1)$.
Graphical Verification
Plotting $4^m$ times the absolute error of the approximation given in $\text{(5d)}$ shows that it is bounded and appears to tend to a limit around $0.6$.